An aerosol-optical module based on Mie scattering theory has been implemented in the Nested Air Quality Prediction Modeling System (NAQPMS), and a new coupler has been developed to deal with the interaction between the mesoscale meteorology model WRF (Weather Research and Forecasting Model) and NAQPMS. The one-way off-line and two-way coupled WRF-NAQPMS models are compared to simulate the severe haze in the Beijing-Tianjin-Hebei area from 27 September to 1 October 2013. The results show that the simulated meteorological elements and PM2.5 concentrations from the two-way coupled model with the aerosol direct radiation effect are more consistent with observations. During the haze period, the boundary layer meteorological elements change significantly because of the aerosol direct radiation effect over the Beijing-Tianjin-Hebei area: Incoming solar radiation is reduced by 25%, the 2-m temperature decreases by 1 ℃, the turbulent kinetic energy is reduced by 25%, the 10-m wind speed decreases by up to 0.2 m/s, and the planetary boundary layer (PBL) height is reduced by 25%. These changes make the atmospheric boundary layer more stable and further exacerbate air pollution over the areas where it is already severe, for example, the PM2.5 concentration increases by up to 30% over Shijiazhuang City. The analysis indicates that there is a positive feedback mechanism between haze and boundary layer meteorology, and the two-way coupled model incorporating this feedback is helpful for accurate simulation and forecasting of haze pollution processes.