A new parameterization of eddy diffusivities for nocturnal boundary-layer modeling

Exchange coefficients and mixing lengths under stable stratification have been studied through measurements of mean wind velocity and temperature in the nocturnal boundary layer. For values of the gradient Richardson number lower than 0.15, our measurements fit well the relation of Delage (1974). Beyond Ri = 0.15, the decrease of mixing length is much slower. So a new parameterization of turbulent exchanges is suggested. When introduced in a model of the nocturnal boundary layer, it results in a thickening of the turbulent and inversion layers.     

Upper-ocean heat budget and ocean eddy transport in the south-east Pacific in a high-resolution coupled model

We present an analysis of the oceanic heat advection and its variability in the upper 500 m in the southeastern tropical Pacific (100W–75W, 25S–10S) as simulated by the global coupled model HiGEM, which has one of the highest resolutions currently used in long-term integrations. The simulated climatology represents a temperature advection field arising from transient small-scale (<450 km) features, with structures and transport that appear consistent with estimates based on available observational data for the mooring at 20S, 85W. The transient structures are very persistent (>4 months), and in specific locations they generate an important contribution to the local upper-ocean heat budget, characterised by scales of a few hundred kilometres, and periods of over a year. The contribution from such structures to the local, long-term oceanic heat budget however can be of either sign, or vanishing, depending on the location; and, although there appears some organisation in preferential areas of activity, the average over the entire region is small. While several different mechanisms may be responsible for the temperature advection by transients, we find that a significant, and possibly dominant, component is associated with vortices embedded in the large-scale, climatological salinity gradient associated with the fresh intrusion of mid-latitude intermediate water which penetrates north-westward beneath the tropical thermocline.     

A simple model for the adjustment of velocity spectra in unstable conditions downstream of an abrupt change in roughness and heat flux

A simple model for the equilibrium spectra of velocity fluctuations in the unstable surface layer is developed. All three component spectra are written as a sum of two spectra For the horizontal spectra, the two parts scale with % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaacIcacaWGUb% GaamOEamaaBaaaleaacaWGPbaabeaakiaac+cacaWG1bGaaiilaiaa% dQhadaWgaaWcbaGaamyAaaqabaGccaGGVaGaamitaiaacMcaaaa!4232!\[(nz_i /u,z_i /L)\]and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaad6gacaWG6b% Gaai4laiaadwhaaaa!3B5E!\[nz/u\], respectively; the vertical spectrum can be written entirely as a function of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaad6gacaWG6b% Gaai4laiaadwhaaaa!3B5E!\[nz/u\] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadQhacaGGVa% Gaamitaaaa!3A42!\[z/L\].The equilibrium spectra are utilized as part of a model describing the development of velocity spectra downwind of a change in surface roughness and heat flux.Results are shown for the streamwise component and compared with hotwire measurements from the RISØ 78 experiment. The model shows excellent agreement with the measurements.on leave of absence from RISØ National Laboratory Roskilde, Denmark.     

Mid-latitude atmospheric responses to heat and vorticity forcing using a linear baroclinic model

Based on diagnostic analysis of reanalysis data for 58-year, the distribution characteristics of decadal variability in diabatic heating, transient eddy heating and transient eddy vorticity forcing related to the sea surface temperature (SST) anomalies over the North Pacific, as well as their relationship with anomalous atmospheric circulation have been investigated in this paper. A linear baroclinic model(LBM) was used to investigate atmospheric responses to idealized and realistic heat and vorticity forcing anomalies, and then to compare relative roles of different kinds of forcing in terms of geopotential height responses. The results illustrate that the responses of atmospheric height fields to the mid-latitude heating can be either baroclinic or barotropic. The response structure is sensitive to the relative horizontal location of heating with respect to the background jet flow, as well as to the vertical profile of heating. The response to the idealized deep heating over the eastern North Pacific, mimicking the observed heating anomaly, is baroclinic. The atmospheric response to the mid-latitude vorticity forcing is always barotropic, resulting in a geopotential low that is in phase with the forcing. The atmospheric responses to the realistic heat and vorticity forcing show the similar results, suggesting that diabatic heating, transient eddy heating and transient eddy vorticity forcing can all cause atmospheric anomalies and that the vorticity forcing plays a relatively more important role in maintaining the equivalent-barotropic structure of geopotential height anomalies.     

一个用于气候模式的简单冻土过程参数化方案的建立和检验

在NCAR／LSM的基础上，发展了一个简单的冻土过程参数化方案，并使用苏联6个站的水气象观测资料考察了耦合了新方案模式的气候效应。在新方案中，加入了对含冰量的求解和在相变过程中的能量变化;并使用Johanson的方案替代了模式中原有的土壤导热率的参数化方案，考虑了含冰量对土壤水热性质的影响。原模式和改进后模式的模拟结果的比较得到，冻土过程方案能够合理的模拟土壤列中的能量收支及水热性质随含冰量的变化。随着入渗的减少和径流的增加，春季的土壤湿度减小。因此，热通量的分配和土壤温度也产生了相应的变化。     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.     

A simple parameterization of sub-grid scale open water for climate models

The effects of small fractions ( < 30%) of open water covering a grid element are currently neglected even in atmospheric general circulation models (AGCMs) which incorporate complex land surface parameterization schemes. Here, a method for simulating sub-grid scale open water is proposed which permits any existing land surface model to be modified to account for open water. This new parameterization is tested as an addition to an advanced land surface scheme and, as expected, is shown to produce general increases in the surface latent heat flux at the expense of the surface sensible heat flux. Small changes in temperature are associated with this change in the partitioning of available energy which is driven by an increase in the wetness of the grid element. The sensitivity of the land surface to increasing amounts of open water is dependent upon the type of vegetation represented. Dense vegetation (with a high leaf area index) is shown to complicate the apparently simple model sensitivity and indicates that more advanced methods of incorporating open water into AGCMs need to be considered and compared against the parameterization suggested here. However, the sensitivity of one land surface model to incorporating open water is large enough to warrant consideration of its incorporation into climate models.     

斜压大气中的特征波动和地转适应过程

本文对斜压大气的小扰动方程作了严格和全面的数学分析:论证了关于垂直边界条件的正确提法以及初始场一些特殊问题,澄清了迄今不正确的和混乱的提法;给出了解的存在、唯一和稳定性的证明。在这基础上,对大气的特征波动作了分类,并证明内波的能量总是向宇宙空间弥散掉,最终强度趋于零。因此,即使扰动的水平范围为无限,或者严格地考虑到地球表面为球面,也有地转适应过程,这是以在未研究过的机理。     

Applying a simple three-dimensional eddy correlation system for latent and sensible heat flux to contrasting forest canopies

Summary A simple eddy correlation system is presented that allows on-line calculation of latent and sensible heat fluxes. The system is composed of a three dimensional propeller anemometer, a thermocouple and a capacitance relative humidity sensor. Results from two contrasting sites demonstrate the capability of the system to measure turbulent fluxes under varying conditions. A dry mixed (dominantly coniferous) forest in hilly terrain in Austria is compared to a well irrigated, heavily transpiring, deciduous pecan orchard in the Southwest of the US. The US site shows insufficient closure of the energy balance that is attributed to non-turbulent fluxes under advective conditions in a stable boundary layer (Blanford et al., 1991) while the Austrian site exhibits almost perfect closure with the use of the very same instruments when the boundary layer is convective and advection is negligible.With 4 Figures     

Validating the turbulence parameterization schemes of a numerical model using eddy dissipation rate and turbulent kinetic energy measurements in terrain-disrupted airflow

A number of turbulence parameterization schemes are available in the latest version (6.0) of the Regional Atmospheric Modelling System (RAMS). Chan in Meteorol Atmos Phys 103:145–157, (2009), studied the performance of these schemes by simulating the eddy dissipation rate (EDR) distribution in the vicinity of the Hong Kong International Airport (HKIA) and comparing with the EDR measurements of remote-sensing instruments at the airport. For the e-l (turbulent kinetic energy ? mixing length) scheme considered in that study, the asymptotic mixing length was assumed to be a constant. This assumption is changed in the present paper, a variable asymptotic mixing length is chosen and simulations of EDR fields are repeated for terrain-disrupted airflow in the vicinity of HKIA. It is found that, with a variable asymptotic mixing length, the performance of the e-l scheme is greatly improved. With suitable choice of the empirical constants in the turbulence closure, the accuracy of the EDR profile (in comparison with LIDAR and wind profiler measurements) is found to be comparable with that predicted by the Deardorff scheme. A study on the sensitivity of the simulation results to these empirical constants has also been performed. Moreover, as a follow-up of the previous study of Chan in Meteorol Atmos Phys 103:145–157, (2009), case studies have been conducted on the following issues of the model simulation of turbulence for aviation application: (a) the effect of vertical gridding on the simulation results, (b) possibility of false alarm (such as over-forecasting of EDR value) in light turbulence cases, and (c) the performance in the simulation of other turbulence intensity metric for aviation purpose, e.g. TKE.     

A semi-analytical model of barotropic and baroclinic flows for an open Panama Gateway

A semi-analytical model of the Panama throughflow is presented. The model expresses the throughflow transport as a function of deep water formation in the North Pacific and in the North Atlantic, and of the Panama Gateway depth. The model is derived from the integral of the momentum equation along a circumpolar path, and can be interpreted from the point of view of the vorticity balance. The important conditions are whether the deep water, whose location is considered to be above the bottom water formed around Antarctica, originates from the North Atlantic or from the North Pacific, and whether the Panama Gateway is shallower than the lower boundary of the deep water. The present model indicates that the barotropic transport through the Panama Gateway is eastward, except for the case where the deep water is formed in the North Pacific and the sill of the Panama Gateway is shallow. The baroclinic structure of the Panama throughflow depends on whether the deep water is formed in the North Pacific or in the North Atlantic. These qualitative implications of the model are consistent with recent numerical studies and proxy-based paleoceanographic studies. Numerical experiments performed in the present study reinforce confidence in the semi-analytical model.     

A baroclinic planetary boundary-layer model,and its application to the Wangara data

Nondimensional parameters characteristic of the outer part of the planetary boundary layer have been determined by fitting a simple, Ekman-type theory to a number of averaged, observed velocity distributions, using the Wangara data of Clarke et al. (1971). The theoretical model is based on constant eddy viscosity in the outer layer and a linear variation of the geostrophic wind with height. At the lower boundary of the outer layer, the condition is applied that stress and velocity are parallel. This yields an equation for the cross-isobar angle as a function of drag coefficient, depth coefficient and nondimensional thermal wind.The data could be sorted into three well-defined, distinct groups, each characterized by a more or less constant value of the depth coefficient. The group with the lowest value of this parameter contains most of the nighttime data, the middle group the remaining nighttime data and most of the daytime ones, and the group with the largest depth, daytime data with cold air advection. The difference between the lowest and highest depth coefficients found here is about a factor of three.Within each group separately, the theoretically derived cross-isobar angle agrees remarkably well with the observed one, as a function of thermal wind.     

Effects of the baroclinic adjustment on the tropopause in the NCEP–NCAR reanalysis

In this work, we study the mean tropopause structure from the National Center for Environmental Prediction–National Center for Atmospheric Research reanalysis in the framework of baroclinic adjustment theories, focusing on the impact of baroclinic eddies on the mean tropopause height. In order to measure the effects of such perturbations, we introduce an appropriate global index that selects events of high baroclinic activity and allows us to distinguish the phases of growth and decay of baroclinic waves. We then composite the tropopause mean structure before and after baroclinic events, finding that baroclinic disturbances cause the zonally averaged midlatitude winter tropopause height to rise. Our results establish the importance of baroclinic adjustment processes for midlatitude tropopause dynamics.     

A numerical world ocean general circulation model: Part II. A baroclinic experiment

A new six-layer world ocean general circulation model based on the primitive system of equations is described in detail and its performance in the case of a homogeneous ocean is described. These test integrations show that the model is capable of reproducing the observed mean barotropic or vertically-integrated transport, as well as the seasonal variability of the major ocean gyres. The surface currents, however, are dominated by the Ekman transport, and such non-linear features as the western boundary currents and the equatorial countercurrents are poorly represented. The abyssal boundary countercurrents are also absent due to the lack of thermohaline forcing. The most conspicuous effect of the bottom topography on a homogeneous ocean is seen in the Southern ocean where the calculated Antarctic circumpolar transport through the Drake passage ( ≈ 10 Sv, with bathymetry included) greatly underestimates the observed transport (≈ 100 Sv).     

A one-dimensional photochemical model of the troposphere with planetary boundary-layer parameterization

The results from a one-dimensional photochemical model of the troposphere representative of summertime conditions at Northern Hemisphere mid-latitudes are presented. A parameterization of mixing processes within the planetary boundary layer (PBL) has been incorporated into the model for both the daytime convective PBL and the formation of the nocturnal PBL. One result of the parameterized PBL is that the concentrations of some trace species in the free troposphere are 20–30% higher than when mixing processes are described by a vertical eddy diffusion coefficient which is held constant with respect to height and time.The calculations indicate that the lifetime of the oxides of nitrogen (NO _x =NO+NO₂) against photochemical conversion to nitric acid (HNO₃) during summertime conditions is on the order of 6 h. This lifetime is short enough to deplete most of the NO _x in the PBL, resulting in the finding that other reactive nitrogen species (HNO₃ and peroxyacetyl nitrate) are more abundant than NO _x throughout the free troposphere, even though NO _x is the most abundant reactive nitrogen species at the surface. The effects of the inclusion of anthropogenic nonmethane hydrocarbon (NMHC) chemistry are also discussed. The inclusion of NMHC chemistry has a pronounced effect on the photochemistry of tropospheric oxone and increases thein situ column production by more than 30%.     

The subsurface transport of heat and moisture and its effect on the environment: A numerical model

A numerical model was developed to study the transport of heat and vapor under the surface of bare soil and soil covered by some materials such as asphalt and concrete under no rainfall conditions. The computational results provide a good match with the experimental data. The results show that the transport of water vapor inside the soil has an important effect on the subsurface distribution of temperature, especially for bare soil. Because of evaporation, the temperature of bare soil is much lower than that under covered surfaces throughout the day and the temperature of the surface covered by asphalt is extremely high-higher than the atmospheric temperature even at night. An increase of thickness of the covering material further increases the temperature and heat stored under surfaces. The stored heat is released to the atmosphere at night, contributing to environmental effects such as the urban heat island.     

A parametric model of vertical eddy fluxes in the atmosphere

A scheme for the representation of the vertical eddy fluxes of heat, momentum and water vapour in a forecast model is presented. An important feature of the scheme is the dependence of the diffusion coefficients on the static stability of the atmosphere. Two tests are presented, using the scheme in a one-dimensional model: the simulation of the diurnal cycle, and the transformation of a polar air mass moving over the warm sea.     

Convective momentum transport in a simple multicloud model for organized convection

从多层原始方程组出发,分析了不同尺度天气系统各层的动能水平梯度与该层地转偏差之间的关系,并对南京2010年7月22日暴雨过程的动能水平梯度模进行了诊断。结果表明：在天气尺度和α中尺度系统中,各层动能水平梯度模与该层地转偏差的模大致成正比。对天气尺度系统,该比例系数为地转参数,与该系统本身无关。对α中尺度系统,该比例系数与该系统本身有关。在β中尺度系统中,各层动能水平梯度模小于等于两项之和,而这两项中第一项与地转偏差的模成正比,另一项则为地转偏差时间导数的模。动能水平梯度大值区运动的非平衡性和爆发性强。在暴雨落区和飑线上,各层动能水平梯度的模均较大。通过对动能水平梯度模的诊断,能够从风场角度来分析各尺度系统,这对中尺度系统的诊断分析尤为重要。

     

Evaluation of six parameterization approaches for the ground heat flux

Summary There are numerous approaches to the parameterization of the ground heat flux that use different input data, are valid for different times of the day, and deliver results of different quality. Six of these approaches are tested in this study: three approaches calculating the ground heat flux from net radiation, one approach using the turbulent sensible heat flux, one simplified in situ measurement approach, and the force-restore method. On the basis of a data set recorded during the LITFASS-2003 experiment, the strengths and weaknesses of the approaches are assessed. The quality of the best approaches (simplified measurement and force-restore) approximates that of the measured data set. An approach calculating the ground heat flux from net radiation and the diurnal amplitude of the soil surface temperature also delivers satisfactory daytime results. The remaining approaches all have such serious drawbacks that they should only be applied with care. Altogether, this study demonstrates that ground heat flux parameterization has the potential to produce results matching measured ones very well, if all conditions and restrictions of the respective approaches are taken into account.