Effects of Mesoscale Surface Thermal Heterogeneity on Low-Level Horizontal Wind Speeds

Using large-eddy simulation, we investigate characteristics of horizontal wind speed at 100 m above the ground, with surface heat-flux variations that are sinusoidal with amplitudes of 0, 50, and 200 W m⁻² and wavelengths of 16, 32, and 128 km, and no background flow. When the amplitude is 200 W m⁻², wind speeds induced by the surface-flux variations on scales of 16 and/or 32 km have multiple temporal oscillations from 0600 to 1800 local standard time. The positive peaks first appear before noon. In contrast, for wind speeds induced by the 128-km surface heterogeneity, a single oscillation occurs in the late afternoon, which is much larger than those generated by the 16- and 32-km surface heterogeneity. In addition, at the oscillation onset the kurtosis of the velocity increment over a distance of 1 km significantly increases, which implies intermittency in the generation of 1-km scale eddies. The spatially intermittent energy cascade generated by surface heterogeneity scaled down to 1-km eddies is analogous to the well-known intermittent energy cascade in the inertial subrange. The kurtosis of the 1-km eddies is much larger with the 128-km surface heterogeneity than with the 16- and 32-km heterogeneities. Thus we conclude that localized rapid changes of low-level horizontal wind speed may be caused by significant local surface heterogeneity on scales between a few tens and a few hundreds of kilometres.     

An Improved Approach for Parameterizing Surface-Layer Turbulent Transfer Coefficients in Numerical Models

Based on classic iterative computation results, new equations to calculate the surface turbulent transfer coefficients are proposed, which allow for large ratios of the momentum and heat roughness lengths. Compared to the Launiainen scheme, our proposed scheme generates results closer to classical iterative computations. Under unstable stratification, the relative error in the Launiainen scheme increases linearly with increasing instability, even exceeding 15%, while the relative error of the present scheme is always less than 8.5%. Under stable stratification, the Launiainen scheme uses two equations, one for 0 < Ri _B ≤ 0.08 and another for 0.08 < Ri _B ≤ 0.2, and does not consider the condition that Ri _B > 0.2, while its relative errors in the region 0 < Ri _B ≤ 0.2 exceed 31 and 24% for momentum and heat transfer coefficients, respectively. In contrast, the present scheme uses only one equation for 0 < Ri _B ≤ 0.2 and another equation for Ri _B > 0.2, and the relative error of the present scheme is always less than 14%.     

A study of ozone in the Colorado mountains

The seasonal and diurnal variations of ozone mixing ratios have been observed at Niwot Ridge. Colorado. The ozone mixing ratios have been correlated with the NO _x (NO+NO₂) mixing ratios measured concurrently at the site. The seasonal and diurnal variations in O₃ can be reasonably well understood by considering photochemistry and transport. In the winter there is no apparent systematic diurnal variation in the O₃ mixing ratio because there is little diurnal change of transport and a slow photochemistry. In the summer, the O₃ levels at the site are suppressed at night due to the presence of a nocturnal inversion layer that isolated ozone near the surface, where it is destroyed. Ozone is observed to increase in the summer during the day. The increases in ozone correlate with increasing NO _x levels, as well as with the levels of other compounds of anthropogenic origin. We interpret this correlation as in-situ or in-transit photochemical production of ozone from these precursors that are transported to our site. The levels of ozone recorded approach 100 ppbv at NO _x mixing ratios of approximately 3 ppbv. Calculations made using a simple clean tropospheric chemical model are consistent with the NO _x -related trend observed for the daytime ozone mixing ratio. However, the chemistry, which does not include nonmethane hydrocarbon photochemistry, underestimates the observed O₃ production.     

Doppler Lidar Measurements of Vertical Velocity Spectra in the Convective Planetary Boundary Layer

We utilized a Doppler lidar to measure spectra of vertical velocity w from 390m above the surface to the top of the daytime convective boundary layer (CBL). The high resolution 2μm wavelength Doppler lidar developed by the NOAA Environmental Technology Laboratory was used to detect the mean radial velocity of aerosol particles. It operated continuously during the daytime in the zenith-pointing mode for several days in summer 1996 during the Lidars-in-Flat-Terrain experiment over level farmland in central Illinois, U.S.A. The temporal resolution of the lidar was about 1 s, and the range-gate resolution was about 30m. The vertical cross-sections were used to calculate spectra as a function of height with unprecedented vertical resolution throughout much of the CBL, and, in general, we find continuity of the spectral peaks throughout the depth of the CBL. We compare the observed spectra with previous formulations based on both measurements and numerical simulations, and discuss the considerable differences, both on an averaged and a case-by-case basis. We fit the observed spectra to a model that takes into account the wavelength of the spectral peak and the curvature of the spectra across the transition from low wavenumbers to the inertial subrange. The curvature generally is as large or larger than the von Kármán spectra. There is large case-to-case variability, some of which can be linked to the mean structure of the CBL, especially the mean wind and the convective instability. We also find a large case-to-case variability in our estimates of normalized turbulent kinetic energy dissipation deduced from the spectra, likely due for the most part to a varying ratio of entrainment flux to surface flux. Finally, we find a relatively larger contribution to the low wavenumber region of the spectra in cases with smaller shear across the capping inversion, and suggest that this may be due partly to gravity waves in the inversion and overlying free atmosphere.     

Parameterizations of the daytime friction velocity,temperature scale,and upslope flow over gently inclined terrain in calm synoptic conditions

A set of new parameterizations for the friction velocity and temperature scale over gently sloped terrain and in calm synoptic conditions are theoretically derived. The friction velocity is found to be proportional to the product of the square root of the total accumulated heating in the boundary layer and the sinusoidal function of the slope angle, while the temperature scale is proportional to the product of the boundary layer depth, the sinusoidal function of the slope angle and the potential temperature gradient in the free atmosphere. Using the new friction velocity parameterization, together with a parameterization of eddy diffusivity and an initial potential temperature profile around sunrise, an improved parameterization for the thermally induced upslope flow profile is derived by solving the Prandtl equations. The upslope flow profile is found to be simply proportional to the friction velocity.     

Atmospheric Disturbances that Generate Intermittent Turbulence in Nocturnal Boundary Layers

Using the unprecedented observational facilities deployed duringthe 1999 Cooperative Atmosphere-Surface Exchange Study (CASES-99),we found three distinct turbulent events on the night of 18October 1999. These events resulted from a density current,solitary wave, and internal gravity wave, respectively. Our studyfocuses on the turbulence intermittency generated by the solitarywave and internal gravity wave, and intermittent turbulenceepisodes associated with pressure change and wind direction shiftsadjacent to the ground. Both the solitary and internal gravitywaves propagated horizontally and downward. During the passage ofboth the solitary and internal gravity waves, local thermal andshear instabilities were generated as cold air was pushed abovewarm air and wind gusts reached to the ground. These thermal andshear instabilities triggered turbulent mixing events. Inaddition, strong vertical acceleration associated with thesolitary wave led to large non-hydrostatic pressure perturbationsthat were positively correlated with temperature. The directionaldifference between the propagation of the internal gravity waveand the ambient flow led to lateral rolls. These episodic studiesdemonstrate that non-local disturbances are responsible for localthermal and shear instabilities, leading to intermittentturbulence in nocturnal boundary layers. The origin of thesenon-local disturbances needs to be understood to improve mesoscalenumerical model performance.     

The spectral velocity tensor for homogeneous boundary-layer turbulence

We have postulated a simple model for the spectral tensor _ij (k) of an anisotropic, but homogeneous turbulent velocity field. It is a simple generalization of the spectral tensor _inf ^{ij piso}(k) for isotropic turbulence and we show how in the limit of isotropy, _ij (k) becomes equal to _inf ^{ij piso}(k). Whereas _inf ^{ij piso}(k) is determined entirely by one scalar function of k = ¦k¦, namely the energy spectrum, we need three independent scalar functions of k to specify _ij (k). We show how it is possible by means of the three stream-wise velocity component spectra to determine the three scalar functions in _ij (k) by solving two uncoupled, ordinary linear differential equations of first and second order. The analytic form of the component spectra each has a set of three parameters: the variance and the integral length scale of the velocity component and a dimensionless parameter, which governs the curvature of the spectrum in the transition domain from the inertial subrange towards lower wave numbers. When the three sets of parameters are the same, the three spectra correspond to isotropic turbulence and they are all interrelated and related to the energy spectrum. We show how it is possible to obtain these spectral forms in the neutral surface layer and in the convective boundary layer from data reported in the literature. The spectral tensor is used to predict the lateral coherences for all three velocity components and these predictions are compared with coherences obtained in two experiments, one using three masts at a horizontally homogeneous site in Denmark and one employing two aircraft flying in formation over eastern Colorado. Comparison shows reasonable agreement although with considerable experimental scatter.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

An analytic formulation for NO and NO2 flux profiles in the atmospheric surface layer

The chemical reactivity of NO and NO₂ is so rapid that their fluxes and concentrations can be considerably modified from that expected for conserved variables in the atmospheric surface layer, even as low as a meter above the surface. Fitzjarrald and Lenschow (1983) have calculated flux and mean concentration profiles for NO, NO₂ and O₃ in the surface layer using numerical techniques. However, their solutions do not approach the photostationary state at large heights. Here we solve a simpler set of equations analytically (i.e. we assume a constant O₃ concentration and neutral hydrodynamic stability), and are able to show how the flux profiles behave at large heights assuming that the concentrations approach their photostationary values. We find, for example, that at large heights the ratio of the flux of NO to that of NO₂ is equal to the ratio of their concentrations. These results are relevant to estimating surface fluxes of NO and NO₂, and are most applicable to nonurban environments where NO and NO₂ concentrations are usually much less than O₃ concentration.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Boundary-Layer Turbulence Over The Nebraska Sandhills

Data from National Centre for Atmospheric Research (NCAR) Queen Air boundary-layer flights over the Nebraska Sandhills are analyzed to investigate the effects of these low hills on boundary-layer turbulence. The Sandhills are an area of anisotropic rolling terrain with characteristic wavelengths of order 2km and rms height variations of order 25m. The biggest impact is found in early morning flight data where horizontal velocity perturbations appear at the same wavelengths as the terrain and variances (normalised by u ² , where u is the local friction velocity) are significantly enhanced relative to standard flat terrain values. By contrast the vertical velocity variance seems less affected and terrain effects are much less evident in data from the afternoon convective boundary layer.