scholarly journals LES study of the impact of moist thermals on the oxidative capacity of the atmosphere in southern West Africa

2018 ◽  
Vol 18 (9) ◽  
pp. 6601-6624 ◽  
Author(s):  
Fabien Brosse ◽  
Maud Leriche ◽  
Céline Mari ◽  
Fleur Couvreux

Abstract. The hydroxyl radical (OH) is a highly reactive species and plays a key role in the oxidative capacity of the atmosphere. We explore the potential impact of a convective boundary layer on reconciling the calculation–measurement differences for OH reactivity (the inverse of OH lifetime) attributable to the segregation of OH and its reactants by thermals and the resulting modification of averaged reaction rates. The large-eddy simulation version of the Meso-NH model is used, coupled on-line with a detailed chemistry mechanism to simulate two contrasted biogenic and urban chemical regimes. In both environments, the top of the boundary layer is the region with the highest calculated segregation intensities but with the opposite sign. In the biogenic environment, the inhomogeneous mixing of isoprene and OH leads to a maximum decrease of 30 % of the mean reaction rate in this zone. In the anthropogenic case, the effective rate constant for OH reacting with aldehydes is 16 % higher than the averaged value. OH reactivity is always higher by 15 to 40 % inside thermals in comparison to their surroundings as a function of the chemical environment and time of the day. Since thermals occupy a small fraction of the simulated domain, the impact of turbulent motions on domain-averaged total OH reactivity reaches a maximum decrease of 9 % for the biogenic case and a maximum increase of 5 % for the anthropogenic case. Accounting for the segregation of air masses by turbulent motions in regional and global models may increase OH reactivity in urban environments but lower OH reactivity in biogenic environments. In both cases, segregation alone is insufficient for resolving the underestimation between observed and modeled OH reactivity.

2012 ◽  
Vol 12 (19) ◽  
pp. 9335-9353 ◽  
Author(s):  
H. G. Ouwersloot ◽  
J. Vilà-Guerau de Arellano ◽  
A. C. Nölscher ◽  
M. C. Krol ◽  
L. N. Ganzeveld ◽  
...  

Abstract. We studied the atmospheric boundary layer (ABL) dynamics and the impact on atmospheric chemistry during the HUMPPA-COPEC-2010 campaign. We used vertical profiles of potential temperature and specific moisture, obtained from 132 radio soundings, to determine the main boundary layer characteristics during the campaign. We propose a classification according to several main ABL prototypes. Further, we performed a case study of a single day, focusing on the convective boundary layer, to analyse the influence of the dynamics on the chemical evolution of the ABL. We used a mixed layer model, initialized and constrained by observations. In particular, we investigated the role of large scale atmospheric dynamics (subsidence and advection) on the ABL development and the evolution of chemical species concentrations. We find that, if the large scale forcings are taken into account, the ABL dynamics are represented satisfactorily. Subsequently, we studied the impact of mixing with a residual layer aloft during the morning transition on atmospheric chemistry. The time evolution of NOx and O3 concentrations, including morning peaks, can be explained and accurately simulated by incorporating the transition of the ABL dynamics from night to day. We demonstrate the importance of the ABL height evolution for the representation of atmospheric chemistry. Our findings underscore the need to couple the dynamics and chemistry at different spatial scales (from turbulence to mesoscale) in chemistry-transport models and in the interpretation of observational data.


2020 ◽  
Author(s):  
Jan Weinkaemmerer ◽  
Ivan Bašták Ďurán ◽  
Jürg Schmidli

<p>In the convective boundary layer over mountainous regions, the mean values and the fluxes of quantities like heat, mass, and momentum are strongly influenced by thermally induced flows. Several studies have pointed out that the enhanced warming of the air inside a valley can be explained by the valley-volume effect whereas the cross-valley circulation leads to a net export of heat to the free atmosphere. We are interested in the influence of an upper-level wind on the local circulations and the boundary-layer properties, both locally and in terms of the horizontal mean, as this aspect has not yet received much attention. LES are carried out over idealized, two-dimensional topographies using the CM1 numerical model. For the analysis, turbulent, mean-circulation, and large-scale contributions are systematically distinguished. Also, budget analyses are performed for the turbulence kinetic energy and the turbulent heat and mass flux. Based on the first results for periodic topographies, no crucial influence on the horizontally averaged heat-flux and temperature profile can be observed, even though the flow pattern of the thermal wind is qualitatively changed. In addition to that, the impact on moisture transport will be evaluated and simulations over different topographies as well as for different atmospheric conditions and surface properties will be presented.</p>


2015 ◽  
Vol 15 (20) ◽  
pp. 11981-11998 ◽  
Author(s):  
M. N. Lang ◽  
A. Gohm ◽  
J. S. Wagner

Abstract. Idealized large-eddy simulations were performed to investigate the impact of different mountain geometries on daytime pollution transport by thermally driven winds. The main objective was to determine interactions between plain-to-mountain and slope wind systems, and their influence on the pollution distribution over complex terrain. For this purpose, tracer analyses were conducted over a quasi-two-dimensional mountain range with embedded valleys bordered by ridges with different crest heights and a flat foreland in cross-mountain direction. The valley depth was varied systematically. It was found that different flow regimes develop dependent on the valley floor height. In the case of elevated valley floors, the plain-to-mountain wind descends into the potentially warmer valley and replaces the opposing upslope wind. This superimposed plain-to-mountain wind increases the pollution transport towards the main ridge by an additional 20 % compared to the regime with a deep valley. Due to mountain and advective venting, the vertical exchange is 3.6 times higher over complex terrain than over a flat plain. However, the calculated vertical exchange is strongly sensitive to the definition of the convective boundary layer height. In summary, the impact of the terrain geometry on the mechanisms of pollution transport confirms the necessity to account for topographic effects in future boundary layer parameterization schemes.


2014 ◽  
Vol 14 (8) ◽  
pp. 4135-4167 ◽  
Author(s):  
K. Toyota ◽  
A. P. Dastoor ◽  
A. Ryzhkov

Abstract. Atmospheric mercury depletion events (AMDEs) refer to a recurring depletion of mercury occurring in the springtime Arctic (and Antarctic) boundary layer, in general, concurrently with ozone depletion events (ODEs). To close some of the knowledge gaps in the physical and chemical mechanisms of AMDEs and ODEs, we have developed a one-dimensional model that simulates multiphase chemistry and transport of trace constituents throughout porous snowpack and in the overlying atmospheric boundary layer (ABL). This paper constitutes Part 2 of the study, describing the mercury component of the model and its application to the simulation of AMDEs. Building on model components reported in Part 1 ("In-snow bromine activation and its impact on ozone"), we have developed a chemical mechanism for the redox reactions of mercury in the gas and aqueous phases with temperature dependent reaction rates and equilibrium constants accounted for wherever possible. Thus the model allows us to study the chemical and physical processes taking place during ODEs and AMDEs within a single framework where two-way interactions between the snowpack and the atmosphere are simulated in a detailed, process-oriented manner. Model runs are conducted for meteorological and chemical conditions that represent the springtime Arctic ABL characterized by the presence of "haze" (sulfate aerosols) and the saline snowpack on sea ice. The oxidation of gaseous elemental mercury (GEM) is initiated via reaction with Br-atom to form HgBr, followed by competitions between its thermal decomposition and further reactions to give thermally stable Hg(II) products. To shed light on uncertain kinetics and mechanisms of this multi-step oxidation process, we have tested different combinations of their rate constants based on published laboratory and quantum mechanical studies. For some combinations of the rate constants, the model simulates roughly linear relationships between the gaseous mercury and ozone concentrations as observed during AMDEs/ODEs by including the reaction HgBr + BrO and assuming its rate constant to be the same as for the reaction HgBr + Br, while for other combinations the results are more realistic by neglecting the reaction HgBr + BrO. Speciation of gaseous oxidized mercury (GOM) changes significantly depending on whether or not BrO is assumed to react with HgBr to form Hg(OBr)Br. Similarly to ozone (reported in Part 1), GEM is depleted via bromine radical chemistry more vigorously in the snowpack interstitial air than in the ambient air. However, the impact of such in-snow sink of GEM is found to be often masked by the re-emissions of GEM from the snow following the photo-reduction of Hg(II) deposited from the atmosphere. GOM formed in the ambient air is found to undergo fast "dry deposition" to the snowpack by being trapped on the snow grains in the top ~1 mm layer. We hypothesize that liquid-like layers on the surface of snow grains are connected to create a network throughout the snowpack, thereby facilitating the vertical diffusion of trace constituents trapped on the snow grains at much greater rates than one would expect inside solid ice crystals. Nonetheless, on the timescale of a week simulated in this study, the signal of atmospheric deposition does not extend notably below the top 1 cm of the snowpack. We propose and show that particulate-bound mercury (PBM) is produced mainly as HgBr42− by taking up GOM into bromide-enriched aerosols after ozone is significantly depleted in the air mass. In the Arctic, "haze" aerosols may thus retain PBM in ozone-depleted air masses, allowing the airborne transport of oxidized mercury from the area of its production farther than in the form of GOM. Temperature dependence of thermodynamic constants calculated in this study for Henry's law and aqueous-phase halide complex formation of Hg(II) species is a critical factor for this proposition, calling for experimental verification. The proposed mechanism may explain observed changes in the GOM–PBM partitioning with seasons, air temperature and the concurrent progress of ozone depletion in the high Arctic. The net deposition of mercury to the surface snow is shown to increase with the thickness of the turbulent ABL and to correspond well with the column amount of BrO in the atmosphere.


2017 ◽  
Vol 17 (3) ◽  
pp. 1623-1640 ◽  
Author(s):  
Metodija M. Shapkalijevski ◽  
Huug G. Ouwersloot ◽  
Arnold F. Moene ◽  
Jordi Vilà-Guerau de Arrellano

Abstract. By characterizing the dynamics of a convective boundary layer above a relatively sparse and uniform orchard canopy, we investigated the impact of the roughness-sublayer (RSL) representation on the predicted diurnal variability of surface fluxes and state variables. Our approach combined numerical experiments, using an atmospheric mixed-layer model including a land-surface-vegetation representation, and measurements from the Canopy Horizontal Array Turbulence Study (CHATS) field experiment near Dixon, California. The RSL is parameterized using an additional factor in the standard Monin–Obukhov similarity theory flux-profile relationships that takes into account the canopy influence on the atmospheric flow. We selected a representative case characterized by southerly wind conditions to ensure well-developed RSL over the orchard canopy. We then investigated the sensitivity of the diurnal variability of the boundary-layer dynamics to the changes in the RSL key scales, the canopy adjustment length scale, Lc, and the β = u*/|U| ratio at the top of the canopy due to their stability and dependence on canopy structure. We found that the inclusion of the RSL parameterization resulted in improved prediction of the diurnal evolution of the near-surface mean quantities (e.g. up to 50 % for the wind velocity) and transfer (drag) coefficients. We found relatively insignificant effects on the modelled surface fluxes (e.g. up to 5 % for the friction velocity, while 3 % for the sensible and latent heat), which is due to the compensating effect between the mean gradients and the drag coefficients, both of which are largely affected by the RSL parameterization. When varying Lc (from 10 to 20 m) and β (from 0.25 to 0.4 m), based on observational evidence, the predicted friction velocity is found to vary by up to 25 % and the modelled surface-energy fluxes (sensible heat, SH, and latent heat of evaporation, LE) vary up to 2 and 9 %. Consequently, the boundary-layer height varies up to 6 %. Furthermore, our analysis indicated that to interpret the CHATS measurements above the canopy, the contributions of non-local effects such as entrainment, subsidence and the advection of heat and moisture over the CHATS site need to be taken into account.


2018 ◽  
Vol 35 (8) ◽  
pp. 1687-1699 ◽  
Author(s):  
Adam L. Houston ◽  
Jason M. Keeler

AbstractThe objective of the research presented is to assess the impact of sensor response and aircraft airspeed on the accuracy of in situ observations collected by small unmanned aircraft systems profiling the convective boundary layer or transecting airmass boundaries. Estimates are made using simulated aircraft flown within large-eddy simulations. Both instantaneous errors (differences between observed temperature, which include the effects of sensor response and airspeed, and actual temperature) and errors in representation (differences between serial observations and representative snapshots of the atmospheric state) are considered. Synthetic data are retrieved assuming a well-aspirated first-order sensor mounted on rotary-wing aircraft operated as profilers in a simulated CBL and fixed-wing aircraft operated through transects across a simulated airmass boundary. Instantaneous errors are found to scale directly with sensor response time and airspeed for both CBL and airmass boundary experiments. Maximum errors tend to be larger for airmass boundary transects compared to the CBL profiles. Instantaneous errors for rotary-wing aircraft profiles in the CBL simulated for this work are attributable to the background lapse rate and not to turbulent temperature perturbations. For airmass boundary flights, representation accuracy is found to degrade with decreasing airspeed. This signal is most pronounced for flights that encounter the density current wake. When representation errors also include instantaneous errors resulting from sensor response, instantaneous errors are found to be dominant for flights that remain below the turbulent wake. However, for flights that encounter the wake, sensor response times generally need to exceed ~5 s before instantaneous errors become larger than errors in representation.


2019 ◽  
Vol 9 (2) ◽  
pp. 357 ◽  
Author(s):  
Mi-Sug Kim

The turbulence kinetics model (TKM) describes an overall reaction rate for microscopic mass transfer phenomenon expressed as separation intensity, Is, in a turbulent reacting flow. This study examines the effects of turbulent mixing in the convective boundary layer (CBL) on essential NOx–O3–Hydrocarbon photochemistry containing sources of NO and a surrogate reactive hydrocarbon. The modeling approach applies for all species except OH with an assumption of a photostationary steady state. The TKM results reveal principal findings as follows: (1) effects of turbulence on reaction rates lead to significant segregations throughout most of the CBL in reaction pairs NO + O3, RH + OH and NO + HO2; (2) segregations permit significantly higher concentrations of NO and RH to build up and endure in the CBL than would occur for a non-turbulent atmosphere; (3) turbulent segregation influences limit and shift the ranges of NO and O3 concentrations compared to the non-turbulent case; (4) while there are differences between the TKM results and those for a published Large Eddy simulation (LES) of the same chemical system, there are also strong similarities. Therefore, a future study remains to compare model results to observations if and when appropriately time-resolved measurements of reacting species are obtained.


2013 ◽  
Vol 6 (1) ◽  
pp. 721-790 ◽  
Author(s):  
S. Strada ◽  
S. R. Freitas ◽  
C. Mari ◽  
K. M. Longo ◽  
R. Paugam

Abstract. Wild-fires release huge amounts of aerosol and hazardous trace gases in the atmosphere. The residence time and the dispersion of fire pollutants in the atmosphere can range from hours to days and from local to continental scales. These various scenarios highly depend on the injection height of smoke plumes. The altitude at which fire products are injected in the atmosphere is controlled by fire characteristics and meteorological conditions. Injection height however is still poorly accounted in chemistry transport models for which fires are sub-grid scale processes which need to be parametrised. Only recently, physically-based approaches for estimating the fire injection heights have been developed which consider both the convective updrafts induced by the release of fire sensible heat and the impact of background meteorological environment on the fire convection dynamics. In this work, two different models are used to simulate fire injection heights in contrasted meteorological scenarios: a Mediterranean arson fire and two Amazonian deforestation fires. A Eddy-Diffusivity/Mass-Flux approach, formerly developed to reproduce convective boundary layer in the non-hydrostatic meteorological model Meso-NH, is compared to the 1-D Plume Rise Model. For both models, radiosonde data and re-analyses from the European Center for Medium-Range Weather Forecasts (ECMWF) have been used as initial conditions to explore the sensitivity of the models responses to different meteorological forcings. The two models predict injection heights for the Mediterranean fire between 1.7 and 3.3 km with the Meso-NH/EDMF model systematically higher than the 1-D PRM model. Both models show a limited sensitivity to the meteorological forcings with a 20–30% difference in the injection height between radiosondes and ECMWF data for this case. Injection heights calculated for the two Amazonian fires ranges from 5 to 6.5 km for the 1-D PRM model and from 2 to 4 km for the Meso-NH/EDMF model. The difference of smoke plume heights between the two models can reach 3–4 km. A large difference is obtained for the windy-wet Amazonian fire by the 1-D PRM model with a injection height 1.5 km higher when ECMWF re-analyses are used compared to the run with the radiosonde forcing. For the Mediterranean case, both models forecast a plume injection height above the boundary layer, although there are evidences that this particular fire propagated near the surface, highlighting the current limitations of the two approaches.


2008 ◽  
Vol 8 (2) ◽  
pp. 6793-6824
Author(s):  
F. Chosson ◽  
R. Paoli ◽  
B. Cuenot

Abstract. Detailed ship plume simulations in various convective boundary layer situations have been performed using a Lagrangian Dispersion Model driven by a Large Eddy Simulation Model. The simulations focus on early stage (1–2 h) of plume dispersion regime and take into account the effects of plume rise on dispersion. Results are presented in an attempt to provide to chemical modellers community a realistic description of the impact of characteristic dispersion on exhaust ship plume chemistry. Plume dispersion simulations are used to derive analytical dilution rate functions. Even though results exhibit striking effects of plume rise parameter on dispersion patterns, it is shown that initial buoyancy fluxes at ship stack have minor effect on plume dilution rate. After initial high dispersion regimes a simple characteristic dilution time scale can be used to parameterize the subgrid plume dilution effects in large scale chemistry models. The results show that this parameter is directly related to the typical turn-over time scale of the convective boundary layer.


Sign in / Sign up

Export Citation Format

Share Document