首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z 0) is a key parameter. We present, for the first time, spatial and temporal variations of z 0 in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z 0 is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z 0≈ 10−4 m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z 0 quickly increases well above 10−2 m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z 0 > 5  ×  10−3 m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z 0 = 10−5 m. At the beginning of winter, we observed that single snow events have the potential to lower z 0 for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.  相似文献   

2.
Vertical profiles of wind speed, temperature and humidity were used to estimate the roughness lengths for momentum (z 0), heat (z H ) and moisture (z Q) over smooth ice and snow surfaces. The profile-measurements were performed in the vicinity of a blue ice field in Queen Maud Land, East Antarctica. The values ofz 0 over ice (3·10–6 m) seem to be the smallest ever obtained over permanent, natural surfaces. The settling of snow on the ice and the loss of momentum at saltating snow particles serve as momentum dissipating processes during snow-drift events, expressed as a strong dependence ofz 0 on u#.The scalar roughness lengths and surface temperature can be evaluated from the temperature and humidity profile measurements if the ratioz H /z Q is specified. This new method circumvents the difficult measurement of surface temperature. The scalar roughness lengths seem to be approximately equal toz0 for a large range of low roughness Reynolds numbers, despite the frequent occurrence of drifting snow. Possible reasons for this agreement with theory of non-saltating flow are discussed.  相似文献   

3.
A common parametrization over snow-covered surfaces that are undergoing saltation is that the aerodynamic roughness length for wind speed (z 0) scales as au*2/g{\alpha u_\ast^2/g}, where u * is the friction velocity, g is the acceleration of gravity, and α is an empirical constant. Data analyses seem to support this scaling: many published plots of z 0 measured over snow demonstrate proportionality to u*2{u_\ast^2 }. In fact, I show similar plots here that are based on two large eddy-covariance datasets: one collected over snow-covered Arctic sea ice; another collected over snow-covered Antarctic sea ice. But in these and in most such plots from the literature, the independent variable, u *, was used to compute z 0 in the first place; the plots thus suffer from fictitious correlation that causes z 0 to unavoidably increase with u * without any intervening physics. For these two datasets, when I plot z 0 against u * derived from a bulk flux algorithm—and thus minimize the fictitious correlation—z 0 is independent of u * in the drifting snow region, u * ≥ 0.30 ms−1. I conclude that the relation z0 = au*2/g{z_0 = \alpha u_\ast^2/g} when snow is drifting is a fallacy fostered by analyses that suffer from fictitious correlation.  相似文献   

4.
We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z 0). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z 0 below a threshold value of approximately 10?3 m, the model performs well and in accord with earlier studies. However, with growing hummock size, z 0 increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z 0 = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z 0  >  10?3 m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z s /z 0)  =  1.5  ? 0.2 ln (Re *)  ? 0.11(ln (Re *))2.  相似文献   

5.
A 4-month deployment on Ice Station Weddell (ISW) in the western Weddell Sea yielded over 2000 h of nearly continuous surface-level meteorological data, including eddy-covariance measurements of the turbulent surface fluxes of momentum, and sensible and latent heat. Those data lead to a new parameterization for the roughness length for wind speed, z0, for snow-covered sea ice that combines three regimes: an aerodynamically smooth regime, a high-wind saltation regime, and an intermediate regime between these two extremes where the macroscale or `permanent' roughness of the snow and ice determines z0. Roughness lengths for temperature, zT, computed from this data set corroborate the theoretical model that Andreas published in 1987. Roughness lengths for humidity,zQ, do not support this model as conclusively but are all, on average, within an order of magnitude of its predictions. Only rarely arezTand zQ equal to z0. These parameterizations have implications for models that treat the atmosphere-ice-ocean system.  相似文献   

6.
  相似文献   

7.
Using the relationship between the bulk Richardson numberR z and the Obukhov stability parameterz/L (L is the Obukhov length), formally obtained from the flux-profile relationships, methods to estimatez/L are discussed. Generally,z/L can not be uniquely solved analytically from flux-profile relationships, and it may be defined using routine observations only by iteration. In this paper, relationships ofz/L in terms ofR z obtained semianalytically were corrected for variable aerodynamic roughnessz 0 and for aerodynamic-to-temperature roughness ratiosz 0/z T, using the flux-profile iteration procedure. Assuming the so-called log-linear profiles to be valid for the nearneutral and moderately stable region (z/L<1), a simple relationship is obtained. For the extension to strong stability, a simple series expansion, based on utilisation of specified universal functions, is derived.For the unstable region, a simple form based on utilisation of the Businger-Dyer type universal functions, is derived. The formulae yield good estimates for surfaces having an aerodynamic roughness of 10–5 to 10–1 m, and an aerodynamic-to-temperature roughness ratio ofz 0/z T=0.5 to 7.3. When applied to the universal functions, the formulae yield transfer coefficients and fluxes which are almost identical with those from the iteration procedure.  相似文献   

8.
It is essential to quantify the background reactivity of smog-chambers, since this might be the major limitation of experiments carried out at low pollutant concentrations typical of the polluted atmosphere. Detailed investigation of three chamber experiments at zero-NO x in the European Photoreactor (EUPHORE) were carried out by means of rate-of-production analysis and two uncertainty analysis tools: local uncertainty analysis and Monte Carlo simulations with Latin hypercube sampling. The chemical mechanism employed was that for methane plus the inorganic subset of the Master Chemical Mechanism (MCMv3.1). Newly installed instruments in EUPHORE allowed the measurement of nitrous acid and formaldehyde at sub-ppb concentrations with high sensitivity. The presence of HONO and HCHO during the experiments could be explained only by processes taking place on the FEP Teflon walls. The HONO production rate can be described by the empirical equation W(HONO)EUPHORE dry = a × j NO 2× exp (− T 0/T) in the low relative humidity region (RH < 2%, a = 7.3×1021 cm−3, T 0 = 8945K), and by the equation W(HONO)EUPHORE humid = W(HONO)EUPHORE dry+ j NO 2× b × RH q in the higher relative humidity region (2% < RH < 15%, b = 5.8×108 cm−3 and q = 0.36, and RH is the relative humidity in percentages). For HCHO the expression W(HCHO)EUPHORE = c × j NO 2exp (− T′0/T) is applicable (c = 3.1×1017 cm−3 and T′0 = 5686 K). In the 0–15% relative humidity range OH production from HONO generated at the wall is about a factor of two higher than that from the photolysis of 100 ppb ozone. Effect of added NO2 was found to be consistent with the dark HONO formation rate coefficient of MCMv3.1.  相似文献   

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

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

11.
Aerodynamic Roughness Length of Fresh Snow   总被引：1，自引：1，他引：0  
This study presents the results from a series of wind-tunnel experiments designed to investigate the aerodynamic roughness length z 0 of fresh snow under no-drift conditions. A two-component hot-film anemometer was employed to obtain vertical profiles of velocity statistics in a zero pressure gradient turbulent boundary layer for flow over naturally deposited snow surfaces. The roughness of these snow surfaces was measured by means of digital photography to capture characteristic length scales that can be related to z 0. Our results show that, under aerodynamically rough conditions, the mean value of the roughness length for fresh snow is \({\langle{z}_{0}\rangle= 0.24}\) mm with a standard deviation σ(z 0) = 0.05 mm. In this study, we show that variations in z 0 are associated with variations in the roughness geometry. The roughness measurements suggest that the estimated values of z 0 are consistent with the presence of irregular roughness structures that develop during snowfalls that mimic ballistic deposition processes.  相似文献   

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

13.
Specific studies about the stable isotope composition (18O/16O and D/H) of atmospheric icy conglomerations are still scarce. The present work offers, for the first time, a very detailed analysis of oxygen and hydrogen isotopic signatures of unusually large ice conglomerations, or “megacryometeors”, that fell to the ground in Spain during January 2000. The hydrochemical analysis is based on the bulk isotopic composition and systematic selective sampling (deuterium isotopic mapping) of eleven selected specimens. δ18O and δD (V-SMOW) of all samples fall into the Meteoric Water Line matching well with typical tropospheric values. The distribution of the samples on Craig's line suggests either a variation in condensation temperature and/or different residual fractions of water vapour (Rayleigh processes). Three of the largest megacryometeors exhibited unequivocally distinctive negative values (δ18O = −17.2%0 and δD = −127 %0 V-SMOW), (δ18O = −15.6%0 and δD = −112%0 V-SMOW) and (δ18O = −14.4%0 and δD = −100%0 V-SMOW), suggesting an atmospheric origin typical of the upper troposphere. Theoretical calculations indicate that the vertical trajectory of growth was lower than 3.2 km. During the period in which the fall of megacryometeors occurred, anomalous atmospheric conditions were observed to exist: a substantial lowering of the tropopause with a deep layer of saturated air below, ozone depression and strong wind shear. Moreover, these large ice conglomerations occurred during non-thunderstorm conditions, suggesting an alternative process of ice growth was responsible for their formation.  相似文献   

14.
A numerical model of airflow in the lowest 50–100 m of the atmosphere above changes in surface roughness and temperature or heat flux has been developed based on boundary layer approximations, the Businger-Dyer hypotheses for the non-dimensional wind shear and heat flux and a mixing length hypothesis.Results have been obtained for several situations, in particular, airflow with neutral upstream conditions encountering a step change in surface temperature or heat flux with no roughness change. In these cases large increases in shear stress at the outer edge of the internal boundary layer are predicted. The case of unstable upstream flow encountering a step change to zero heat flux is also considered.Two situations that may be encountered near the shores of the Great Lakes are considered.Notation B Businger-Dyer constant (= 16.0) in form for M, H - c p Specific heat at constant pressure - g Acceleration due to gravity - H Upward vertical heat flux - H 0 , H 1 Surface heat fluxes for x < 0, x 0 - k von Kármán's constant ( = 0.4) - l Mixing length - L Monin-Obukhov length - L 0 Upstream value of L - m Ratio of roughness lengths (= z 1/z 0) - RL * Non-dimensional parameter, see Equations (20, 22 and 24) - RL 1 * Same as RL * but with z 1 scaling (= mRL *) - T Scaled temperature - T 0 (z) Upstream temperature profile - u 0, u 1(x) Surface friction velocities for x < 0, x 0 - U, W Horizontal and vertical mean velocities - U 0 (z) Upstream velocity profile - x, z Horizontal and vertical coordinates - z i Local roughness length  相似文献   

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

16.
 The effect of a snow cover on sea ice accretion and ablation is estimated based on the ‘zero-layer’ version sea ice model of Semtner, and is examined using a coupled atmosphere-sea ice model including feedbacks and ice dynamics effects. When snow is disregarded in the coupled model the averaged Antarctic sea ice becomes thicker. When only half of the snowfall predicted by the atmospheric model is allowed to land on the ice surface sea ice gets thicker in most of the Weddell and Ross Seas but thinner in East Antarctic in winter, with the average slightly thicker. When twice as much snowfall as predicted by the atmospheric model is assumed to land on the ice surface sea ice also gets much thicker due to the large increase of snow-ice formation. These results indicate the importance of the correct simulation of the snow cover over sea ice and snow-ice formation in the Antarctic. Our results also illustrate the complex feedback effects of the snow cover in global climate models. In this study we have also tested the use of a mean value of 0.16 Wm-1 K-1 instead of 0.31 for the thermal conductivity of snow in the coupled model, based on the most recent observations in the eastern Antarctic and Bellingshausen and Amundsen Seas, and have found that the sea ice distribution changes greatly, with the ice becoming much thinner by about 0.2 m in the Antarctic and about 0.4 m in the Arctic on average. This implies that the magnitude of the thermal conductivity of snow is of considerable importance for the simulation of the sea ice distribution. An appropriate value of the thermal conductivity of snow is as crucial as the depth of the snow layer and the snowfall rate in a sea ice model. The coupled climate models require accurate values of the effective thermal conductivity of snow from observations for validating the simulated sea ice distribution under the present climate conditions. Received: 20 November 1997/Accepted: 27 July 1998  相似文献   

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

18.
A coupled atmosphere-ocean model developed at the Institute for Space Studies at NASA Goddard Space Flight Center (Russell et al., 1995) was used to verify the validity of Haney-type surface thermal boundary condition, which linearly connects net downward surface heat flux Q to air / sea temperature difference △T by a relaxation coefficient k. The model was initiated from the National Centers for Environmental Prediction (NCEP) atmospheric observations for 1 December 1977, and from the National Ocean Data Center (NODC) global climatological mean December temperature and salinity fields at 1° ×1° resolution. The time step is 7.5 minutes. We integrated the model for 450 days and obtained a complete model-generated global data set of daily mean downward net surface flux Q, surface air temperature TA,and sea surface temperature To. Then, we calculated the cross-correlation coefficients (CCC) between Q and △T. The ensemble mean CCC fields show (a) no correlation between Q and △T in the equatorial regions, and (b) evident correlation (CCC≥ 0.7) between Q and △T in the middle and high latitudes.Additionally, we did the variance analysis and found that when k= 120 W m-2K-1, the two standard deviations, σQ and σk△T, are quite close in the middle and high latitudes. These results agree quite well with a previous research (Chu et al., 1998) on analyzing the NCEP re-analyzed surface data, except that a smaller value of k (80 W m-2K-1) was found in the previous study.  相似文献   

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

20.
The goal of this paper is to quantitatively formulate some necessary conditions for the development of intense atmospheric vortices. Specifically, these criteria are discussed for tropical cyclones (TC) and polar lows (PL) by using bulk formulas for fluxes of momentum, sensible heating, and latent heating between the ocean and the atmosphere. The velocity scale is used in two forms: (1) as expressed through the buoyancy flux b and the Coriolis parameter lc for rotating fluids convection, and (2) as expressed with the cube of velocity times the drag coefficient through the formula for total kinetic energy dissipation in the atmospheric boundary layer. In the quasistationary case the dissipation equals the generation of the energy. In both cases the velocity scale can be expressed through temperature and humidity differences between the ocean and the atmosphere in terms of the reduced gravity, and both forms produce quite comparable velocity scales. Using parameters b and lc, we can form scales of the area and, by adding the mass of a unit air column, a scale of the total kinetic energy as well. These scales nicely explain the much smaller size of a PL, as compared to a TC, and the total kinetic energy of a TC is of the order 1018-1019 J. It will be shown that wind of 33 m s-1 is produced when the total enthalpy fluxes between the ocean and the atmosphere are about 700 W m-2 for a TC and 1700 W m-2 for a PL, in association with the much larger role of the latent heat in the first case and the stricter geostrophic constraints and larger static stability in the second case. This replaces the mystical role of 26oC as a criterion for TC origin. The buoyancy flux, a product of the reduced gravity and the wind speed, together with the atmospheric static stability, determines the rate of the penetrating convection. It is known from the observations that the formation time for a PL reaching an altitude of 5--6 km can be only a few hours, and a day, or even half a day, for a TC reaching 15--18 km. These two facts allow us to construct curves on the plane of Ts and ΔT=Ts-Ta to determine possibilities for forming an intense vortex. Here, Ta is the atmospheric temperature at the height z=10 m. A PL should have ΔT>20oC in accordance with the observations and numerical simulations. The conditions for a TC are not so straightforward but our diagram shows that the temperature difference of a few degrees, or possibly even a fraction of a degree, might be sufficient for TC development for a range of static stabilities and development times.  相似文献   

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

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