Application of the Priestley-Taylor evaporation model to assess the influence of soil moisture on the evaporation from a large weighing lysimeter and class a pan

The influence of soil moisture on evaporation from a 6-m grass-covered lysimeter and from Class A pans was assessed for one summer using the -parameter of the Priestley-Taylor evaporation model appropriate for the individual surfaces computed on a daily basis. Net radiation over the pan was estimated from above-grass measurements using a correlation established between the two, using measurements made in the previous two summers. Changes in heat storage of the water were considered in the derivation of for the pan. A unique relationship for the particular conditions of the site was determined between the for the lysimeter and soil moisture, approaching 1.29 at soil moisture near field capacity, but decreasing to as low as 0.5 for dry soil. The corresponding relationship for the pan showed more scatter, but this was improved by using 5-day running means of evaporation and stratifying the data in terms of wind speed to yield a family of curves. Values for at wet soil conditions varied from 1.07 for 100 km day^–1 wind run to 1.17 for 250 km day^–1 wind run. For each curve, values of increased by about 20%; as the soil dried. The relationships may be used to reduce observed Class A pan evaporation to equivalent values for wet-soil conditions and to estimate near-surface soil moisture and actual evapotranspiration for this particular site. Extension of the technique to other areas requires derivation of similar relationships appropriate for those other locations     

Interaction between the atmospheric and oceanic boundary layers

The two-layer system of an atmosphere over water bodies is reduced to a single-layer problem. Values of the interfacial quantities, such as the friction velocity, the surface velocity, the angles, and , between the surface shear stress and the geostrophic wind velocity and the surface wind velocity, respectively, and the surface roughness, all of which depend upon external parameters, such as the geostrophic wind and stratifications, are obtained. The geostrophic drag coefficient C _d, the geostrophic wind coefficient C _f, and the angles , and , of the turbulent flow at the sea-air interface are functions of a dimensionless number, mfG/kg, with S ₁ and S ₂ as two free stratification parameters. The surface roughness is uniquely determined from the geostrophic wind rather than from the wind profile in the boundary layer.Formerly Visiting Research Associate, Applied Physics Branch, Earth Observations Division, NASA-Manned Spacecraft Center, Houston, Texas.     

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.     

The effect of synoptic-scale advection on the performance of the Priestley-Taylor evaporation formula

Measurements of the appropriate parameters for the calculation of the latent heat flux over a black spruce forest in northern Quebec were carried out in August, 1980. Values of the Priestley-Taylor parameter, , were derived by exploiting the Bowen-ratio-energy-balance (BREB) technique. Derived values of are then related to synoptic-scale warm air advection, derived from surface synoptic charts and tephigrams of the planetary boundary layer. It is found that when warm air advection is present, > 1.26, especially when the surface is wet. When advection enhancement is removed, however, values of approach unity. A new approach to calculating the latent heat flux, when warm air advection is present, is therefore proposed.     

The behaviour of aerosol optical depth during clear summer days in Athens,Greece

Summary The integral aerosol optical depths (k ) at the hour of 08:20 Local Standard Time (LST), are compared with those calculated previously at 11:20 and 14:20 LST, for clear days during summer in Athens over the period 1962–1988. The mean values at 08:20 LST were consistently lower than the values at 11:20 and 14:20 LST. The influence of the vertical wind profile on the values ofk was also investigated. A comparison was made of the wind profiles at 02:00 and 14:00 LST, for days in which the 11:20 and 14:20 LST values ofk were 0.200 andk 0.350, respectively. The corresponding bulk wind shear s was also found for the period 1980–1988. The most significant results occurred with the first category of days. The resultant wind velocities from the surface to the 900 hPa level, in each hour were higher by 2–4 m·s^–1 with respect to the corresponding values for the second category. At 02:00 LST the bulk wind shear showed a considerable difference (1.8) between the two categories of days in the surface to 700 hPa layer at 02:00 LST. Finally, the associated weather conditions that appear to initiate a period of low values ofk (k 0.200) at 11:20 and 14:20 LST were examined for the period 1980–1988. Fifteen such cases were identified and it was found that they all occurred after the passage of weak cold fronts.With 6 Figures     

Evaporation from a blanket bog in a foggy coastal environment

Frequent fog severely restricts evaporation from blanket bogs in Newfoundland because it more than halves the radiant energy input, and it eliminates the vapor pressure deficit, resulting in evaporation at the equilibrium rate (average = 0.99 during fog). During these periods, there is no surface resistance to evaporation because the bog has been wetted by fog drip, and although the latent heat flux dominates over sensible heat (average = 0.8), both are small. In contrast, the surface dries during clear periods, increasing the surface resistance to evaporation so that sensible heat becomes more important ( = 1.05). When the mosses are dry, evaporation is below the equilibrium rate ( = 0.87), although the higher available energy ensures that actual evaporation is higher. During clear periods, daily evaporation averaged 2.5 mm, compared to 1.1 and 0.7 mm for fog and rain, respectively. The suppressed evaporation at this site is important in maintaining appropriate hydrological conditions for blanket bog development.     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

Comments on ’a Numerical Simulation of Boundary-Layer Flows Near Shelterbelts‘ by H. Wang and E. Takle

We implemented the windbreak model of Wang and Takle to investigate whythose authors could have obtained better agreement with the Bradley andMulhearn wind reduction measurements in the far lee of a fence, than wasobtained by Wilson using a very similar model. According to our experiencethe improvement of the Wang-Takle simulations (relative to Wilson's)largely arises from their having used a too-shallow computational domain(8H, versus Wilson's 47H; H being the windbreak height).     

A generalization of the Kolmogorov-von Kármán relationship and some further implications on the values of the constants

A general relation between the Kolmogorov (₁) and von Kármán (k) constants appropriate to flows in the surface boundary layer is obtained by including stability and divergence dependent terms in the development. The result is used to show that the author's earlier inference of the valuek = 0.36 drawn from the observed value ,k ^4/3=0.141 cannot have been underestimated due to the neglect of stability or divergence effects. It is concluded that if k is to be 0.4, then ₁ must be 0.5.Work performed under the auspices of the U.S. Energy Research and Development Administration.     

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.     

A note on the analogy between momentum transfer across a rough solid surface and the air-sea interface

The aerodynamic classification of the resistance laws above solid surfaces is based on the use of a so-called Reynolds roughness number Re _s =h _s u _*/, whereh _s is the effective roughness height, -viscosity,u _*-friction velocity. The recent experimental studies reported by Toba and Ebuchi (1991), demonstrated that the observed variability of the sea roughness cannot be explained only on the basis of the classification of aerodynamic conditions of the sea surface proposed by Kitaigorodskii and Volkov (1965) and Kitaigorodskii (1968) even though the latter approach gains some support from recent experimental studies (see for example Geernaertet al. 1986). In this paper, an attempt is made to explain some of the recently observed features of the variability of surface roughness (Toba and Ebuchi, 1991; Donelanet al., 1993). The fluctuating regime of the sea surface roughness is also described. It is shown that the contribution from the dissipation subrange to the variability of the sea surface can be very important and by itself can explain Charnock's (1955) regime.     

Shear stress results from a mixed spectral finite-difference model: Application to the Askervein Hill project data

Beljaars et al. (1987) developed a model for neutrally stratified boundary-layer flow over roughness changes and topography. It has been discovered that a constant parameter, , was missing in the algebraic-stress closure equations of their paper. This omission also occurred in the MSFD model code but only affects the Askervein Hill shear-stress results for the E-- turbulence closure in Beljaars et al. It also changes the stress results of Karpik (1988), but not his conclusions regarding the robustness of his improved numerical scheme. The present paper demonstrates the effect of the omission of the parameter, , and tests the sensitivity of the model to variations in its value. The new results are compared with the data and model results of Zeman and Jensen (1987).     

The behaviour of a mixed-layer model of the convective boundary layer coupled to a big leaf model of surface energy partitioning

Daily mean values of the Priestley-Taylor coefficient, ¯, are derived from a simple model of the daily growth of a convective boundary layer. For a particular control set of driving environmental variables, ¯ is related to the prescribed bulk surface resistance, r_S by 1/¯ = 1/₀ + mr_S for parameters ₀ and m. The dependence of the parameters ₀ and m on weather is explored and a potential use of this linear relation to provide information about regional values of r_S is indicated.     

A THEORETICAL BASIS FOR THE PRIESTLEY-TAYLOR COEFFICIENT

The relationship between potential evaporation and arealevaporation is assessed using a closed-box model of the convectiveboundary layer (CBL). Potential evaporation is defined as theevaporation that would occur from a hypothetical saturated surface,with radiative properties similar to those of the whole area, and smallenough that the excess moisture flux does not modify thecharacteristics of the CBL. It is shown that the equilibrium rate ofpotential evaporation is given by Ep0=E0,where E0 is the equilibrium evaporation (radiative termof the Penman formula), and is a coefficient similar to thePriestley-Taylor coefficient. Its expression is , where is the areal surface resistance, r_a is the localaerodynamic resistance, and is the dimensionless slope of thesaturation specific humidity at the temperature of the air. Itscalculated value is around 1 for any saturated surface surrounded bywater, about 1.3 for saturated grass surrounded by well-watered grassand can be greater than 3 over saturated forest surrounded by forest.The formulation obtained provides a theoretical basis to the overallmean value of 1.26, empirically found by Priestley and Taylor for thecoefficient . Examining, at the light of this formulation, thecomplementary relationship between potential and actual evaporation(as proposed by Bouchet and Morton), it appears that the sum ofthese two magnitudes is not a constant at equilibrium, but depends onthe value of the areal surface resistance.     

Monoterpene concentrations in and above a forest of scots pine

Diurnal and vertical ambient air measurements of the monoterpenes have been made in and above a Scots pine (Pinus sylvestris) forest of central Sweden, within the boreal northern coniferous biome. Sampling was done with Tenax TA, and analysis by GC and ion trap detection. Daytime mixing ratios were on the order of tenths of a ppbv from the forest floor to the top of the forest, and a factor of 2 or 3 lower above the forest. Mixing ratios at night were at the ppbv level, highest near the forest floor and the crown, and decreased with height above the forest. The highest total concentration observed was 8 ppbv inside the forest at 3 am (GMT). The average terpene composition was 3-carene 32%, -pinene 29%, limonene 18%, -pinene 10%, -phellandrene 7%, camphene 5%, and sabinene at less than 2%. The 3-carene/-pinene ratio varied with wind direction and speed, relative humidity, and wet/dry vegetation, but not with ozone or NO₂ concentration, solar radiation, or temperature. Variations in the observed terpene composition at the sampling site are mainly caused by the influence of other vegetation in the vicinity of the site. It would seem that wet Scots pine emits more 3-carene relative to -pinene than does dry pine.     

Analytical investigations of horizontal meso-scale momentum equations with Newtonian nudging

Summary This study presents an analytical investigation of the local behaviour of the solution to a mesoscale model with Newtonian nudging when observed winds are time varying. The analysis examines each Fourier component of the time series of observed winds. Unlike the case with a constant observed wind, the nudged wind vector does not asymptotically approach the observed wind. In response to sinusoidal oscillation of the observed wind, the nudged wind vector is always on a half circle connecting the vector ends of the observed and un-nudged modelled winds. When nudging parameter 0, the nudged wind vector approaches the un-nudged wind; when , the nudged wind vector approaches the observed wind. For commonly used values of nudging parameter , the modelled wind field always carries errors.A target nudging scheme is devised in this study in order to ensure the model result is identical to observed winds with sinusoidal oscillation. Investigation shows that such a target wind exists for a finite value of , and the magnitude of the target-nudging term is about the same as that of a normal nudging term if f, wheref is the Coriolis parameter and is the frequency of the wind oscillation.With 7 Figures     

The effects of orography on precipitation

In this paper, a model simulating the effects of topography and altitude on precipitation is presented. Topography has its maximum effect on precipitation when the angle which the wind makes with the slope direction approaches zero and the inclination of the slope is near 45°. The smaller the angle , the greater the influence of slope on precipitation. When <45°, the larger the inclination, the greater the influence of slope on precipitation and the less the difference in precipitation between the windward and the leeward slopes. When <45°, the reverse holds. But for in the range of 0°–45° and in the range 30°–60°, differences in precipitation on both the windward and leeward slopes are not so well marked and can be neglected in general. In condition of uniform slope inclination, precipitation on the windward slope increases with altitude at first and then decreases after attaining a height (H _m ) of maximum precipitation; alsoH _m is greater, the drier the air mass. When the terrain on the windward side is stepped in shape, it is possible that more than one height of maximum precipitation will occur.     

Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation

The use of analytical solutions of the diffusion equation for footprint prediction is explored. Quantitative information about the footprint, i.e., the upwind area most likely to affect a downwind flux measurement at a given height z, is essential when flux measurements from different platforms, particularly airborne ones, are compared. Analytical predictions are evaluated against numerical Lagrangian trajectory simulations which are detailed in a companion paper (Leclerc and Thurtell, 1990). For neutral stability, the structurally simple solutions proposed by Gash (1986) are shown to be capable of satisfactory approximation to numerical simulations over a wide range of heights, zero displacements and roughness lengths. Until more sophisticated practical solutions become available, it is suggested that apparent limitations in the validity of some assumptions underlying the Gash solutions for the case of very large surface roughness (forests) and tentative application of the solutions to cases of small thermal instability be dealt with by semi-empirical adjustment of the ratio of horizontal wind to friction velocity. An upper limit of validity of these solutions for z has yet to be established.     

The use of σ T as a temperature scale in the atmospheric surface layer

The standard deviation of temperature _T is proposed as a temperature scale and as a velocity scale to describe the behaviour of turbulent flows in the Atmospheric Surface Layer (ASL), instead of _* andu _* of the Monin—Obukhov similarity theory, and ofT _f andU _f used for free convection stability conditions. On the basis of experimental evidence reported in the literature, it is shown that _T T _f andv _* U _f in the free convection region, and _{T *} andv _* U _* in nearneutral and stable conditions. This implies that the proposed scales can be applied for all stabilities. Furthermore, a new length scale is proposed and its relation with Obukhov length is given. Also, a simple semi-empirical expression is presented with which _T andv _* can be evaluated in a rather simple way. Some examples of practical applications are given, e.g., a stability classification for unstable conditions.     

On the enhanced smoothing over topography in some mesometeorological models

An equation is derived for the components of the horizontal (turbulent) frictional force in the -coordinate system with special attention to mesometeorological flow models. The starting point is the horizontal equation of motion in its flux-form in the -system in which we replace (following Reynolds' procedure) the velocity components u,v and % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbai % aaaaa!37B8! \[ \dot \sigma \] aswell as other relevant quantities by terms of the form u = + u,..., = ± + % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbai % Gbauaaaaa!37C3! \[ \dot \sigma ' \] , etc. ( = time average of u; u = fluctuating part of u.) Next, the equation is averaged with respect to time and terms which we believe are small in mesometeorological flows, are neglected. On expressing by an appropriate expression that involves w, the result shows the appearance of two new terms which, have not been considered previously in the published literature. While the expression earlier used in the literature involved the -derivative of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG3bGbauaaaaaaaa!380B!\[\overline {u'w'} \] alone, the new terms add the -derivatives of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaadaahaaWcbeqaaiaaikdaaaaaaaaa!37EC!\[\overline {u'^2 } \] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG2bGbauaaaaaaaa!380A!\[\overline {u'v'} \] for the x-component of the force, and the -derivatives of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG2bGbauaadaahaaWcbeqaaiaaikdaaaaaaaaa!37ED!\[\overline {v'^2 } \]} and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG2bGbauaaaaaaaa!380A!\[\overline {u'v'} \] for the y-component, where and are the slopes of the -surfaces in the x- and y-directions, respectively. Further, a few numerical simulations of the sea-breeze over topography are carried out with and without the correction terms. It is shown that when corrections terms are not included the effective smoothing is stronger above the sloping regions and may amount to as high as 50 percent of the convergence with slopes of ~.04. The ìnclusìon of the new terms does not lead to any special computational difficulties and for that reason there is no compelling reason to neglect them, all the more so because, as is shown, the addition of the new terms results in a consistent apportioning of the degree of horizontal diffusion.On leave from CIMMS, Norman, OK.Now visiting Dept. of Met., Helsinki, Finland.