Wave-Dependence of Friction Velocity, Roughness Length, and Drag Coefficient over Coastal and Open Water Surfaces by Using Three Databases

The parameterization of friction velocity, roughness length, and the drag coefficient over coastal zones and open water surfaces enables us to better understand the physical processes of air-water interaction. In context of measurements from the Humidity Exchange over the Sea Main Experiment （HEXMAX）, we recently proposed wave-parameter dependent approaches to sea surface friction velocity and the aerodynamic roughness by using the dimensional analysis method. To extend the application of these approaches to a range of natural surface conditions, the present study is to assess this approach by using both coastal shallow （RASEX） and open water surface measurements （Lake Ontario and Grand Banks ERS-1 SAR） where wind speeds were greater than 6.44 m s-1. Friction velocities, the surface aerodynamic roughness, and the neutral drag coefficient estimated by these approaches under moderate wind conditions were compared with the measurements mentioned above. Results showed that the coefficients in these approaches for coastal shallow water surface differ from those for open water surfaces, and that the aerodynamic roughness length in terms of wave age or significant wave height should be treated differently for coastal shallow and open water surfaces.

Test of the Louis scheme and COARE algorithm for friction velocity in different wind-sea/swell regimes

The Louis scheme and the COARE algorithm （version 3.0） are tested against eddy covariance and inertial dissipation methods for friction velocity estimates in different wind-sea/swell regimes. Atmospheric forcing data, tabulated by Donelan et al. （1997.J Phys Oceanog, 27：2087-2099）, were collected from a mast on the foredeck of a SWATH （small water-plane area, twin hull） ship in deep sea off the State of Virginia during the surface wave dynamics experiment. These data are representative of low to moderate wind regimes. The aerodynamic roughness length is determined by using the Charnock relationship. The intercomparison shows that the Louis scheme and the COARE algorithm underestimate the friction velocity by 6% and 3% respectively under pure wind sea conditions, 15% and 13% respectively under cross swell conditions, and 21% and 17% respectively under counter swell conditions. The analysis shows that these underestimations were caused by the method chosen to determine the aerodynamic roughness length because it significantly underestimates the aerodynamic roughness length. It is especially true under the cross swell and counter swell conditions.

Parameterizations of the Daytime Friction Velocity, Temperature Scale, and Upslope Flow over Gently Inclined Terrain in Calm Synoptic Conditions

A set of new parameterizations for the friction velocity and temperature scale over gently sloped terrain and in calm synoptic conditions are theoretically derived. The friction velocity is found to be proportional to the product of the square root of the total accumulated heating in the boundary layer and the sinusoidal function of the slope angle, while the temperature scale is proportional to the product of the boundary layer depth, the sinusoidal function of the slope angle and the potential temperature gradient in the free atmosphere. Using the new friction velocity parameterization, together with a parameterization of eddy diffusivity and an initial potential temperature profile around sunrise, an improved parameterization for the thermally induced upslope flow profile is derived by solving the Prandtl equations. The upslope flow profile is found to be simply proportional to the friction velocity.

Internal structures and high-velocity frictional properties of a bedding-parallel carbonate fault at Xiaojiaqiao outcrop activated by the 2008 Wenchuan earthquake

This paper reports internal structures of a bedding-parallel fault in Permian limestone at Xiaoji-aqiao outcrop that was moved by about 0.5 m during the 2008 MW7.9 Wenchuan earthquake. The fault is located about 3 km to the south from the middle part of Yingxiu-Beichuan fault, a major fault in the Longmenshan fault system that was moved during the earthquake. The outcrop is also located at Anxian transfer zone between the northern and central segments of Yingxiu-Beichuan fault where fault system is complex. Thus the fault is an example of subsidiary faults activated by Wenchuan earthquake. The fault has a strike of 243°or N63°E and a dip of 38°NW and is nearly optimally oriented for thrust motion, in contrast to high-angle coseismic faults at most places. Surface outcrop and two shallow drilling studies reveal that the fault zone is several centimeters wide at most and that the coseismic slip zone during Wenchuan earthquake is about 1 mm thick. Fault zone contains foliated cataclasite, fault breccia, black gouge and yellowish gouge. Many clasts of foliated cataclasite and black gouge contained in fault breccia indicate multiple slip events along this fault. But fossils on both sides of fault do not indicate clear age difference and overall displacement along this fault should not be large. We also report results from high-velocity friction experiments conducted on yellowish gouge from the fault zone using a rotary shear low to high-velocity frictional testing apparatus. Dry experiments at normal stresses of 0.4 to 1.8 MPa and at slip rates of 0.08 to 1.35 m/s reveal dramatic slip weakening from the peak friction coeffcient of around 0.6 to very low steady-state friction coeffcient of 0.1–0.2. Slip weakening parameters of this carbonate fault zone are similar to those of clayey fault gouge from Yingxiu-Beichuan fault at Hongkou outcrop and from Pingxi fault zone. Our experimental result will provide a condition for triggering movement of subsidiary faults or off-fault damage during a large earthquake.

Study on the Log-Linear Velocity Profile of Near-Bed Tidal Currentin Estuarine and Coastal Waters

Many observed data show that the near-bed tidal velocity profile deviates from the usual logarithmic law. The amount of deviation may not be large, but it results in large errors when the logarithmic velocity profile is used to calculate the bed roughness height and friction velocity （or shear stress）. Based on their investigation, Kuo et al. （1996） indicate that the deviation amplitude may exceed 100%. On the basis of fluid dynamic principle, the profile of the near-bed tidal velocity in estuarine and coastal waters is established by introducing Prandtl＇ s mixing length theory and Von Kannan selfsimilarity theory. By the fitting and calculation of the near-bed velocity profde data observed in the west Solent, England, the results are compared with those of the usual logarithmic model, and it is shown that the present near-bed tidal velocity profile model has such advantages as higher fitting precision, and better inner consistency between the roughness height and friction velocity. The calculated roughness height and friction velocity are closer to reality. The conclusions are validated that the logarithmic model underestimates the roughness height and friction velocity during tidal acceleration and overestimates them during tidal deceleration.

Study on the double-logarithmic profile of tidal flow velocity in the near-bed layers

Tidal current velocity profile in the near-bed layers has been widely studied. The results showed that velocity profile in the near-bed layer obviously departure from the traditional logarithmic profile, due to the acceleration or deceleration. Although the logarithmic linear profile can reduce the rate of deviation from this, only it is a lower-order approximate solution. In this paper, considering the unsteady and non-linear features of tidal motion, the double logarithmic profile near-bed layers in estuarine and coastal waters is established on the assumption that the turbulent shear stress along the water depth was parabolic distribution, and on the basis of Prandtl＇s mixing length theory and yon Karman＇s self-similar theory. Having been verified the data observed at the West Solent in the south of England, and comparison of the logarithmic linear profile, it found that the double logarithmic profile is more precious than the latter. At last, the discussed results showed that： （1） The parabolic distribution of the tidal shear stresses verified good by the field data and experimental data, can be better reflected the basic features of the tidal shear stress deviating from linear distribution that is downward when to accelerate, upward when to decelerate. （2） The traditional logarithmic velocity profile is the zero-order approximation solution of the double logarithmic profile, the logarithmic linear profile is the first order, and the logarithmic parabolic profile is the second order. （3） Ignoring the conditions of diffusion and convection in the tida movement, the double logarithmic profile can reflect the tidal properties of acceleration or deceleration, so that the calculation of the friction velocity and roughness length are more reasonable. When the acceleration or the deceleration is about zero, the double logarithmic profile becomes the logarithmic profile.

An approach to determine coeffi cients of logarithmic velocity vertical profile in the bottom boundary layer

Velocity vertical profiles in the bottom boundary layer are important to understand the oceanic circulation.The logarithmic vertical profile,u=A ln z+B,is the universal profile for the horizontal velocity in the boundary layer,in which two coefficients(A and B)need to be determined.The two coefficients are the functions of the friction velocity(u_(*))and the roughness length(z_(0)),and they are calculated using u_(*)and z_(0).However,the measurement of u_(*)and z_(0) is a challenge.In the present study,an approach is developed to estimate the two coefficients(A and B)by using a series of fl ume laboratory experiments with fl at boundary and regularly distributed cylinders as the rough boundaries.An acoustic doppler velocimeter(ADV)is used to measure the velocity vertical profiles of the steady flow.Using the measured velocity data,the regressed logarithmic profiles are obtained.Based on the series of the A and B values,the mathematical formula for A and B are statistically established as the function of the cylinder height,inflow velocity,and the water depth,which avoids the measurement of the friction velocity and the roughness length.

Preparation and Properties of Friction Materials by Using Two Kinds of Fibrous Industrial Minerals

The basic technology and properties of the brake blocks made of modified needle-like wollastonite and fibrous sepiolite were intensively researched.The impact strengthes and fixed velocity friction of the brake blocks prepared by different recipes were tested. The testing results show that it is feasible for needle-like wollastonite and fibrous sepiolite to take the place of asbestos as the reinforced materials of friction materials.The braking effect of the brake blocks is the best when the ratio of the needle-like wollastonite to the fibrous sepiolite was 1∶6.

Rolling velocity and relative motion of particle detector in local granular flow

The velocity of a particle detector in granular flow can be regarded as the combination of rolling and sliding velocities.The study of the contribution of rolling velocity and sliding velocity provides a new explanation to the relative motion between the detector and the local granular flow.In this study,a spherical detector using embedded inertial navigation technology is placed in the chute granular flow to study the movement of the detector relative to the granular flow.It is shown by particle image velocimetry(PIV)that the velocity of chute granular flow conforms to Silbert’s formula.And the velocity of the detector is greater than that of the granular flow around it.By decomposing the velocity into sliding and rolling velocity,it is indicated that the movement of the detector relative to the granular flow is mainly caused by rolling.The rolling detail shown by DEM simulation leads to two potential mechanisms based on the position and drive of the detector.

Sand flux estimation during a sand-dust storm at Tazhong area of Taklimakan Desert, China

In this paper, the sand transport during a sand-dust storm in the Tazhong area of the central Taklimakan Desert from 11：29 to 23：56 on July 19, 2008 was observed and measured in real time. The sand flux at Tazhong was estimated using sand transport empirical formulas. The critical friction velocity at Tazhong was 0.24 m/s and the functional relation between the wind speed and sediment discharge at the height of 2 m was established. It was also found that the calculated values by Lettau＇s sediment discharge formula were close to those of the instrument measurements. The horizontal sand flux and the vertical sand flux during this sand-dust storm at Tazhong were respectively 258.67×10-4 kg/（m·s） and 40.07×10-7 kg/（m2·s）.

Wind speed scaling and the drag coefficient

Wind speed scaling in similarity law in wind-generated waves and the drag coefficient are studied. In analyzing the data in the wind wave channel, it is found that the u＊ scaling greatly reduces the scatter in the U10 scaling. The u＊ scaling has much less scatter than the scaling using other wind speeds. The friction velocity seems to play a distinctive role in wave growth. The result is important in the applications of the similarity law and in wave modeling. In theory it gives an insight into the mechanism of wind wave interaction. It is found that wave steepness is important in influencing the drag coefficient. The variability of the coefficients in the currently widely used drag form can be explained by the differences in wave steepness in the observations. A drag coefficient model with wind speed and wave steepness as parameters is proposed. An explanation for Kahma＇s result that the u, scaling does not reduce the scatter in the U10 scaling is given.

Size distribution of PM_(20) observed to the north of the Tibetan Plateau

The size distribution of airborne dust particles is an important parameter in the measurement of dust emissions due to wind erosion,and a quantitative and accurate description is necessary.Observations regarding the size distribution of airborne dust particles are currently lacking in Tibetan Plateau(TP).This study aims to obtain a dust particle size distribution and compare the difference in spatial distribution at a field site(FS)and two urban observation sites(Minfeng and Hetian),in the north of the TP,under the condition of the dust-days and clear-days.The observation data was collected from 5 July to 4 August,2019.The mass concentration of dust(PM20)was measured with a 10-stage quartz crystal microbalance cascade impactor.The results suggested that 91.9%of the particles were less than 1.4μm on clear-days,and particles of d≥1.4μm increased to 27.2%with the occurrence of the dust events.More than 80%of the airborne dust particles were less than 1.4μm for each friction velocity.The proportion of d<0.7μm was the greatest at 3.0 m(above the surface,similarly hereinafter)on dust-days,while the proportion of d>0.7μm was the greatest at 0.5 m on clear-days.During urban observation,with the occurrence of dust events,the concentration of a particle size greater than 2.5μm increased most significantly in Minfeng,while the concentration of a particle size less than 0.7μm increased most significantly in Hetian.Moreover,the proportion of particles with d<0.7μm had the smallest difference,while the proportion of particles with d>1.4μm had the most obvious difference among these three sites.

Effect of stones on the sand saltation threshold during natural sand and dust storms in a stony desert in Tsogt-Ovoo in the Gobi Desert,Mongolia

Non-erodible elements such as stones and vegetation are key to controlling wind erosion and dust emission in drylands.Stony deserts are widely distributed in the Gobi Desert,but the effect of stones on wind erosion and dust emission have not been well studied,except under artificial conditions.In this study,we evaluated the effect of stones on wind erosion and dust emission by measuring the sand saltation threshold in a stony desert in Tsogt-Ovoo in the Gobi Desert,Mongolia,under natural surface conditions during sand and dust storms.We quantified the amount of stones by measuring the roughness density,and determined the threshold friction velocity for sand saltation by measuring wind speed and sand saltation count.Our results showed that the threshold friction velocity increased with the roughness density of stones.In the northern part of the study area,where neither a surface crust nor vegetation was observed,the roughness density of stones was 0.000 in a topographic depression(TD),0.050 on a northern slope(N.SL),and 0.160 on the northern mountain(N.MT).The mean threshold friction velocity values were 0.23,0.41,and 0.57 m/s at the TD,N.SL,and N.MT sites,respectively.In the southern part of the study area,the roughness density values of stones were 0.000 and 0.070-0.320 at the TD and southern slope sites,respectively,and the mean threshold friction velocities were 0.23 and 0.45-0.71 m/s,respectively.We further compared the observed threshold friction velocities with simulated threshold friction velocities using Raupach's theoretical roughness correction and the measured roughness density values,and found that Raupach's roughness correction worked very well in the simulation of threshold friction velocity in the stony desert.This means that the results of our stone measurement can be applied to a numerical dust model.

Numerical analysis of supersonic jet flow and dust transport induced by air ingress in a fusion reactor

During a loss of vacuum accident(LOVA),the air ingress into a vacuum vessel(VV)may lead to radioactive dust resuspension,migration,and even explosion,thereby posing a great threat to the safe operation of future fusion reactors;thus,it is crucial to understand the flow characteristics and radioactive dust transport behavior induced by LOVA.However,only a few studies have identified the characteristics of the highly under-expanded jet flow at a scale of milliseconds during LOVA.Particularly,the occurrence and behavior of a Mach disk is yet to be captured in existing studies.In this study,we used a more advanced model with a finer mesh and adaptive mesh strategies to capture the Mach disk in a VV during LOVA.In detail,a computational fluid dynamics–discrete phase model one-way coupled multiphase approach was established using the computational fluid dynamics code ANSYS FLUENT and applied to the analysis during the first seconds of LOVA.The results showed that air ingress into the VV behaved like a highly free under-expanded jet at the initial stage and Mach disk was formed at~6 ms.Moreover,the flow field dramatically changed at the position of the Mach disk.The jet core before the Mach disk had a maximum velocity of~8 Mach with the corresponding lowest static pressure(~100 Pa)and temperature(few tens of K).The friction velocities in the lower part of the VV,which is an area of concern due to dust deposition,were generally larger than 15 m/s near the inlet region.Lastly,the crude prediction of the particle trajectories demonstrated the spiral trajectories of the dust following the air motion.Therefore,this study provided a basis for further safety analysis and accident prevention related to dust transport and explosion in future fusion reactors.

Dependence of Estimating Whitecap Coverage on Currents and Swells

The shipboard measurements of whitecap coverage(W)and the meteorological and oceanographic information from two cruises in the South China Sea and Western Pacific are explored for estimating W.This study aims at evaluating how to im-prove the parameterizations of W while considering the effects of currents and swells on wave breaking.Currents indeed affect W in a way that winds with following currents can decrease W,whereas winds with opposing currents can increase W.Then,10-m wind speed over sea surface(U_(10))is calibrated by subtracting the current velocity from U_(10) when the propagating directions of winds and currents are aligned.By contrast,when the direction is opposite,U_(10) is calibrated by adding the parallel velocity com-ponent of the opposing current to U_(10).The power fits of W dependence on the U_(10)-related parameters of U_(10),friction velocity,wind sea Reynolds number in terms of this calibrated-U_(10) obtain better results than those directly fitted to U_(10).Different from the effect of currents on W,wind blowing along the crest line of swells may contribute to the increase in W.The conclusions suggest that U_(10) should be calibrated first before parameterizing W in areas with a strong current or some swell-dominant areas.

Parameterization of Wave Breaking Probability and Whitecap Coverage

The model for whitecap coverage and wave breaking probability are parameterized by the dimensionless wind fetch X^-. This paper aims at replacing X^- with other parameters such as the average wave period T^-, wind speed U10 or wave age ξ in order to improve the suitability and convenience of the model for application. First, W and B are expressed in terms of T^- and U10, which are relatively easy to measure in the field. Further, U10 is replaced with the friction velocity U. by use of the empirical relationship. As wave age has been widely used to parameterize spectral models of ocean waves and air-sea fluxes, W and B are then expressed as a simple function of wave age, respectively. The new forms of the model obtained are W= 1 - Ф（3.02ξ0＂76） and B = exp（ - 4.54ξ^1.52） . The two forms are mere applicable in pracrice, since ξ is relatively easy to measure or determine from wave and wind records. Comparisons between these expressions and data collected from published literature are made and agreement is fairly good.

Characteristic Study of the Boundary Layer Parameters over the Arabian Sea and the Bay of Bengal Using the QuikSCAT Dataset

The marine atmospheric boundary layer （MABL） plays a vital role in the transport of momentum and heat from the surface of the ocean into the atmosphere. A detailed study on the MABL characteristics was carried out using high-resolution surface-wind data as measured by the QuikSCAT （Quick scatterometer） satellite. Spatial variations in the surface wind, frictional velocity, roughness parameter and drag coefficient for the different seasons were studied. The surface wind was strong during the southwest monsoon season due to the modulation induced by the Low Level Jetstream. The drag coefficient was larger during this season, due to the strong winds and was lower during the winter months. The spatial variations in the frictional velocity over the seas was small during the post-monsoon season （-0.2 m s^-1）. The maximum spatial variation in the frictional velocity was found over the south Arabian Sea （0.3 to 0.5 m s^-1） during the southwest monsoon period, followed by the pre-monsoon over the Bay of Bengal （0.1 to 0.25 m s^-1）. The mean wind-stress curl during the winter was positive over the equatorial region, with a maximum value of 1.5×10^-7 N m^-3, but on either side of the equatorial belt, a negative wind-stress curl dominated. The area average of the frictional velocity and drag coefficient over the Arabian Sea and Bay of Bengal were also studied. The values of frictional velocity shows a variability that is similar to the intraseasonal oscillation （ISO） and this was confirmed via wavelet analysis. In the case of the drag coefficient, the prominent oscillations were ISO and quasi-biweekly mode （QBM）. The interrelationship between the drag coefficient and the frictional velocity with wind speed in both the Arabian Sea and the Bay of Bengal was also studied.

Dynamics analysis for flux of carbon dioxide across seasurface

The CO2 exchange between atmosphere and seawater has ho discussed. The physical process of exchange at seasurface is related closely with the atmospheric forcing. The result of a dynamics analysis of .the CO2 flux across the seasurface, and a computed example over the shelf area of the East China Sea are reported. The analysis and calculation show that the CO2 vertical flux at seawater is mainly dependent on the CO2 concentration difference between seawater and atmosphere. Its exchange velocity is affected by the friction velocity, solubility and resistivity. The agreement between this exchange velocity and recent wind tunnel experimental results is good with wind speed over the range4～ 10m/s. This flux can increase very fast with the increase of wind speed, but very small with the increase oftemperature. During the three cruise expeditions in October 1993, and April and October 1994, the CO2 vertical flux at seasurface was directed downward over the shelf area of the East China Sea. Because the CO2 Partial pressure difference between seawater and atmosphere was negative value, the flux value was - 45. 5μg/(m2.h) in average.

Numerical analysis of submarine landslides using a smoothed particle hydrodynamics depth integral model

Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Currently, commercial calculation programs such as BING have limitations in simulating underwater soil movements. All of these processes can be consistently simulated through a smoothed particle hydrodynamics（SPH） depth integrated model. The basis of the model is a control equation that was developed to take into account the effects of soil consolidation and erosion. In this work, the frictional rheological mode has been used to perform a simulation study of submarine landslides. Time-history curves of the sliding body＇s velocity, height,and length under various conditions of water depth, slope gradient, contact friction coefficient, and erosion rate are compared; the maximum sliding distance and velocity are calculated; and patterns of variation are discussed.The findings of this study can provide a reference for disaster warnings and pipeline route selection.

AN APPROACH TO SEPARATING THE CURRENT VELOCITY IN ESTUARINE AND COASTAL WATERS

In the near-shore waters, the actual flow is mainly induced by tide, wind and salinity, and the river water runoff should also be included as a component in the estuary waters. The interactions among these major components are very complicated. Many approaches were proposed to study isolated tide and wind-driven currents or run-off based on the measured velocity, with all its components taken as a whole. In this article, firstly, based on the actual hydrodynamic characteristics of estuarine and coastal waters, an approach is proposed to separate the measured velocity by considering the theoretical current velocity profiles and using the least squares method. The vertical structures of tidal, wind-driven currents, density current and runoff can be obtained as well as their proportions in the measured velocity. Then, this approach is applied to the analysis of velocity data obtained in the North Branch of Yangtze River estuary and of laboratory test data. The results are found to be satisfactory. Finally, this approach is used to separate the measured velocity in the South Branch of Yangtze River estuary, to determine not only the bed friction velocity and roughness height, but also the surface wind stress, and to estimate the wind velocity data above the water surface. The results show that this method is simple in principle, practical in use, and reasonable in obtained results. So it can be used to effectively analyze the field data.