首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
The standard E – model generates aplanetary boundary layerthat appears to be much too deep. The cause of theproblem is traced to the equation for the dissipationrate () of turbulent kinetic energy (E), specifically theparameterization of dissipation production anddestruction. In the context of atmosphericboundary-layer modelling, we argue that a part of thedissipation production should be modelled as the inputto the spectral cascade from the energy-containingpart of the spectrum, with a characteristic length , while the equilibrium imbalancebetween local production and destruction ofdissipation is modelled as proportional toE2/E, as in the standard model. Wepropose an E – – turbulence closurescheme, in which both the mixing length, m, and are prescribed. The importance ofthe equation is diminished, though itstill determines the dissipation rate in the Eequation.  相似文献   

2.
Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both SO and AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m3. In the case of aerosol carbon, we obtained a more gradual increase of AC up to an LWC of 0.5 g/m3. At higher LWCs, _ remained relatively constant at 0.60. As the differences between SO and A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.  相似文献   

3.
An E- turbulence model is used to study air-sea interaction characteristics and turbulence structure using a coupled model for air-sea boundary layers. The E- turbulence model consists of equations for the turbulent kinetic energy, the energy-dissipation, and for the turbulent exchange coefficient expressed in terms of turbulent kinetic energy and energy-dissipation. The energy-dissipation equations for the air-sea interface are solved analytically to obtain boundary conditions for energy-dissipation at the interface. The air-sea interaction and turbulence characteristics of the E- model are compared with those of the mixing-length model and with available observations.The simulations demonstrate that the air-sea interaction parameters obtained by the E- model agree well with observations. The numerical studies also show that the E- turbulence model with appropriate constants can give good results in modeling coupled air-sea boundary-layer flows.  相似文献   

4.
Several commonly used turbulence closure schemes forthe atmospheric boundarylayer (ABL) are applied to simulate neutral, nocturnal and diurnal cyclesituations in a one-dimensional ABL. Results obtained with the differentschemes, E-, E- and its modified versions, and twoversions ofthe q 2 Level 2.5, are compared and discussed.  相似文献   

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

6.
A numerical case study with a second-order turbulence closure model is proposed to study the role of urban canopy layer (UCL) for the formation of the nocturnal urban boundary layer (UBL). The turbulent diffusion coefficient was determined from an algebraic stress model. The concept of urban building surface area density is proposed to represent the UCL. Calculated results were also compared with field observation data. The height of the elevated inversion above an urban center was simulated and found to be approximately twice the average building height. The turbulent kinetic energy k, energy dissipation rate , and turbulence intensities u 2 and w 2 increase rapidly at the upwind edge of the urban area. The Reynolds stress uw displayed a nearly uniform profile inside the UBL, and the vertical sensible heat flux w had a negative value at the inversion base height. This indicates that the downward transport of sensible heat from the inversion base may play an important role in the formation of the nocturnal UBL.  相似文献   

7.
The study is focused on the retrieval and validation of the turbulent kinetic energy dissipation rate from spectral width measurements made by a UHF wind profiler in the convective atmospheric boundary layer. The possibility to deduce , which is one of the most important parameters for boundary-layer study and monitoring, from Doppler radar spectral width measurements has a firm theoretical basis established by numerous earlier works. However, the major drawback of this approach lies in various meteorological and instrumental sources of non-turbulent spectral width broadening which have to be recognised and accounted for. In the first part of the study, the theoretical background of the retrieval is detailed and all possible causes of spectral broadening are listed and evaluated. In the second part, the method is applied to four days of UHF diurnal boundary-layer observations, collected during the TRAC-98 experiment, for which in situ aircraft measurements were available. Comparison between radar-retrieved and in situ aircraft measurements yields a fairly good agreement with a linear correlation coefficient of about 0.9 and a residual bias less than 2 × 10-4 m2 s-3. The analysis of derived from vertical and off-zenith observation leads to the recommendation, that in the boundary layer, where the wind is usually moderate, data collected by off-zenith beams should be used. Indeed, the measurement of the vertical spectral width, less affected by the large-scale broadening factor, can still be altered by the ground clutter removal.  相似文献   

8.
A mesoscale Planetary Boundary Layer (PBL) model with a simple turbulence closure scheme based on the turbulence kinetic energy (TKE) equation and the dissipation () equation is used to simulate atmospheric flow over mesoscale topography. Comparative studies with different parameterizations suggest that with a proper closure assumption for turbulence dissipation, the E-model can simulate the circulation induced by the mesoscale topography with results similar to those obtained using the E- model. On the other hand, the first-order closure using O'Brien's cubic interpolation for eddy diffusivities (K) generally produces much larger K profiles in the stable and the unstable regions, which is believed to be due to the overprediction of the height of the PBL. All models with the TKE equation yield quite similar ensemble mean fields, which are found to be little sensitive to the closure assumption for turbulence dissipation, though their predicted magnitudes of TKE and K may differ appreciably. A discussion on the diurnal evolution of the mesoscale topography-induced circulation and the spatial variations of the turbulence fluxes in the surface layer is also given based on the E- model results.  相似文献   

9.
The common representation of frequency spectra in meteorology and climatology is discussed. It is pointed out that this representation is misleading since spectral peaks and spectral gaps are obtained even when the spectrum density is monotonously decreasing in the whole frequency range. A plea is made for using the spectrum distribution function, F() which gives an unambiguous picture of the distribution of variance with frequency.  相似文献   

10.
Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km2 Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.  相似文献   

11.
A mass-flux approach is applied to observational data obtained in a convective boundary layer topped with stratocumulus clouds. The observational data were obtained from aircraft measurements during the Atlantic Stratocumulus Transition Experiment (ASTEX). A conditional sampling method is used to calculate average updraft and downdraft values. The vertical fluxes calculated with the mass-flux approach are found to be proportional to the real (measured) fluxes, with a proportionality factor being about 0.6. This value is predicted by theory for two variables having a joint Gaussian distribution function; proportionality factor = 2-1 0.637. The horizontal fractional entrainment and detrainment rates calculated from the data ( 1–2 × 10-2 m-1) are an order of magnitude higher than the rates obtained by large eddy simulations for cumulus convection ( 2–3 × 10-3 m-1) and two orders of magnitude higher than those used in modelling cumulus convection with a mass-flux scheme in an operational weather forecast model ( 3 × 10-4 m-1). A numerical mass-flux model for the thermodynamics was developed and showed that results are in good agreement when compared with measured profiles of the liquid water content.  相似文献   

12.
In this paper, we lay the foundations of a systematic mathematical formulation for the governing equations for flow through an urban canopy (e.g., coarse-scaled building array) where the effects of the unresolved obstacles on the flow are represented through a distributed mean-momentum sink. This, in turn, implies additional corresponding terms in the transport equations for the turbulence quantities. More specifically, a modified k-- model is derived for the simulation of the mean wind speed and turbulence for a neutrally stratified flow through and over a building array, where groups of buildings in the array are aggregated and treated as a porous medium. This model is based on time averaging the spatially averaged Navier--Stokes equations, in which the effects of the obstacle--atmosphere interaction are included through the introduction of a volumetric momentum sink (representing drag on the unresolved buildings in the array).The k-- turbulence closure model requires two additional prognostic equations, namely one for the time-averaged resolved-scale kinetic energy of turbulence,, and another for the dissipation rate, , of . The transport equation for is derived directly from the transport equation for the spatially averaged velocity, and explicitly includes additional sources and sinks that arise from time averaging the product of the spatially averaged velocity fluctuations and the distributed drag force fluctuations. We show how these additional source/sink terms in the transport equation for can be obtained in a self-consistent manner from a parameterization of the sink term in the spatially averaged momentum equation. Towards this objective, the time-averaged product of the spatially averaged velocity fluctuations and the distributed drag force fluctuations can be approximated systematically using a Taylor series expansion. A high-order approximation is derived to represent this source/sink term in the transport equation for . The dissipation rate () equation is simply obtained as a dimensionally consistent analogue of the equation. The relationship between the proposed mathematical formulation of the equations for turbulent flow within an urban canopy (where the latter is treated as a porous medium) and an earlier heuristic two-band spectral decomposition for parameterizing turbulence in a plant canopy is explored in detail.  相似文献   

13.
Jackson and Hunt's (1975) equation for the depth of the inner layer of flows over low hills does not depend on any closure assumption as contrarily supposed in literature. This equation contains a constant which can arbitrarily be specified. It is suggested that this inner-layer constant should be determined from experimental data. A preliminary check with some data from the Askervein experiment suggests that Jackson and Hunt's equation fits these data almost as well as Jensen's equation provided that fitted inner-layer constants are used.  相似文献   

14.
Summary A radiative transfer model has been used to determine the large scale effective 6.6 GHz and 37 GHz optical depths of the vegetation cover. Knowledge of the vegetation optical depth is important for satellite-based large scale soil moisture monitoring using microwave radiometry. The study is based on actual observed large scale surface soil moisture data and observed dual polarization 6.6 and 37 GHz Nimbus/SMMR brightness temperatures over a 3-year period. The derived optical depths have been compared with microwave polarization differences and polarization ratios in both frequencies and with Normalized Difference Vegetation Index (NDVI) values from NOAA/AVHRR. A synergistic approach to derive surface soil emissivity from satellite observed brightness temperatures by inverse modelling is described. This approach improves the relationship between satellite derived surface emissivity and large scale top soil moisture fromR 2=0.45 (no correction for vegetation) toR 2=0.72 (after correction for vegetation). This study also confirms the relationship between the microwave-based MPDI and NDVI earlier described and explained in the literature.List of Symbols f frequency [Hz] - f i(p) fractional absorption at polarizationp - h surface roughness - h h cos2 - H horizontal polarization - n i complex index of refraction - p polarization (H orV) - R s microwave surface reflectivity - T B(p) brightness temperature at polarizationp - T * normalized brightness temperature - T polarization difference (T v-T H) - T s temperature of soil surface - T c temperature of canopy - T max daily maximum air temperature - T min daily minimum air temperature - V vertical polarization - soil moisture distribution factor; also used for the constant to partition the influence of bound and free water components to the dielectric constant of the mixture - empirical complex constant related to soil texture - microwave transmissivity of vegetation (=e –) - * effective transmissivity of vegetation (assuming =0) - microwave emissivity - s emissivity of smooth soil surface - rs emissivity of rough soil surface - vs emissivity of vegetated surface - soil moisture content (% vol.) - K dielectric constant [F·m–1] - K fw dielectric constant of free water [F·m–1] - K ss dielectric constant of soil solids [F·m–1] - K m dielectric constant of mixture [F·m–1] - K o permittivity of free space [8.854·10–12 F·m–1] - high frequency limit ofK wf [F·m–1] - wavelength [m] - incidence angle [degrees from nadir] - polarization ratio (T H/T V) - b soil bulk density [gr·cm–3] - s soil particle density [gr·cm–3] - R surface reflectivity in red portion of spectrum - NIR surface reflectivity in near infrared portion of spectrum - eff effective conductivity of soil extract [mS·cm–1] - vegetation optical depth - 6.6 vegetation optical depth at 6.6 GHz - 37 vegetation optical depth at 37 GHz - * effective vegetation optical depth (assuming =0) - single scattering albedo of vegetation With 12 Figures  相似文献   

15.
Local Similarity Relationships In The Urban Boundary Layer   总被引：5，自引：3，他引：2  
To investigate turbulent structures in an urban boundary layer (UBL) with many tallbuildings, a number of non-dimensional variable groups based on turbulent observationsfrom a 325-m meteorological tower in the urban area of Beijing, China, are analyzedin the framework of local similarity. The extension of surface-layer similarity to localsimilarity in the stable and unstable boundary layer is also discussed. According to localsimilarity, dimensionless quantities of variables: e.g., velocity and temperature standarddeviations i/u*l (i=u,v,w) andT/T*l,correlation coefficients of uw and wT covariance, gradients of wind and temperaturem and h, and dissipation rates of turbulent kinetic energy (TKE) andtemperature variance and N can be represented as a functiononly of a local stability parameter z/, where is the local Obukhovlength and z is the height above ground. The average dissipation rates of TKE andtemperature variance are computed by using the u spectrum, and the uw and wTcospectra in the inertial subrange. The functions above were found to be in a goodagreement with observational behaviour of turbulence under unstable conditions, butthere were obvious differences in the stable air.  相似文献   

16.
Analytical solutions for the Ekman layer   总被引：1，自引：0，他引：1  
The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U *z (0.4,U *=friction velocity,z=elevation) to the geostrophic asymptoteK GU *d forzd is solved through an expansion in fd/U *1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU *d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U *d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U * (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U * and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U *d[1–exp(–z/d)] are developed in appendices.  相似文献   

17.
This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu *(1–e –(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u *, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.  相似文献   

18.
E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c 2 /e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c 2 = 1.92, = 1.3, e = 1). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and e to 1.1 and 1.6, respectively, giving, with c 2 = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.  相似文献   

19.
From measurements in the atmospheric surface layer over a paddy field, the Kolmogorov constants for CO2 and longitudinal wind velocity were obtained. In this study, the nondimensional dissipation rate nc = (1–16 v )-1/2 for CO2 variance and = (1–16 v )-1/4 – v for turbulent energy were used, assuming the equality of the local production term and the local dissipation term, and neglecting the divergence flux term in the budget equation. The value of the constant for CO2 was consistent with recent determinations for temperature and humidity. The constant for longitudinal wind velocity showed good agreement with other recent observations.  相似文献   

20.
Flux Footprints Within and Over Forest Canopies   总被引：8，自引：1，他引：7  
The characteristics of turbulence within a forest arespatially heterogeneous and distinct from thoseassociated with the surface boundary layer. Consequently, the size and probability distribution offlux footprints emanating from sources below aforest canopy have the potential to differ from thoseobserved above forests.A Lagrangian random walk model was used toinvestigate this problem since no analytical solutionof the diffusion equation exists. Model calculations suggest that spatialcharacteristics of flux footprint probabilitydistributions under forest canopies are muchcontracted, compared to those evaluated in the surfaceboundary layer. The key factors affecting thestatistical spread of the flux footprint, and theposition of the peak of its probability distribution,are horizontal wind velocity and the standarddeviations of vertical and horizontal velocityfluctuations. Consequently, canopies, which attenuatemean horizontal wind speed, or atmospheric conditions,which enhance vertical velocity fluctuations, willcontract flux footprint distributions mostly near thefloor of a forest. It was also found that theprobability distributions of the flux footprint arenarrower when horizontal wind velocity fluctuationsare considered, instead of the simpler case that considers only vertical velocity fluctuations and meanhorizontal wind velocity.  相似文献   

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

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