# The Near-Surface Small-Scale Spatial and Temporal Variability of Sensible and Latent Heat Exchange in the Svalbard Region: A Case Study

Authors: Georg Jocher; Florian Karner; Christoph Ritter; Roland Neuber; Klaus Dethloff; Friedrich Obleitner; Joachim Reuder; Thomas Foken
Journal: ISRN Meteorology
 (2012)

---

## Abstract

We present data from two eddy covariance systems for determining the sensible and latent heat flux and the shear stress near the earth's surface. One measures continuously since September 2010 close to the village Ny-Ålesund, Svalbard, the other one was located on Kongsvegen glacier for a short period in April 2011. Two examples for small-scale variability are discussed: near surface external gravity waves associated with katabatic wind from the Broggerbreen glacier located a few kilometres southwestern of Ny-Ålesund, and an episode when the two eddy systems at the different measurement sites captured very different conditions at the same time. In case of gravity wave motion the eddy covariance method results in fictitious positive fluxes due to the strong correlation between temperature and vertical wind, which has to be considered carefully.  The comparison between the two sites showed that generally the dynamical processes near the Earth’s surface for Ny-Ålesund and the Kongsvegen glacier are different and local. But there are also cases of synchronization due to synoptic influences, and then the same processes are visible at the two measurement sites. Both examples show that the boundary layer in Ny-Ålesund is not only affected by the main glaciers in the east of the village, but also by other orographical characteristics and synoptic issues. Therefore, the only meaningful way to deal with point measurements is to consider them in context with the surrounding orography and the general meteorological conditions.


---

## Body

## 1. Introduction

Dynamical processes near the Earth’s surface are important and affect the whole atmospheric boundary layer. In Arctic environments, near-surface temperature inversions are prominent phenomena which essentially influence surface energy exchanges during stable stratification and drive important local and regional atmospheric circulations. A proper treatment of such phenomena in atmospheric models is still a problem, which in part is due to improper horizontal resolution (about 25–50 km). On the other hand, existing subgrid-scale parameterizations have still problems to describe these processes in a correct way [1]. Further, there is also a need of relevant data for input and validation of model results, which are still sparse in Arctic environments particularly. So, it is essential to get more data of the near surface processes to complement and amplify the existing work on this field at Svalbard (e.g., [2] or [3]). This is of most importance for a better understanding of the underlying physics, for the appreciation of the interaction of these processes with the whole Arctic atmospheric boundary layer, and finally, to take account of these processes in numerical weather prediction and climate models as much as possible. Another approach to solve the problem of the linkage between small-scale and large-scale processes is, for example, the Large-Eddy simulation [4].Boundary layer measurements started in Ny-Ålesund several years ago. Beine et al. [7] analysed the wind profile and found prevailing winds from south-east which were interpreted as katabatic outflows from the glaciers to the east, namely, the Kongsvegen and the Kronebreen glacier. Argentini et al. [8] presented data on wind and thermal structure of the boundary layer and found stably stratified conditions with gravity waves excited at low wind speeds. Forced convection due to mechanical turbulence was found at higher wind speeds. Recently, Vihma et al. [9] investigated air temperature and humidity inversions and low-level jets over two Svalbard fjords, one was the Kongsfjorden, where Ny-Ålesund is located. These examples and also the recent construction of the 32 m high Italian Amundsen Nobile Climate Change tower [10] near Ny-Ålesund demonstrate the significance of boundary layer measurements to understand the processes in the whole atmosphere as much as possible.The mentioned investigations already showed a high temporal and spatial variability of the near-surface exchange processes. The reason for this high variability can probably be found in the complex orography around the site, including the coast. Thus, Ny-Ålesund is indirectly influenced by the Gulf Stream and resembles more tundra climate, hence, results concerning the boundary layer processes derived at that place cannot be simply generalized for the Arctic as a whole. Therefore, the site of Ny-Ålesund is in actual fact not well suited for Arctic boundary layer measurements, in spite of its easy accessibility, the disparities in space and time are too pronounced. A more detailed understanding on how the synoptic meteorological situation interacts with orography and small-scale processes is, therefore, urgently needed.In this work data from two eddy covariance stations, one close to the village and a second on Kongsvegen glacier are analysed. A case of gravity wave motion at the measurement site near Ny-Ålesund, triggered by katabatic outflows of the Broggerbreen glacier and the orography was captured as an example for the high temporal near-surface process variability. Moreover, cases of similar and different conditions in comparison of both measurement sites will be discussed to analyse the small-scale spatial variability in this region.The presented periods (November 2010 for discussing the temporal variability, April 2011 for the spatial variability) were chosen due to their representativeness of the mentioned issues.


## 2. Sites and Methods

### 2.1. Eddy Covariance System Ny-Ålesund

Ny-Ålesund on Svalbard (Figure1) is a centre for different polar research institutions, also for the Alfred Wegener Institute AWI (AWIPEV research station). The eddy covariance system  (N 78° 55.287′, E 011° 54.851′)  is located about 300 m south of the meteorological observatory of Ny-Ålesund about half-way between the observatory and the Zeppelin mountain on the rather flat bare tundra. The position has mainly been chosen because of the quite homogeneous footprint area around.The Kongsfjord region on Svalbard (upper left side), including both measurement sites Ny-Ålesund and Kongsvegen glacier (marked with arrows). The overview picture is taken from [5]. The upper right picture [6] shows the Kongsvegen measurement site. Beneath a topographical map of the Kongsfjord region (http://www.toposvalbard.npolar.no/) is shown. The arrows in the map mark the dominant wind directions in Ny-Ålesund, most often are 1 and 2 (eastern directions). 3 marks the direction downward the Broggerbreen glacier, 4 the free ocean, and the yellow dot the eddy covariance system on the Kongsvegen glacier.
(a)
(b)The snow height during the relevant periods was 0.4 m in April 2011 (effective measurement height 1.7 m) and 0.05 m in November 2010 (effective measurement height 2.05 m).


### 2.2. Eddy Covariance System Kongsvegen Glacier

The distance between the eddy covariance complex near Ny-Ålesund and the measurement system on the Kongsvegen glacier is about 15 km. Kongsvegen is one of the few Svalbard glaciers where long-term meteorological and mass balance measurements are performed. The glacier drains into Kongsfjorden and originates at a saddle connecting to the south flowing Sveabreen (800 m a.s.l.). It covers an area of 101 km², and has a length of 26 km, slopes are gentle and have a north-westerly aspect. The turbulence measurements were performed at the tongue of the glacier at an elevation of 162 m a.s.l. (N 78° 50.725′, E 012° 40.106′). The local slope of the surface is oriented towards NW, that is, towards the Kongsfjorden coast. The surface conditions at this site are characterized by snow during winter and bare ice during summer. Transitions, the changes from snow to ice and vice versa, usually occur in May and September. At the time of the turbulence measurements considered here (April 2011), the surface was still covered by snow. The measured snow depth was 67 cm (referring to the underlying glacier ice) and the average density of the snow pack was 331 kg m−3. The surface roughness may be characterized as slightly undulating sastrugi structures upon which small ripples were present.So far, micrometeorological or turbulence measurements have not been performed at Kongsvegen glacier except of a study performed about 150 m higher up at the glacier during summer 2008 [6]. These data revealed temperature inversions in excess of 10°C within the lowest 15 m above the surface and a katabatic wind maximum occurring frequently 5–10 m above the surface. The katabatic wind speeds, therefore, reached a maximum velocity of 4 ms−1, the absolute wind maximum downward the glacier as combination of katabatic and synoptic forcing went up to 15 ms−1. Turbulence characteristics were not yet evaluated.Turbulence measurements cannot routinely be performed at the glacier due to problems concerning for example appropriate power supply, proper alignment of instruments at the moving glacier surface, accretion of snow, and rime or data retrieval. The data considered in this study are based on measurements during the period from 6 April to 1 May 2011, contemporary to the Pan-Arctic Measurements and Arctic Regional Climate Model Simulation Project (PAMARCMiP) 2011, when the station was visited at least twice a week for proper maintenance. The equipment was mounted at a tripod which was frozen into the ice and stabilized by guy wires. To minimize disturbances by the mast and supports, the sensors were mounted towards north-east that is, perpendicular with respect to the prevailing wind direction (down slope).


### 2.3. Instrumentation and Data Processing

The measurements employ a sonic anemometer CSAT 3 for determining the high-frequency fluctuations of the wind components and the temperature and an LiCor 7500 for the humidity fluctuations for both stations. All the values were sampled with 20 Hz frequency. The sensors were mounted at a height of 2.1 m above ground, the data postprocessing was made with the internationally compared eddy covariance software TK 3 [11, 12]. This software includes the option to use spike detection, the conversion from the buoyancy flux to the “real” sensible heat flux and the rotation of the measured data in the actual streaming field (double rotation was used here) which was applied for all calculations in this context.Most important for a meaningful use of the eddy covariance method [13] is the exact application of all corrections and quality tests according to Foken et al. [14] and Foken et al. [15]. The eddy covariance method needs horizontally homogeneous surfaces, steady-state conditions, and ideally well-developed turbulence. Obviously these mentioned requirements pose severe limitations especially in Arctic regions frequently. Svalbard is not horizontally homogeneous, well-developed turbulence is quite rare in the Arctic and steady-state conditions not the normal case. This has to be considered carefully in the data handling, in Lüers and Bareiss [2], for example, are some issues of that problem discussed.As background, the following quantities were collected with a lower frequency (1 Hz, second mast) at the Ny-Ålesund site: wind velocity, wind direction, temperature, air pressure, absolute and relative humidity all at 2 m height, as well as the surface temperature and snow height during winter. Further, there is an additional temperature profile measurement on this mast (3 instruments at 1.5, 1, and 0.5 m height) and a small soil profile next to the mast with measurements of temperature in 1 and 10 cm soil depth, a heat flux plate in 10 cm depth and a TDR (Time Domain Reflection) instrument for measuring the soil humidity.


### 2.4. The Meteorological Background

#### 2.4.1. General Meteorological Conditions at the Two Measurements Sites

Standard meteorological measurements were performed using automatic weather stations at the two turbulence measurement sites including, for example, the radiation components and surface height changes. These data allow the documentation of the general conditions during the investigation period. Further, these studies are valuable with respect to the interpretation and spatial representativeness of the turbulence data.Considering mean annual values (May 2010 to April 2011), air temperature is about 2 degrees lower at the Kongsvegen site in comparison with Ny-Ålesund (−4.3 versus −6.4°C). This reflects the different elevation as well as different surface conditions. The latter mainly evolves during summer, when there is permanent snow and ice at the glacier, while bare tundra is present at Ny-Ålesund. Mean annual albedo is correspondingly higher at the glacier in the annual mean (0.3 versus 0.6). The atmospheric boundary layer at the glacier is characterized by lower humidity and higher wind speeds than the site Ny-Ålesund.Figure2(b) shows the typical frequency distribution of wind directions observed at the glacier from May 2010 to April 2011. The pattern is characterized by frequent and strong south-easterly winds. Their nature is primarily associated to topographically modified synoptic winds and katabatic forcing. The SE-NW alignment of the surrounding mountain ranges and the presence of Kongsvegpasset play a major role in this context. The secondary peak from north-westerly directions is associated to cold air advection in the rear of low-pressure systems passing the area.Average wind direction frequency distribution at the Ny-Ålesund (a) and Kongsvegen measurement site (b) from May 2010 until April 2011.
(a)
(b)Modelling studies by, for example, Skeie and Gronas [16] or Sandvik and Furevik [17] suggest that glacier winds also influence the wind conditions at Ny-Ålesund located about 15 km downstream of Kongsvegen glacier. This is corroborated by our observational data (Figure 2(a)) which reflect the general orientation of Kongsfjorden orography and the prevalence of winds from the glacierized areas to the southeast. Beine et al. [7] also stated that the predominant flow at Ny-Ålesund is from the east-southeast due to katabatic flow from the Kongsvegen glacier around 15 km to the east of Ny-Ålesund. But this is not always valid, as will be shown later in Section 3. The comparatively frequent winds from south-west are explained as down-slope katabatic winds draining along the slopes of Zeppelin mountain range south of Ny-Ålesund. Svendsen  et al. [18], Argentini et al. [19] and Hartmann et al. [20] showed that the occurrence of these winds is strongly related to the position of the low pressure systems in the area.


#### 2.4.2. Selection of the Periods for the Case Studies

From 25 March to 6 May 2011, the international PAMARCMiP (Polar Airborne Measurements and Arctic Regional Climate Model Simulation Project) 2011 took place. This campaign aimed to measure sea ice thickness and extent next to various key meteorological quantities, including aerosol from a Basler BT 67 aircraft. From 30 April to 5 May, the plane operated from Longyearbyen. Two overflights over Ny-Ålesund could be realised for aerosol, radiation, and boundary layer measurements. The turbulence measurements at Kongsvegen glacier were performed to further support these activities focusing on glacier environments, which are prominent in the Kongsfjord area and Svalbard. Moreover, these data may be used for validation of energy balance calculations at the glacier based on standard parameterizations, and to investigate the glacier’s atmospheric boundary layer and the skill of numerical weather prediction and regional climate model output, respectively.The preliminary analysis of the available turbulence dataset revealed 2 interesting typical features, representing the temporal and spatial variability of micrometeorological processes in this region, which will be investigated in detail in this study. The evaluated data stem from two periods, that is, 6 November 2010 to 9 November 2010 and 8 April 2011 to 9 April 2011. The former addresses gravity wave phenomena observed at Ny-Ålesund and related methodical issues. The latter focuses on the spatial variability of the exchange characteristics comparing contemporary data collected at Ny-Ålesund and Kongsvegen. The two periods were chosen due to their representativeness of the mentioned issues.


#### 2.4.3. Meteorological Background of the Investigated Days in November 2010

(a) 6 November 2010 and 7 November 2010
Figure3 shows a summary of relevant variables in the course of the investigated days. On the left panel, a period with gravity waves is depicted, clearly to see are the inversion and the quite small wind velocities coming from about 210°.
The radiation budget is negative because of the long-wave radiation loss at the surface. For the investigated period with gravity waves, the Brunt-Vaisala frequency is quite high (shown later). This reflects the extreme difference between the surface temperature and the temperature at the measurement height caused by radiative cooling at this time.z/L is negative because of the fictitious positive sensible heat flux, which was measured due to the existence of the gravity waves.Temperature at 2 m and at the surface, humidity, wind velocity (line) and wind direction (dots), radiation budget and the stability parameterz/L for the time period 6 November 2010–9 November 2010 (from the top to the bottom).
(a)
(b)(b) 8 November 2010 and 9 November 2010
On the right panel of Figure3, a period without gravity waves is shown, the inversion is almost gone, the wind speed is larger, and the wind direction is different. The radiation budget is nearly around zero because of clouds (less longwave radiation loss), the stability parameter is zero or positive.


#### 2.4.4. Meteorological Background of the Investigated Days in April 2011

(a) 8 April 2011
On this day it was quite cold (between −13 and −7°C) and dry, the sun was shining and the atmosphere was clear. The wind speed was always lower than 5 ms−1, the direction was oscillating around south. In the evening clouds began to appear. To visualize this, Figure 4 with the wind speed and wind direction is added, produced with data from the BSRN weather mast, a 10 m height mast with different meteorological instruments. Used were the data from the 2 and 10 m wind velocity and direction measurement.
Because of its height, this mast has a larger footprint area, and is, therefore, ideal for reference measurements. Furthermore, it is located nearby Ny-Ålesund, about 200 m in the east of the formerly introduced eddy covariance system.Wind speed and wind direction for 8 April 2011 at the BSRN weather mast, plotted are the values for 2 (green) and 10 (black) m height.
(a)
(b)(b) 9 April 2011
It was a little bit warmer (between −6 and −3°C) with more humidity on this day in comparison to the 8 April 2011. The whole day there was a cloud cover, the wind velocity was still quite low, the wind direction was oscillating. Quite suddenly at about 17 UTC the wind freshed up, the wind direction changed to west-northwest and the temperature went down to about −10°C quite quickly (Figure5).
As an overview, diverse meteorological values for both measurement sites and both days are summarized in Figure6. In the left panel are the values shown for the Ny-Ålesund site, in the right panel are the values shown for the Kongsvegen site. z/L is fluctuating at the Ny-Ålesund site because of the small sensible heat flux. Is the sensible heat flux fluctuating around zero, the Obukhov length is also fluctuating between positive and negative values, and thereby z/L. Detailed analyses are following in Section 3.Wind speed and wind direction for 9 April 2011 at the BSRN weather mast, plotted are the values for 2 (green) and 10 (black) m height.
(a)
(b)Temperature in 2 m and at the surface, relative humidity, wind direction (dots) and wind velocity (line), radiation budget and the stability parameterz/L for the period 8 April and 9 April 2011 (from the top to the bottom).
(a)
(b)


## 2.1. Eddy Covariance System Ny-Ålesund

Ny-Ålesund on Svalbard (Figure1) is a centre for different polar research institutions, also for the Alfred Wegener Institute AWI (AWIPEV research station). The eddy covariance system  (N 78° 55.287′, E 011° 54.851′)  is located about 300 m south of the meteorological observatory of Ny-Ålesund about half-way between the observatory and the Zeppelin mountain on the rather flat bare tundra. The position has mainly been chosen because of the quite homogeneous footprint area around.The Kongsfjord region on Svalbard (upper left side), including both measurement sites Ny-Ålesund and Kongsvegen glacier (marked with arrows). The overview picture is taken from [5]. The upper right picture [6] shows the Kongsvegen measurement site. Beneath a topographical map of the Kongsfjord region (http://www.toposvalbard.npolar.no/) is shown. The arrows in the map mark the dominant wind directions in Ny-Ålesund, most often are 1 and 2 (eastern directions). 3 marks the direction downward the Broggerbreen glacier, 4 the free ocean, and the yellow dot the eddy covariance system on the Kongsvegen glacier.
(a)
(b)The snow height during the relevant periods was 0.4 m in April 2011 (effective measurement height 1.7 m) and 0.05 m in November 2010 (effective measurement height 2.05 m).


## 2.2. Eddy Covariance System Kongsvegen Glacier

The distance between the eddy covariance complex near Ny-Ålesund and the measurement system on the Kongsvegen glacier is about 15 km. Kongsvegen is one of the few Svalbard glaciers where long-term meteorological and mass balance measurements are performed. The glacier drains into Kongsfjorden and originates at a saddle connecting to the south flowing Sveabreen (800 m a.s.l.). It covers an area of 101 km², and has a length of 26 km, slopes are gentle and have a north-westerly aspect. The turbulence measurements were performed at the tongue of the glacier at an elevation of 162 m a.s.l. (N 78° 50.725′, E 012° 40.106′). The local slope of the surface is oriented towards NW, that is, towards the Kongsfjorden coast. The surface conditions at this site are characterized by snow during winter and bare ice during summer. Transitions, the changes from snow to ice and vice versa, usually occur in May and September. At the time of the turbulence measurements considered here (April 2011), the surface was still covered by snow. The measured snow depth was 67 cm (referring to the underlying glacier ice) and the average density of the snow pack was 331 kg m−3. The surface roughness may be characterized as slightly undulating sastrugi structures upon which small ripples were present.So far, micrometeorological or turbulence measurements have not been performed at Kongsvegen glacier except of a study performed about 150 m higher up at the glacier during summer 2008 [6]. These data revealed temperature inversions in excess of 10°C within the lowest 15 m above the surface and a katabatic wind maximum occurring frequently 5–10 m above the surface. The katabatic wind speeds, therefore, reached a maximum velocity of 4 ms−1, the absolute wind maximum downward the glacier as combination of katabatic and synoptic forcing went up to 15 ms−1. Turbulence characteristics were not yet evaluated.Turbulence measurements cannot routinely be performed at the glacier due to problems concerning for example appropriate power supply, proper alignment of instruments at the moving glacier surface, accretion of snow, and rime or data retrieval. The data considered in this study are based on measurements during the period from 6 April to 1 May 2011, contemporary to the Pan-Arctic Measurements and Arctic Regional Climate Model Simulation Project (PAMARCMiP) 2011, when the station was visited at least twice a week for proper maintenance. The equipment was mounted at a tripod which was frozen into the ice and stabilized by guy wires. To minimize disturbances by the mast and supports, the sensors were mounted towards north-east that is, perpendicular with respect to the prevailing wind direction (down slope).


## 2.3. Instrumentation and Data Processing

The measurements employ a sonic anemometer CSAT 3 for determining the high-frequency fluctuations of the wind components and the temperature and an LiCor 7500 for the humidity fluctuations for both stations. All the values were sampled with 20 Hz frequency. The sensors were mounted at a height of 2.1 m above ground, the data postprocessing was made with the internationally compared eddy covariance software TK 3 [11, 12]. This software includes the option to use spike detection, the conversion from the buoyancy flux to the “real” sensible heat flux and the rotation of the measured data in the actual streaming field (double rotation was used here) which was applied for all calculations in this context.Most important for a meaningful use of the eddy covariance method [13] is the exact application of all corrections and quality tests according to Foken et al. [14] and Foken et al. [15]. The eddy covariance method needs horizontally homogeneous surfaces, steady-state conditions, and ideally well-developed turbulence. Obviously these mentioned requirements pose severe limitations especially in Arctic regions frequently. Svalbard is not horizontally homogeneous, well-developed turbulence is quite rare in the Arctic and steady-state conditions not the normal case. This has to be considered carefully in the data handling, in Lüers and Bareiss [2], for example, are some issues of that problem discussed.As background, the following quantities were collected with a lower frequency (1 Hz, second mast) at the Ny-Ålesund site: wind velocity, wind direction, temperature, air pressure, absolute and relative humidity all at 2 m height, as well as the surface temperature and snow height during winter. Further, there is an additional temperature profile measurement on this mast (3 instruments at 1.5, 1, and 0.5 m height) and a small soil profile next to the mast with measurements of temperature in 1 and 10 cm soil depth, a heat flux plate in 10 cm depth and a TDR (Time Domain Reflection) instrument for measuring the soil humidity.


## 2.4. The Meteorological Background

### 2.4.1. General Meteorological Conditions at the Two Measurements Sites

Standard meteorological measurements were performed using automatic weather stations at the two turbulence measurement sites including, for example, the radiation components and surface height changes. These data allow the documentation of the general conditions during the investigation period. Further, these studies are valuable with respect to the interpretation and spatial representativeness of the turbulence data.Considering mean annual values (May 2010 to April 2011), air temperature is about 2 degrees lower at the Kongsvegen site in comparison with Ny-Ålesund (−4.3 versus −6.4°C). This reflects the different elevation as well as different surface conditions. The latter mainly evolves during summer, when there is permanent snow and ice at the glacier, while bare tundra is present at Ny-Ålesund. Mean annual albedo is correspondingly higher at the glacier in the annual mean (0.3 versus 0.6). The atmospheric boundary layer at the glacier is characterized by lower humidity and higher wind speeds than the site Ny-Ålesund.Figure2(b) shows the typical frequency distribution of wind directions observed at the glacier from May 2010 to April 2011. The pattern is characterized by frequent and strong south-easterly winds. Their nature is primarily associated to topographically modified synoptic winds and katabatic forcing. The SE-NW alignment of the surrounding mountain ranges and the presence of Kongsvegpasset play a major role in this context. The secondary peak from north-westerly directions is associated to cold air advection in the rear of low-pressure systems passing the area.Average wind direction frequency distribution at the Ny-Ålesund (a) and Kongsvegen measurement site (b) from May 2010 until April 2011.
(a)
(b)Modelling studies by, for example, Skeie and Gronas [16] or Sandvik and Furevik [17] suggest that glacier winds also influence the wind conditions at Ny-Ålesund located about 15 km downstream of Kongsvegen glacier. This is corroborated by our observational data (Figure 2(a)) which reflect the general orientation of Kongsfjorden orography and the prevalence of winds from the glacierized areas to the southeast. Beine et al. [7] also stated that the predominant flow at Ny-Ålesund is from the east-southeast due to katabatic flow from the Kongsvegen glacier around 15 km to the east of Ny-Ålesund. But this is not always valid, as will be shown later in Section 3. The comparatively frequent winds from south-west are explained as down-slope katabatic winds draining along the slopes of Zeppelin mountain range south of Ny-Ålesund. Svendsen  et al. [18], Argentini et al. [19] and Hartmann et al. [20] showed that the occurrence of these winds is strongly related to the position of the low pressure systems in the area.


### 2.4.2. Selection of the Periods for the Case Studies

From 25 March to 6 May 2011, the international PAMARCMiP (Polar Airborne Measurements and Arctic Regional Climate Model Simulation Project) 2011 took place. This campaign aimed to measure sea ice thickness and extent next to various key meteorological quantities, including aerosol from a Basler BT 67 aircraft. From 30 April to 5 May, the plane operated from Longyearbyen. Two overflights over Ny-Ålesund could be realised for aerosol, radiation, and boundary layer measurements. The turbulence measurements at Kongsvegen glacier were performed to further support these activities focusing on glacier environments, which are prominent in the Kongsfjord area and Svalbard. Moreover, these data may be used for validation of energy balance calculations at the glacier based on standard parameterizations, and to investigate the glacier’s atmospheric boundary layer and the skill of numerical weather prediction and regional climate model output, respectively.The preliminary analysis of the available turbulence dataset revealed 2 interesting typical features, representing the temporal and spatial variability of micrometeorological processes in this region, which will be investigated in detail in this study. The evaluated data stem from two periods, that is, 6 November 2010 to 9 November 2010 and 8 April 2011 to 9 April 2011. The former addresses gravity wave phenomena observed at Ny-Ålesund and related methodical issues. The latter focuses on the spatial variability of the exchange characteristics comparing contemporary data collected at Ny-Ålesund and Kongsvegen. The two periods were chosen due to their representativeness of the mentioned issues.


### 2.4.3. Meteorological Background of the Investigated Days in November 2010

(a) 6 November 2010 and 7 November 2010
Figure3 shows a summary of relevant variables in the course of the investigated days. On the left panel, a period with gravity waves is depicted, clearly to see are the inversion and the quite small wind velocities coming from about 210°.
The radiation budget is negative because of the long-wave radiation loss at the surface. For the investigated period with gravity waves, the Brunt-Vaisala frequency is quite high (shown later). This reflects the extreme difference between the surface temperature and the temperature at the measurement height caused by radiative cooling at this time.z/L is negative because of the fictitious positive sensible heat flux, which was measured due to the existence of the gravity waves.Temperature at 2 m and at the surface, humidity, wind velocity (line) and wind direction (dots), radiation budget and the stability parameterz/L for the time period 6 November 2010–9 November 2010 (from the top to the bottom).
(a)
(b)(b) 8 November 2010 and 9 November 2010
On the right panel of Figure3, a period without gravity waves is shown, the inversion is almost gone, the wind speed is larger, and the wind direction is different. The radiation budget is nearly around zero because of clouds (less longwave radiation loss), the stability parameter is zero or positive.


### 2.4.4. Meteorological Background of the Investigated Days in April 2011

(a) 8 April 2011
On this day it was quite cold (between −13 and −7°C) and dry, the sun was shining and the atmosphere was clear. The wind speed was always lower than 5 ms−1, the direction was oscillating around south. In the evening clouds began to appear. To visualize this, Figure 4 with the wind speed and wind direction is added, produced with data from the BSRN weather mast, a 10 m height mast with different meteorological instruments. Used were the data from the 2 and 10 m wind velocity and direction measurement.
Because of its height, this mast has a larger footprint area, and is, therefore, ideal for reference measurements. Furthermore, it is located nearby Ny-Ålesund, about 200 m in the east of the formerly introduced eddy covariance system.Wind speed and wind direction for 8 April 2011 at the BSRN weather mast, plotted are the values for 2 (green) and 10 (black) m height.
(a)
(b)(b) 9 April 2011
It was a little bit warmer (between −6 and −3°C) with more humidity on this day in comparison to the 8 April 2011. The whole day there was a cloud cover, the wind velocity was still quite low, the wind direction was oscillating. Quite suddenly at about 17 UTC the wind freshed up, the wind direction changed to west-northwest and the temperature went down to about −10°C quite quickly (Figure5).
As an overview, diverse meteorological values for both measurement sites and both days are summarized in Figure6. In the left panel are the values shown for the Ny-Ålesund site, in the right panel are the values shown for the Kongsvegen site. z/L is fluctuating at the Ny-Ålesund site because of the small sensible heat flux. Is the sensible heat flux fluctuating around zero, the Obukhov length is also fluctuating between positive and negative values, and thereby z/L. Detailed analyses are following in Section 3.Wind speed and wind direction for 9 April 2011 at the BSRN weather mast, plotted are the values for 2 (green) and 10 (black) m height.
(a)
(b)Temperature in 2 m and at the surface, relative humidity, wind direction (dots) and wind velocity (line), radiation budget and the stability parameterz/L for the period 8 April and 9 April 2011 (from the top to the bottom).
(a)
(b)


## 2.4.1. General Meteorological Conditions at the Two Measurements Sites

Standard meteorological measurements were performed using automatic weather stations at the two turbulence measurement sites including, for example, the radiation components and surface height changes. These data allow the documentation of the general conditions during the investigation period. Further, these studies are valuable with respect to the interpretation and spatial representativeness of the turbulence data.Considering mean annual values (May 2010 to April 2011), air temperature is about 2 degrees lower at the Kongsvegen site in comparison with Ny-Ålesund (−4.3 versus −6.4°C). This reflects the different elevation as well as different surface conditions. The latter mainly evolves during summer, when there is permanent snow and ice at the glacier, while bare tundra is present at Ny-Ålesund. Mean annual albedo is correspondingly higher at the glacier in the annual mean (0.3 versus 0.6). The atmospheric boundary layer at the glacier is characterized by lower humidity and higher wind speeds than the site Ny-Ålesund.Figure2(b) shows the typical frequency distribution of wind directions observed at the glacier from May 2010 to April 2011. The pattern is characterized by frequent and strong south-easterly winds. Their nature is primarily associated to topographically modified synoptic winds and katabatic forcing. The SE-NW alignment of the surrounding mountain ranges and the presence of Kongsvegpasset play a major role in this context. The secondary peak from north-westerly directions is associated to cold air advection in the rear of low-pressure systems passing the area.Average wind direction frequency distribution at the Ny-Ålesund (a) and Kongsvegen measurement site (b) from May 2010 until April 2011.
(a)
(b)Modelling studies by, for example, Skeie and Gronas [16] or Sandvik and Furevik [17] suggest that glacier winds also influence the wind conditions at Ny-Ålesund located about 15 km downstream of Kongsvegen glacier. This is corroborated by our observational data (Figure 2(a)) which reflect the general orientation of Kongsfjorden orography and the prevalence of winds from the glacierized areas to the southeast. Beine et al. [7] also stated that the predominant flow at Ny-Ålesund is from the east-southeast due to katabatic flow from the Kongsvegen glacier around 15 km to the east of Ny-Ålesund. But this is not always valid, as will be shown later in Section 3. The comparatively frequent winds from south-west are explained as down-slope katabatic winds draining along the slopes of Zeppelin mountain range south of Ny-Ålesund. Svendsen  et al. [18], Argentini et al. [19] and Hartmann et al. [20] showed that the occurrence of these winds is strongly related to the position of the low pressure systems in the area.


## 2.4.2. Selection of the Periods for the Case Studies

From 25 March to 6 May 2011, the international PAMARCMiP (Polar Airborne Measurements and Arctic Regional Climate Model Simulation Project) 2011 took place. This campaign aimed to measure sea ice thickness and extent next to various key meteorological quantities, including aerosol from a Basler BT 67 aircraft. From 30 April to 5 May, the plane operated from Longyearbyen. Two overflights over Ny-Ålesund could be realised for aerosol, radiation, and boundary layer measurements. The turbulence measurements at Kongsvegen glacier were performed to further support these activities focusing on glacier environments, which are prominent in the Kongsfjord area and Svalbard. Moreover, these data may be used for validation of energy balance calculations at the glacier based on standard parameterizations, and to investigate the glacier’s atmospheric boundary layer and the skill of numerical weather prediction and regional climate model output, respectively.The preliminary analysis of the available turbulence dataset revealed 2 interesting typical features, representing the temporal and spatial variability of micrometeorological processes in this region, which will be investigated in detail in this study. The evaluated data stem from two periods, that is, 6 November 2010 to 9 November 2010 and 8 April 2011 to 9 April 2011. The former addresses gravity wave phenomena observed at Ny-Ålesund and related methodical issues. The latter focuses on the spatial variability of the exchange characteristics comparing contemporary data collected at Ny-Ålesund and Kongsvegen. The two periods were chosen due to their representativeness of the mentioned issues.


## 2.4.3. Meteorological Background of the Investigated Days in November 2010

(a) 6 November 2010 and 7 November 2010
Figure3 shows a summary of relevant variables in the course of the investigated days. On the left panel, a period with gravity waves is depicted, clearly to see are the inversion and the quite small wind velocities coming from about 210°.
The radiation budget is negative because of the long-wave radiation loss at the surface. For the investigated period with gravity waves, the Brunt-Vaisala frequency is quite high (shown later). This reflects the extreme difference between the surface temperature and the temperature at the measurement height caused by radiative cooling at this time.z/L is negative because of the fictitious positive sensible heat flux, which was measured due to the existence of the gravity waves.Temperature at 2 m and at the surface, humidity, wind velocity (line) and wind direction (dots), radiation budget and the stability parameterz/L for the time period 6 November 2010–9 November 2010 (from the top to the bottom).
(a)
(b)(b) 8 November 2010 and 9 November 2010
On the right panel of Figure3, a period without gravity waves is shown, the inversion is almost gone, the wind speed is larger, and the wind direction is different. The radiation budget is nearly around zero because of clouds (less longwave radiation loss), the stability parameter is zero or positive.


## 2.4.4. Meteorological Background of the Investigated Days in April 2011

(a) 8 April 2011
On this day it was quite cold (between −13 and −7°C) and dry, the sun was shining and the atmosphere was clear. The wind speed was always lower than 5 ms−1, the direction was oscillating around south. In the evening clouds began to appear. To visualize this, Figure 4 with the wind speed and wind direction is added, produced with data from the BSRN weather mast, a 10 m height mast with different meteorological instruments. Used were the data from the 2 and 10 m wind velocity and direction measurement.
Because of its height, this mast has a larger footprint area, and is, therefore, ideal for reference measurements. Furthermore, it is located nearby Ny-Ålesund, about 200 m in the east of the formerly introduced eddy covariance system.Wind speed and wind direction for 8 April 2011 at the BSRN weather mast, plotted are the values for 2 (green) and 10 (black) m height.
(a)
(b)(b) 9 April 2011
It was a little bit warmer (between −6 and −3°C) with more humidity on this day in comparison to the 8 April 2011. The whole day there was a cloud cover, the wind velocity was still quite low, the wind direction was oscillating. Quite suddenly at about 17 UTC the wind freshed up, the wind direction changed to west-northwest and the temperature went down to about −10°C quite quickly (Figure5).
As an overview, diverse meteorological values for both measurement sites and both days are summarized in Figure6. In the left panel are the values shown for the Ny-Ålesund site, in the right panel are the values shown for the Kongsvegen site. z/L is fluctuating at the Ny-Ålesund site because of the small sensible heat flux. Is the sensible heat flux fluctuating around zero, the Obukhov length is also fluctuating between positive and negative values, and thereby z/L. Detailed analyses are following in Section 3.Wind speed and wind direction for 9 April 2011 at the BSRN weather mast, plotted are the values for 2 (green) and 10 (black) m height.
(a)
(b)Temperature in 2 m and at the surface, relative humidity, wind direction (dots) and wind velocity (line), radiation budget and the stability parameterz/L for the period 8 April and 9 April 2011 (from the top to the bottom).
(a)
(b)


## 3. Results and Discussion

### 3.1. Temporal Variability of Sensible and Latent Heat Exchange in Arctic Regions: Near-Surface Gravity Waves

The occurrence of a shallow inversion layer above the surface, snow, or ice in Polar Regions in about 1 to 2 m height is nowadays known (e.g., [21]), a disturbed temperature profile in the lowest meters of the atmosphere can be caused by strong cooling and by strong heating at the mentioned interface.In both cases of a narrow inversion the layers below and above the change of temperature gradient are mostly decoupled, a problem for the measurement of heat fluxes, because the determined fluxes do not represent the heat fluxes at the surface (Lüers and Bareiss  2010 [22]). For exact conclusions about the vertical energy transfer, the different layers have to be observed separately.If the inversion is caused by radiative cooling (for external gravity waves, the density of the lower layer has to be larger than the density of the layer above the change of temperature gradient), at the interface between the layers below and above the inversion, external gravity waves can develop [23]. External gravity waves are oscillating vertically and propagating horizontally. The range of gravity wave motions is given by the Brunt-Vaisala frequency, which describes the maximum oscillation frequency of a vertically displaced mass element in the vertical column [24].Figure7 shows the potential temperatures at the surface and the measurement height (a) and the Brunt-Vaisala frequencies (b) for the period 6 November 2010 00 UTC to 7 November 2010 24 UTC on the left side, a period with gravity wave developing. As contrast on the right side in the same graphic, the period 8 November 2010 00 UTC to 9 November 2010 24 UTC is plotted, a period mostly without gravity wave developing. Only a short distance in this period shows the preconditions for gravity waves.Potential temperatures ((a), green is the potential temperature at the measurement height, black at the surface) and Brunt-Vaisala frequency (in 1/s, (b)) for the period 6 November 2010, 00–10 November 2010, 24 UTC.
(a)
(b)The discrepancy in the temperature and frequency behaviour is obvious. Per definition gravity waves should be seen in the air pressure developing, the wavelet plot Figure8(a), for which the Mexican hat wavelet was used, approves this clearly for our measurement period. For example, the air pressure at the measurement height for a period from 6 November 2010 (08 to 16 UTC) is pictured, a period where gravity waves were found. The powerful processes in the frequency range of the gravity waves (at about scale range 4) are evident. In contrast, in Figure 8(b), the air pressure for a period from 9 November 2010 (16 to 24 UTC) is plotted, a period where we did not find gravity waves. Figure 8 at the bottom shows as additional information the frequency distribution of the air pressure signal from the investigated period with gravity waves, therefore, the signal mean was eliminated. The external gravity waves frequency range in this case begins at about 2.8*10−4 per second (marked by a line), which means nearly 1 per hour, and continues obviously in the lower frequency range. So, it is to assume that there are probably waves with different frequencies.Wavelet coefficients of the air pressure in 2 m height for the period 6 November 2010, 08–16 UTC (a). Wavelet coefficients of the air pressure in 2 m height for the period 9 November 2010, 16–24 UTC (b). Scale means in this context delta (signal) *c (c is a continuous numbering vector with the length of the investigated signal range). The more powerful the detected signal is the more light it is in the graphic. Note: on the left side the resolution of the pressure axis is 0.1 hPa per line, on the right side 0.5 hPa per line (in the upper part of the graphics). The graphic at the bottom shows the frequency distribution of the investigated air pressure period with gravity waves.
(a)
(b)
(c)Our data suggest the existence of several criteria for the development of gravity waves.First of all, clear sky conditions are needed, the outgoing longwave radiation has to be considerably larger than the incoming longwave radiation that a near surface inversion can develop. The inversion is the main precondition for the development of this kind of external gravity waves. In contrast, clouds prevent the strong radiative cooling of the atmosphere-surface-interface, there is no near-surface inversion and so the formation of near surface external gravity is suppressed. Furthermore, calm air is essential, the wind velocity has to be low (<5 m/s), such that turbulent mixing by wind shear is suppressed. In this local case, the wind has a favoured direction (about 180–260 degrees) during the existence of gravity waves, a hint for orographical stimulated gravity waves at this site. In this direction the Broggerbreen glacier is located a few kilometers away. It is highly probable that katabatic outflow from this glacier triggers the development of the observed gravity waves, further investigations are ongoing and required.It can be assumed that the detected gravity waves are external gravity waves without a significant vertical component, because in 10 m height, the gravity wave motion is almost invisible (Figure11(a)). Internal gravity waves, which are quite well known in the stably stratified (polar) boundary layer (e.g., [25, 26]), could develop components in all space directions.For estimating the synoptic influence on the discussed processes and to prove the potential linkage between the large-scale flow and processes at small scales, Figure12 shows the wind distribution in 850 hPa pressure level for the periods 6 November 2010-7 November 2010, 12-12 UTC (a) and 8 November 2010-9 November 2010, 12-12 UTC (b), taken from the NCEP/NCAR online tool. The pressure level 850 h Pa was chosen, because it represents quite well the common synoptic situation, orographical effects can be probably neglected in this height for this region. Such estimates are important to improve the general knowledge about the dominant processes and more or less detailed done before for other places, like, for example, the southwest of Norway [27].The wind is coming from northerly directions for the whole pictured space of time, a little bit more from west for the period 8 November 2010-9 November 2010, 12-12 UTC. It can be assumed that the wind distribution is not so important in this case. Fundamental in this relation is the humidity as clouds are the dominating factor in the matter of longwave radiation loss.Another fact has to be considered in this case: during periods with gravity waves positive sensible heat fluxes are measured with the eddy covariance method (Figure13(c)). The figure shows that exact at the times with high Brunt Vaisala frequencies and high potential temperature gradient between surface and measurement height, thus in the times when gravity waves occurs, the sensible heat flux shows positive values. It is most important to point out that these positive flux components are no turbulent fluxes, only a fictitious result of the calculation method due to the strong correlation between vertical wind and temperature in case of gravity waves [28]. The best way to get rid of the fictitious fluxes in case of gravity waves is to filter the data before flux calculation considering the Brunt Vaisala frequency to be sure to have the “pure” turbulent signal [29] for calculating the turbulent fluxes.For detailed statements to the wavelength and phase speed of the gravity waves, different point measurements would be needed [25], which were not available for this case study.


### 3.2. Spatial Variability of Sensible and Latent Heat Exchange in Arctic regions: Contemporary Turbulence Measurements at Ny-Ålesund and the Kongsvegen Glacier

#### 3.2.1. A Period with Different Conditions

Figure14 shows the sensible and latent heat fluxes for the period 8 April 2011, 12–24 UTC at the two measurement sites. On the Kongsvegen glacier strong katabatic winds at this period, coming downward the glacier from the east occurred. Probably the general synoptic situation forced the katabatic wind, Figure 9(a) shows that the wind in 850 hPa pressure level came from south-easterly directions for the period 8 April 2011, 18–24 UTC. The shear stress caused a negative sensible heat flux and positive latent heat flux on the glacier. While on the Kongsvegen glacier the fluxes followed this typical behaviour, these values were fluctuating around zero at the eddy covariance measurement site near Ny-Ålesund. A view on the wind direction and wind speed at the two sites makes this different behaviour understandable (Figure 14 at the bottom). In contrast to the Kongsvegen site, the wind velocity at the Ny-Ålesund site is quite low at this period; the wind direction is oscillating in the range between south and west. Apparently, the general synoptic situation plays no role for the flux behaviour near Ny-Ålesund at this time. The mountains in the east of Ny-Ålesund may avoid this. The processes at the two measurement sites are completely different at this period; the katabatic outflow of the Kongsvegen obviously did not reach Ny-Ålesund.Wind distribution for the period 8 April 2011, 18–24 UTC in 850 hPa pressure level (a). Wind distribution for the period 9 April 2011, 18–24 UTC in 850 hPa pressure level (b). Source for both is the online tool of NCEP/NCAR.
(a)
(b)


#### 3.2.2. A period with Quite the Same Conditions

One day later, the processes on the two measurement sites changed completely. Figure10 shows that after a short period of changeover, the fluxes at the two measurement sites showed quite the same behaviour. Looking at the wind direction and wind speed at the relevant time makes this changing behaviour (Figure 10 at the bottom) explainable. If the katabatic wind on the glacier disappeared, both sites are quite comparable at this time.  With the change of the wind direction on 9 April 2011 and, therefore, the change of the flux behaviour at both sites, the flux values are getting more and more equal. Figure 9(b) shows the general synoptic situation for the period 9 April 2011, 18–24 UTC in 850 hPa pressure level. Well to see is the strong wind from westerly directions. The wind can blow upward the fjord and reaches the measurement site Ny-Ålesund and the measurement site on the Kongsvegen glacier in the same way. So, there is a synchronization of the observed processes, the general synoptic situation seems to play an important role.Sensible heat flux, latent heat flux, wind speed, and wind direction for the period 9 April 2011 15:30–24:00 UTC. The green line are values from the Kongsvegen site, the black line from the site near Ny-Ålesund.
(a)
(b)
(c)
(d)Temperature in 2 (green) and 10 (black) m height for 6 November-7 November 2010 (clear conditions, gravity waves can develop; (a)) and for 8 November-9 November 2010 (mostly cloudy, gravity waves are only sporadic able to develop; (b)).
(a)
(b)Wind distribution for the periods 6 November 2010, 12 UTC-7 November 2010, 12 UTC (a) and 8 November 2010, 12 UTC-9 November 2010, 12 UTC (b), both for 850 hPa pressure level. Source is the online tool of NCEP/NCAR.
(a)
(b)Potential temperatures (a), Brunt Vaisala frequencies (b), and sensible heat flux (c) for the period 6 November 2010, 00 UTC–9 November 2010, 24 UTC.
(a)
(b)
(c)Sensible heat flux, latent heat flux, wind speed and wind direction for the period 8 April 2011 12–24 UTC. The green line are values from the Kongsvegen site, the black line from the site near Ny-Ålesund.
(a)
(b)
(c)
(d)For Sections3.2.1 and 3.2.2 differences in the surface conditions were neglected, because there was a closed snow layer at this time also near Ny-Ålesund, for both places the same roughness (0,24 mm; [30]) was assumed.


## 3.1. Temporal Variability of Sensible and Latent Heat Exchange in Arctic Regions: Near-Surface Gravity Waves

The occurrence of a shallow inversion layer above the surface, snow, or ice in Polar Regions in about 1 to 2 m height is nowadays known (e.g., [21]), a disturbed temperature profile in the lowest meters of the atmosphere can be caused by strong cooling and by strong heating at the mentioned interface.In both cases of a narrow inversion the layers below and above the change of temperature gradient are mostly decoupled, a problem for the measurement of heat fluxes, because the determined fluxes do not represent the heat fluxes at the surface (Lüers and Bareiss  2010 [22]). For exact conclusions about the vertical energy transfer, the different layers have to be observed separately.If the inversion is caused by radiative cooling (for external gravity waves, the density of the lower layer has to be larger than the density of the layer above the change of temperature gradient), at the interface between the layers below and above the inversion, external gravity waves can develop [23]. External gravity waves are oscillating vertically and propagating horizontally. The range of gravity wave motions is given by the Brunt-Vaisala frequency, which describes the maximum oscillation frequency of a vertically displaced mass element in the vertical column [24].Figure7 shows the potential temperatures at the surface and the measurement height (a) and the Brunt-Vaisala frequencies (b) for the period 6 November 2010 00 UTC to 7 November 2010 24 UTC on the left side, a period with gravity wave developing. As contrast on the right side in the same graphic, the period 8 November 2010 00 UTC to 9 November 2010 24 UTC is plotted, a period mostly without gravity wave developing. Only a short distance in this period shows the preconditions for gravity waves.Potential temperatures ((a), green is the potential temperature at the measurement height, black at the surface) and Brunt-Vaisala frequency (in 1/s, (b)) for the period 6 November 2010, 00–10 November 2010, 24 UTC.
(a)
(b)The discrepancy in the temperature and frequency behaviour is obvious. Per definition gravity waves should be seen in the air pressure developing, the wavelet plot Figure8(a), for which the Mexican hat wavelet was used, approves this clearly for our measurement period. For example, the air pressure at the measurement height for a period from 6 November 2010 (08 to 16 UTC) is pictured, a period where gravity waves were found. The powerful processes in the frequency range of the gravity waves (at about scale range 4) are evident. In contrast, in Figure 8(b), the air pressure for a period from 9 November 2010 (16 to 24 UTC) is plotted, a period where we did not find gravity waves. Figure 8 at the bottom shows as additional information the frequency distribution of the air pressure signal from the investigated period with gravity waves, therefore, the signal mean was eliminated. The external gravity waves frequency range in this case begins at about 2.8*10−4 per second (marked by a line), which means nearly 1 per hour, and continues obviously in the lower frequency range. So, it is to assume that there are probably waves with different frequencies.Wavelet coefficients of the air pressure in 2 m height for the period 6 November 2010, 08–16 UTC (a). Wavelet coefficients of the air pressure in 2 m height for the period 9 November 2010, 16–24 UTC (b). Scale means in this context delta (signal) *c (c is a continuous numbering vector with the length of the investigated signal range). The more powerful the detected signal is the more light it is in the graphic. Note: on the left side the resolution of the pressure axis is 0.1 hPa per line, on the right side 0.5 hPa per line (in the upper part of the graphics). The graphic at the bottom shows the frequency distribution of the investigated air pressure period with gravity waves.
(a)
(b)
(c)Our data suggest the existence of several criteria for the development of gravity waves.First of all, clear sky conditions are needed, the outgoing longwave radiation has to be considerably larger than the incoming longwave radiation that a near surface inversion can develop. The inversion is the main precondition for the development of this kind of external gravity waves. In contrast, clouds prevent the strong radiative cooling of the atmosphere-surface-interface, there is no near-surface inversion and so the formation of near surface external gravity is suppressed. Furthermore, calm air is essential, the wind velocity has to be low (<5 m/s), such that turbulent mixing by wind shear is suppressed. In this local case, the wind has a favoured direction (about 180–260 degrees) during the existence of gravity waves, a hint for orographical stimulated gravity waves at this site. In this direction the Broggerbreen glacier is located a few kilometers away. It is highly probable that katabatic outflow from this glacier triggers the development of the observed gravity waves, further investigations are ongoing and required.It can be assumed that the detected gravity waves are external gravity waves without a significant vertical component, because in 10 m height, the gravity wave motion is almost invisible (Figure11(a)). Internal gravity waves, which are quite well known in the stably stratified (polar) boundary layer (e.g., [25, 26]), could develop components in all space directions.For estimating the synoptic influence on the discussed processes and to prove the potential linkage between the large-scale flow and processes at small scales, Figure12 shows the wind distribution in 850 hPa pressure level for the periods 6 November 2010-7 November 2010, 12-12 UTC (a) and 8 November 2010-9 November 2010, 12-12 UTC (b), taken from the NCEP/NCAR online tool. The pressure level 850 h Pa was chosen, because it represents quite well the common synoptic situation, orographical effects can be probably neglected in this height for this region. Such estimates are important to improve the general knowledge about the dominant processes and more or less detailed done before for other places, like, for example, the southwest of Norway [27].The wind is coming from northerly directions for the whole pictured space of time, a little bit more from west for the period 8 November 2010-9 November 2010, 12-12 UTC. It can be assumed that the wind distribution is not so important in this case. Fundamental in this relation is the humidity as clouds are the dominating factor in the matter of longwave radiation loss.Another fact has to be considered in this case: during periods with gravity waves positive sensible heat fluxes are measured with the eddy covariance method (Figure13(c)). The figure shows that exact at the times with high Brunt Vaisala frequencies and high potential temperature gradient between surface and measurement height, thus in the times when gravity waves occurs, the sensible heat flux shows positive values. It is most important to point out that these positive flux components are no turbulent fluxes, only a fictitious result of the calculation method due to the strong correlation between vertical wind and temperature in case of gravity waves [28]. The best way to get rid of the fictitious fluxes in case of gravity waves is to filter the data before flux calculation considering the Brunt Vaisala frequency to be sure to have the “pure” turbulent signal [29] for calculating the turbulent fluxes.For detailed statements to the wavelength and phase speed of the gravity waves, different point measurements would be needed [25], which were not available for this case study.


## 3.2. Spatial Variability of Sensible and Latent Heat Exchange in Arctic regions: Contemporary Turbulence Measurements at Ny-Ålesund and the Kongsvegen Glacier

### 3.2.1. A Period with Different Conditions

Figure14 shows the sensible and latent heat fluxes for the period 8 April 2011, 12–24 UTC at the two measurement sites. On the Kongsvegen glacier strong katabatic winds at this period, coming downward the glacier from the east occurred. Probably the general synoptic situation forced the katabatic wind, Figure 9(a) shows that the wind in 850 hPa pressure level came from south-easterly directions for the period 8 April 2011, 18–24 UTC. The shear stress caused a negative sensible heat flux and positive latent heat flux on the glacier. While on the Kongsvegen glacier the fluxes followed this typical behaviour, these values were fluctuating around zero at the eddy covariance measurement site near Ny-Ålesund. A view on the wind direction and wind speed at the two sites makes this different behaviour understandable (Figure 14 at the bottom). In contrast to the Kongsvegen site, the wind velocity at the Ny-Ålesund site is quite low at this period; the wind direction is oscillating in the range between south and west. Apparently, the general synoptic situation plays no role for the flux behaviour near Ny-Ålesund at this time. The mountains in the east of Ny-Ålesund may avoid this. The processes at the two measurement sites are completely different at this period; the katabatic outflow of the Kongsvegen obviously did not reach Ny-Ålesund.Wind distribution for the period 8 April 2011, 18–24 UTC in 850 hPa pressure level (a). Wind distribution for the period 9 April 2011, 18–24 UTC in 850 hPa pressure level (b). Source for both is the online tool of NCEP/NCAR.
(a)
(b)


### 3.2.2. A period with Quite the Same Conditions

One day later, the processes on the two measurement sites changed completely. Figure10 shows that after a short period of changeover, the fluxes at the two measurement sites showed quite the same behaviour. Looking at the wind direction and wind speed at the relevant time makes this changing behaviour (Figure 10 at the bottom) explainable. If the katabatic wind on the glacier disappeared, both sites are quite comparable at this time.  With the change of the wind direction on 9 April 2011 and, therefore, the change of the flux behaviour at both sites, the flux values are getting more and more equal. Figure 9(b) shows the general synoptic situation for the period 9 April 2011, 18–24 UTC in 850 hPa pressure level. Well to see is the strong wind from westerly directions. The wind can blow upward the fjord and reaches the measurement site Ny-Ålesund and the measurement site on the Kongsvegen glacier in the same way. So, there is a synchronization of the observed processes, the general synoptic situation seems to play an important role.Sensible heat flux, latent heat flux, wind speed, and wind direction for the period 9 April 2011 15:30–24:00 UTC. The green line are values from the Kongsvegen site, the black line from the site near Ny-Ålesund.
(a)
(b)
(c)
(d)Temperature in 2 (green) and 10 (black) m height for 6 November-7 November 2010 (clear conditions, gravity waves can develop; (a)) and for 8 November-9 November 2010 (mostly cloudy, gravity waves are only sporadic able to develop; (b)).
(a)
(b)Wind distribution for the periods 6 November 2010, 12 UTC-7 November 2010, 12 UTC (a) and 8 November 2010, 12 UTC-9 November 2010, 12 UTC (b), both for 850 hPa pressure level. Source is the online tool of NCEP/NCAR.
(a)
(b)Potential temperatures (a), Brunt Vaisala frequencies (b), and sensible heat flux (c) for the period 6 November 2010, 00 UTC–9 November 2010, 24 UTC.
(a)
(b)
(c)Sensible heat flux, latent heat flux, wind speed and wind direction for the period 8 April 2011 12–24 UTC. The green line are values from the Kongsvegen site, the black line from the site near Ny-Ålesund.
(a)
(b)
(c)
(d)For Sections3.2.1 and 3.2.2 differences in the surface conditions were neglected, because there was a closed snow layer at this time also near Ny-Ålesund, for both places the same roughness (0,24 mm; [30]) was assumed.


## 3.2.1. A Period with Different Conditions

Figure14 shows the sensible and latent heat fluxes for the period 8 April 2011, 12–24 UTC at the two measurement sites. On the Kongsvegen glacier strong katabatic winds at this period, coming downward the glacier from the east occurred. Probably the general synoptic situation forced the katabatic wind, Figure 9(a) shows that the wind in 850 hPa pressure level came from south-easterly directions for the period 8 April 2011, 18–24 UTC. The shear stress caused a negative sensible heat flux and positive latent heat flux on the glacier. While on the Kongsvegen glacier the fluxes followed this typical behaviour, these values were fluctuating around zero at the eddy covariance measurement site near Ny-Ålesund. A view on the wind direction and wind speed at the two sites makes this different behaviour understandable (Figure 14 at the bottom). In contrast to the Kongsvegen site, the wind velocity at the Ny-Ålesund site is quite low at this period; the wind direction is oscillating in the range between south and west. Apparently, the general synoptic situation plays no role for the flux behaviour near Ny-Ålesund at this time. The mountains in the east of Ny-Ålesund may avoid this. The processes at the two measurement sites are completely different at this period; the katabatic outflow of the Kongsvegen obviously did not reach Ny-Ålesund.Wind distribution for the period 8 April 2011, 18–24 UTC in 850 hPa pressure level (a). Wind distribution for the period 9 April 2011, 18–24 UTC in 850 hPa pressure level (b). Source for both is the online tool of NCEP/NCAR.
(a)
(b)


## 3.2.2. A period with Quite the Same Conditions

One day later, the processes on the two measurement sites changed completely. Figure10 shows that after a short period of changeover, the fluxes at the two measurement sites showed quite the same behaviour. Looking at the wind direction and wind speed at the relevant time makes this changing behaviour (Figure 10 at the bottom) explainable. If the katabatic wind on the glacier disappeared, both sites are quite comparable at this time.  With the change of the wind direction on 9 April 2011 and, therefore, the change of the flux behaviour at both sites, the flux values are getting more and more equal. Figure 9(b) shows the general synoptic situation for the period 9 April 2011, 18–24 UTC in 850 hPa pressure level. Well to see is the strong wind from westerly directions. The wind can blow upward the fjord and reaches the measurement site Ny-Ålesund and the measurement site on the Kongsvegen glacier in the same way. So, there is a synchronization of the observed processes, the general synoptic situation seems to play an important role.Sensible heat flux, latent heat flux, wind speed, and wind direction for the period 9 April 2011 15:30–24:00 UTC. The green line are values from the Kongsvegen site, the black line from the site near Ny-Ålesund.
(a)
(b)
(c)
(d)Temperature in 2 (green) and 10 (black) m height for 6 November-7 November 2010 (clear conditions, gravity waves can develop; (a)) and for 8 November-9 November 2010 (mostly cloudy, gravity waves are only sporadic able to develop; (b)).
(a)
(b)Wind distribution for the periods 6 November 2010, 12 UTC-7 November 2010, 12 UTC (a) and 8 November 2010, 12 UTC-9 November 2010, 12 UTC (b), both for 850 hPa pressure level. Source is the online tool of NCEP/NCAR.
(a)
(b)Potential temperatures (a), Brunt Vaisala frequencies (b), and sensible heat flux (c) for the period 6 November 2010, 00 UTC–9 November 2010, 24 UTC.
(a)
(b)
(c)Sensible heat flux, latent heat flux, wind speed and wind direction for the period 8 April 2011 12–24 UTC. The green line are values from the Kongsvegen site, the black line from the site near Ny-Ålesund.
(a)
(b)
(c)
(d)For Sections3.2.1 and 3.2.2 differences in the surface conditions were neglected, because there was a closed snow layer at this time also near Ny-Ålesund, for both places the same roughness (0,24 mm; [30]) was assumed.


## 4. Conclusions

In this work, two clear proofs for small-scale processes of boundary layer dynamics around Ny-Ålesund, derived from two eddy covariance measurement systems, one on the Kongsvegen glacier, other on a smooth plane close to Ny-Ålesund, are presented.On the one hand the formation of gravity waves during polar night conditions as example for high temporal variability of flux behaviour was detected: strong temperature inversion conditions due to longwave radiation loss and calm wind are the main preconditions for this phenomenon, which is triggered by katabatic outflow from the Broggerbreen glacier in the southwest of Ny-Ålesund and winds from this direction flowing over the orography in this direction. To minimize the influence of temporal variability and to get the pure turbulent signal, it is recommendable for prospective work on this field to filter the raw turbulence data got by an eddy covariance system. In this way, one can get rid of the longwave components caused by gravity waves and orographical effects like secondary circulations and thereby avoid wrong conclusions while interpreting turbulence data.On the other hand drastic differences in flux measurements on different sites at the same time were observed, among other things depending on katabatic outflows from glaciers, as example for high spatial variability of flux behaviour: the best way to deal with the spatial context in future is to provide as much measurement locations as possible and then compare.Based on this work, we propose that Ny-Ålesund is not only influenced by katabatic outflows from Kongsvegen glacier, as was previously suggested [7], but additionally also from the smaller Broggerbreen and the synoptic situation. Depending on the actual main wind direction, the Kongsvegen glacier (wind from south easterly direction) or the Broggerbreen glacier (wind from south to west direction) influence the dynamics in Ny-Ålesund, both in different ways. Therefore, in future all glaciers in Kongsfjord region and the synoptic scale blocking effect by the mountains surrounding Ny-Ålesund should be considered for understanding causes of small-scale dynamics. The small-scale variations in the near-surface behaviour of momentum, sensible and latent heat exchange measured around Ny-Ålesund are obviously an interaction of many different local effects and the synoptic situation. Because of this issue, it is very difficult to make universal statements about the Arctic atmospheric boundary layer at this place.So, the eddy covariance measurements will be continued on the described measurement site near Ny-Ålesund. A radiometer (for measurements of temperature and humidity profiles) and a wind Lidar (for measurements of 3D wind profiles) will support the eddy covariance measurements to achieve a complete as possible view on the conditions of the atmospheric boundary layer and to observe the entrainment behaviour between the atmospheric boundary layer and the free troposphere.To get more knowledge about the near surface temperature and humidity profiles, additional sensors to determine these quantities were installed in 4, 7, and 10 m height on the already established BSRN weather mast, which is not far away from the eddy covariance system. Finally, for 2012, the additional usage of a tethered balloon system is planned.



---

*Source: 357925-2012-07-19.xml*