Flux-profile Relationships for Wind Speed and Temperature in the Stable Atmospheric Boundary Layer

Wind and temperature profiles in the stable boundary layer were analyzed in the context of MoninObukhov similarity. The measurements were made on a 60-m tower in Kansas during October 1999 (CASES-99). Fluxprofile relationships, obtained from these measurements in their integral forms, were established for wind speed and temperature. Use of the integral forms eliminates the uncertainty and accuracy issues resulting from gradient computations. The corresponding stability functions, which were nearly the same for momentum and virtual sensible heat, were found to exhibit different features under weakly stable conditions compared to those under strongly stable conditions. The gradient stability functions were found to be linear, namely _m = 1+ 5.8 and _h = 1 + 5.4 up to a limit of the MoninObukhov stability parameter = 0.8; this is consistent with earlier findings. However, for stronger stabilities beyond a transition range, both functions were observed gradually to approach a constant, with a value of approximately 7. To link these two distinct regimes, a general but pliable functional form with only two parameters is proposed for the stability functions, covering the entire stability range from neutral to very stable conditions.     

On the similarity in atmospheric fluctuations of carbon dioxide,water vapor and temperature over vegetated fields

This paper describes the similarity between atmospheric fluctuations of carbon dioxide, water vapor and temperature using data which cover a wide range of instability (0.02 < < 10). The is the Monin-Obukhov stability parameter including the humidity effect.The spectral analysis shows that the coherency between fluctuations of carbon dioxide and water vapor or temperature is very close to unity, and the phase difference is basically out of phase for whole frequency ranges analyzed. The stability dependence of the normalized standard deviation of carbon dioxide is very similar to those of water vapor and temperature. The normalized standard deviation is about 2.5 under near neutral conditions, and it decreases with increasing instability following the -1/3; power law as (-)^-1/3. The skewness factors of carbon dioxide, water vapor and temperature show a systematic departure with increasing instabilities for 0.02 < s- < 1, and level off at high instabilities for 1 < -\s < 10. The stability dependence of the flatness factors is not so clear as that noted in the standrard deviation and skewness factors. Dissipation rates of carbon dioxide, water vapor and temperature variance are well related to the spectral peak wavelength. This seems to be real since the local production and local dissipation rates are the main terms, almost balancing one another in the variance budget equations for scalar entities.     

The Value Of The β Coefficient In The Relaxed Eddy Accumulation Method In Terms Of Fourth-Order Moments

The commonly measured value of in the relaxed eddy accumulationmethod of about 0.56is shown to arise from the non-Gaussiannature of turbulence. Fourth-orderGram–Charlier functions forthe two-dimensional probability distributionsof variation in the horizontal component of wind velocityand concentrations of water vapour, carbondioxide and methane with respect to thevertical component of wind velocity are used to examinethe value of .An analytical solution for ispresented in terms of fourth-order moments.Under mean conditions, this solution givesa value for of0.557. Variation of is shown to be controlledprimarily by the ratio of the mean ofc'w3 (where c'is relevant to the entity of interest andw' is vertical component of windvelocity) to the correlationcoefficient between the entity concentrationand vertical component of wind velocity.     

Markov-chain simulation of particle dispersion in inhomogeneous flows: The mean drift velocity induced by a gradient in Eulerian velocity variance

The Langevin equation is used to derive the Markov equation for the vertical velocity of a fluid particle moving in turbulent flow. It is shown that if the Eulerian velocity variance _wE is not constant with height, there is an associated vertical pressure gradient which appears as a force-like term in the Markov equation. The correct form of the Markov equation is: w(t + t) = aw(t) + b _wE + (1 – a)T _L ( _wE ²)/z, where w(t) is the vertical velocity at time t, a random number from a Gaussian distribution with zero mean and unit variance, T _L the Lagrangian integral time scale for vertical velocity, a = exp(–t/T _L), and b = (1 – a ²)^1/2. This equation can be used for inhomogeneous turbulence in which the mean wind speed, _wE and T _L vary with height. A two-dimensional numerical simulation shows that when this equation is used, an initially uniform distribution of tracer remains uniform.     

The budget of carbon dioxide variance in the surface layer over vegetated fields

The budget equation for carbon dioxide variance can be represented by production, dissipation and flux divergence terms. Each term is measured under near neutral to moderately unstable conditions over vegetated fields. The flux divergence term is about an order of magnitude smaller than production and dissipation terms, though it shows a loss for 0.006 < _v < 1 and a gain for 1 < - _v < 10. Here, _v is the Monin-Obukhov stability parameter including humidity effect. As expected from a closure of the budget, the nondimensional production and dissipation terms are basically identical and represented by the same functional form: (1–16 _v )^–1/2.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

Determination of similarity functions of the resistance laws for the planetary boundary layer using surface-layer similarity functions

Functional forms of the universal similarity functions A, B (for wind components parallel and normal to the surface stress), and C (for potential temperature difference) are determined based on the generalized theory of the resistance laws for the Planetary Boundary Layer (PBL). The similarity-profile functions for the surface layer are matched with the velocity and temperature-defect profiles that are assumed to have shapes modified by certain powers of nondimensional height z/h, where h is the PBL height. The powers of the outer-layer profile functions are determined, so that the functions become negligible in the surface layer. To close the temperature defect law, an assumption that the temperature gradient across the top of the PBL is continuous with the stratification of the overlying atmosphere is used. The result of this assumption is that nondimensional momentum and temperature profiles in the PBL can be described in terms of four basic ratios: (1) roughness ratio = /h (2) scale-height ratio =|f|h/u_*, (3) ambient stratification parameter =h/_*, and (4) stability parameter =h/L, where L is the Monin-Obukhov length, z₀ is the surface roughness, is the upper-air stratification, u _* is the friction velocity, and _* is the temperature scale at the surface. For stable conditions, the scale-height ratio can be related to the atmospheric stability and the upperair stratification, and the generalized similarity and Rossby number similarity theories become identical. Under appropriate boundary conditions, function A is explicitly dependent on the stability parameter , while B is a function of scale-height ratio , which in turn depends on the stability. Function C is shown to be dependent on the stability and the upper-air stratification, due to the closure assumption used for the temperature profile.The suggested functional forms are compared with other empirical approximations by several authors. The general framework used to determine the functional forms needs to be tested against good boundary-layer measurements.     

Spectral Characteristics and Correction of Long-Term Eddy-Covariance Measurements Over Two Mixed Hardwood Forests in Non-Flat Terrain

We present turbulence spectra and cospectra derived from long-term eddy-covariancemeasurements (nearly 40,000 hourly data over three to four years) and the transferfunctions of closed-path infrared gas analyzers over two mixed hardwood forests inthe mid-western U.S.A. The measurement heights ranged from 1.3 to 2.1 times themean tree height, and peak vegetation area index (VAI) was 3.5 to 4.7; the topographyat both sites deviates from ideal flat terrain. The analysis follows the approach ofKaimal et al. (Quart. J. Roy. Meteorol. Soc. 98, 563–589, 1972) whose results were based upon 15 hours of measurements atthree heights in the Kansas experiment over flatter and smoother terrain. Both thespectral and cospectral constants and stability functions for normalizing and collapsingspectra and cospectra in the inertial subrange were found to be different from those ofKaimal et al. In unstable conditions, we found that an appropriate stabilityfunction for the non-dimensional dissipation of turbulent kinetic energy is of the form ₍₎ = (1 - b^-)^-1/4 - c^-, where representsthe non-dimensional stability parameter. In stable conditions, a non-linear functionG_xy() = 1 + b_xy^c _xy (c_xy < 1) was found to benecessary to collapse cospectra in the inertial subrange. The empirical cospectralmodels of Kaimal et al. were modified to fit the somewhat more (neutraland unstable) or less (stable) sharply peaked scalar cospectra observed over forestsusing the appropriate cospectral constants and non-linear stability functions. Theempirical coefficients in the stability functions and in the cospectral models varywith measurement height and seasonal changes in VAI. The seasonal differencesare generally larger at the Morgan Monroe State Forest site (greater peak VAI) andcloser to the canopy.The characteristics of transfer functions of the closed-path infrared gas analysersthrough long-tubes for CO₂ and water vapour fluxes were studied empirically. This was done by fitting the ratio between normalized cospectra of CO₂ or watervapour fluxes and those of sensible heat to the transfer function of a first-order sensor.The characteristic time constant for CO₂ is much smaller than that for water vapour. The time constant for water vapour increases greatly with aging tubes. Three methods were used to estimate the flux attenuations and corrections; from June through August, the attenuations of CO₂ fluxes are about 3–4% during the daytime and 6–10% at night on average. For the daytime latent heat flux (Q_E), the attenuations are foundto vary from less than 10% for newer tubes to over 20% for aged tubes. Correctionsto Q_E led to increases in the ratio (Q_H + Q_E)/(Q^* - Q_G) by about 0.05 to0.19 (Q_H is sensible heat flux, Q^* is net radiation and Q_G is soil heat flux),and thus are expected to have an important impact on the assessment of energy balanceclosure.     

Numerical Simulation of Strongly Stratified Flow Over a Three-Dimensional Hill

A numerical study of stably stratified flow over a three-dimensional hill is presented. Large-eddy simulation is used here to examine in detail the laboratory experimental flows described in the landmark work of Hunt and Snyder about stratified flow over a hill. The flow is linearly stratified and U/Nh is varied from 0.2 to 1.0. Here N and U are the buoyancy frequency and freestream velocity respectively, and h is the height of the hill. The Reynolds number based on the hill height is varied from 365 to 2968. The characteristic flow patterns at various values of U/Nh have been obtained and they are in good agreement with earlier theoretical and experimental results. It is shown that the flow field cannot be predicted by Drazin's theory when recirculation exists at the leeside of the hill even at UNh 1. The wake structure agrees well with a two-dimensional wake assumption when U/Nh 1 but lee waves start to influence the wake structure as U/Nh increases. The dividing-streamline heights obtained in the simulation are in accordance with experimental results and Sheppard's formula. The energy loss along the dividing streamline due to friction/turbulence approximately offsets the energy gained from pressure field. When lee waves are present, linear theory always underestimates the amplitude and overestimates the wavelength of three-dimensional lee waves. The simulated variations of drag coefficients with the parameterK (=ND/ U) are qualitatively consistent with experimental data and linear theory. Here D is the depth of the tank.     

Non-dimensional wind and temperature profiles in the atmospheric surface layer: A re-evaluation

Previous results of non-dimensional wind and temperature profiles as functions of ( = z/L) show systematic deviations between different experiments. These discrepancies are generally believed not to reflect real differences but rather instrumental shortcomings. In particular, it is clear that flow distortion has not been adequately treated in most previous experiments. In the present paper, results are presented from a surface-layer field experiment where great care was taken to remove any effects from this kind of error and also to minimize other measuring errors. Data from about 90 30-min runs with turbulence measurements at three levels (3, 6, and 14 m) and simultaneous profile data have been analysed to yield information on flux-gradient relationships for wind and temperature.The flux measurements themselves show that the fluxes of momentum and sensible heat are constant within ± 7% on average for the entire 14 m layer in daytime conditions and when the stratification is slightly stable. For more stable conditions, the flux starts to decrease systematically somewhere in the layer 6 to 14 m. From a large body of data for near-neutral conditions (¦¦ 0.1), values are derived for von Kármán's constant: 0.40 ± 0.01 and for _h at neutrally, 0.95 ± 0.04. The range of uncertainty indicated here is meant to include statistical uncertainty as well as the effect of possible systematic errors.Data for _m and _h for an extended stability range (1 > > – 3) are presented. Several formulas for _m and _h appearing in the literature have been used in a comparative study. But first all the formulas have been modified in accordance with the following assumptions: = 0.40 and ( _h ) _{= 0} = 0.95; deviations from this result in the various studies are due to incomplete correction for flow distortion. After new corrections are introduced, the various formulas were compared with the present measurements and with each other. It is found that after this modification, the most generally used formulas for _m and _h for unstable conditions, i.e., those of Businger et al. (1971) and Dyer (1974) agree with each other to within ± 10% and with the present data. For stable conditions, the various formulas still disagree to some extent. The conclusion in relation to the present data is not as clear as for the unstable runs, because of increased scatter. It is, however, found that the modified curve of Businger et al. (1971) for _h fits the data well, whereas for _m , Dyer's (1974) curve appears to give slightly better agreement.     

Bistability of the thermohaline circulation identified through comprehensive 2-parameter sweeps of an efficient climate model

The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (F_a) and atmospheric moisture diffusivity (K_q). In a large ensemble of spinup experiments, the values of F_a and K_q are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of F_a and K_q, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.     

Synchronous Scan and Excitation-Emission Matrix Fluorescence Spectroscopy of Water-Soluble Organic Compounds in Atmospheric Aerosols

Three-dimensional excitation–emission matrix (EEM) fluorescence spectra of water-soluble organic compounds (WSOC) from aerosol samples were measured and compared with those reported in the literature for natural dissolved organic matter. The EEM profiles of the WSOC presented three characteristic excitation/emission (_Exc/_Em) peaks: 240/405 nm, 310/405 nm and 280/340 nm. The fluorescence intensities at _Exc/_Em240/405 nm and _Exc/_Em310/405 nm are located at wavelengths shorter than those reported for aquatic humic substances, indicating a smaller content of both aromatic structures and condensed unsaturated bond systems in the WSOC fraction. The EEM profiles of fractions obtained by the isolation procedure of the WSOC by the XAD resins showed that a fractionation has occurred and the XAD-8 eluate is highly representative of the total WSOC of collected aerosol. Synchronous scan spectra were more detailed than conventional fluorescence emission spectra, appearing more suitable for studying multicomponent samples such as the WSOC from atmospheric aerosols.     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.     

An improved calibration for tethered balloon wind measurements

A previously published technique for using tethered spherical balloons as anemometers for measuring light low-level winds has been further developed. Earlier data on the relationship between the aerodynamic drag coefficient and the Reynolds number of spherical rubber balloons were combined with a large number of new data and re-analysed; and the errors in the relationship were estimated. The results allowed a more accurate calculation of wind speed from the deflection of a tethered balloon from the vertical. When combined with a new technique for calculating the effects of the tether, this enabled light to moderate low-level winds at fixed heights up to 600 m or more to be measured with simple, cheap, and readily mobile equipment; and a slight modification of the technique allowed measurement of winds in and above fog. Wind speeds measured by the ballon technique showed reasonably good agreement with measurements by an anemometer carried beneath the balloon.Glossary of Symbols a, b, c Coefficients in the relationship between lnC _d and lnR - A Quantity under square root in solution for lnV whena0 - C _d Wind drag coefficient for balloon - C _dc Value ofC _d given by calibration curve of Table I - D Dynamic wind pressure force on balloon - F Buoyant free lift of balloon with load - Re Reynold's number of balloon (sphere) - R = Re/10⁵ - r Radius of sphere - T Tension in tether - V Wind speed - ₈₃() =(lnC _dc -lnC _d ) when 83° , or 0 for other - Error in lnC _d - Elevation of tether where attached to balloon - Elevation of balloon from ground tether point - Molecular viscosity of air - Ratio of circumference to diameter of circle - Density of air     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

Analysis and trends of atmospheric turbidity parameters over Cairo

Summary This work deals with the Linke turbidity factor, based on total spectrum observations of the direct solar beam and aerosol turbidity parametersa _a , , and based on observations in broad spectral bands. Diurnal and seasonal variations of these turbidity parameters were analyzed for the period 1975 to 1991.Annual variations of these parameters show low values in winter and high values in both spring and summer. The extinction coefficients decrease with increase of both wavelength and optical airmass. Trend analysis shows an increase in aerosol extinction coefficient below 0.63 m, and a slight decrease for longer wave-lengths.Linear regression relations are also constructed to estimate botha _a and whenT _L is available. The relations show thata _a can be estimated with errors below 20%. The relation with the parameter, may give better results when it is estimated by assigning a fixed value of .Nomenclature AV Monthly and total average of the measured parameter - a Atmospheric optical thicknes - a _a Aerosol optical thickness - a _r Mean value of optical thickness of. Rayleigh atmosphere over all wavelengths - a _o Ozone optical thickness - a _o Ozone absorption coefficient - a _w Water vapor optical thickness - COR Correlation coefficient of the linear relation in percentage - Ex₁, Ex₂, Ex₃ Aerosol extinction coefficients in the bands .2–.53, .53–.63, .63–.695, respectively - I _(o) Normal incident direct solar radiation, under clear sky condition - I _o Extraterrestrial insolation at normal incidence - m _r Relative (optical) air mass - NO Number of the observations used in either making the relation or the verification - RMSE Root mean square error of the linear relation - RMSE% Percentage value of the root mean square error relative to the average measured value AV - T _L Linke turbidity factor - T Dry bulb temperature in °C - u _o Ozone layer thickness, cm - Z Zenith angle - Ångström wavelength exponent - Ångström turbidity coefficient - Wavelength - Y The year number after 1975 With 5 Figures     

Equilibrium evaporation beneath a growing convective boundary layer

Expressions for the equilibrium surface Bowen ratio ( _s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer _i and the entrainment parameter A _R . If A_R is put equal to zero, the solution for _s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA _R and _i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.     

Turbulence in a Katabatic Flow

Turbulence measurements performed in a stable boundary layer over the sloping ice surface of the Vatnajökull in Iceland are described. The boundary layer, in which katabatic forces are stronger than the large-scale forces, has a structure that closely resembles that of a stable boundary layer overlying a flat land surface, although there are some important differences. In order to compare the two situations the set-up of the instruments on an ice cap in Iceland was reproduced on a flat grass surface at Cabauw, the Netherlands. Wind speed and temperature gradients were calculated and combined with flux measurements made with a sonic anemometer in order to obtain the local stability functions _m and _h as a function of the local stability parameter z/L. Unlike the situation at Cabauw, where _m was linear as a function of z/L, in the katabatically forced boundary layer, the dependence of _m on stability was found to be non-linear and related to the height of the wind maximum. Thermal stratification and the depth of the stable boundary layer however seem to be rather similar under these two different forcing conditions.Furthermore, measurements on the ice were used to construct the energy balance. These showed good agreement between observed melt and components contributing to the energy balance: net radiation (supplying 55% of the energy), sensible heat flux (30%) and latent heat flux (15%).Local sources and sinks in the turbulent kinetic energy budget are summed and indicate a reasonable balance in near-neutral conditions but not in more stable situations. The standard deviation of the velocity fluctuations _u, _v, and _w, can be scaled satisfactorily with the local friction velocity u_* and the standard deviation of the temperature fluctuation with the local temperature scale _*.