Planetary Boundary Layer Mixing Research Articles

Modeling diurnal variation of surface PM<sub>2.5</sub> concentrations over East China with WRF-Chem: impacts from boundary-layer mixing and anthropogenic emission

Abstract. Diurnal variation of surface PM2.5 concentration (diurnal PM2.5) could dramatically affect aerosol radiative and health impacts and can also well reflect the physical and chemical mechanisms of air pollution formation and evolution. So far, diurnal PM2.5 and its modeling capability over East China have not been investigated and therefore are examined in this study. Based on the observations, the normalized diurnal amplitude of surface PM2.5 concentrations averaged over East China is weakest (∼1.2) in winter and reaches ∼1.5 in other seasons. The diurnal PM2.5 shows the peak concentration during the night in spring and fall and during the daytime in summer. The simulated diurnal PM2.5 with WRF-Chem and its contributions from multiple physical and chemical processes are examined in the four seasons. The simulated diurnal PM2.5 with WRF-Chem is primarily controlled by planetary boundary layer (PBL) mixing and emission variations and is significantly overestimated against the observation during the night. This modeling bias is likely primarily due to the inefficient PBL mixing of primary PM2.5 during the night. The simulated diurnal PM2.5 is sensitive to the PBL schemes and vertical-layer configurations with WRF-Chem. Besides the PBL height, the PBL mixing coefficient is also found to be the critical factor determining the PBL mixing of pollutants in WRF-Chem. With reasonable PBL height, the increase in the lower limit of the PBL mixing coefficient during the night can significantly reduce the modeling biases in diurnal PM2.5 and also the mean concentrations, particularly in the major cities of East China. It can also reduce the modeling sensitivity to the PBL vertical-layer configurations. The diurnal variation and injection height of anthropogenic emissions also play roles in simulating diurnal PM2.5, but the impact is relatively smaller than that from the PBL mixing. This study underscores that more efforts are needed to improve the boundary mixing process of pollutants in models with observations of PBL structure and mixing fluxes in addition to PBL height, in order to simulate reasonably the diurnal PM2.5 over East China. The diurnal variation and injection height of anthropogenic emissions must also be included to simulate the diurnal PM2.5 over East China.

Influences of planetary boundary layer mixing parameterization on summertime surface ozone concentration and dry deposition over North China

We present a regional modeling study that analyzes how planetary boundary layer (PBL) and surface layer parameterizations influence surface ozone concentrations and dry deposition fluxes over the Beijing-Tianjin-Hebei region in summer (July 2015). We use the Weather Research and Forecasting Model coupled to Chemistry (WRF-Chem) to simulate surface ozone concentration and dry deposition, and examine three PBL schemes: the Yonsei University (YSU), Mellor–Yamada–Janjić (MYJ), and Asymmetric Convective Model version 2 (ACM2) schemes. The model sensitivity to surface layer schemes is also tested by coupling the ACM2 PBL scheme with either the revised MM5-similarity scheme or the Pleim-Xiu surface layer scheme. Key physical and chemical factors for ozone dry deposition parameterization are analyzed to explore the root causes of model discrepancies. We find that all simulations overestimate the summertime daily mean ozone concentrations over North China (42 ppbv in observations vs. 43–50 ppbv in model results for July 2015) that are caused by high biases in daytime ozone and partly compensated by low biases in nighttime ozone. The YSU scheme has the largest overestimate in daily mean ozone concentration, but best reproduces the ozone diurnal cycle. The ACM2 scheme shows the largest underestimates of surface ozone over North China during nighttime, which can be explained by its weakest vertical mixing leading to high NOx concentrations and strong ozone titration near surface. The choices of PBL and surface layer schemes lead to over 20% differences in ozone dry deposition fluxes due to differences in simulated surface ozone concentrations and dry deposition velocities. We find the differences in dry deposition velocity are mainly caused by differences in Monin-Obukhov length during nighttime and surface temperature during daytime. Our study emphasizes the needs to better understand these key PBL and surface factors for reducing the uncertainties in model simulation of surface ozone concentration and dry deposition.

The Implications for Radiative Cloud Forcing via the Link Between Shallow Convection and Planetary Boundary Layer Mixing

AbstractIn this study, the role of shallow convection and its links with planetary boundary layer (PBL) mixing is analyzed using a global atmospheric model. All of the simulations are conducted with the Korean Integrated Model (KIM), which is a next‐generation operational global model for Korea Meteorological Administration. The simulation results show that the direct effect of shallow convection enhances additional vertical mixing inside the shallow convective layer, whereas its indirect effect reduces vertical mixing inside the PBL by changing the PBL height, thus vertical mixing strength. It is also found that the impact of shallow convection on low‐level clouds is tightly related to the overlapping between PBL and shallowing convection. When the PBL height is located in the higher (lower) part of the shallow convective layer, the mixing across the top of the PBL becomes smaller (larger). Therefore, the diffusivity profile of KIM shallow convection scheme is modified according to the overlapping pattern of PBL mixing and shallow convection. The systematic bias of downward solar radiation at the surface of KIM is improved with the revised diffusivity profile of the shallow convection scheme proposed in this study.

Impact of Turbulent Mixing in the Stratocumulus-Topped Boundary Layer on Numerical Weather Prediction

The impact of enhanced turbulent mixing induced by radiative cooling at the top of the stratocumulus-topped boundary layer (STBL) on numerical weather prediction is examined. An additional term involving top-down turbulent mixing via in-cloud radiative cooling is applied to the Yonsei University (YSU) planetary boundary layer (PBL) parameterization scheme using a top-down diffusivity profile and cloud-top entrainment. The modified scheme is evaluated in an advection fog case over the Yellow Sea of Korea using the Weather Research and Forecasting (WRF) model and in global medium-range forecasts using the Global/Regional Integrated Model system (GRIMs). In the fog case simulation, consideration of the additional top-down mixing parameterization in the YSU PBL simulates less formation and more rapid dispersion of the fog. As a result, the modified scheme simulates a drier and warmer boundary layer and a moister and cooler layer above the PBL. The modified algorithm also improves surface temperature prediction over the Yellow Sea accompanying early dissipation of the fog. In the global medium-range forecast experiment, the modified scheme simulates overall enhanced PBL mixing over the STBL in the tropics and subtropical ocean, showing drier and warmer regions near the surface and moister and cooler regions above the PBL, resulting in prediction of reduced low level cloud amount and increased downward shortwave radiation at the surface. The modified scheme appears to improve systematic bias in temperature and humidity in the lower troposphere compared to the control simulation.

Diurnal, seasonal, and annual trends in tropospheric CO in Southwest London during 2000–2015: Wind sector analysis and comparisons with urban and remote sites

Ambient carbon monoxide (CO) and meteorological parameters measured at the Egham (EGH) semi-rural site in SW London during 2000–2015 have permitted wind sector analysis of diurnal and seasonal cycles, and interpretation of long-term trends. CO daily amplitudes are used as a proxy for anthropogenic emissions. At EGH, morning and evening peaks in CO arise from the dominant contribution of road transport sources. Smaller amplitudes are observed during weekends than weekdays due to lower combustion emissions, and for mornings compared to evenings due to the timing of the development and break-up of the nocturnal inversion layer or planetary boundary layer (PBL). A wavelet transform revealed that the dominant mode of CO variability is the annual cycle, with apparent winter maxima likely due to increased CO emissions from domestic heating with summer minima ascribed to enhanced dispersion and dilution during the annual maximum of PBL mixing heights.Over the last two decades, both mitigation measures to reduce CO emissions and also a major switch to diesel cars, have accompanied a change at EGH from the dominance of local diurnal sources to a site measuring close to Atlantic background levels in summer months. CO observed in the S and SW wind sectors has declined by 4.7 and 5.9 ppb yr−1 respectively. The EGH CO record shows the highest levels in the early 2000s, with levels in E and calm winds comparable to those recorded at background stations in Greater London. However, since 2012, levels in S-SW sector have become more comparable with Mace Head background except during rush-hour periods. Marked declines in CO are observed during 2000–2008 for the NE, E, SE (London) and calm wind sectors, with the smallest declines observed for the S, SW and W (background) sectors. For the majority of wind sectors, the decline in CO is less noticeable since 2008, with an apparent stabilisation for NE, E and SE after 2009. The EGH CO data record exhibits a similar but slower exponential decay, but from a much lower starting concentration, than do CO data recorded at selected monitoring sites in urban areas in SE England. CO/CO2 residuals determined using hourly averaged datain the diurnal cycle demonstrate a clear decline in CO from 2000 to 2015 during daily periods of increased vehicle traffic, which is consistent with a sustained reduction in CO emissions from the road transport sector.

Simulated precipitation diurnal cycles over East Asia using different CAPE-based convective closure schemes in WRF model

Closure assumption in convection parameterization is critical for reasonably modeling the precipitation diurnal variation in climate models. This study evaluates the precipitation diurnal cycles over East Asia during the summer of 2008 simulated with three convective available potential energy (CAPE) based closure assumptions, i.e. CAPE-relaxing (CR), quasi-equilibrium (QE), and free-troposphere QE (FTQE) and investigates the impacts of planetary boundary layer (PBL) mixing, advection, and radiation on the simulation by using the weather research and forecasting model. The sensitivity of precipitation diurnal cycle to PBL vertical resolution is also examined. Results show that the precipitation diurnal cycles simulated with different closures all exhibit large biases over land and the simulation with FTQE closure agrees best with observation. In the simulation with QE closure, the intensified PBL mixing after sunrise is responsible for the late-morning peak of convective precipitation, while in the simulation with FTQE closure, convective precipitation is mainly controlled by advection cooling. The relative contributions of different processes to precipitation formation are functions of rainfall intensity. In the simulation with CR closure, the dynamical equilibrium in the free troposphere still can be reached, implying the complex cause-effect relationship between atmospheric motion and convection. For simulations in which total CAPE is consumed for the closures, daytime precipitation decreases with increased PBL resolution because thinner model layer produces lower convection starting layer, leading to stronger downdraft cooling and CAPE consumption. The sensitivity of the diurnal peak time of precipitation to closure assumption can also be modulated by changes in PBL vertical resolution. The results of this study help us better understand the impacts of various processes on the precipitation diurnal cycle simulation.

Evaluation of simulated climatological diurnal temperature range in CMIP5 models from the perspective of planetary boundary layer turbulent mixing

This study examines the effects of modeled planetary boundary layer (PBL) mixing on the simulated temperature diurnal cycle climatology over land in 20 CMIP5 models with AMIP simulations. When compared with observations, the magnitude of diurnal temperature range (DTR) is systematically underestimated over almost all land areas due to a widespread warm bias of daily minimum temperature (Tmin) and mostly a cold bias of daily maximum temperature (Tmax). Analyses of the CMIP5 multi-model ensemble means suggest that the biases of the simulated PBL mixing could very likely contribute to the temperature biases. For the regions with the cold bias in Tmax, the daytime PBL mixing is generally underestimated. The consequent more dry air entrainment from the free atmosphere could help maintain the surface humidity gradient, and thus produce more surface evaporation and potentially lower the Tmax. The opposite situation holds true for the regions with the warm bias of Tmax. This mechanism could be particularly applicable to the regions with moderate and wet climate conditions where surface evaporation depends more on the surface humidity gradient, but less on the available soil moisture. For the widespread warm bias of Tmin, the widely-recognized overestimated PBL mixing at nighttime should play a dominant role by transferring more heat from the atmosphere to the near-surface to warm the Tmin. Further analyses using the high resolution CFMIP2 output also support the CMIP5 results about the connections of the biases between the simulated turbulent mixing and the temperature diurnal cycle. The large inter-model variations of the simulated temperature diurnal cycle primarily appear over the arid and semi-arid regions and boreal arctic regions where the model differences in the PBL turbulence mixing could make equally significant contributions to the inter-model variations of DTR, Tmax and Tmin compared to the model differences in surface radiative processes. These results highlight the importance and need for accurate descriptions of the PBL processes with respect to the turbulent mixing in order to improve the temperature diurnal cycle simulations in climate models.

Atmospheric transport and chemistry of trace gases in LMDz5B: evaluation and implications for inverse modelling

Abstract. Representation of atmospheric transport is a major source of error in the estimation of greenhouse gas sources and sinks by inverse modelling. Here we assess the impact on trace gas mole fractions of the new physical parameterizations recently implemented in the atmospheric global climate model LMDz to improve vertical diffusion, mesoscale mixing by thermal plumes in the planetary boundary layer (PBL), and deep convection in the troposphere. At the same time, the horizontal and vertical resolution of the model used in the inverse system has been increased. The aim of this paper is to evaluate the impact of these developments on the representation of trace gas transport and chemistry, and to anticipate the implications for inversions of greenhouse gas emissions using such an updated model. Comparison of a one-dimensional version of LMDz with large eddy simulations shows that the thermal scheme simulates shallow convective tracer transport in the PBL over land very efficiently, and much better than previous versions of the model. This result is confirmed in three-dimensional simulations, by a much improved reproduction of the radon-222 diurnal cycle. However, the enhanced dynamics of tracer concentrations induces a stronger sensitivity of the new LMDz configuration to external meteorological forcings. At larger scales, the inter-hemispheric exchange is slightly slower when using the new version of the model, bringing them closer to observations. The increase in the vertical resolution (from 19 to 39 layers) significantly improves the representation of stratosphere/troposphere exchange. Furthermore, changes in atmospheric thermodynamic variables, such as temperature, due to changes in the PBL mixing modify chemical reaction rates, which perturb chemical equilibriums of reactive trace gases. One implication of LMDz model developments for future inversions of greenhouse gas emissions is the ability of the updated system to assimilate a larger amount of high-frequency data sampled at high-variability stations. Others implications are discussed at the end of the paper.

A thermal plume model for the Martian convective boundary layer

The Martian planetary boundary layer (PBL) is a crucial component of the Martian climate system. Global climate models (GCMs) and mesoscale models (MMs) lack the resolution to predict PBL mixing which is therefore parameterized. Here we propose to adapt the “thermal plume” model, recently developed for Earth climate modeling, to Martian GCMs, MMs, and single‐column models. The aim of this physically based parameterization is to represent the effect of organized turbulent structures (updrafts and downdrafts) on the daytime PBL transport, as it is resolved in large‐eddy simulations (LESs). We find that the terrestrial thermal plume model needs to be modified to satisfyingly account for deep turbulent plumes found in the Martian convective PBL. Our Martian thermal plume model qualitatively and quantitatively reproduces the thermal structure of the daytime PBL on Mars: superadiabatic near‐surface layer, mixing layer, and overshoot region at PBL top. This model is coupled to surface layer parameterizations taking into account stability and turbulent gustiness to calculate surface‐atmosphere fluxes. Those new parameterizations for the surface and mixed layers are validated against near‐surface lander measurements. Using a thermal plume model moreover enables a first‐order estimation of key turbulent quantities (e.g., PBL height and convective plume velocity) in Martian GCMs and MMs without having to run costly LESs.

Effects of Diurnal Cycle on a Simulated Asian Summer Monsoon

Abstract This study investigates the effects of the diurnal cycle on monsoonal circulations over Asia in summer with a focus on precipitation. To this end, two sets of experiments are designed in a regional climate modeling framework forced by reanalysis data. The control experiment is a normal integration in which radiation is computed hourly, whereas the no-diurnal experiment is an experimental integration in which the daily averaged solar flux is computed once a day. Analysis of the results from the two experiments reveals that the diurnal cycle enhances the daily averaged sensible heat flux over land and the latent flux over oceans, which means that daytime net solar heating exceeds nighttime cooling in terms of the effects in surface climate and monsoonal circulations. Seasonal precipitation increased by about 3% over land and 11% over oceans. The surface hydroclimate over land is strongly influenced by the interaction between land and the atmosphere, and results in cooler surface temperatures except over the Tibetan Plateau. Over oceans, a robust increase in precipitation results from enhanced planetary boundary layer mixing. The diurnal cycle over the Tibetan Plateau region is found to decrease surface albedo by melting snow during the daytime, which contributes to the formation of the thermal low near the surface and the Tibetan high in the upper troposphere. The resultant monsoonal precipitation is modulated by an increase (decrease) in precipitation over northern (southern) India. This modulation results in the summer monsoon over East Asia being shifted northward.

Intercomparison of Planetary Boundary-Layer Parametrizations in the WRF Model for a Single Day from CASES-99

This study compares five planetary boundary-layer (PBL) parametrizations in the Weather Research and Forecasting (WRF) numerical model for a single day from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) field program. The five schemes include two first-order closure schemes—the Yonsei University (YSU) PBL and Asymmetric Convective Model version 2 (ACM2), and three turbulent kinetic energy (TKE) closure schemes—the Mellor–Yamada–Janjic (MYJ), quasi-normal scale elimination (QNSE), and Bougeault–Lacarrere (BouLac) PBL. The comparison results reveal that discrepancies among thermodynamic surface variables from different schemes are large at daytime, while the variables converge at nighttime with large deviations from those observed. On the other hand, wind components are more divergent at nighttime with significant biases. Regarding PBL structures, a non-local scheme with the entrainment flux proportional to the surface flux is favourable in unstable conditions. In stable conditions, the local TKE closure schemes show better performance. The sensitivity of simulated variables to surface-layer parametrizations is also investigated to assess relative contributions of the surface-layer parametrizations to typical features of each PBL scheme. In the surface layer, temperature and moisture are more strongly influenced by surface-layer formulations than by PBL mixing algorithms in both convective and stable regimes, while wind speed depends on vertical diffusion formulations in the convective regime. Regarding PBL structures, surface-layer formulations only contribute to near-surface variability and then PBL mean properties, whereas shapes of the profiles are determined by PBL mixing algorithms.

Impacts of boundary layer mixing on pollutant vertical profiles in the lower troposphere: Implications to satellite remote sensing

Mixing in the planetary boundary layer (PBL) affects vertical distributions of air tracers in the lower troposphere. An accurate representation of PBL mixing is critical for chemical-transport models (CTMs) for applications sensitive to simulations of the vertical profiles of tracers. The full mixing assumption in the widely used global CTM GEOS-Chem has recently been supplemented with a non-local PBL scheme. This study analyzes the impact of the non-local scheme on model representation of PBL mixing, consequences for simulations of vertical profiles of air tracers and surface air pollution, and implications for model applications to the interpretation of data retrieved from satellite remote sensing. The non-local scheme significantly improves simulations of the vertical distributions for NO 2 and O 3, as evaluated using aircraft measurements in summer 2004. It also reduces model biases over the U.S. by more than 10 ppb for surface ozone concentrations at night and by 2–5 ppb for peak ozone in the afternoon, as evaluated using ground observations. The application to inverse modeling of anthropogenic NO x emissions for East China using satellite retrievals of NO 2 from OMI and GOME-2 suggests that the full mixing assumption results in 3–14% differences in top–down emission budgets as compared to the non-local scheme. The top–down estimate combining the non-local scheme and the Lin et al. inverse modeling approach suggests a magnitude of 6.6 TgN yr −1 for emissions of NO x over East China in July 2008 and 8.0 TgN yr −1 for January 2009, with the magnitude and seasonality in good agreement with bottom–up estimates.

Seasonal CO2 rectifier effect and large‐scale extratropical atmospheric transport

In atmospheric transport models, the covariation of the atmospheric transport and annually neutral biospheric CO2 flux is usually evident as the annual zonal mean surface CO2 concentration gradient. Using the NIES transport model and CO2 flux from the Biome‐BGC model, the covariations of different transport mechanisms and CO2 flux were examined and quantified. Including the covariation of the total transport (processes included in the NIES model) and CO2 flux, the annual average pole to pole CO2 concentration gradient is 3.5 ppm and interhemispheric difference of the average extratropical surface concentration is 2.5 ppm. The conventional covariation mechanism of the seasonal variation of planetary boundary layer mixing height and CO2 flux accounts for approximately 45% of the CO2 concentration gradient. Another important contribution to the CO2 concentration gradient in this model is the covariation of the extratropical anomaly transport (mainly by cyclones and anticyclones) and the biospheric flux, which accounts for about 55%. This alternate physical mechanism is the association of stronger meridional (north–south) anomaly transport (under strong baroclinic instability condition) with higher CO2 concentration from soil respiration in the winter and weaker anomaly transport (weak baroclinic instability condition) with lower CO2 concentration from photosynthetic uptake in the summer. The net result of the meridional transport and flux covariation is a north–south annual zonal mean CO2 concentration gradient.

Global model simulation of summertime U.S. ozone diurnal cycle and its sensitivity to PBL mixing, spatial resolution, and emissions

Simulation of summertime U.S. surface ozone diurnal cycle is influenced by the model representation of planetary boundary layer (PBL) mixing, spatial resolution, and precursor emissions. These factors are investigated here for five major regions (Northeast, Midwest, Southeast, California, and Southwest) by using the Model for Ozone And Related chemical Tracers version 2.4 (MOZART-2.4), with important modifications, to conduct sensitivity experiments for summer 1999 with three PBL mixing schemes, two horizontal resolutions and two emissions datasets. Among these factors, the PBL mixing is dominant. The default non-local scheme well reproduces the observed ozone diurnal variation, where the timing for the afternoon maximum and the morning minimum is within 1 h of the observed; biases for the minimum are less than 5 ppb except over the Southeast; and the ozone maximum–minimum contrast (OMMC) is within 10 ppb of observations except for the overprediction by 18.9 ppb over the Northeast. In contrast, the local scheme significantly overestimates the OMMC by 10–34 ppb over all regions as ozone and precursors are trapped too close to the ground. On the other hand, the full-mixing assumption underestimates the OMMC by 0–25 ppb, except over the Northeast, as the nighttime ozone decline is greatly underpredicted. As compared to PBL mixing, the effects of horizontal resolutions and precursor emissions being used are smaller but non-negligible. Overall, with the non-local mixing scheme, relatively high horizontal resolution (∼1.1°) and updated emissions data, the modified MOZART is capable of simulating the main features of the observed ozone diurnal cycle.

Measurements of carbon dioxide on very tall towers: results of the NOAA/CMDL program

Measurements of carbon dioxide (CO2) mixing ratios have been carried out since 1992 on a 610-m tall communications tower in North Carolina and since 1994 on a 447-m tall tower in Wisconsin. The data provide insights into the influence of pollution (fossil fuel combustion), biological exchange, boundary layer dynamics, and advective transport on CO2 mixing ratios over the continents. In this paper, we provide an overview of the data, describe access to the data, and suggest ways in which these results could be used to improve simulations of the global carbon cycle. In particular, the data will be very useful to constrain model estimates of covariance between terrestrial ecosystem fluxes of CO2 and diurnal and seasonal variations of planetary boundary layer mixing.

Use of a land-surface-transfer scheme (LSX) in a global climate model: the response to doubling stomatal resistance

One response of vegetation to future increases in atmospheric CO 2 may be a widespread increase in stomatal resistance. Such a response would increase plant water usage efficiency while still allowing CO 2 assimilation at current rates. The associated reduction in transpiration rates has the potential of causing significant modifications in climate on regional and global scales. This paper describes the effects of a uniform doubling of the stomatal resistance parameterization in a global climate model (GENESIS). The model includes a land-surface transfer scheme (LSX) that accounts for the physical effects of vegetation, including stomatal resistance and transpiration, which is described in detail in an appendix. The atmospheric general circulation model is a heavily modified version of the NCAR Community Climate Model version 1 with new treatments of clouds, penetrative convection, planetary boundary layer mixing, solar radiation, the diurnal cycle, and semi-Lagrangian transport of water vapor. The other surface models include multi-layer models of soil, snow and sea ice, and a 50-m slab ocean mixed layer. The effects of doubling the stomatal resistance parameterization are largest in heavily forested regions: tropical South America, and parts of the Northern Hemispheric boreal forests in Canada, Russia and Siberia in summer. The primary surface changes are a decrease in evapotranspiration, an increase in upward sensible heat flux, and a surface-air warming. Secondary effects include shifts in the ITCZ which cause large increases in precipitation, soil moisture and runoff in western tropical South America, and decreases in these quantities in northern subtropical Africa. Noticeable changes in relative humidity, cloudiness and meridional circulation occur throughout the troposphere. The global effects on atmospheric temperature and specific humidity are small fractions of those found in other doubled CO 2 experiments. However, unlike doubled CO 2 the signs of those changes combine to give relatively large reductions in relative humidity and cloudiness. It is suggested that the stomatal-resistance effect and other plant responses to large-scale environmental perturbations should be included in models of future climate.