ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 1 The surface layer observed by a high- resolution sodar at DOME C, Antarctica STEFANIA ARGENTINI1, IGOR PETENKO1, ANGELO VIOLA1, GIANGIUSEPPE MASTRANTONIO1, ILARIA PIETRONI1, GIAMPIETRO CASASANTA1,*, ERIC ARISTIDI2, CHRISTOPHE GHENTON3 1 Istituto di Scienze dell'Atmosfera e del Clima, CNR, Rome, Italy 2Laboratoire Lagrange, UNS/CNRS/OCA, Nice Cedex 2, France 3 Laboratoire de Glaciologie et Géophysique de l’Environnement, UJF-Grenoble 1/CNRS, Gre- noble, France Abstract A one-year field experiment started on December 2011 at the French - Italian station of Concordia at Dome C, East Antarctic Plateau. The objective of the experiment was the study of the surface layer turbu- lent processes under stable/very stable stratifications, and the mechanisms leading to the formation of the warming events. A sodar was improved to achieve the vertical/temporal resolution needed to study these processes. The system, named surface layer sodar (SL-sodar), may operate both in high vertical resolution (low range) and low vertical resolution (high range) modes. SL-sodar observations were complemented with in situ turbulence and radiation measurements. A few preliminary results, concerning the standard summer diurnal cycle, a summer warming event, and unusually high frequency boundary layer atmos- pheric gravity waves are presented. I. INTRODUCTION t Dome C, light wind and clear sky favor weak turbulence and mixing, and strong temperature gradients near the surface. The boundary layer height varies depending on the relative contribution of the mechanical and thermal generation of turbulence. Because of the extremely low temperature and humidity, and the high elevation, Dome C is a potentially ideal site for astronomical observa- tions. For this reason, the optical turbulence over the Antarctic plateau has been a subject of studies by astronomers [Lawrence et al. 2004, Aristidi et al. 2005, Agabi et al. 2006, Lascaux et al. 2009]. At Dome C, Ricaud et al. [2012] made an ex- periment to monitor the vertical evolution of the planetary boundary layer (PBL) tempera- ture and humidity in the transition from win- ter to summer, by using a microwave radiome- ter operating at 60 GHz and 183 GHz. Due to the instrument low vertical resolution, the fine structure of the thermal turbulence could not be evidenced. * Corresponding author: Giampietro Casasanta, g.casasanta@isac.cnr.it A ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 2 Ghenton et al. [2010] analyzed the Dome C 45 m meteo tower measurements (tempera- ture, humidity, wind speed and direction measurements) for a three-week period in summer 2008. The main task of their work was to compare these measurements with the 6- hourly European Center for Medium-Range Forecasts (ECMWF) analyses and the daily ra- diosoundings. Pietroni et al. [2012], using the temperature profiles measured with a passive microwave radiometer [Kadygrov and Pick 1998, Argen- tini et al. 2004], characterized the behavior of the surface-based temperature inversions over the course of a year. They found that during the winter and the summer “nights” strong temperature inversions allow for a mixing depth of a few tens of meters with a quiescent layer above, decoupled from the surface layer. During the summer, despite the low surface temperatures, weak convection generates the development of a mixed layer characterized by a maximum depth of 200-400 m [Argentini et al. 2005]. The diurnal behavior of this mixed layer, monitored with a sodar, was described by Mastrantonio et al. [1999], Argentini et al. [2005], and King et al. [2006]. The sodar meas- urements, because of the membrane ringing just after the tone burst emission, allowed for a first echo recording starting at 20-30 m, de- pending on the membrane ringing time. Be- cause of this limitation, those measurements did not allow to study the surface turbulent layer under stable conditions, neither in sum- mer nor in winter. An advanced high-resolution sodar named surface-layer sodar (hereafter SL-sodar), allow- ing for the lowest observation height at ≈ 2 m and a vertical resolution of ≈ 2 m, was devel- oped by the ISAC-CNR [Argentini et al. 2011]. The SL-sodar was deployed at Concordia sta- tion after a preliminary test period at the ISAC-CNR research centre of Rome [Argentini et al. 2011]. In this paper, a few preliminary re- sults from the summer season are shown. II. SITE AND INSTRUMENTATION Concordia is a permanent station located at Dome C (75.1° S, 123.3° E, 3233 m a.s.l.), on the East Antarctic plateau, at approximately 1000 km from the nearest coast. One-year in situ turbulence and radiation measurements, as well as SL-sodar observations, were carried out at Concordia station from December 2011 up to December 2012. Table 1. SL-sodar setting parameters. Mode 1 Mode 2 Carrier Frequency 2000 Hz 4850 Hz Pulse duration 50 ms 10 ms Repetition rate 3 s 2 s Maximum range 430 m 280 m Lowest height 8 m 2 m Vertical resolution 8 m 2 m ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 3 The SL-sodar [Argentini et al. 2011] is an im- proved version of the sodar described by Mastrantonio et al. [1999] and Argentini and Pietroni [2010], with the possibility of zooming in the atmospheric surface-layer thermal tur- bulent structure. The SL-sodar consists of 3 horn-type antennas, placed symmetrically around a 1.2 m diameter parabolic receiving antenna, emitting simultaneously acoustic pulses at the same frequency. The receiving antenna is noise-protected by a shielding struc- ture of 1.5 m L x 1.5 m W x 2.0 m H in size. The transmitting and receiving circuits are kept separated to minimize the “cross-talk” be- tween channels. The carrier frequency, the pulse duration, and the pulse repetition rate Figure 1. Sodargram for 28 December 2011 with mixing height estimate (straight line) (a), downwelling longwave radiation (𝐿𝑊 ↓ ) and sonic temperature ( 𝑇𝑠) (b), wind speed and direction (c), momentum (𝑢∗) and heat fluxes (𝐻0) (d). ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 4 can change according to the two modes listed in Table 1. “Mode 1” allows to monitor the convective mixed layer, while “Mode 2”, with higher vertical resolution, is used to investi- gate the near-surface stable layer. Measurements of turbulence were made with a Metek USA-1, a three-axes sonic thermo- anemometer (sampling frequency of 10 Hz) installed on a 3.5 m mast. The heat and mo- mentum fluxes are estimated using the eddy covariance method [Lee et al. 2004]. The longwave and shortwave radiation compo- nents (up and down) are measured using Kipp & Zonen CNR1 pyrgeometers and pyranome- ters, installed at 1.5 m above the snow surface. In this paper, unless told otherwise, the local standard time (LST) is used. III. RESULTS Summertime ABL diurnal behavior The facsimile recording of the verti- cal/temporal variation of the acoustic backscattering (sodargram) “depicts” the thermal structure of the atmosphere. The scat- tering elements producing the change of echo intensity are the small-scale temperature in- homogeneities due to thermal turbulence. Temperature fluctuations are usually associat- ed with the convective plumes originating from the surface, or with potential temperature gradients and wind shear usually occurring in the inversion layers. During the summer, the boundary layer at Dome C can reach a depth of 200-400 m. Therefore, the “Mode 1” setting (see Table 1) was used to catch the whole ver- tical evolution during the daily cycle. Figure 1a shows the typical boundary layer sodargram during a clear summer day (28 De- cember 2011), with the superimposed estimate of the mixing height (MH, straight line). A sta- ble boundary layer occurs between 0000 and 0900 LST; from 0900 to 1630 LST the MH in- creases because of the convective activity, and then it drops to 50 m because of the surface ra- diative cooling. The MH was estimated using a method origi- nally proposed by Beyrich and Weill [1993], which uses the backscattered range corrected signal (RCS). Under stable nocturnal stratifica- tion, the MH was determined either from the minimum of the first derivative or from the maximum curvature of the RCS, depending on the stage of the planetary boundary layer evo- lution, and on the shape of the sodar profile. Under convective conditions, the MH was es- Figure 2. Sodargram for 5 February 2012, the full dots represent the mixing height estimate. High ampli- tude waves are observed between 0800–0930 LST. ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 5 timated as the height at which an elevated sec- ondary maximum occurs (i.e. the height of the turbulent zone characterizing the top of the mixing layer). For the same day, Figure 1b shows the downward longwave radiation 𝐿𝑊 ↓ and the sonic temperature 𝑇𝑠, Figure 1c the wind speed and direction. The heat turbu- lent flux 𝐻0 and friction velocity 𝑢∗ are plotted in Figure 1d. 𝐿𝑊 ↓ ranges between 100 and 150 W m-2, while 𝑇𝑠 reaches its minimum (- 35°C) at 0300 LST, and its maximum (-25°C) between 1200 and 1500 LST. The direction in- dicates a wind from the continent persisting the whole day. The maximum wind speed (6 ms-1) occurs because of the momentum transfer from the free atmosphere to the sur- face layer, as a consequence of the turbulent mixing (confirmed by the positive and increas- ing values for 𝐻0 and 𝑢∗) during convective hours. Gravity waves in the ABL Between 2 and 5 February 2012, waves with periods of a few minutes are observed under stable stratification for more than 35% of the time. The resolution achieved by the SL-sodar with the “Mode2” setting (Table 1) allowed to visualize the fine structure of these wave pat- terns. At the transition time from the stable to the unstable boundary layer (between 0800 and 1000 LST) the capping inversion layer os- cillates with an amplitude that reaches 70 m (Figure 2). The apparent period of these oscil- lating structures was estimated through the spectral analysis of the sonic temperature and the wind components of the sonic anemome- ters, installed at 3 different levels (7.0, 22.8, and 37.5 m) on a 45-m meteorological tower [Gen- thon et al. 2009] located at ≈ 1 km from the SL- sodar. In Figures 3a and 3b, the sodargram and the power spectral density of the temperature and horizontal wind components measured at 37.5 m (between 0945 and 1000 LST on 5 February) are shown. A peak occurs simultaneously in the temperature and wind components Power Spectral Density (PSD) at 0.0098 Hz, corre- sponding to a period of 102 s. The analysis of the sonic anemometers measurements at the three levels, limited to the time interval 0800- Figure 3. Sodar echogram for February 5, 2012 (0930 to 1000 LST) (a); Power Spectral Density (PSD) of the horizontal wind components (Ux, Uy) and sonic temperature (𝑇𝑠) at 39 m (b) for the indicated 10- minutes interval (the scale is given in arbitrary units). ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 6 0930 LST, gives an apparent period ranging between 90 and 120 s. The apparent period remains approximately the same also when the inversion strength and height change. This be- havior indicates that the origin of these waves might be a disturbance (probably the wind shear) originating between the inversion layer and the free atmosphere. A similar behavior was observed during other days. A summer warming event Warming events of particular intensity were regularly observed at Dome C during the win- ter [Argentini et al. 2001, Petenko et al. 2007, Ghenton et al. 2013]. During these events the surface temperature sometimes has a sharp in- crease of 20-40 °C [Argentini et al. 2001], reach- ing then the typical summer values. Studies carried out at South Pole [Carroll 1982, Stone et al. 1990, Stone and Kahl 1991, Stone 1993] have evidenced that these warming events are generally observed in presence of clouds. Neff [1999], analyzing the particles tra- jectories across Antarctica, found that these warming processes are mostly due to warm and moist air intrusion and to the condensa- tion of nuclei originating from the Weddell Sea, producing a wide variety of cloud types. Carroll [1982] suggested two possible mecha- nisms of this phenomenon: the advection of warm air, and/or the vertical mixing of air from different layers. Schwerdtfeger and Weller [1977] related the surface warming to the variation of long-wave radiation emitted by the clouds associated to the moist air in the upper part of the atmosphere. Figure 4. Downwelling longwave radiation (𝐿𝑊 ↓) and sonic temperature (𝑇𝑠) with superimposed the linear trends (a), the wind speed and direction (b), the heat (𝐻0) and momentum (𝑢∗) fluxes (c) for days 8- 14 January 2012. ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 7 The measurements collected during a summer warming event observed between 8 and 17 January 2012 have been analyzed. Figure 4 shows the time series of the downwelling longwave radiation 𝐿𝑊 ↓ and the sonic tem- perature 𝑇𝑠 (Figure 4a), the wind speed and di- rection (Figure 4b), the heat flux 𝐻0 and the friction velocity 𝑢∗ (Figure 4c) during the se- lected period. Starting from 9 January 2012, wet and warm coastal air masses are advected from the coast toward the Dome C area. Until the end of Jan- uary 10 the 𝐿𝑊 ↓ and 𝑇𝑠 (Figure 4a) do not show the typical diurnal behavior observed during the previous and the following days. An increasing trend is evident starting on 9 January. The wind direction changes, rotating from S to NE-NW, and the wind speed is low up to 1200 LST of 10 January. Due to the pres- ence of the clouds, the downwelling long wave radiation and the surface temperature increase, initiating the convective activity shown in the sodargram of Figure 5. The vertical mixing re- duces the decoupling between the boundary layer and the free atmosphere. As a conse- quence, the wind speed increases: at 2400 LST of 11 January the wind speed is ≈ 6 m s-1 (Fig- ure 4b). 𝐻0 and 𝑢∗ (Figure 4c) confirm this anomalous behavior. At the end of 11 January Figure 5. Sodargram for 10 January 2012 with the mixing height estimate (dots). ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 8 the behavior is again the typical summer one, with a peak in the wind speed at 1200 LST. The sodargram for 10 January, with the mixing height estimate superimposed (dots), confirms this hypothesis. A clear and intense convective activity is observed during the whole day even in the nighttime between 9 and 10 January. IV. SUMMARY The main results of this observational study can be summarized as follows: • during the summer, under steady weather conditions, the atmospheric boundary layer thermal structure is characterized by the alter- nation of a stable stratified layer with a con- vective boundary layer, following a behavior similar to that observed at mid-latitudes. • A regular wave activity was observed within the inversion layer. The time period of these waves ranges between 90 s and 120 s, and their origin can be attributed to the wind shear across the inversion layer. • The summer warming events take origin from the presence of clouds advected from the coast toward the Dome C area. Clouds modify the surface radiation budget by increasing the downward longwave radiation, which in turn produces an increase of surface temperature, leading to convection. Due to the decrease of the temperature inversion strength, the vertical mixing, combined with the wind shear, allows the transport of warm air from the upper parts of the atmosphere towards the surface, further contributing to the surface warming. Acknowledgements. This research is sup- ported by the Italian National Research Programme (PNRA) in the frame of French-Italian projects for Dome C. The authors thank Mr. A. Conidi, Mr. N. Ferrara, and Mr. S. Ciampichetti for the technical support provided in the realization of the SL-sodar antennas, and the logistics staff of the Concordia station during the field work REFERENCES [Agabi et al. 2009] Agabi, A., E. Aristidi, M. Azouit, E. Fossat, F. Martin, T. Sadibekova, J. Vernin, and A. Ziad (2006). First Whole At- mosphere Nighttime Seeing Measurements At Dome C, Antarctica, Publ. Astron. Soc. Pac., 118, 840, 344-348. [Argentini et al. 2001] Argentini, S., I. Petenko, G. Mastrantonio, V. Bezverkhnii and A. Viola (2001). Spectral characteristics of East Antarcti- ca meteorological parameters during 1994, J. Geophys. Res., 106, 12463-12476. [Argentini et al. 2004] Argentini S, A. Conidi, A. Viola, G. Mastrantonio, N. Ferrara, I. Peten- ko, E.N. Kadygrov, A.V. Koldaev and A. Via- zankin (2004). Temperature measurements at Dome C using a new microwave temperature profiler, X Workshop on Antarctic Atmos- phere, 22-24 October 2003. [Argentini et al. 2005] Argentini S., A. Viola, A. Sempreviva and I. Petenko (2005). Summer boundary-layer height at the plateau site of Dome C, Antarctica, Boundary-Layer Meteor- ol., 115, 409-422. [Argentini et al. 2011] Argentini, S., I. Pietroni, G. Mastrantonio, I. Petenko and A. Viola (2011). Use of a high-resolution sodar to study surface-layer turbulence at night. Boundary- Layer Meteorol., 143, 177-188. [Aristidi et al. 2005] Aristidi, E., A. Agabi, E. Fossat, M. Azouit, F. Martin, T. Sadibekova, T. Travouillon, J. Vernin, And A. Ziad (2005), Site Testing In Summer At Dome C, Antarctica, As- tron. Astrophys., 444, 2, 651-659. [Beyrich and Weill 1993] Beyrich, F. and A. Weill (1993). Some aspects of determining the ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 9 stable boundary layer depth from sodar data, Boundary-Layer Meteorol., 63, 97-116. [Carrol 1982] Carroll, J.J. (1982). Long-term means and short-term variability of the surface energy balance components at the South Pole, J. Geophys. Res., 87, 4277-4286. [Genthon et al. 2010] Genthon C. , D. Six, V. Favier, S. Argentini. A. Pellegrini (2010). Me- teorological Atmospheric Boundary Layer Measurements and ECMWF Analyses During Summer At Dome C, Antarctica, J. Geophys. Res. Atmos., 115, D05104. [Genthon et al. 2013] Genthon C., H. Gallee, D. Six, P. Grigioni, A. Pellegrini (2013). Two years of atmospheric boundary layer observation on a 45-m tower at Dome C on the Antarctic plat- eau, J. Geophys. Res. Atmos. 118, 3218-3232. [Kadygrov and Pick 1988] Kadygrov E.N., D.R. Pick (1998). The potential for temperature re- trieval from an angular-scanning single- channel microwave radiometer and some comparisons with in situ observations, Mete- orol. Appl., 5, 393-404. [King et al. 2006] King, J.C., S. Argentini and P.S. Anderson (2006). Contrasts between the summertime surface energy balance and boundary layer structure at Dome C and Hal- ley stations, Antarctica J. Geophys. Res., 111, D02105. [Lascaux et al 2009] Lascaux, F., E. Masciadri, S. Hagelin, J. Stoesz (2009). Mesoscale Optical Turbulence Simulations At Dome C, Mon. Not. R. Astron. Soc., 398(3), 1093-1104. [Lawrence et al. 2004] Lawrence J.S., M.C.B. Ashley, A. Tokovinin, T. Trouvillon (2004). Ex- ceptional Astronomical Seeing Conditions Above Dome C In Antarctica, Nature, 431, 278- 281. [Lee et al. 2004] Lee, X., W. Massman and B. Law (2004). Handbook of Micrometeorology, Kluwer Academic Publishers, Dordrecht. [Mastrantonio et al. 1999] Mastrantonio, G., V. Malvestuto, S. Argentini, T. Georgiadis and A. Viola (1999). Evidence of a convective bounda- ry layer developing on the Antarctic Plateau during the summer, Meteorol. Atmos. Phys., 71, 127-132. [Neff 1999] Neff, W.D. (1999). Decadal time scale trends and variability in the tropospheric circulation over the South Pole, J. Geophys. Res., 104, 217-251. [Petenko et al. 2007] Petenko I., S. Argentini, I. Pietroni, G. Mastrantonio, A. Viola (2007). Warming events during winter in Antarctica. In Calacino, M. and Rafanelli, C. (eds.), Proc. 11th Workshop Italian Research on Antarctic Atmosphere, Rome, 10-12 April 2007, p.119- 132. [Pietroni et al. 2012] Pietroni, I., S. Argentini, I. Petenko and R. Sozzi (2012). Measurements and parameterizations of the atmospheric boundary layer height at Dome C, Antarctica, Boundary-Layer Meteorol., 143, 189-206. [Ricaud et al. 2012] Ricaud P., C. Ghenton, P. Durand, J. L. Attié, F. Carminati, G. Canut, J. F. Vanacker, L. Moggio, Y. Courcoux, A. pelle- grini, T. Rose (2012). Summer to winter diurnal variabilities of temperature and water vapour in the lowermost troposphere as observed by HAMSTRAD over Dome C, Antarctica. Boundary-Layer Meteorol., 143, 227-259. [Schwerdtfeger and Weller 1977] Schwerdtfeger, P., G. Weller (1977). Radiative heat transfer processes in snow and ice, in Businger , J.A. (Ed.), Meteorological Studies at Plateau Station, Antarctic Research Series, 25, Amer. Geophys. Union, 35-39. [Stone et al. 1990] Stone, R.S., G.E. Dutton and J.J. DeLuisi (1990). Surface radiation and tem- perature variations associated with Cloudiness at the South Pole, Antarct. J. U. S, 24, 230-232. ANNALS OF GEOPHYSICS, 56, 5, 2013; 10.4401/ag-6347 10 [Stone and Kahl 1991] Stone, R.S., and J.D. Kahl (1991). Variations in Boundary layer Properties associated with Clouds and Transi- ent weather disturbance at the South Pole dur- ing winter, J. Geophys. Res., 96, 5137-5144. [Stone 1993] Stone R.S. (1993). Properties of austral winter clouds derived from radiometric profiles at the South Pole, J. Geophys. Res., 98, 12961-1297. The surface layer observed by a high-resolution sodar at DOME C, Antarctica Stefania Argentini1, Igor Petenko1, Angelo Viola1, Giangiuseppe Mastrantonio1, Ilaria Pietroni1, Giampietro Casasanta1,*, Eric Aristidi2, Christophe Ghenton3 1 Istituto di Scienze dell'Atmosfera e del Clima, CNR, Rome, Italy 2Laboratoire Lagrange, UNS/CNRS/OCA, Nice Cedex 2, France 3 Laboratoire de Glaciologie et Géophysique de l’Environnement, UJF-Grenoble 1/CNRS, Grenoble, France Abstract Table 1. SL-sodar setting parameters.