Dry Deposition Methods Based on Turbulence Kinetic Energy: Part 2. Extension to Particle Deposition Using a Single-Point Model.

Magnitude of atmospheric turbulence, a key driver of several processes that contribute to aerosol (i.e., particle) deposition, is underrepresented in current models. Various formulations have been developed to model particle dry deposition; all these formulations typically rely on friction velocity and some use additional ad hoc factors to represent enhanced impacts of turbulence. However, none were formally linked with the three-dimensional (3-D) turbulence. Here, we propose a set of 3-D turbulence-dependent resistance formulations for particle dry deposition simulation and intercompare the performance of new resistance formulations with that obtained from using the existing formulations and measured dry deposition velocity. Turbulence parameters such as turbulence velocity scale, turbulence factor, intensity of turbulence, effective sedimentation velocity, and effective Stokes number are newly introduced into two different particle deposition schemes to improve turbulence strength representation. For an assumed particle size distribution, the newly proposed schemes predict stronger diurnal variation of particle dry deposition velocity and are comparable to corresponding measurements while existing formulations indicate large underpredictions. We also find that the incorporation of new turbulence parameters either introduced or added stronger diurnal variability to sedimentation velocity and collection efficiencies values, making the new schemes predict higher deposition values during daytime and nighttime when compared to existing schemes. The findings from this research may help improve the capability of dry deposition schemes in regional and global models.

Dry Deposition Methods Based on Turbulence Kinetic Energy: Part 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single-Point Model.

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. To address these issues with dry deposition, firstly we take advantage of three-dimensional (3-D) turbulence information in estimating resistances by proposing and validating a 3-D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3-D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (1) present a set of resistance formulations for ozone (O3) based on the 3-D turbulence velocity scale; (2) intercompare estimations of such resistances with those obtained using existing formulations; and, (3) evaluate simulated O3 fluxes using a single-point dry deposition model against long-term observations of O3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

Impact of scale-aware deep convection on the cloud liquid and ice water paths and precipitation using the Model for Prediction Across Scales (MPAS-v5.2).

The cloud liquid water path (LWP), ice water path (IWP), and precipitation simulated with uniform- and variable-resolution numerical experiments using the Model for Prediction Across Scales (MPAS) are compared against Clouds and the Earth's Radiant Energy System (CERES) and Tropical Rainfall Measuring Mission data. Our comparison between monthly-mean model diagnostics and satellite data focuses on the convective activity regions of the tropical Pacific Ocean, extending from the Tropical Eastern Pacific Basin where trade wind boundary layer clouds develop to the Western Pacific Warm Pool characterized by deep convective updrafts capped with extended upper-tropospheric ice clouds. Using the scale-aware Grell-Freitas (GF) and Multiscale Kain-Fritsch (MSKF) convection schemes in conjunction with the Thompson cloud microphysics, uniform-resolution experiments produce large biases between simulated and satellite-retrieved LWP, IWP, and precipitation. Differences in the treatment of shallow convection lead the LWP to be strongly overestimated when using GF, while being in relatively good agreement when using MSKF compared to CERES data. Over areas of deep convection, uniform- and variable-resolution experiments overestimate the IWP with both MSKF and GF, leading to strong biases in the top-of-the-atmosphere longwave and shortwave radiation relative to satellite-retrieved data. Mesh refinement over the Western Pacific Warm Pool does not lead to significant improvement in the LWP, IWP, and precipitation due to increased grid-scale condensation and upward vertical motions. Results underscore the importance of evaluating clouds, their optical properties, and the top-of-the-atmosphere radiation budget in addition to precipitation when performing mesh refinement global simulations.

Studying Scale Dependency of Aerosol Cloud Interactions using Multi-Scale Cloud Formulations.

The WRF-ACI model configuration is used to investigate the scale dependency of aerosol-cloud interactions (ACI) across the "grey zone" scales for grid and subgrid-scale clouds. The impacts of ACI on weather are examined across regions in the eastern and western U. S. at 36, 12, 4, and 1 km grid spacing for short-term periods during the summer of 2006. ACI impacts are determined by comparing simulations with current climatological aerosol levels to simulations with aerosol levels reduced by 90%. The aerosol-cloud lifetime effect is found to be the dominant process leading to suppressed precipitation in regions of the eastern U.S., while regions in the western U. S. experience offsetting impacts on precipitation from the cloud lifetime effect and other effects that enhance precipitation. Generally, the cloud lifetime effect weakens with decreasing grid spacing due to a decrease in relative importance of autoconversion compared to accretion. Subgrid-scale ACI are dominant at 36 km, while grid-scale ACI are dominant at 4 and 1 km. At 12 km grid spacing, grid-scale and subgrid-scale ACI processes are comparable in magnitude and spatial coverage, but random perturbations in grid-scale-ACI impacts make the overall grid-scale-ACI impact appear muted. This competing behavior of grid and subgrid-scale clouds complicate the understanding of ACI at 12 km within the current WRF modeling framework. The work implies including subgrid-scale-cloud microphysics and ice/mixed phase cloud ACI processes may be necessary in weather and climate models to study ACI effectively.

The Weather Research and Forecasting Model with Aerosol-Cloud Interactions (WRF-ACI): Development, Evaluation, and Initial Application.

The Weather Research and Forecasting Model with Aerosol-Cloud Interactions (WRF-ACI) is developed for studying aerosol effects on gridscale and subgrid-scale clouds using common aerosol activation and ice nucleation formulations and double-moment cloud microphysics in a scale-aware subgrid-scale parameterization scheme. Comparisons of both the standard WRF and WRF-ACI models' results for a summer season against satellite and reanalysis estimates show that the WRF-ACI system improves the simulation of cloud liquid and ice water paths. Correlation coefficients for nearly all evaluated parameters are improved, while other variables show slight degradation. Results indicate a strong cloud lifetime effect from current climatological aerosols increasing domain average cloud liquid water path and reducing domain average precipitation as compared to a simulation with aerosols reduced by 90%. Increased cloud-top heights indicate a thermodynamic invigoration effect, but the impact of thermodynamic invigoration on precipitation is overwhelmed by the cloud lifetime effect. A combination of cloud lifetime and cloud albedo effects increases domain average shortwave cloud forcing by ~3.0 W m-2. Subgrid-scale clouds experience a stronger response to aerosol levels, while gridscale clouds are subject to thermodynamic feedbacks because of the design of the WRF modeling framework. The magnitude of aerosol indirect effects is shown to be sensitive to the choice of autoconversion parameterization used in both the gridscale and subgrid-scale cloud microphysics, but spatial patterns remain qualitatively similar. These results indicate that the WRF-ACI model provides the community with a computationally efficient tool for exploring aerosol-cloud interactions.

Mitigation of severe urban haze pollution by a precision air pollution control approach.

Severe and persistent haze pollution involving fine particulate matter (PM2.5) concentrations reaching unprecedentedly high levels across many cities in China poses a serious threat to human health. Although mandatory temporary cessation of most urban and surrounding emission sources is an effective, but costly, short-term measure to abate air pollution, development of long-term crisis response measures remains a challenge, especially for curbing severe urban haze events on a regular basis. Here we introduce and evaluate a novel precision air pollution control approach (PAPCA) to mitigate severe urban haze events. The approach involves combining predictions of high PM2.5 concentrations, with a hybrid trajectory-receptor model and a comprehensive 3-D atmospheric model, to pinpoint the origins of emissions leading to such events and to optimize emission controls. Results of the PAPCA application to five severe haze episodes in major urban areas in China suggest that this strategy has the potential to significantly mitigate severe urban haze by decreasing PM2.5 peak concentrations by more than 60% from above 300 µg m-3 to below 100 µg m-3, while requiring ~30% to 70% less emission controls as compared to complete emission reductions. The PAPCA strategy has the potential to tackle effectively severe urban haze pollution events with economic efficiency.

Sensitivities of Summertime Mesoscale Circulations in the Coastal Carolinas to Modifications of the Kain-Fritsch Cumulus Parameterization.

Two mesoscale circulations, the Sandhills circulation and the sea breeze, influence the initiation of deep convection over the Sandhills and the coast in the Carolinas during the summer months. The interaction of these two circulations causes additional convection in this coastal region. Accurate representation of mesoscale convection is difficult as numerical models have problems with the prediction of the timing, amount, and location of precipitation. To address this issue, the authors have incorporated modifications to the Kain-Fritsch (KF) convective parameterization scheme and evaluated these mesoscale interactions using a high-resolution numerical model. The modifications include changes to the subgrid-scale cloud formulation, the convective turnover time scale, and the formulation of the updraft entrainment rates. The use of a grid-scaling adjustment parameter modulates the impact of the KF scheme as a function of the horizontal grid spacing used in a simulation. Results indicate that the impact of this modified cumulus parameterization scheme is more effective on domains with coarser grid sizes. Other results include a decrease in surface and near-surface temperatures in areas of deep convection (due to the inclusion of the effects of subgrid-scale clouds on the radiation), improvement in the timing of convection, and an increase in the strength of deep convection.

High reduction of ozone and particulate matter during the 2016 G-20 summit in Hangzhou by forced emission controls of industry and traffic.

Many regions in China experience air pollution episodes because of the rapid urbanization and industrialization over the past decades. Here we analyzed the effect of emission controls implemented during the G-20 2016 Hangzhou summit on air quality. Emission controls included a forced closure of highly polluting industries, and limiting traffic and construction emissions in the cities and surroundings. Particles with aerodynamic diameter lower than 2.5 µm (PM2.5) and ozone (O3) were measured. We also simulated air quality using a forecast system consisting of the two-way coupled Weather Research and Forecast and Community Multi-scale Air Quality (WRF-CMAQ) model. Results show PM2.5 and ozone levels in Hangzhou during the G-20 Summit were considerably lower than previous to the G-20 Summit. The predicted concentrations of ozone were reduced by 25.4%, whereas the predicted concentrations of PM2.5 were reduced by 56%.

Does temperature nudging overwhelm aerosol radiative effects in regional integrated climate models?

Nudging (data assimilation) is used in many regional integrated meteorology-air quality models to reduce biases in simulated climatology. However, in such modeling systems, temperature changes due to nudging could compete with temperature changes induced by radiatively active and hygroscopic short-lived tracers leading to two interesting dilemmas: when nudging is continuously applied, what are the relative sizes of these two radiative forces at regional and local scales? How do these two forces present in the free atmosphere differ from those present at the surface? This work studies these two issues by converting temperature changes due to nudging into pseudo radiative effects (PRE) at the surface (PRE_sfc), in troposphere (PRE_atm), and at the top of atmosphere (PRE_toa), and comparing PRE with the reported aerosol radiative effects (ARE). Results show that the domain-averaged PRE_sfc is smaller than ARE_sfc estimated in previous studies and this work, but could be significantly larger than ARE_sfc at local scales. PRE_atm is also much smaller than ARE_atm. These results indicate that appropriate nudging methodology could be applied to the integrated models to study aerosol radiative effects at continental/regional scales, but it should be treated with caution for local scale applications.

Attribution of the United States "warming hole": aerosol indirect effect and precipitable water vapor.

Aerosols can influence the climate indirectly by acting as cloud condensation nuclei and/or ice nuclei, thereby modifying cloud optical properties. In contrast to the widespread global warming, the central and south central United States display a noteworthy overall cooling trend during the 20(th) century, with an especially striking cooling trend in summertime daily maximum temperature (Tmax) (termed the U.S. "warming hole"). Here we used observations of temperature, shortwave cloud forcing (SWCF), longwave cloud forcing (LWCF), aerosol optical depth and precipitable water vapor as well as global coupled climate models to explore the attribution of the "warming hole". We find that the observed cooling trend in summer Tmax can be attributed mainly to SWCF due to aerosols with offset from the greenhouse effect of precipitable water vapor. A global coupled climate model reveals that the observed "warming hole" can be produced only when the aerosol fields are simulated with a reasonable degree of accuracy as this is necessary for accurate simulation of SWCF over the region. These results provide compelling evidence of the role of the aerosol indirect effect in cooling regional climate on the Earth. Our results reaffirm that LWCF can warm both winter Tmax and Tmin.

Chemical transport models: the combined non-local diffusion and mixing schemes, and calculation of in-canopy resistance for dry deposition fluxes.

BACKGROUND, AIM, AND SCOPE: Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. MATERIALS AND METHODS: The upward mixing rate forming the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. RESULTS: This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. DISCUSSION: In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants calculated with chemical and air-quality models. In addition, this scheme was also compared with a commonly used, local, eddy-diffusivity scheme. Simulated concentrations of NO2 by the COM scheme and new parameterization of the in-canopy resistance are closer to the observations when compared to those obtained from using the local eddy-diffusivity scheme. CONCLUSIONS: Concentrations calculated with the COM scheme and new parameterization of in-canopy resistance, are in general higher and closer to the observations than those obtained by the local, eddy-diffusivity scheme (on the order of 15-22%). RECOMMENDATIONS AND PERSPECTIVES: To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO2) were compared for the years 1999 and 2002. The comparison was made for the entire domain used in simulations performed by the chemical European Monitoring and Evaluation Program Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.

Development of a nonlocal convective mixing scheme with varying upward mixing rates for use in air quality and chemical transport models.

BACKGROUND, AIM, AND SCOPE: Asymmetrical convective non-local scheme (CON) with varying upward mixing rates is developed for simulation of vertical turbulent mixing in the convective boundary layer in air quality and chemical transport models. MATERIALS AND METHODS: The upward mixing rate form the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. RESULTS: This scheme provides a less rapid mass transport out of surface layer into other layers than other asymmetrical convective mixing schemes. DISCUSSION: In this paper, we studied the performance of a nonlocal convective mixing scheme with varying upward mixing in the atmospheric boundary layer and its impact on the concentration of pollutants calculated with chemical and air-quality models. This scheme was additionally compared versus a local eddy-diffusivity scheme (KSC). Simulated concentrations of NO(2) and the nitrate wet deposition by the CON scheme are closer to the observations when compared to those obtained from using the KSC scheme. CONCLUSIONS: Concentrations calculated with the CON scheme are in general higher and closer to the observations than those obtained by the KSC scheme (of the order of 15-20%). Nitrate wet deposition calculated with the CON scheme are in general higher and closer to the observations than those obtained by the KSC scheme. RECOMMENDATIONS AND PERSPECTIVES: To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO(2)) and nitrate wet deposition was compared for the year 2002. The comparison was made for the whole domain used in simulations performed by the chemical European Monitoring and Evaluation Programme Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.