首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
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.  相似文献   

2.
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.  相似文献   

3.
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.  相似文献   

4.
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.  相似文献   

5.
Wind profile and eddy-correlation data obtained at two sites on a melting glacier surface in Iceland during the summer of 1996 are presented. Throughout the experiment the surface roughness increased rapidly from smooth to very rough, with the largest roughness element height obtained being about 1.7 m. In a layer close to the rough surface we find that the wind speed profiles were disturbed showing horizontal inhomogeneities as in a roughness sublayer. Its height was approximately two times the height of the main roughness elements (h) at both sites throughout the experiment. From the wind profiles and eddy-correlation data we calculated corrections for the displaced zero plane as a function of time and compared these with results obtained from a drag partitioning model. In general, the agreement was reasonable considering the ranges of uncertainty but the results indicate that the increasing horizontal anisotropy of the surface probably limits the use of the model. The values obtained for the roughness lengths are in good agreement with those calculated from a simple linear model, i.e., z0/h = 0.5 with the frontal area index. Above the roughness sublayer the wind profiles, normalised standard deviations of wind speed, and the balance of the turbulence kinetic energy budget behaved as over an ideal homogeneous surface thereby confirming similarity of the flow.  相似文献   

6.
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.  相似文献   

7.
Effects of Wall Heating on Flow Characteristics in a Street Canyon   总被引：5，自引：4，他引：1  
We develop a large-eddy simulation (LES) model based on a meteorological numerical model for a real scale street-canyon flow with rough building facets heated by a given temperature. The model is applied to a canyon with the aspect ratio of unity for two idealized heating scenarios: (1) the roof and the entire upstream wall are heated, named as ‘assisting cases’, and (2) the roof and the entire downstream wall are heated, named as ‘opposing cases’. These facets were heated up to 15 K above the air temperature. A wall function for temperature is proposed for a rough facet with an assumption that the thermal roughness length, z 0T, is much smaller than the aerodynamic roughness length, z 0. It is demonstrated that the sensible heat flux and canyon-air temperature are significantly influenced by the near-facet process that is parametrized by z 0T as the primary factor; other processes such as in-canyon mixing and roof-level exchange are secondary. This new finding strongly suggests that it is vital to choose an appropriate value of z 0T in a numerical simulation of street-canyon flows with the facet-air exchange processes of heat or any scalar. The finding also raises an awareness of the demand for carefully designed laboratory or field experiments of quantifying z 0T values for various urban surfaces. For the opposing cases, an unsteady penetrating narrow updraft zone appears occasionally along the heated wall and this feature is consistent field observations. The unique result indicates the superior capability of LES. The results of this study can be used to guide the parametrization of turbulent processes inside the urban canopy layer.  相似文献   

8.
Roughness length for heat over an urban canopy   总被引：1，自引：0，他引：1  
The roughness length for heat zT was evaluated over an urban canopy, using the measured sensible heat flux and radiometric temperature. To overcome thermal heterogeneity in the urban area, the measured radiometric temperature was transformed into the equivalent temperature of an upward longwave radiation flux. The equivalent temperature was found to provide an effective parameterization of the radiometric temperature. The daytime average of the resulting ln(zT/z0) was 10, where z0 is the aerodynamic roughness length. This result generally agrees with previous studies; however, the anthropogenic heat is a large uncertainty, which could cause an error at least 240% in zT.  相似文献   

9.
Abstract

A physically‐based numerical model was developed to estimate the temporal course of the surface energy flux densities and the soil temperatures in dry and wet bare soils. Aerodynamic heat, vapour and momentum transfer theory was used to calculate the sensible and latent heat flux densities at the surface under diabatic and adiabatic conditions. A finite‐difference solution of the differential equation describing one‐dimensional heat transfer was used to calculate the surface soil heat flux density and soil profile temperatures. The surface temperature was determined iteratively by the simultaneous solution of equations describing radiative, heat and momentum transfer at the surface. The model was tested with measurements from energy balance studies conducted on a dry, sandy soil and a wet, silt loam soil, and was found to predict accurately the surface energy fluxes and soil temperatures over three‐day periods under conditions of potential and negligible evaporation. The sensitivity of the model to uncertainties in the aerodynamic roughness lengths for momentum (z0) and heat (zT) is reported. Values for z0 and Z0/ZT of 0.5 mm and 3.0, respectively, resulted in the best agreement between modelled and measured values of the fluxes and temperatures for both soils.  相似文献   

10.
The roughness length for momentum (z0m), zero-plane displacementheight (d), and roughness length for heat (z0h) are importantparameters used to estimate land-atmosphere energy exchange. Although many different approaches have been developed to parameterizemomentum and heat transfer, existing parameterizations generally utilizehighly simplified representations of vegetation structure. Further, a mismatch exists between the treatments used for momentum and heat exchange and those used for radiative energy exchanges. In this paper, parameterizations are developed to estimate z0m, d, and z0h for forested regimes using information related to tree crown density and structure. The parameterizations provide realistic representationfor the vertical distribution of foliage within canopies, and include explicit treatment for the effects of the canopy roughness sublayer and leaf drag on momentum exchange. The proposed parameterizationsare able to realistically account for site-to-site differences in roughness lengths that arise from canopy structural properties.Comparisons between model predictions and field measurements show good agreement, suggesting that the proposed parameterizations capture the most important factors influencing turbulent exchange of momentumand heat over forests.  相似文献   

11.
A parameterisation ofthe effective roughness length is presented for an arbitrary givenroughness distribution z0(x,y) over flat terrainat neutral stratification.Beyond pure averaging, it takes into account the spatial structure of the distribution, especially the influence of length scales, and inflow direction.To allow for interactions between different rough areas, Boussinesq-approximated equations with a turbulence closure of first order are considered and solved using perturbation theory.As a result, the logarithm of the effective roughness length isrepresented as a sum over the product of the Fourier transformation of log z0 and a so-called dynamic function, which describesthe response of the flow field to a single wavelength of z0.Although the numerical expenditure is larger than for simple averaging formulae,this method could be used by large-scale models to calculateeffective roughness lengths in every close-to-surface grid cell.  相似文献   

12.
The surface flux calculation method proposed by Louis (1979) is extended by allowing momentum and heat to have different roughness lengths (z 0 andz T respectively). Our approach is to extend the traditional analysis by a more careful consideration of potential temperature structure near the surface. This consideration leads to a redefinition of the bulk Richardson number. For bulk transfer coefficients, our method shows a significant improvement over the original Louis method when compared with the theoretical surface-layer model applied to rough surfaces. Numerical experiments simulating sea breezes in 2D show that our extension is crucially important in simulating light wind and low humidity conditions.  相似文献   

13.
Estimation of Roughness Parameters Within Sparse Urban-Like Obstacle Arrays   总被引：1，自引：0，他引：1  
We conduct wind-tunnel experiments on three different uniform roughness arrays composed of sparsely distributed rectangular cylinders for the estimation of surface parameters. Roughness parameters such as the roughness length z 0 and zero-plane displacement d are extracted using a best-fit approximation of the measured wind velocity. We also perform a large-eddy simulation (LES) to confirm that four sampling points are sufficient to surrogate a space average above the canopy layer of the sparse roughness arrays. We propose a new morphological model from a systematic analysis of experimental data on the arrays. The friction velocity predicted by the proposed model agrees well with the peak value of the measured Reynolds shear stress ${(-\left<\overline{u'w'}\right>)^{0.5}}${(-\left<\overline{u'w'}\right>)^{0.5}}. The proposed model is further validated in an additional wind-tunnel experiment conducted on a scaled configuration of a real urban area exposed to four wind directions. The proposed model is found to perform very well particularly in the estimation of the friction velocity, readily leading to a better estimation of turbulence, which is essential for an accurate prediction of pollutant dispersion.  相似文献   

14.
Although the bulk aerodynamic transfer coefficients for sensible (C H ) and latent (C E ) heat over snow and sea ice surfaces are necessary for accurately modeling the surface energy budget, they have been measured rarely. This paper, therefore, presents a theoretical model that predicts neutral-stability values of C H and C E as functions of the wind speed and a surface roughness parameter. The crux of the model is establishing the interfacial sublayer profiles of the scalars, temperature and water vapor, over aerodynamically smooth and rough surfaces on the basis of a surface-renewal model in which turbulent eddies continually scour the surface, transferring scalar contaminants across the interface by molecular diffusion. Matching these interfacial sublayer profiles with the semi-logarithmic inertial sublayer profiles yields the roughness lengths for temperature and water vapor. When coupled with a model for the drag coefficient over snow and sea ice based on actual measurements, these roughness lengths lead to the transfer coefficients. C E is always a few percent larger than CH. Both decrease monotonically with increasing wind speed for speeds above 1 m s–1, and both increase at all wind speeds as the surface gets rougher. Both, nevertheless, are almost always between 1.0 × 10–3 and 1.5 × 10–3.  相似文献   

15.
A comparison between various parameterizations for the bulk transfercoefficient for heat and momentum is carried out for a wide range ofatmospheric stability and values of the roughness lengths for momentum and heat,z0m and z0h respectively.It is confirmed that the parameterization of Launiainen compares wellto a numerical iterative solution for the Obukhov length L as function of the gradients of wind speed and temperatureover a limited range of z0m/z0h and stability conditions.For 0m/z0h > 500, an alternativeinterpolation formulation of Holtslag and Ek,in combination with the formulation of Launiainen, provides a better approximation.  相似文献   

16.
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.  相似文献   

17.
A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z0) and heat (zT). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z0 that took intoaccount the form drag effect of sea-ice ridges together withzT = 0.1z0. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity Kmin. A simple ad hocparameterization seems applicable, where Kmin is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.  相似文献   

18.
不同形状冰晶权重假定对冰云光学和辐射特性的影响   总被引：1，自引：0，他引：1  
陈琪  张华 《气象学报》2018,76(2):279-288
在BCC_RAD辐射传输模式和包含多形状冰晶粒子的冰云光学性质参数化方案的基础上,详细分析了不同冰晶粒子权重选取对冰云光学和辐射特性的影响。结果显示,不同形状冰晶粒子权重的选取对长波带平均消光系数、单次散射比、不对称因子和短波带平均不对称因子均有较大的影响。冰晶粒子权重选取对长波辐射通量有很大影响:对长波向下辐射通量,权重选择不同可在云底处造成高达10.50 W/m2的差别;对长波向上辐射通量,权重选择不同可在云顶处造成高达15.05 W/m2的差别。冰晶粒子权重选择对短波辐射通量也存在较大影响:对短波向下辐射通量,权重选择不同可在云底处造成高达12.48 W/m2的差别;对短波向上辐射通量,权重选择不同可在云顶处造成高达10.23 W/m2的差别。冰晶粒子权重选择对长波加热率影响较大,在云顶处和云底处分别可达1.31和-2.06 K/d。研究表明,不同形状冰晶粒子权重的选取对冰云光学性质和辐射计算均有较大的影响,在长波区间尤其明显。   相似文献   

19.
The roughness lengths for momentum and temperature are calculated using the profile method on amelting glacier surface. Data from a 5-level 9-m meteorological mastpositioned near the edge of Breidamerkurjökull, an outlet glacier of the Vatnajökull ice cap Iceland, are used for the calculations. The data are selected to avoid the presence of the katabatic wind speedmaximum which would otherwise alter the scaling laws of the surface layer. The surface roughness length for momentum is determined to be 1.0 mm, similar to other estimates made on flat melting ice surfaces. The surface roughness length for temperature is found to be in good agreement with previously proposed surface renewal theories for the observed roughness Reynolds number range of 30 * 70.  相似文献   

20.
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.  相似文献   

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

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