Effect of rainfall-induced diabatic heating over southern China on the formation of wintertime haze on the North China Plain

Abstract. During the winters (December–February) between 1985 and
2015, the North China Plain (NCP, 30–40.5∘ N,
112–121.5∘ E) suffered many periods of heavy haze, and these
episodes were contemporaneous with extreme rainfall over southern China;
i.e., south rainfall–north haze events. The formation of such haze events
depends on meteorological conditions which are related to the atmospheric
circulation associated with rainfall over southern China, but the underlying
physical mechanism remains unclear. This study uses observations and model
simulations to demonstrate that haze over the NCP is modulated by anomalous
anticyclonic circulation caused by the two Rossby wave trains, in
conjunction with the north–south circulation system, which ascends over
southern China, moves north into northern China near 200–250 hPa, and then
descends in the study area. Moreover, in response to rainfall heating,
southern China is an obvious Rossby wave source, supporting waves along the
subtropical westerly jet waveguide and finally strengthening anticyclonic
circulation over the NCP. Composite analysis indicates that these changes
lead to a stronger descending motion, higher relative humidity, and a weaker
northerly wind, which favors the production and accumulation of haze over
the NCP. A linear baroclinic model simulation reproduced the observed
north–south circulation system reasonably well and supports the diagnostic
analysis. Quasi-geostrophic vertical pressure velocity diagnostics were used
to quantify the contributions to the north–south circulation system made
by large-scale adiabatic forcing and diabatic heating (Q). The results
indicated that the north–south circulation system is induced mainly by
diabatic heating related to precipitation over southern China, and the
effect of large-scale circulation is negligible. These results provide the
basis for a more comprehensive understanding of the mechanisms that drive
the formation of haze over the NCP.



Introduction
Extensive heavy haze on the North China Plain (NCP, • E) has a detrimental effect on both human health and social activities (Chen et al., 2017;Hughes et al., 2018;Lelieveld et al., 2019;Li et al., 2019). These haze events are caused by emissions of pollutants combined with unfavorable meteorological conditions (Yang et al., 2016;Cai et al., 2017;Zhang et al., 2021). Although emissions play an important role in the generation of haze, numerous studies have suggested that meteorology is also a significant factor in the occurrence of extreme haze events (Quan et al., 2011;Wang et al., 2015;Gao et al., 2016;Stirnberg et al., 2021). For instance, Dang and Liao (2019) found that large interannual variations in the frequency and intensity of severe winter haze days were driven mainly by changes in meteorology. Chen and Wang (2015) also showed that the occurrence of severe haze events over northern China during the winter generally correlates with meteorological factors. Y. Q.  provided evidence that the accumulation of pollutants caused by unfavorable meteorological conditions have offset the decreases caused by emissions reductions during the COVID-19 lockdown, leading to the high aerosol concentrations over Beijing-Tianjin-Hebei that developed between 7 and 14 February 2020. Unfortunately, continued global warming will further increase the incidence of haze days in China by reducing the wind strength (Cai et al., 2017;Xu et al., 2019). Callahan and Mankin (2020) found that climate change will lead to hazefavorable conditions over Beijing becoming more frequent, but that internal variability can generate large uncertainties in these projections. There is little doubt that developing an improved understanding of the factors and mechanisms that cause haze is one of the greatest challenges facing researchers over the coming decades.
Overall, the role of meteorology modulated by the largescale circulation in the generation of haze is crucial but uncertain. Large-scale circulations regulate meteorological conditions like reducing dispersion by atmospheric teleconnection and facilitating the accumulation of haze pollutants (Y. . Hence, the predictability of wintertime heavy haze over the NCP stems from the underlying meteorological factors. These include, but are not limited to, Eurasian snow cover and central Siberian soil moisture ; the sea surface temperature (SST) anomalies related to El Niño-Southern Oscillation (ENSO) in the tropical Pacific (Feng et al., 2019;Yu et al., 2020; and the Atlantic oceans ; the Pacific Decadal Oscillation (Zhao et al., 2016); the Arctic Oscillation (Cai et al., 2017;; the preceding Antarctic Oscillation (i.e., August-September-October) (Z. ; and the North Atlantic Oscillation (Feng et al., 2019;M. Li et al., 2021), as well as Arctic sea ice changes (Zou et al., 2020), especially in the Beaufort Sea Figure 1. Schematic illustration of the local north-south circulation system. The circular red (blue) arrow indicates the anomalous cyclone (anticyclone), the vertical red (blue) arrow indicates ascending (descending) motion, the thick black arrow indicates northward movement in the upper-level troposphere, H indicates the anticyclonic anomaly, and the white cloud indicates rainfall (see An et al., 2020, for caption details). (Yin et al., 2019a;Li and Yin, 2020) and Chukchi Sea (Yin et al., 2019b).
Recent studies have further documented that heavy haze over the NCP can be attributed to an anticyclonic anomaly over northeastern Asia (AANA) caused by westerly jet wave trains in the middle-to-upper troposphere over the Eurasian continent Wang et al., 2019;. As a synoptic-scale circulation, this AANA is accompanied by anomalous southeasterly winds near the surface, as well as a temperature inversion layer and anomalous vertical motion in the surrounding areas, which encourages the development of severe haze (Zhong et al., 2019;. It is worth mentioning that  also found that the atmospheric circulation related to rainfall over southern China further supports the maintenance of heavy haze in the AANA background over the NCP via the local northsouth circulation system (Fig. 1). However, their interpretation was based only on isolated extreme heavy haze events during November and December 2015 . To determine whether this finding was simply a special case, further investigations will be required that are based on a greater number of haze events with the same background of atmospheric circulation as rainfall in southern China over a longer period. Consequently, in this study, we aim to investigate the mechanisms associated with haze formation over the NCP between 1985 and 2015, with an emphasis on the role of circulation related to rainfall-induced diabatic heating in southern China under the same background of atmospheric circulation. In addition, given the potential role of circulation related to rainfall over southern China in maintaining heavy haze over the NCP, we will also explore diabatic heating using a linear baroclinic model (LBM) run under heating forcing in southern China.
The remainder of the paper is organized as follows. The second section describes the datasets and methods used in this work. This is followed by a more extensive section describing the haze events and associated weather patterns. The fourth section introduces the role of diabatic heating caused by rainfall on haze formation over the NCP, which we simulated using the LBM. The paper concludes with a brief summary and discussion of the formation of haze over the NCP.

Data
We obtained quality-controlled in situ daily rainfall data from 194 stations in China covering the period January 1985 to February 2015 from the Chinese Meteorological Administration (http://www.nmic.cn/, last access: 12 November 2017) to determine the distribution of rainfall over southern China. Daily visibility data from meteorological stations in the region bounded by 15-55 • N, 105-135 • E, from January 1985 to December 2015 were also obtained from the Chinese Meteorological Administration. Such observed rainfall (Li and Sun, 2015;Ding and Li, 2017) and visibility (e.g., Liu et al., 2017; data have been widely used in previous research on extreme rainfall and haze events. In addition, the daily PM 2.5 concentration dataset for China for the period 1980-2019 was collated by Yang (2020) using data available at https://zenodo.org/ record/4293239#.YJI3J8DiuUn (last access: 12 April 2021). These daily PM 2.5 concentrations were in excellent agreement with ground measurements, with a coefficient of determination of 0.95 and a mean relative error of 12 % (H. . The atmospheric reanalysis data, including geopotential, zonal and meridional wind, relative humidity, air temperature, and vertical velocity, were obtained from the National Centers for Environmental Prediction -National Center for Atmospheric Research (NCEP-NCAR) NCEP-DOE AMIP-II Reanalysis (R-2) dataset (Kanamitsu et al., 2002), provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA (http://www.esrl.noaa.gov/psd/, last access: 20 December 2020). The data used in this study cover a 31-year period from 1985 to 2015, with a 1 d resolution, a horizontal spatial resolution of 2.5 • longitude × 2.5 • latitude, and vertical levels from 1000 to 100 hPa.

Methods
In accordance with the standards set by the Chinese Meteorological Administration (2010), we defined haze as a day on which the daily mean visibility and relative humidity were less than 10 km and 80 %, respectively, and when no rain, snow, or sand and dust storms occurred. This identification method has been applied to the forecast of haze-fog by Chi-nese Meteorological Administration and is particularly useful for studying haze (e.g., Liu et al., 2017;. To analyze the distribution of heating associated with rainfall, we calculated the atmospheric apparent heat source (Q 1 ) according to the equation obtained by Yanai et al. (1973): where c p denotes the specific heat at constant pressure, T is the air temperature, t is the time, ω is the vertical pressure velocity, the static stability σ = (RT/c p p) − (∂T /∂p), R is the gas constant, p is the pressure, V is the horizontal wind vector, ∇ is the horizontal gradient operator, L is the latent heat of condensation, and q is the specific humidity. Here, Q 1 represents the total diabatic heating (including radiation, latent heating, and surface heat flux) and the subgrid-scale heat flux convergences (Yanai et al., 1973). According to Nie et al. (2020), the quasi-geostrophic (QG) vertical pressure velocity (ω) diagnostics is can be used to decompose ω in extreme precipitation into one part (ω D ) due to large-scale adiabatic forcing (F ) and another part (ω Q ) due to diabatic heating (Q). The QG ω equation reads where σ = − RT p ∂ p lnθ is the dry static stability, and f is the Coriolis parameter. the terms Adv ζ = −V g ·∇ζ and Adv T = −V g · ∇T are the horizontal advection of geostrophic absolute vorticity (ζ ) and temperature (T ), respectively, by the geostrophic winds.
To explore the propagation of anomalous Rossby waves along the subtropical westerly jet waveguide in the Northern Hemisphere, the horizontal stationary wave activity flux can be calculated as follows (Takaya and Nakamura, 2001): where W is the wave activity flux with unit of m 2 s −2 , ψ(= /f ) is the geostrophic stream function, (gpm) is the geopotential height, and f (= 2 sin φ) is the Coriolis parameter. U(= (U, V ) T ; m s −2 ) is the basic flow. We used the daily reanalysis data, i.e., the daily zonal wind, meridional wind, and anomalous geopotential height (for the stream function), to calculate the vector W . Panels (a-n) represent events 1-13, respectively, as described in Table 1.
We calculated the linearized Rossby wave source according to Sardeshmukh and Hoskins (1988), which can be expressed as follows: Here, u = (u, v) denotes the horizontal wind velocity, ∇ H is the horizontal gradient, and the subscript χ(ψ) represents the divergent (rotational) component. Overbars indicate the climatological mean and primes signify anomalies. The LBM (Watanabe and Kimoto, 2000), a simple dry model, has been widely used to examine the steady linear response to idealized diabatic heating (Lu and Lin, 2009;Sampe and Xie, 2010;Xu et al., 2020;Hu et al., 2020). This model consists of basic equations linearized with respect to the mean state of the December (0)-February (1) (DJF) climatology from the NCEP-DOE reanalysis for 1981-2010. The version used here has a horizontal resolution of T42 (roughly equivalent to 2.8 • ) and 20 vertical sigma levels. Using the time integration methods, the LBM was run up to 30 d, and the variable (i.e., zonal wind, meridional wind, and vertical velocity) on the last day (i.e., day 30) was taken as the steady response to the prescribed diabatic heating over southern China.
All composite maps were obtained from the average of each individual case.

Weather patterns related to heavy haze events on the NCP
Ding and Li (2017) (Ding and Li, 2017). At the same time, significant haze was trapped on the NCP, with poor visibility (<10 km) over large areas, which is similar to the conditions described by .
In particular, more intense and widespread haze events occurred over the NCP on 22-24 December 2007, 1-5 January 1992, and 12-17 January 2012. A natural question that arises is what kind of atmospheric circulation regulates the weather phenomena associated with SR-NH events. Figure 3 shows the composite anomalous air temperature at 1000 hPa, the absolute values of relative humidity at 925 hPa, and the anomalous wind vectors at 1000 hPa for the 13 SR-NH events. The positive air temperature anomaly over the NCP, which is caused by warmer and humid airflow from the low-level western North Pacific and transported by the easterly wind anomaly (Fig. 3a and c), creates favorable moisture conditions for the hygroscopic growth of haze particles, and further deteriorates haze in the NCP . Additionally, relative humidity over China shows a triode pattern, with low values over the NCP (≤ 80 %), meaning that these haze cases are haze rather than haze-fog . However, the relative humidity over southern China is close to 100 %, and this is the result of the southwesterly airflow along the coast of southern China and the rainfall (Figs. 2, 3b, c, and 4b). The anomalous easterlies on the eastern side of the NCP not only transport warm and moist air that promotes the development of haze over the NCP, but also weaken the East Asian winter monsoon (An Figure 5. Composite map of the Rossby wave source (shading, unit: s −2 ) at 200 hPa for the 13 SR-NH events. According to Eq. (4), the variables should be composited first and then used to calculate the Rossby wave source; hence, there is no t test here. et al., 2020), which is not conducive to the horizontal diffusion of haze (Fig. 3c).
To illustrate the reasons for the above changes in the meteorological factors, we next consider large-scale circulation. The composite map of the meridional wind anomaly reveals a substantial wave train with alternate positive and negative meridional wind anomalies at 200 hPa over the midlatitudes of the Northern Hemisphere (Fig. 4a). This circulation pattern is suggestive of a Rossby wave train emanating from the North Atlantic (Li and Sun, 2015;Ding and Li, 2017;Huang et al., 2020). According to the Rossby wave theory (Hoskins and Ambrizzi,  1993), the subtropical westerly jet, as a waveguide, allows the Rossby wave to arrive at the NCP (Fig. 4a). The wave activity fluxes show this wave train extending eastwards from the North Atlantic to eastern Europe, then southeastwards to the Arabian Peninsula, the Tibetan Plateau, and finally to the NCP (Fig. 4a), indicating that atmospheric circulation over the NCP is regulated by this wave train. As a result, there is a prominent anticyclonic anomaly over the NCP (Fig. 4b), which agrees with the analysis of . At 500 hPa, the anomalies are similar to the negative phase of the Eurasian teleconnection (EU), although they are shifted south and west of the canonical position of the three centers of the EU pattern (Wallace and Gutzler, 1981) (Fig. 4b).
This negative EU-like pattern is not conducive to the development of the East Asian winter monsoon , and this further encourages the development of heavy haze over the NCP. In addition, the AANA related to EU-like pattern also supports haze development over the NCP via the anomalous descending motion and southerly wind along the east coast of China (Figs. 3c and 4b). In addition, the anoma-lous southwesterly airflow along the coast of southern China transports large amounts of water vapor to this area and is one of the conditions that leads to the heavy rainfall over southern China (Li and Sun, 2015;Ding and Li, 2017). Previous studies have demonstrated that the convergence and divergence anomalies in the upper troposphere can be regarded as an effective Rossby wave source for the stationary Rossby wave (Hoskins and Ambrizzi, 1993;Branstator, 2002;Watanabe, 2000;Chen et al., 2020). The strong divergence causes rainfall in southern China and is located at 200 hPa (Li and Sun, 2015;Ding and Li, 2017;An et al., 2020), but does it also act as the Rossby wave source to strengthen the Rossby wave along the subtropical westerly jet waveguide? The composite map of the Rossby wave source for SR-NH events reveals a significant negative Rossby wave source in the upper troposphere over southern China around 25-30 • N, 110-120 • E (Fig. 5), where positive precipitation and strong ascending motion anomalies are located (Figs. 2 and 6b). Therefore, the divergence anomalies in the upper troposphere (Fig. 7a) associated with positive precipitation and anomalous ascending motion (Figs. 2 and 6a, b) play a crucial role in the formation of the strong negative Rossby wave source over southern China, which excites the Rossby wave that propagates to the downstream regions. This means that there is a positive feedback process involving subtropical westerly jet waveguide rainfall over southern China. The rainfall over southern China is largely caused by the subtropical westerly jet waveguide (Li and Sun, 2015;Ding and Li, 2017;, which would in turn strengthen vertical motion and high-level divergence by an intensification of latent heat release, generating a Rossby wave source, and hence strengthen the wave train along the subtropical westerly jet. This process reinforces the Rossby wave that originated from the North Atlantic. Wang et al. (2019) and  found that such waves can account for haze over the NCP. The ascending motion over southern China that is related to the diabatic heating caused by the rainfall will be examined further below.

Possible physical mechanisms driving SR-NH events: the role of rainfall-induced diabatic heating
In the previous section, we presented a diagnostic analysis of the atmospheric circulation during periods of haze over the NCP. We now focus on the influence of circulation related to diabatic heating on haze formation over the NCP. The composite sections show when precipitation occurred in southern China, and that it was accompanied by a strong vertical ascending motion (Fig. 6a, b). At the same time, the NCP was controlled by an obvious descending motion (Fig. 6a,  c), which was related to the ascending motion over southern China and forms the local north-south circulation system together with the ascending motion in southern China (Fig. 6a). The descending motion over the NCP develops with the AANA as shown in Fig. 4b. The divergent wind, as well as the velocity potential, also confirmed that there was strong divergence near 200 hPa over southern China and strong convergence in northern China (Fig. 7a). The latitude-height cross section of the stream function and relative vorticity shows that there was strong convergence in the middle-to upper levels of the troposphere, but weak divergence (strong convergence) at lower (higher) levels over the NCP (Fig. 7b), and this was responsible for the descending motion and not conducive to the diffusion of haze. The centers of the convergence and divergence over southern China were the opposite to those over the NCP (Fig. 7b), which is the typical circulation pattern that accompanied rainfall (Li and Sun, 2015). The ascending motion over southern China is not only related to the subtropical westerly jet waveguide (Ding and Li, 2017) but may be also related to the diabatic heating caused by the rainfall. The diabatic heating in the atmosphere caused by precipitation may intensify the local ascending motion Xu et. al., 2020). To estimate the diabatic heating produced by rainfall extremes in southern China, we calculated Q 1 quantities using Eq. (3) as described above. Figure 8 shows the composite Q 1 during the period of heavy rainfall in southern China. An obviously positive Q 1 appears in southern China when rainfall occurs, with the maximum value of Q 1 near 26 • N, 115 • E (Fig. 8a). Shortwave radiation from the Sun rarely reaches the lower troposphere during rainfall because it is blocked by the clouds; therefore, a positive value for Q 1 indicates the rainfall process may release large amounts of heat. To investigate the vertical distribution of Q 1 , we calculated the average value of Q 1 value between 1000 and 100 hPa over the domain 20-30 • N, 110-120 • E. Figure 9b shows that the maximum value of Q 1 (2.6 K d −1 ) occurred at 500 hPa (Fig. 8b). The results shown in Fig. 9a also confirm that the rainfall acts as a heat source. Analyzing the distribution of Q 1 may help to identify the specific heating processes that are occurring in the atmosphere (Yanai et al., 1973). The above analysis corroborates that the rainfall process does release a lot of heat, which may further encourage ascending motion.
To further validate the above-mentioned rainfall-heating circulation, a numerical experiment was conducted. The results outlined above indicate that the local north-south circulation system seems to be closely associated with the diabatic heating over southern China. Using the LBM, we performed a numerical simulation to complement the observational results and test the plausibility of the proposed hazerainfall link. This numerical experiment simulated the atmospheric response to heat forcing induced by heavy rainfall in southern China (Fig. 9a). The experiment was run with diabatic heating centered over southern China (26 • N, 115 • E), which essentially matched the maximum value and variance of Q 1 , and the heavy rainfall located as shown in Fig. 8a. The maximum heating, with an amplitude of 2 K d −1 , was set at 500 hPa (Fig. 9b). Figure 9 illustrates the 300 hPa wind response to the diabatic heating over southern China. An obvious divergence in the wind vectors occurs over eastern China. The whole layer is divergent, which means that diabatic heating is conducive to upward motion, and the airflow diffuses northwards at 200 hPa.
The local north-south circulation system simulated by the LBM is similar to the observed spatial pattern of the local north-south circulation system (Figs. 6 and 10), both of which are broadly similar to the results obtained by An et al. (2020). As depicted in Fig. 11b, the upper-level negative and low-level positive relative vorticity over southern China (20-30 • N) indicates that the atmosphere is strongly baroclinic, with an anticyclonic circulation in the upper levels and a cyclonic circulation at lower levels over southern China. The direct product of this circulation is strong ascending motion over southern China. Meanwhile, upper-level convergence and low-level weak divergence over the NCP support descending motion there, which is conducive to haze. In addition, at 200 hPa, the positive relative vorticity located over the NCP reinforces the AANA (Fig. 11a). In summary, this LBM experiment further confirms that the circulation related to the rainfall over southern China plays an obvious assisting role in maintaining haze over the NCP.
From our earlier discussion, it is clear that the diabatic heat released by rainfall over southern China has a significant effect on the ascending motion, except for large-scale circulation background; i.e., the subtropical westerly jet waveguide and south branch trough suggested by Li and Sun (2015). However, a crucial issue that remains to be addressed is which contributes more to the vertical movement. In reality, in a moist atmosphere, the vertical motion depends not only on dynamic forcing by large-scale perturbations but also on the driving by the latent heating released from convection (Nie et al., 2020). To quantify the contribution of vertical motion due to rainfall heating, we calculated the components of the ω equation as described in Eq. (4). First, we find that negative (positive) ω values occur mainly over southern China (the NCP), which indicates that southern China is controlled by ascending (descending) motion (Fig. 12). Both the ascending motion over southern China and the descending motion over the NCP are strongest at 500 hPa (Fig. 12). The top panel in Fig. 12 presents the diabatic heating term, which is in agreement with that of ω, especially in the mid-high levels of the troposphere (i.e., 300 and 500 hPa), whereas the dry dynamic forcing term (middle panel in Fig. 12) clearly differs from ω, and its values in the lower level of troposphere (below 700 hPa) are very weak and unlikely to be responsible for enhancing vertical motion (Fig. 12). This means that the strengthened ascending motion caused by the subtropical westerly jet waveguide enhances the production of rainfall over southern China, resulting in increased diabatic heating, which in turn reinforces the ascending motion and hence forms a positive feedback that links the diabatic heating to the ascending motion when there is a plentiful supply of moisture. The above analysis again confirms that diabatic heating related to rainfall plays a crucial role in the ascending motion over southern China, which in turn maintains the local north-south circulation system, leading to a shallow atmospheric boundary layer height, finally aggravates the haze pollution over the NCP (not shown).  show the diabatic heating term R pf 2 ∇ 2 Q , the dry forcing term (F = 1 f ∂ p Adv ζ + R pf 2 ∇ 2 Adv T ), where Q is calculated using the Eq. (1), and ω (reanalysis data), respectively. The unit of ω is Pa s −1 . The unit of the F and Q terms is Pa −1 s −1 . For ease of viewing, and to reduce the differences in scale, the F , Q terms, and ω were multiplied by 4 × 10 7 , 0.8 × 10 7 , and 4, respectively. Dotted white regions indicate areas at the 95 % confidence level based on the two-tailed Student t test.   Figure 14. As Fig. 4 but for the six south rainfall-no north haze (SR-noNH) events in Table 2.

Discussion and conclusions
Our study investigated the mechanisms associated with production over the NCP related to rainfall heating over southern China based on a combination of observational composites and model results. We found that the appearance of some haze events over the NCP is the product of the circulation associated with rainfall over southern China in conjunction with two wave trains along the subtropical westerly jet and polar front jet waveguides. The extreme rainfall events over southern China analyzed in this study were caused mainly by the subtropical westerly jet waveguide (Li and Sun, 2015;Ding and Li, 2017) and the diabatic heat related to it induces the secondary circulation (referred to as the local northsouth circulation system in this paper). On the one hand, this heating strengthens the ascending motion over southern China and the descending motion over the NCP (Fig. 13).
On the other hand, this heating stimulates the wave train as the Rossby wave source and strengthens the background subtropical westerly jet wave train. These changes eventually lead to the strengthening of the anomalous anticyclone in northeast Asia, resulting in the strong descending movement, weak northerly wind, and warm and moisture-laden flow over the NCP. As a consequence, heavy haze is maintained over the NCP. When the AANA over the NCP moves a little to the east (Fig. 14b), heavy rainfall also occurs in southern China and triggers a similar the local north-south circulation system (not shown), but there is then no haze (i.e., visibility was >10 km) over the NCP (i.e., a SR-noNH event). This implies that the AANA is one of the important factors in the occurrence of severe haze events over the NCP (Fig. 4), and this agrees with the results of Chen et al. (2020) and . In addition, PM 2.5 concentrations were lower (most areas below 120 µg m −3 ) in the first 7 d during nine SR-noNH events (Fig. 15). This implies that emissions are the source of haze over the NCP, and meteorological factors worsen haze (e.g., Wang et al., 2015;Stirnberg et al., 2021). The shortcoming of this study is that there are not many samples (i.e., only 19 samples in present results), and further analysis might be needed from more samples in the future. A comprehensive understanding of why there is no haze in the NCP despite some rainfall in southern China (six SR-noNH events) also needs to be added in the future.
In addition to rainfall over southern China associated with the subtropical westerly jet waveguide, rainfall is also associated with ENSO (Ma et al., 2018). Therefore, ENSO may also be, except for the westerly jet waveguide, another potential factor in haze production over the NCP and rainfall over southern China. For instance, there has been rainfall over southern China and haze over the NCP during El Niño years (i.e., 1988, 1992, 2007, and 2015). The present study does not consider the role of ENSO with respect to rainfall over southern China and haze over the NCP, but this will be the focus of future work. The complex interconnections among atmospheric systems make it difficult to determine whether any single factor is the dominant control on haze production over the NCP, and it is possible that the synergistic effects of many influencing factors may play a more important role in the occurrence of haze. As mentioned in the introduction, there are many meteorological factors that play an important role in haze production in northern China. Previous indepth studies have considered various aspects of the formation mechanisms of haze in northern China; i.e., ENSO (Yu et al., 2020;, the Arctic Oscillation (Cai et al., 2017), and Arctic sea ice (Yin et al., 2019a, b;Li and Yin, 2020). However, whether such factors have a synergistic effect on haze production (e.g., Arctic sea ice and the tropical ocean) has not been considered here but would be a worthwhile focus for a future study. Data availability. The visibility observational data are available from the China Meteorological Administration (http://data.cma.cn/, CMA, 2017), the reanalysis dataset is available at NCEP/NCAR (https://www.esrl.noaa.gov/psd/data/gridded/, NCEP/NCAR, 2020), and the daily PM 2.5 concentration dataset for China for the period 1980-2019 is available at https://doi.org/10.5281/zenodo.4293239 (Yang, 2020).
Author contributions. XA, LS, and LC designed the experiments and carried them out. XD and YL developed the model code and performed the simulations. XA downloaded and analyzed the reanalysis data and prepared all the figures. XD prepared the manuscript with contributions from all co-authors. LS, WC, and JH revised the manuscript.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements. This research was supported by the National Natural Science Foundation of China (grant no. 41975008), Fundamental Research Funds for the Central Universities (grant no. 201961004), and National Key R&D Program of China (grant no. 2019YFA0607002). All the authors are very grateful to the China Meteorological Administration (http://data.cma.cn/, last access: 12 November 2017) and NCEP/NCAR (https://www.esrl. noaa.gov/psd/data/gridded/, last access: 20 December 2020) for data used in this study. Review statement. This paper was edited by Peter Haynes and reviewed by two anonymous referees.