The performance of RAMS in representing the convective boundary layer structure in a very steep valley

Data from a comprehensive field study in the Riviera Valley of Southern Switzerland are used to investigate convective boundary layer structure in a steep valley and to evaluate wind and temperature fields, convective boundary layer height, and surface sensible heat fluxes as predicted by the mesoscale model RAMS. Current parameterizations of surface and boundary layer processes in RAMS, as well as in other mesoscale models, are based on scaling laws strictly valid only for flat topography and uniform land cover. Model evaluation is required to investigate whether this limits the applicability of RAMS in steep, inhomogeneous terrain. One clear-sky day with light synoptic winds is selected from the field study. Observed temperature structure across and along the valley is nearly homogeneous while wind structure is complex with a wind speed maximum on one side of the valley. Upvalley flows are not purely thermally driven and mechanical effects near the valley entrance also affect the wind structure. RAMS captured many of the observed boundary layer characteristics within the steep valley. The wind field, temperature structure, and convective boundary layer height in the valley are qualitatively simulated by RAMS, but the horizontal temperature structure across and along the valley is less homogeneous in the model than in the observations. The model reproduced the observed net radiation, except around sunset and sunrise when RAMS does not take into account the shadows cast by the surrounding topography. The observed sensible heat fluxes fall within the range of simulated values at grid points surrounding the measurement sites. Some of the scatter between observed and simulated turbulent sensible heat fluxes are due to sub-grid scale effects related to local topography.     

A coupled numerical model to investigate the air–sea interaction at the coastal upwelling area of Cabo Frio,Brazil

This work investigates the capability of an oceanic numerical model dynamic and thermodynamically coupled to a three-dimensional mesoscale atmospheric numerical model to simulate the basic features of the air–sea interaction in the coastal upwelling area of Cabo Frio (RJ, Brazil). The upwelling/downwelling regime is an important feature in the oceanic circulation of Cabo Frio and determines the sustainability of local ecosystems. This regime is predominantly driven by the atmospheric circulation and is well documented, being suitable to be used as test reference for atmospheric and oceanic coupled and uncoupled models. The oceanic boundary conditions, coastline shape and coupling effect have been tested. The uncoupled oceanic model forced by a NE (SW) wind field generates a realistic upwelling (downwelling) phenomenon regardless of the proximity of the lateral boundary and how realistic is the shape of the coastline. The atmospheric-oceanic coupled model generates an upwelling location and intensity similar to the uncoupled simulation, but the upwelling is gradually enhanced by the sea-breeze circulation. It also generates vertical profiles of mixing ratio that compare better to the observations than the uncoupled simulation and air potential temperature and wind vertical profiles that represent particular features of the atmospheric circulation at Cabo Frio.     

A Numerical and Field Investigation of Surface Heat Fluxes from Small Wind-Sheltered Waterbodies in Semi-Arid Western Australia

Shelterbelts are used for a variety of purposes in agricultural environments, primarily because of their ability to improve the downwind microclimate. Excessive evaporative losses from small, agricultural water supply reservoirs in semi-arid Western Australia motivated a combined numerical modelling and field investigation into the potential for using shelterbelts to reduce evaporation from these open waterbodies. A numerical model of the disturbed momentum and turbulence fields in the region modified by the wind-shelter was employed and accounted for the presence of a waterbody downwind. The model was coupled with conservation equations for heat and moisture and sensible and latent heat fluxes were estimated from the simulated momentum, temperature and humidity fields. The numerical simulations were tested against four days of field data from two experiments conducted in the agricultural districts of southwest Western Australia that measured boundary-layer evolution over a variety of small waterbodies protected by artifical and natural wind-shelters. The model provided good predictions of windspeed during neutral conditions, but inadequate specification of the upwind boundary during non-neutral stabilities resulted in the model failing to capture any sensitivity to atmospheric stability as seen in the field data. Despite this limitation, the temperature and humidity fields were adequately captured by the model, and evaporative mass flux predictions also agreed well with estimates taken from water-balance measurements. It is concluded that well-designed wind-shelters can reduce evaporation from open waterbodies by 20–30% as a result of reductions in the velocity scales responsible for removing moisture from the water surface. The model can be used to estimate the values of various shelterbelt design parameters (e.g., porosity, height) that could be applied in the field to provide optimum evaporation reductions.     

Weak-wind mesoscale meandering in the nocturnal boundary layer

This study examines the strength and statistical behavior of mesoscale motions on time scales up to 1 h using eight data sets over different surface types. The mesoscale motions include internal gravity waves, microfront-like structures, horizontal modes, and a complex variety of other signatures, perhaps resulting from superposition of different modes. With weak large-scale flow, the mesoscale motions lead to meandering of the wind direction, as found in previous studies. However, the meandering often takes the form of sudden wind shifts rather than oscillation of wind direction. The relative strength and impact of such mesoscale motions are examined in terms of the constancy of the wind vector, the within-record standard deviation of the wind direction and the ratio of a meso-velocity scale to the speed of the large-scale flow. The strength of the mesoscale flow varies by an order of magnitude between nights at a given site and varies systematically between sites. The statistics of the vertical structure of such motions are examined for two of the data sets, both with sonic anemometers at seven levels.     

Impact of anthropogenic heat emissions on meteorological parameters and air quality in Beijing using a high-resolution model simulation

• The Large scale Urban Consumption of energY model was updated and coupled with WRF. • Anthropogenic heat emissions altered the precipitation and its spatial distribution. • A reasonable AHE scheme could improve the performance of simulated PM_2.5. • AHE aggravated the O₃ pollution in urban areas. Anthropogenic heat emissions (AHE) play an important role in modulating the atmospheric thermodynamic and kinetic properties within the urban planetary boundary layer, particularly in densely populated megacities like Beijing. In this study, we estimate the AHE by using a Large-scale Urban Consumption of energY (LUCY) model and further couple LUCY with a high-resolution regional chemical transport model to evaluate the impact of AHE on atmospheric environment in Beijing. In areas with high AHE, the 2-m temperature (T₂) increased to varying degrees and showed distinct diurnal and seasonal variations with maxima in night and winter. The increase in 10-m wind speed (WS₁₀) and planetary boundary layer height (PBLH) exhibited slight diurnal variations but showed significant seasonal variations. Further, the systematic continuous precipitation increased by 2.1 mm due to the increase in PBLH and water vapor in upper air. In contrast, the precipitation in local thermal convective showers increased little because of the limited water vapor. Meanwhile, the PM_2.5 reduced in areas with high AHE because of the increase in WS₁₀ and PBLH and continued to reduce as the pollution levels increased. In contrast, in areas where prevailing wind direction was opposite to that of thermal circulation caused by AHE, the WS₁₀ reduced, leading to increased PM_2.5. The changes of PM_2.5 illustrated that a reasonable AHE scheme might be an effective means to improve the performance of PM_2.5 simulation. Besides, high AHE aggravated the O₃ pollution in urban areas due to the reduction in NO_x.     

An evaluation of the dissimilarity in heat and momentum transport through quadrant analysis for an unstable atmospheric surface layer flow

Environmental Fluid Mechanics - To elucidate the role of each fluid motion in the transport of momentum and heat fluxes in an unstable atmospheric surface layer (ASL) flow from single point...     

A mesoscale vortex in a small stratified lake

A mesoscale vortex structure in the small stratified Lake Stechlin has been revealed by field experiments with satellite-tracked quasi-lagrangian drifters. The vortex with a radius of about 200 m drifted at 300 m/day along the western bight of the lake with nearly constant rotation speed of 3 cpd. Analysis of kinematical properties of the vortex motion demonstrates solid body character of rotation. Extrapolation of the vortex drift trajectory over the period preceding the observations combined with data on local winds and seiche dynamics has allowed tracing the vortex fate from its generation point. The normal modes analysis of the internal seiching in the lake reveals the vortex generation mechanism to be the interaction of certain seiche modes with local bottom topography and suggests generation of the mesoscale vortices to be the a regular feature of the lake circulation. Analysis of vorticity suggests additional energy supply to rotational flow, possibly from inverse cascading of small-scale turbulent motions—a feature typical for quasi-2D turbulence. The vortices can play an important role in the energy transport from basin-scale motions to small-scale boundary mixing. They can also contribute significantly to the horizontal heterogeneity of phyto- and zooplankton distribution as well as to the transport of dissolved matter such as nutrients between littoral and profundal areas. The topographically generated traveling vortices represent an analog of the synoptic eddies in the Ocean and in the Atmosphere, whereas their role in the lake hydrodynamics is practically unknown.     

Numerical experiments with RAMS model in highly complex terrain

RAMS mesoscale model has been applied to simulate the atmospheric circulation in highly complex terrain, in the Hindu-Kush Karakorum and Himalaya regions. The goal of the work is to assess the sensitivity of the model to the grid spacing and related resolution, the smoothing of the orography, the microphysics parameterizations, in such heterogeneous and challenging conditions. As a follow up, the capability of the model in correctly capturing the atmospheric processes is tested. Two main cases, where some measured data were available, have been considered: a period during which a flood occurred and a period characterised by high-pollution episodes. RAMS model provided sensible results, the predictions reasonably reproduce the observed data at a regional scale, but the most local characteristics cannot be definitely described even at the typical resolution of 1 km order.     

Numerical simulation of sand dune erosion

Erosion of sand or other granular material is a subject of utmost importance in several fields of practical interest, including industrial processes or environmental issues. Resulting from intricate interaction between the incident flow field and localized body forces responsible for the granular material cohesion, erosion is a particularly complex phenomenon. The present work addresses this problem, proposing a numerical method to compute the time evolution of a sand dune subjected to aeolian erosion, along with the associated entrainment and deposition fluxes. Turbulent fluid flow is computed through a three-dimensional Navier-Stokes solver based on a generalized coordinate system. A Lagrangian approach is adopted for tracking the trajectories of particles entrained in the saltation regime, thus allowing prediction of the corresponding deposition locations. Different models for saltation fluxes are tested, along with several formulations for the creeping-to-saltation flux ratio, creeping threshold and creeping distance. Comparison with results from wind tunnel experiments is very encouraging, stressing the relative importance of creeping in the erosion process for the presently studied conditions.     

Numerical modeling of minor gas constituents and aerosols in the atmosphere

In this paper, the structure of meso- and regional-scale atmospheric models is briefly described. The models incorporate transport of multicomponent pollutant transport/transformation and aerosol formation through nucleation, condensation–evaporation, and coagulation mechanisms. Also, an optimization problem is considered aimed at minimizing the emission source magnitudes. The results of numerical experiments are presented for a number of applications to demonstrate the model capabilities in solving environmental problems for different areas. Analysis of the numerical results shows that the models developed can be successfully used for solving regional environmental problems.     

A regional forest ecosystem carbon budget model: impacts of forest age structure and landuse history

This study investigated the impacts of landuse history and forest age structure on regional carbon fluxes for the forests in the Pacific Northwest of the United States based on a two-stage modeling strategy. In the first stage, an individual-based forest ecosystem carbon flux model (IntCarb) at stand scale is developed. IntCarb combines components from the ZELIG and CENTURY models to simulate forest development and heterotrophic respiration, respectively. Stand scale carbon fluxes simulated by IntCarb strongly depend on stand age. A forest stand can be a carbon sink for up to 200 years old with a peak at 30–40 years old. Old-growth stands are carbon neutral to the atmosphere in the long term. For any particular year, an old-growth stand can be either a carbon sink or source. The interannual variation of Net Ecosystem Productivity (NEP) for an old-growth stand is primarily determined by heterotrophic respiration. Due to the high spatial variability of stand ages, forest age structure needs to be taken into account to improve estimation of carbon budgets of forest ecosystems over large areas. In the stand stage, a regional carbon budget model (RegCarb) is developed to estimate regional carbon fluxes over large areas based on forest age structure, adjusting for the nonrespiratory carbon losses (timber harvesting). Our initial estimate with RegCarb for the Pacific Northwest of the United States indicates that this region was a tremendous carbon source to the atmosphere from 1890 to 1990 due to extensive logging of old-growth forest. Projection for the role of forests in this region in global carbon cycle in the future strongly depend on the amount of timber to be harvested, i.e. how the age structure of forests in this region is to be altered.     

A Numerical Study of the Dynamics of the Riverine Thermal Bar in a Deep Lake

A numerical model based on a Finite Volume formulation of the Navier–Stokes equations is used to simulate a range of scenarios leading to a thermal bar formed by a river inflow to an idealised deep lake. The results presented here show that small riverine salinity increases have a profound effect on the dynamics of the thermal bar, suppressing horizontal propagation of the plume and raising the possibility of a thermal bar which is capable of sinking to great depths. This finding is particularly relevant to Lake Baikal in Siberia, where the vigorous deep-water renewal is still not fully understood. An analysis of the buoyancy forces governing the depth of penetration of the thermal bar plume shows that realistic salinity gradients are an important factor in determining the circulation of Baikal waters. Observations of the saline curtailment of the thermal bar's horizontal propagation also reveal a potential for reduced productivity in the ecosystem of any temperate river delta during the Spring renewal period.     

Shallow flows: 2D or not 2D?

It is commonly assumed that shallow flows are in good approximation two-dimensional (2D) or quasi-2D. We will provide evidence that this is not always the case, and that the simple scaling argument based on the continuity equation does not always apply. Laboratory experiments on vortex flows in shallow fluid layers have revealed that locally significant three-dimensional (3D) effects and substantial vertical motions may occur, clearly destroying the assumed 2D character of the flow. For example, in the case of a dipolar vortex structure, an oscillatory vertical motion is observed in the vortex cores, while a spanwise circulation roll is present in front of the travelling dipole. These laboratory observations are confirmed by 3D numerical flow simulations. Attention will be given to a correct scaling analysis, in which both the aspect ratio of the fluid depth and a typical horizontal scale and the Reynolds number play a role.     

Subgrid-Scale Modeling of Reacting Scalar Fluxes in Large-Eddy Simulations of Atmospheric Boundary Layers

In large-eddy simulations of atmospheric boundary layer turbulence, the lumped coefficient in the eddy-diffusion subgrid-scale (SGS) model is known to depend on scale for the case of inert scalars. This scale dependence is predominant near the surface. In this paper, a scale-dependent dynamic SGS model for the turbulent transport of reacting scalars is implemented in large-eddy simulations of a neutral boundary layer. Since the model coefficient is computed dynamically from the dynamics of the resolved scales, the simulations are free from any parameter tuning. A set of chemical cases representative of various turbulent reacting flow regimes is examined. The reactants are involved in a first-order reaction and are injected in the atmospheric boundary layer with a constant and uniform surface flux. Emphasis is placed on studying the combined effects of resolution and chemical regime on the performance of the SGS model. Simulations with the scale-dependent dynamic model yield the expected trends of the coefficients as function of resolution, position in the flow and chemical regime, leading to resolution-independent turbulent reactant fluxes.     

Clustering of aerosols in atmospheric turbulent flow

A mechanism of formation of small-scale inhomogeneities in spatial distributions of aerosols and droplets associated with clustering instability in the atmospheric turbulent flow is discussed. The particle clustering is a consequence of a spontaneous breakdown of their homogeneous space distribution due to the clustering instability, and is caused by a combined effect of the particle inertia and a finite correlation time of the turbulent velocity field. In this paper a theoretical approach proposed in Elperin et al. (2002) Phys Rev E 66:036302 is further developed and applied to investigate the mechanisms of formation of small-scale aerosol inhomogeneities in the atmospheric turbulent flow. The theory of the particle clustering instability is extended to the case when the particle Stokes time is larger than the Kolmogorov time scale, but is much smaller than the correlation time at the integral scale of turbulence. We determined the criterion of the clustering instability for the Stokes number larger than 1. We discussed applications of the analyzed effects to the dynamics of aerosols and droplets in the atmospheric turbulent flow.     

High Spatial and Temporal Variability of Dry Deposition in a Coastal Region

A real meteorological situation characterized by strong spatial and temporal variability of the meteorological fields in a coastal region of eastern Denmark is examined in view of the transport of a passive tracer and dry deposition. Model simulations using a full mesoscale NWP model (COAMPS^TM) at different horizontal resolutions are performed. A realistic simulation showed that the differences in the amount of dry deposited matter can reach one order of magnitude and larger, during the period of one afternoon, depending on the model horizontal resolution. In addition, the results of an idealized experiment with straight coastline indicate that the horizontal model resolution alone is responsible for most of the differences. This study confirms the importance of high spatial and temporal resolution modelling for environmental applications.     

疏浚对大连湾生物可利用磷的水体浓度和沉积物-水界面交换通量的影响

利用FeO/CAM复合膜对大连湾已疏浚区与疏浚影响区海水水体中生物可利用磷(Bio-available phosphorus,BAP)浓度和水质指标进行调查,并通过实验室培养法、间隙水浓度梯度估算法和薄膜扩散梯度技术(Diffusivegradients in thin films technique,DGT)法分别测定沉积物-水界面间BAP交换通量.结果表明,由于受到调查海区附近正在进行的疏浚工程的影响,疏浚影响区水体中浊度和悬浮物含量平均值高于已疏浚区,且疏浚影响区水体中BAP的平均浓度(22.2μg L-1)高于已疏浚水体(19.5μg L-1);而疏浚影响区叶绿素a和溶解氧的浓度明显低于已疏浚区.沉积物-水界面间BAP交换通量研究结果显示,BAP主要由沉积物向上覆水中扩散,沉积物是BAP交换的源;间隙水浓度梯度估算法和DGT法测定的BAP交换通量相似,但远低于实验室培养法测定值;BAP在已疏浚区沉积物-水界面间的交换通量小于疏浚影响区,说明在一定时间内疏浚可以减弱沉积物内源磷的释放.     

The energy balance closure problem: an overview

This paper gives an overview of 20 years of research on the energy balance closure problem. It will be shown that former assumptions that measuring errors or storage terms are the reason for the unclosed energy balance do not stand up because even turbulent fluxes derived from documented methods and calibrated sensors, net radiation, and ground heat fluxes cannot close the energy balance. Instead, exchange processes on larger scales of the heterogeneous landscape have a significant influence. By including these fluxes, the energy balance can be approximately closed. Therefore, the problem is a scale problem and has important consequences to the measurement and modeling of turbulent fluxes.     

Daytime air pollution transport mechanisms in stable atmospheres of narrow versus wide urban valleys

Transport of air pollutants emitted from urban valleys can be strongly restricted by interactions between static and dynamic factors including topographic forcing, low-level atmospheric stability related to temperature inversions, and urban heat island-induced circulations. Interplay between these processes has a complex and dynamic nature, and is determinant for the evolution of different ventilation mechanisms and the associated impacts on air quality. Here we investigate these transport mechanisms through large eddy simulations using EULAG, an established model for multiscale flows, to simulate an idealized atmospheric environment in narrow versus wide urban valleys during critical conditions for air quality (high atmospheric stability). Our results show how the ventilation of valleys depends on a dynamic (variable during the daytime) balance between interacting and sometimes competing processes related to thermally-driven slope flows, urban heat island-induced flows, and the trapping effect of atmospheric stability; and how valley width affects this balance. Particularly important is that the time-space distribution of pollutants (a passive tracer) varies greatly between both valleys despite having the same urban area and emission rates. These variations lead to pollutants being mostly concentrated in different areas of the narrow and wide valleys. We discuss the mechanisms behind these results and their potential implications for real urban valleys. Further understanding of these mechanisms is crucial for explaining the occurrence of severe air pollution episodes and informing related decision-making processes in urban valleys.