Parameterizing turbulent exchange over sea ice: the ice station weddell results

A 4-month deployment on Ice Station Weddell (ISW) in the western Weddell Sea yielded over 2000 h of nearly continuous surface-level meteorological data, including eddy-covariance measurements of the turbulent surface fluxes of momentum, and sensible and latent heat. Those data lead to a new parameterization for the roughness length for wind speed, z₀, for snow-covered sea ice that combines three regimes: an aerodynamically smooth regime, a high-wind saltation regime, and an intermediate regime between these two extremes where the macroscale or `permanent' roughness of the snow and ice determines z₀. Roughness lengths for temperature, z_T, computed from this data set corroborate the theoretical model that Andreas published in 1987. Roughness lengths for humidity,z_Q, do not support this model as conclusively but are all, on average, within an order of magnitude of its predictions. Only rarely arez_Tand z_Q equal to z₀. These parameterizations have implications for models that treat the atmosphere-ice-ocean system.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

Characteristics of Boundary Layer Structure during a Persistent Haze Event in the Central Liaoning City Cluster,Northeast China

The characteristics of boundary layer structure during a persistent regional haze event over the central Liaoning city cluster of Northeast China from 16 to 21 December 2016 were investigated based on the measurements of particulate matter (PM) concentration and the meteorological data within the atmospheric boundary layer (ABL). During the observational period, the maximum hourly mean PM2.5 and PM10 concentrations in Shenyang, Anshan, Fushun, and Benxi ranged from 276 to 355 μg m^–3 and from 378 to 442 μg m^–3, respectively, and the lowest hourly mean atmospheric visibility (VIS) in different cities ranged from 0.14 to 0.64 km. The central Liaoning city cluster was located in the front of a slowly moving high pressure and was mainly controlled by southerly winds. Wind speed (WS) within the ABL (< 2 km) decreased significantly and WS at 10-m height mostly remained below 2 m s^–1 during the hazy episodes, which was favorable for the accumulation of air pollutants. A potential temperature inversion layer existed throughout the entire ABL during the earlier hazy episode [from 0500 Local Time (LT) 18 December to 1100 LT 19 December], and then a potential temperature inversion layer developed with the bottom gradually decreased from 900 m to 300 m. Such a stable atmospheric stratification further weakened pollutant dispersion. The atmospheric boundary layer height (ABLH) estimated based on potential temperature profiles was mostly lower than 400 m and varied oppositely with PM_2.5 in Shenyang. In summary, weak winds due to calm synoptic conditions, strong thermal inversion layer, and shallow atmospheric boundary layer contributed to the formation and development of this haze event. The backward trajectory analysis revealed the sources of air masses and explained the different characteristics of the haze episodes in the four cities.     

A Large-Eddy Simulation Study of Water Vapour and Carbon Dioxide Isotopes in the Atmospheric Boundary Layer

A large-eddy simulation model developed at the National Center for Atmospheric Research (NCAR) is extended to simulate the transport and diffusion of C¹⁸OO, H₂¹⁸O and ¹³CO₂ in the atmospheric boundary layer (ABL). The simulation results show that the ¹⁸O compositions of leaf water and the ABL CO₂ are moderately sensitive to wind speed. The variations in the ¹⁸O composition of water vapour are an order of magnitude greater than those in the ¹³C and ¹⁸O compositions of CO₂ both at turbulent eddy scales and across the capping inversion. In a fully-developed convective ABL, these isotopic compositions are well mixed as with other conserved atmospheric quantities. The Keeling intercepts determined with the simulated high-frequency turbulence time series do not give a reliable estimate of the ¹⁸O composition of the surface water vapour flux and may be a reasonable approximation to the ¹³C and ¹⁸O compositions of the surface CO₂ flux in the late afternoon only after a deep convective ABL has developed. We suggest that our isotopic large-eddy simulation (ISOLES) model should be a useful tool for testing and formulating research hypotheses on land–air isotopic exchanges.     

Airborne surveys of the atmospheric boundary layer above the marginal ice zone on the Newfoundland shelf

Abstract

Airborne measurements in the atmospheric boundary layer (ABL) above the marginal ice zone (MIZ) on the Newfoundland Shelf reveal strong lateral variations in mean wind, temperature and the vertical fluxes of heat and momentum under conditions of cold, off‐ice wind. Flux measurements in (and near) the surface layer indicate that the neutral 10‐m drag coefficient depends on ice concentration, ranging from 2 × 10^‐3 at 10% coverage to 5 × 10^‐3 at 90%. Furthermore, cross‐ice‐edge transects consistently show increasing wind speed, temperature and heat flux in the off‐ice direction, but the momentum flux may either increase or decrease, depending on the relative importance of surface buoyancy flux and roughness. For the conditions encountered in this experiment, it appears surface wave maturity does not have a significant influence on the drag coefficient in fetch‐limited regimes near the ice edge.     

‘A Case Study on a Polar Cold Air Outbreak over Fram Strait using a Mesoscale Weather Prediction Model’

The mesoscale weather prediction model ’Lokal-Modell’ (LM) of the Deutscher Wetterdienst is applied to the situation of an Arctic cold air outbreak in the Fram Strait region in April 1998. Observations are available from a flight along 5⁰E carried out during the ARTIST campaign. Initial and time-dependent boundary data for the simulation are taken from a larger scale operational model system. Using the standard configuration of LM, the simulation reproduced the propagation of cold air and the characteristic structure of the atmospheric boundary layer (ABL) in fair agreement with the observations. However, a detailed comparison revealed three basic problems. Firstly, there is evidence that the available data on sea-ice conditions were insufficient approximations to the true state for several reasons. A modification of the sea-ice data towards observations revealed that parts of the discrepancies were due to the original sea-ice data. Secondly, a control run with the model in its standard configuration shows an insufficient warming of the ABL downstream of the ice edge due to underestimation of surface heat fluxes. A simple modification of the approach for the scalar roughness length resulted in the strongest benefit, while comparative studies showed only a slight sensitivity to different types of parametrisation of turbulent mixing or the inclusion of an additional moist convection parametrisation. Thirdly, in all the simulations the deepening of the convective ABL downstream of the ice edge is weaker than observed. This may be partly due to the thermal stratification above the ABL in the analysis data, which is more stable than observed; but it may also be a hint to the fact that processes near the inversion are insufficiently parametrised in mesoscale models with resolutions as used in LM. The simulated cloud layer in the convective ABL is similar to that observed with respect to condensate content, a sharply defined cloud top, a diffuse lower bound, and continuous light precipitation.     

Boundary-layer observations over water and arctic sea-ice during on-ice air flow

Observations made on 8 and 9 May 1988 by aircraft and two ships in and around the marginal ice zone of the Fram Strait during on-ice air flow under cloudy and cloud-free conditions are presented.The thermodynamic modification of the air mass moving from the open water to the ice over horizontal distances of 100–300 km is only a few tenth of a degree for temperature and a few tenth of a gram per kilogram for specific humidity. This is due to the small temperature differences between sea and ice surfaces. During the day, the ice surface is even warmer than the sea surface. The stably stratified 200–400 m deep boundary layer is often topped by a moisture inversion leading to downward fluxes of sensible as well as latent heat.The radiation and energy balance at the surface are measured as functions of ice cover, cloud cover and sun elevation angle. The net radiationR _Nis the dominating term of the energy budget. During the day, the difference ofR _Nbetween clear and overcast sky is only a few W/m² over ice, but 100–200 W/m² over water. During the night,R _Nover ice is more sensitive to cloud cover.The kinematic structure is characterized by strong shears of the longitudinal and the transversal wind component. The profile of the latter one shows an inflection point near the top of the boundary layer. Dynamically-driven roll circulations are numerically separated from the mean flow. The secondary flow patterns have wavelengths of about 1 km and contribute substantially to the total variances and covariances.     

The effect of snow on Antarctic sea ice simulations in a coupled atmosphere-sea ice model

 The effect of a snow cover on sea ice accretion and ablation is estimated based on the ‘zero-layer’ version sea ice model of Semtner, and is examined using a coupled atmosphere-sea ice model including feedbacks and ice dynamics effects. When snow is disregarded in the coupled model the averaged Antarctic sea ice becomes thicker. When only half of the snowfall predicted by the atmospheric model is allowed to land on the ice surface sea ice gets thicker in most of the Weddell and Ross Seas but thinner in East Antarctic in winter, with the average slightly thicker. When twice as much snowfall as predicted by the atmospheric model is assumed to land on the ice surface sea ice also gets much thicker due to the large increase of snow-ice formation. These results indicate the importance of the correct simulation of the snow cover over sea ice and snow-ice formation in the Antarctic. Our results also illustrate the complex feedback effects of the snow cover in global climate models. In this study we have also tested the use of a mean value of 0.16 Wm^-1 K^-1 instead of 0.31 for the thermal conductivity of snow in the coupled model, based on the most recent observations in the eastern Antarctic and Bellingshausen and Amundsen Seas, and have found that the sea ice distribution changes greatly, with the ice becoming much thinner by about 0.2 m in the Antarctic and about 0.4 m in the Arctic on average. This implies that the magnitude of the thermal conductivity of snow is of considerable importance for the simulation of the sea ice distribution. An appropriate value of the thermal conductivity of snow is as crucial as the depth of the snow layer and the snowfall rate in a sea ice model. The coupled climate models require accurate values of the effective thermal conductivity of snow from observations for validating the simulated sea ice distribution under the present climate conditions. Received: 20 November 1997/Accepted: 27 July 1998     

北极低空急流和低层逆温特征观测分析

利用北冰洋冰表面热量平衡计划1997年10月中旬至1998年10月上旬的探空气球探测结果,分析了北极地区近地层逆温和低空急流特征.结果表明,96％的观测时次(11:15和23:15,协调世界时)出现近地层逆温,其中22％的逆温为贴地逆温,70％的逆温厚度在250～850 m之间,冬半年贴地逆温发生频率、逆温层厚度和逆温层内的温度变化都明显要大于夏半年.全年间低空急流出现频率为41％,平均高度为520 m,最大频率出现在150 m附近,70％的急流出现在600m高度以下.急流平均风速为10.6m·s-1,风速在4～13 m·s-1范围内的急流约占总数的75％,东和东南方向为全年急流的主导风向.根据对急流核和地面风速之间转换角分布的分析,惯性震荡可能是北极低空急流的主要成因.     

General Characteristics of the Atmospheric Boundary Layer Over a Flaw Lead Polynya Region in Winter and Spring

A time series of microwave radiometric profiles over Arctic Canada’s Cape Bathurst (70°N, 124.5°W) flaw lead polynya region from 1 January to 30 June, 2008 was examined to determine the general characteristics of the atmospheric boundary layer in winter and spring. A surface based or elevated inversion was present on 97% of winter (January–March) days, and on 77% of spring (April–June) days. The inversion was the deepest in the first week of March (≈1100 m), and the shallowest in June (≈250 m). The mean temperature and absolute humidity from the surface to the top of the inversion averaged 250.1 K (−23.1°C), and 0.56 × 10⁻³ kg m⁻³ in winter, and in spring averaged 267.5 K (−5.6°C), and 2.77 × 10⁻³ kg m⁻³. The median winter atmospheric boundary-layer (ABL) potential temperature profile provided evidence of a shallow, weakly stable internal boundary layer (surface to 350 m) topped by an inversion (350–1,000 m). The median spring profile showed a shallow, near-neutral internal boundary layer (surface to 350 m) under an elevated inversion (600–800 m). The median ABL absolute humidity profiles were weakly positive in winter and negative in spring. Estimates of the convergence of sensible heat and water vapour from the surface that could have produced the turbulent internal boundary layers of the median profiles were 0.67 MJ m⁻² and 13.1 × 10⁻³ kg m⁻² for the winter season, and 0.66 MJ m⁻² and 33.4 × 10⁻³ kg m⁻² for the spring season. With fetches of 10–100 km, these accumulations may have resulted from a surface sensible heat flux of 15–185 W m⁻², plus a surface moisture flux of 0.001–0.013 mm h⁻¹ (or a latent heat flux of 0.7–8.8 W m⁻²) in winter, and 0.003–0.033 mm h⁻¹ (or a latent heat flux of 2–22 W m⁻²) in spring.     

Generation and atmospheric heat exchange of coastal polynyas in the Weddell Sea

The forcing mechanisms for Antarctic coastal polynyas and the thermodynamic effects of existing polynyas are studied by means of an air-sea-ice interaction experiment in the Weddell Sea in October and November 1986.Coastal polynyas develop in close relationship to the ice motion and form most rapidly with offshore ice motion. Narrow polynyas occur frequently on the lee side of headlands and with strong curvature of the coastline. From the momentum balance of drifting sea ice, a forcing diagram is constructed, which relates ice motion to the surface-layer wind vector v _z and to the geostrophic ocean current vector c _g . In agreement with the data, wind forcing dominates when the wind speed at a height of 3 m exceeds the geostrophic current velocity by a factor of at least 33. This condition within the ocean regime of the Antarctic coastal current usually is fulfilled for wind speeds above 5 m/s at a height of 3 m.Based on a nonlinear parameter estimation technique, optimum parameters for free ice drift are calculated. Including a drift dependent geostrophic current in the ice/water drag yields a maximum of explained variance (91%) of ice velocity.The turbulent heat exchange between sea ice and polynya surfaces is derived from surface-layer wind and temperature data, from temperature changes of the air mass along its trajectory and from an application of the resistance laws for the atmospheric PBL. The turbulent heat flux averaged over all randomly distributed observations in coastal polynyas is 143 W/m². This value is significantly different over pack ice and shelf ice surfaces, where downward fluxes prevail. The large variances of turbulent fluxes can be explained by variable wind speeds and air temperatures. The heat fluxes are also affected by cloud feedback processes and vary in time due to the formation of new ice at the polynya surface.Maximum turbulent fluxes of more than 400 W/m² result from strong winds and low air temperatures. The heat exchange is similarly intense in a narrow zone close to the ice front, when under weak wind conditions, a local circulation develops and cold air associated with strong surface inversions over the shelf ice is heated above the open water.     

Analysis of Atmospheric Boundary Layer Height Characteristics over the Arctic Ocean Using the Aircraft and GPS Soundings

Utilizing aircraft sounding data collected from the Surface Heat Budget of the Arctic Ocean (SHEBA, 1998) campaign, the authors evaluated commonly-used profile methods for Arctic ABL height estimation by validating against the’true’ABL height from aircraft sounding profiles, where ABL height is defined as the top of the layer with significant turbulence. Furthermore, the best performing method was used to estimate ABL height from the one-year GPS soundings obtained during SHEBA (October 1997-October 1998). It was found that the temperature gradient method produces a reliable estimate of ABL height. Additionally, the authors determined optimal threshold values of temperature gradient for stable boundary layer (SBL) and convective boundary layer (CBL) to be 6.5 K/100 m and 1.0 K/100 m, respectively. The maximum ABL height during the year was 1150 m occurred in May. Median values of Arctic ABL height in May, June, July, and August were 400 m, 430 m, 180 m, and 320 m, respectively. Arctic ABL heights are clearly higher in the spring than in the summer.     

Fuzzy logic based melting layer recognition from 3 GHz dual polarization radar: appraisal with NWP model and radio sounding observations

The advent of polarimetry makes it possible to categorize hydrometeor inferences more accurately by providing detailed information of the scattering properties. In light of this, the authors have developed a fuzzy logic based system for the recognition of melting layer in the atmosphere. The fuzzy system is based on characterizing melting layer scatterers from non-melting scatterers using five crisp inputs, namely, horizontal reflectivity (Z _H), differential reflectivity (Z _DR), co-polar correlation coefficient (ρ _HV), linear depolarization ratio (LDR) and height of radar measurements (H). For the implementation of melting layer recognition, the study employs the dual polarized signatures from the 3 GHz Chilbolton Advanced Meteorological Radar (CAMRA). Furthermore, a simple but effective averaging procedure for melting level estimation from a volume RHI scan is proposed. The proposed scheme has been evaluated with Weather Research and Forecasting (WRF) model simulated and radio soundings retrieved melting level height over a total of 84 RHI scan-based bright band cases. The results confirm that the estimated melting level heights from the proposed method are in good agreement with the WRF model and radio sounding observations. The 3 GHz radar melting level height estimates correspond with the R ² and RMSE values of 0.92 and 0.24 km, respectively, when compared to the radio soundings, and 0.93 and 0.21 km, respectively, when compared to the WRF model results. Moreover, the related R ² and RMSE values are reported as 0.93 and 0.22 km respectively between the WRF and radio soundings retrievals. This implies that the downscaled WRF modelled melting level height may also be used for operational or research needs.     

An Unprecedented Record Low Antarctic Sea-ice Extent during Austral Summer 2022

Seasonal minimum Antarctic sea ice extent (SIE) in 2022 hit a new record low since recordkeeping began in 1978 of 1.9 million km² on 25 February, 0.17 million km² lower than the previous record low set in 2017. Significant negative anomalies in the Bellingshausen/Amundsen Seas, the Weddell Sea, and the western Indian Ocean sector led to the new record minimum. The sea ice budget analysis presented here shows that thermodynamic processes dominate sea ice loss in summer through enhanced poleward heat transport and albedo–temperature feedback. In spring, both dynamic and thermodynamic processes contribute to negative sea ice anomalies. Specifically, dynamic ice loss dominates in the Amundsen Sea as evidenced by sea ice thickness (SIT) change, while positive surface heat fluxes contribute most to sea ice melt in the Weddell Sea.     

A Large-Eddy Simulation Study of the Influence of Subsidence on the Stably Stratified Atmospheric Boundary Layer

The influence of the large-scale subsidence rate, S, on the stably stratified atmospheric boundary layer (ABL) over the Arctic Ocean snow/ice pack during clear-sky, winter conditions is investigated using a large-eddy simulation model. Simulations of two 24-h periods are conducted while varying S between 0, 0.001 and 0.002 ms⁻¹, and the resulting quasi-equilibrium ABL structures and evolutions are examined. Simulations conducted with S = 0 yield a boundary layer that is deeper, more strongly mixed and cools more rapidly than the observations. Simulations conducted with S > 0 yield improved agreement with the observations in the ABL height, potential temperature gradients and bulk heating rates. We also demonstrate that S > 0 limits the continuous growth of the ABL observed during quasi-steady conditions, leading to the formation of a nearly steady ABL of approximately uniform depth and temperature. Subsidence reduces the magnitudes of the stresses, as well as the implied eddy-diffusivity coefficients for momentum and heat, while increasing the vertical heat fluxes considerably. Subsidence is also observed to increases the Richardson number to values in excess of unity well below the ABL top.     

Climate fluctuations of the Weddell Sea and its surroundings in a transient climate change scenario

The response of the Weddell Sea and Antarctic Peninsula to anthropogenic forcing simulated by a global climate model is analyzed. The model, despite its low resolution, is able to capture several aspects of the observed regional pattern of climate change. A strong warming and depletion of the sea ice cover in the western Weddell Sea contrasts with a slight cooling and a sea-ice extension in the eastern Weddell Sea. This simulated long-term climate change is modulated by interdecadal variability but also affected by some abrupt regional changes in the oceanic mixed layer depth. Between 1960 and 2030 a reorganization of the deep convection activity in the Weddell Sea sustains the opposition between the eastern and western Weddell Sea. The deep convection collapses in the western Weddell Sea in the 2030s. The sea ice retreat trend is then followed by an increase of the sea ice cover in the western Weddell Sea. In the eastern Weddell Sea another abrupt collapse of the deep convection activity occurs around 2080. This event is followed by a rapid cooling and sea ice extension during the next 20 years. Most of the surface changes are associated with large-scale atmospheric circulation changes that project on the dominant mode of natural variability but also with oceanic convection and circulation changes.     

A Simple Parametrization for the Concentration Variance Dissipation in a Lagrangian Single-Particle Model

A new quasi-analytical mixed-layer model is formulated describing the evolution of the convective atmospheric boundary layer (ABL) during cold-air outbreaks (CAO) over polar oceans downstream of the marginal sea-ice zones. The new model is superior to previous ones since it predicts not only temperature and mixed-layer height but also the height-averaged horizontal wind components. Results of the mixed-layer model are compared with dropsonde and aircraft observations carried out during several CAOs over the Fram Strait and also with results of a 3D non-hydrostatic (NH3D) model. It is shown that the mixed-layer model reproduces well the observed ABL height, temperature, low-level baroclinicity and its influence on the ABL wind speed. The mixed-layer model underestimates the observed ABL temperature only by about 10 %, most likely due to the neglect of condensation and subsidence. The comparison of the mixed-layer and NH3D model results shows good agreement with respect to wind speed including the formation of wind-speed maxima close to the ice edge. It is concluded that baroclinicity within the ABL governs the structure of the wind field while the baroclinicity above the ABL is important in reproducing the wind speed. It is shown that the baroclinicity in the ABL is strongest close to the ice edge and slowly decays further downwind. Analytical solutions demonstrate that the \(\mathrm{e}\)-folding distance of this decay is the same as for the decay of the difference between the surface temperature of open water and of the mixed-layer temperature. This distance characterizing cold-air mass transformation ranges from 450 to 850 km for high-latitude CAOs.     

The Parameterisation of Scalar Transfer over Rough Ice

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

Temperature stratification of the atmospheric boundary layer over the Deccan Plateau,India, during the summer monsoon

The vertical and horizontal temperature structure of the atmospheric boundary layer (ABL) were studied using aircraft observations made in the lowest 2.4 km above ground level during the summer monsoon.The vertical temperature structure of the ABL in the region may be classified into the following four categories.Category The ABL consisted of two layers of thickness 700–900 m separated by a thin transition layer. The lapse rates in the former two layers were dry adiabatic.Category The lowest layer of the ABL of thickness 400–600 m was adiabatically stratified and the overlying layer was stable with gradients of potential temperature 4–5°C km^–1. The stable layer contained a thin adiabatic stratified layer of 200–300 m thickness at a height of 1.5 km.Category The lowest 200–400 m layer of the ABL was adiabatically stratified and the overlying layer was stable with potential temperature gradients of 5–6 °C km¹.Category The ABL was mainly stable with potential temperature gradients of 6 °C km^–1 or greater. Occasionally thin layers with adiabatic stratification were found embedded in the ABL.The temperature distribution of the horizontal temperature at 900 m was mainly normal. The high-frequency portion of the spectra lying between 0.05 and 0.16 Hz (corresponding to wave length 1 km to 300 m) oscillated around the –\2/3 power law line. The spectral curve showed a significant peak at 0.011 Hz having a wave-length of 5 km.Department of Geoscience, North Carolina State University, Raleigh, NC, 27650, U.S.A.