Effective parameters of surface energy balance in heterogeneous landscape

This paper addresses the problem of estimating surface fluxes at large scale over heterogeneous terrain, and the corresponding determination of effective surface parameters. Two kinds of formulation are used to calculate the fluxes of sensible and latent heat: the basic diffusion equations (Ohm's law type) and the Penman-Monteith equations. The strategy explored is based upon the principle of flux conservation, which stipulates that the average flux over a large area is simply the area-weighted mean of the contributions from the different patches making up the area. We show that the application of this strategy leads to different averaging schemes for the surface parameters, depending on the type of flux (latent heat, sensible heat) and on the type of formulation used to express the flux. It appears that the effective value of a given parameter must be appraised for each individual application, because it is not unique, but differs according to the magnitude being conserved and the equation used to express this magnitude. Numerical simulations are carried out to test over contrasted areas the aggregation procedures obtained. The areal fluxes estimated from these effective parameters, together with the areal fluxes calculated by means of a simple areal averaging of the parameters, are compared to the true average fluxes, calculated as area-weighted means of the elementary fluxes. The aggregation procedures obtained prove to be much more accurate for estimating areal fluxes and for closing the energy balance equation than those based upon simple areal averaging of the parameters.     

Design of land-air parameterization scheme (LAPS) for modelling boundary layer surface processes

Summary A land-air parametrization scheme (LAPS) describes mass, energy and momentum transfer between the land surface and the atmosphere. The scheme is designed as a software package which can be run as part of an atmospheric model or a stand-alone scheme. A single layer approach is chosen for the physical and biophysical scheme background. The scheme has six prognostic variables: two temperatures (one for the canopy vegetation and one for soil surface), one interception storage, and three soil moisture storage variables. The scheme's upper boundary conditions are: air temperature, water vapour pressure, wind speed, radiation and precipitation at some reference level within the atmospheric boundary layer. The sensible and latent heat are calculated using resistance representation. The evaporation from the bare soil is parametrized using the scheme. The soil part is designed as a three-layer model which is used to describe the vertical transfer of water in the soil.The performances of the LAPS scheme were tested using the results of meteorological measurements over a maize field at the experimental site De Sinderhoeve (The Netherlands). The predicted partitioning of the absorbed radiation into sensible and latent heat fluxes is in good agreement with observations. Also, the predicted leaf temperature agrees quite well with the observed values.With 9 Figures     

Determining the optimal soil temperature scheme for atmospheric modelling applications

The treatment of the land surface can have a significant impact on the performance of atmospheric models, influencing the surface energy balance and near surface atmospheric variables. In numerical weather prediction models it is especially important to reproduce the observed diurnal cycle in screen-level temperature, which requires an accurate representation of the surface temperature, and therefore an accurate and computationally efficient representation of soil heat storage and transfer. We present a technique for analysing the accuracy of numerical soil temperature schemes, and a methodology for choosing the optimal layer thicknesses for schemes with a given number of layers. Furthermore, the analysis suggests that first generation land surface schemes, which typically used a layer-average surface temperature, may be more accurate in this respect than the latest land surface schemes, which tend to use a skin surface temperature boundary condition.The British Crowns right to retain a non-exclusive royalty-free license in and to any copyright is acknowledged.     

Simple indices of global climate variability and change Part II: attribution of climate change during the twentieth century

Five simple indices of surface temperature are used to investigate the influence of anthropogenic and natural (solar irradiance and volcanic aerosol) forcing on observed climate change during the twentieth century. These indices are based on spatial fingerprints of climate change and include the global-mean surface temperature, the land-ocean temperature contrast, the magnitude of the annual cycle in surface temperature over land, the Northern Hemisphere meridional temperature gradient and the hemispheric temperature contrast. The indices contain information independent of variations in global-mean temperature for unforced climate variations and hence, considered collectively, they are more useful in an attribution study than global mean surface temperature alone. Observed linear trends over 1950–1999 in all the indices except the hemispheric temperature contrast are significantly larger than simulated changes due to internal variability or natural (solar and volcanic aerosol) forcings and are consistent with simulated changes due to anthropogenic (greenhouse gas and sulfate aerosol) forcing. The combined, relative influence of these different forcings on observed trends during the twentieth century is investigated using linear regression of the observed and simulated responses of the indices. It is found that anthropogenic forcing accounts for almost all of the observed changes in surface temperature during 1946–1995. We found that early twentieth century changes (1896–1945) in global mean temperature can be explained by a combination of anthropogenic and natural forcing, as well as internal climate variability. Estimates of scaling factors that weight the amplitude of model simulated signals to corresponding observed changes using a combined normalized index are similar to those calculated using more complex, optimal fingerprint techniques.     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.     

Heat fluxes and roll circulations over the western Gulf Stream during an intense cold-air outbreak

Turbulence and heat fluxes in the marine atmospheric boundary layer (MABL) for the roll vortex regime, observed during the Genesis of Atlantic Lows Experiment (GALE) over the western Gulf Stream, have been studied. The spectral analysis suggests that cloud streets (roll vortices) are vertically organized convection in the MABL having the same roll scale for both the cloud layer and subcloud layer, and that the roll spacing is about three times the MABL depth. The roll circulations contribute significantly to the sensible (temperature) and latent heat (moisture) fluxes with importance increasing upward. Near the MABL top, these fluxes are primarily due to roll vortices which transfer both sensible heat and moisture upward in the lower half of the convective MABL. Near the MABL top, the roll circulations transfer sensible heat downward and moisture upward in the clear thermal-street region, but roll vortices influenced by evaporative cooling can transfer sensible heat upward and moisture downward in the cloud-street region. Near the cloud-top, the upward buoyancy flux due to evaporative cooling is highly related to the roll circulations near the inversion.For the lower half of the MABL, the normalized temperature flux decreases upward more rapidly than the humidity flux, which is mainly because potential temperature () increases slightly upward while humidity (q) decreases slightly upward above the unstable surface layer. The gradient production (associated with the gradient) is a source for the temperature flux in the unstable surface layer but changes to a sink in the mixed layer, while the gradient production (associated with the q gradient) acts as a source for the humidity flux in both the unstable surface and mixed layers. The results suggest that the entrainment at the MABL top might affect the budgets of temperature and humidity fluxes in the lower MABL, but not in the unstable surface layer.Caelum Research Corporation, Silver Spring, MD, 20901, U.S.A.     

A Spectrum-Independent Procedure for Correcting Eddy Fluxes Measured with Separated Sensors

We investigate flux underestimates in eddy correlation measurements that are caused by horizontal separation of the sensors. A common eddy correlation setup consists of a sonic anemometer and a humidity sensor which, because of its bulk, must be placed some distance away from the sonic path, leading to a flux loss (of latent heat). Utilizing an additional fast temperature sensor placed near the humidity sensor, we develop a procedure for correcting for this loss. The procedure simultaneously corrects the sensible heat flux for the difference between true temperature and sonic temperature. Our correction procedure, which does not depend on the shape of the cospectrum, is then compared to the widely-used procedure following Moore (1986), which assumes a cospectral model ('Kansas Model). Both correction methods are applied to data collected within the internal boundary layer over a rice paddy, downwind of arid land. Under conditions of good fetch, they were found to agree well. Under poor fetch conditions, the model-based correction tended to be too small, while the spectrum-independent combined correction was robust. The latter is thus recommended for situations where the cospectral shape can be expected to deviate from the 'Kansas shape.     

Interannual variation of Gulf Stream heat transport in a high-resolution model forced by reanalysis data

The variability present in a 1/6th degree Atlantic ocean simulation forced by analysed wind stress and heat flux over a 20-year period is investigated by means of heat transport diagnostics. A section is defined which follows the Gulf Stream and its seaward extension, and transport of heat across this section is analysed to reveal the physical mechanisms responsible for intergyre heat exchanges on a variety of time scales. Heat transport across another section that crosses the Gulf Stream is also diagnosed to reveal the temporal behaviour of the gyre circulation. The Ekman response to wind stress variations accounts for the annual cycle and much of the interannual variability in both measures. For the intergyre heat transports, cancellation by transient-mean flow terms leads to a very weak annual cycle. Transient eddies account for approximately half the total intergyre transport of 0.7 Petawatts. They also account for a significant fraction of the interannual variability, but separate experiments with repeated-annual-cycle forcing indicate that the transient eddy component of the heat transport variability is internally generated. Links between the intergyre transport, the wind-driven gyre circulation, the surface heat budget and the atmospheric North Atlantic Oscillation are discussed.     

Evaluation of a land-surface scheme at Cabauw

Summary We study the response of the land-surface to prescribed atmospheric forcing for 31 May 1978 at Cabauw, Netherlands, using the land-surface scheme from the Coupled Atmospheric boundary layer-Plant-Soil (CAPS) model. Results from model runs show realistic daytime surface fluxes are produced using a canopy conductance formulation derived from Cabauw data (for 1987, a different year), and un-tuned parameterizations of root density (near-uniform with depth) and soil heat flux (reduced thermal conductivity through vegetation). Sensitivity of model-calculated surface heat fluxes to initial values of soil moisture is also examined. Results of this study provide the land-surface base state for a coupled land–atmosphere modeling study.     

Stability functions correct at the free convection limit and consistent for for both the surface and Ekman layers

For the thermal stability function _h used to calculate heat and moisture fluxes in the surface layer, we choose a formulation which has the theoretically correct free convection limit % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeikaiabgk% HiTGqaciaa-PhacaqGVaGaamitaiaabMcadaahaaWcbeqaaiabgkHi% TiaaigdacaGGVaGaaG4maaaaaaa!3DFE!\[{\rm{(}} - z{\rm{/}}L{\rm{)}}^{ - 1/3} \]. We then use the experimental result that z/L Ri to deduce a formulation with an exponent -1/6 for the momentum stability function _m. This formulation also resolves the matching problem at the interface between the surface and Ekman layers. The proposed functions are found to remain reasonably close to another formulation that is well supported by observations and has exponents -1/2 for _h and -1/4 for _m. The intent of the proposals is mainly to clarify and simplify the parameterization of the convective boundary layer in present day atmospheric models, without significantly altering the results.     

Atmospheric water vapor flux,bifurcation of the thermohaline circulation,and climate change

Latitudinal heat transport in the ocean and atmosphere represents a fundamental process of the Earth's climate system. The ocean component of heat transport is effected by the thermohaline circulation. Changes in this circulation, and hence latitudinal heat transport, would have a significant effect on global climate. Paleoclimate evidence from the Greenland ice cores and deep sea sediment cores suggests that during much of glacial time the climate system oscillated between two different states. Bimodal equilibrium states of the thermohaline circulation have been demonstrated in climate models. We address the question of the role of the atmospheric hydrological cycle on the global thermohaline circulation and the feedback to the climate system through changes in the ocean's latitudinal heat transport, with a simple coupled ocean-atmosphere energy-salt balance model. Two components of the atmospheric hydrological cycle, i.e., latitudinal water vapor transport and the net flux of water vapor from the Atlantic to the Pacific Ocean appear to play separate roles. If the inter-basin transport is sufficiently large, small changes in water vapor transport over the North Atlantic can effect bifurcation or a rapid transition between two different equilibria in the global thermohaline circulation; maximum difference between the modes occurs in the North Atlantic. If the inter-basin transport is from the Pacific to the Atlantic and sufficiently large, latitudinal vapor transport in the North Pacific controls the bifurcations, with maximum changes occurring in the North Pacific. For intermediate values of inter-basin transport, no rapid transitions occur in either basin. In the regime with vapor flux from the Atlantic to the Pacific, the on mode has strong production of deep water in the North Atlantic and a large flux of heat to the atmosphere from the high latitude North Atlantic. The off mode has strong deep water production in the Southern Ocean and weak production in the North Pacific. Heat transport into the high latitude North Atlantic by the ocean is reduced to about 20% of the on mode value. For estimated values of water vapor transport for the present climate the model asserts that while water vapor transport from the Atlantic to the Pacific Ocean is sufficiently large to make the North Atlantic the dominant region for deep water production, latitudinal water vapor transport is sufficiently low that the thermohaline circulation appears stable, i.e., far from a bifurcation point. This conclusion is supported to some extent by the fact that the high latitude temperature of the atmosphere as recorded in the Greenland ice cores has changed little over the last 9000 years.     

Surface renewal analysis for sensible and latent heat flux density

High frequency temperature measurements were recorded at five heights and surface renewal (SR) analysis was used to estimate sensible heat flux density (H) over 0.1 m tall grass. Traces of the temperature data showed ramp-like structures, and the mean amplitude and duration of these ramps were used to calculate H using structure functions. Data were compared with H values measured with a sonic anemometer. Latent heat flux density (E) was calculated using an energy balance and the results were compared with E computed from the sonic anemometer data. SR analysis provided good estimates of H for data recorded at all heights but the canopy top and at the highest measurement level, which was above the fully adjusted boundary layer.     

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.     

Atmospheric correction of infrared measurements of sea surface temperature using channels at 3.7, 11 and 12 Μm

Atmospheric effects upon the radiometric determination of surface temperature were studied for channels centered at 3.7, 11 and 12 m. The error due to the atmosphere is least for the channel centered at 3.7 m, which is a real advantage. The use of a linear combination of two or all three of these channels allows one to eliminate most of the atmospheric effect. If instrumental noise of from 0.1 to 0.2 K is accounted for in each channel, the best results are obtained by a combination of the two channels at 3.7 and 12 m.     

Changing thermal topography of the Baltimore-Washington corridor: 1950–1979

A study of the temperature field for the Baltimore-Washington region reveals that since 1950 there has been the development of an urban heat corridor. Trend surface analysis shows that there has been an inversion of the thermal topography in the region as a saddle of rising temperatures has emerged, replacing a trough of lowered temperatures through three decades. Steady population growth throughout the area is seen to be the most important contributor. The strengthening of the heat corridor is best expressed in the summer months and weakest during the winter months. As population continues to grow in this region, the thermal topography is sure to be modified even further.It is recommended that further study be devoted to temperature and precipitation changes in regions experiencing urban growth.     

Application of an infrared carbon dioxide and humidity instrument to studies of turbulent transport

A calibration equation and some results of the field performance of an infrared instrument, which is designed to measure simultaneous fluctuations of atmospheric carbon dioxide and water vapor, are described. Field observations show that the instrument is suitable for simultaneous measurement of turbulent fluxes of carbon dioxide and water vapor in conjunction with a sonic anemometer. Measured values of carbon dioxide and water vapor fluxes show diurnal variations characterized by crop activity with respect to assimilation, respiration and evapotranspiration. Carbon dioxide is transferred downward during the daytime and upward at night, while latent heat and sensible heat are transferred in the opposite sense. The non-dimensional gradient of carbon dioxide is expressed in the following form under weak unstable conditions: _c = (1 – 16 _v )^-1/2. Here, _v is the Monin-Obukhov stability parameter including the humidity effect. This relation was originally proposed for temperature and humidity. Thus, the results indicate that the turbulent mechanisms of carbon dioxide fluctuations are similar to those of other scalar entities. This is strongly supported by the high correlation coefficient found between fluctuations of carbon dioxide and temperature or humidity in the air layer over crop fields.     

Estimation of surface heat and momentum fluxes using the flux-variance method above uniform and non-uniform terrain

Eddy-correlation measurements above an uneven-aged forest, a uniform-irrigated bare soil field, and within a grass-covered forest clearing were used to investigate the usefulness of the fluxvariance method above uniform and non-uniform terrain. For this purpose, the Monin and Obukhov (1954) variance similarity functions were compared with direct measurements. Such comparisons were in close agreement for momentum and heat but not for water vapor. Deviations between measured and predicted similarity functions for water vapor were attributed to three factors: 1) the active role of temperature in surface-layer turbulence, 2) dissimilarity between sources and sinks of heat and water vapor at the ground surface, and 3) the non-uniformity in water vapor sources and sinks. It was demonstrated that the latter non-uniformity contributed to horizontal gradients that do not scale with the vertical flux. These three factors resulted in a turbulence regime that appeared more efficient in transporting heat than water vapor for the dynamic convective sublayer but not for the dynamic sublayer. The agreement between eddy-correlation measured and flux-variance predicted sensible heat flux was better than that for latent heat flux at all three sites. The flux-variance method systematically overestimated the latent heat flux when compared to eddy-correlation measurements. It was demonstrated that the non-uniformity in water vapor sources reduced the surface flux when compared to an equivalent uniform terrain subjected to identical shear stress, sensible heat flux, and atmospheric water vapor variance. Finally, the correlation between the temperature and water vapor fluctuations was related to the relative efficiency of surface-layer turbulence in removing heat and water vapor. These relations were used to assess critical assumptions in the derivation of the flux-variance formulation.     

Nearly steady convection and the boundary-layer budgets of water vapor and sensible heat

The budgets of water vapor and sensible heat in the convective atmospheric boundary (mixed) layer are analyzed by means of a simple slab approach adapted to steady large-scale advective conditions with radiation and cloud activity. The entrainment flux for sensible heat is assumed to be a linear function of the surface flux. The flux of water vapor at the top of the mixed layer is parameterized by extending the first-order Betts-Deardorff approach, i.e., by adopting linear changes for both the specific humidity and the flux across the mixed layer and across the inversion layer of finite thickness. In this way the dissimilarity of sensible heat and water vapor transport in the mixed layer can be taken into account. The experimental data were obtained from the Air Mass Transformation Experiment (AMTEX). The entrainment constant for sensible heat at the top of the mixed layer was found to have values similar to those observed in other weakly convective situations, i.e., around 0.4 to 0.6. This appears to indicate that the effect of mechanical turbulence was not negligible; however, the inclusion of this effect in the formulation did not improve the correlation. In contrast to the first-order approach, the zero-order approach, i. e., the jump equation commonly used for the flux of a scalar at the inversion, (ovwc ) _h = w_e c (where w _e is the entrainment velocity and c the concentration jump across the inversion), was found to be invalid and incapable of describing the data.     

Potential of continental-scale circulation for the determination of local daily surface variables

Summary This study deals with one aspect of statistical downscaling, viz. links between the continental-scale upper-air circulation and surface weather variables on a daily scale. The circulation-to-weather links are expressed in terms of multiple regression between either grid point values or intensities of circulation variability modes and weather elements, including temperature variables, relative humidity, cloudiness, sunshine duration, zonal and meridional wind components, precipitation, and atmospheric pressure. The upper-air circulation influences the surface pressure, temperature and zonal wind most strongly. The relative humidity, cloudiness and sunshine duration appear not to be connected with circulation in winter at all. The low-frequency part of circulation (i.e. processes with periods longer than 10 days) is most efficient in specifying surface weather variables. The circulation-to-weather links manifest considerable intra- and interdecadal variations, posing doubts on the applicability of the downscaling-from-circulation method of constructing climate change scenarios.With 11 Figures     

Statistical properties of optical refractive index fluctuations in the marine boundary layer

The optical refractive index fluctuation has been determined from measurements of turbulent temperature and humidity fluctuations in the atmospheric surface layer over the ocean. Probability density, spectral density and even- and odd-order structure functions of the refractive index fluctuation are presented. The spectral density exhibits a significant -5/3 inertial subrange as a result of the existence of a -5/3 subrange in the spectra of temperature and humidity fluctuations and in the temperature-humidity co-spectrum. The behaviour in the inertial subrange of fourth- and sixth-order structure functions of the refractive index is in reasonable agreement with that predicted by the analysis of Antonia and Van Atta (1975, 1978). The third-order structure functions of the refractive index exhibit an approximately linear dependence on separation in the inertial subrange, in agreement with a ramp model for the large coherent structure of the motion.