Parameterizing turbulent exchange over sea ice: the ice station weddell results

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, z₀, 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 z₀. Roughness lengths for temperature, z_T, computed from this data set corroborate the theoretical model that Andreas published in 1987. Roughness lengths for humidity,z_Q, do not support this model as conclusively but are all, on average, within an order of magnitude of its predictions. Only rarely arez_Tand z_Q equal to z₀. These parameterizations have implications for models that treat the atmosphere-ice-ocean system.     

Critical Evaluation of Scalar Roughness Length Parametrizations Over a Melting Valley Glacier

We present a field investigation over a melting valley glacier on the Tibetan Plateau. In the ablation zone, aerodynamic roughness lengths (z _0M ) vary on the order of 10⁻⁴–10⁻² m, whose evolution corresponds to three melt phases with distinct surface cover and moisture exchange: snow (sublimation/evaporation), bare ice (deposition/condensation), and ice hummocks (sublimation/evaporation). Bowen-ratio similarity is validated in the stably stratified katabatic winds, which suggests a useful means for data quality check. A roughness sublayer is regarded as irrelevant to the present ablation season, because selected characteristics of scalar turbulence over smooth snow are quite similar to those over hummocky ice. We evaluate three parametrizations of the scalar roughness lengths (z _0T for temperature and z _0q for humidity), viz. key factors for the accurate estimation of sensible heat and latent heat fluxes using the bulk aerodynamic method. The first approach is based on surface-renewal models and has been widely applied in glaciated areas; the second has never received application over an ice/snow surface, despite its validity in (semi-)arid regions; the third, a derivative of the first, is proposed for use specifically over rough ice defined as z _0M > 10⁻³ m or so. This empirical z _0M threshold value is deemed of general relevance to glaciated areas (e.g. ice sheet/cap and valley/outlet glaciers), above which the first approach gives notably underestimated z _0T,q . The first and the third approaches tend to underestimate and overestimate turbulent heat/moisture exchange, respectively, frequently leading to relative errors higher than 30%. Comparatively, the second approach produces fairly low errors in energy flux estimates both in individual melt phases and over the whole ablation season; it thus emerges as a practically useful choice to parametrize z _0T,q in glaciated areas. Moreover, we find all three candidate parametrizations unable to predict diurnal variations in the excess resistances to humidity transfer, thus encouraging more efforts for improvement.     

The Parameterisation of Scalar Transfer over Rough Ice

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 ₀). 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 ₀ 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 ₀ increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z ₀ = 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 ₀  >  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 ₀)  =  1.5  ? 0.2 ln (Re _*)  ? 0.11(ln (Re _*))².     

Aerodynamic Roughness Length of Fresh Snow

This study presents the results from a series of wind-tunnel experiments designed to investigate the aerodynamic roughness length z ₀ 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 ₀. 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.05 mm. In this study, we show that variations in z ₀ are associated with variations in the roughness geometry. The roughness measurements suggest that the estimated values of z ₀ are consistent with the presence of irregular roughness structures that develop during snowfalls that mimic ballistic deposition processes.     

A simple extension of the Louis method for rough surface layer modelling

The surface flux calculation method proposed by Louis (1979) is extended by allowing momentum and heat to have different roughness lengths (z ₀ andz T respectively). Our approach is to extend the traditional analysis by a more careful consideration of potential temperature structure near the surface. This consideration leads to a redefinition of the bulk Richardson number. For bulk transfer coefficients, our method shows a significant improvement over the original Louis method when compared with the theoretical surface-layer model applied to rough surfaces. Numerical experiments simulating sea breezes in 2D show that our extension is crucially important in simulating light wind and low humidity conditions.     

Heat and mass transfer to and from surfaces with dense vegetation or similar permeable roughness

A differential equation is obtained to describe the concentration of passive admixtures (water vapor, sensible heat, pollutants, CO₂, etc.) of turbulent flow inside a dense and uniform vegetational canopy. The profiles of eddy diffusivity, wind speed and shear stress are assumed to be exponential decay functions of depth below the top of the canopy. This equation is solved for the case of a vegetation with constant concentration of the admixture at the foliage surfaces. The solution is used to formulate bulk mass or heat transfer coefficients, which can be applied to practical problems involving surfaces covered with a vegetation or with similar porous or fibrous roughness elements. The results are shown to be consistent with experimental data presented by Chamberlain (1966), Garratt and Hicks (1973) and Garratt (1978). Calculations with the model illustrate that, as compared to its behavior over surfaces with bluff roughness elements, ln(z ₀/z _0c ) (wherez ₀ is the momentum roughness andz _0c , the scalar roughness) for permeable roughness elements is relatively insensitive tou _* and practically independent ofz ₀.     

A theory for the scalar roughness and the scalar transfer coefficients over snow and sea ice

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.     

The estimation of the surface-layer parameters from wind velocity,temperature and humidity profiles by least-squares methods

The estimation of the surface-layer parameters u _* (friction velocity), _* and q _* (temperature and humidity scales), _r and q _r (temperature and humidity reference values), z _o (roughness length) and d (zero-displacement) from vertical profiles of wind velocity, temperature and humidity by least-squares methods is described. The estimation is based on the flux-gradient relationships and the constant flux assumption for the transfer of momentum, sensible heat and matter near the Earth's surface.Test calculations were carried out with the vertical profile data from the GREIV I 1974 experiment and the Great Plains Turbulence Project.     

The Fallacy of Drifting Snow

A common parametrization over snow-covered surfaces that are undergoing saltation is that the aerodynamic roughness length for wind speed (z ₀) scales as au_*²/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 ₀ measured over snow demonstrate proportionality to u_*²{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 ₀ in the first place; the plots thus suffer from fictitious correlation that causes z ₀ to unavoidably increase with u _* without any intervening physics. For these two datasets, when I plot z ₀ against u _* derived from a bulk flux algorithm—and thus minimize the fictitious correlation—z ₀ is independent of u _* in the drifting snow region, u _* ≥ 0.30 ms⁻¹. I conclude that the relation z₀ = au_*²/g{z_0 = \alpha u_\ast^2/g} when snow is drifting is a fallacy fostered by analyses that suffer from fictitious correlation.     

Turbulence measurements above a pine forest

Eddy fluxes of momentum, sensible and latent heat, and turbulence spectra measured over the Thetford Forest during 10 days in the Spring of 1973 are described. The measured total heat flux (H + E) for 122 20-min periods agreed closely on average with independent estimates from an energy balance method. There was evidence that the energy balance data gave small systematic overestimates of available energy during the hours before noon, compensated by slight underestimates for the remainder of the day. A comparison of measured wind speeds and friction velocities in neutral stability confirmed the validity of the aerodynamic method for estimating momentum fluxes at heights of a few roughness lengths above the canopy. In stable conditions the log-linear wind profileU = (u _*/k)(ln ((z -d)/z _o) + (z -d -z _o)/L) with = 3.4 ± 0.4 provided a good fit to the data. Spectra in unstable conditions were generally more sharply peaked than those measured by other workers over smoother terrain: differences were less marked in the case of vertical velocity in stable conditions. Temperature spectra in these stable conditions showed high energy at relatively low wavenumbers, andwT cospectra showed a cospectral gap; both of these results were associated with an intermittent sawtooth structure in the temperature fluctuations.Now at the Meteorological Office, Bracknell     

Radiometrically determined skin temperature and scalar roughness to estimate surface heat flux part II: Performance of parameterized scalar roughness for the determination of sensible heat

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.     

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet

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 ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ 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 ₀ 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 ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ 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 ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ 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.     

A comparison of surface renewal theory with the observed roughness length for temperature on a melting glacier surface

The roughness lengths for momentum and temperature are calculated using the profile method on amelting glacier surface. Data from a 5-level 9-m meteorological mastpositioned near the edge of Breidamerkurjökull, an outlet glacier of the Vatnajökull ice cap Iceland, are used for the calculations. The data are selected to avoid the presence of the katabatic wind speedmaximum which would otherwise alter the scaling laws of the surface layer. The surface roughness length for momentum is determined to be 1.0 mm, similar to other estimates made on flat melting ice surfaces. The surface roughness length for temperature is found to be in good agreement with previously proposed surface renewal theories for the observed roughness Reynolds number range of 30 _* 70.     

Heat and Mass Transfer to and from Surfaces with Dense Vegetation or Similar Permeable Roughness

A differential equation is obtained to describe the concentration of passive admixtures (water vapor, sensible heat, pollutants, CO₂, etc.) of turbulent flow inside a dense and uniform vegetational canopy. The profiles of eddy diffusivity, wind speed and shear stress are assumed to be exponential decay functions of depth below the top of the canopy. This equation is solved for the case of a vegetation with constant concentration of the admixture at the foliage surfaces. The solution is used to formulate bulk mass or heat transfer coefficients, which can be applied to practical problems involving surfaces covered with a vegetation or with similar porous or fibrous roughness elements. The results are shown to be consistent with experimental data presented by Chamberlain (1966), Garratt and Hicks (1973) and Garratt (1978). Calculations with the model illustrate that, as compared to its behavior over surfaces with bluff roughness elements, ln(z ₀/ z _oc ) (where z ₀ is the momentum roughness and Z _oc the scalar roughness) for permeable roughness elements is relatively insensitive to u _* and practically independent of z ₀.     

Estimation of the sensible heat flux of a semi-arid area using surface radiative temperature measurements

In the framework of an international field program for the study of semi-arid areas, observations were done in the region called La Crau in southern France. In this paper, the use of the surface radiative temperature for the determination of the sensible heat flux is addressed. We found that, once proper values of the roughness length of momentum (z ₀) and heat (z ₀h) are set, the sensible heat flux can be reliably predicted with a one-layer resistance model using standard observations of wind speed and air temperature, together with the surface temperature. The latter quantity has to be known with a precision better than ±2°C. From our observations, the value of the parameterB ^–1k ^–1 In (z ₀ z ₀h) was found to be 9.2, which falls between values quoted by Brutsaert (1982) for grass and bluff bodies.     

Momentum and heat transfers in the surface layer over a frozen sea

Wind and temperature profiles from a 10-m mast were measured over the frozen Baltic Sea. The Monin-Obukhov similarity theory applies well to runs carefully selected according to stationarity criteria. This provides relatively low-scattered results with the roughness length z ₀ 0.04 cm, the drag coefficient C _D 1.5 × 10^-3 and the Stanton number C _H 1.00 × 10^-3 under near-neutral conditions. The roughness length, however, behaves in a peculiar way under intense stratification conditions. The reasons proposed for this could lead to an extension of the theory. The mechanisms for momentum and heat transfers are also examined, revealing that there are two different regimes, the smooth and the rough, for the wall friction but not for the heat transfer. Further, the scatter of the data for momentum and heat transfer under the aerodynamically rough regime can be explained by the distinction between the type of thermal stability conditions.     

An Improvement of the Louis Scheme for the Surface Layer in an Atmospheric Modelling System

The Louis scheme for calculating the vertical eddy fluxes within the atmospheric surface layer is improved by broadening the original assumptions. In our approach, the momentum and heat transfer roughness lengths (z₀ and z_T respectively) can be different, and z₀ need not be negligibly small compared with the lowest height (z) in modelling. For these conditions, we choose more consistent wind and potential temperature profile forms, then derive new algorithms for calculating fluxes. Improvement is demonstrated for a wide range of z/L (L is the Obukhov length), z/z₀ and z₀ z_T, by comparing these fluxes with those derived from a theoretical surface-layer model. The improved algorithms can be used in atmospheric modelling systems for more varied surfaces and a wide range of atmospheric stability.     

On the role of roughness lengths in flux parameterizations of boundary-layer models

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.     

Flux-gradient profiles for momentum and heat over an urban surface

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

Numerical simulation of turbulent convective flow over wavy terrain

By means of a large-eddy simulation, the convective boundary layer is investigated for flows over wavy terrain. The lower surface varies sinusoidally in the downstream direction while remaining constant in the other. Several cases are considered with amplitude up to 0.15H and wavelength ofH to 8H, whereH is the mean fluid-layer height. At the lower surface, the vertical heat flux is prescribed to be constant and the momentum flux is determined locally from the Monin-Obukhov relationship with a roughness lengthz _o=10^–4 H. The mean wind is varied between zero and 5w _*, wherew _* is the convective velocity scale. After rather long times, the flow structure shows horizontal scales up to 4H, with a pattern similar to that over flat surfaces at corresponding shear friction. Weak mean wind destroys regular spatial structures induced by the surface undulation at zero mean wind. The surface heating suppresses mean-flow recirculation-regions even for steep surface waves. Short surface waves cause strong drag due to hydrostatic and dynamic pressure forces in addition to frictional drag. The pressure drag increases slowly with the mean velocity, and strongly with /H. The turbulence variances increase mainly in the lower half of the mixed layer forU/w _*>2.