Turbulence Characteristics of the Stable Boundary Layer Over a Mid-Latitude Glacier. Part I: A Combination of Katabatic and Large-Scale Forcing

Profile and eddy-correlation (heights of 4 and 10 m) measurements performed on the Pasterze glacier (Austria) are used to study the characteristics of the stable boundary layer under conditions of katabatic and large-scale forcing. We consider cases where large-scale forcing results in a downslope (or following) ambient wind. The analysis of averaged spectra and cospectra reveals low frequency perturbations that have a large influence on the variances of temperature and horizontal wind components and also alter the cospectra of momentum and sensible heat flux. Only the spectrum of the vertical wind speed is comparable to universal spectra. The low frequency perturbations occur as brief intermittent events and result in downward entrainment of ambient air thereby producing enhanced downward sensible heat fluxes and downward as well as upward momentum fluxes with various magnitudes and timescales. After the variances were high pass filtered, the normalised standard deviations of wind speed and temperature compare favourably to findings in the literature within the range 0>z/L>0.5. For larger z/L they deviate as a result of an increased influence from low frequency perturbations and thus non-stationarity. In line with this, the turbulent kinetic energy budget (at 4 m height) indicates that production (shear) is in balance with destruction (buoyancy and dissipation) within the range 0>z/L>0.3. Non-dimensional gradients of wind speed within the range 0>z/L>0.3 have a slope of about 3.5. The scatter for the dimensionless temperature gradient is quite large, and the slope is comparable to that for wind speed gradients. For z/L>0.3 the imbalance in the turbulent kinetic energy budget grows and non-dimensional gradients for wind speed and temperature deviate considerably from accepted values as a result of increased non-stationarity. Average roughness lengths for momentum and sensible heat flux derived from wind speed and temperature profiles are respectively 1 × 10^-3 m and 6 × 10^-5 m, consistent with the literature. The ratio (z_0h/z_0m) compares to those predicted by surface renewal models. A variation of this ratio with the roughness Reynolds number is not indicated by our data.     

Influence of the boundary layer height on the global air–sea surface fluxes

Results from large-eddy simulations and field measurements have previously shown that the velocity field is influenced by the boundary layer height, z _i , during close to neutral, slightly unstable, atmospheric stratification. During such conditions the non-dimensional wind profile, φ _m , has been found to be a function of both z/L and z _i /L. At constant z/L, φ _m decreases with decreasing boundary layer height. Since φ _m is directly related to the parameterizations of the air–sea surface fluxes, these results will have an influence when calculating the surface fluxes in weather and climate models. The global impact of this was estimated using re-analysis data from 1979 to 2001 and bulk parameterizations. The results show that the sum of the global latent and sensible mean heat fluxes increase by 0.77 W m⁻² or about 1% and the mean surface stress increase by 1.4 mN m⁻² or 1.8% when including the effects of the boundary layer height in the parameterizations. However, some regions show a larger response. The greatest impact is found over the tropical oceans between 30°S and 30°N. In this region the boundary layer height influences the non-dimensional wind profile during extended periods of time. In the mid Indian Ocean this results in an increase of the mean annual heat fluxes by 2.0 W m⁻² and an increase of the mean annual surface stress by 2.6 mN m⁻².     

An Improvement of the Louis Scheme for the Surface Layer in an Atmospheric Modelling System

The Louis scheme for calculating the vertical eddy fluxes within the atmospheric surface layer is improved by broadening the original assumptions. In our approach, the momentum and heat transfer roughness lengths (z₀ and z_T respectively) can be different, and z₀ need not be negligibly small compared with the lowest height (z) in modelling. For these conditions, we choose more consistent wind and potential temperature profile forms, then derive new algorithms for calculating fluxes. Improvement is demonstrated for a wide range of z/L (L is the Obukhov length), z/z₀ and z₀ z_T, by comparing these fluxes with those derived from a theoretical surface-layer model. The improved algorithms can be used in atmospheric modelling systems for more varied surfaces and a wide range of atmospheric stability.     

Wind profiles at the Boulder Tower

Analysis of wind profiles at the Boulder Tower (BAO) leads to these conclusions:

1. The variation of roughness with wind direction found earlier is confirmed. Roughness lengths measured on the tower are larger than those measured close to the surface.
2. The profiles and measurements of Reynolds stress are consistent with a von-Karman constant of 0.35.
3. The form φ_m=(1?15z/L)^-1/3 fits best in the range -0.6 < z/L < 0. In the range 0 < z/L < 0.5, θ _m ～ 1 + 4.7z/L provides a good fit to the observations. For z/L < 0.1, φ _m also depends on h, the thickness of the PBL. For z/L < -0.6, Φ _m approaches the constant 0.5, in contrast to all previous suggestions. For larger stabilities, the upper level is usually not in the surface layer, and wind ratios become independent of z/L.
4. With snow cover, the effective roughness diminishes to about 1 cm, even for directions for which the roughness length without snow is large.
5. Estimation of winds at 100 or 150 m from information near the surface is best for similarity theory provided that the ratio of height to Monin-Obukhov L is less than 0.1. For larger z/L, simple power laws seem more appropriate.

     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

Structure of the marine atmospheric boundary layer during two cold air outbreaks of varying intensities: GALE 86

Aircraft, surface, upper air and satellite measurements have been used to observe the evolution and growth of the convective Marine Atmospheric Boundary Layer (MABL) offshore of North Carolina in close proximity to the Gulf Stream, during the intense cold air outbreak of 28 January 1986 and the moderate event of 12 February 1986, as part of the Genesis of Atlantic Lows Experiment (GALE). Air mass modification processes, driven primarily by the ocean-atmosphere exchanges of surface turbulent sensible and latent heat fluxes, caused the overlying air mass to warm and moisten as it advected over the warmer waters of the eastern United States continental shelf. Maximum observed near-surface total heat fluxes were 1045 and 811 W·m^–2 over the core of the Gulf Stream, for 28 January and 12 February 1986, respectively. The observed changes in the overlying air mass occurred almost instantaneously as the ambient flow traversed different underlying SST conditions.The turbulent structure showed a buoyancy-dominated MABL below approximately 0.8z/h. However, shear was also observed to be an important production term above 0.8z/h and below 0.1z/h for the 28 January 1986 event. Dissipation of turbulent kinetic energy was the dominant destruction term in the budgets, but vertical transport of energy was a strong contributor below 0.5z/h, above which this term became a source of turbulent energy. Additionally, the normalized standard deviations of the horizontal velocity components showed a near-equal contribution to the turbulence, while the vertical velocity components displayed the characteristic mid-layer maximum profile observed for a convective, well-mixed boundary layer.     

Meteorological results of MONSOON-88 expedition (pre-monsoon period)

Mean atmospheric circulation, moisture budget and net heat exchange were studied during a pre-monsoon period (18th March to 3rd May, 1988), making use of the data collected on board Akademik Korolev in the central equatorial and southern Arabian Sea region. The net heat exchange (R _n ) is found to be about 20 W m^–2 for a small area (0–4° N; 55–60° E), 50% less than the dimatological value. The mean value of net radiation (140 W m^–2) is less than the climatological value, which was due to higher cloud amount. The higher SST enhanced both the latent and sensible heat fluxes.The mean atmospheric circulation obtained from the upper air data is quite convincing. The mean exchange coefficient (C _e ) estimated from the moisture budget is about 1.0 × 10^–3 for a wind speed of 4 m s^–1. This value is slightly lower than that obtained by the usual methods.National Institute of Oceanography, RC, 52-Kirlampudi layout, Visakhapatnam — 530 023.India Meteorological Department, Gauhati.     

Transfer coefficients of sensible heat on a snowmelt surface

Summary Field observations were carried out in order to determine the transfer coefficeient of sensible heat flux above a melting snow surface at the Moshiri experimental site. The coefficient is calculated as the ratio of the sensible heat flux determined by the eddy correlation method using a sonic anemometer to the product of the wind speed and the temperature difference between the air and snow surface. The sensible heat fluxes are also compared with the result of the precise heat balance observations. The nondimensional transfer coefficienth shows a good correlation with the atmospheric stability,Ri. The value ofh=2.3×10^–3 is obtained in the range of 0<Ri<0.1, however, it is smaller and scattered under stronger stable condition (Ri>0.1). The dimensional transfer coefficent for sensible heat flux is calculated, and a linear relationship is obtained as a function of the logarithm of atmospheric stability.With 11 Figures     

Coherent Structures near the Surface in a Strongly Sheared Convective Boundary Layer Generated by Large-Eddy Simulation

Well-developed low speed and high temperature streaks in association with the alignment of convection cells are observed in a large-eddy-simulation (LES) generated strongly sheared convective boundary-layer flow, which is driven by a geostrophic wind speed of 15 m s^-1 and a surface kinematic heat flux of 0.05 K m s^-1. Vortices that drive streaky structures are identified through an eigenvalue method (lambda;₂method) near the surface. These vortices are highly elongated along the quasi-streamwise direction alternating sign of the x-component of vorticity (_x). By conditional sampling of fully developed vortices, a statistically significant coherent structure is educed. The educed vortex is elongated to the streamwise direction with the elevation angle of about 17° above the horizontal surface. However, the horizontal tilting is not clearly demonstrated in the present simulation. Fluctuation fields in the domain of the educed vortex show the existence of a low speed and high temperature streak as a direct consequence of momentum and heat transport by vortical motions. The strong ejection(upward transport of low momentum or high temperature)occurring at the higher level than that of the strong sweep (downward transport of high momentum and low temperature) can be explained from the spatial distribution of the fluctuationfields of velocity and temperature. The contribution of ejection to the Reynolds stress at z/h₁ = 0.18 is about 75%, which is slightly greater than that (70% at z/h₁ = 0.173) for the neutrally stratified atmospheric boundary layer. Ejection is also found to be dominant for the turbulent heat flux.     

On the response of the oceanic wind-driven circulation to atmospheric CO2 increase

A global, flux-corrected climate model is employed to predict the surface wind stress and associated wind-driven oceanic circulation for climate states corresponding to a doubling and quadrupling of the atmospheric CO₂ concentration in a simple 1% per year CO₂ increase scenario. The model indicates that in response to CO₂ increase, the position of zero wind stress curl in the mid-latitudes of the Southern Hemisphere shifts poleward. In addition, the wind stress intensifies significantly in the mid-latitudes of the Southern Hemisphere. As a result, the rate of water circulation in the subpolar meridional overturning cell in the Southern Ocean increases by about 6 Sv (1 Sv=10⁶ m³ s⁻¹) for doubled CO₂ and by 12 Sv for quadrupled CO₂, implying an increase of deep water upwelling south of the circumpolar flow and an increase of Ekman pumping north of it. In addition, the changes in the wind stress and wind stress curl translate into changes in the horizontal mass transport, leading to a poleward expansion of the subtropical gyres in both hemispheres, and to strengthening of the Antarctic Circumpolar Current. Finally, the intensified near-surface winds over the Southern Ocean result in a substantial increase of mechanical energy supply to the ocean general circulation.     

Impact of waves on air-sea exchange of sensible heat and momentum

The impact of sea waves on sensible heat and momentum fluxes is described. The approach is based on the conservation of heat and momentum in the marine atmospheric surface layer. The experimental fact that the drag coefficient above the sea increases considerably with increasing wind speed, while the exchange coefficient for sensible heat (Stanton number) remains virtually independent of wind speed, is explained by a different balance of the turbulent and the wave-induced parts in the total fluxes of momentum and sensible heat.Organised motions induced by waves support the wave-induced stress which dominates the surface momentum flux. These organised motions do not contribute to the vertical flux of heat. The heat flux above waves is determined, in part, by the influence of waves upon the turbulence diffusivity.The turbulence diffusivity is altered by waves in an indirect way. The wave-induced stress dominates the surface flux and decays rapidly with height. Therefore the turbulent stress above waves is no longer constant with height. That changes the balance of the turbulent kinetic energy and of the dissipation rate and, hence the diffusivity.The dependence of the exchange coefficient for heat on wind speed is usually parameterized in terms of a constant Stanton number. However, an increase of the exchange coefficient with wind speed is not ruled out by field measurements and could be parametrized in terms of a constant temperature roughness length. Because of the large scatter, field data do not allow us to establish the actual dependence. The exchange coefficient for sensible heat, calculated from the model, is virtually independent of wind speed in the range of 3–10 ms^-1. For wind speeds above 10 ms^-1 an increase of 10% is obtained, which is smaller than that following from the constant roughness length parameterization.The investigation was in part supported by the Netherlands Geosciences Foundation (GOA) with financial aid from the Netherlands Organization for Scientific Research (NWO).     

Turbulent transport of carbon dioxide over a paddy field

Turbulent fluctuations in CO₂ concentrations over a paddy field are measured by a fastresponse device with an open sensing path. This IR device coupled with a sonic anemometer constitutes an eddy correlation instrument to measure CO₂ fluxes. Three experiments were conducted in the surface layer over paddy 90 cm high. The stability (z – d)/L ranged from -0.14 to 0.20, where L denotes the Monin-Obukhov length.CO₂ power spectra show the range of applicability of the -2/3 power law to be between f = 0.2 and f = 2, where f is the frequency normalized by wind speed and height. The cospectral estimate between CO₂ and vertical component of wind speed ranging from f = 0.005 to f = 2 shows a peak at about f = 0.15 under near-neutral stratification.Hourly means of CO₂ flux measured by the eddy correlation method increase with intensity of net radiation. The maximum value of downward flux of CO₂ rises to 0.6 mg cm^-2 hr^-1 over the paddy field at the stage of ear emergence.Some turbulence statistics relating to the CO₂ transport are evaluated: the correlation coefficient between CO₂ and vertical velocity is about -0.3, and that between CO₂ and humidity attains -0.7 ~ -0.8 under unstable stratification; nondimensional gradients _c for CO₂ and _m for wind speed are 0.89 and 0.99, respectively.     

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.     

On the bulk parameterization of surface fluxes for various conditions and parameter ranges

A comparison between various parameterizations for the bulk transfercoefficient for heat and momentum is carried out for a wide range ofatmospheric stability and values of the roughness lengths for momentum and heat,z_0m and z_0h respectively.It is confirmed that the parameterization of Launiainen compares wellto a numerical iterative solution for the Obukhov length L as function of the gradients of wind speed and temperatureover a limited range of z_0m/z_0h and stability conditions.For _0m/z_0h > 500, an alternativeinterpolation formulation of Holtslag and Ek,in combination with the formulation of Launiainen, provides a better approximation.     

Present and future near-surface wind climate of Greenland from high resolution regional climate modelling

The present and twenty-first century near-surface wind climate of Greenland is presented using output from the regional atmospheric climate model RACMO2. The modelled wind variability and wind distribution compare favourably to observations from three automatic weather stations in the ablation zone of southwest Greenland. The Weibull shape parameter is used to classify the wind climate. High values (κ > 4) are found in northern Greenland, indicative of uniform winds and a dominant katabatic forcing, while lower values (κ < 3) are found over the ocean and southern Greenland, where the synoptic forcing dominates. Very high values of the shape parameter are found over concave topography where confluence strengthens the katabatic circulation, while very low values are found in a narrow band along the coast due to barrier winds. To simulate the future (2081–2098) wind climate RACMO2 was forced with the HadGEM2-ES general circulation model using a scenario of mid-range radiative forcing of +4.5 W m^?2 by 2100. For the future simulated climate, the near-surface potential temperature deficit reduces in all seasons in regions where the surface temperature is below the freezing point, indicating a reduction in strength of the near-surface temperature inversion layer. This leads to a wind speed reduction over the central ice sheet where katabatic forcing dominates, and a wind speed increase over steep coastal topography due to counteracting effects of thermal and katabatic forcing. Thermally forced winds over the seasonally sea ice covered region of the Greenland Sea are reduced by up to 2.5 m s^?1.     

Charnock’s Roughness Length Model and Non-dimensional Wind Profiles Over the Sea

An analysis tool for the study of wind speed profiles over the water has been developed. The profiles are analysed using a modified dimensionless wind speed and dimensionless height, assuming that the sea surface roughness can be predicted by Charnock’s roughness length model. In this form, the roughness dependency on wind speed is extracted and the variations on the wind profile are due solely to atmospheric stability. The use of the Charnock’s non-dimensional wind profile is illustrated using data collected from a meteorological mast installed in the Danish North Sea. The best fit with the observed mean non-dimensional wind profile under neutral atmospheric conditions is found using a value of 1.2 × 10⁻² for Charnock’s parameter. The stability correction on the neutral wind profile suggested by the Businger-Dyer relations was found to perform well over the sea.     

Surface-Layer Characteristics in the Stable Boundary Layer with Strong and Weak Winds

An extensive meteorological dataset obtained from the plumevalidation experiment conducted by the Electric Power Research Institute (EPRI) atKincaid during 1980–1981 is analysed for studying the characteristic differences in thesurface-layer parameters in strong and weak wind stable conditions. The surface-layerparameters are computed using the similarity functions _m and _h proposed byBeljaars and Holtslag. The weak winds are characterized using the geostrophic wind speedas well as the wind speed at the 10-m level. The surface fluxes are found to be finitein weak wind conditions.Empirical formulations for the eddy diffusivities of momentum(K_M) and heat (K_H), and drag (C_D) and heat exchange (C_H) coefficients, as powerlaw functions of the bulk Richardson number (Ri_B), are proposed under both strong andweak wind conditions. Results are close to those based on observations taken from the IndianInstitute of Technology low wind diffusion experiment, the Land surface processes experiment,the Hanford diffusion experiment, the Cabauw field experiment and the Cooperative Atmospheric SurfaceExchange Study 1999 (CASES-99) experiment. In addition, the fluxes obtained fromthe proposed empirical relations are in good agreement with those based on similarity theory as wellas the turbulence measurements taken from the CASES-99 experiment.     

Turbulent exchange above a pine forest,I: Fluxes and gradients

Measurements of gradients of wind, temperature and humidity and of the corresponding turbulent fluxes have been carried out over a sparse pine forest at Jädra»s in Sweden. In order to ascertain that correct gradient estimates were obtained, two independent measuring systems were employed: one system with sensors at 10 fixed levels on a 51 m tower and another with reversing sensors for temperature and humidity, covering the height interval 23 to 32 m. Turbulent fluxes were measured at three levels simultaneously. Data from three field campaigns: in June 1985, June 1987 and September 1987 have been analyzed. The momentum flux is found on the average to be virtually constant from tree top level, at 20 to 50 m. The average fluxes of sensible and latent heat are not so well behaved. The ratio of the non-dimensional gradients of wind and temperature to their corresponding values under ideal conditions (low vegetation) are both found to be small immediately above the canopy (about 0.3 for temperature and 0.4 for wind). With increasing height, the ratios increase, but the values vary substantially with wind direction. The ratios are not found to vary systematically with stability (unstable stratification only studied). The ratio of the non-dimensional humidity gradient to the corresponding non-dimensional potential temperature gradient (equivalent to k _h /k _w ) is found to be unity for (z – d)/L _v less than about –0.1 and about 1.4 for near neutral stratification, but the scatter of the data is very large.     

Forest-Air Fluxes Of Carbon, Water And Energy Over Non-Flat Terrain

A field study of surface-air exchange of carbon, water, and energy was conducted at a mid-latitude, mixed forest on non-flat terrain to investigate how to best interpret biological signals from the eddy flux data that may be subject to advective influences. It is shown that during periods of Southwest winds (sector with mild topography), the eddy fluxes are well-behaved in terms of energy balance closure, the existence of a constant flux layer, consistency with chamber observations and the expected abiotic controls on the fluxes. Advective influences are evident during periods with wind from a steep (15%) slope to the Northeast of the tower. These influences appear more severe on CO₂ flux, particularly in stable air, than on the energy fluxes. Large positive flux of CO₂ (> 23 mol m^-2 s^-1) occurs frequently at night. The annual sum of the carbon flux is positive, but the issue about whether the forest is a source of atmospheric carbon remains inconclusive.Attempts are made to assess vertical advectionusing the data collected on a single tower. Over the Southwestsector, vertical advection makes a statistically significant but small contribution to the 30-min energy imbalance and CO₂ flux variations. Contributions by horizontal advection may be larger but cannot be verified directly by the current experimental method.     

Variations of atomic mercury concentration in atmospheric surface layer over the Ussuri Bay of the Sea of Japan during the typhoon Bolaven passage in 2012

Presented are the results of measuring the concentration of atomic mercury Hg⁰ on August 28–30, 2012 in the atmospheric surface layer over the Ussuri Bay in the Sea of Japan during the typhoon Bolaven passage. It is revealed that the Hg⁰ concentration during this period varied from 1.7 to 3.3 ng/m³. The maximum values were observed at the coming of mercury-enriched air masses from the Yellow Sea area at the simultaneous decrease in the air pressure and increase in the wind speed.