首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
The Louis scheme for calculating the vertical eddy fluxes within the atmospheric surface layer is improved by broadening the original assumptions. In our approach, the momentum and heat transfer roughness lengths (z0 and zT respectively) can be different, and z0 need not be negligibly small compared with the lowest height (z) in modelling. For these conditions, we choose more consistent wind and potential temperature profile forms, then derive new algorithms for calculating fluxes. Improvement is demonstrated for a wide range of z/L (L is the Obukhov length), z/z0 and z0 zT, by comparing these fluxes with those derived from a theoretical surface-layer model. The improved algorithms can be used in atmospheric modelling systems for more varied surfaces and a wide range of atmospheric stability.  相似文献   

2.
To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z 0) is a key parameter. We present, for the first time, spatial and temporal variations of z 0 in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z 0 is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z 0≈ 10−4 m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z 0 quickly increases well above 10−2 m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z 0 > 5  ×  10−3 m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z 0 = 10−5 m. At the beginning of winter, we observed that single snow events have the potential to lower z 0 for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.  相似文献   

3.
We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z 0). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z 0 below a threshold value of approximately 10?3 m, the model performs well and in accord with earlier studies. However, with growing hummock size, z 0 increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z 0 = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z 0  >  10?3 m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z s /z 0)  =  1.5  ? 0.2 ln (Re *)  ? 0.11(ln (Re *))2.  相似文献   

4.
采用北京325 m铁塔2008—2012年的单层超声观测资料,基于莫宁-奥布霍夫相似理论(Monin-Obukhov similarity theory)和前人提出的最小误差分析方法,计算了铁塔周边下垫面的零平面位移高度和动力粗糙度长度。结果表明,由于铁塔位于北京市区,其周边下垫面呈现极其复杂的非均匀性,所以对应铁塔周边不同的扇区,零平面位移高度和动力粗糙度长度各有不同。平均而言,在2008—2012年间,铁塔周边下垫面的零平面位移高度为34.4 m,动力粗糙度长度为1.16 m。此外,综合前人的计算结果发现,铁塔周边的零平面位移高度和动力粗糙度长度在2001年之前呈显著增加的趋势,而在2001年以后并未增长,这一现象与铁塔周边的城市化进程相对应。  相似文献   

5.
This paper examines the practical importance of stability, baroclinicity, and acceleration in the bulk ABL similarity formulations, in light of the random errors inherent in field measurements. This is done by propagating the measurement uncertainties through a theoretical model for the bulk ABL similarity functionsA 0 andB 0, under a range of assumed (but always unstable) conditions. It is shown that random measurement errors and acceleration effects may overwhelm most effects of baroclinicity and stability, once conditions are at least slightly unstable. Because of this, it is hard to discern a clear functional dependence ofA 0 andB 0 on degree of instability. Thus, for a given value ofh i/z0, whereh i is the inversion height andz 0 is the surface roughness, the geostrophic drag coefficient, which depends onA 0 andB 0, and weakly onh i/z0, may also be taken to be nearly independent of degree of instability.  相似文献   

6.
Mean wind velocity profiles were measured by means of radio-windsondes over the Landes region in southwestern France, which consists primarily of pine forests with scattered villages and clearings with various crops. Analysis of neutral profiles indicated the existence of a logarithmic layer between approximately z – d 0 = 67(±18)z 0 and 128(+-32)z 0 (z is the height above the ground, z 0 the surface roughness and d 0 the displacement height). The upper limit can also be given as z – d 0 = 0.33 (±0.18)h, where h is the height of the bottom of the inversion. The profiles showed that the surface roughness of this terrain is around 1.2 m and the displacement height 6.0 m. Shear stresses derived from the profiles were in good agreement with those obtained just above the forest canopy at a nearby location with the eddy correlation method by a team from the Institute of Hydrology (Wallingford, England).  相似文献   

7.
We present a field investigation over a melting valley glacier on the Tibetan Plateau. In the ablation zone, aerodynamic roughness lengths (z 0M ) vary on the order of 10−4–10−2 m, whose evolution corresponds to three melt phases with distinct surface cover and moisture exchange: snow (sublimation/evaporation), bare ice (deposition/condensation), and ice hummocks (sublimation/evaporation). Bowen-ratio similarity is validated in the stably stratified katabatic winds, which suggests a useful means for data quality check. A roughness sublayer is regarded as irrelevant to the present ablation season, because selected characteristics of scalar turbulence over smooth snow are quite similar to those over hummocky ice. We evaluate three parametrizations of the scalar roughness lengths (z 0T for temperature and z 0q for humidity), viz. key factors for the accurate estimation of sensible heat and latent heat fluxes using the bulk aerodynamic method. The first approach is based on surface-renewal models and has been widely applied in glaciated areas; the second has never received application over an ice/snow surface, despite its validity in (semi-)arid regions; the third, a derivative of the first, is proposed for use specifically over rough ice defined as z 0M > 10−3 m or so. This empirical z 0M threshold value is deemed of general relevance to glaciated areas (e.g. ice sheet/cap and valley/outlet glaciers), above which the first approach gives notably underestimated z 0T,q . The first and the third approaches tend to underestimate and overestimate turbulent heat/moisture exchange, respectively, frequently leading to relative errors higher than 30%. Comparatively, the second approach produces fairly low errors in energy flux estimates both in individual melt phases and over the whole ablation season; it thus emerges as a practically useful choice to parametrize z 0T,q in glaciated areas. Moreover, we find all three candidate parametrizations unable to predict diurnal variations in the excess resistances to humidity transfer, thus encouraging more efforts for improvement.  相似文献   

8.
By non-dimensionalizing a trajectory-simulation (TS) model of turbulent dispersion, it is shown that the dimensionless concentration z 0cu*/kQ (cu */kQ) due to a continuous line (area) source of strength Q in the atmospheric surface layer depends only on z/z 0, x/z 0, z 0/L and z s/z0, where z s is the source height. The TS model is used to tabulate concentration profiles due to ground-level line and area sources. Concentration profiles generated by the TS model for elevated sources are shown to be inconsistent with the Reciprocal Theorems of Smith (1957) and it is suggested that this is because the flux-mean gradient closure scheme inherent in the Reciprocal Theorem is invalid for an elevated source.  相似文献   

9.
Improving and validating land surface models based on integrated observations in deserts is one of the challenges in land modeling. Particularly, key parameters and parameterization schemes in desert regions need to be evaluated in-situ to improve the models. In this study, we calibrated the land-surface key parameters and evaluated several formulations or schemes for thermal roughness length (z 0h ) in the common land model (CoLM). Our parameter calibration and scheme evaluation were based on the observed data during a torrid summer (29 July to 11 September 2009) over the Taklimakan Desert hinterland. First, the importance of the key parameters in the experiment was evaluated based on their physics principles and the significance of these key parameters were further validated using sensitivity test. Second, difference schemes (or physics-based formulas) of z 0h were adopted to simulate the variations of energy-related variables (e.g., sensible heat flux and surface skin temperature) and the simulated variations were then compared with the observed data. Third, the z 0h scheme that performed best (i.e., Y07) was then selected to replace the defaulted one (i.e., Z98); the revised scheme and the superiority of Y07 over Z98 was further demonstrated by comparing the simulated results with the observed data. Admittedly, the revised model did a relatively poor job of simulating the diurnal variations of surface soil heat flux, and nighttime soil temperature was also underestimated, calling for further improvement of the model for desert regions.  相似文献   

10.
The spray content in the surface boundary layer above an air—water interface was determined by a series of measurements at various feteches and wind speeds in a laboratory facility. The droplet flux density N(z) can be described in terms of the scaling flux density N* and von Karman constant K throguh the equation, N(z)/N* = −(1/K) ln(z/z0d) where z is height above the mean water level and z0d is the droplet boundary layer thickness. N* is given by a unique relationship in terms of the roughness Reynolds number u*σ/ν where σ is the root-mean-square surface displacement. Spray inception occurred for u* 0.3. The dominant mode of spray generation in the present and most other laboratory tests, as well as in available field data, appears to be bubble bursting.  相似文献   

11.
Vertical profiles of wind speed, temperature and humidity were used to estimate the roughness lengths for momentum (z 0), heat (z H ) and moisture (z Q) over smooth ice and snow surfaces. The profile-measurements were performed in the vicinity of a blue ice field in Queen Maud Land, East Antarctica. The values ofz 0 over ice (3·10–6 m) seem to be the smallest ever obtained over permanent, natural surfaces. The settling of snow on the ice and the loss of momentum at saltating snow particles serve as momentum dissipating processes during snow-drift events, expressed as a strong dependence ofz 0 on u#.The scalar roughness lengths and surface temperature can be evaluated from the temperature and humidity profile measurements if the ratioz H /z Q is specified. This new method circumvents the difficult measurement of surface temperature. The scalar roughness lengths seem to be approximately equal toz0 for a large range of low roughness Reynolds numbers, despite the frequent occurrence of drifting snow. Possible reasons for this agreement with theory of non-saltating flow are discussed.  相似文献   

12.
A differential equation is obtained to describe the concentration of passive admixtures (water vapor, sensible heat, pollutants, CO2, etc.) of turbulent flow inside a dense and uniform vegetational canopy. The profiles of eddy diffusivity, wind speed and shear stress are assumed to be exponential decay functions of depth below the top of the canopy. This equation is solved for the case of a vegetation with constant concentration of the admixture at the foliage surfaces. The solution is used to formulate bulk mass or heat transfer coefficients, which can be applied to practical problems involving surfaces covered with a vegetation or with similar porous or fibrous roughness elements. The results are shown to be consistent with experimental data presented by Chamberlain (1966), Garratt and Hicks (1973) and Garratt (1978). Calculations with the model illustrate that, as compared to its behavior over surfaces with bluff roughness elements, ln(z 0/z 0c ) (wherez 0 is the momentum roughness andz 0c , the scalar roughness) for permeable roughness elements is relatively insensitive tou * and practically independent ofz 0.  相似文献   

13.
Using the relationship between the bulk Richardson numberR z and the Obukhov stability parameterz/L (L is the Obukhov length), formally obtained from the flux-profile relationships, methods to estimatez/L are discussed. Generally,z/L can not be uniquely solved analytically from flux-profile relationships, and it may be defined using routine observations only by iteration. In this paper, relationships ofz/L in terms ofR z obtained semianalytically were corrected for variable aerodynamic roughnessz 0 and for aerodynamic-to-temperature roughness ratiosz 0/z T, using the flux-profile iteration procedure. Assuming the so-called log-linear profiles to be valid for the nearneutral and moderately stable region (z/L<1), a simple relationship is obtained. For the extension to strong stability, a simple series expansion, based on utilisation of specified universal functions, is derived.For the unstable region, a simple form based on utilisation of the Businger-Dyer type universal functions, is derived. The formulae yield good estimates for surfaces having an aerodynamic roughness of 10–5 to 10–1 m, and an aerodynamic-to-temperature roughness ratio ofz 0/z T=0.5 to 7.3. When applied to the universal functions, the formulae yield transfer coefficients and fluxes which are almost identical with those from the iteration procedure.  相似文献   

14.
Temperature variance and temperature power spectra in the unstable surface layer have always presented a problem to the standard Monin-Obukhov similarity model. Recently that problem has intensified with the demonstration by Smedman et al. (2007, Q J Roy Meteorol Soc 133: 37–51) that temperature spectra and heat-flux cospectra can have two distinct peaks in slightly unstable conditions, and by McNaughton et al. (2007, Nonlinear Process Geophys 14: 257–271) who showed that the wavenumber of the peak of temperature spectra in a convective boundary layer (CBL), closely above the surface friction layer (SFL), can be sensitive to the CBL depth, z i. Neither the two-peak form at slight instability nor the dependence of peak position on z i at large instability is compatible with the Monin-Obukhov model. Here we examine the properties of temperature spectra and heat-flux cospectra from between these extremes, i.e. from within the unstable SFL, in two experiments. The analysis is based on McNaughton’s model of the turbulence structure in the SFL. According to this model, heat is transported through most of the SFL by sheet plumes, created by the action of impinging outer eddies. The smallest and most effective of these outer eddies have sizes that scale on SFL depth, z s. The z s-scale eddies and plumes are organised within the overall convection pattern in the CBL, and in turn they organise the motion of smaller eddies within the SFL, whose sizes scale on height, z. The main experimental results are: (1) the peak amplitudes of the temperature spectra in the SFL are collapsed with a scaling factor (zsz)1/3eo2/3{(z_{\rm s}z)^{1/3}\varepsilon_{\rm o}^{2/3}} divided by the square of the surface temperature flux, where eo{\varepsilon_{\rm o}} is the dissipation rate of turbulent energy in the outer CBL (above the SFL); (2) the peak wavenumbers of the temperature spectra are collapsed with the mixed length scale (z i z s)1/2; (3) the peak wavenumbers of the heat-flux cospectra are collapsed with the doubly-mixed length scale (z i z s)1/4 z 1/2; (4) for z/z s < 0.03, the peak in the cospectrum is replaced by another peak at a wavenumber about a magnitude larger. This peak’s position scales on z; (5) all these findings are consistent with the observations of Smedman et al.  相似文献   

15.
It is proposed that the ratios of the standard deviations of the horizontal velocity components to the friction velocity in the surface layer under convective conditions depend only onz i /L wherez i is the height of the lowest inversion andL is the Monin-Obukhov length. This hypothesis is tested by using observations from several data sets over uniform surfaces and appears to fit the data well. Empirical curves are fitted to the observations which have the property that at largez i /-L, the standard deviations become proportional tow *, the convective scaling velocity.Fluctuations of vertical velocity obtained from the same experiments scale withz/L, wherez is the height above the surface, in good agreement with Monin-Obukhov theory.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

16.
An approximate method for calculating the relationship between z/L(z = reference height, L = Obukhov length) and the bulk Richardsonnumber is presented. If this relationship is known, the momentum andheat fluxes can be computed easily without any iteration. The avoidance of iteration can speed up computationsin large-scale models considerably (up to 10 times) and cases which do not converge or converge very slowly cannot occur. The proposed formulae take into account the difference between momentum (z0M) and heat roughnesslengths (z0H). Because the roughness lengths are not neglected at any step of the derivation, the resulting analytical formulae can be used not only between the surface and the reference height but also between two finite levels z1 andz2 (by replacing z0M and z0H by z1 and z by z2). Theequations remain correct even in the limit z1 z2.The formulae are based upon the (partially modified) Businger–Dyer flux–profile relationships and,consequently, they are restricted to predominantly homogeneous terrain.These new approximations are an improvement over the existing solutions because they are simpler than most of the formulae in the literature and are able to match the numerical exact solution for different parameter sets (Businger, Dyer, Högström) with an maximum error of about 2% for a wide range of z/L, z/z0M and z0M/z0H.Furthermore, in stable conditions, schemes with and without a finitecritical bulk Richardson number can be approximated. The possibleambiguity of the exact solution =f(RIB) in (moderately) stable conditions is discussed briefly. The performance of the new formulae is compared to the exact numerical solution and to different formulae proposed in the literature.  相似文献   

17.
An easily-evaluated expression for the dimensionless concentration profile (z/z 0,/z 0, z 0/L) = = cu */kQ (or z 0cu*/kQ) downwind of a continuous ground-level area (or line) source in the stable surface layer is obtained by integrating the diffusion equation using the Shwetz approximation method (c = concentration, Q = source strength, k = von Kárman's constant). The analytical solution compares closely with concentration profiles obtained using a trajectory-simulation model over a useful range of heights, the important discrepancies occurring at the upper edge of the plume. The analytical solution is used to generate predictions of ground-level concentration for the Project Prairie Grass experiments; good agreement with the observations is obtained at all downwind distances (50 to 800 m).  相似文献   

18.
In the framework of an international field program for the study of semi-arid areas, observations were done in the region called La Crau in southern France. In this paper, the use of the surface radiative temperature for the determination of the sensible heat flux is addressed. We found that, once proper values of the roughness length of momentum (z 0) and heat (z 0h) are set, the sensible heat flux can be reliably predicted with a one-layer resistance model using standard observations of wind speed and air temperature, together with the surface temperature. The latter quantity has to be known with a precision better than ±2°C. From our observations, the value of the parameterB –1k –1 In (z 0 z 0h) was found to be 9.2, which falls between values quoted by Brutsaert (1982) for grass and bluff bodies.  相似文献   

19.
A simple new model is proposed to predict the distribution of wind velocity and surface shear stress downwind of a rough-to-smooth surface transition. The wind velocity is estimated as a weighted average between two limiting logarithmic profiles: the first log law, which is recovered above the internal boundary-layer height, corresponds to the upwind velocity profile; the second log law is adjusted to the downwind aerodynamic roughness and local surface shear stress, and it is recovered near the surface, in the equilibrium sublayer. The proposed non-linear form of the weighting factor is equal to ln(z/z 01)/ln(δ i /z 01), where z, δ i and z 01 are the elevation of the prediction location, the internal boundary-layer height at that downwind distance, and the upwind surface roughness, respectively. Unlike other simple analytical models, the new model does not rely on the assumption of a constant or linear distribution for the turbulent shear stress within the internal boundary layer. The performance of the new model is tested with wind-tunnel measurements and also with the field data of Bradley. Compared with other existing analytical models, the proposed model shows improved predictions of both surface shear stress and velocity distributions at different positions downwind of the transition.  相似文献   

20.
The flux footprint is the contribution, per unit emission, of each element of a surface area source to the vertical scalar flux measured at height z m ; it is equal to the vertical flux from a unit surface point source. The dependence of the flux footprint on crosswind location is shown to be identical to the crosswind concentration distribution for a unit surface point source; an analytic dispersion model is used to estimate the crosswind-integrated flux footprint. Based on the analytic dispersion model, a normalized crosswind-integrated footprint is proposed that principally depends on the single variable z/z m , where z is a measure of vertical dispersion from a surface source. The explicit dependence of the crosswind-integrated flux footprint on downwind distance, thermal stability and surface roughness is contained in the dependence of z on these variables. By also calculating the flux footprint with a Lagrangian stochastic dispersion model, it is shown that the normalized flux footprint is insensitive to the analytic model assumption of a self-similar vertical concentration profile.The National Center for Atmospheric Research is funded by the National Science Foundation.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号