On the bottom relief effect on wind-driven currents in a homogeneous ocean

The three-dimensional model of stationary wind-driven currents in a homogeneous ocean of a variable depth is investigated. The model is linear but includes horizontal and vertical turbulent mixings. Two cases of the behaviour of the isolines of the function $\text{?/H}$ are considered, namely: (1) all isolines $\text{?/H}$ start at one part of the coastline and end in another part of it, and (2) a certain isoline $\text{?/H}$ exists which is tangential to the coastline. Here ? is the Coriolis parameter, and H is the depth of the ocean. The first case is the simplest one; it arises in particular if H = constant and the coasts are meridional. The second case is marked by the boundary current separation from the coast. The paper deals with the boundary layers which arise at the surface, bottom, side boundary and inside the ocean.     

Forward-in-time and Backward-in-time Dispersion in the Convective Boundary Layer: the Concentration Footprint

The turbulence field obtained using a large-eddy simulation model is used to simulate particle dispersion in the convective boundary layer with both forward-in-time and backward-in-time modes. A Lagrangian stochastic model is used to treat subgrid-scale turbulence. Results of forward dispersion match both laboratory experiments and previous numerical studies for different release heights in the convective boundary layer. Results obtained from backward dispersion show obvious asymmetry when directly compared to results from forward dispersion. However, a direct comparison of forward and backward dispersion has no apparent physical meaning and might be misleading. Results of backward dispersion can be interpreted as three-dimensional or generalized concentration footprints, which indicate that sources in the entire boundary layer, not only sources at the surface, may influence a concentration measurement at a point. Footprints at four source heights in the convective boundary layer corresponding to four receptors are derived using forward and backward dispersion methods. The agreement among footprints derived with forward and backward methods illustrates the equivalence between both approaches. The paper shows explicitly that Lagrangian simulations can yield identical footprints using forward and backward methods in horizontally homogeneous turbulence.     

城市边界层数值模式研究以及在香港地区复杂地形下的应用

由于城市表面复杂的下垫面的影响,使得城市边界层风温场结构较其他下垫面有很大不同.作者通过将城市中500 m网格内的各种地表类型,按照各自在网格中所占的百分比及各自的地表参数加权平均,得到此网格的平均的地表参数,以此准确反映下垫面的情况,建立分辨率为500 m的城市边界层能量平衡模式,将此城市边界层能量平衡模式嵌入动力学框架,并用中尺度模式MM5作为初始条件和边界条件,建立一个既考虑中尺度背景场又详细考虑城市下垫面复杂性与多样性的城市边界层模式系统.将模式系统运用于香港复杂地形下的边界层特点的模拟研究.通过与观测值的比较,模式能够较准确的模拟出海陆风、城市热岛等热力过程,及气流过山引起的绕流等动力过程,并且通过对边界层高度的模拟预测污染扩散的条件等.说明模式系统具有模拟在中尺度的背景场的控制下海陆风环流、过山堆积和绕流及城市热力影响的能力.     

A model for the height of the internal boundary layer over an area with an irregular coastline

A model for the time and space variation of the internal boundary-layer height over a land area with an irregular coastline is presented. It is based on the analytical model of the boundary-layer height proposed by Gryning and Batchvarova (1990) and Batchvarova and Gryning (1991), The model accounts for the temperature jump and the mean vertical air motion at the top of the internal boundary-layer. Four cases from experiments in Nanticoke and Vancouver are used for model validation. The agreement between the calculated and measured internal boundary layer height at the observational sites is fairly good. The input information for the model consist of wind speed and direction, friction velocity and kinematic heat flux in time and space for the area, and the potential temperature gradient and the mean vertical air motion above the internal boundary layer. For the experiments used in the validation the effect of subsidence is relatively important in the afternoon under low wind speed high pressure conditions, lowering the height of the internal boundary layer by up to 10%, and it is negligible in the morning hours. The effect of the mixing height over the sea is found to be negligible.     

On the transient adjustment of a mid-latitude abyssal ocean basin with realistic geometry: the constant depth limit

The early stages in the adjustment of a mid-latitude abyssal basin with realistic geometry are studied using an inverted one and one-half layer model of the Eastern Mediterranean Sea as a natural test basin. The model is forced with a localized sidewall mass source and a compensating distributed mass sink. A flat bottom basin is investigated for comparison with existing theories on abyssal gyral spin-up, and as a precursor to a study with realistic topography. As in existing theories, the early adjustment is dominated by sub-inertial Kelvin and Rossby waves. Obstacles and the varying coastal geometry do not impede the passage of the Kelvin wave, though the circuit time of the main Kelvin wave signal is reduced by an aggregate 6% for the abyssal Eastern Mediterranean basin. The scattering of the Kelvin wave due to small-scale variations in the coastline is also shown not to be significant to the adjustment. The relatively short period of time needed to reach a statistical steady state is attributed to western boundary current formation in response to local Kelvin wave dynamics. Upon cessation of the sidewall forcing, sub-inertial motion controls the spin-down adjustment with basin-scale Rossby waves becoming the most pronounced feature of the flow. Two dynamical issues of particular interest emerge in these simulations: the retardation of Kelvin wave propagation around the abyssal basin and the roles of detrainment and sidewall forcing in the interior vorticity balance. An idealized simulation using an elliptical basin is used to illustrate that the mechanism for Kelvin wave retardation is a geometrically induced dispersion due to large-scale variations in the coastline. A dynamical analysis of the interior circulation shows that detrainment alone does not develop a Sverdrup response. Both the localized sidewall injection and the detrainment are needed to describe the interior dynamics, with both poleward and equatorward flows developing during the adjustment.     

A Homogeneous Langevin Equation Model, Part ii: Simulation of Dispersion in the Convective Boundary Layer

We present a Lagrangian stochastic model of vertical dispersion in the convective boundary layer (CBL). This model is based on a generalized Langevin equation that uses the simplifying assumption that the skewed vertical velocity probability distribution is spatially homogeneous. This approach has been shown to account for two key properties of CBL turbulence associated with large-scale coherent turbulent structures: skewed vertical velocity distributions and long velocity correlation time. A 'linear-skewed' form of the generalized Langevin equation is used, which has a linear (in velocity) deterministic acceleration and a skewed random acceleration. 'Reflection' boundary conditions for selecting a new velocity for a particle that encounters a boundary were investigated, including alternatives to the standard assumption that the magnitudes of the particle incident and reflected velocities are positively correlated. Model simulations were tested using cases for which exact, analytic statistical properties of particle velocity and position are known, i.e., well-mixed spatial and velocity distributions. Simulations of laboratory experiments of CBL dispersion show that (1) the homogeneous linear-skewed Langevin equation model (as well as an alternative 'nonlinear-Gaussian' Langevin equation model) can simulate the important aspects of dispersion in the CBL, and (2) a negatively-correlated-speed reflection boundary condition simulates the observed dispersion of material near the surface in the CBL significantly better than alternative reflection boundary conditions. The homogeneous linear-skewed Langevin equation model has the advantage that it is computationally more efficient than the homogeneous nonlinear-Gaussian Langevin equation model, and considerably more efficient than inhomogeneous Langevin equation models.     

Large-Eddy Simulation of air Pollution Dispersion in the Nocturnal Cloud-Topped Atmospheric Boundary Layer

Effects of stratocumulus clouds on the dispersion of contaminants are studied in the nocturnal atmospheric boundary layer. The study is based on a large-eddy simulation (LES) model with a bulk parametrization of clouds. Computations include Lagrangian calculations of atmospheric dispersion of a passive tracer released from point sources at various heights above the ground. The results obtained show that the vertical diffusion is non-Gaussian and depends on the location of a source in the boundary layer.     

Polar coordinates shallow water storm surge model for the coast of Bangladesh

A coastal shallow water model generally involves at least one coastal boundary which is curvilinear in nature. Representation of this curvilinear boundary needs special attention when a finite difference scheme is used. To incorporate the curvilinear coastal boundary and off-shore islands properly in the numerical scheme, we need fine resolution which may not be necessary away from the coast. In order to tackle this problem, a nested numerical scheme may be used where fine resolution for the coastal boundary region is nested into a course resolution scheme covering the whole analysis area. But this involves comparatively more memory and CPU time in the solution process. The present study aims at developing a storm surge model that may subsequently be modified for operational forecasting purpose for the coast of Bangladesh. A vertically integrated model is developed in a cylindrical polar coordinate system capable of incorporating bending of the coastline and off-shore islands with considerable accuracy without using any nesting. A comparison of CPU time and memory in the numerical computations between two types of modelling is described and it is found that the cylindrical polar coordinates model is economic. The model is applied for estimating the water levels at different coastal and island stations associated with a few storms that hit the coast of Bangladesh. It is found that the computed and observed water levels are in good agreement. Considering the computing time and other factors, it is found that the Polar Coordinates model is more suitable for operational forecasting purpose along the coast of Bangladesh.     

基于离海距和GIS技术的福建低温精细监测

利用福建省1:250000的DEM资料和67个气象站气温观测数据,在建立最低气温与经纬度、海拔高度相关推算方程的基础上, 融合离海距因子, 对2008—2010年冬季3个冷空气过程的最低气温 (数值及分布状态) 进行精细模拟,同时总结出利用逐步回归及综合残差平方和选取适宜离海距的方法。结果表明:融合离海距因子后,对冷空气过程最低气温的模拟效果更好。随着过程平均降温幅度的增大,离海距对过程最低气温模拟值的贡献率有减小趋势。不同冷空气过程的离海距大小存在差异,总体上以50 km为标准,再进一步得出适宜离海距。离海距以外区域最低气温模拟适用经度、纬度、海拔高度3个因子确定的地理气候方程进行,以内区域则适用在上述最低气温模拟方法的基础上融合离海距因子进行,以达到提高低温监测模拟精度和体现海洋对陆地温度调节能力的目的。经检验,模拟结果与实际情况基本相符。     

城市环境对重气扩散的影响及其验证

开发了适合城市环境下的应急重气扩散模型——SLAB_URBAN模型,该模型能够对城市环境下重气的传输扩散过程进行模拟。模型的原理基于重气扩散浅层理论,采用了新的城市边界层和扩散参数的参数化方案。该方案考虑了城市冠层内特有的风和湍流扩散的特征,能够体现城市边界层和湍流对重气扩散的影响。对美国盐湖城Urban2000的城市扩散试验进行模拟,主要验证下风方向观测弧所观测到的气体最大小时平均浓度与源释放速率的比值。结果表明,模型能够比较好地模拟出下风方向上浓度的分布特征。另外,与国外同类城市扩散模型的比较来看,SLAB_URBAN模型的模拟能力居于前列。     

一个预测沿岸陆上熏烟扩散的随机游动模式

本文建立了一个对流边界层中的随机游动扩散模式,并用KNRC试验No.64的资料作了验证.然后考虑混合层顶不规则夹卷界面的作用,应用于楠蒂科克沿岸陆上熏烟扩散问题,将模拟的地面轴线浓度与观测资料以及三个熏烟扩散预测模式的结果作了比较.结果表明:模式的模拟效能良好.所有数值计算均可在微机上实现,便于推广应用,更具实用性.     

A Two-Dimensional Numerical Study of the Impact of a City on Atmospheric Circulation and Pollutant Dispersion in a Coastal Environment

The urban impact on the sea breeze is studied by means of a mesoscale model with a detailed urban parameterisation. Four simulations are carried out on an idealised two-dimensional flat domain. In the base case, half of the domain is characterised by seaand the other half by rural land. In the urban case, an urban area 10 km wide is added near the shoreline. Simulations are performed for a moist rural soil (weak sea breeze) and for a dry rural soil (strong sea breeze). Results are analysed in order to evaluate the impact of the city on the wind, temperature and turbulent kinetic energy fields. The dispersion of a passive tracer emitted near the coastline is, also, used in the comparison. Results show that the city accelerates the sea-breeze formation in the morning (combinations of urban circulation and sea breeze), but it slows thesea-breeze front penetration. Moreover, the presence of the city enhances the recirculation processes and strongly modifies the pollutant dispersion. These effects are enhanced for a moist rural soil.     

Turbulent dispersion in coastal atmospheric boundary layers: An application of a Lagrangian model

A Lagrangian model is applied to simulate the dispersion of passive tracers (in particular, water vapour) in coastal atmospheric boundary layers under onshore wind conditions. When applied to convective boundary layers over uniform surfaces, the model gives results in agreement with those of similar studies. Numerical simulation of turbulent dispersion in coastal areas also reproduces the basic features known from experimental studies. Under onshore wind conditions, the humidity field is plume-shaped with the maximum vertical transport being over land downstream of the coast line. The model shows that the surface sensible heat flux over land, the static stability of the onshore air flow and the onshore wind speed are the most important factors determining the basic features of turbulent dispersion in coastal areas.     

Dispersion of passive tracer in the atmospheric convective boundary layer with wind shears: a review of laboratory and numerical model studies

Summary Paper reviews recent laboratory and numerical model studies of passive gaseous tracer dispersion in the atmospheric convective boundary layer (CBL) with surface and elevated wind shears. Atmospheric measurement data used for validation of these two model techniques are briefly discussed as well. A historical overview is given of laboratory studies of dispersion in the atmospheric CBL. Model studies of tracer dispersion in two CBL types, the (i) non-steady, horizontally homogeneous CBL and (ii) quasi-stationary, horizontally heterogeneous CBL, are reviewed. The discussion is focused on the dispersion of non-buoyant plume emitted from a point source located at different elevations within the CBL. Approaches towards CBL modeling employed in different laboratory facilities (water tanks and wind tunnels) are described. The reviewed numerical techniques include Large Eddy Simulation (LES) and Lagrangian modeling. Numerical data on dispersion in the sheared CBL is analyzed in conjunction with experimental results from wind-tunnel CBLs.     

Lagrangian Stochastic Modeling of Dispersion in the Stable Boundary Layer

Lagrangian stochastic models are well-suited for modeling dispersion in the stable boundary layer, especially in complex terrain. This note briefly describes the formulations and application of a Lagrangian stochastic model to predict dispersion of tracers released within nocturnal drainage flows.     

Centennial oscillations in an ocean-ice coupled model

Irregular centennial oscillations, with a spectral peak at 106 years, were obtained from an ocean-ice coupled model for the North Atlantic with realistic coastline and bottom topography. The model’s thermohaline circulation is forced by mixed boundary conditions, i.e., a Haney-type relaxation condition for temperature, but an equivalent virtual salt flux condition for salinity. All forcing fields are taken from the observed monthly mean climatological wind stress and buoyancy fluxes. The oscillations appeared in the form of a surface-intensified tripole in both the sea surface temperature and salinity fields located in the vicinity of the Labrador Sea. The oscillations involve a delicate interplay between heat and fresh water advection by meridional overturning circulation, horizontal gyres, vertical convection, and the seasonal cycle. The oscillations are primarily control?led by the salinity component of the circulation; however, sea ice plays a minor role in driving the oscillations observed in the model. On the other hand, a regular seasonal cycle in the forcing fields is an important ingredient for the centennial oscillations.     

Boundary-Layer Wind Structure in a Landfalling Tropical Cyclone

In this study, a slab boundary layer model with a constant depth is used to analyze the boundary-layer wind structure in a landfalling tropical cyclone. Asymmetry is found in both the tangential and radial components of horizontal wind in the tropical cyclone boundary layer at landfall. For a steady tropical cyclone on a straight coastline at landfall, the magnitude of the radial component is greater in the offshoreflow side and the tangential component is greater over the sea, slightly offshore, therefore the greater total wind speed occurs in the offshore-flow side over the sea. The budget analysis suggests that： （1） a greater surface friction over land produces a greater inflow and the nonlinear effect advects the maximum inflow downstream, and （2） a smaller surface friction over the sea makes the decrease of the tangential wind component less than that over land. Moreover, the boundary layer wind structures in a tropical cyclone are related to the locations of the tropical cyclone relative to the coastline due to the different surface frictions. During tropical cyclone landfall, the impact of rough terrain on the cyclone increases, so the magnitude of the radial component of wind speed increases in the offshore-flow side and the tangential component outside the radius of maximum wind speed decreases gradually.     

Impact of Land Surface Heterogeneity on Mesoscale Atmospheric Dispersion

Prior numerical modelling studies show that atmospheric dispersion is sensitive to surface heterogeneities, but past studies do not consider the impact of a realistic distribution of surface heterogeneities on mesoscale atmospheric dispersion. While these focussed on dispersion in the convective boundary layer, the present work also considers dispersion in the nocturnal boundary layer and above. Using a Lagrangian particle dispersion model (LPDM) coupled to the Eulerian Regional Atmospheric Modeling System (RAMS), the impact of topographic, vegetation, and soil moisture heterogeneities on daytime and nighttime atmospheric dispersion is examined. In addition, the sensitivity to the use of Moderate Resolution Imaging Spectroradiometer (MODIS)-derived spatial distributions of vegetation characteristics on atmospheric dispersion is also studied. The impact of vegetation and terrain heterogeneities on atmospheric dispersion is strongly modulated by soil moisture, with the nature of dispersion switching from non-Gaussian to near-Gaussian behaviour for wetter soils (fraction of saturation soil moisture content exceeding 40%). For drier soil moisture conditions, vegetation heterogeneity produces differential heating and the formation of mesoscale circulation patterns that are primarily responsible for non-Gaussian dispersion patterns. Nighttime dispersion is very sensitive to topographic, vegetation, soil moisture, and soil type heterogeneity and is distinctly non-Gaussian for heterogeneous land-surface conditions. Sensitivity studies show that soil type and vegetation heterogeneities have the most dramatic impact on atmospheric dispersion. To provide more skilful dispersion calculations, we recommend the utilisation of satellite-derived vegetation characteristics coupled with data assimilation techniques that constrain soil-vegetation-atmosphere transfer (SVAT) models to generate realistic spatial distributions of surface energy fluxes.     

Mathematical model for hermitized atmospheric dispersion in low winds with eddy diffusivities as linear functions of down wind distance

Summary In the present paper, an attempt is made for generalized the atmospheric diffusion operator. This can be accomplished by employing the realizability procedure, to identify a surface operator, that ensures self-adjointness’ of the atmospheric diffusion operator. The dispersion modeling in low wind speeds assumes importance because of the high frequency of occurrence and episodic nature of these poor diffusion conditions. A steady-state mathematical model for hermitized model has been calculated for the dispersion of air pollutants in low winds by taking into account the diffusion in the three coordinate directions and advection along the mean wind. The eddy diffusivities have been parameterized in terms of downwind distance for near source dispersion (Arya, 1995). The constants involved in this parameterization are the squares of intensities of turbulence. An analytical solution for resulting advection-diffusion equation with the physically relevant boundary conditions has been obtained. The solution has been used to simulate the field tracer data collected at IIT Delhi in low wind convective conditions.     

The Influence Of Vertical Wind Direction Shear On Dispersion In The Convective Boundary Layer, And Its Incorporation In Coastal Fumigation Models

The mean concentration distributionwithin a plume released from a point source in the atmosphericboundary layer can be greatly influenced by the systematic turningof wind with height (i.e. vertical wind direction shear). Such aninfluence includes a deflection of the plume centroid, with anassociated shearing of the vertical plume cross-section, and anenhancement of dispersion, in the horizontal plane. Wind directionshear is normally not accounted for in coastal fumigation models,although dispersion observations with shear acting as acontrolling parameter are not uncommon. A three-dimensionalLagrangian stochastic model is used to investigate the influenceof uniform wind direction shear on the diffusion of a point-sourceplume within the horizontally homogeneous convective boundarylayer, with the source located at the top of the boundary layer.Parameterisations are developed for the plume deflection andenhanced dispersion due to shear within the framework of aprobability density function (PDF) approach, and compared with theLagrangian model results. These parameterisations are thenincorporated into two applied coastal fumigation models: a PDFmodel, and a commonly used model that assumes uniform andinstantaneous mixing in the vertical direction. The PDF modelrepresents the vertical mixing process more realistically. A moreefficient version of the PDF model, which assumes a well-mixedconcentration distribution in the vertical at large times, isapplied to simulate sulfur dioxide data from the Kwinana CoastalFumigation Study. A comparison between the model results and thedata show that the model performs much better when the wind-sheareffects are included.