The modeling of the atmospheric boundary layer during convective conditions has long been a major source of uncertainty in the numerical modeling of meteorological conditions and air quality. Much of ...the difficulty stems from the large range of turbulent scales that are effective in the convective boundary layer (CBL). Both small-scale turbulence that is subgrid in most mesoscale grid models and large-scale turbulence extending to the depth of the CBL are important for the vertical transport of atmospheric properties and chemical species. Eddy diffusion schemes assume that all of the turbulence is subgrid and therefore cannot realistically simulate convective conditions. Simple nonlocal closure PBL models, such as the Blackadar convective model that has been a mainstay PBL option in the fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model (MM5) for many years and the original asymmetric convective model (ACM), also an option in MM5, represent large-scale transport driven by convective plumes but neglect small-scale, subgrid turbulent mixing. A new version of the ACM (ACM2) has been developed that includes the nonlocal scheme of the original ACM combined with an eddy diffusion scheme. Thus, the ACM2 is able to represent both the supergrid- and subgrid-scale components of turbulent transport in the convective boundary layer. Testing the ACM2 in one-dimensional form and comparing it with large-eddy simulations and field data from the 1999 Cooperative Atmosphere–Surface Exchange Study demonstrates that the new scheme accurately simulates PBL heights, profiles of fluxes and mean quantities, and surface-level values. The ACM2 performs equally well for both meteorological parameters (e.g., potential temperature, moisture variables, and winds) and trace chemical concentrations, which is an advantage over eddy diffusion models that include a nonlocal term in the form of a gradient adjustment.
A new combined local and nonlocal closure atmospheric boundary layer model called the Asymmetric Convective Model, version 2, (ACM2) was described and tested in one-dimensional form and was compared ...with large-eddy simulations and field data in Part I. Herein, the incorporation of the ACM2 into the fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5) is described. Model simulations using the MM5 with the ACM2 are made for the summer of 2004 and evaluated through comparison with surface meteorological measurements, rawinsonde profile measurements, and observed planetary boundary layer (PBL) heights derived from radar wind profilers. Overall model performance is as good as or better than similar MM5 evaluation studies. The MM5 simulations with the ACM2 compare particularly well to PBL heights derived from radar wind profilers during the afternoon hours. The ACM2 is designed to simulate the vertical mixing of any modeled quantity realistically for both meteorological models and air quality models. The next step, to be described in a subsequent article, is to incorporate the ACM2 into the Community Multiscale Air Quality (CMAQ) model for testing and evaluation.
In January 2013, a severe regional haze occurred over the North China Plain. An online-coupled meteorology-chemistry model was employed to simulate the impacts of aerosol-meteorology interactions on ...fine particles (PM2.5) pollution during this haze episode. The response of PM2.5 to meteorology change constituted a feedback loop whereby planetary boundary layer (PBL) dynamics amplified the initial perturbation of PM2.5. High PM2.5 concentrations caused a decrease of surface solar radiation. The maximal decrease in daily average solar radiation reached 53% in Beijing, thereby leading to a more stable PBL. The peak PBL height in Beijing decreased from 690 m to 590 m when the aerosol extinction was considered. Enhanced PBL stability suppressed the dispersion of air pollutants, and resulted in higher PM2.5 concentrations. The maximal increase of PM2.5 concentrations reached 140 g m−3 in Beijing. During most PM2.5 episodes, primary and secondary particles increased simultaneously. These results imply that the aerosol-radiation interactions played an important role in the haze episode in January 2013.
The Pleim–Xiu land surface model, Pleim surface layer scheme, and Asymmetric Convective Model (version 2) are now options in version 3.0 of the Weather Research and Forecasting model (WRF) Advanced ...Research WRF (ARW) core. These physics parameterizations were developed for the fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model (MM5) and have been used extensively by the air quality modeling community, so there was a need based on several factors to extend these parameterizations to WRF. Simulations executed with the new WRF physics are compared with simulations produced with the MM5 and another WRF configuration with a focus on the replication of near-surface meteorological conditions and key planetary boundary layer features. The new physics in WRF is recommended for retrospective simulations, in particular, those used to drive air quality simulations. In the summer, the error of all variables analyzed was slightly lower across the domain in the WRF simulation that used the new physics than in the similar MM5 configuration. This simulation had an even lower error than the other more common WRF configuration. For the cold season case, the model simulation was not as accurate as the other simulations overall, but did well in terms of lower 2-m temperature error in the western part of the model domain (plains and Rocky Mountains) and most of the Northeast. Both MM5 and the other WRF configuration had lower errors across much of the southern and eastern United States in the winter. The 2-m water vapor mixing ratio and 10-m wind were generally well simulated by the new physics suite in WRF when contrasted with the other simulations and modeling studies. Simulated planetary boundary layer features were compared with both wind profiler and aircraft observations, and the new WRF physics results in a more precise wind and temperature structure not only in the stable boundary layer, but also within most of the convective boundary layer. These results suggest that the WRFperformance is now at or above the level of MM5. It is thus recommended to drive future air quality applications.
Concentrations of both fine particulate matter (PM2.5) and ozone (O3) in the United States (US) have decreased significantly since 1990, mainly because of air quality regulations. Exposure to these ...air pollutants is associated with premature death. Here we quantify the annual mortality burdens from PM2.5 and O3 in the US from 1990 to 2010, estimate trends and inter-annual variability, and evaluate the contributions to those trends from changes in pollutant concentrations, population, and baseline mortality rates. We use a fine-resolution (36 km) self-consistent 21-year simulation of air pollutant concentrations in the US from 1990 to 2010, a health impact function, and annual county-level population and baseline mortality rate estimates. From 1990 to 2010, the modeled population-weighted annual PM2.5 decreased by 39 %, and summertime (April to September) 1 h average daily maximum O3 decreased by 9 % from 1990 to 2010. The PM2.5-related mortality burden from ischemic heart disease, chronic obstructive pulmonary disease, lung cancer, and stroke steadily decreased by 54 % from 123 700 deaths year−1 (95 % confidence interval, 70 800–178 100) in 1990 to 58 600 deaths year−1 (24 900–98 500) in 2010. The PM2.5-related mortality burden would have decreased by only 24 % from 1990 to 2010 if the PM2.5 concentrations had stayed at the 1990 level, due to decreases in baseline mortality rates for major diseases affected by PM2.5. The mortality burden associated with O3 from chronic respiratory disease increased by 13 % from 10 900 deaths year−1 (3700–17 500) in 1990 to 12 300 deaths year−1 (4100–19 800) in 2010, mainly caused by increases in the baseline mortality rates and population, despite decreases in O3 concentration. The O3-related mortality burden would have increased by 55 % from 1990 to 2010 if the O3 concentrations had stayed at the 1990 level. The detrended annual O3 mortality burden has larger inter-annual variability (coefficient of variation of 12 %) than the PM2.5-related burden (4 %), mainly from the inter-annual variation of O3 concentration. We conclude that air quality improvements have significantly decreased the mortality burden, avoiding roughly 35 800 (38 %) PM2.5-related deaths and 4600 (27 %) O3-related deaths in 2010, compared to the case if air quality had stayed at 1990 levels (at 2010 baseline mortality rates and population).
Air quality across the northern hemisphere over the past two decades has witnessed dramatic changes, with continuous improvement in developed countries in North America and Europe, but a contrasting ...sharp deterioration in developing regions of Asia.
This study investigates the historical trend in the long-term exposure to PM
and PM
-related premature mortality (PM
-mortality) and its response to changes in emission that occurred during 1990-2010 across the northern hemisphere. Implications for future trends in human exposure to air pollution in both developed and developing regions of the world are discussed.
We employed the integrated exposure-response model developed by Health Effects Institute to estimate the PM
-mortality. The 1990-2010 annual average PM
concentrations were obtained from the simulations using the WRF-CMAQ model. Emission mitigation efficiencies of sulfur dioxide (SO
), nitrogen oxides (NO
), ammonia (NH
), and primary PM are estimated from the PM
-mortality responses to the emission variations.
Estimated PM
-mortalities in East Asia and South Asia increased by 21% and 85% respectively, from 866,000 and 578,000 in 1990, to 1,048,000 and 1,068,000 in 2010. PM
-mortalities in developed regions (i.e., Europe and high-income North America) decreased substantially by 67% and 58% respectively.
Over the past two decades, correlations between population and PM
have become weaker in Europe and North America due to air pollution controls but stronger in East Asia due to deteriorating air quality. Mitigation of primary PM appears to be the most efficient way for increasing health benefits (i.e., providing the largest mortality reduction per unit emissions). However, reductions in emissions of NH
are needed to maximize the effectiveness of NO
emission controls. Citation: Wang J, Xing J, Mathur R, Pleim JE, Wang S, Hogrefe C, Gan CM, Wong DC, Hao J. 2017. Historical trends in PM
-related premature mortality during 1990-2010 across the northern hemisphere. Environ Health Perspect 125:400-408; http://dx.doi.org/10.1289/EHP298.
This note describes a simple scheme for analytical estimation of the surface-layer similarity functions from state variables. What distinguishes this note from the many previous papers on this topic ...is that this method is specifically targeted for numerical models in which simplicity and economic execution are critical. In addition, it has been in use in a mesoscale meteorological model for several years. For stable conditions, a very simple scheme is presented that compares well to the iterative solution. The stable scheme includes a very stable regime in which the slope of the stability functions is reduced to permit significant fluxes to occur, which is particularly important for numerical models in which decoupling from the surface can be an important problem. For unstable conditions, simple schemes generalized for varying ratios of aerodynamic roughness to thermal roughness (z₀/z0h
) are less satisfactory. Therefore, a simple scheme has been empirically derived for a fixedz₀/z0h
ratio, which represents quasi-laminar sublayer resistance.
We describe and analyze the results of the third global energy and water cycle experiment atmospheric boundary layer Study intercomparison and evaluation study for single-column models. Each of the ...nineteen participating models was operated with its own physics package, including land-surface, radiation and turbulent mixing schemes, for a full diurnal cycle selected from the Cabauw observatory archive. By carefully prescribing the temporal evolution of the forcings on the vertical column, the models could be evaluated against observations. We focus on the gross features of the stable boundary layer (SBL), such as the onset of evening momentum decoupling, the 2-m minimum temperature, the evolution of the inertial oscillation and the morning transition. New process diagrams are introduced to interpret the variety of model results and the relative importance of processes in the SBL; the diagrams include the results of a number of sensitivity runs performed with one of the models. The models are characterized in terms of thermal coupling to the soil, longwave radiation and turbulent mixing. It is shown that differences in longwave radiation schemes among the models have only a small effect on the simulations; however, there are significant variations in downward radiation due to different boundary-layer profiles of temperature and humidity. The differences in modelled thermal coupling to the land surface are large and explain most of the variations in 2-m air temperature and longwave incoming radiation among models. Models with strong turbulent mixing overestimate the boundary-layer height, underestimate the wind speed at 200 m, and give a relatively large downward sensible heat flux. The result is that 2-m air temperature is relatively insensitive to turbulent mixing intensity. Evening transition times spread 1.5 h around the observed time of transition, with later transitions for models with coarse resolution. Time of onset in the morning transition spreads 2 h around the observed transition time. With this case, the morning transition appeared to be difficult to study, no relation could be found between the studied processes, and the variation in the time of the morning transition among the models.
The Pleim–Xiu land surface model (PX LSM) has been improved by the addition of a second indirect data assimilation scheme. The first, which was described previously, is a technique in which soil ...moisture is nudged according to the biases in 2-m air temperature and relative humidity between the model- and observationbased analyses. The new technique involves nudging the deep soil temperature in the soil temperature force-restore (FR) model according to model bias in 2-m air temperature only during nighttime. While the FR technique is computationally efficient and very accurate for the special conditions for which it was derived, it is very dependent on the deep soil temperature that drives the restoration term of the surface soil temperature equation. Thus, adjustment of the deep soil temperature to optimize the 2-m air temperature during the night, when surface forcing is minimal, provides significant advantages over other methods of deep soil moisture initialization. Simulations of the Weather Research and Forecasting Model (WRF) using the PX LSM with and without the new deep soil temperature nudging scheme demonstrate substantial benefits of the new scheme for reducing error and bias of the 2-m air temperature. The effects of the new nudging scheme are most pronounced in the winter (January 2006) during which the model’s cold bias is greatly reduced. Air temperature error and bias are also reduced in a summer simulation (August 2006) with the greatest benefits in less vegetated and more arid regions. Thus, the deeptemperaturenudging scheme complements the soilmoisturenudging scheme because it is most effective for conditions in which the soil moisture scheme is least effective, that is, when evapotranspiration is not important (winter and arid climates).
The fourth phase of the Air Quality Model Evaluation
International Initiative (AQMEII4) is conducting a diagnostic
intercomparison and evaluation of deposition simulated by regional-scale air
quality ...models over North America and Europe. In this study, we analyze
annual AQMEII4 simulations performed with the Community Multiscale Air
Quality Model (CMAQ) version 5.3.1 over North America. These simulations
were configured with both the M3Dry and Surface Tiled Aerosol and Gas
Exchange (STAGE) dry deposition schemes available in CMAQ. A comparison of
observed and modeled concentrations and wet deposition fluxes shows that the
AQMEII4 CMAQ simulations perform similarly to other contemporary
regional-scale modeling studies. During summer, M3Dry has higher ozone
(O3) deposition velocities (Vd) and lower mixing ratios than STAGE
for much of the eastern US, while the reverse is the case over eastern
Canada and along the US West Coast. In contrast, during winter STAGE has higher
O3 Vd and lower mixing ratios than M3Dry over most of the
southern half of the modeling domain, while the reverse is the case for much
of the northern US and southern Canada. Analysis of the diagnostic
variables defined for the AQMEII4 project, i.e., grid-scale and land-use-specific effective conductances and deposition fluxes for the major dry deposition pathways, reveals generally higher summertime stomatal and
wintertime cuticular grid-scale effective conductances for M3Dry and
generally higher soil grid-scale effective conductances (for both vegetated
and bare soil) for STAGE in both summer and winter. On a domain-wide basis,
the stomatal grid-scale effective conductances account for about half of the
total O3 Vd during daytime hours in summer for both schemes.
Employing land-use-specific diagnostics, results show that daytime Vd varies
by a factor of 2 between land use (LU) categories. Furthermore, M3Dry vs. STAGE
differences are most pronounced for the stomatal and vegetated soil pathway
for the forest LU categories, with M3Dry estimating larger effective
conductances for the stomatal pathway and STAGE estimating larger effective
conductances for the vegetated soil pathway for these LU categories. Annual
domain total O3 deposition fluxes differ only slightly between M3Dry
(74.4 Tg yr−1) and STAGE (76.2 Tg yr−1), but pathway-specific fluxes to
individual LU types can vary more substantially on both annual and seasonal
scales, which would affect estimates of O3 damage to sensitive
vegetation. A comparison of two simulations differing only in their LU
classification scheme shows that the differences in LU cause seasonal mean
O3 mixing ratio differences on the order of 1 ppb across large portions
of the domain, with the differences generally being largest during summer and in
areas characterized by the largest differences in the fractional coverages
of the forest, planted and cultivated, and grassland LU categories. These
differences are generally smaller than the M3Dry vs. STAGE differences
outside the summer season but have a similar magnitude during summer.
Results indicate that the deposition impacts of LU differences are caused
by differences in the fractional coverages and spatial distributions of
different LU categories and the characterization of these categories
through variables like surface roughness and vegetation fraction in lookup
tables used in the land surface model and deposition schemes. Overall, the
analyses and results presented in this study illustrate how the diagnostic
grid-scale and LU-specific dry deposition variables adopted for AQMEII4 can
provide insights into similarities and differences between the CMAQ M3Dry
and STAGE dry deposition schemes that affect simulated pollutant budgets and
ecosystem impacts from atmospheric pollution.