Water vapour flux profiles in the convective boundary layer

Summary Water vapour flux profiles in the atmospheric boundary layer have been derived from measurements of water vapour density fluctuations by a ground-based Differential Absorption Lidar (DIAL) and of vertical wind fluctuations by a ground-based Doppler lidar. The data were collected during the field experiment LITFASS-2003 in May/June 2003 in the area of Lindenberg, Germany. The eddy-correlation method was applied, and error estimates of ±50 W/m² for latent heat flux were found. Since the sampling error dominates the overall measurement accuracy, time intervals between 60 and 120 min were required for a reliable flux calculation, depending on wind speed. Rather large errors may occur with low wind speed because the diurnal cycle restricts the useful interval length. In the lower height range, these measurements are compared with DIAL/radar-RASS fluxes. The agreement is good when comparing covariance and error values. The lidar flux profiles are well complemented by tower measurements at 50 and 90 m above ground and by area-averaged near surface fluxes from a network of micrometeorological stations. Water vapour flux profiles in the convective boundary layer exhibit different structures mainly depending on the magnitude of the entrainment flux. In situations with dry air above the boundary layer a positive entrainment flux is observed which can even exceed the surface flux. Flux profiles which linearly increase from the surface to the top of the boundary layer are observed as well as profiles which decrease in the lower part and increase in the upper part of the boundary layer. In situations with humid air above the boundary layer the entrainment flux is about zero in the upper part of the boundary layer and the profiles in most cases show a linear decrease.     

江苏省人为热排放现状及未来变化趋势分析

依据1990~2013年江苏省能源消费统计资料和1990、1995、2000、2005、2010年的人口空间分布数据,估计江苏省人为热排放量、探讨其时空分布特征并预测其未来变化趋势。研究结果表明,20多年来江苏省的人为热排放量持续增长,从1990年全省平均的0.59 W/m~2增加到2013年的2.85 W/m~2。排放的空间分布也不均匀,总体上江苏南部地区的排放通量高于江苏北部。近些年江苏省的人为热排放已成为区域问题,高值区已连接成片,2010年江苏省大部分地区的人为热排放高于2.5 W/m~2,南部主要城市排放通量大于10 W/m~2、北部城市也多大于5 W/m~2。江苏各城市的城区人为热排放通量稳步增长,2002年后增长加速、增量最高可达到2 W m~(-2) a~(-1)。在2030年和2050年江苏省的平均人为热排放通量将分别达到5.7 W/m~2和9.1 W/m~2,可能对中国东部的气候和大气环境产生重要影响。     

The Turbulent Lagrangian Time Scale in Forest Canopies Constrained by Fluxes,Concentrations and Source Distributions

One-dimensional Lagrangian dispersion models, frequently used to relate in-canopy source/sink distributions of energy, water and trace gases to vertical concentration profiles, require estimates of the standard deviation of the vertical wind speed, which can be measured, and the Lagrangian time scale, T _L , which cannot. In this work we use non-linear parameter estimation to determine the vertical profile of the Lagrangian time scale that simultaneously optimises agreement between modelled and measured vertical profiles of temperature, water vapour and carbon dioxide concentrations within a 40-m tall temperate Eucalyptus forest in south-eastern Australia. Modelled temperature and concentration profiles are generated using Lagrangian dispersion theory combined with source/sink distributions of sensible heat, water vapour and CO₂. These distributions are derived from a multilayer Soil-Vegetation-Atmospheric-Transfer model subject to multiple constraints: (1) daytime eddy flux measurements of sensible heat, latent heat, and CO₂ above the canopy, (2) in-canopy lidar measurements of leaf area density distribution, and (3) chamber measurements of CO₂ ground fluxes. The resulting estimate of Lagrangian time scale within the canopy under near-neutral conditions is about 1.7 times higher than previous estimates and decreases towards zero at the ground. It represents an advance over previous estimates of T _L , which are largely unconstrained by measurements.     

Organised Motion and Radiative Perturbations in the Nocturnal Canopy Sublayer above an Even-Aged Pine Forest

Using time series measurements of velocity, carbon dioxideand water vapour concentration, and temperature collected justabove a 15 m tall even-aged pine forest, we quantify the roleof organized motion on scalar and momentum transport withinthe nocturnal canopy sublayer (CSL). We propose a frameworkin which the nocturnal CSL has two end-members, bothdominated by organised motion. These end-members representfully developed turbulent flows at near-neutral or slightly stablestratification and no turbulence for very stable stratification.Our analysis suggests that ramps dominate scalar transport fornear-neutral and slightly stable conditions, while linear canopywaves dominate the flow dynamics for very stable conditions.For intermediate stability, the turbulence is highly damped andoften dominated by fine scale motions. Co-spectral analysissuggests that ramps are the most efficient net scalar mass-transportingagent while linear canopy waves contribute little to net scalartransport between the canopy and atmosphere for averagingintervals that include complete wave cycles. However, canopywaves significantly contribute to the spectral properties of thescalar time series. Ramps are the most frequently occurringorganised motion in the nocturnal CSL for this site.Numerous night-time runs, however, resided between thesetwo end-members. Our analysis suggests that whenradiative perturbations are sufficient large (>20 W m^-2 innet radiation), the flow can switch from being highly dampedfine-scale turbulence to being organized with ramp-like properties. We also found that when ramps are already the dominant eddymotion in the nocturnal CSL, radiative perturbations have aminor impact on scalar transport. Finally, in agreement withprevious studies, we found that ramps and canopy waves havecomparable length scales of about 30–60 metres. Consequencesto night-time flux averaging are also discussed.     

Evaluation of Mediterranean Sea water and heat budgets simulated by an ensemble of high resolution regional climate models

Air-sea heat and freshwater water fluxes in the Mediterranean Sea play a crucial role in dense water formation. Here, we compare estimates of Mediterranean Sea heat and water budgets from a range of observational datasets and discuss the main differences between them. Taking into account the closure hypothesis at the Gibraltar Strait, we have built several observational estimates of water and heat budgets by combination of their different observational components. We provide then three estimates for water budget and one for heat budget that satisfy the closure hypothesis. We then use these observational estimates to assess the ability of an ensemble of ERA40-driven high resolution (25 km) Regional Climate Models (RCMs) from the FP6-EU ENSEMBLES database, to simulate the various components, and net values, of the water and heat budgets. Most of the RCM Mediterranean basin means are within the range spanned by the observational estimates of the different budget components, though in some cases the RCMs have a tendency to overestimate the latent heat flux (or evaporation) with respect to observations. The RCMs do not show significant improvements of the total water budget estimates comparing to ERA40. Moreover, given the large spread found in observational estimates of precipitation over the sea, it is difficult to draw conclusions on the performance of RCM for the freshwater budget and this underlines the need for better precipitation observations. The original ERA40 value for the basin mean net heat flux is ?15 W/m² which is 10 W/m² less than the value of ?5 W/m² inferred from the transport measurements at Gibraltar Strait. The ensemble of heat budget values estimated from the models show that most of RCMs do not achieve heat budget closure. However, the ensemble mean value for the net heat flux is ?7 ± 21 W/m², which is close to the Gibraltar value, although the spread between the RCMs is large. Since the RCMs are forced by the same boundary conditions (ERA40 and sea surface temperatures) and have the same horizontal resolution and spatial domain, the reason for the large spread must reside in the physical parameterizations. To conclude, improvements are urgently required to physical parameterizations in state-of-the-art regional climate models, to reduce the large spread found in our analysis and to obtain better water and heat budget estimates over the Mediterranean Sea.     

Observations of atmospheric downward radiation in the Tokyo area

The size and variation of the urban increase of atmospheric downward radiation have been observed in the Tokyo area. Downward radiation flux, air temperature and humidity were observed by automobile traverses along the expressways across the greater Tokyo area from the Kanagawa (southwest) to the Chiba (east) sides. The observations were carried out on clear and rather calm nights in the winters of 1975–1978. The results show that Tokyo receives increased downward radiation in comparison with the amount received in the surrounding rural areas. The area of increased downward flux extends about 40 km at its widest along the observation route when the winds are weak, and the rural-urban increase of downward radiation becomes about 8% on the Kanagawa side and about 10% on the Chiba side. The area of increased downward flux is close to that of the urban heat island, where an urban decrease of water vapour content is also observed. The present analysis shows that the large urban increase of radiation flux is due to the difference between the vertical temperature profiles over the urban and rural regions. The development of a surface temperature inversion over the rural area and the urban decrease in water vapour content provide a rather smaller rural-urban difference in downward radiation compared with the development of a heat island.     

Observations of Atmospheric Downward Radiation in the Tokyo Area

The size and variation of the urban increase of atmospheric downward radiation have been observed in the Tokyo area. Downward radiation flux, air temperature and humidity were observed by automobile traverses along the expressways across the greater Tokyo area from the Kanagawa (southwest) to the Chiba (east) sides. The observations were carried out on clear and rather calm nights in the winters of 1975–1978. The results show that Tokyo receives increased downward radiation in comparison with the amount received in the surrounding rural areas. The area of increased downward flux extends about 40 km at its widest along the observation route when the winds are weak, and the rural-urban increase of downward radiation becomes about 8% on the Kanagawa side and about 10% on the Chiba side. The area of increased downward flux is close to that of the urban heat island, where an urban decrease of water vapour content is also observed. The present analysis shows that the large urban increase of radiation flux is due to the difference between the vertical temperature profiles over the urban and rural regions. The development of a surface temperature inversion over the rural area and the urban decrease in water vapour content provide a rather smaller rural-urban difference in downward radiation compared with the development of a heat island.     

从FY-4静止气象卫星估算晴空地表下行长波辐射通量的反演模式

通过红外辐射传输模拟,计算了446183条全球晴空大气廓线的地表下行长波辐射通量以及FY-4成像仪通道亮温,应用统计回归分析模拟结果,建立了通道亮温与低层大气有效辐射温度的回归关系、大气柱总可降水量与低层大气比辐射率的回归关系,依据黑体辐射定律建立了由卫星观测估算地表下行长波辐射通量的反演模式。初步地,模式应用于FY-4成像仪代理资料Meteosat-8卫星的SEVIRI（Spring Enhanced Visible and Infared Imager）仪器观测数据,估算了2006年8月1日00:00（协调世界时,下同）、06:00、12:00、18:00地理范围在（45°S～45°N,45°W～45°E）的地表下行长波辐射通量,结果与用ECMWF 6 h预报场资料经验计算的通量相比,系统均方根误差（RMS）依次为12.1、13.0、20.7、12.5 W/m2,相关系数分别为0.9256、0.9291、0.9042、0.9325,相比于GOES-R同类产品（RMS=13.7 W/m2）,这一直接把卫星通道亮温与低层大气温度相联系的反演模式,其精度性能达到了应用研究对其反演产品的质量要求。     

Radiative surface temperature and energy balance of a wheat canopy

The sensible heat loss from a stand of winter wheat was calculated from radiometric measurements of crop surface temperature, measurements of air temperature, and an atmospheric resistance to momentum transfer; corresponding latent heat flux was obtained through the energy balance equation. These estimates of sensible and latent heat were compared with fluxes from the Bowen Ratio method. When radiative temperature was derived using a measured canopy emissivity of 0.98, calculations of sensible heat flux were systematically 50–100 W m^-2 less than Bowen Ratio values. The two techniques agreed more closely when an apparent emissivity of 0.96 was used with an apparent reflectivity of 0.03. The mean difference between the estimates of latent heat flux was then -16 ± 32 W m^-2.The surface temperature method showed less systematic error in comparison with the Bowen Ratio values than did estimates using the aerodynamic method.On leave from: University of Nottingham, School of Agriculture, Loughborough LE12 5RD.     

An error analysis of profile flux,stability, and roughness length measurements made in the marine atmospheric surface layer

An analytical error analysis of profile-derived fluxes of heat, moisture, and momentum, along with stability and roughness length, is performed using the accuracies of the constituent temperature, humidity, and wind speed measurements. Five experiments, representing more than two thirds of the existing marine profile data presently contained in the literature, are compared. Much of the profile data examined was used to develop the transfer coefficients presently employed by a large number of competing bulk aerodynamic flux schemes. Depending upon the experiment, typical profile-method measurement errors were found to range from 15 to 35% for a sensible heat flux of ± 10 W m^-2, from 15 to 105% for a latent heat flux of ± 100 W m^-2, from 10 to 40% for a stress of 0.05 N m^-2, from 15 to 60% for a Monin-Obukhov stability of ± 0.05, and from 25 to 100% for a roughness length of 2 × 10^-4 m. Smaller magnitude flux values were found to contain typical errors as large as 100% for sensible heat flux, 300% for latent heat flux, and 60% for stress.     

Boundedness of Turbulent Temperature Probability Distributions,and their Relation to the Vertical Profile in the Convective Boundary Layer

Higher-order moments, minima and maxima of turbulent temperature and water vapour mixing ratio probability density functions measured with an eddy-covariance system near the ground were related to each other and to vertical boundary-layer profiles of the same scalars obtained through airborne soundings. The dependence of kurtosis on squared skewness showed a kurtosis intercept below the Gaussian expectation, suggesting a compression of the probability density function by the presence of natural boundaries. This hypothesis was corroborated by comparing actual minima and maxima of turbulent fluctuations to estimates obtained from the first four sample moments by fitting a four-parameter beta distribution. The most sharply defined boundaries were found for the minima of temperature datasets during the day, indicating that negative temperature fluctuations at the sensor are limited by the availability of lower temperatures in the boundary layer. By comparison to vertical profiles, it could be verified that the turbulent minimum of temperature near the ground is close to the minimum of potential temperature in the boundary layer. The turbulent minimum of water vapour mixing ratio was found to be equal to the mixing ratio at a height above the minimum of the temperature profile. This height roughly agrees with the top of the non-local unstable domain according to bulk Richardson number profiles. We conclude that turbulence statistics measured near the surface cannot be solely explained by local effects, but contain information about the whole boundary layer including the entrainment zone.     

A Lagrangian Model to Predict the Modification of Near-Surface Scalar Mixing Ratios and Air–Water Exchange Fluxes in Offshore Flow

A model was developed to predict the modification with fetch in offshore flow of mixing ratio, air–water exchange flux, and near-surface vertical gradients in mixing ratio of a scalar due to air–water exchange. The model was developed for planning and interpretation of air–water exchange flux measurements in the coastal zone. The Lagrangian model applies a mass balance over the internal boundary layer (IBL) using the integral depth scale approach, previously applied to development of the nocturnal boundary layer overland. Surface fluxes and vertical profiles in the surface layer were calculated using the NOAA COARE bulk algorithm and gas transfer model (e.g., Blomquist et al. 2006, Geophys Res Lett 33:1–4). IBL height was assumed proportional to the square root of fetch, and estimates of the IBL growth rate coefficient, α, were obtained by three methods: (1) calibration of the model to a large dataset of air temperature and humidity modification over Lake Ontario in 1973, (2) atmospheric soundings from the 2004 New England Air Quality Study and (3) solution of a simplified diffusion equation and an estimate of eddy diffusivity from Monin–Obukhov similarity theory (MOST). Reasonable agreement was obtained between the calibrated and MOST values of α for stable, neutral, and unstable conditions, and estimates of α agreed with previously published parametrizations that were valid for the stable IBL only. The parametrization of α provides estimates of IBL height, and the model estimates modification of scalar mixing ratio, fluxes, and near-surface gradients, under conditions of coastal offshore flow (0–50 km) over a wide range in stability.     

昆明城市热岛效应的数值模拟研究

本文利用WRF（V3.9.1）模式中耦合Noah/SLUCM方案作为Control试验,研究了土地利用类型（Md04试验）、陆面过程（NoUCM试验）和湖泊（Nolake试验）对城市热岛强度及昆明城市气象要素水平、垂直的时空分布影响。主要结论如下:（1）四个试验城市热岛强度的平均日变化趋势相似,白天城市热岛强度较弱、夜间较强,在20时（北京时,下同）左右达到最大值。城市冠层（湖泊）对城市热岛有较明显的减（增）温,Control-NoUCM（Nolake）试验中,平均日最大差值为?0.79°C（+1.07°C）。（2）从能量平衡方程分析Control-Md04试验,感热（潜热）通量的差值为+46.18（?79.71）W m?2,潜热通量释放大于感热通量的绝对值。Control-NoUCM试验中,感热（潜热）通量的差值为?40.88（+29.60）W m?2;因NoUCM试验未考虑几何建筑物储热与遮挡,太阳辐射大部分被地表所吸收,导致感热通量偏大。（3）四种试验中,15（07）时边界层高度达到最大（小）值。NoUCM（Nolake）试验中城市边界层高度分别降低103 m（32 m）左右,而Md04试验中城市边界层高度增加102 m左右。（4）湖泊（滇池）对城市热岛环流影响的试验表明,湖泊上空垂直运动较弱,但水平方向湖陆风较大,这有利于向城市输送水汽,增加干空气湿度,使城市中空气的水汽含量增加,同时增大潜热能量释放,降低感热通量,减小了垂直温度梯度。     

Airborne Investigation Of Arctic Boundary-Layer Fronts Over The Marginal Ice Zone Of The Davis Strait

Three aircraft-based studies of boundary-layer fronts (BLFs) werecarried out during the experiment KABEG in April 1997near the sea-ice edge over the Davis Strait. The zone of strongestcross-front horizontal gradients showed a typical length scaleof 20 km, while the along-front scale was observed to beseveral hundreds of kilometres.The observed BLFs were stronger than the few previously reportedcases. Horizontal gradients of potential temperature and specifichumidity ranged up to 3 K or 0.25 g kg^-1over 20 km, respectively.Low-level winds were around 15 m s^-1 parallel to the ice edge.The capping inversion sloped from between250 and 400 m over sea ice to between 400 and 700 m over ocean.For two BLF cases turbulent fluxes and energy budgets are calculated.Turbulent energy fluxes show a factor 2 to 3 contrast acrossthe ice edge and range from 15 to 50 W m^-2 over sea iceand from 50 to 100 W m^-1 over open ocean.The mean boundary-layer energy budgets are dominated bycold, dry horizontal advection, which is exceededby vertical heat flux convergence.The momentum budgets are dominated by pressure gradient force,Coriolis force and momentum flux divergence.     

临近绿洲的沙漠上空近地面层内水汽输送特征

1991年7月在HEIFE实验区沙漠站进行了一次水汽输送观测实验。观测发现近地面层内比湿廓线常常出现极小值，该极小值高度为零水汽通量面的高度。这个高度以上水汽向下输送，廓线呈逆湿特征，这个高度以下水汽向上输送，廓线呈蒸发特征。沙漠上空近地面层内水汽不是常通量，在一个固定高度上观测到的水汽通量和“潜热通量”不能代表地面的蒸发量和潜热通量。     

Comparisons of xylem sap flow and water vapour flux at the stand level and derivation of canopy conductance for Scots pine

Summary Simultaneous measurements of xylem sap flow and water vapour flux over a Scots pine (Pinus sylvestris) forest (Hartheim, Germany), were carried out during the Hartheim Experiment (HartX), an intensive observation campaign of the international programme REKLIP. Sap flow was measured every 30 min using both radial constant heating (Granier, 1985) and two types of Cermak sap flowmeters installed on 24 trees selected to cover a wide range of the diameter classes of the stand (min 8 cm; max 17.5 cm). Available energy was high during the observation period (5.5 to 6.9 mm.day^–1), and daily cumulated sap flow on a ground area basis varied between 2.0 and 2.7 mm day^–1 depending on climate conditions. Maximum hourly values of sap flow reached 0.33 mm h^–1, i.e., 230 W m^–2.Comparisons of sap flow with water vapour flux as measured with two OPEC (One Propeller Eddy Correlation, University of Arizona) systems showed a time lag between the two methods, sap flow lagging about 90 min behind vapour flux. After taking into account this time lag in the sap flow data set, a good agreement was found between both methods: sap flow = 0.745^* vapour flux,r ² = 0.86. The difference between the two estimates was due to understory transpiration.Canopy conductance (g _c ) was calculated from sap flow measurements using the reverse form of Penman-Monteith equation and climatic data measured 4 m above the canopy. Variations ofg _c were well correlated (r ² = 0.85) with global radiation (R) and vapour pressure deficit (vpd). The quantitative expression forg _c =f (R, vpd) was very similar to that previously found with maritime pine (Pinus pinaster) in the forest of Les Landes, South Western France.With 6 Figures     

Aircraft measurements of boundary-layer turbulence over the central Equatorial Pacific Ocean

We analyze gust probe measurements obtained in the convective boundary layer over the central equatorial Pacific as part of the Equatorial Pacific Ocean Climate Studies (EPOCS) program. From the lowest level flights, the bulk transfer coefficients are found to be 1.1 × 10^-3 and 1.4 × 10^-3 for latent and sensible heat fluxes, respectively. Vertical profiles of water vapor density, potential temperature and wind velocity are obtained as are the profiles of the fluxes of latent and sensible heat and momentum. From the extrapolated profiles, we obtain surface fluxes of 120 W m^-2 and 13 W m^-2 for latent and sensible heat, respectively, and 0.11 N m^-2 for momentum. The 10 m drag coefficient is 1.5 × 10^-3. Two convergence boxes, north and south of the ITCZ, are analyzed. Enhanced convergence is found in the northern trades relative to the southern trades. The advective acceleration is found to be comparable in magnitude to the other terms in the horizontal equations of motion.     

Vertical eddy diffusivity estimations in Swan river estuary

A field experiment is described in which diapycnal diffusivity is estimated by direct and indirect methods in Swan river estuary, Perth, Western Australia. The microstructure profiles were collected in a narrow and straight part of this estuary using a portable flux profiler (PFP). The profiles were segmented into stationary parts and the rate of vertical eddy diffusivity was estimated for the segments within the pycnocline. The direct measurement showed that instantaneous flux could be positive or negative with a low net rate of vertical eddy diffusivity for mass of about 6.5 × 10⁻⁸ m² s⁻¹. All the indirect methods overestimated this rate. However, within the indirect methods, the method of Osborn yielded the poorest estimate while the method of Osborn and Cox gave the closest estimate.     

Convective characteristics of the nocturnal urban boundary layer as observed with Doppler sodar and Raman lidar

Convective plume patterns, characteristic of clear sky and light wind daytime boundary layers over land, were observed for two nights with a tri-axial Doppler sodar operated in the central area of Rome during the summer of 1994. An urban heat island effect, combined with a continuation of a breeze from the sea late into night during both days, is believed to be responsible for the observed nocturnal convection. Estimates of the surface heat flux and the vertical velocity scaling parameter are obtained from profiles of vertical velocity variance, and the Raman lidar water vapor measurements are used to obtain the humidity scaling parameter. Convective scaling results for vertical wind and humidity fairly agree with the results of other experiments and models. On the basis of available measurements, it appears that mixed-layer similarity formulations used to characterize the daytime convective boundary layer over horizontally homogeneous surfaces can also be applied to the nocturnal urban boundary layer during periods of reasonable convective activity.     

Mean Profiles of Moisture Fluxes in Snow-Filled Boundary Layers

Profiles of moisture fluxes have been examined for convective boundary layers containing clouds and snow, using data derived from aircraft measurements taken on four dates during the 1983/1984 University of Chicago lake-effect snow project. Flux profiles were derived from vertical stacks of aircraft cross-wind flight legs taken at various heights over Lake Michigan near the downwind shore. It was found that, if ice processes are taken into account, profiles of potential temperature and water content were very similar to those presented in past studies of convective boundary layers strongly heated from below. Profiles of total water content and equivalent potential temperature adjusted for ice were nearly invariant with height, except very near the top of the boundary layer, suggesting that internal boundary-layer mixing processes were rapid relative to the rates at which heat and vapour were transported into the boundary layer through entrainment and surface fluxes. Ice was found to play a significant, measurable role in boundary-layer moisture fluxes. It was estimated that 40 to 57% of the upward vapour flux was returned to the surface in the form of snow, converting about 45 to 64% of the surface latent heat flux into sensible heat in the snow-producing process. Assuming advective fluxes are relatively small (thought to be appropriate after the first few tens of km over the lake as suggested by past studies), the boundary layer was found to warm at a rate faster than could be explained by surface heat fluxes and latent heat releases alone, the remainder of the heating presumably coming from radiational processes and entrainment. Discussions of moisture phase change processes throughout the boundary layer and estimates of errors of these flux measurements are presented.