<p><strong>Abstract.</strong> Organic aerosol (OA) concentrations are simulated over the Beijing-Tianjin-Hebei (BTH) region from 9 to 26 January, 2014 using the Weather Research and Forecasting model coupled with chemistry (WRF-CHEM), with the goal of examining the impact of heterogeneous HONO sources on SOA formation and the SOA formation from different pathways during wintertime haze days. The model generally performs well in simulating air pollutants and organic aerosols against measurements in BTH. Model results show that heterogeneous HONO sources substantially enhance the near-surface SOA formation, increasing regional average near-surface SOA concentration by about 46.3<span class="thinspace"></span>% during the episode. Oxidation and partitioning of primary organic aerosols treated as semi-volatile dominate the SOA formation, contributing 58.9<span class="thinspace"></span>% of the near-surface SOA mass in BTH. Irreversible uptake of glyoxal and methylglyoxal on aerosol surfaces constitutes the second most important SOA formation pathway during the episode, with SOA contribution increasing from 8.5<span class="thinspace"></span>% in non-haze conditions to 30.2<span class="thinspace"></span>% in haze conditions. Additionally, direct emissions of glyoxal and methylglyoxal from residential living sources contribute about 25.5<span class="thinspace"></span>% to the total SOA mass on average in BTH. Our study highlights the importance of heterogeneous HONO sources and primary residential emissions of glyoxal and methylglyoxal to SOA formation in winter over BTH.</p>