Model experiments on free-convection heat and mass transfer of leaves and plant elements

Effects of orientation, shape, and surface structure on the free convection mass transfer of plates, leaf- and plant-models were studied in an electrolytic system. The results were extrapolated to the transfer of heat and mass in air, in an effort to predict realistic boundary-layer transfer coefficients for leaves and plant-like surfaces. Flow visualization complemented the investigation.The results indicate that increases in transfer of 50% over that of a vertical plate are possible for non-vertical rough leaves, but that a 25% increase is a more reasonable estimate for an average plant. Evidence is also presented that more reliable caluclations of the transfer properties of pine needles are obtained if the needles are approximated by vertical cylinders, rather than by horizontal ones.     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

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.     

Observed variations of σθ and σΦ in the nocturnal drainage flow in a deep valley

The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms^-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and     

Recurrence statistics of concentration fluctuations in plumes within a near-neutral atmospheric surface layer

This study examines the statistical properties of the concentration derivative, , for a dispersing plume in a near-neutrally stratified atmospheric surface layer. Towards this goal, the probability density function (pdf) of , and the conditional pdf of given a fixed concentration level, , have been measured. These pdfs are found to be modeled well by a generalizedq-Gaussian (gqG) distribution with intermittency exponent,q, equal to 0.3 and 3/4, respectively. These results highlight the strong intermittency effect (patchiness) of the small-scale concentration eddy structures in the plume. The distribution of time intervals between successive high peaks in the squared derivative process, x², is found to be well approximated by a power-law distribution, implying that occurrences of these high peaks are much more clustered than would be predicted by a Poisson or shot-noise process. The results are used to improve models for the joint pdf of and , and for the expected number of upcrossings per unit time interval of a fixed concentration level that have been proposed by Kristensenet al. (1989). The predictions of the improved models are in accord with observations, and suggest that the intercorrelation between and must be explicitly incorporated if good estimates of the upcrossing intensity are to be obtained.     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

AN E – ε – ℓ TURBULENCE CLOSURE SCHEME FOR PLANETARY BOUNDARY-LAYER MODELS: THE NEUTRALLY STRATIFIED CASE

The standard E – model generates aplanetary boundary layerthat appears to be much too deep. The cause of theproblem is traced to the equation for the dissipationrate () of turbulent kinetic energy (E), specifically theparameterization of dissipation production anddestruction. In the context of atmosphericboundary-layer modelling, we argue that a part of thedissipation production should be modelled as the inputto the spectral cascade from the energy-containingpart of the spectrum, with a characteristic length , while the equilibrium imbalancebetween local production and destruction ofdissipation is modelled as proportional toE2/E, as in the standard model. Wepropose an E – – turbulence closurescheme, in which both the mixing length, m, and are prescribed. The importance ofthe equation is diminished, though itstill determines the dissipation rate in the Eequation.     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

Properties of spatial temperature derivatives in the atmospheric surface layer

A pair of parallel cold wires separated in either the vertical or lateral direction was used to obtain the three components _x, _y, _z of the temperature derivative in the streamwise, lateral and vertical directions, respectively. The average absolute skewness values of _x and _z are nonzero and approximately equal, while the skewness of _y is approximately zero. These results appear to be consistent with the presence of a large, three-dimensional organised structure in the surface layer. There is an apparent low-frequency contamination in the spectral density of _y and _z due mainly to small errors in estimating the sensitivity of the cold wires. The temperature derivatives were high-pass filtered, the filter being set to remove possible contributions from the large structure and to minimise low-frequency sensitivity contamination. The filtered rms ratios \~_x/\~_y and \~_x/\~_z were in the range 0.7 to 0.9, a result in qualitative agreement with that obtained in the laboratory boundary layer by Sreenivasan et al. (1977). The skewness of filtered _x or _z is negligible, consistent with local isotropy of small-scale temperature fluctuations and in support of the high wavenumber spectral isotropy discussed in Antonia and Chambers (1978).     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.     

The Effect of Atmospheric Stability on the Surface-Layer Characteristics in a Low-Wind Area of Tropical West Africa

This study details the observed effects ofatmospheric stability on characteristics of thesurface layer in a low wind speed (U = 1.5 m s-1)regime of tropical West Africa. Theaerodynamic roughness length, z0, anddisplacement height, d, obtained from profilewind-speed data at our bush land site (height 2 m)have values of 0.24 ± 0.10 m and 1.54 ± 0.04 mrespectively. In the unstable range (-2.5 < Ri < -0.1; Riis gradient Richardson number), thestandard deviation in wind speed fluctuations, u, increased from 0.57 ± 0.19 m s-1 toa maximum of 0.7 ± 0.2 m s-1 in near-neutralconditions, and in the stable range, the parameterdecreased rapidly to 0.41 ± 0.15 m s-1 at Ri 0.2.In the same stability range, the horizontal winddispersion, , decreased withincreasing stability from 19 ± 8 deg. to 13 ± 5 deg.The surface-layer integral quantity, u/u*, when plottedas a function of stability, is in agreement with theempirical results. The ratio ofsensible heat flux (estimated) to the net radiationranged between 0.1 and 0.2 at nighttime,increasing to about 0.5 during the daytime, and showeda strong dependency on season.     

Gas and Particulate Products Distribution from the Photooxidation of α-Humulene in the Presence of NOx, Natural Atmospheric Air and Sunlight

The photooxidation of -humulene in the presence of NO_x, natural sunlight, and rural background air was investigated using a combination of gas chromatography-mass spectrometry (GC-MS) and high performance liquid chromatography (HPLC). Identification and quantification of gas and particulate reaction products were reported over the course of the reaction. The daytime photooxidation was carried out in a large outdoor smog chamber (190 m³). A wide range of ring retaining and ring opening products in the gas and particle phase are reported. On average, measured gas and particle phase products accounted for 44% of the reacted -humulene carbon. Measurements show that a number of reaction products with low vapor pressures (e.g. 3-seco--humulone aldehyde, 7-seco--humulone aldehyde, -humulal aldehyde, -humulene 3-oxide or -humulene 7-oxide, -humulaic/alic acid isomers, and 3-seco--14-hydroxyhumulone aldehyde) were found in the early stage of the reaction and may play an important role in the early formation of secondary organic aerosol. A detailed mechanism is proposed to account for most products observed in this investigation.     

On the cumulative frequency distributions ofk-day period amounts of precipitation for any station in The Netherlands,while k≧30

Summary A simple method of representing cumulative frequency distributions ofk-day period amounts of precipitation (30k1800), deduced from observations at Hoofddorp during 87 years and from observations at Winterswijk during 73 years, by means of cumulative Poisson distributions, had been developed. Poisson's constant appeared to depend linearly onk only, while the Poisson unitq appeared to be only proportional to the annual variation.In this paper it is shown that the representation of observed cumulative frequency distributions with the aid of Poisson distributions is true for any station in The Netherlands in the interval 30k1800. A map with isopleths of -values fork=360 presents a simple picture, so that for any station an -value fork=360 can be estimated. For the process of representation the meank-day period amount of precipitation is also necessary. This value may be obtained from a map for each of the 12 months, each map covered with isohyetal lines, which give a distribution of the monthly precipitation.

---

Zusammenfassung Es wurde eine Methode entwickelt, um in einfacher Weise Summenfunktionen vonk-tägigen Niederschlagsmengen (30k1800), welche aus Beobachtungen in Hoofddorp während 87 Jahren und in Winterswijk während 73 Jahren angestellt worden sind, durch kumulative Poissonsche Häufigkeitsverteilungen darzustellen. Dabei ergab sich, daß die Poissonsche Konstante nur linear vonk abhängig ist, während die Poissonsche Einheitq sich als proportional zum Jahresgang erwies.In der vorliegenden Untersuchung wird nun gezeigt, daß die Darstellung der beobachteten kumulativen Häufigkeitsverteilungen mit Hilfe von Poisson-Verteilungen für irgendeine Station in den Niederlanden im Bereich von 30k1800 gültig ist. Eine Karte mit Isoplethen der -Werte fürk=360 ergibt ein einfaches Bild, aus dem für jede beliebige Station ein -Wert fürk=360 abgelesen werden kann. Für eine solche Darstellung muß auch die Niederschlagsmenge derk-tägigen Perioden bekannt sein; diese Werte erhält man aus monatlichen Isohyeten-Karten, die die monatliche Niederschlagsverteilung wiedergeben.

---

Résumé On avait développé une méthode pour représenter d'une façon simple les distributions des sommes cumulées dek jours (30k1800) des quantités de précipitations déduites des observations faites à Hoofddorp pendant 87 ans et à Winterswijk pendant 73 ans. On avait démontré que la constante de Poisson ne dépend dek que linéairement et que l'unité de Poissonq est proportionnelle à la variation annuelle.De la présente étude ressort que, dans l'intervalle de 30k1800, la représentation des distributions des fréquences cumulées observées au moyen des formules de Poisson est valable pour une station quelconque aux Pays-Bas. Une carte d'isoplèthes de la valeur pourk=360 donne une image suffisamment simple pour que la valeur d'une station quelconque aux Pays-Bas puisse être déterminée sans difficulté. Pour une telle représentation il faut connaître également les moyennes des quantités des périodes dek jours de la précipitation. Ces valeurs peuvent se déduire de cartes mensuelles des isohyètes moyennes.

---

---

With 1 Figure

---

Dedicated to Dr.Anders K. Ångström on the occasion of his 70th birthday.     

Mass transfer to bean leaves

A vapour of radio-lead (²¹²Pb) has been used to measure the Sherwood number, Sh, of model leaves at various angles of incidence,, to the airstream in a wind tunnel. The results for=0 are compared with Pohlhausen's formula and the results for 0, with Powell's experiments. The local values of Sh on the upwind and downwind sides of discs have been obtained. For leaves in the canopy, Sh was found to be about 25% greater than would be predicted by applying Pohlhausen's equation without correction for orientation.     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

The dynamical range of global circulations — II

The dynamical range of global atmospheric circulations is extended to specialized parameter regions by evaluating the influence of the rotation rate () on axisymmetric, oblique, and diurnally heated moist models. In Part I, we derived the basic range of circulations by altering for moist and dry atmospheres with regular and modified surfaces. Again we find the circulations to be composed of only a few elementary forms. In axisymmetric atmospheres, the circulations consist of a single jet in the rotational midrange (^*=1/2–1) and of double jets in the high range (^*=2–4), together with one or two pairs of Hadley and Ferrel cells; where (^*=/ _E ) is the rotation rate normalized by the terrestrial value. These circulations differ from those predicted by firstorder symmetric-Hadley (SH₁) theory because the moist inviscid atmosphere allows a greater nonlinearity and prefers a higher-order meridional mode. The axisymmetric circulations do, however, resemble the mean flows of the natural system — but only in low latitudes, where they underlie the quasi-Hadley (QH) element of the MOIST flows. In midlatitudes, the axisymmetric jets are stronger than the natural jets but can be reduced to them by barotropic and baroclinic instablities. Oblique atmospheres with moderate to high tilts ( _P =25°–90°) have the equator-straddling Hadley cell and the four basic zonal winds predicted by the geometric theory for the solstitial-symmetric-Hadley (SSH) state: an easterly jet and a westerly tradewind in the summer hemisphere, and a westerly jet and an easterly tradewind in the winter hemisphere. The nonlinear baroclinic instability of the winter westerly produces a Ferrel cell and the same eddy fluxes as the quasi-geostrophic QG element, while the instability of the summer easterly jet produces a QG-Hadley (QGH) element with a unique, vertically bimodal eddy momentum flux. At high _P and low ^*, the oblique atmospheres reach a limiting state having global easterlies, a pole-to-pole Hadley cell, and a warm winter pole. At low tilts _P <10°, the oblique circulations have a mix of solstitial and equinoctial features. Diurnal heating variations exert a fundamental influence on the natural-Hadley (NH) circulations of slowly rotating systems, especially in the singular range where the zonal winds approach extinction. The diurnality just modifies the NH element in the upper singular range (1/45^*1/16), but completely transforms it into a subsolar-antisolar Halley circulation in the lower singular range (0^*<1/45). In the modified NH flows, the diurnality acts through the convection to enhance the generation of the momentum-transferring planetary waves and, thereby, changes the narrow polar jets of the nondiurnal states into broad, super-rotating currents. Circulation theory for these specialized flows remains rudimentary. It does not explain fully how the double jets and the multiple cells arise in the axisymmetric atmospheres, how the QGH element forms in the oblique atmospheres, or how waves propagate in the slowly rotating diurnal atmospheres. But eventually all theories could, in principle, be compared against planetary observation: with Mars testing the QGH elements; Jupiter, the high-range elements; Titan, the equinoctial and solstitital axisymmetric states; and Venus, the diurnally modified NH flows.     

Estimation of the Monin-Obukhov similarity functions from a spectral model

From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function _m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - C_i constants i=u, v, w, (see Equation (5) - k_me,k_mT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - k_mi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - k_sb, k_c characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - _m nondimensional mean shear - _H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - _D nondimensional divergence of vertical transports of turbulent linetic energy     

On Modelling the One-Dimensional Atmospheric Boundary Layer

Several commonly used turbulence closure schemes forthe atmospheric boundarylayer (ABL) are applied to simulate neutral, nocturnal and diurnal cyclesituations in a one-dimensional ABL. Results obtained with the differentschemes, E-, E- and its modified versions, and twoversions ofthe q ² Level 2.5, are compared and discussed.     

A note on free convection

The limiting condition of µ_*, 0 in the unstable atmospheric boundary layer is usually referred to as free convection. Some of the similarity laws that are proposed for this condition do not agree with experiment. A mechanism is proposed in this paper to show why the asymptotic free convection condition is never completely reached near the surface. It is shown that local turbulent shear production plays an important role in shaping the temperature profile even when the average velocity is zero.Contribution No. 267, Dept. of Atmospheric Sciences, University of Washington.     

Diurnal variation of horizontal wind direction fluctuations in complex terrain at Geysers,Cal.

Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km² Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s^-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.