Turbulence Characteristics of the Stable Boundary Layer Over a Mid-Latitude Glacier. Part I: A Combination of Katabatic and Large-Scale Forcing

Profile and eddy-correlation (heights of 4 and 10 m) measurements performed on the Pasterze glacier (Austria) are used to study the characteristics of the stable boundary layer under conditions of katabatic and large-scale forcing. We consider cases where large-scale forcing results in a downslope (or following) ambient wind. The analysis of averaged spectra and cospectra reveals low frequency perturbations that have a large influence on the variances of temperature and horizontal wind components and also alter the cospectra of momentum and sensible heat flux. Only the spectrum of the vertical wind speed is comparable to universal spectra. The low frequency perturbations occur as brief intermittent events and result in downward entrainment of ambient air thereby producing enhanced downward sensible heat fluxes and downward as well as upward momentum fluxes with various magnitudes and timescales. After the variances were high pass filtered, the normalised standard deviations of wind speed and temperature compare favourably to findings in the literature within the range 0>z/L>0.5. For larger z/L they deviate as a result of an increased influence from low frequency perturbations and thus non-stationarity. In line with this, the turbulent kinetic energy budget (at 4 m height) indicates that production (shear) is in balance with destruction (buoyancy and dissipation) within the range 0>z/L>0.3. Non-dimensional gradients of wind speed within the range 0>z/L>0.3 have a slope of about 3.5. The scatter for the dimensionless temperature gradient is quite large, and the slope is comparable to that for wind speed gradients. For z/L>0.3 the imbalance in the turbulent kinetic energy budget grows and non-dimensional gradients for wind speed and temperature deviate considerably from accepted values as a result of increased non-stationarity. Average roughness lengths for momentum and sensible heat flux derived from wind speed and temperature profiles are respectively 1 × 10^-3 m and 6 × 10^-5 m, consistent with the literature. The ratio (z_0h/z_0m) compares to those predicted by surface renewal models. A variation of this ratio with the roughness Reynolds number is not indicated by our data.     

Seafloor sediment classification from single beam echo sounder data using LVQ network

Seafloor sediment classification based on echo characteristics obtained from single-beam echosounder is very useful in remote and instant sediment classification. Results of different classification techniques using such data provide robust results when the acoustic beam has a normal incidence with the seabottom. This may not always be true and show poor classification, with the data acquired during rough sea periods corresponding to both oblique and normal incidence of the acoustic pulse, due to roll and pitch motion of the ship. In the present study, an attempt is made to exploit the artificial neural network (ANN) techniques for better classification with such data. Learning Vector Quantisation (LVQ) is a supervised learning algorithm of ANN that is found to be an effective tool and show good performance. The input data to the network include the roughness index (E1) and hardness index (E2) derived from echo characteristics. The network utilizes the competitive learning, a distance function in the first layer and a linear function in the second layer. The network was tried with a different size of hidden neurons and training data size to see the influence on classification. It is found that with ten neurons in the first layer and four neurons in the second layer good performance in classification for the data was achieved.     

Palaeovelocity: A parsimonious proposal

Using data collected in her field studies, Maizels (1983) compared five methods of calculating palaeovelocity from coarse gravel deposits. Her paper indicated that the five methods provide estimates which differ by up to 2 x. There are several errors in the published equations, however. Upon correction, the smallest estimate is within 3% per cent of the largest. Other published data sets suggest that the distinction among a variety of palaeovelocity equations is not worth pursuing in light of the amount of scatter in typical data. This paper recommends use of the Thompson and Campbell (1979) modification of the basic Keulegan flow resistance equation because it fits the data as well as any other, whilst covering a more comprehensive range of flow conditions. Problems of appropriate roughness length and entrainmnet function remain.     

Momentum Absorption in Rough-Wall Boundary Layers with Sparse Roughness Elements in Random and Clustered Distributions

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 ₀/h = f(λ) between the normalised roughness length z ₀/ h and the roughness density λ (where z ₀ 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 ₀/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 ₀/h = f(λ) relationship, through increased local horizontal heterogeneity; (3) this scatter obscures any trend in the z ₀/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.     

Sensitivity of high-resolution tropical cyclone intensity forecasts to surface flux parameterization

Surface flux parameterization schemes used in current dynamic models are primarily based upon measurements at low and moderate wind speeds. Recent studies show that these parameterization schemes may be incorrect at high wind speeds (e.g., tropical cyclone forecasts). Five high-resolution numerical model experiments are designed to assess the sensitivity of tropical cyclone intensity forecasts to changes in the surface flux parameterization. The sensitivity experiments are conducted by running 48 h forecasts of the Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) for six selected tropical cyclones with individual modifications to surface flux calculation that include: (1) limiting the surface stress for wind speeds greater than 33 m s⁻¹, or 64 knots (kt); (2) computing the stress at the top of the model bottom grid layer (MBGL) by averaging results from surface layer similarity and turbulence mixing parameterization for wind speeds greater than 33 m s⁻¹; (3) increasing the roughness lengths for heat and moisture transfer by a factor of ten; (4) setting the roughness lengths for heat and moisture transfer to 1/10 of the momentum roughness length; and (5) cooling the sea surface temperature (SST) by a prescribed rate at high winds. Averaged responses for the six storms to these sensitivity tests show that: (i) the limit on surface stress at high winds significantly increases the cyclone intensity in 48 h forecasts; (ii) the averaged surface layer stress at high winds increases the cyclone intensity but to a much lesser degree than limiting the surface stress; (iii) large increases in the roughness lengths for heat and moisture transfer are needed to significantly impact the intensity forecast; (iv) the different roughness length formula for surface transfer coefficients notably increases C _h/C _d ratio from 0.59 to 0.79 for 25 m s⁻¹ and 0.41 to 0.75 for 50 m s⁻¹ that significantly increases the predicted cyclone intensity; and (v) cooling of the SST by −5.8°C in 48 h reduces the maximum surface wind speed by −32 kt, or 16.5 m s⁻¹, at 48 h forecast. These results suggest that a surface flux parameterization scheme suitable for tropical cyclone intensity forecast must correctly model the leveling-off character of surface stress and C _h/C _d ratio at high winds. All modifications to surface flux calculation have little influence on 48 h track forecasts, even though they may significantly impact the intensity forecasts.

A Roughness Sublayer Wind Profile Above A Non-Uniform Surface

In atmospheric models for different scales the underlying surface consists of patches of bare soil and plant communities with different morphological parameters. Experimental evidence indicates that there is a significant departure of the wind profile above a vegetative surface from that predicted by the logarithmic relationship, which gives values that are greater than those observed. This situation can seriously disturb the physical picture concerning the transfers of momentum, heat and water vapour from the surface into the atmosphere.The intention of this paper is to generalise the calculation of exchange of momentum between the atmosphere and a non-homogenous vegetative surface, and to derive a general equation for the wind speed profile in a roughness sublayer under neutral conditions. Furthermore, these results are extended to non-neutral cases. The suggested expression for the wind profile is compared with some earlier approaches and the observations obtained above a broad range of plant communities.     

Statistics and spectra of turbulence under different roughness in the near-neutral atmospheric surface layer

Accurate knowledge of the contacts between surface roughness and the resultant wind speed are important for climatic models, wind power meteorology, agriculture and erosion hazards especially on sand saltation in arid and semi-arid environments, where vegetation cover is scarce. In this study, synchronous measurements of three-dimensional wind speed below 5 m are carried out in three different surface roughness conditions in Minqin, China, and the difference in the turbulence statistics and the structure of the very large-scale motions (VLSMs) were revealed. The results show that the slope of the mean wind profile (MVP), the turbulent kinetic energy (TKE) and Reynolds stress increase with the surface roughness. The roughness seems to suppress the ejection events and the surface roughness will not only weaken the energy of the VLSMs, but also reduce the scale values of VLSMs near the wall. These influences may cause some changes regarding the dust transportation in streamwise and vertical directions during the sand and dust storm (SDS). That is, the decrease of the mean velocity near the ground will reduce the dust transportation in the streamwise direction and influence of the roughness on the ejection and sweep events will change the dust transportation in the vertical direction. Furthermore, the increase of roughness will weaken the scale and energy of VLSMs, which will lead to the decrease of the capacity of dust transportation. © 2019 John Wiley & Sons, Ltd.     

近岸陆地与海面风速关系分析研究

分析了近岸陆地与海上风速同步观测资料(1989年5—12月),发现海上风速与陆上风速之比并非总大于1。海、陆风速之比值不仅有明显的日变化,也随季节不同而异。进一步研究指出,海、陆风速比值是陆地地面粗糙度、大气稳定度、测站离岸线距离以及海面风速等因素的函数。利用大气边界层的理论和已有实验结果解释了海、陆风速比值的观测事实,并给出了该比值变化规律的一些定性结论。     

Introduction of a new friction routine into the SWAN model that evaluates roughness due to bedform and sediment size changes

The significant loss of wave energy due to seabed interaction in finite depths is a known effect and bottom friction terms are used in the wave models to account for this dissipation. In this paper, a new bottom-interaction function is tested by means of the SWAN model, based on measurements at two field sites, Lake George and Lakes Entrance, both in Australia. The function accounts for dependence of the friction on the formation process of bottom ripples and on the grain size of the sediment. The overall improvement of the model prediction both for the wave height and wave period is demonstrated.     

Terrestrial Laser Scanning of grain roughness in a gravel-bed river

This paper demonstrates the application of Terrestrial Laser Scanning (TLS) to determine the full population of grain roughness in gravel-bed rivers. The technique has the potential to completely replace the need for complex, time-consuming manual sampling methods. Using TLS, a total of 3.8 million data points (mean spacing 0.01 m) were retrieved from a gravel bar surface at Lambley on the River South Tyne, UK. Grain roughness was extracted through determination of twice the local standard deviation (2σ_z) of all the elevations in a 0.15 m radius moving window over the data cloud. 2σ_z values were then designated to each node on a 5 cm regular grid, allowing fine resolution DEMs to be produced, where the elevation is equivalent to the grain roughness height. Comparisons are made between TLS-derived grain roughness and grid-by-number sampling for eight 2 m² patches on the bar surface. Strong relationships exist between percentiles from the population of 2σ_z heights with measured a-, b-, and c-axes, with the closest matches appearing for the c-axis. Although strong relationships exist between TLS-derived grain roughness (2σ_z), variations in the degree of burial, packing and imbrication, results in very different slope and intercept exponents. This highlights that conventional roughness measurement using gravel axis length should be used with caution as measured axes do not necessarily represent the actual extent to which the grain protrudes into the flow. The sampling error inherent in conventional sampling is also highlighted through undertaking Monte Carlo simulation on a population of 2000 clasts measured using the grid-by-number method and comparing this with the TLS-derived population of grain roughness heights. Underestimates of up to − 23% and overestimates of up to + 50% were found to occur when considering the D₈₄, and − 20% and overestimates of up to + 36% were found to occur when considering the D_50.