Estimation of the turbulent fluxes in the surface layer using the inertial dissipative method: A Monte Carlo error analysis

Inertial dissipation methods are based on measurements of the structure parameters of momentum, temperature and humidity. The most serious problems arise when experimental errors and uncertainties in the meteorological constants produce uncertainties in the calculated fluxes.In order to design an experimental layout to use an optical technique to measure fluxes, we need to perform an accurate error analysis. A Monte Carlo analysis is presented here to simulate the results and to analyze the propagation of statistical errors affecting the measurements. Results show that the coupling of inertial dissipation methods with optical methods is promising and a field experiment, based on the joint techniques, is feasible.This work was partially supported by M.P.I. (60%) grants of the Italian Ministry of Education.     

Regional surface fluxes from satellite-derived surface temperatures (AVHRR) and radiosonde profiles

Radiometric surface temperatures, derived from measurements by the AVHRR instrument aboard the NOAA-9 and the NOAA-11 polar orbiting satellites, were used in combination with wind velocity and temperature profiles measured by radiosondes, to calculate surface fluxes of sensible heat. The measurements were made during FIFE, the First ISLSCP (International Satellite Land Surface Climatology Project) Field Experiment, in a hilly tall grass prairie area of northeastern Kansas. The method of calculation was based on turbulent similarity formulations for the atmospheric boundary layer. Good agreement (r = 0.7) was obtained with reference values of sensible heat flux, taken as arithmetic means of measurements with the Bowen ratio method at six ground stations. The values of evaporation (latent heat fluxes), derived from these sensible heat fluxes by means of the energy budget, were also in good agreement (r = 0.94) with the corresponding reference values from the ground stations.     

The aerosol module in the INM RAS climate model

The aerosol module is included into the INM RAS climate model. The module computes the evolution of main aerosols: sea salt, mineral dust, sulfate aerosol, and black and organic carbon. Aerosol surface fluxes, advection, gravitational sedimentation, surface absorption, and scavenging by precipitation are taken into account to compute aerosol concentration variations. Model aerosol distribution is used to compute radiation fluxes. The ten-year run of the climate model is performed. The climatology of model aerosol is considered. The aerosol mass, integral source values, optical thickness, and radiative forcing are presented. The results are compared with the data of other models and observations.     

Estimation of turbulent sensible heat and momentum fluxes over a heterogeneous urban area using a large aperture scintillometer

The accurate determination of surface-layer turbulent fluxes over urban areas is critical to understanding urban boundary layer (UBL) evolution. In this study, a remote-sensing technique using a large aperture scintillometer (LAS) was investigated to estimate surface-layer turbulent fluxes over a highly heterogeneous urban area. The LAS system, with an optical path length of 2.1 km, was deployed in an urban area characterized by a complicated land-use mix (residential houses, water body, bare ground, etc.). The turbulent sensible heat (Q H) and momentum fluxes (τ) were estimated from the scintillation measurements obtained from the LAS system during the cold season. Three-dimensional LAS footprint modeling was introduced to identify the source areas ("footprint") of the estimated turbulent fluxes. The analysis results showed that the LAS-derived turbulent fluxes for the highly heterogeneous urban area revealed reasonable temporal variation during daytime on clear days, in comparison to the land-surface process-resolving numerical modeling. A series of sensitivity tests indicated that the overall uncertainty in the LAS-derived daytime Q H was within 20%-30% in terms of the influence of input parameters and the non-dimensional similarity function for the temperature structure function parameter, while the estimation errors in τ were less sensitive to the factors of influence, except aerodynamic roughness length. The 3D LAS footprint modeling characterized the source areas of the LAS-derived turbulent fluxes in the heterogeneous urban area, revealing that the representative spatial scales of the LAS system deployed with the 2.1 km optical path distance ranged from 0.2 to 2 km2 (a "micro-α scale"), depending on local meteorological conditions.     

Comparison of areally averaged turbulent fluxes over non-homogeneous terrain: Results from the EFEDA-field experiment

In June 1991 the EFEDA-field experiment (ECHIVAL Field Experiment in a Desertification-Threatened Area) was carried out in the Spanish province Castilla-La Mancha, to improve the understanding of the interactions between the soil, the vegetation and the atmosphere.Here results of energy balance studies at the Barrax site are given, one of the three intensively studied experimental sites within Castilla-La Mancha. This area is characterized by a large fraction of irrigated fields (40%) while the remaining 60% was fallow land at the end of June 1991. The energy balances over these two characteristic land-use classes totally differ. While for the irrigated fields the evapotranspiration is dominant, for the non-irrigated fields the sensible and the soil heat fluxes dominate and the latent heat flux is nearly negligible.In order to achieve areally averaged turbulent fluxes, surface, SODAR and aircraft data have been used. Comparing the surface fluxes from all three facilities, it can be found that:The sensible heat flux estimation from SODAR (_w-method) gives reliable results when a calibration of _w is done with another independent system (e.g. aircraft).Aircraft measurements in conjunction with energy budget methods yield surface fluxes of sensible heat, which are about 20% lower than the areally averaged values calculated by the surface measurements. The areally averaged latent heat fluxes from aircraft and surface measurements agree better than the sensible heat fluxes.     

Heat and moisture fluxes on the time scale of 20 to 60 days over the Indian monsoon area

Summary In this paper we have studied the low frequency variability of the sensible and latent heat flux over the Indian monsoon area. We have used an atmospheric energy budget (vertical integrated heat sources and moisture sinks), as well as the similarity theory in order to compute the surface fluxes on a darly basis. Mainly, the three following data sets were used: the First GARP Global Experiment analyzed data, the TIROS-N outgoing longwave radiation data and the Monsoon Experiment precipitation data.Our three main findings are the following. First, the variability of the temperature and the specific humidity at the surface is more important over the land than over the sea on the intraseasonal time scale (30% over land, but 20% over sea). For the wind an energy peak appears clearly around 30–40 days. The surface fluxes show an uneven variance percentage field (10% to 40%); the energy peaks stretch from 10 to 40 days. Second, the wind has a significant influence on the surface fluxes, except at some locations exclusively over the land areas. Of the temperature and the specific humidity, the temperature is the one which influences the fluxes the most. (This influence may be very strong over land.) The specific humidity may have a significant influence, over the land and sea, at the same time. Thus, one cannot neglect the influence of temperature and specific humidity over land on the intraseasonal time scale. Third, we have found a close relation between the propagation of low frequency waves and the propagation of surface flux patterns. This may suggest a feedback mechanism which relates surface processes to the northward propagation of these waves over India.With 17 FiguresOn leave from Etablissement d'études et de recherches méteorologiques Paris, France     

Optical attenuation in fog at a wavelength of 1.55 micrometers

As part of a series of studies on laser propagation for terrestrial free space optical (FSO) telecommunications or laser telecommunications, an experiment was conducted to determine the relationship between visibility in fog and optical attenuation (dB/km) at a laser wavelength of 1.55 μm. In the telecommunications industry, a semi-empirical equation, called the Kruse formula [Kruse, P.W., McGlauchlin, L., McQuistan, R.B., 1962. Elements of infrared technology: generation, transmission, and detection. John Wiley and Sons, New York] is typically used to calculate expected attenuation for a given meteorological visibility. The Kruse formula, however, was developed to relate meteorological visibility to optical attenuation over wavelengths from the visible to the near infrared (IR), and for dust and small particle aerosols with dimensions much smaller than the wavelength. Typically, suspended small aerosols have diameters that average about 0.1 μm while fog droplets have diameters that range upward from 2.5 μm with mean diameters that exceed 10 μm in some fogs. Therefore, application of the Kruse formula to attenuation in fog is not appropriate since fogs consist mainly of particles much larger than the laser wavelength. As part of the experiment, a transmissometer with an 85-m baseline and a dynamic range of 60 dB operated for thirteen months in an area prone to radiation fog. A commercial visibility sensor, similar to those used at airports, was located near the middle of the optical path of the transmissometer and operated over the same period. The largest attenuation measured at this site was just over 300 dB/km, corresponding to a visibility of 32 m. The key finding of the study is that the generally accepted Kruse formula relating visibility and optical attenuation may be too pessimistic at low visibilities, and actual attenuation values for a given visibility may be more than 20% lower than previously thought. At visibilities exceeding about 650 m, the Kruse formula gives a good estimate of optical attenuation.     

Lidar remote sounding of cirrus clouds and comparison of simulated fluxes with surface and METEOSAT observations

Three cirrus cloud cases have been remotely sounded near Paris by a ground-based backscatter lidar and broadband radiometers. Some cirrus properties (optical depth, emissivity, height) are derived from these measurements and used to compare radiative transfer calculations to surface and METEOSAT observations of broadband irradiances.For a useful comparison, the three cirrus cases were selected to have different morphologies and optical properties: June 29, 1993—thin cirrus cloud (thickness 1.5 km, optical depth 0.22); September 6, 1993—thick cirrus cloud (thickness 5 km, optical depth 2.7); and November 16, 1993—inhomogeneous and geometrically thick cirrus cloud (thickness 3.5–6.5 km) but optically thin (optical depth 0.82).At surface, the differences between measurements and model range from 1.5 to 4 Wm⁻² for longwave fluxes, and from 20 to 70 Wm⁻² for shortwave fluxes.At the top of the atmosphere, the differences between METEOSAT measurements and model are in fair agreement for longwave fluxes (up to 50 Wm⁻²). However, unexpected high differences are found for shortwave fluxes (up to 144 Wm⁻²) due to cirrus clouds heterogeneities and uncertainties in their microphysical properties and especially the occurrence of high reflectivity due to horizontally oriented ice crystals at the cloud top, which are not taken into account by the Model presently.     

Consequences of Incomplete Surface Energy Balance Closure for CO2 Fluxes from Open-Path CO2/H2O Infrared Gas Analysers

We present an approach for assessing the impact of systematic biases in measured energy fluxes on CO₂ flux estimates obtained from open-path eddy-covariance systems. In our analysis, we present equations to analyse the propagation of errors through the Webb, Pearman, and Leuning (WPL) algorithm [Quart. J. Roy. Meteorol. Soc. 106, 85–100, 1980] that is widely used to account for density fluctuations on CO₂ flux measurements. Our results suggest that incomplete energy balance closure does not necessarily lead to an underestimation of CO₂ fluxes despite the existence of surface energy imbalance; either an overestimation or underestimation of CO₂ fluxes is possible depending on local atmospheric conditions and measurement errors in the sensible heat, latent heat, and CO₂ fluxes. We use open-path eddy-covariance fluxes measured over a black spruce forest in interior Alaska to explore several energy imbalance scenarios and their consequences for CO₂ fluxes.     

Problems of determining snow albedo for the high arctic: Research note

Abstract

The measurement of radiation fluxes suffers from inaccuracies at low solar elevations and this poses a problem for determining the snow albedo at high latitudes. From the data of Resolute, NWT, three situations were observed: (1) an often‐reported situation when albedo decreases with increasing solar elevation, (2) an inverse situation when albedo increases with increasing solar elevation and (3) no obvious relationship. There were also cases when albedo exceeded 100%. The possible causes for such anomalous conditions or for erroneous albedos include instrument response deviating from the cosine law, instrument tilt, sensing of the sun by the inverted pyranometer and change in the spectral quality of incoming radiation with changing solar height. However, omission of the radiation values measured during the period of low solar elevation will not seriously affect the prediction of snowmelt. In this note, we have identified the anomalies and suggested possible causes; but further investigations are required to verify the causative mechanisms.     

Analysis of three-dimensional Eliassen-Palm fluxes in the lower stratosphere

Using the monthly mean NCEP/NCAR reanalysis dataset, the three-dimensional Eliassen-Palm (EP) fluxes of quasi-stationary wave propagation in the lower stratosphere were computed for each month from November to March for the period from 1958 to 2007. It is shown that the upward planetary wave propagation from the troposphere to the stratosphere generally occurs over the northern Eurasia, while their weak downward propagation is observed in Labrador and southern Greenland regions in the lower stratosphere. Interannual variations of the vertical EP fluxes also have the dipole-like spatial pattern with the opposite anomalies in the West and East hemispheres which are most prominent in January–February. Significant differences in the interaction of the zonal circulation of the stratosphere in the beginning of winter (November–December) and mid-to-late winter (January–March) are revealed. Intensification of the planetary waves’ penetration into the stratosphere in December causes changes in the stratospheric dynamics, creating the “preconditions” for the stratospheric warming appearances in January, but such a mechanism is not detected in February. In the years with the cold polar vortex, the “stratospheric bridge” is formed with the strengthening of the upward EP flux over the northern Eurasia and downward EP flux over the North Atlantic.     

Boreal summer intraseasonal variability simulated in the NCEP climate forecast system: insights from moist static energy budget and sensitivity to convective moistening

The NCEP Climate Forecast System (CFS) with the relaxed Arakawa Schubert (RAS, hereafter referred to as CTRL) convection scheme of Moorthi and Suarez exhibits better performance in representing boreal summer tropical intraseasonal variability as compared with a simulation using simplified Arakawa–Schubert scheme. The intraseasonal moist static energy (MSE) budget is analyzed in this version of the CFS model (CTRL), which produces realistic eastward and northward propagation characteristics. The moist and thermodynamic processes involved in the maintenance and propagation of the poleward moving intraseasonal oscillation (ISO) disturbances are examined here. Budget diagnostics show that horizontal MSE advection is the principal component of the budget, contributing to the poleward movement of the convection. The injection of MSE moistens the atmosphere north of the convective area causing the poleward movement of convection by destabilization of the atmosphere. The moistening process is mainly contributed by the climatological wind acting on the anomalous moisture gradient as confirmed from the examination of moisture advection equation. While surface enthalpy fluxes (consisting of radiative and surface turbulent heat fluxes) maintain the ISO anomalies, they oppose the MSE tendency due to horizontal advection thus regulating the poleward propagation characteristics. In addition, the model results show that wind–evaporation feedback dominates over cloud–radiation feedback for ISO propagation; this is in contrast to our estimates using the newly available European Centre for Medium Range Weather Forecasts Interim reanalysis. Sensitivity experiments suggest that intraseasonal variability in the CFS model with the RAS scheme is highly sensitive to the parameterization of both the shallow convection and the convective rain evaporation and downdrafts. Removal of these components adversely affects the propagation characteristics and greatly reduces the amplitude of intraseasonal variability. Our results support the primary importance of the moisture preconditioning ahead of the ISO and the physical relationship between moisture and precipitation. For realistic ISO simulations, models need to represent these features appropriately.     

The measurement of turbulent surface-layer fluxes by use of bichromatic scintillation

Scintillation measurements with a HeNe and a CO₂ laser were used to derive turbulent fluxes of heat and momentum in the surface layer. This was achieved by the structure constant or dissipation technique, i.e., by relating the measured structure constants and inner scales of refractive index fluctuations to structure constants of temperature fluctuations and dissipation rates of turbulent kinetic energy, respectively, and then assuming Monin-Obukhov similarity.The resulting heat fluxes agree well with measurements using the eddy correlation technique but for averaging periods of 10 min, the optical data show a much smoother and physically more plausible behaviour. The optically derived friction velocities are in good agreement with estimates derived from wind velocity and surface roughness. It was also observed that for stationary conditions, 1-min averaged optical measurements already provide good estimates for longer averaged heat and momentum fluxes.Even though some uncertainty remains about the empirical constants and Monin-Obukhov similarity expressions used, the method clearly proves to be of great value for monitoring surface-layer turbulence.     

Contrasting characteristics of the surface energy balance between the urban and rural areas of Beijing

A direct comparison of urban and rural surface energy balances, as well as a variety of other variables including incoming shortwave/longwave radiation and aerosol optical depth, is conducted for the Beijing metropolitan area. The results indicate that, overall, the urban area receives a smaller amount of incoming shortwave radiation but a larger amount of incoming longwave radiation. However, comparisons in the aerosol optical depth and cloud fraction at the two locations suggest that neither aerosol optical depth nor cloud fraction alone can explain the difference in the incoming shortwave radiation. The urban–rural differences in the incoming longwave radiation are unlikely to be caused by the presence of more abundant greenhouse gases over the urban area, as suggested by some previous studies, given that water vapor is the most dominant greenhouse gas and precipitable water is found to be less in urban areas. The higher incoming longwave radiation observed over the urban area is mostly likely due to the higher temperatures of the ambient air. The urban area is also found to always produce higher sensible heat fluxes and lower latent heat fluxes in the growing season. Furthermore, the urban area is associated with a larger amount of available energy(the sum of sensible and latent heat fluxes) than the rural area, except in May and October when evapotranspiration in the rural area significantly exceeds that in the urban area. This study provides observational evidence of urban–rural contrasts in relevant energy-balance components that plausibly arise from urban–rural differences in atmospheric and land-surface conditions.     

Scintillometer Intercomparison Study—Continued

An earlier study by one of the authors reported significant differences of up to 21% in linear regression slopes between six Kipp & Zonen large-aperture scintillometers. In this note, the consistency of this increasingly popular instrument for measuring sensible heat fluxes at the km scale was quantified by comparing measurements from four Scintec boundary-layer scintillometers and one large-aperture scintillometer over nearly identical transects. The Kipp & Zonen instrument’s sensible heat fluxes were more than 20% larger than those from the Scintec instruments, while the difference in regression slopes amongst the Scintec instruments was 3% or less.     

前向散射能见度仪校准技术

基于前向散射能见度仪的工作原理、结构特点及使用环境,分析了前向散射能见度仪的测量误差来源、产生测量误差原因及误差模型,归纳了前向散射能见度仪出厂标定方法、标定参数及标定参数修正方法。为确保能见度观测数据具有准确性,并根据将来量值传递需要,在分析了能见度量值传递参数基础上,提出了前向散射能见度仪的实验室标定条件和方法,梳理了前向散射能见度仪现场光学装置校准、校准板校准及人工校准方法,为今后建立能见度量值传递体系和编写相关校准规范提供参考。     

Monin–Obukhov Similarity Functions for the Structure Parameters of Temperature and Humidity

Monin–Obukhov similarity functions for the structure parameters of temperature and humidity are needed to derive surface heat and water vapour fluxes from scintillometer measurements and it is often assumed that the two functions are identical in the atmospheric surface layer. Nevertheless, this assumption has not yet been verified experimentally. This study investigates the dissimilarity between the turbulent transport of sensible heat and water vapour, with a specific focus on the difference between the Monin–Obukhov similarity functions for the structure parameters. Using two datasets collected over homogeneous surfaces where the surface sources of sensible heat and water vapour are well correlated, we observe that under stable and very unstable conditions, the two functions are similar. This similarity however breaks down under weakly unstable conditions; in that regime, the absolute values of the correlations between temperature and humidity are also observed to be low, most likely due to large-scale eddies that transport unsteadiness, advection or entrainment effects from the outer layer. We analyze and demonstrate how this reduction in the correlation leads to dissimilarity between the turbulent transport of these two scalars and the corresponding Monin–Obukhov similarity functions for their structure parameters. A model to derive sensible and latent heat fluxes from structure parameters without measuring the friction velocity is tested and found to work very well under moderately to strongly unstable conditions (−z/L > 0.5). Finally, we discuss the modelling of the cross-structure parameter over wet surfaces, which is crucial for correcting water vapour effects on optical scintillometer measurements and also for obtaining surface sensible and latent heat fluxes from the two-wavelength scintillometry.     

Sensitivity of Offshore Surface Fluxes and Sea Breezes to the Spatial Distribution of Sea-Surface Temperature

A series of numerical sensitivity experiments is performed to quantify the impact of sea-surface temperature (SST) distribution on offshore surface fluxes and simulated sea-breeze dynamics. The SST simulations of two mid-latitude sea-breeze events over coastal New England are performed using a spatially-uniform SST, as well as spatially-varying SST datasets of 32- and 1-km horizontal resolutions. Offshore surface heat and buoyancy fluxes vary in response to the SST distribution. Local sea-breeze circulations are relatively insensitive, with minimal differences in vertical structure and propagation speed among the experiments. The largest thermal perturbations are confined to the lowest 10% of the sea-breeze column due to the relatively high stability of the mid-Atlantic marine atmospheric boundary layer (ABL) suppressing vertical mixing, resulting in the depth of the marine layer remaining unchanged. Minimal impacts on the column-averaged virtual potential temperature and sea-breeze depth translates to small changes in sea-breeze propagation speed. This indicates that the use of datasets with a fine-scale SST may not produce more accurate sea-breeze simulations in highly stable marine ABL regimes, though may prove more beneficial in less stable sub-tropical environments.     

Thermal radiative fluxes through inhomogeneous cloud fields: a sensitivity study using a new stochastic cloud generator

We analyze the effects of flat and bumpy top, fractional and internally inhomogeneous cloud layers on large area-averaged thermal radiative fluxes. Inhomogeneous clouds are generated by a new stochastic model: the tree-driven mass accumulation process (tdMAP). This model is able to provide stratocumulus and cumulus cloud fields with properties close to those observed in real clouds. A sensitivity study of cloud parameters is done by analyzing differences between 3D fluxes simulated by the spherical harmonic discrete ordinate method and three “standard” models likely to be used in general circulation models: plane-parallel homogeneous cloud model (PPH), PPH with fractional cloud coverage model (FCPPH) and independent pixel approximation model (IPA). We show that thermal fluxes are strong functions of fractional cloud coverage, mean optical depth, mean geometrical thickness and cloud base altitude. Fluctuations of “in-cloud” horizontal variability in optical depth and cloud-top bumps have negligible effects in the whole. We also showed that PPH, FCPPH and IPA models are not suitable to compute thermal fluxes of flat top fractional inhomogeneous cloud layer, except for completely overcast cloud. This implies that horizontal transport of photon at thermal wavelengths is important when cloudy cells are separated by optically thin regions.