首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Based on data obtained during the Hualhe River Basin Experiment （HUBEX） in 1999, this study intends to detect the quantitative discrepancies in the momentum （τ0）, sensible heat （H0） and latent heat （E0） fluxes among six sets of similarity functions with the aerodynamic method. It also aims to clarify the applicability of the functions under stable conditions. The relative discrepancy was studied with the normalized transfer coefficients for τ0, H0 and E0, namely CD, CH and CQ, respectively. Except for one set of functions that adopted a rather small von Kármán＇s constant （0.365）, the relative discrepancy in τ0 among the other functions was less than 10%, while that in H0（E0） sometimes reached 25% when the bulk Richardson number （R/B） was less than 0.07. The absolute discrepancy in the fluxes was studied with statistical computations. Among the six sets of functions, the discrepancy in τ0, H0 and E0 sometimes reached 0.03 kg m^-1 s^-2, 4 W m^-2 and 10 W m^-2, respectively, and the discrepancy in the energy balance ratio sometimes exceeded 0.1. Furthermore, when RiB exceeded the critical value （Ric） for a specific set of functions, no fluxes could be derived with the functions. It is therefore suggested that RiB be compared with Ric before computing the fluxes if RiB is less than Ric. Finally, two sets of nonlinear similarity functions are recommended, due to their unlimited applicability in terms of RiB.  相似文献   

2.
Data collected during the SHEBA and CASES-99 field programs are employed to examine the flux–gradient relationship for wind speed and temperature in the stably stratified boundary layer. The gradient-based and flux-based similarity functions are assessed in terms of the Richardson number Ri and the stability parameter z/Λ*, z being height and Λ* the local Obukhov length. The resulting functions are expressed in an analytical form, which is essentially unaffected by self-correlation, when thermal stratification is strong. Turbulence within the stably stratified boundary layer is classified into four regimes: “nearly-neutral” (0 < z/Λ* < 0.02), “weakly-stable” (0.02 < z/Λ* < 0.6), “very-stable” (0.6 < z/Λ* < 50), and “extremely-stable” (z/Λ* > 50). The flux-based similarity functions for gradients are constant in “nearly-neutral” conditions. In the “very-stable” regime, the dimensionless gradients are exponential, and proportional to (z/Λ*)3/5. The existence of scaling laws in “extremely-stable” conditions is doubtful. The Prandtl number Pr decreases from 0.9 in nearly-neutral conditions and to about 0.7 in the very-stable regime. The necessary condition for the presence of steady-state turbulence is Ri < 0.7.  相似文献   

3.
Air temperature T a , specific humidity q, CO2 mole fraction χ c , and three-dimensional winds were measured in mountainous terrain from five tall towers within a 1 km region encompassing a wide range of canopy densities. The measurements were sorted by a bulk Richardson number Ri b . For stable conditions, we found vertical scalar differences developed over a “transition” region between 0.05 < Ri b < 0.5. For strongly stable conditions (Ri b > 1), the vertical scalar differences reached a maximum and remained fairly constant with increasing stability. The relationships q and χ c have with Ri b are explained by considering their sources and sinks. For winds, the strong momentum absorption in the upper canopy allows the canopy sublayer to be influenced by pressure gradient forces and terrain effects that lead to complex subcanopy flow patterns. At the dense-canopy sites, soil respiration coupled with wind-sheltering resulted in CO2 near the ground being 5–7 μmol mol−1 larger than aloft, even with strong above-canopy winds (near-neutral conditions). We found Ri b -binning to be a useful tool for evaluating vertical scalar mixing; however, additional information (e.g., pressure gradients, detailed vegetation/topography, etc.) is needed to fully explain the subcanopy wind patterns. Implications of our results for CO2 advection over heterogenous, complex terrain are discussed.  相似文献   

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

5.
The empirical dependence of turbulence Prandtl number (Pr) on gradient Richardson number (Ri) is presented, derived so as to avoid the effects of self-correlation from common variables. Linear power relationships between the underlying variables that constitute both Pr and Ri are derived empirically from flux and profile observations. Pr and Ri are then reconstructed from these power laws, to indicate their interdependence whilst avoiding self-correlation. Data are selected according to the stability range prior to regression, and the process is iterated from neutral to higher stability until error analysis indicates the method is no longer valid. A Butterworth function is fitted to the resulting Pr −1(Ri) regression to give an empirical summary of the analysis. The form suggests that asymptotically Pr −1 decreases as Ri 3/2. Scatter in the data increases above Ri ~ 1, however, indicating additional constraints to Pr are not captured by Ri alone in this high stability regime. The Butterworth function is analytic for all Ri > 0, and may be included in suitable boundary-layer parameterisation schemes where the turbulent diffusivity for heat is derived from the turbulent diffusivity for momentum.  相似文献   

6.
We use a reduced complexity climate model with a three-dimensional ocean component and realistic topography to investigate the effect of stratification-dependent mixing on the sensitivity of the North Atlantic subpolar gyre (SPG), and the Atlantic meridional overturning circulation (AMOC), to idealized CO2 increase and peaking scenarios. The vertical diffusivity of the ocean interior is parameterized as κ ∼ N −α, where N is the local buoyancy frequency. For all parameter values 0 ≤ α ≤ 3, we find the SPG, and subsequently the AMOC, to weaken in response to increasing CO2 concentrations. The weakening is significantly stronger for α ≥ αcr ≈ 1.5. Depending on the value of α, two separate model states develop. These states remain different after the CO2 concentration is stabilized, and in some cases even after the CO2 concentration has been decreased again to the pre-industrial level. This behaviour is explained by a positive feedback between stratification and mixing anomalies in the Nordic Seas, causing a persistent weakening of the SPG.  相似文献   

7.
In unstable conditions, the set of equations defining the aerodynamic resistance to sensible heat transfer, r a , cannot be solved analytically. An iterative technique must be used to obtain r a exactly, but this is cumbersome and time consuming. In this paper, a new, empirical equation is presented relating the ratio, Q, of the aerodynamic resistances in neutral and unstable conditions, to the bulk Richardson number, Ri B . The equation takes the form Q = a + b(–Ri) c , where a, b and c are empirical functions of (z – d)/z om . This model is shown to predict r awith a mean absolute error of 0.06 s m–1 over the ranges -15 < Ri B < 0 and 10 < (z – d)/z om < 2300. Statistical comparison with other equations that have been proposed for r a in unstable conditions indicates the superior precision of the model presented here.  相似文献   

8.
The parameterization of the dimensionless entrainment rate (w e /w *) versus the convective Richardson number (Ri δθ ) is discussed in the framework of a first-order jump model (FOM). A theoretical estimation for the proportionality coefficient in this parameterization, namely, the total entrainment flux ratio, is derived. This states that the total entrainment flux ratio in FOM can be estimated as the ratio of the entrainment zone thickness to the mixed-layer depth, a relationship that is supported by earlier tank experiments, and suggesting that the total entrainment flux ratio should be treated as a variable. Analyses show that the variability of the total entrainment flux ratio is actually the effect of stratification in the free atmosphere on the entrainment process, which should be taken into account in the parameterization. Further examination of data from tank experiments and large-eddy simulations demonstrate that the different power laws for w e /w * versus Ri δθ can be interpreted as the variability of the total entrainment flux ratio. These results indicate that the dimensionless entrainment rate depends not only on the convective Richardson number but also upon the total entrainment flux ratio.  相似文献   

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

10.
A non-iterative analytical scheme is developed for unstable stratification that parametrizes the Monin–Obukhov stability parameter \(\zeta \) (\({=}z{/}L\), where z is the height above the ground and L is the Obukhov length) in terms of bulk Richardson number (\(Ri_B\)) within the framework of Businger–Dyer type similarity functions. The proposed scheme is valid for a wide range of roughness lengths of heat and momentum. The absolute relative error in the transfer coefficients of heat and momentum is found to be less than 1.5% as compared to those obtained from an iterative scheme for Businger–Dyer type similarity functions. An attempt has been made to extend this scheme to incorporate the similarity functions having a theoretically consistent free convection limit. Further, the performance of the scheme is evaluated using observational data from two different sites. The proposed scheme is simple, non-iterative and relatively more accurate compared to the schemes reported in the literature and can be used as a potential alternative to iterative schemes used in numerical models of the atmosphere.  相似文献   

11.
Ambient respirable particles (PM10; aerodynamic diameter ≤10 μm) collected in a tropical urban environment (Delhi, India) during December 2008-November 2009 were characterized with respect to 16 US EPA priority polycyclic aromatic hydrocarbons (PAHs) and 8 major and trace metals (Fe, Mn, Cd, Cu, Ni, Pb, Zn and Cr). Concentrations of Σ16PAHs (annual mean: 74.7 ± 50.7 ng m−3, range 22.1–258.4 ng m−3) and most metallic species were at least an order of magnitude greater than values reported from similar locations worldwide. Seasonal variations in Σ16PAHs were significant (p < 0.001) with highest levels in winter while crustal and anthropogenic metals showed significant but mutually opposite seasonal dependence. Statistically significant associations were observed between chemical species and various meteorological parameters. The PAH profile was dominated by combustion-derived large-ring species (~85%) that were essentially local in origin. Principal component analysis–multiple linear regression (PCA-MLR) apportioned four sources: crustal dust (73%), vehicular emission (21%), coal combustion (4%) and industrial emission (2%) that was further validated by hierarchical cluster analysis (HCA). Temporal trend analysis showed that crustal sources were predominant in summer (p < 0.05) while the remaining sources were most active in winter. Summertime intrusions of Saharan dust were identified with the help of aerosol maps and air parcel backward trajectories. Inhalation cancer risk assessment showed that up to 3,907 excess cancer cases (357 for PAHs, 122 for Cd, 2040 for Cr (VI) and 1387 for Ni) are likely in Delhi considering lifetime inhalation exposure to these chemicals at their current concentrations.  相似文献   

12.
In stably stratified flows vertical movement of eddies is limited by the fact that kinetic energy is converted into potential energy, leading to a buoyancy displacement scale z B . Our new mixing-length concept for turbulent transport in the stable boundary layer follows a rigid-wall analogy, in the sense that we assume that the buoyancy length scale is similar to neutral length scaling. This implies that the buoyancy length scale is: ℓ B  = κ B z B , with κ B ≈ κ, the von Karman constant. With this concept it is shown that the physical relevance of the local scaling parameter z/Λ naturally appears, and that the α coefficient of the log-linear similarity functions is equal to c/κ 2, where c is a constant close to unity. The predicted value α ≈ 1/κ 2 = 6.25 lies within the range found in observational studies. Finally, it is shown that the traditionally used inverse linear interpolation between the mixing length in the neutral and buoyancy limits is inconsistent with the classical log-linear stability functions. As an alternative, a log-linear consistent interpolation method is proposed.  相似文献   

13.
An approximate expression of the sky radiance in almucantar is confirmed in this paper. Under the conditions with molecular optical depth<0.3, aerosol (or cloud) optical depth<0.6, solar zenith angle< 75° and scattering angle<30°, the errors of the sky radiance and the aerosol phase function computed from the expression are respectively less than 8 % and 10 %, but the phase functions computed from Box-Deepak formula and single-scattering formula can be 84 % and 260 % more than the exact value, respectively. Furthermore, by use of the expression we have inferred the volume scattering functions of 1°≤θ≤30° for retrieving aerosol size distribution with the conclusion that the distribution information in the range of 0.1相似文献   

14.
We analyze climate change in a cost–benefit framework, using the emission and concentration profiles of Wigley et al. (Nature 379(6562):240–243, 1996). They present five scenarios that cover the period 1990–2300 and are designed to reach stabilized concentration levels of 350, 450, 550, 650 and 750 ppmv, respectively. We assume that the damage cost in each year t is proportional to the corresponding gross world product and the square of the atmospheric temperature increase (ΔT(t)). The latter is estimated with a simple two-box model (representing the atmosphere and deep ocean). Coupling the damage cost with the abatement cost, we interpolate between the five scenarios to find the one that is optimal in the sense of minimizing the sum of discounted annual (abatement plus damage) costs over a time horizon of N years. Our method is simpler than ‘traditional’ models with the same purpose, and thus allows for a more transparent sensitivity study with respect to the uncertainties of all parameters involved. We report our central result in terms of the stabilized emission level E o and concentration level p o (i.e. their values at t = 300 years) of the optimal scenario. For the central parameter values (that is, N = 150 years, a discount rate r dis = 2%/year and a growth rate r gro = 1%/year of gross world product) we find E o  = 8.0 GtCO2/year and p o = 496 ppmv. Varying the parameters over a wide range, we find that the optimal emission level remains within a remarkably narrow range, from about 6.0 to 12 GtCO2/year for all plausible parameter values. To assess the significance of the uncertainties we focus on the social cost penalty, defined as the extra cost incurred by society relative to the optimum if one makes the wrong choice of the emission level as a result of erroneous damage and abatement cost estimates. In relative terms the cost penalty turns out to be remarkably insensitive to errors. For example, if the true damage costs are three times larger or smaller than the estimate, the total social cost of global climate change increases by less than 20% above its minimum at the true optimal emission level. Because of the enormous magnitude of the total costs involved with climate change (mitigation), however, even a small relative error implies large additional expenses in absolute terms. To evaluate the benefit of reducing cost uncertainties, we plot the cost penalty as function of the uncertainty in relative damage and abatement costs, expressed as geometric standard deviation and standard deviation respectively. If continued externality analysis reduces the geometric standard deviation of relative damage cost estimates from 5 to 4, the benefit is 0.05% of the present value G tot of total gross word product over 150 years (about $3.9 × 1015), and if further research reduces the standard deviation of relative abatement costs from 1 to 0.5, the benefit is 0.03% of G tot .  相似文献   

15.
For 390 ten-minute samples of turbulent flux, made with a trivane above a lake, the vertical alignment is determined within 0.1 ° through azimuth-dependent averaging. One degree of instrumental misalignment is found to produce an average tilt error of 9 ± 4% for momentum flux, and 4 ± 2% for heat flux. The tilt error in the vertical momentum flux depends mainly ons u/u*, and cannot be much diminished with impunity by high-pass pre-filtering of the turbulence signals. The effects of rain on trivane measurements of vertical velocity are shown to be negligible at high wind speeds, and adaptable to correction in any case.The normalized vertical velocity variance,s w/u*, appears to be proportional to the square root ofz/L for unstable stratification. For a wind speed range of 2 to 15 m s–1, the eddy correlation stresses measured at 4- and 8-m heights can be reasonably well estimated by using a constant drag coefficientC d=1.3 X 10-3, while cup anemometer profile measurements give an overestimate of eddy stress at high wind speeds. A good stress estimate is also obtained from the elevation variance; it is suggested that trivane measurement of this variance might be made from a mobile platform, e.g., a moderately stabilized spar buoy.  相似文献   

16.
Studies have shown a systematic error in the Kipp & Zonen large-aperture scintillometer (K&ZLAS) measurements of the sensible heat flux, H. We improved on these studies and compared four K&ZLASs with a Wageningen large-aperture scintillometer at the Chilbolton Observatory. The scintillometers were installed such that their footprints were the same and independent flux measurements were made along the measurement path. This allowed us to compare H and the direct scintillometer output, the refractive index structure parameter, Cn2{C_{n}^{2}} . Furthermore, spectral analysis was performed on the raw scintillometer signal to investigate the characteristics of the error. Firstly, correlation coefficients ≥ 0.99 confirm the robustness of the scintillometer method, and secondly we discovered two systematic errors: the low-Cn2{C_{n}^{2}} error and the high-Cn2{C_{n}^{2}} error. The low-Cn2{C_{n}^{2}} error is a non-linear error that is caused by high-frequency noise, and we suspect the error to be caused by the calibration circuit in the receiver. It varies between each K&ZLAS, is significant for H ≤ 50 W m−2, and we propose a solution to remove this error using the demodulated signal. The high-Cn2{C_{n}^{2}} error identified by us is the systematic error found in previous studies. We suspect this error to be caused by poor focal alignment of the receiver detector and the transmitter light-emitting diode that causes ineffective use of the Fresnel lens in the current Kipp & Zonen design. It varies between each K&ZLAS (35% up to 240%) and can only be removed by comparing with a reference scintillometer in the field.  相似文献   

17.
Summary Sublayer-Stanton numbers, Bi, of heat and matter for the interfacial sublayer over aerodynamically smooth surfaces determined for forced convective conditions by elementary and numerical integration are reviewed and evaluated. The results are based on Roth’s modified Heisenberg model for the spectral energy transfer in the equilibrium range under locally isotropic conditions and the approaches of Reichardt, Elser, Deissler, van Driest, Rannie, Sheppard, and Spalding for the normalized eddy diffusivity Km/v. The results substantiate that with the exception of Sheppard’s Km/v-approach all formulations are appropriate to provide sublayer-Stanton numbers with a sufficient degree of accuracy. From a theoretical point of view the Km/v-relationships of Roth, Reichardt, van Driest, and Spalding are to be preferred, when a turbulent Prandtl number Prt = 1 is presumed. Since within the framework of mesoscale meteorological modelling numerical integration techniques would consume too much CPU-time because of the large number of near-wall grid points, a parameterization formula for the sublayer-Stanton number is proposed and evaluated. Compared to the Bi −1-results obtained by numerical integration, this kind of parameterization leads to a relative error of less than 5 percent for roughness Reynolds numbers, ηr, ranging from 30 to 600. Received January 2001 Revised November 5, 2001  相似文献   

18.
Observations obtained over a glacier surface in a predominantlykatabatic flow and with a distinctwind maximum below 13-m height are presented. The data werecollected using a 13-m high profilemast and two sonic anemometers (at about 2.5-m and 10-m heights).The spectra at frequencies belowthat of the turbulence range appear to deviate considerably fromthe curves obtained by Kaimal andco-workers during the 1968 Kansas experiment. The characteristicsof these deviations are compared tothe observations of others in surface-layers disturbed by anykind of large-scale outer-layer (orinactive) turbulence. In our case the disturbances arelikely to be induced by the highmountain ridges that surround the glacier. Moreover, the deviationsobserved in the cospectra seemto result from an, as yet, unspecified interaction between theinactive outer-layer turbulenceand the local surface-layer turbulence. Near the distinctwind maximum turbulence production ceasedwhile turbulence itself did not, probably the result ofturbulence transport from other levels. Consequently, we studied thelocal similarity relations using w instead of u* as an alternative velocity scale. Wellbelow the wind maximum, and for relatively low stability(0< Rig <0.2), the flow behaves accordingto well established local-scaling similarity relationshipsin the stable boundary layer. For higherstability (Rig > 0.2), and near or above the wind maximum, the boundary-layer structure conforms tothat of z-less stratification suggesting that the eddy sizeis restricted by the local stability ofthe flow. In line with this we observed that the sensibleheat fluxes relate remarkably well to thelocal flow parameters.  相似文献   

19.
Results from large-eddy simulations and field measurements have previously shown that the velocity field is influenced by the boundary layer height, z i , during close to neutral, slightly unstable, atmospheric stratification. During such conditions the non-dimensional wind profile, φ m , has been found to be a function of both z/L and z i /L. At constant z/L, φ m decreases with decreasing boundary layer height. Since φ m is directly related to the parameterizations of the air–sea surface fluxes, these results will have an influence when calculating the surface fluxes in weather and climate models. The global impact of this was estimated using re-analysis data from 1979 to 2001 and bulk parameterizations. The results show that the sum of the global latent and sensible mean heat fluxes increase by 0.77 W m−2 or about 1% and the mean surface stress increase by 1.4 mN m−2 or 1.8% when including the effects of the boundary layer height in the parameterizations. However, some regions show a larger response. The greatest impact is found over the tropical oceans between 30°S and 30°N. In this region the boundary layer height influences the non-dimensional wind profile during extended periods of time. In the mid Indian Ocean this results in an increase of the mean annual heat fluxes by 2.0 W m−2 and an increase of the mean annual surface stress by 2.6 mN m−2.  相似文献   

20.
Aerosol and rain samples were collected between 48°N and 55°S during the KH-08-2 and MR08-06 cruises conducted over the North and South Pacific Ocean in 2008 and 2009, to estimate dry and wet deposition fluxes of atmospheric inorganic nitrogen (N). Inorganic N in aerosols was composed of ~68% NH4+ and ~32% NO3– (median values for all data), with ~81% and ~45% of each species being present on fine mode aerosol, respectively. Concentrations of NH4+ and NO3− in rainwater ranged from 1.7–55 μmol L−1 and 0.16–18 μmol L−1, respectively, accounting for ~87% by NH4+ and ~13% by NO3− of total inorganic N (median values for all data). A significant correlation (r = 0.74, p < 0.05, n = 10) between NH4+ and methanesulfonic acid (MSA) was found in rainwater samples collected over the South Pacific, whereas no significant correlations were found between NH4+ and MSA in rainwater collected over the subarctic (r = 0.42, p > 0.1, n = 6) and subtropical (r = 0.33, p > 0.5, n = 6) western North Pacific, suggesting that emissions of ammonia (NH3) by marine biological activity from the ocean could become a significant source of NH4+ over the South Pacific. While NO3− was the dominant inorganic N species in dry deposition, inorganic N supplied to surface waters by wet deposition was predominantly by NH4+ (42–99% of the wet deposition fluxes for total inorganic N). We estimated mean total (dry + wet) deposition fluxes of atmospheric total inorganic N in the Pacific Ocean to be 32–64 μmol m−2 d−1, with 66–99% of this by wet deposition, indicating that wet deposition plays a more important role in the supply of atmospheric inorganic N than dry deposition.  相似文献   

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

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