Sensitivity of the dynamics of volcanic eruption columns to their shape

This paper presents a one-dimensional steady-state model to investigate the sensitivity of the dynamics of sustained eruption columns to radius variations with height due to thermal expansion of the entrained air, and decreases in atmospheric pressure with height. In contrast to a number of previous models using an equation known as the entrainment assumption, the new model is based on similarity arguments to derive an equation set equivalent to the model proposed by Woods [Bull Volcanol 50:169–193, 1988]. This approach allows investigation of the effect of gas compressibility on the entrainment rate of ambient air, which has been little examined for a system in which a decrease in pressure significantly affects the density stratification of a compressible fluid. The new model provides results that include two end members: one in which the volume change within the eruption columns affects only the radial expansion without changing the vertical motion, and the other is the converse. The Woods [Bull Volcanol 50:169–193, 1988] model can be regarded as being between those two end members. The range of uncertainty arises because the extremely high temperature of discharged materials from a volcanic vent, and the exceptional terminal height of the eruption columns, allow significant expansion of the gas component in the eruption columns, making them behave differently from common turbulent plumes. This study indicates that the maximum height of the eruption columns is affected considerably by this uncertainty, particularly when the eruption columns extend above a height of 10 km, at which the pressure is about one-fourth the pressure at the ground surface. Column collapse may also be suppressed in wider parameter ranges than previously estimated. However, the uncertainty can be reduced by measuring column radii through a vertical profile during actual volcanic eruptions. Accordingly, this paper suggests that appropriate observation of eruption column shapes is essential for improving our understanding of the dynamics of eruption columns.     

Transport of Temperature and Humidity Variance and Covariance in the Marine Surface Layer

In this paper we address the budget of potential temperature T and moisture mixing ratio q variances as well as the q-T covariance budget. We focus on the vertical transport and study the quantities contained in these terms. Estimates of transport terms are rare and to the best of our knowledge they have never been shown together in the literature. For this purpose an extensive set of humidity and temperature turbulence data from two heights in a marine environment has been analysed. We discuss and compare our results with results from earlier studies for a range of stability from near-neutral to unstable. Similarity and differences between T and q behaviour in the surface layer over water are presented.     

A lidar study of the atmospheric entrainment zone and mixed layer over Hong Kong

In this paper we present the lidar study of the atmospheric boundary layer using the City University of Hong Kong lidar system. Three cases of determination of the entrainment zone thickness and the mixed layer depth in the atmospheric boundary layer over Hong Kong were selected for detailed study. The data collected have been analysed using the visual inspection method and the Steyn et al. [J. Atmos. Ocean. Technol. 16 (1999) 953] detection method. During the cold front passage, the mixed layer depth stayed at the level almost constant over the transition period, and the entrainment zone was thickening at a steady entrainment rate.     

Boundary-Layer Entrainment Estimation Through Assimilation of Radiosonde and Micrometeorological Data into a Mixed-Layer Model

In this study we estimate boundary-layer growth and entrainment bycombining radiosonde and micrometeorological observations with asimple coupled boundary-layer and land surface model. A variational(smoothing) approach is used to find the optimal estimate ofentrainment over the daytime window. This method is appealingbecause it accounts for the uncertainty in the various data streams,while enforcing the dynamics of the model (i.e., water and energybudgets in the boundary layer, mixed-layer growth, etc.). Thetraditional variational framework was modified in this study toinclude an ensemble approach, which not only yields a (mean) estimateof entrainment, but a measure of its uncertainty as well. Themethodology is applied to a field experiment site in Kansas. Resultsfrom this study indicate a much larger ratio of entrainment to surfacefluxes compared to early literature values from other sites. However,our results are consistent with recent estimates at the site usingindependent estimation methods. In tests where radiosonde data werewithheld, reasonable skill in entrainment estimation was still shown,suggesting the potential for more widespread applications where onlymicrometeorological data are available. Finally, the data assimilationframework presented here has not traditionally been used inatmospheric boundary-layer studies, and may provide a useful approachfor studying other aspects of the boundary layer in the future.     

Wind tunnel investigation of the influence of surface moisture content on the entrainment and erosion of beach sand by wind using sands from tropical humid coastal southern China

Wind erosion has major impacts on dune growth, desertification, and architecture on sea coasts. The deflation threshold shear velocity is a crucial parameter in predicting erosion, and surface moisture greatly affects this threshold and thus sand stability. Wind tunnel studies have shown that reduced moisture contents decrease entrainment thresholds and increase wind erosion, but field and wind tunnel test data is lacking for tropical humid coastal areas. In this study, we investigated the influence of surface moisture contents (at 1 mm depth) on sand entrainment and erosion using tropical humid coastal sands from southern China. Shear velocities were deduced from velocity profiles above the sand. The threshold shear velocity increased linearly with increasing ln100M (M, gravimetric moisture content). The increase was steepest below a moisture content of 0.0124 (i.e., at M_1.5, the moisture content in the sand at a matric potential of − 1.5 MPa). We compared several popular models that predict threshold shear velocity of moisture sediment, and found substantial differences between their predicted results. At a surface moisture content of 0.0124, the predicted increase in the wet threshold shear velocity compared with the dry threshold shear velocity ranged from 34% to 195%. The empirical model of Chepil and Selah simulated the data well for M < 0.0062 (i.e., 0.5M_1.5), whereas Belly's empirical model simulated the data best for M > 0.0062. Wind erosion modulus increased with increasing effective wind velocity following a power function with a positive exponent at all moisture contents, but decreased with increasing surface moisture content following a power function with a negative exponent. When wind speed and moisture content varied simultaneously, wind erosion modulus was proportional to the 0.73 power of effective wind velocity, but inversely proportional to the 1.48 power of M. The increase in resistance to erosion at low moisture contents probably results from cohesive forces in the water films surrounding the sand particles. At a moisture content near M_1.5, wind erosion ceases nearly for all wind velocities that we tested.     

Vertical mixing of oil droplets by breaking waves

Oil spilled on a sea surface can be dispersed by a variety of natural processes, of which the influence of breaking waves is dominant. Breaking waves are able to split the slick into small droplets, facilitating oil mixing in the water column. Vertical dynamics of the droplets plays a major role in the oil mass exchange between the slick and the water column. In this paper a mathematical model of oil droplet mixing by breaking waves is developed. The model uses a kinetic approach to describe the vertical exchange of the droplets at the interface between the slick and the water column. The majority of the coefficients and parameters are conveniently combined into a single "mixing factor". The model is verified using sensitivity analysis and empirical formulae of other authors. The model permits a rapid estimation of the amount of dispersed oil under the breaking waves. The ultimate goal of the research is to parameterise influence of breaking waves on vertical mixing of oil droplets to be used in a general 3-D oil spill model.     

Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-eddy Simulations. Part 1: Free Convection

Convection in a quasi-steady, cloud-free, shear-free atmospheric boundary layer is investigated based on a large-eddy simulation model. The performed tests indicate that the characteristic (peak) values of statistical moments at the top of the mixed layer are proportional to the interfacial scales (from gradients of scalars in the interfacial layer). Based on this finding a parameterization is proposed for profiles of scalar variances. The parameterization employs two, semi-empirical similarity functions F_m(z/z_i) andF_i(z/z_i), multiplied by a combination of the mixed-layer scales and the interfacial scales.     

Dynamics of medium-intensity dense water plumes in the Arkona Basin, Western Baltic Sea

In this study, the dynamics of medium-intensity inflow events over Drogden Sill into the Arkona Sea are investigated. Idealised model simulations carried out with the General Estuarine Transport Model suggest that most of the salt transport during such inflow events occur north of Kriegers Flak, a shoal with less than 20 m water depth surrounded by water depths of more than 40 m. This assumption about the pathway is supported by recent ship-based observations in the Arkona Sea during a medium-intensity inflow event. The propagation of a saline bottom plume could be observed during several days after having passed Drogden Sill. In the area north of Kriegers Flak the plume was about 10 m thick, and propagated with more than 0.5 m s⁻¹ and a salinity of up to 20 psu (with ambient water salinity being 8 psu) eastwards. Although the model simulations were idealised, the structural agreement between the observation and model result was good. The structure and pathways of these medium-intensity inflow events are of specific interest due to the plans for erecting extensive offshore wind farms in the Arkona Sea which may under certain circumstances lead to increased entrainment of ambient water into the bottom plumes.     

Diurnal cycle of precipitation over the Maritime Continent simulated by a spectral cumulus parameterization

The validity of a spectral cumulus parameterization (spectral scheme) for simulating a diurnal cycle of precipitation over the Maritime Continent (MC) was examined using a regional atmospheric model. The impacts of entrainment parameterization and each type of convective closure, i.e., non-equilibrium (or equilibrium) closure for deep convection, mid-level, and shallow convective closures, were also examined. When vertically variable entrainment and appropriate convective closures were employed, the model adequately simulated a diurnal cycle of precipitation over both land and ocean as compared to the observation. Analysis regarding the entrainment parameterization revealed that variable entrainment parameterization was needed not only for simulating better mean patterns of precipitation, but also for more realistic phases of diurnal cycles. The impacts of convective closures appeared in the differences in the precipitation amplitude. Analysis on diurnal cycles of convective properties and tendencies revealed that the cycles between boundary layer forcing and convective heating determined convection strength and were affected by each type of convective closure. It can be concluded that the spectral scheme with appropriate convective closures is able to simulate a realistic diurnal cycle over the MC.     

风切变对边界层对流影响的大涡模拟研究

利用"西北干旱区陆-气相互作用野外观测实验"加密观测期间在敦煌站的观测资料以及大涡模式,模拟了对流边界层的发展,以及示踪物从混合层向残留层传输的时空变化。模拟的对流边界层的结构及演变特征与实测结果基本一致。进一步通过有风切变和无风切变的敏感性数值试验,研究了风切变对垂直速度、位温和示踪物浓度的水平分布以及示踪物传输高度的影响。研究结果表明,在有风切变的试验中(甚至风切变仅存在于近地层中),对流边界层的增长加强,而且示踪物被传输的高度也较高。与浮力驱动的对流边界层相比,由浮力和风切变共同驱动的边界层中上升气流较弱而下沉气流较强,但前者的上升气流与下沉气流的分布在垂直方向上更为倾斜。由于夹卷作用的增强,浮力和风切变共同驱动的对流边界层较浮力驱动的对流边界层暖。在夹卷层,浮力和风切变共同驱动的边界层对流的上升气流和下沉气流都比浮力驱动的边界层对流中的强,而且垂直速度的概率密度函数分布也较对称,其位温和示踪物浓度的概率密度函数分布也比浮力驱动的边界层中的平直。对湍流动能收支的分析也表明风切变对湍流动能有重要影响,尤其对夹卷层中的湍流动能切变产生项影响较大。示踪物浓度的概率密度函数垂直分布显示,浮力驱动的边界层中示踪物浓度随高度变化较小,而浮力和风切变共同驱动的边界层中示踪物浓度随高度递减,但是示踪物传输的高度比较高。