首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 109 毫秒
1.
Summary This paper investigates the influence of the planetary boundary-layer (PBL) parameterization and the vertical distribution of model layers on simulations of an Alpine foehn case that was observed during the Mesoscale Alpine Programme (MAP) in autumn 1999. The study is based on the PSU/NCAR MM5 modelling system and combines five different PBL schemes with three model layer settings, which mainly differ in the height above ground of the lowest model level (z 1). Specifically, z 1 takes values of about 7 m, 22 m and 36 m, and the experiments with z 1 = 7 m are set up such that the second model level is located at z = 36 m. To assess if the different model setups have a systematic impact on the model performance, the simulation results are compared against wind lidar, radiosonde and surface measurements gathered along the Austrian Wipp Valley. Moreover, the dependence of the simulated wind and temperature fields at a given height (36 m above ground) on z 1 is examined for several different regions. Our validation results show that at least over the Wipp Valley, the dependence of the model skill on z 1 tends to be larger and more systematic than the impact of the PBL scheme. The agreement of the simulated wind field with observations tends to benefit from moving the lowest model layer closer to the ground, which appears to be related to the dependence of lee-side flow separation on z 1. However, the simulated 2 m-temperatures are closest to observations for the intermediate z 1 of 22 m. This is mainly related to the fact that the simulated low-level temperatures decrease systematically with decreasing z 1 for all PBL schemes, turning a positive bias at z 1 = 36 m into a negative bias at z 1 = 7 m. The systematic z 1-dependence is also observed for the temperatures at a fixed height of 36 m, indicating a deficiency in the self-consistency of the model results that is not related to a specific PBL formulation. Possible reasons for this deficiency are discussed in the paper. On the other hand, a systematic z 1-dependence of the 36-m wind speed is encountered only for one out of the five PBL schemes. This turns out to be related to an unrealistic profile of the vertical mixing coefficient. Correspondence: Günther Z?ngl, Meteorologisches Institut der Universitat München, 80333 München, Germany  相似文献   

2.
A wind-tunnel experiment has been used to investigate momentum absorption by rough surfaces with sparse random and clustered distributions of roughness elements. An unusual (though longstanding) method was used to measure the boundary-layer depth δ and friction velocity u * and thence to infer the functional relationship z 0/h = f(λ) between the normalised roughness length z 0/ h and the roughness density λ (where z 0 is the roughness length and h the mean height of the roughness elements). The method for finding u * is based on fitting the velocity defect in the outer layer to a functional form for the dimensionless velocity-defect profile in a canonical zero-pressure-gradient boundary layer. For the conditions investigated here, involving boundary layers over sparse roughness with strong local heterogeneity, this velocity-defect-law method is found to be more robust than several alternative methods for finding u * (uw covariance, momentum integral and slope of the logarithmic velocity profile).The experimental results show that, (1) there is general agreement in the relationship z 0/h = f(λ) between the present experiment with random arrays and other wind-tunnel experiments with regular arrays; (2) the main effect of clustering is to increase the scatter in the z 0/h = f(λ) relationship, through increased local horizontal heterogeneity; (3) this scatter obscures any trend in the z 0/h = f(λ) relationship in response to clustering; and (4) the agreement between the body of wind-tunnel data (taken as a whole) and field data is good, though with scatter for which it is likely that a major contribution stems from local horizontal heterogeneity in the field.  相似文献   

3.
Summary In this paper, we evaluate the applicability of flux-gradient relationships for momentum and heat for urban boundary layers within the Monin-Obukhov similarity (MOS) theory framework. Although the theory is widely used for smooth wall boundary layers, it is not known how well the theory works for urban layers. To address this problem, we measured the vertical profiles of wind velocity, air temperature, and fluxes of heat and momentum over a residential area and compared the results to theory. The measurements were done above an urban canopy whose mean height zh is 7.3 m. 3-D sonic anemometers and fine wire thermocouples were installed at 4 heights in the region 1.5zh < z < 4zh. We found the following: (1) The non-dimensional horizontal wind speed has good agreement with the stratified logarithmic profile predicted using the semi-empirical Monin-Obukov similarity (MOS) function, when it was scaled by the surface friction velocity that is derived from the shear stress extrapolated to the roof-top level. (2) The scaled gradient of horizontal wind speed followed a conventional semi-empirical function for a flat surface at a level (z/zh = 2.9), whereas, in the vicinity of the canopy height was larger than the commonly-used empirical relationship. (3) The potential temperature profile above the canopy shows dependency on the atmospheric stability and the scaled gradient of temperature is in good agreement with a conventional shear function for heat. In the case of heat, the dependency on height was not found. (4) The flux-gradient relationship for momentum and heat in the region 1.5zh < z < 4zh was rather similar to that for flat surfaces than that for vegetated canopies.  相似文献   

4.
Flux-profile relationships based on surface-layer similarity theory are used to derive relationships between the Monin-Obukhov stability parameter = z/L and the bulk Richardson number Ri b . In contrast to previous studies, the roughness length for heat, z 0h ,is assumed unequal to the roughness length for momentum, z 0m .For the stable case, an analytic expression of in terms of Ri b can be derived and in the unstable case, the solution is obtained through a simple iterative process.Errors introduced from the simplification of z 0h = z 0m are evaluated and are shown to be very significant in most cases. Thus, this error in many practical applications may invalidate the intended solution.  相似文献   

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

6.
It is well known that if turbulent mass convection is modelled as diffusion, errors result unless trajectories from the source (ath) to the point of observation (z p ) comprise many statistically-independent segments (Taylor, 1921). We show that this is not guaranteed merely by the Lagrangian timescale () at the source being small (e.g., source at ground), but that a better criterion istmax[(h), (z p )], wheret is a typical travel time toz p .  相似文献   

7.
The roughness height z 0 and the zero-plane displacement height d 0 were determined for a region of complex terrain in the Pre-Alps of Switzerland. This region is characterized by hills of the order of 100 m above the valley elevations, and by distances between ridges of the order of 1 km; it lies about 20 to 30 km north from the Alps. The experimental data were obtained from radiosonde observations under near neutral conditions. The analysis was based on the assumption of a logarithmic profile for the mean horizontal wind existing over one half of the boundary layer. The resulting (z 0/h) and (d 0/h) (where h is the mean height of the obstacles) were found to be in reasonable agreement with available relationships in terms of placement density and shape factor of the obstacles, which were obtained in previous experiments with h-scales 2 to 4 orders of magnitude smaller than the present ones.  相似文献   

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

9.
Impacts of different closure schemes in the Mixed Spectral Finite-Difference model (Beljaarset al., 1987) for neutrally stratified atmospheric surface-layer flow over complex terrain are studied. Six different closure schemes, (Z+z 0), Mixing Length,E–(Z+z 0),E–,E–– andq 2 l are compared. Model results for flow over an infinite series of sinusoidal ridges are examined in the context of the inner and outer layers defined by Jackson and Hunt (1975). Results are compared with rapid distortion estimates of the changes in normal stresses. The effects of streamline curvature are also examined in a qualitative sense.  相似文献   

10.
Summary In this paper the results of an urban measurement campaign are presented. The experiment took place from July 1995 to February 1996 in Basel, Switzerland. A total of more than 2000 undisturbed 30-minute runs of simultaneous measurements of the fluctuations of the wind vector u′, v′, w′ and the sonic temperature θ s ′ at three different heights (z=36, 50 and 76 m a.g.l.) are analysed with respect to the integral statistics and their spectral behaviour. Estimates of the zero plane displacement height d calculated by the temperature variance method yield a value of 22 m for the two lower levels, which corresponds to 0.92 h (the mean height of the roughness elements). At all three measurement heights the dimensionless standard deviation σ w /u * is systematically smaller than the Monin-Obukhov similarity function for the inertial sublayer, however, deviations are smaller compared to other urban turbulence studies. The σθ/θ* values follow the inertial sublayer prediction very close for the two lowest levels, while at the uppermost level significant deviations are observed. Profiles of normalized velocity and temperature variances show a clear dependence on stability. The profile of friction velocity u * is similar to the profiles reported in other urban studies with a maximum around z/h=2.1. Spectral characteristics of the wind components in general show a clear dependence on stability and dimensionless measurement height z/h with a shift of the spectral peak to lower frequencies as thermal stability changes from stable to unstable conditions and as z/h decreases. Velocity spectra follow the −2/3 slope in the inertial subrange region and the ratios of spectral energy densities S w (f)/S u (f) approach the value of 4/3 required for local isotropy in the inertial subrange. Velocity spectra and spectral peaks fit best to the well established surface layer spectra from Kaimal et al. (1972) at the uppermost level at z/h=3.2. Received September 26, 1997 Revised February 15, 1998  相似文献   

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

12.
Aerodynamic Scaling for Estimating the Mean Height of Dense Canopies   总被引：1，自引：1，他引：0  
We used an aerodynamic method to objectively determine a representative canopy height, using standard meteorological measurements. The canopy height may change if the tree height is used to represent the actual canopy, but little work to date has focused on creating a standard for determining the representative canopy height. Here we propose the ‘aerodynamic canopy height’ h a as the most effective means of resolving the representative canopy height for all forests. We determined h a by simple linear regression between zero-plane displacement d and roughness length z 0, without the need for stand inventory data. The applicability of h a was confirmed in five different forests, including a forest with a complex canopy structure. Comparison with stand inventory data showed that h a was almost equivalent to the representative height of trees composing the crown surface if the forest had a simple structure, or to the representative height of taller trees composing the upper canopy in forests with a complex canopy structure. The linear relationship between d and z 0 was explained by assuming that the logarithmic wind profile above the canopy and the exponential wind profile within the canopy were continuous and smooth at canopy height. This was supported by observations, which showed that h a was essentially the same as the height defined by the inflection point of the vertical profile of wind speed. The applicability of h a was also verified using data from several previous studies.  相似文献   

13.
The iteration scheme of Robinson (1962) for the determination ofz 0 andd under neutral atmospheric stability has been generalized in a relatively simple manner. It now accounts for nonneutral conditions and permits simultaneous calculation ofd,z 0, and flux densities of latent and sensible heat. Using profile data sets given in the literature, the generalized scheme yields results which agree quite well with those obtained by Kramm (1989).  相似文献   

14.
The steady-state, homogeneous and barotropic equations of motion within the planetary boundary layer are solved with the assumption that the coefficient of eddy viscosity varies as K(Z) = K O(1–Z/h) p , where h is the height of the bounday layer and p is a parameter which depends on atmospheric stability. The solutions compare favourably with observed velocity profiles based on the Wangara data.  相似文献   

15.
It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u */w *, where u * and w * are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z i , where z i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.  相似文献   

16.
This work presents the results of a comparison between the global storm activity rate IRS and electric field intensity E0Z. The permanent analysis of the IRS may become an important tool for testing Global Electric Circuit models. IRS is determined by a new method that uses the background component of the first 7 Schumann resonances (SR). The rate calculations are based on ELF observations carried out in 2005 and 2006 in the observatory station “Hylaty” of the Jagiellonian University in the Eastern Carpathians (Kułak, A., Zięba, S., Micek, S., Nieckarz, Z., 2003. Solar variations in extremely low frequency propagation parameters: I. A two-dimensional telegraph equation (TDTE) model of ELF propagation and fundamental parameters of Schumann resonances, J. Geophys. Res., 108, 1270, doi:10.1029/2002JA009304). Diurnal runs of the IRS rate were compared with diurnal runs of E0Z amplitudes registered at the Earth's surface in the Geophysical Observatory of the Polish Academy of Sciences in Świder (Kubicki, M., 2005. Results of Atmospheric Electricity and Meteorological Observations, S. Kalinowski Geophysical Observatory at Świder 2004, Pub. Inst. Geophysics Polish Academy of Sciences, D-68 (383), Warszawa.). The days with the highest values of the correlation coefficient (R) between amplitudes of both observed parameters characterizing atmosphere electric activity are shown. The seasonal changes of R, IRS and E0Z are also presented.  相似文献   

17.
Measurements of atmospheric turbulence made during the Surface Heat Budget of the Arctic Ocean Experiment (SHEBA) are used to examine the profile stability functions of momentum, φ m , and sensible heat, φ h , in the stably stratified boundary layer over the Arctic pack ice. Turbulent fluxes and mean meteorological data that cover different surface conditions and a wide range of stability conditions were continuously measured and reported hourly at five levels on a 20-m main tower for 11 months. The comprehensive dataset collected during SHEBA allows studying φ m and φ h in detail and includes ample data for the very stable case. New parameterizations for φ m (ζ) and φ h (ζ) in stable conditions are proposed to describe the SHEBA data; these cover the entire range of the stability parameter ζ = z/L from neutral to very stable conditions, where L is the Obukhov length and z is the measurement height. In the limit of very strong stability, φ m follows a ζ 1/3 dependence, whereas φ h initially increases with increasing ζ, reaches a maximum at ζ ≈ 10, and then tends to level off with increasing ζ. The effects of self-correlation, which occur in plots of φ m and φ h versus ζ, are reduced by using an independent bin-averaging method instead of conventional averaging.  相似文献   

18.
A method based on Giorgi (1997a, 1997b) and referred to as ‘combined approach’, which is a combination of mosaic approach and analytical-statistical-dynamical approach, is proposed. Compared with those of other approaches, the main advantage of the combined approach is that it not only can represent both interpatch and intrapatch variability, but also cost less computational time when the land surface heterogeneity is considered. Because the independent variable of probability density function (PDF) is extended to the single valued function of basic meteorological characteristic quantities, which is much more universal, the analytical expressions of the characteristic quantities, (e.g., drag coefficient, snow coverage, leaf surface aerodynamical resistance) affected by roughness length are derived, when the roughness length (and/ or the zero plane displacement) heterogeneity has been mainly taken into account with the approach. On the basis of the rule which the PDF parameters should follow, we choose a function y of the roughness length z0 as the PDF independent variable, and set different values of the two parameters width ratio αn and height ratio γ of PDF (here a linear, symmetric PDF is applied) for sensitivity experiments, from which some conclusions can be drawn, e.g., relevant characteristic terms show different sensitivities to the heterogeneous characteristic (i.e., roughness length), which suggests that we should consider the heterogeneities of the more sensitive terms in our model instead of the heterogeneities of the rest, and which also implies that when the land surface scheme is coupled into the global or regional atmospheric model, sensitivity tests against the distribution of the heterogeneous characteristic are very necessary; when the parameterαn is close to zero, little heterogeneity is represented, andαn differs with cases, which have an upper limit of about 0.6; in the reasonable range ofαn, a peak-like distribution of roughness length can be depicted by a small value ofγ, etc..  相似文献   

19.
Momentum and Heat Transfer over Urban-like Surfaces   总被引：2，自引：2，他引：0  
Momentum and heat transfer was examined for the urban-like surfaces used within the Comprehensive Outdoor Scale MOdel (COSMO) experiments. Simultaneous and comparative meteorological measurements were made over a pair of scale models with different block geometries. These data were analyzed to investigate the influence of height variations, obstacle elongation, and packing density, λ p , of blocks on the aerodynamic properties. In addition, the robustness of theoretical expressions of bulk transfer coefficients for momentum and heat with respect to geometric parameters was examined. Our analyses showed: (1) the theoretical framework for the bulk transfer coefficient for momentum, C m , and that for heat, C h , was applicable for homogeneous building arrays, (2) the sensitivity of C h to the surface geometry was smaller than that of C m , (3) the transfer coefficients were increased by variations of block heights, but not by elongation of blocks, (4) first-order approximations of C m and C h for an array of blocks with two different heights can be made by applying simple theoretical assumptions to include the effects of height variation, and (5) variations of block heights increased the momentum flux significantly, but caused little change in the sensible heat flux. This can be explained by the feedback mechanism of aerodynamic– thermal interaction; aerodynamic mixing decreased both the advective velocity and the vertical temperature gradient.  相似文献   

20.
The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height zi, which depend on the surface friction velocity u*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h* = |f|zi/u* > 0.15.The simple analytical model performed quite well over all values of h*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of zi and a strongly stable background stratification (h* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.  相似文献   

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

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