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.     

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.     

Radiative Heat Transfer and Hydrostatic Stability in Nocturnal Fog

We have performed a one-dimensional and transient radiative heat transfer analysis in order to investigate interaction between atmospheric radiation and convective instability within a nocturnal fog. The radiation element method using the Ray Emission Model (REM²), which is a generalized numerical method, in conjunction with a line-by-line (LBL) method, is employed to attain high spectral resolution calculations for anisotropically scattering fog. The results show that the convective instability has a strong dependence on radiative properties of the fog. For the condition of a 20-m droplet diameter and liquid water content of 0.1 × 10^–3 kg m^–3;, the temperature profile within the fog becomes S shaped, and a convective instability layer forms in the middle or lower level of the fog. However, for the same water content and a 40-m diameter droplet, no strong convective instability layer forms, whereas for a 10-m diameter droplet a strong convective instability is observed.     

A study of a fair weather boundary layer in TOGA-COARE: Parameterization of surface fluxes in large scale and regional models for light wind conditions

The study focuses on a way to parameterize the effect of subgrid scale convective motions on surface fluxes in large scale and regional models for the case of light surface winds. As previously proposed, these subgrid effects are assumed to scale with the convection intensity through the relationship: where is the mean velocity of the wind, U₀ the velocity of the mean wind, w_* the free convection velocity, and an empirical coefficient to be determined. Both observations and numerical simulation are presently used to determine the free convection coefficient .Large eddy simulation of a fair weather convective boundary layer case observed during TOGA-COARE is performed. Comparisons between observations and the simulation of surface properties and vertical profiles in the planetary boundary layer are presented. The simulated vertical turbulent fluxes of heat, moisture and buoyancy range well within estimates from aircraft measurements.The most important result is that the true free convection coefficient , directly estimated from simulation, leads to a value of 0.65, smaller than the ones estimated from temporal and spatial variances. Using observations and simulation, estimates of from temporal and spatial variances are obtained with similar values 0.8. From both theoretical derivations and numerical computations, it is shown that estimates of the true from variances are possible but only after applying a correction factor equal to 0.8. If this correction is not used, is overestimated by about 25%. The time and space sampling problem is also addressed in using numerical simulations.     

Plötzlicher Herztod und Wetter

Zusammenfassung An einem Material von 100 Fällen eindeutig gesicherten Herztodes, die im Institut für Gerichtliche und Soziale Medizin in Frankfurt am Main aus einem bestimmten Bereich — Raum Frankfurt am Main — und einer bestimmten zeit—Kalenderjahre 1948 bis 1950—zur Autopsie gelangten, einerseits und den bioklimatischen Daten für den gleichen Raum und die gleiche Zeit anderseits, wurde der Zusammenhang von Herztodesfällen und Wettervorgängen mit der Wahrscheinlichkeitsüberlegung, dem Wahrscheinlichkeitsintegral und formalstatistischen Methoden überprüft. Die Auswertung ergab, daß eine starke Abhängigkeit des Herztodes vom Wettergeschehen besteht. Die Häufung der Herztodesfälle konnte für die Sammelgruppe Wetterstörungen, für labiles Aufgleiten, Warmfront mit Aufgleiten und Kaltfront mit Turbulenz ebenso eindeutig gesichert werden, wie die Verminderung der Herztodesfälle in störungsfreien Zeiten. Die hohe Wahrscheinlichkeit einer Beziehung besteht ferner für Kaltfront mit Aufgleiten und in geringerem Maße auch für Föhndurchbruch. Während bei Okklusion das Material zu einem abschließenden Urteil zu gering war, konnte bei allen übrigen Wetterlagen eine Beziehung zum plötzlichen Herztod nicht nachgewiesen werden. Auch für die Gruppen, bei denen eine enge Beziehung statistisch gesichert werden konnte, sei betont, daß die Wettervorgänge nicht die alleinige Ursache des plötzlichen Herztodes, sondern nur ein auslösendes Moment bei schon vorher schwer geschädigten Herzen darstellen.

---

Summary On the basis of 100 ascertained cases of cardiac death from a certain region—province Frankfort on-the-Main—and a certain period—calendar years 1948 to 1950—, submitted to autopsy in the Institute for Legal and Social Medicine at Frankfort on-the-Main, and of bioclimatic data from the same region and period, the relation between cardiac death and meteorological phenomena was examined by means of probability considerations, the probability integral, and methods of formal statistics. Evaluation of this material showed cardiac death to be strongly dependent on weather conditions. Increase of frequency of cardiac death could be established without any doubt for the collective groups weather disturbances, unstable ascending motion, warm front with ascending motion, cold front with turbulence, whereas decrease of frequency was found in undisturbed periods. Furthermore, high probability of an existing relation was found for cold front with ascending motion and, to a lesser degree, also for foehn invasion. As to occlusion no definite conclusions could be drawn for scarcity of observations. For all the other weather types no relation to sudden cardiac death could be proved. Meteorological phenomena are by no means the only cause of sudden cardiac death, not even for the groups in which a close relation could be established statistically. They have but an initiatory effect in case of prior severe heart diseases.

---

Résumé Disposant d'un matériel médical de 100 cas de mort cardiaque qui ont fait l'objet d'autopsies à l'Institut de médecine légale et sociale de Francfort sur le Mein (région de Francfort et période de 1948 à 1950), et disposant d'autre part pour la même région et pour la même période d'observations bioclimatologiques, l'auteur a étudié la relation de ces décès avec les phénomènes météorologiques à l'aide des méthodes modernes de statistique, de l'intégrale de probabilité etc. Le résultat montre une étroite dépendance centre les cas de mort cardiaque et les conditions météorologiques. On a pu vérifier sans doute possible l'augmentation des décès lors de «perturbations», d'«ascendances labiles», de «fronts chauds avec glissement ascendant», de «fronts froids avec turbulence», de même que leur diminution lors de périodes non troublées. La forte corrélation subsiste avec les «fronts froids avec ascendance» et à un moindre degré avec «l'invasion du foehn». Tandis que pour les «occlusions» le matériel disponible était trop restreint pour aboutir à une conclusion, l'examen des situations autres que les précédentes n'a pas conduit à une relation intime. Il convient de remarquer que pour tous les groupes de météores ayant donné une corrélation statistique très nette, les phénomènes météorologiques n'ont joué qu'un rôle déclancheur chez des malades déjà gravement atteints.

---

---

Mit 2 Textabbildungen.     

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.     

Interactions between internal gravity waves and their traumatic effect on a continuous stratification

An investigation is made of the detailed kinematic conditions which characterize breaking of internal waves in a continuous density field. It has been observed that breaking is associated with localized, abruptly appearing intensifications of the density gradient. Turbulence subsequently develops from these intensifications or traumata.Interactions between two or more disparate waves produce distortions more extreme than would be deduced from linear combination. This is exploited experimentally by crossing two weak internal wave beams from separate sources.Lagrangian integration of the motion field defined by the beams gives extrema of isopycnic slope, Richardson number defined by shear across isopycnics and dynamic vertical density gradient, and dynamic density gradient, under conditions for the appearance of traumata.     

Thermal comfort of man in different urban environments

Summary On July 29, 1985, a hot summer day, biometeorological measurements were performed simultaneously in three different urban structures within the city of Munich and in the trunk space of a nearby tall spruce forest. Based on the results of these experiments the following thermophysiologically relevant biometeorological indices were calculated: Predicted mean vote, skin wettedness and physiologically equivalent temperature. These three indices are derived from different models for the human energy balance. They allow the assessment of the thermal components of the microclimates at the selected sites with regard to application in urban planning. The results quantitatively show the great heat stress in the urban structure street canyon, exposed to south, whereas in the trunk space of the tall spruce forest there is nearly an optimal climate even on hot summer days. Between these extremes the results for street canyon, exposed to north show a little higher heat load than for backyard with trees.

---

Zusammenfassung An einem heißen Sommertag, dem 29. Juli 1985, wurden in drei Stadtstrukturen in München und im Stammraum eines nahegelegenen Fichtenhochwaldes zeitgleich biometeorologische Messungen durchgeführt. Mit den Meßergebnissen wurden folgende thermophysiologisch relevante biometeorologische Indizes berechnet: Predicted mean vote, Hautbenetzungsgrad und physiologisch äquivalente Temperatur. Diese drei Indizes beruhen auf verschiedenen Modellen zur menschlichen Energiebilanz. Mit den drei Indizes wurden die thermischen Komponenten der Mikroklimate an den ausgewählten Meßplätzen im Hinblick auf Stadtplanungsaufgaben bewertet. Die Ergebnisse zeigen quantitativ die relativ große Hitzebelastung bei der Stadtstruktur Straßenschlucht, nach Süd exponiert, während im Stammraum des Fichtenhochwaldes selbst an heißen Sommertagen nahezu optimale Bedingungen herrschen. Zwischen diesen Extremen liegen die Ergebnisse für die anderen Meßplätze, wobei für Straßenschlucht, nach Nord exponiert die Wärmebelastung etwas höher als für Innenhof mit Bäumen ist.

---

---

With 6 Figures     

Carbon dioxide emissions in a methane economy

Increasing reliance on natural gas (methane) to meet global energy demands holds implications for atmospheric CO₂ concentrations. Analysis of these implications is presented, based on a logistic substitution model viewing energy technologies like biological species invading an econiche and substituting in case of superiority for existing species. This model suggests gas will become the dominant energy source and remain so for 50 years, peaking near 70 percent of world supply. Two scenarios of energy demand are explored, one holding per capita consumption at current levels, the second raising the global average in the year 2100 to the current U.S. level. In the first (efficiency) scenario concentrations peak about 450 ppm, while in the second (long wave) they near 600 ppm. Although projected CO₂ concentrations in a methane economy are low in relation to other scenarios, the projections confirm that global climate warming is likely to be a major planetary concern throughout the twenty-first century. A second finding is that data on past growth of world per capita energy consumption group neatly into two pulses consistent with longwave theories in economics.     

Stable Isotope Ratios: Hurricane Olivia

The oxygen and hydrogen isotopic compositions of rains from HurricaneOlivia (1994) in the eastern Pacific were measured. The rains werecollected on 24 and 25 September during airplane flights conducted at anelevation of 3 km. Hurricane Olivia peaked in intensity to a category-4storm between the two dates. Isotope ratios of rains from HurricaneOlivia were markedly lower ( ¹⁸O = –13.9to –28.8) than that of rain collected from a thunderstormat an elevation of 2.3 km outside the influence of Olivia (¹⁸O = –3.8). A distinct decrease in isotoperatios from the first day to the next ( ¹⁸O =–18.4 to –21.9) in Hurricane Olivia wasattributed to decreased updraft velocities and outflow aloft. Thisshifted the isotopic water mass balance so that fewer hydrometeors werelifted and more ice descended to flight level. A decrease in the averagedeuterium excess from the first day to the next (d = 15.5 to 7.1)was attributed to an increase in the relative humidity of the watervapor `source' area. We hypothesize that the `source' region for therain was in the boundary layer near the storm center and that becausethe hurricane was at peak intensity prior to the second day the relative humidity was higher.     

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.     

Fog chemistry during EUMAC Joint Cases: Analysis of routine measurements in southern Germany and model calculations

Summary The role of fog events for acid production and moist deposition in southern Germany during three EUMAC Joint Cases has been investigated by analysis of routine measurements and simulations with a one-dimensional fog-chemistry model. To identify the chemical and meteorological conditions, routine measurements by different institutions have been interpreted. The periods under consideration include a smog episode with low photooxidant concentrations during the Winter Case in February 1982. The Wet Case in spring 1986 represents a period with higher photooxidant concentrations. Conditions which are mostly characterized by low SO₂ and oxidant concentrations and comparatively high pH-values in the fog are given during the SANA 1 case in autumn 1990. Fog mostly occurs as a subscale phenomenon, but sometimes it can also cover large areas and it can contribute significantly to moist deposition.The model results indicate that the liquid phase sulfate production in the fog layer may even exceed the gas phase production during 24 hours within a layer of the same height occasionally. On the other hand, during the SANA 1 case the sulfate production in the fog was extremely low at night due to lack of oxidants and SO₂. Depending on the dissipation time of the fog a remarkable effect on the photolysis rates is possible. Since a significant amount of particulate mass is lost by moist deposition during fog, it is evident that fog events can have a noticeable effect on some of the gas phase constituents which are easily soluble.With 17 Figures     

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     

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     

Northern hemisphere atmospheric response to changes of atlantic ocean SST on decadal time scales: a GCM experiment

Analyses indicate that the Atlantic Ocean seasurface temperature (SST) was considerably colder at the beginning than in the middle of the century. In parallel, a systematic change in the North Atlantic sea-level pressure (SLP) pattern was observed. To find out whether the SST and SLP changes analyzed are consistent, which would indicate that the SST change was real and not an instrumental artifact, a response experiment with a low-resolution (T21) atmospheric GCM was performed. Two perpetual January simulations were conducted, which differ solely in the Atlantic Ocean (40° S-60° N) SST: the cold simulation utilizes the SSTs for the period 1904–1913; the warm simulation uses the SSTs for the period 1951–1960. Also, a control run with the model's standard SST somewhat between the cold and warm SST was made. For the response analysis, a rigorous statistical approach was taken. First, the null hypothesis of identical horizontal distributions was subjected to a multivariate significance test. Second, the level of recurrence was estimated. The multivariate statistical approaches are based on hierarchies of test models. We examined three different hierarchies: a scale-dependent hierarchy based on spherical harmonics (S), and two physically motivated ones, one based on the barotropic normal modes of the mean 300 hPa flow (B) and one based on the eigenmodes of the advection diffusion operator at 1000 hPa (A). The intercomparison of the cold and warm experiments indicates a signal in the geostrophic stream function that in the S-hierarchy is significantly nonzero and highly recurrent. In the A-hierarchy, the low level temperature field is identified as being significantly and recurrently affected by the altered SST distribution. The SLP signal is reasonably similar to the SLP change observed. Unexpectedly, the upper level stream-function signal does not appear to be significantly nonzero in the B-hierarchy. If, however, the pairs of experiments warm versus control and cold versus control are examined in the B-hierarchy, a highly significant and recurrent signal emerges. We conclude that the cold versus warm response is not a small disturbance that would allow the signal to be described by eigenmodes of the linear system. An analysis of the three-dimensional structure of the signal leads to the hypothesis that two different mechanisms are acting to modify the model's mean state. At low levels, local heating and advection are dominant, but at upper levels the extratropical signal is a remote responce to modifications of the tropical convection.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dilmenil.AWI Publication no. 254     

Measurements of solar radiation and estimation of optical depth in the high Andes of Peru

Summary During an expedition to the high Andes of Southern Peru in June–July 1977, measurements of direct solar radiation in four spectral bands (0.270–0.530–0.630–0.695–2.900 ) were conducted at six sites in elevations ranging from sea level to 5645 m. These measurements were evaluated in Langley plots to determine total optical depths () and irradiances at the top of the atmosphere. In addition, water vapor optical depths (_wv) were calculated from the mean radiosounding over Lima during the expedition, and Rayleigh (_ray) and ozone (_oz) optical depths were obtained from published tabulations. Subtracting _ray, _oz, and _wv from yielded estimates of aerosol optical depth _aer. The components _ray and _oz decrease from the shorter towards the longer wavelength bands and from the lower towards the higher elevation sites; _aer also decreases towards the higher elevations. Particularly pronounced is the decrease of _aer and from the lowlands of the Pacific coast to the highlands of the interior, reflecting the effect of a persistent lower-tropospheric inversion and the contrast from the marine boundary layer to the clear atmosphere of the high Andes.With 4 Figures     

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     

An alternative analysis of flux-gradient relationships at the 1976 ITCE

An extensive micrometeorological data set from the 1976 International Turbulence Comparison Experiment (ITCE) is analysed to determine flux-gradient relationships in an unstable atmosphere for momentum, sensible heat and water vapour transfers. The data are first analysed for internal consistency, resulting in the rejection of some data. Following a least-square fit to the remaining data in the form /k = (1 – z/L)^-/k, rounded-off values of k, , and are selected for each form of transfer consistent with the statistical accuracy of the measurements. The equations finally adopted are _M = (1 – 28z/L)^-1/4 and _{H, W} = (1 – 14z/L)^-1/2 with k _M = k_H = k_W = 0.40.These expressions fit the averaged observations to within a few per cent in the stability range of the experiment (-4 < z/L < -0.004).     

Stable boundary layers: Observations,models and variability part I: Modelling and measurements

Selected field measurements of evening stable boundary layers are presented in detail comparable with published Large Eddy Simulation results. Such models appear to match idealized theories more closely than do some boundary-layer observations. Any attempt to compare detailed observations with idealized models therefore highlights the variability of the real boundary layer.Here direct turbulence measurements across the stable boundary layer from a heterogeneous and an ideal site are contrasted. Recommendations are made for the information needed to distinguish heterogeneous and ideal cases.The companion paper (Part II) discusses further the issues of data, analysis in the presence of variability, and the effects of averaging over heterogeneous terrain.Part of UK Meteorological Office Atmospheric Process Research Division.     

A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.