Uncoupled multi-layer model for the transfer of sensible and latent heat flux densities from vegetation

The sensible heat flux density C and the latent heat flux density E are coupled in the case of a multi-layer model of vegetation. Therefore two linearly independent combinations of C and E, the enthalpy flux density H and the saturation heat flux density J, are introduced. Two electrical analogues, for H and J, are designed. They are equivalent to the resistance scheme for C and E, but uncoupled. Penman's formulas for C and E, which are applicable only to single-layer models, can be expressed equivalently in terms of H and J. This version of Penman's formulas can be extended easily to multi-layer canopies.     

Modelling radiation quantities and photolysis frequencies in the troposphere

STAR (System for Transfer of Atmospheric Radiation) was developed to calculate accurately and efficiently the irradiance, the actinic flux, and the radiance in the troposphere. Additionally a very efficient calculation scheme to computer photolysis frequencies for 21 different gases was evolved. STAR includes representative data bases for atmospheric constituents, especially aerosol particles. With this model package a sensitivity study of the influence of different parameter on photolysis frequencies in particular of O₃ to Singlet D oxygen atoms, of NO₂, and of HCHO was performed. The results show the quantitative effects of the influence of the solar zenith angle, the ozone concentration and vertical profile, the aerosol particles, the surface albedo, the temperature, the pressure, the concentration of NO₂, and different types of clouds on the photolysis frequencies.Notation I _A(, ) actinic flux - I _H(, ) irradiance - L(, , , ) radiance - wavelength - azimuth angle - cosine of zenith angle - _s cosine of solar zenith angle - optical depth - _s scattering coefficient - _c extinction coefficient - _o single scattering albedo - p _mix mixed phase function - g _mix mixed asymmetry factor - J _gas photolysis frequency     

In-water and remote measurements of ocean color

Spectral measurements of downwelling irradiance, E _d(), above the surface, and of upwelling irradiance just below the surface, E _u(), allow computation of spectral values of the diffuse reflectance R() = E _u()/E _d(); this yields full information about the true color and brightness of the ocean. Typical results are presented and interpreted for waters very different in turbidity and phytoplankton content. Conversely, the possibility of infering the water content from R() data at selected wavelengths is examined in terms of the respective number of equations and unknowns. The necessary use of assumptions and of empirical laws is emphasized.The magnitude of the useful signal emerging from the water, and the magnitude of the additional signals due to specular reflexion at the interface and to atmospheric scattering, are compared on the basis of spectroradiometric measurements performed within and above the sea, from different altitudes. These unwanted signals are dominant, causing a drastic change in the spectral composition of the light received by a remote sensor. The evaluation of the atmospheric effect must be very precise in order to recover the marine signal with a sufficient accuracy for a meaningful application of any kind of algorithm.This work was supported in part by the Centre National d'Exploitation des Océans (under contract CNEXO 77/1695) and in part by the Centre National de la Recherche Scientifique (RCP 247 & ERA 278).     

Heat and momentum transport in an oak forest canopy

With the use of a sonic anemometer, vertical heat and momentum fluxes were measured at three different levels in an oak forest canopy. A quadrant analysis of the resulting data shows that approximately half of the transport occurs in extreme events lasting about 5 to 10% of the time. The partition of transport into momentum sweeps, bursts and interactions shows good agreement with existing data. The heat flux is analysed by observing the fluxes during the different momentum events and considering concurrent momentum and heat flux intensities by means of conditional probabilities. While low intensity (normal) events show similar probability distributions throughout the canopy, different structures appear at the three measurement heights for high intensity (extreme) events that can tentatively be explained by taking the temperature profile into account. This supports the idea that these events are coherent motions with scales comparable to the canopy height.     

A new theory of the energy spectrum

In a recent paper, the author introduced a new viscous boundary layer, called the mesolayer, in turbulent shear flow. Its importance stems from its location between the inner and outer regions which are controlled by the law of the wall and Reynolds number similarity, respectively. This intrusion prevents the classical overlap assumption which appears to be fundamental in the derivation of the classical logarithmic behavior. The mesolayer has a thickness proportional to Taylor's microscale . This, and the analogy between the energy equation for the spectrum function of isotropic turbulence and the momentum equation for shear flow, suggest the existence of a similar region in wavenumber space with wavenumber k ~ ^-1. This mesoregion separates the inner region k ~ k _s(where k _s^-1 and is the Kolmogorov length) and the outer region k k _e(where k _e ^-1 and l is the energy-containing eddy size) and again invalidates the overlap assumption which appears to be fundamental in the derivation of the classical k ^-5/3-behavior of the inertial subrange.Incorporation of the mesoregion into the argument leads to a new theory with k ^-5/3-behavior in two regions (^-1 k k _s) and (k _e k ^-1) although with two different coefficients of proportionality (Kolmogorov constants). This leads to a wandering of the spectrum curve about the classical k ^-5/3 line similar to a wandering in turbulent shear flow about the logarithmic curve. This is clearly indicated by the data for the variation of the Kolmogorov constant.Other data support the new theory. In particular, the location of the point k _mwhere the curve of the nonlinear energy-transfer function goes through zero shows agreement with the theory, i.e., k _m^-1.     

The Askervein Hill Project: Mean wind variations at fixed heights above ground

This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.     

Stability Effects on Heat and Moisture Fluxes at Sea

During the 1996 ASGAMAGE experiment we measured windspeed, air temperature T_a, watertemperature T_s, humidity and the momentum,heat and moisture fluxes at a research platform offthe Dutch coast. For each quantity we used several(sets of) instruments simultaneously. This allowed usto make an extensive assessment of the quality of themeasurements and to find optimal values for thevarious quantities for each run. From these values wecalculated C_H and C_E, theStanton and Dalton numbers, and reduced them to 10-mheight and neutral conditions. For this reductionwe made a separate analysis for the effect ofinclusion or non-inclusion of the assumption that theroughness length for heat or moisture is the same forthe neutral and non-neutral cases. Differences inthe reduced data due to this assumption turned out tobe well within the measurement error.For C_H we distinguished three separategroups of data: stable (A), unstable with_{a s} (B) and unstablewith thetas;_a > _s (C), with indicating the potential temperature.The stable data separate into two groups, depending onwater temperature and/or the wave field. The data ofgroup B showed a relation with wave age. The data ofgroup C consistently gave negative values forC_H, a result that might indicate conversion oflatent heat into sensible heat through condensation ofwater vapour just above the water surface. An attemptto re-analyse the data in terms of density fluxes,combining the effects of heat and moisture, still gavenegative transfer coefficients for group C.For the moisture flux we found the more conventionalresult of a separation in stable and unstable values;these categories showed a clear difference, but notrends with, for example, wind speed.We conclude that standard Monin-Obukhov similaritytheory cannot explain our data.     

On the correlation between turbulent velocity and temperature derivatives in the atmospheric surface layer

Velocity and temperature derivatives were obtained at a height of 4 m in the atmospheric surface layer above land. With the assumption of local isotropy, these measurements are used to obtain some statistics of the turbulent energy and scalar dissipation fields. These statistics include the variances of the logarithms of the scalar and velocity dissipation fields and the correlation between these logarithms. When used in conjunction with the hypotheses for fluctuations in turbulent dissipation rates of Obukhov and Kolmogorov, the statistics suggest that the dependence of the flatness factor of temperature derivative on the turbulent Rynolds number R is not as large as that which had been previously reported in the literature. The experimental data indicate a R ^0.5 dependence for the kurtosis of the temperature derivative and a R ^0.15 variation for the strain rate-scalar dissipation correlation.     

Turbulence measurements above a pine forest

Eddy fluxes of momentum, sensible and latent heat, and turbulence spectra measured over the Thetford Forest during 10 days in the Spring of 1973 are described. The measured total heat flux (H + E) for 122 20-min periods agreed closely on average with independent estimates from an energy balance method. There was evidence that the energy balance data gave small systematic overestimates of available energy during the hours before noon, compensated by slight underestimates for the remainder of the day. A comparison of measured wind speeds and friction velocities in neutral stability confirmed the validity of the aerodynamic method for estimating momentum fluxes at heights of a few roughness lengths above the canopy. In stable conditions the log-linear wind profileU = (u _*/k)(ln ((z -d)/z _o) + (z -d -z _o)/L) with = 3.4 ± 0.4 provided a good fit to the data. Spectra in unstable conditions were generally more sharply peaked than those measured by other workers over smoother terrain: differences were less marked in the case of vertical velocity in stable conditions. Temperature spectra in these stable conditions showed high energy at relatively low wavenumbers, andwT cospectra showed a cospectral gap; both of these results were associated with an intermittent sawtooth structure in the temperature fluctuations.Now at the Meteorological Office, Bracknell     

Numerical integration errors in calculated tropospheric photodissociation rate coefficients

Tropospheric photodissociation rate coefficients (J values) were calculated for NO₂, O₃, HNO₂, CH₂O, and CH₃CHO using high spectral resolution (0.1 mm wavelength increments), and compared to the J values obtained with numerically degraded resolution (=1, 2, 4, 6, 8, and 10 nm, and several commonly used nonuniform grids). Depending on the molecule, substantial errors can be introduced by the larger increments. Thus for =10 nm, errors are less than 1% for NO₂, less than 2% for HNO₂, +6.5% to -16% for CH₂O, -6.9% to +24% for CH₃CHO, and -24% to +110% for O₃. The errors for CH₂O arise from the fine structure of its absorption spectrum, and are prevalently negative (underestimate of J). The errors for O₃, and to a lesser extent for CH₃CHO, arise mainly from under-resolving the overlap of the molecular action spectrum and the tropospheric actinic flux in the wavelength region of stratospheric ozone attenuation. The sign of those errors depends on whether the actinic flux is averaged onto the grid before or after the radiative transfer calculation. In all cases studied, grids with 2 nm produced errors no larger than 5%.     

A comparison of bowen ratio and eddy correlation sensible and latent heat flux measurements above deciduous forest

Sensible and latent heat flux densities (H and E) were measured above a mature, 18 m deciduous forest during July and August, 1988, using the Bowen ratio-energy balance (BREB) and eddy correlation (EC) methods. EC estimates ofH and E underestimated day-time surface available energy by 11%. EC also partitioned available energy differently than BREB. for/L<0.0, EC favouredH and BREB favoured E. Practical and theoretical limitations of the BREB and EC methods above forests are discussed. The most plausible causes for the failure of EC to close the surface energy balance are a low frequency loss of flux and the failure of a single point measurement to account for the spatial dispersive flux. The most plausible causes of the EC-BREB energy partitioning anomaly are the invalidity of the BREB similarity assumption and the violation of flux-gradient diffusion assumptions in the near-field diffusion region.     

Numerical simulation of turbulent convective flow over wavy terrain

By means of a large-eddy simulation, the convective boundary layer is investigated for flows over wavy terrain. The lower surface varies sinusoidally in the downstream direction while remaining constant in the other. Several cases are considered with amplitude up to 0.15H and wavelength ofH to 8H, whereH is the mean fluid-layer height. At the lower surface, the vertical heat flux is prescribed to be constant and the momentum flux is determined locally from the Monin-Obukhov relationship with a roughness lengthz _o=10^–4 H. The mean wind is varied between zero and 5w _*, wherew _* is the convective velocity scale. After rather long times, the flow structure shows horizontal scales up to 4H, with a pattern similar to that over flat surfaces at corresponding shear friction. Weak mean wind destroys regular spatial structures induced by the surface undulation at zero mean wind. The surface heating suppresses mean-flow recirculation-regions even for steep surface waves. Short surface waves cause strong drag due to hydrostatic and dynamic pressure forces in addition to frictional drag. The pressure drag increases slowly with the mean velocity, and strongly with /H. The turbulence variances increase mainly in the lower half of the mixed layer forU/w _*>2.     

Field measurements of snow-drift threshold and mass fluxes,and related model simulations

Field measurements were carried out to calculate the threshold friction velocity for snow saltation, and mass fluxes during snow drift. The wind was measured in three components by an ultrasonic anemometer, and the mass fluxes were determined using an optical sensor (snow particle counter), acoustic sensors (Flowcapt) and mechanical traps. The threshold friction velocity was found to be correlated to the grain size (R²=0.75). The mass flux measurements were compared with numerical simulations of snow drift, and it was demonstrated that the maximum snow transport takes place at shear stress values of roughly two times the average shear stress over 20 min. By implementing a probability distribution for the shear stress the mass flux was simulated with only the mean measured value of the shear stress as input. This procedure enables the future use of the numerical model for operational applications.     

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.     

Growth Of The Summer Daytime Convective Boundary Layer At Anand

The heights of the daytime convective boundary layer (CBL), computed by a one-dimensional model for a bare soil surface at a semi-arid station,Anand, during the dry and hot summer month of May 1997, are presented. As input, the model requires surface heat flux, friction velocity and air temperature as functions of time. Temperature data at the one-metre level from a tower and sonic anemometer data at 9.5 m collected during the period 13–17 May 1997 in the Land Surface Processes Experiment (LASPEX-97) are used to compute hourly values of surface heat flux, friction velocity and Obukhov length following the operational method suggested by Holtslag and Van Ulden [J. Climate Appl. Meteorol. 22,517–529 (1983)]. The model has been tested with different values for the potential temperature gradient ( ) above the inversion. The model-estimated CBL heights comparefavourably with observed heights obtained from radiosonde ascents.     

A SIMPLE FORMULA FOR ATTENUATION OF EDDY FLUXES MEASURED WITH FIRST-ORDER-RESPONSE SCALAR SENSORS

A simple formula, (1 + (2fmc))-1,is proposed to estimate the attenuation of a scalar flux measurement made by eddy-correlation using a fast-response anemometer and a linear, first-order-response scalar sensor with a characteristic time constant c.In this formula, =7/8 for neutral and unstable stratification within the surface-flux layer and =1 both within the convective boundary layer (CBL) and for stable stratification in the surface layer.fm is the frequency of the peak of the logarithmic cospectrum and can be estimated from fm = nm /z, where z is the measurement height and is thewind speed at that height. The dimensionless frequency at the cospectral maximum nm is estimated here from observations of its behavioras a function of atmospheric stability, z/L within the surface layeror z/zi within the CBL, where L is the Obukhov stability length and zi is the depth of the CBL. The predicted dependence of flux attenuation on measurement height is discussed.     

Turbulence Reynolds number and the turbulent kinetic energy balance in the atmospheric surface layer

The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .     

Der Umkehreffekt von Götz Messungen des Zenithimmelslichtes mit dem Multiplier

Zusammenfassung Bisher wurde der Umkehreffekt vonGötz nur mittels Photozelle oder Zählrohr untersucht. Die Anwendung des Multipliers als Meßinstrument bietet wegen dessen wesentlich größerer Empfindlichkeit viele Vorteile. So kann damit noch bei Zenitdistanzen der Sonnez>90° bei kleineren Wellenlängen als bisher und bei engerem Spalt des Spektrometers (bessere Monochromasie) gemessen werden.Resultate: 1. Die Verschiebung der Minima der Umkehrkurven,i _n /i _m =f(z), erfolgtgleichsinnig mit abnehmender Wellenlänge _n nach kleineren Zenitdistanzen (entgegen den Angaben vonSchein undStoll, aber in Übereinstimmung mit der Theorie nachGötz undDobson). 2. Zum ersten Mal beziehen sich hier die Angaben aufeffektive Wellenlängen, die ihrerseits noch eine Funktion vonz sind. 3. Aus den Messungen bei den kleinsten, bisher noch nicht zugänglichen Wellenlängen folgt, daß deren Streuniveaus bei tiefen Sonnenständen oberhalb der Schichten größten Ozongehaltes liegen müssen.

---

Summary Until now the Umkehr-effect of Götz (inversion of intensity-ratio in zenith-scattered sunlight) has been investigated only by means of photocells or Geiger counters. The application of the multiplier as measuring instrument is of greater advantage because of its much higher sensitivity. Measurements can be made with it even with the sun's zenith distancesz>90° at shorter wave-lengths and with a narrower slit of the spectrometer (better monochromasy).Results: (1) The shift of the minima of the Umkehr-curvesi _n /i _m =f(z), erfolgtgleichsinnig mit abnehmender Wellenlänge _n shows the same sense with decreasing wave-length _n to shorter zenith distances (in opposition to the assertions ofSchein andStoll, but in agreement with the theory ofGötz andDobson). (2) For the first time the statements refer here to effective wave-lengths which in their turn are a function ofz. (3) From measurements of the shortest wave-lengths which until now have not yet been accessible, it can be concluded that at low sun heights the scattering level is situated above the layers with the greatest ozone content.

---

Résumé Jusqu'à présent le «Umkehreffect» deGötz, c. à d. l'effet d'inversion des relations d'intensité dans le rayonnement solaire zénithal diffus, n'a été étudié qu'au moyen de cellules photo-électriques ou de compteurs Geiger. En utilisant le multiplicateur comme instrument de mesure on obtient l'avantage d'une bien plus grande sensibilité. Ainsi il est possible de mesurer encore à des distances zénithales du soleil dez>90° par longueurs d'ondes plus courtes que jusqu'ici et avec une fente du spectromètre plus étroite, ce qui nous garantit une meilleure monochromasie.Résultats: 1) Le déplacement des minima des courbes d'inversioni _n /i _m =f(z), erfolgtgleichsinnig mit abnehmender Wellenlänge _n s'effectue dans le même sens que la diminution des longuers d'ondes _n , par distances zénithales plus petites (contrairement aux assertions deSchein et deStoll, mais en concordance avec la théorie selonGötz etDobson). — 2) Pour la première fois les données se rapportent ici à des longueurs d'ondes effectives qui, de leur côté, sont encore fonctions dez. — 3) Il ressort des mesures des plus courtes longueurs d'ondes, qui, jusqu'ici, restaient hors d'atteinte, que pour des hauteurs du soleil basses, leur niveau de fiffusion se trouve au-dessus des couches ayant le plus grand contenu d'ozone.

---

---

Mit 9 Textabbildungen     

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.     

Triggering of atmospheric circulations by moisture inhomogeneities of the earth's surface

A two-dimensional mesoscale soil-atmosphere model is used to simulate the triggering of atmospheric convection by horizontally varying soil water content. The variation is periodic with a wavelength between 4 and 40 km, which is considered a realistic scale for the variation of land surface characteristics. Three stages of convection can be clearly discerned: a short initial stage when convection sets in and where the size of the conective cells is determined by , a mature stage with well developed cells whose size is still determined by , and a decay/transformation stage, characterized by the formation of narrow regions of strong updrafts and wide regions of moderate downdrafts, independent of . Parameters relevant for the transition are given, and the importance of the feedback between soil and atmosphere is demonstrated. The dependence of convective parameters, e.g., height of the convective layer, vertical velocity and fluxes of heat and moisture on is investigated. The calculations of the mature stage are compared with the predictions of a linear model.