Observed and modelled “chemical weather” during ESCOMPTE

The new MOdèle de Chimie Atmosphérique à Grande Echelle (MOCAGE) three-dimensional multiscale chemistry and transport model (CTM) has been applied to study heavy pollution episodes observed during the ESCOMPTE experiment. The model considers the troposphere and lower stratosphere, and allows the possibility of zooming from the planetary scale down to the regional scale over limited area subdomains. Like this, it generates its own time-dependent chemical boundary conditions in the vertical and in the horizontal. This paper focuses on the evaluation and quantification of uncertainties related to chemical and transport modelling during two intensive observing periods, IOP2 and IOP4 (June 20–26 and July 10–14, 2001, respectively). Simulations are compared to the database of four-dimensional observations, which includes ground-based sites and aircraft measurements, radiosoundings, and quasi-continuous measurements of ozone by LIDARs. Thereby, the observed and modelled day-to-day variabilities in air composition both at the surface and in the vertical have been assessed. Then, three sensitivity studies are conducted concerning boundary conditions, accuracy of the emission dataset, and representation of chemistry. Firstly, to go further in the analysis of chemical boundary conditions, results from the standard grid nesting set-up and altered configurations, relying on climatologies, are compared. Along with other recent studies, this work advocates the systematic coupling of limited-area models with global CTMs, even for regional air quality studies or forecasts. Next, we evaluate the benefits of using the detailed high-resolution emissions inventory of ESCOMPTE: improvements are noticeable both on ozone reactivity and on the concentrations of various species of the ozone photochemical cycle especially primary ones. Finally, we provide some insights on the comparison of two simulations differing only by the parameterisation of chemistry and using two state-of-the-art chemical schemes for regional photochemical modelling. Regional air quality modelling is found to be highly sensitive to the emission inventory dataset and also to the vertical and horizontal boundary conditions and detailed representation of chemistry. Interestingly enough, they infer the same range of errors compared to total model errors.     

The Kuroshio Current System as a jet and twin “relative” recirculation gyres embedded in the Sverdrup circulation

By analyzing the results of a realistic ocean general circulation model (OGCM) and conducting a series of idealized OGCM experiments, the dynamics of the Kuroshio Current System is examined. In the realistic configuration, the Kuroshio Current System is successfully simulated when the horizontal resolution of OGCMs is increased from 1/2° to 1/10°. The difference between the two experiments shows a jet, the model’s Kuroshio Extension, and a pair of cyclonic and anticyclonic, “relative,” recirculation gyres (RRGs) on the northern and southern flanks of the jet. We call them recirculation gyres because they share some features with ordinary recirculation gyres in previous studies, and we add the adjective “relative” to emphasize that they may not be apparent in the total field. Similar zonal jet and RRGs are obtained also in the idealized model with a rectangular basin and a flat bottom with a horizontal resolution of 1/6°. The northern RRG is generated by the injection of high potential vorticity (PV) created in the viscous sublayer of the western boundary current, indicating the importance of a no-slip boundary condition. Since there is no streamline with such high PV in the Sverdrup interior, the eastward current in the northern RRG region has to lose its PV anomaly by viscosity before connecting to the interior. In the setup stage this injection of high PV is carried out by many eddies generated from the instability of the western boundary current. This high PV generates the northern RRG, which induces the separation of the western boundary current and the formation of the zonal jet. In the equilibrium state, the anomalous high PV values created in the viscous sublayer are carried eastward in the northern flank of the zonal jet. The southern RRG is due to the classical Rhines–Young mechanism, where low PV values are advected northward within the western boundary inertial sublayer, and closed, PV-conserving streamlines form to the south of the Kuroshio Extension, allowing slow homogenization of the low PV anomalies. The westward-flowing southern branch of this southern RRG stabilizes the inertial western boundary current and prevents its separation in the northern half of the Sverdrup subtropical gyre, where the western boundary current is unstable without the stabilizing effect of the southern RRG. Therefore, in the equilibrium state, the southern RRG should be located just to the north of the center of the Sverdrup subtropical gyre, which is defined as the latitude of the Sverdrup streamfunction maximum. The zonal jet (the Kuroshio Extension) and the northern RRG gyre are formed to the north of the southern RRG. This is our central result. This hypothesis is confirmed by a series of sensitivity experiments where the location of the center of the Sverdrup subtropical gyre is changed without changing the boundaries of the subtropical gyre. The locations of the zonal jets in the observed Kuroshio Current System and Gulf Stream are consistent as well. Sensitivities of the model Kuroshio Current System are also discussed with regard to the horizontal viscosity, strength of the wind stress, and coastline.     

The stability of buoyancy-driven coastal currents

Buoyancy-driven boundary currents were generated in the laboratory by releasing buoyant fluid from a source adjacent to a vertical boundary in a rotating container. The boundary removed the Coriolis force parallel to it, allowing the buoyant fluid to spread in a current along the boundary. Ise of a cylindrical boundary and a line source that released fluid uniformly around the circumference enabled an axisymmetric (zonal) current to be produced. With the continuous release of fluid from the source, the current grew in width and depth until it became unstable to non-axisymmetric disturbances. The wavelength and phase velocities of the disturbances were consistent with a model of baroclinic instability of two-layer flow when frictional dissipation due to Ekman layers is included. However, when the current only occupied a small fraction of the total depth, barotropic processes were also thought to be important, with the growing waves gaining energy from the horizontal shear.In other experiments, gravity currents were produced by a point source adjacent to either a zonal (circular) or a meridional (radial) vertical boundary. The currents were also observed to become unstable to the same upstream breaking waves as those on the continuous zonal current. Finally, some comparisons are made with oceanic coastal currents.     

Overturning and wind-driven circulation in a low-order ocean–atmosphere model

A low-order ocean–atmosphere model is presented which combines coupling through heat exchange at the interface and wind stress forcing. The coupling terms are derived from the boundary conditions and the forcing terms of the constituents. Both the ocean and the atmosphere model are based on Galerkin truncations of the basic fluid dynamical equations. Hence, the coupled model can readily be extended to include more physics and more detail. The model presented here is the simplest of a hierarchy of low-order ocean–atmosphere models. The behaviour of the coupled model is investigated by means of geometric singular perturbation theory and bifurcation analysis. Two ways are found in which the slow time scales can play a role in the coupled dynamics. In the first scenario, a limit cycle on the overturning time scale is created. The associated oscillatory behaviour is governed by internal ocean dynamics. In the second scenario, intermittent behaviour occurs between periodic and chaotic regimes in parameter space.     

A simulation of heat flow in the lower troposphere by a laboratory model

A laboratory model of the troposphere has been constructed in which the heat flow above the ground in the lowest kilometer may be simulated. The scaling criteria, derived from the governing equations with a number of approximations, are 104 for length, 200 for velocity and 0.2 for temperature. Water is the simulating fluid, which moves through a closed-circuit water-tunnel. A set of electrically heated aluminum plates simulate the heating of the ground by insolation. Temperature and velocity fine structure are measured. The fluid isotherms are inferred from shadowgraph photography. All observations are confined to conditions of natural and free convection.A comparison between the temperature and velocity profiles within the prototype and the model, for identical boundary conditions, indicates that the model is a good representation of the prototype to heights of 1000 meters.It is found that the vertical transfer of warm air is along cylindrical columns, which take on different forms in natural and in free convection. In the former, the columns carry heat aloft into caps which rise at some 1.2 meter per second, or which may merge with neighbouring columns to continue upwards to heights of several hundred meters. If an elevated temperature-inversion layer is present, these convective columns perturb it, initiating horizontally-travelling waves. In free convection, the columns are inclined downwind, with their upper extremities extending almost horizontally; the cap is absent and the wave-like tail may separate and move across the fluid independently of its parent.     

Initializing a Mesoscale Boundary-Layer Model with Radiosonde Observations

A mesoscale boundary-layer model is used to simulate low-level regional wind fields over the La Plata River of South America, a region characterized by a strong daily cycle of land–river surface-temperature contrast and low-level circulations of sea–land breeze type. The initial and boundary conditions are defined from a limited number of local observations and the upper boundary condition is taken from the only radiosonde observations available in the region. The study considers 14 different upper boundary conditions defined from the radiosonde data at standard levels, significant levels, level of the inversion base and interpolated levels at fixed heights, all of them within the first 1500 m. The period of analysis is 1994–2008 during which eight daily observations from 13 weather stations of the region are used to validate the 24-h surface-wind forecast. The model errors are defined as the root-mean-square of relative error in wind-direction frequency distribution and mean wind speed per wind sector. Wind-direction errors are greater than wind-speed errors and show significant dispersion among the different upper boundary conditions, not present in wind speed, revealing a sensitivity to the initialization method. The wind-direction errors show a well-defined daily cycle, not evident in wind speed, with the minimum at noon and the maximum at dusk, but no systematic deterioration with time. The errors grow with the height of the upper boundary condition level, in particular wind direction, and double the errors obtained when the upper boundary condition is defined from the lower levels. The conclusion is that defining the model upper boundary condition from radiosonde data closer to the ground minimizes the low-level wind-field errors throughout the region.     

A multi-limit formulation for the equilibrium depth of a stably stratified boundary layer

Currently no expression for the equilibrium depth of the turbulent stably-stratified boundary layer is available that accounts for the combined effects of rotation, surface buoyancy flux and static stability in the free flow. Various expressions proposed to date are reviewed in the light of what is meant by the stable boundary layer. Two major definitions are thoroughly discussed. The first emphasises turbulence and specifies the boundary layer as a continuously and vigorously turbulent layer adjacent to the surface. The second specifies the boundary layer in terms of the mean velocity profile, e.g. by the proximity of the actual velocity to the geostrophic velocity. It is shown that the expressions based on the second definition are relevant to the Ekman layer and portray the depth of the turbulence in the intermediate regimes, when the effects of static stability and rotation essentially interfere. Limiting asymptotic regimes dominated by either stratification or rotation are examined using the energy considerations. As a result, a simple equation for the depth of the equilibrium stable boundary layer is developed. It is valid throughout the range of stability conditions and remains in force in the limits of a perfectly neutral layer subjected to rotation and a rotation-free boundary layer dominated by surface buoyancy flux or stable density stratification at its outer edge. Dimensionless coefficients are estimated using data from observations and large-eddy simulations. Well-known and widely used formulae proposed earlier by Zilitinkevich and by Pollard, Rhines and Thompson are shown to be characteristic of the above interference regimes, when the effects of rotation and static stability (due to either surface buoyancy flux, or stratification at the outer edge of the boundary layer) are roughly equally important.     

Understanding and simulating the link between African easterly waves and Atlantic tropical cyclones using a regional climate model: the role of domain size and lateral boundary conditions

Using a suite of lateral boundary conditions, we investigate the impact of domain size and boundary conditions on the Atlantic tropical cyclone and african easterly Wave activity simulated by a regional climate model. Irrespective of boundary conditions, simulations closest to observed climatology are obtained using a domain covering both the entire tropical Atlantic and northern African region. There is a clear degradation when the high-resolution model domain is diminished to cover only part of the African continent or only the tropical Atlantic. This is found to be the result of biases in the boundary data, which for the smaller domains, have a large impact on TC activity. In this series of simulations, the large-scale Atlantic atmospheric environment appears to be the primary control on simulated TC activity. Weaker wave activity is usually accompanied by a shift in cyclogenesis location, from the MDR to the subtropics. All ERA40-driven integrations manage to capture the observed interannual variability and to reproduce most of the upward trend in tropical cyclone activity observed during that period. When driven by low-resolution global climate model (GCM) integrations, the regional climate model captures interannual variability (albeit with lower correlation coefficients) only if tropical cyclones form in sufficient numbers in the main development region. However, all GCM-driven integrations fail to capture the upward trend in Atlantic tropical cyclone activity. In most integrations, variations in Atlantic tropical cyclone activity appear uncorrelated with variations in African easterly wave activity.     

An innovative artificial fog production device improved in the European project “FOG”

Transport safety is a major goal in the European Union. Low visibility conditions, especially due to fog, increase the risk of major accidents (chain collision). Innovative products have been developed by the automotive industry, including equipment manufacturers, to increase the level of safety of car passengers and drivers. Testing of these products requires the simulation or artificial reproduction of low visibility (fog) conditions with good stability and reproducibility characteristics. We report on the results of the European Union funded “FOG” project to improve road transport safety through fog production in an experimental test chamber located at the Clermont-Ferrand laboratory for research on road safety and visibility. The project developed a prototype of a small-scale climatic chamber, an improved fog production spraying device, a laser-based visibility measurement device, a reduced scale transmissometer, and a combined indoor climate-fog production simulation software. The ability of the fog chamber to test for driver reaction was also investigated. Recent developments include a device able to produce stable visibility levels and homogeneous fog, representative of various types of natural water droplet distribution. The fog characteristics were determined and compared to natural fog. Results are presented for a selection of conditions including stabilized visibility levels for dense fog and two kinds of droplet distributions.     

Inertial Effects on the Vertical Transport of Suspended Particles in a Turbulent Boundary Layer

In many atmospheric flows, a dispersed phase is actively suspended by turbulence, whose competition with gravitational settling ultimately dictates its vertical distribution. Examples of dispersed phases include snow, sea-spray droplets, dust, or sand, where individual elements of much larger density than the surrounding air are carried by turbulent motions after emission from the surface. In cases where the particle is assumed to deviate from local fluid motions only by its gravitational settling (i.e., they are inertialess), traditional flux balances predict a power-law dependence of particle concentration with height. It is unclear, however, how particle inertia influences this relationship, and this question is the focus of this work. Direct numerical simulations are conducted of turbulent open-channel flow, laden with Lagrangian particles of specified inertia; in this way the study focuses on the turbulent transport which occurs in the lowest few meters of the planetary boundary layer, in regions critical for connecting emission fluxes to the fluxes felt by the full-scale boundary layer. Simulations over a wide range of particle Stokes number, while holding the dimensionless settling velocity constant, are performed to understand the role of particle inertia on vertical dispersion. It is found that particles deviate from their inertialess behaviour in ways that are not easily captured by traditional theory; concentrations are reduced with increasing Stokes number. Furthermore, a similarity-based eddy diffusivity for particle concentration fails as particles experience inertial acceleration, precluding a closed-form solution for particle concentration as in the case of inertialess particles. The primary consequence of this result is that typical flux parametrizations connecting surface emission models (e.g., saltation models or sea-spray generation functions) to elevated boundary conditions may overestimate particle concentrations due to the reduced vertical transport caused by inertia in between; likewise particle emission may be underestimated if inferred from concentration measurements aloft.     

旋转扰动流体中壁面湍流不稳定及其拟序结构研究I：实验

重点介绍和讨论了中性条件下旋转扰动流体中边界层强迫不稳定及其相关的一些问题,阐述了旋转体系中切变驱动边界层不稳定的动力学特征。这些不稳定状态的研究在大气物理学、流体动力学、海洋学等多个领域中引起科学家极大的兴趣,近年来在实验和理论研究中都得到了不断的发展。意大利都灵大学基础物理系地球科学实验组通过水槽旋转实验方法,不断改变水槽启动或结束时的旋转运动速度,以及底部壁面粗糙度等要素,所得到的实验结果与SDBL理论非常一致。     

Stationary Rossby wave propagation in a shear flow along a reflective boundary

Summary This study investigates the impact of lateral boundary conditions on the propagation and dispersion of locally excited Rossby waves in a zonally periodic, barotropic, quasigeostrophic channel model on the β-plane. We use basic flows with either a linear meridional shear or a jet-like profile. On the southern boundary of the channel we impose either a rigid wall or a radiation condition, whereas the northern sidewall is permeable for Rossby waves. We compare the numerical solutions found for a reflecting southern boundary in a weakly dissipative flow to the solutions obtained from a WKB-analysis for the corresponding unforced nondissipative situation. Furthermore, we compare the generalized Eliassen-Palm flux vectors to the ray paths of Rossby wave packets, obtained from WKB ray tracing. In particular, we focus our investigation on the two-dimensional structure of trapped modal waves and wavetrains in a simple linear numerical model. Summarizing our results, we find that along the reflective wall, trapped modal wave structures as well as reflected wavetrains occur with characteristics (e.g., wavenumbers, turning latitudes) similar to the ones computed using asymptotic methods. In a linear sheared flow wave packets are trapped for all zonal wave numbers in contrast to a jet-like mean flow which has a selective effect on the waves; i.e., a turning latitudes phenomenon between the coast and the flow maximum occurs for short waves, while long waves can propagate freely across the zonal mean flow. This comes out clearly when studying the stream lines of the Eliassen-Palm flux vectors of the numerical model simulations. Furthermore, due to the reflected wave activity, the dispersion of Rossby waves is influenced by the southern boundary condition not only in the vicinity of the border but also in regions away from the boundary. These results appear to be important on the one hand for the existence of trapped Rossby waves in large-scale oceanic shear flows along a zonally oriented coast. And, on the other hand for large-scale boundary waves in conceptional atmospheric channel models which can lead to unwanted resonance effects. Received July 18, 2000/Revised June 9, 2001     

Effects of Time-Dependent Inflow Perturbations on Turbulent Flow in a Street Canyon

Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.     

区域模式的嵌套边界条件研究

在对现今常用的几种嵌套边界条件进行对比研究的基础上，提出一种普遍边界条件，它可以在不同的假定转变成各种常用的嵌套侧边界条件。由此可清晰地分析对比各种常用嵌套边界条件作用时的优缺点以及应注意的问题。     

Formation and Characteristics of Coastal Internal Boundary Layers During Onshore Flows

The development and characteristics of coastal internal boundary layers were investigated in 28 tests. These were made at all seasons and in both gradient and sea-breeze flows but only during mid-day periods. Measurements of turbulence and temperature were taken from a light aircraft which flew traverses across Long Island at successive altitudes parallel to the wind direction. These were used to locate the boundary between modified and unmodified air as a function of height and distance from the coast. The same measurements plus tower measurements of wind, turbulence and temperature, pilot balloon soundings and measurements of land and water surface temperatures by a remote sensing IR thermometer were used to quantify the characteristics of the modified and unmodified air. The boundary layer slope was steep close to the land-water interface and became shallower with downwind distance. Growth of the boundary layer was initially slower with stable lapse rates upwind than with neutral or unstable conditions over the water. An equilibrium height was found in many tests except under conditions of free convection when the internal boundary layer merged into the mixed layer inland and with sea-breeze conditions. The equilibrium height depended on downwind conditions and was greater with low wind speeds and strong land surface heating than with stronger winds and small land-water temperature differences. Current theoretical models are not adequate to predict the height of the boundary layer at the altitudes and distances studied but reasonably good predictions were given by an empirical model developed earlier. Wind speed in the modified air averaged about 70% of that at the coast but turbulence levels were several times higher both near the surface and aloft. These findings have important implications for diffusion from coastal sites.     

Formation and characteristics of coastal internal boundary layers during onshore flows

The development and characteristics of coastal internal boundary layers were investigated in 28 tests. These were made at all seasons and in both gradient and sea-breeze flows but only during mid-day periods. Measurements of turbulence and temperature were taken from a light aircraft which flew traverses across Long Island at successive altitudes parallel to the wind direction. These were used to locate the boundary between modified and unmodified air as a function of height and distance from the coast. The same measurements plus tower measurements of wind, turbulence and temperature, pilot balloon soundings and measurements of land and water surface temperatures by a remote sensing IR thermometer were used to quantify the characteristics of the modified and unmodified air. The boundary layer slope was steep close to the land-water interface and became shallower with downwind distance. Growth of the boundary layer was initially slower with stable lapse rates upwind than with neutral or unstable conditions over the water. An equilibrium height was found in many tests except under conditions of free convection when the internal boundary layer merged into the mixed layer inland and with sea-breeze conditions. The equilibrium height depended on downwind conditions and was greater with low wind speeds and strong land surface heating than with stronger winds and small land-water temperature differences. Current theoretical models are not adequate to predict the height of the boundary layer at the altitudes and distances studied but reasonably good predictions were given by an empirical model developed earlier. Wind speed in the modified air averaged about 70% of that at the coast but turbulence levels were several times higher both near the surface and aloft. These findings have important implications for diffusion from coastal sites.     

Comparison of Large-Eddy Simulation Data with Spatially Averaged Measurements Obtained by Acoustic Tomography – Presuppositions and First Results

An attempt is made to compare results oflarge-eddy simulation (LES) in a convective boundarylayer using the model PALM with experimental data obtained from acoustic travel time tomography.This method provides two-dimensional data arrays, which are considered as more suitable forLES-validation than classical local orline-integrated measurements, because the tomographic data are area- or volume-averaged.For a quantitative comparison with experimental data in general, some prerequisites have to be considered: First of all, the initial and boundary conditions of the LES model have to be provided correctly by the experiment. Considering measurement errors, a sensitivity study was performed to investigate the influence of inaccurate initial and boundary conditions on the simulation results.This showed that for determining some boundary conditions, such as the surface temperature and the roughness length, high measurement accuracies are necessary, which are difficult to reach or which at least require considerable extra measurement efforts.The initial and boundary conditions provided by the Lindenberg experiment in 1999 turned out to be of insufficient accuracy to allow quantitative comparisons.However, a qualitative comparison was performed instead to investigate if the acoustic tomography method is a proper method for comparisons with LES models in general.It showed a good qualitative agreement with some quantitative differences. These differences can partly be explained by the sensitivity of the LES to initial and boundary conditions and by the limitations of the acoustic tomography.     

A solution for simultaneous turbulent heat and vapor transfer between a water surface and the atmosphere

Interaction between sensible heat and water vapor diffusion in the lower atmosphere leads to the necessity of solving two simultaneous turbulent diffusion equations. This solution is obtained by the construction of Green's function which when incorporated in the boundary conditions produces two integral equations. These are solved by transformation into two algebraic equations by means of the Laplace Transformation. The results show how for a simple steady-state case, sensible heat and water vapor transfer and also the water surface temperature depend on the meteorological conditions and the rate of change of energy content of the water body. Due to advection, the water surface temperature and the turbulent fluxes vary in the downwind direction. However, for practical calculations of the mean evaporation or heat transfer, the error introduced by the use of an average temperature is usually quite small and negligible.     

Treatment of LBCs in 2D simulation of convection over hills

A series of idealized model simulations are analyzed to determine the sensitivity of model results to different configurations of the lateral boundary conditions (LBCs) in simulating mesoscale shallow convection over hilly terrain, In the simulations with steady thermal forcing at the model surface, a radiation condition at both boundaries is the best choice under high wind conditions, and the best results are produced when both the normal velocities and the temperature are treated with the radiation scheme in which the phase speed is the same for different variables, When the background wind speed is reasonably small, the LBC configuration with either the radiation or the zero gradient condition at both boundaries tends to make the numerical solution unstable. The choice of a constant condition at the inflow boundary and a radiation outflow boundary condition is appropriate in most cases, In the simulations with diurnal thermal forcing at the model surface, different LBC schemes are combined together to reduce spurious signals induced by the outflow boundary, A specification inflow boundary condition, in which the velocity fields at the inflow boundary are provided using the time-dependent results of a simulation with periodic LBCs over a flat domain, is tested and the results indicate that the specification condition at the inflow boundary makes it possible to use a smaller model domain to obtain reasonable results. The model horizontal domain length should be greater than a critical length, which depends on the domain depth H and the angle between gravity wave phase lines and the vertical, An estimate of minimum domain length is given by [(H-zi)/πU]√N^2L2x-4π^2u^2, where N and U are the background stability and wind speed,respectively, Lx is the typical gravity wavelength scale, and zi is the convective boundary layer (CBL)depth.