首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
A Forest SO2 Absorption Model (ForSAM) was developed to simulate (1) SO2 plume dispersion from an emission source, (2) subsequent SO2 absorption by coniferous forests growing downwind from the source. There are three modules: (1) a buoyancy module, (2) a dispersion module, and (3) a foliar absorption module. These modules were used to calculate hourly abovecanopy SO2 concentrations and in-canopy deposition velocities, as well as daily amounts of SO2 absorbed by the forest canopy for downwind distances to 42 km. Model performance testing was done with meteorological data (including ambient SO2 concentrations) collected at various locations downwind from a coal-burning power generator at Grand Lake in central New Brunswick, Canada. Annual SO2 emissions from this facility amounted to about 30,000 tonnes. Calculated SO2 concentrations were similar to those obtained in the field. Calculated SO2 deposition velocities generally agreed with published values.Notation c air parcel cooling parameter (non-dimensional) - E foliar absorption quotient (non-dimensional) - f areal fraction of foliage free from water (non-dimensional) - f w SO2 content of air parcel - h height of the surface layer (m) - H height of the convective mixing layer (m) - H stack stack height (m) - k time level - k drag coefficient of drag on the air parcel (non-dimensional) - K z eddy viscosity coefficient for SO2 (m2·s–1) - L Monin-Obukhov length scale (m) - L A single-sided leaf area index (LAI) - n degree-of-sky cloudiness (non-dimensional) - N number of parcels released with every puff (non-dimensional) - PAR photosynthetically active radiation (W m–2) - Q emission rate (kg s–2) - r b diffusive boundary-layer resistance (s m–1) - r c canopy resistance (s m–1) - r cuticle cuticular resistance (s m–1) - r m mesophyllic resistance (s m–1) - r s stomatal resistance (s m–1) - r exit smokestack exit radius (m) - R normally distributed random variable with mean of zero and variance of t (s) - u * frictional velocity scale, (m s–1) - v lateral wind vector (m s–1) - v d SO2 dry deposition velocity (m s–1) - VCD water vapour deficit (mb) - z can mean tree height (m) - Z zenith position of the sun (deg) - environmental lapse rate (°C m–1) - dry adiabatic lapse rate (0.00986°C m–1) - von Kármán's constant (0.04) - B vertical velocities initiated by buoyancy (m s–1) - canopy extinction coefficient (non-dimensional) - ()a denotes ambient conditions - ()can denotes conditions at the top of the forest canopy - ()h denotes conditions at the top of the surface layer - ()H denotes conditions at the top of the mixed layer - ()s denotes conditions at the canopy surface - ()p denotes conditions of the air parcels  相似文献   

2.
A model is developed to simulate the potential temperature and the height of the mixed layer under advection conditions. It includes analytic expressions for the effects of mixed-layer conditions upwind of the interface between two different surfaces on the development of the mixed layer downwind from the interface. Model performance is evaluated against tethersonde data obtained on two summer days during sea breeze flow in Vancouver, Canada. It is found that the mixed-layer height and temperature over the ocean has a small but noticeable effect on the development of the mixed layer observed 10 km inland from the coast. For these two clear days, the subsidence velocity at the inversion base capping the mixed layer is estimated to be about 30 mm s–1 from late morning to late afternoon. When the effects of subsidence are included in the model, the mixed-layer height is considerably underpredicted, while the prediction for the mean potential temperature in the mixed layer is considerably improved. Good predictions for both height and temperature can be obtained when values for the heat entrainment ratio,c, 0.44 and 0.68 for these two days respectively for the period from 1000 to 1300 LAT, were used. These values are estimated using an equation including the additional effects on heat entrainment due to the mechanical mixing caused by wind shear at the top of the mixed layer and surface friction. The contribution of wind shear to entrainment was equal to, or greater than, that from buoyant convection resulting from the surface heat flux. Strong wind shear occurred near the top of the mixed layer between the lower level inland flow and the return flow aloft in the sea breeze circulation.Symbols c entrainment parameter for sensible heat - c p specific heat of air at constant pressure, 1010 J kg–1 K–1 - d 1 the thickness of velocity shear at the mixed-layer top, m - Q H surface sensible heat flux, W m–2 - u m mean mixed-layer wind speed, m s–1 - u * friction velocity at the surface, m s–1 - w subsidence velocity, m s–1 - W subsidence warming,oC s–1 - w e entrainment velocity, m s–1 - w * convection velocity in the mixed layer, m s–1 - x downwind horizontal distance from the water-land interface, m - y dummy variable forx, m - Z height above the surface, m - Z i height of capping inversion, m - Z m mixed-layer depth, i.e.,Z i–Zs, m - Z s height of the surface layer, m - lapse rate of potential temperature aboveZ i, K m–1 - potential temperature step atZ i, K - u h velocity step change at the mixed-layer top - m mean mixed-layer potential temperature, K  相似文献   

3.
Summary A zonally averaged global energy balance model with feedback mechanisms was constructed to simulate (i) the poleward limits of ITCZ over the continent and over the ocean and (ii) a simple monsoon system as a result of differential heating between the continent and the ocean. Three numerical experiments were performed with lower boundary as (1) global continent, (2) global ocean and (3) continent-ocean, with freezing latitudes near the poles. Over the continent, midlatitude deserts were found and the ITCZ migrates 25° north and south with seasons. Over a global swamp ocean results do not show migration of ITCZ with time but once the ocean currents are introduced the ITCZ migrates 5° north and south with seasons. It was found that the seasonal migration of ITCZ strongly depends on the meridional distribution of the surface temperature. It was also found that continent influences the location of the oceanic ITCZ. In the tropics northward progression of quasi-periodic oscillations called events are found during the pre- and post-monsoon periods with a period of 8 to 15 days. This result is consistent with the observed quasi-periodic oscillations in the tropical region. Northward propagation of the surface temperature perturbation appears to cause changes in the sensible heat flux which in turn causes perturbations in vertical velocity and latent heat flux fields.List of Symbols vertical average - 0 zonal average - vertical mean of the zonal average - 0s zonal average at the surface - 0a zonal average at 500 mb level - latitude We now define the various symbols used in the model rate of atmospheric heating due to convective cloud formation (K/sec) - dp/dt (N/m2/sec) - density - potential temperature (K) - rate of rotation of the earth (rad/sec) - empirical constant - humidity mixing ratio - * saturated humidity mixing ratio - opacity of the atmosphere - 1,2 factors for downward and upward effective black body long wave radiation from the atmosphere - Stefan-Boltzmann constant - emissivity of the surface - D subsurface temperature (K) - a specific volume - 0xs ,0ys eastward and northward components of surface frictional stress - * vertical velocity at the top of the boundary layer (N/m2/sec) - P Thickness of the boundary layer (mb) - nondimensional function of pressure - P pressure - P a pressure of the model atmosphere (N/m2) - P s pressure at the surface (N/m2) - t time (sec) - U eastward wind speed (m/sec) - V northward wind speed (m/sec) - surface water availability - T absolute temperature (K) - heat addition due to water phase changes - g acceleration due to gravity (m2/sec) - a radius of the earth (m) - R gas constant for dry air (J/Kg/K) - C p specific heat of air at constant pressure (J/Kg/K) - k R/C p - L latent heat of condensation (J/Kg) - f coriolis parameter (rad/sec) - H s H 0s (1) +H 0s (2) +H 0s (3) +H 0s (4) +H 0s (5) (J/m2/Sec)=sum of the rates of vertical heat fluxes per unit surface area, directed toward the surface - H a H 0a (1) +H 0a (2) +H 0a (3) +H 0a (4) (J/m2/Sec)=sum of the rates of heat additions to the atmospheric column per unit horizontal area by all processes - H 0s (1) ,H 0a (1) heat flux due to short wave radiation - H 0s (2) ,H 0a (2) heat flux due to long wave radiation - H 0s (3) ,H 0a (3) heat flux due to small scale convection - H 0s (4) heat flux due to evaporation - H 0a (4) heat flux due to condensation - H 0s (5) heat flux due to subsurface conduction and convection - e * saturation vapor pressure - R solar constant (W/m2) - r a albedo of the atmosphere - r s albedo of the surface - b 2 empirical constant (J/m2/sec) - c 2 empirical constant (J/m2/sec) - e 2 nondimensional empirical constant - f 2 empirical constant (J/m2/sec) - factor proportional to the conductive capacity of the surface medium - a s constant used in Sellers model - b s positive constant of proportionality used in the Sellers model (kg m2/J/sec2) - K HT coefficient for eddy diffusivity of heat (m2/sec) - K HE exchange coefficient for water vapor (m2/sec) - h depth of the water column (m) - z height (m) - V 0ws meridional component of surface current (m/sec) - n cloud amount - G 0,n long wave radiation form the atmosphere for cloud amount n (W/m2) - B 0 long wave radiation from the surface (W/m2) - S 0,n short wave radiation from the atmosphere for cloud amount n (W/m2) - A n albedo factor for a cloud amount n - R f1 large scale rainfall (mm/day) - R f2 small scale rainfall (mm/day) With 22 Figures  相似文献   

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

5.
A set of semi-continuous measurements of temperature, wind and moisture gradients as well as of net radiation and ground heat flux covering a period of about one and a half years has been analysed to give a corresponding set of complete surface energy balance data on an hourly basis. An analysis of the evaporation data so obtained is given.It is shown that surface resistance r S exhibits a diurnal trend: values are smallest (ca. 150 s m-1) a few hours before noon and increase to as much as 800 s m-1 towards dusk. The minimum values tend to be higher during dry periods when the soil moisture is low. There is also some indication that r S decreases rapidly soon after rainfall.An exponential relation is found between the fraction of available energy used as evaporative flux, , and r S for values of r I/rS <0.70, where r I is the climatological resistance. On the other hand, the ratio of r I to r S is linearly correlated with , implying that an equilibrium state is established between the grass surface and the atmosphere, at least from mid-morning to mid-afternoon when the leaves are dry. Near-noon values calculated by Stewart and Thom for Thetford Forest also follow a linear trend.The above two regression results (In (r S) versus r I/rS versus ) are combined to obtain an empirical relation of the form r I=m exp (a-b) which is used to estimate evaporative flux. The estimates are found to be within 20% of calculated values.Guest Scientist from Department of Physics, University of Cape Coast, Cape Coast, Ghana.  相似文献   

6.
A two-dimensional mesoscale model has been developed to simulate the air flow over the Gulf Stream area where typically large gradients in surface temperature exist in the winter. Numerical simulations show that the magnitude and the maximum height of the mesoscale circulation that develops downwind of the Gulf Stream depends on both the initial geostrophic wind and the large-scale moisture. As expected, a highly convective Planetary Boundary Layer (PBL) develops over this area and it was found that the Gulf Stream plays an important role in generating the strong upward heat fluxes causing a farther seaward penetration as cold air advection takes place. Numerical results agree well with the observed surface fluxes of momentum and heat and the mesoscale variation of vertical velocities obtained using Doppler Radars for a typical cold air outbreak. Precipitation pattern predicted by the numerical model is also in agreement with the observations during the Genesis of Atlantic Lows Experiment (GALE).List of Symbols u east-west velocity [m s–1] - v north-south velocity [m s–1] - vertical velocity in coordinate [m s–1] - w vertical velocity inz coordinate [m s–1] - gq potential temperature [K] - q moisture [kg kg–1] - scaled pressure [J kg–1 K–1] - U g the east-south component of geostrophic wind [m s–1] - V g the north-south component of geostrophic wind [m s–1] - vertical coordinate following terrain - x east-west spatial coordinate [m] - y north-south spatial coordinate [m] - z vertical spatial coordinate [m] - t time coordinate [s] - g gravity [m2 s–1] - E terrain height [m] - H total height considered in the model [m] - q s saturated moisture [kg kg–1] - p pressure [mb] - p 00 reference pressure [mb] - P precipitation [kg m–2] - vertical lapse rate for potential temperature [K km–1] - L latent heat of condensation [J kg–1] - C p specific heat at constant pressure [J kg–1 K–1] - R gas constant for dry air [J kg–1 K–1] - R v gas constant for water vapor [J kg–1 K–1] - f Coriolis parameter (2 sin ) [s–1] - angular velocity of the earth [s–1] - latitude [o] - K H horizontal eddy exchange coefficient [m2 s–1] - t integration time interval [s] - x grid interval distance inx coordinate [m] - y grid interval distance iny coordinate [m] - adjustable coefficient inK H - subgrid momentum flux [m2 s–2] - subgrid potential temperature flux [m K s–1] - subgrid moisture flux [m kg kg–1 s–1] - u * friction velocity [m s–1] - * subgrid flux temperature [K] - q * subgrid flux moisture [kg kg–1] - w * subgrid convective velocity [m s–1] - z 0 surface roughness [m] - L Monin stability length [m] - s surface potential temperature [K] - k von Karman's constant (0.4) - v air kinematic viscosity coefficient [m2 s–1] - K M subgrid vertical eddy exchange coefficient for momentum [m2 s–1] - K subgrid vertical eddy exchange coefficient for heat [m2 s–1] - K q subgrid vertical eddy exchange coefficient for moisture [m2 s–1] - z i the height of PBL [m] - h s the height of surface layer [m]  相似文献   

7.
The energy balance was measured for the dry canopy of narrow-leaved snow tussock (Chionochloa rigida), and measurements of transpiration were obtained from a large weighing lysimeter.Typical maximum summer transpiration rates of 0.21–0.43 mmhr-1 (140–290 W m-2) were recorded. The latent heat flux accounted for less than 40% of net radiation. The estimated value of the bulk stomatal resistance (r ST) for 29 days was 158 s m-1, and the decoupling parameter () was 0.17. Transpiration rates were not driven directly by net radiation, but were closely linked to the size of the regional saturation deficit imposed at the level of the canopy by efficient overhead mixing, and were constrained by a large bulk stomatal resistance. A linear relationship between r ST and the saturation deficit is proposed as a realistic method for estimating transpiration for water yield studies of tussock catchments.  相似文献   

8.
The present study involved determination of the experimental energy receipt partitioning over the tropical Amazon forest. Diurnal variation of net radiation (Q *), sensible heat flux (Q H) and latent heat flux (Q E) is presented. The daytimeQ E is in phase withQ * and it is always an important term in the energy balance. The daily averagedQ E is 59 to 100% of the dailyQ * whereasQ H is 5 to 28% at the Amazon forest site (2° 57 S; 59° 57 W) for the sample periods. The results present evidence thatQ E over the Amazon forest is greater thanQ * in the afternoon hours. The role of sensible heat advection in maintaining largeQ E over the forest surface is discussed. Hourly Bowen ratio () values for two campaigns of the Amazon forest micrometeorological experiment are presented. During daylight hours, the differences in are not significant, and exhibit a systematic pattern. The only time that the variation in Bowen ratio increased significantly was at sunrise and sunset when the thermal structure of the air was changing from a strong inversion to lapse and vice versa. The diurnal values changed from –3.50 to 0.85. The mean hourly calculated from values from 07.00 to 16.00 h, varied from 0.05 to 0.85.Diese Studie beschäftigt sich mit der Aufteilung der empfangenen Energie über dem tropischen Amazonasurwald. Es wird der Tagesgang der Strahlungsbilanz (Q *), des fühlbaren (Q H) und des latenten Wärmestromes (Q E) vorgestellt. Während der Tagesstunden istQ E in Phase mitQ * und ist immer ein wichtiger Term der Energiebilanz. Das Tagesmittel vonQ E beträgt 59 bis 100%,Q H 5 bis 28% des täglichenQ * an den Meßtagen bei der Amazonasurwaldstation (2° 57 S; 59° 57 W). Die Ergebnisse legen nahe, daß in den NachmittagsstundenQ E über dem Amazonasurwald größer ist alsQ *. Die Rolle der Advektion von fühlbarer Wärme zur Aufrechterhaltung des großenQ E über der Waldoberfläche wird diskutiert. Für zwei Meßkampagnen wurden die stündlichen Bowenverhältnisse () vorgestellt. Während der Tagesstunden ergaben sich keine signifikanten Änderungen von, während bei Sonnenaufgang und -untergang, wenn der thermische Aufbau der Luft von einer starken Inversion zu neutral und umgekehrt wechselt, die Unterschiede deutlich anstiegen. Die Tageswerte von lagen zwischen –3.50 und 0.85. Die Stundenmittel von 7.00 bis 16.00 Uhr schwankten zwischen 0.05 und 0.85.
With 3 Figures  相似文献   

9.
We have devised a partial differential equation for the prediction of dust concentration in a thin layer near the ground. In this equation, erosion (detachment), transport, deposition and source are parameterised in terms of known quantities. The interaction between a wind prediction model in the boundary layer and this equation affects the evolution of the dust concentration at the top of the surface layer. Numerical integrations are carried out for various values of source strength, ambient wind and particle size. Comparison with available data shows that the results appear very reasonable and that the model should be subjected to further development and testing.Notation (x, y, z, t) space co-ordinates and time (cm,t) - u, v components of horizontal wind speed (cm s–1) - u g, vg components of the geostrophic wind (cm s–1) - V=(u2+v2)1/2 (cm s–1) - (û v)= 1/(h – k) k h(u, v)dz(cm s–1) - V * friction velocity (cm s–1) - z 0 roughness length (cm) - k 1 von Karman constant =0.4 - V d deposition velocity (cm s–1) - V g gravitational settling velocity (cm s–1) - h height of inversion (cm) - k height of surface layer (cm) - potential temperature (°K) - gr potential temperature at ground (°K) - K potential temperature at top of surface layer (°K) - P pressure (mb) - P 0 sfc pressure (mb) - C p/Cv - (t)= /z lapse rate of potential temperature (°K cm–1) - A(z) variation of wind with height in transition layer - B(z) variation of wind with height in transition layer - Cd drag coefficient - C HO transfer coefficient for sensible heat - C dust concentration (g m–3) - C K dust concentration at top of surface layer (g m–3) - D(z) variation with height of dust concentration - u, v, w turbulent fluctuations of the three velocity components (cm s–1) - A 1 constant coefficient of proportionality for heat flux =0.2 - Ri Richardson number - g gravitational acceleration =980 cm s–2 - Re Reynolds number = - D s thickness of laminar sub-layer (cm) - v molecular kinematic viscosity of air - coefficient of proportionality in source term - dummy variable - t time step (sec) - n time index in numerical equations On sabbatical leave at University of Aberdeen, Department of Engineering, September 1989–February 1990.  相似文献   

10.
We present turbulence spectra and cospectra derived from long-term eddy-covariancemeasurements (nearly 40,000 hourly data over three to four years) and the transferfunctions of closed-path infrared gas analyzers over two mixed hardwood forests inthe mid-western U.S.A. The measurement heights ranged from 1.3 to 2.1 times themean tree height, and peak vegetation area index (VAI) was 3.5 to 4.7; the topographyat both sites deviates from ideal flat terrain. The analysis follows the approach ofKaimal et al. (Quart. J. Roy. Meteorol. Soc. 98, 563–589, 1972) whose results were based upon 15 hours of measurements atthree heights in the Kansas experiment over flatter and smoother terrain. Both thespectral and cospectral constants and stability functions for normalizing and collapsingspectra and cospectra in the inertial subrange were found to be different from those ofKaimal et al. In unstable conditions, we found that an appropriate stabilityfunction for the non-dimensional dissipation of turbulent kinetic energy is of the form () = (1 - b-)-1/4 - c-, where representsthe non-dimensional stability parameter. In stable conditions, a non-linear functionGxy() = 1 + bxyc xy (cxy < 1) was found to benecessary to collapse cospectra in the inertial subrange. The empirical cospectralmodels of Kaimal et al. were modified to fit the somewhat more (neutraland unstable) or less (stable) sharply peaked scalar cospectra observed over forestsusing the appropriate cospectral constants and non-linear stability functions. Theempirical coefficients in the stability functions and in the cospectral models varywith measurement height and seasonal changes in VAI. The seasonal differencesare generally larger at the Morgan Monroe State Forest site (greater peak VAI) andcloser to the canopy.The characteristics of transfer functions of the closed-path infrared gas analysersthrough long-tubes for CO2 and water vapour fluxes were studied empirically. This was done by fitting the ratio between normalized cospectra of CO2 or watervapour fluxes and those of sensible heat to the transfer function of a first-order sensor.The characteristic time constant for CO2 is much smaller than that for water vapour. The time constant for water vapour increases greatly with aging tubes. Three methods were used to estimate the flux attenuations and corrections; from June through August, the attenuations of CO2 fluxes are about 3–4% during the daytime and 6–10% at night on average. For the daytime latent heat flux (QE), the attenuations are foundto vary from less than 10% for newer tubes to over 20% for aged tubes. Correctionsto QE led to increases in the ratio (QH + QE)/(Q* - QG) by about 0.05 to0.19 (QH is sensible heat flux, Q* is net radiation and QG is soil heat flux),and thus are expected to have an important impact on the assessment of energy balanceclosure.  相似文献   

11.
Comparisons between sensible heat flux measured using eddy correlation instrumentation and estimated using the temperature fluctuation method are presented for four types of surface in West Africa. Agreement between measured and estimated values is good. Regression of estimated on measured sensible heat flux gave a mean slope of 0.98 with a mean r 2 of 0.94 for bare soil, mature millet, fallow savannah and tiger bush. Estimates of heat flux from temperature fluctuations measured by an instrument mounted beneath a tethered balloon are also shown to be in close agreement with eddy correlation measurements made at the surface (regression slope = 0.98, r 2 = 0.84). The results provide evidence that the ratio /×is indeed a universal function of z/L for all the surface types considered.  相似文献   

12.
The dependence of radiometric surface temperature (s) on view angle and the unclear definition of the aerodynamic temperature, which is the temperature that gives the correct sensible heat flux estimate at a specified roughness length, bring about a challenge in estimating sensible heat flux from s. An analytical-land-atmosphere-radiometer model (ALARM) has been developed to convert s taken at any zenith view angle to a clearly defined equivalent isothermal surface temperature, i, at a defined scalar roughness length. ALARM is an analytical model based on K-theory that links the foliage temperature profile to the radiometric surface temperature and the temperature felt by the turbulent lower atmosphere. ALARM has previously been applied with slightly different values of its parameters to several grassland sites of varying canopy density. Our objective in this study was to apply ALARM to these and to one additional dataset with a single parameterization. When compared to the reference (measured) values of sensible heat flux H, ALARM estimates of H had root mean square errors of about 35 W m-2. These results were comparable to those from two other simple canopy models also tested with these datasets.  相似文献   

13.
The temperature drop T between the ocean surface and the 5-cm depth was recorded during GATE, Phase III. With measured values of the total heat flux Q and an assumption about the thickness of the viscous boundary layer of the ocean, the wind-speed dependence of the factor of proportionality between T and Q is determined. This factor depends on the deviations of the thickness of the conductive layer from the thickness of the viscous layer and possibly partially on the wind stress. A further assumption about the thickness of the conductive layer leads to a wind-speed dependence of the ratio between total wind stress and its wave supporting part of it. This ratio increases from a value 1.5 at 1 m s–1 to 9 at 10 m s–1, which is in agreement with existing estimates.  相似文献   

14.
Several formulations and proposals to determine the value of the radiometric scalar roughness for sensible heatz 0h,r are tested with respect to their performance in the estimation of the sensible heat flux by means of the profile equations derived from Monin-Obukhov similarity theory. The equations are applied to the data set of spatially averaged surface skin temperature and profiles of wind speed and temperature observed in a pasture field during a growing season. The use of a physical model developed for a dense canopy to estimate scalar roughness for sensible heatz 0h,r produced sensible heat fluxH with a correlation coefficientr=0.884, the ratio of means being H s /H=1.19 in a comparison with reference values ofH s . In comparison, a proposal for a fixed value ofz 0h yieldedr=0.887, H s /H=0.879. In both cases, the validity ofz 0h =z 0h,r was assumed. All expressions derived to estimatez 0h,r from a multiple linear regression with such predictors as leaf area index, solar radiation and the ratio of solar radiation to extraterrestrial radiation, were found to produce a better result, withr better than 0.90 and H s /H around 1.0. However, when the constantsc andf of a linear regression equationHs=cH+f are used to evaluate the equations, a marked difference in performance of each formulation appeared. In general, equations with smaller numbers of predictors tend to produce a biased result, i.e., an overestimation ofH at largeH s . These values ofH are used in conjunction with the energy balance equation to derive values of the latent heat fluxLE, which are shown to be in good agreement with the reference valuesLE s , withr greater than 0.97.  相似文献   

15.
The relationship between satellite-derived low-level cloud motion, surface wind and geostrophic wind vectors is examined using GATE data. In the trades, surface wind speeds can be derived from cloud motion vectors by the linear relation: V = 0.62 V s + 1.9 m s–1 with a mean scatter of ±1.3 m s–1. The correlation coefficient between surface and satellite wind speed is 0.25. Considering baroclinicity, i.e., the influence of the thermal wind, the correlation coefficient does not increase, because of the uncertainty of the thermal wind vectors. The ratios of surface to geostrophic wind speed and surface to satellite wind speed are 0.7 and 0.8, respectively, with a statistical uncertainty of ±0.3. Calculations of the ratio of surface to geostrophic wind speed on the basis of the resistance law yield V/V g = 0.8 ± 0.2, in agreement with experimental results. The mean angle difference between the surface and the satellite wind vectors amounts to - 18 °, taking into account baroclinicity. This value is in good agreement with the mean ageostrophic angle - 25 °.  相似文献   

16.
Previous measurements of urban energy balances generally have been limited to densely built, central city sites and older suburban locations with mature tree canopies that are higher than the height of the buildings. In contrast, few data are available for the extensive, open vegetated types typical of low-density residential areas that have been newly converted from rural land use. We made direct measurements of surface energy fluxes using the eddy-covariance technique at Greenwood, a recently developed exurban neighbourhood near Kansas City, Missouri, USA, during an intensive field campaign in August 2004. Energy partitioning was dominated by the latent heat flux under both cloudy and near clear-sky conditions. The mean daytime Bowen ratio (β) values were 0.46, 0.48, and 0.47 respectively for the cloudy, near clear-sky and all-sky conditions. Net radiation (R n ) increased rapidly from dawn (−34 and −58W m−2) during the night to reach a maximum (423 and 630W m−2) after midday for cloudy and near clear-sky conditions respectively. Mean daytime values were 253 and 370W m−2, respectively for the cloudy and near clear-sky conditions, while mean daily values were 114 for cloudy and 171W m−2 for near clear-sky conditions, respectively. Midday surface albedo values were 0.25 and 0.24 for the cloudy and near clear-sky conditions, respectively. The site exhibited an angular dependence on the solar elevation angle, in contrast to previous observations over urban and suburban areas, but similar to vegetated surfaces. The latent heat flux (Q E ), sensible heat flux (Q H ), and the residual heat storage ΔQ s terms accounted for between 46–58%, 21–23%, and 18–31% of R n , respectively, for all-sky conditions and time averages. The observed albedo, R n , and Q E values are higher than the values that have been reported for suburban areas with high summer evapotranspiration rates in North America. These results suggest that the rapidly growing residential areas at the exurban fringe of large metropolitan areas have a surface energy balance that is more similar to the rural areas from which they were developed than it is to the older suburbs and city centres that make up the urban fabric to which they are being joined.  相似文献   

17.
An attempt is made to construct a model, coupling land surface and atmospheric processes in the planetary boundary layer (PBL). A grassland strip in a semi-desert (hereinafter called desert) is presupposed, so as to simulate the case of heterogeneous vegetation cover.Modeling results indicate that every term in the equation of the surface energy balance changes as the air flows over the grassland. The striking contrast of water and energy conditions between the grassland and the desert means that the air over the grassland is cooler and wetter than that over the desert. Consequently, in the heating and dynamic forcing of the air by the underlying surface, heterogeneities arise and are then transferred upward by the turbulent motions. Horizontal differences thus develop in the PBL, resulting in a local circulation. Meanwhile, the horizontal differences affect the free atmosphere through vertical motion at the top of the PBL.List of symbols d 1,d 2,d 3 depths of surface, middle and lower layers of soil - T c ,T 1,T 2,T 3 temperatures of canopy, surface, middle and lower layers of soil - R nc net radiation of canopy layer - c shielding factor of vegetation - Ew, Etc evaporation from wet fraction of foliage and transpiration from dry fraction of foliage - Et 1,Et 2 transpiration of foliage water absorbed by the root in the upper and lower soil, respectively - H c sensible heat of canopy - P c ,D c precipitation rate and drainage of canopy - C s ,C c ,C w heat capacity of soil, canopy and water - w , s density of water and air near the surface - D hydraulic permeability of soil - s saturated value of the ratio of volumetric soil moisture - S g , g solar radiation and surface reflection - H g ,R L g turbulent heat flux and long wave radiation of surface - P g ,E g precipitation rate and evaporation of soil surface - K s soil thermal diffusivity - K (m),K (H),K (q) eddy coefficients of momentum, heat and moisture - u, v, w components of wind speed in three directions - air potential temperature - e turbulent kinetic energy - p atmospheric pressure - C p specific heat of air under constant pressure - R d gas constant - u * friction velocity - * feature temperature - h height of the PBL - f Coriolis parameter - L 0 Monin-Obukhov length - latent heat of vaporization - q specific humidity - M c ,M cm interception water storage of canopy and its maximum - 0 Exner number of largescale background field - perturbation Exner number - u g ,v g components of the geostrophic wind speed Sponsored by the National Natural Science Foundation of China.  相似文献   

18.
Summary The evaporation of deep crops such as forests is usually considered in terms of the two-resistance Penman-Monteith model, though this conflates two of the three resistances actually involved, i.e. the canopy resistancer c between the transpiring leaves and the top of the canopy, and the resistancer s due to the stomates of the leaves. A review of the literature on these and the aerodynamic resistancer a (between the crop and the atmosphere) shows how distinctly different they are, and therefore how inappropriate it is to lump any two together.Once the soil has dried substantially,r s depends approximately onM –2, whereM is the fractional available soil moisture.As regards grassed surfaces,r a is 300/u s m–1, whereu is the wind speed at 2 m.With 2 Figures  相似文献   

19.
A two dimensional model has been set up to investigate the circulation induced by an urban heat island in the absence of synoptic winds. The boundary conditions need to be formulated carefully and due to difficulties arising here, we restrict our attention to cases of initially stable thermal stratification. Heat island circulations are allowed to develop from rest and prior to the appearance of the final symmetric double cell pattern, a transitional multi-cell pattern is observed in some cases. The influence on the steady state circulation of various parameters is studied, among which are eddy transfer coefficients, the heat island intensity, the initial temperature stratification and the heat island size. Some results are presented for a case in which differential surface cooling beneath an initially stable atmosphere produces a circulation and an unstable layer capped by an elevated inversion over the city. It is hoped that this case is vaguely representative of the night-time heat island with no geostrophic wind.Notation cp Specific heat at constant pressure - g Acceleration due to gravity - H Top of integration region - Kz Vertical eddy transfer coefficient - Kx, KxH, Kxm Horizontal eddy transfer coefficients for heat and momentum - l ixing length - p Pressure - p0 Reference surface pressure (1000 mb) - PH (x, t) Pressure at z = H - R Specific gas constant for dry air - t Time - u, w Horizontal and vertical velocities - x, z Horizontal and vertical coordinates - x1, x2 Positions of discontinuities in surface temperature field (see Figure 2) - xa Heat island half-width - xb Boundary of integration region - Parameter in formula for eddy coefficients (variable-K case) = 18.0 - s Intensity of heat island - Potential temperature field - Reference absolute temperature (variable-K case) - r Reference temperature (° C) - s Surface temperature - Q Air density  相似文献   

20.
Summary A numerical model was used to study the behaviour of prototype cold fronts as they approach the Alps. Two fronts with different orientations relative to the Alpine range have been considered. One front approaches from west, a second one from northwest. The first front is connected with southwesterly large-scale air-flow producing pre-frontal foehn, whereas the second front is associated with westerly largescale flow leading to weak blocking north of the Alps.Model simulations with fully represented orography and parameterized water phase conversions have been compared with control runs where either the orography was cut off or the phase conversions were omitted. The results show a strong orographic influence in case of pre-frontal foehn which warms the pre-frontal air and increases the cross-frontal temperature contrast leading to an acceleration of the front along the northern Alpine rim. The latent heat effect was found to depend much on the position of precipitation relative to the surface front line. In case of pre-frontal foehn precipitation only falls behind the surface front line into the intruding cold air where it partly evaporates. In contrary, precipitation already appears ahead of the front in the case of blocking. Thus, the cooling effect of evaporating rain increases the cross-frontal temperature difference only in the first case causing an additional acceleration of the front.List of symbols C pd specific heat capacity of dry air at constant pressure (C pd =1004.71 J kg–1 K–1) - C pv specific heat capacity of water vapour at constant pressure (C pv =1845.96 J kg–1 K–1) - C f propagation speed of a front - x, y horizontal grid spacing (cartesian system) - , horizontal grid spacing (geographic system) - t time step - E turbulent kinetic energy - f Coriolis parameter - g gravity acceleration (g=9.81 ms–1) - h terrain elevation - H height of model lid (H=9000 m) - k Karman constant (k=0.4) - K Mh horizontal exchange coefficient of momentum - K Hh horizontal exchange coefficient of heat and moisture - K Mz vertical exchange coefficient of momentum - K Hz vertical exchange coefficient of heat and moisture - l mixing length - l c specific condensation heat (l c =2500.61 kJ kg–1) - l f specific freezing heat (l f =333.56 kJ kg–1) - l s specific sublimation heat (l s =2834.17 kJ kg–1) - longitude - m 1,m 2,m 3 metric coefficients - p pressure - Exner function - Pr Prandtl number - latitude - M profile function - q v specific humidity - q c specific content of cloud droplets - q i specific content of cloud ice particles - q R specific content of rain drops - q S specific content of snow - R d gas constant of dry air (R d =287.06 J kg–1 K–1) - R v gas constant of water vapour (R v =461.51 J kg–1 K–1) - r E radius of earth (r E =6371 km) - Ri F flux Richardson number - density of dry air - t time - T temperature - dia period of diastrophy - potential temperature - v virtual potential temperature - e equivalent potential temperature - U relative humidity - u, v, w cartesian wind components - u F ,v F front-normal and front-parallel wind components - x, y, z cartesian coordinates - w * transformed vertical wind component - W R speed of falling rain - W S speed of falling snow - z * transformed vertical coordinate Abbreviations GND (above) ground level - MSL (above) mean sea level With 12 Figures  相似文献   

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

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