Yang Zhou, Joshua S. North, Alan M. Rhoades, Jing Tao, William Rudisill, Mark D. Risser, William D. Collins
This study investigates the factors influencing runoff response to atmospheric rivers (ARs) over the U.S. West Coast. We focused on runoff time series variations impacted by AR characteristics (e.g., category and frequency) and land preconditions during Northern Hemisphere cool seasons in the period of 1940–2023. Results show that high-category ARs significantly increase local runoff with higher hourly precipitation rates leading to a greater incremental rate and peak runoff. Extreme runoff increases greatly with the AR category with an increase rate up to 12.5 times stronger than non-extreme runoff. Besides the AR category, land preconditions such as soil moisture and snowpack also play crucial roles in modulating runoff response. We found that runoff induced by weak-category ARs is more sensitive to land preconditions than high-category ARs, with high peak runoff occurring when soil is nearly saturated. Additionally, more than 50% of high-peak-runoff events in snow-covered grid cells are associated with rain-on-snow events particularly for the events associated with weaker ARs. Regression analysis reveals that AR precipitation and land preconditions jointly influence runoff, emphasizing the importance of including soil moisture and snowpack levels in AR impact assessments. The study also highlights the intensified runoff response to back-to-back ARs with short intervals, which may become more frequent with climate warming, posing increased flood risks via facilitating wet soil conditions. Our findings have significant implications for AR risk predictions and the development of prediction models for AR-induced runoff.
{"title":"Atmospheric River Frequency-Category Characteristics Shape U.S. West Coast Runoff","authors":"Yang Zhou, Joshua S. North, Alan M. Rhoades, Jing Tao, William Rudisill, Mark D. Risser, William D. Collins","doi":"10.1029/2024JD041805","DOIUrl":"https://doi.org/10.1029/2024JD041805","url":null,"abstract":"<p>This study investigates the factors influencing runoff response to atmospheric rivers (ARs) over the U.S. West Coast. We focused on runoff time series variations impacted by AR characteristics (e.g., category and frequency) and land preconditions during Northern Hemisphere cool seasons in the period of 1940–2023. Results show that high-category ARs significantly increase local runoff with higher hourly precipitation rates leading to a greater incremental rate and peak runoff. Extreme runoff increases greatly with the AR category with an increase rate up to 12.5 times stronger than non-extreme runoff. Besides the AR category, land preconditions such as soil moisture and snowpack also play crucial roles in modulating runoff response. We found that runoff induced by weak-category ARs is more sensitive to land preconditions than high-category ARs, with high peak runoff occurring when soil is nearly saturated. Additionally, more than 50% of high-peak-runoff events in snow-covered grid cells are associated with rain-on-snow events particularly for the events associated with weaker ARs. Regression analysis reveals that AR precipitation and land preconditions jointly influence runoff, emphasizing the importance of including soil moisture and snowpack levels in AR impact assessments. The study also highlights the intensified runoff response to back-to-back ARs with short intervals, which may become more frequent with climate warming, posing increased flood risks via facilitating wet soil conditions. Our findings have significant implications for AR risk predictions and the development of prediction models for AR-induced runoff.</p>","PeriodicalId":15986,"journal":{"name":"Journal of Geophysical Research: Atmospheres","volume":"130 2","pages":""},"PeriodicalIF":3.8,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2024JD041805","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115763","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Heatwave preconditioned-heavy rainfall (HW_HR), a preconditioned compound event, can cause more damage than a single heatwave or rainstorm. Both heatwaves and rainstorms can be exacerbated by the presence of cities, but the response of their compounding to urbanization remains unclear especially at the hourly scale. Here, we investigate the spatial and temporal responses of hourly HW_HR to a typical urban agglomeration, the Pearl River Delta, using observations and scenario-based numerical simulations. Compared to rural areas, the observations show that HW_HR in urban areas has a higher probability of occurrence and mean intensity, and its diurnal cycle of frequency is narrower, peaking in the afternoon. The temporal and spatial response of HW_HR intensity to urbanization effects is the most significant, with the urbanization-induced increase in HW_HR being five times that in non-heatwave preconditioned-heavy rainfall (noHW_HR). Our simulations support the observations and suggest that urbanization-induced changes are intense and spatially heterogeneous in HW_HR but relatively weak and spatially continuous in noHW_HR. The simulations also suggest that heatwave preconditioning not only amplifies urbanization-induced changes in key variables that alter atmospheric conditions but also provides a pre-storm unstable environment for the urban-induced warm-dry surface to trigger and enhance convection. The sub-daily pre-storm environment suggests that the preconditioning-induced thermodynamic changes gradually decrease, whereas the dynamic changes gradually increase as the event approaches. Our study highlights the importance of understanding urbanization effects on preconditioned events, providing new insights into the role of preconditions in the urban water cycle.
{"title":"Urbanization Intensifies Heavy Hourly Rainfall Preconditioned by Heatwaves","authors":"Zifeng Deng, Gabriele Villarini, Zhaoli Wang, Xushu Wu, Zhaoyang Zeng, Chengguang Lai","doi":"10.1029/2024JD041184","DOIUrl":"https://doi.org/10.1029/2024JD041184","url":null,"abstract":"<p>Heatwave preconditioned-heavy rainfall (HW_HR), a preconditioned compound event, can cause more damage than a single heatwave or rainstorm. Both heatwaves and rainstorms can be exacerbated by the presence of cities, but the response of their compounding to urbanization remains unclear especially at the hourly scale. Here, we investigate the spatial and temporal responses of hourly HW_HR to a typical urban agglomeration, the Pearl River Delta, using observations and scenario-based numerical simulations. Compared to rural areas, the observations show that HW_HR in urban areas has a higher probability of occurrence and mean intensity, and its diurnal cycle of frequency is narrower, peaking in the afternoon. The temporal and spatial response of HW_HR intensity to urbanization effects is the most significant, with the urbanization-induced increase in HW_HR being five times that in non-heatwave preconditioned-heavy rainfall (noHW_HR). Our simulations support the observations and suggest that urbanization-induced changes are intense and spatially heterogeneous in HW_HR but relatively weak and spatially continuous in noHW_HR. The simulations also suggest that heatwave preconditioning not only amplifies urbanization-induced changes in key variables that alter atmospheric conditions but also provides a pre-storm unstable environment for the urban-induced warm-dry surface to trigger and enhance convection. The sub-daily pre-storm environment suggests that the preconditioning-induced thermodynamic changes gradually decrease, whereas the dynamic changes gradually increase as the event approaches. Our study highlights the importance of understanding urbanization effects on preconditioned events, providing new insights into the role of preconditions in the urban water cycle.</p>","PeriodicalId":15986,"journal":{"name":"Journal of Geophysical Research: Atmospheres","volume":"130 2","pages":""},"PeriodicalIF":3.8,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115593","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Using long-term reanalysis data, extreme cold events (ECEs) over East Asia are categorized into long-duration ECEs (L-ECEs, lasting more than 7 days) and short-duration ECEs (S-ECEs, lasting 3–7 days) according to the duration of ECEs during boreal winters from 1959/1960 to 2019/2020. Our results show that L-ECEs over East Asia exhibit greater intensity, higher frequency, and longer persistence than S-ECEs. On annual average, L-ECEs are 1.8 K colder, contribute over 60% of total ECE days, and dominate in longer durations, persisting beyond day 5 and lasting up to 12 days. L-ECEs in East Asia are attributed to the extension of the stratospheric polar vortex (SPV) toward the Atlantic-Euro region. This SPV extension tends to enhance the Atlantic-Euro trough, which is coupled with a strengthened and persistent Eurasian teleconnection pattern and a strengthened Ural ridge and East Asia trough during the decay of ECEs (after day 3). The strengthened East Asia trough, accompanied with cyclonic circulation, facilitates the transport of cold air from higher to lower latitudes. This process makes the ECEs more persistent and longer in duration, potentially developing into L-ECEs. The findings of this study provide valuable insights into the linkages between the extension of the SPV and ECEs over East Asia from a long-term climate statistics perspective instead of relying solely on case-by-case analysis. This study may offer a valuable new predictor for forecasting long-lasting ECE events over East Asia.
{"title":"Longer Duration of Extreme Cold Events Over East Asia Affected by the Stratospheric Polar Vortex Extension","authors":"Dingzhu Hu, Meichen Liu, Zhaoyong Guan, Wenshou Tian, Zhicong Yin","doi":"10.1029/2024JD042299","DOIUrl":"https://doi.org/10.1029/2024JD042299","url":null,"abstract":"<p>Using long-term reanalysis data, extreme cold events (ECEs) over East Asia are categorized into long-duration ECEs (L-ECEs, lasting more than 7 days) and short-duration ECEs (S-ECEs, lasting 3–7 days) according to the duration of ECEs during boreal winters from 1959/1960 to 2019/2020. Our results show that L-ECEs over East Asia exhibit greater intensity, higher frequency, and longer persistence than S-ECEs. On annual average, L-ECEs are 1.8 K colder, contribute over 60% of total ECE days, and dominate in longer durations, persisting beyond day 5 and lasting up to 12 days. L-ECEs in East Asia are attributed to the extension of the stratospheric polar vortex (SPV) toward the Atlantic-Euro region. This SPV extension tends to enhance the Atlantic-Euro trough, which is coupled with a strengthened and persistent Eurasian teleconnection pattern and a strengthened Ural ridge and East Asia trough during the decay of ECEs (after day 3). The strengthened East Asia trough, accompanied with cyclonic circulation, facilitates the transport of cold air from higher to lower latitudes. This process makes the ECEs more persistent and longer in duration, potentially developing into L-ECEs. The findings of this study provide valuable insights into the linkages between the extension of the SPV and ECEs over East Asia from a long-term climate statistics perspective instead of relying solely on case-by-case analysis. This study may offer a valuable new predictor for forecasting long-lasting ECE events over East Asia.</p>","PeriodicalId":15986,"journal":{"name":"Journal of Geophysical Research: Atmospheres","volume":"130 2","pages":""},"PeriodicalIF":3.8,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115594","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The summer precipitation in Northeast China (NEC) exerts vital effects on the local crop yield and food security of our country, however, the prediction skill of it is still limited as yet. This study investigated the impact of snow cover anomaly over the Russian Far East during spring (March–April) on the succeeding early summer (May–June) precipitation variability in NEC as well as the underlying physical mechanisms during the years 1961–2020 based on diagnostic analyses and model simulation. The results showed that, on the interannual timescale, the decreased (increased) spring snow cover over the Russian Far East is conducive to the increased (decreased) NEC early summer precipitation. Further analyses indicated that the negative anomaly of spring snow cover can persist into early summer, which then enhances the diabatic heating and thereby warms the overlying atmosphere via the snow-albedo effect. The warming atmosphere favors the intensified Okhotsk High and an associated prominent anticyclonic circulation over northeast Asia. The anomalous southeast winds on the southwest flank of the anticyclone could bring abundant water vapor to the NEC, induce moisture convergence and ascending motion, and eventually contribute to the above-than-normal precipitation. These physical processes operate in the opposite manner in the case of the positive anomaly of spring snow cover. Our findings are beneficial to further improve the understanding and prediction ability of the summer precipitation variability over NEC.
{"title":"Impact of Spring Snow Cover Anomaly Over the Russian Far East on the Early Summer Precipitation Variability in Northeast China","authors":"Taotao Zhang, Xiaoyi Wang, Yaoming Song, Haishan Chen","doi":"10.1029/2024JD042118","DOIUrl":"https://doi.org/10.1029/2024JD042118","url":null,"abstract":"<p>The summer precipitation in Northeast China (NEC) exerts vital effects on the local crop yield and food security of our country, however, the prediction skill of it is still limited as yet. This study investigated the impact of snow cover anomaly over the Russian Far East during spring (March–April) on the succeeding early summer (May–June) precipitation variability in NEC as well as the underlying physical mechanisms during the years 1961–2020 based on diagnostic analyses and model simulation. The results showed that, on the interannual timescale, the decreased (increased) spring snow cover over the Russian Far East is conducive to the increased (decreased) NEC early summer precipitation. Further analyses indicated that the negative anomaly of spring snow cover can persist into early summer, which then enhances the diabatic heating and thereby warms the overlying atmosphere via the snow-albedo effect. The warming atmosphere favors the intensified Okhotsk High and an associated prominent anticyclonic circulation over northeast Asia. The anomalous southeast winds on the southwest flank of the anticyclone could bring abundant water vapor to the NEC, induce moisture convergence and ascending motion, and eventually contribute to the above-than-normal precipitation. These physical processes operate in the opposite manner in the case of the positive anomaly of spring snow cover. Our findings are beneficial to further improve the understanding and prediction ability of the summer precipitation variability over NEC.</p>","PeriodicalId":15986,"journal":{"name":"Journal of Geophysical Research: Atmospheres","volume":"130 2","pages":""},"PeriodicalIF":3.8,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115117","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Inorganic nitrate (NO3−), a crucial component of fine particulate matter (PM2.5), has not shown a consistent decrease, despite an obvious decrease of nitrogen oxide (NOx) and PM2.5. The atmospheric oxidation process for nitrate formation has been deemed a key factor in pollution; however, the changes of sources and oxidation pathways during a particular haze episode require further investigation. Here, daily dual isotopes (δ15N and δ18O) were used to quantify the sources and oxidation pathways of nitrate formation in Qingdao, a port city in Northern China, from September 2017 to February 2018. This study also includes a detailed introduction to two haze episodes. δ15N and δ18O results show that both fractions of nocturnal oxidation for nitrate formation and NOx from coal combustion were lower in warmer season and higher in colder season. The fractions increased with increasing PM2.5 under low PM2.5 concentration while the fractions were not significantly changed under higher PM2.5 concentration, dominated by nocturnal oxidation (70.6% ± 9.7%) and coal combustion (66.1% ± 18.2%), respectively. The haze episode 1 was attributed to smoke transported over long distances, which provided a large amount of aerosol particles to absorb more locally formed gaseous HNO3 or N2O5. In this episode, meteorological and air quality factors, nitrate sources, and formation mechanism did not obviously change. Haze episode 2 was caused by unfavorable meteorological factors that enhanced local nitrate accumulation. As pollution worsened, the oxidation pathway shifted from OH oxidation to N2O5 hydrolysis, and the primary source changed from coal combustion to more vehicle exhaust.
{"title":"Changes in the Dominant Contributions of Nitrate Formation and Sources During Haze Episodes: Insights From Dual Isotopic Evidence","authors":"Hong-Wei Xiao, Tian-Shu Chen, Qi-Jie Zhang, Rong Wang, Hao Xiao, Yu Xu, Wen-Kai Guan, Ai-Min Long, Hua-Yun Xiao","doi":"10.1029/2024JD042175","DOIUrl":"https://doi.org/10.1029/2024JD042175","url":null,"abstract":"<p>Inorganic nitrate (NO<sub>3</sub><sup>−</sup>), a crucial component of fine particulate matter (PM<sub>2.5</sub>), has not shown a consistent decrease, despite an obvious decrease of nitrogen oxide (NOx) and PM<sub>2.5</sub>. The atmospheric oxidation process for nitrate formation has been deemed a key factor in pollution; however, the changes of sources and oxidation pathways during a particular haze episode require further investigation. Here, daily dual isotopes (δ<sup>15</sup>N and δ<sup>18</sup>O) were used to quantify the sources and oxidation pathways of nitrate formation in Qingdao, a port city in Northern China, from September 2017 to February 2018. This study also includes a detailed introduction to two haze episodes. δ<sup>15</sup>N and δ<sup>18</sup>O results show that both fractions of nocturnal oxidation for nitrate formation and NOx from coal combustion were lower in warmer season and higher in colder season. The fractions increased with increasing PM<sub>2.5</sub> under low PM<sub>2.5</sub> concentration while the fractions were not significantly changed under higher PM<sub>2.5</sub> concentration, dominated by nocturnal oxidation (70.6% ± 9.7%) and coal combustion (66.1% ± 18.2%), respectively. The haze episode 1 was attributed to smoke transported over long distances, which provided a large amount of aerosol particles to absorb more locally formed gaseous HNO<sub>3</sub> or N<sub>2</sub>O<sub>5</sub>. In this episode, meteorological and air quality factors, nitrate sources, and formation mechanism did not obviously change. Haze episode 2 was caused by unfavorable meteorological factors that enhanced local nitrate accumulation. As pollution worsened, the oxidation pathway shifted from OH oxidation to N<sub>2</sub>O<sub>5</sub> hydrolysis, and the primary source changed from coal combustion to more vehicle exhaust.</p>","PeriodicalId":15986,"journal":{"name":"Journal of Geophysical Research: Atmospheres","volume":"130 2","pages":""},"PeriodicalIF":3.8,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115116","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Amit Singh Chandel, Chandan Sarangi, Karl Rittger, Rakesh K. Hooda, Antti-Pekka Hyvärinen
In this study, we used satellite observations to identify 10 typical dust-loading events over the Indian Himalayas. Next, the aerosol microphysical and optical properties during these identified dust storms are characterized using cotemporal in situ measurements over Mukteshwar, a representative site in Indian Himalayas. Relative to the background values, the mass of coarse particles (size range between 2.5 and 10 μm) and the extinction coefficient were found to be enhanced by 400% (from 24 ± 15 to 98 ± 40 μg/m3) and 175% (from 89 ± 57 Mm−1 to 156 ± 79 Mm−1), respectively, during these premonsoonal dust-loading events. Moreover, based on the air mass trajectory, these dust storms can be categorized into two categories: (a) mineral dust events (MDEs), which involve long-range transported dust plumes traversing through the lower troposphere to reach the Himalayas and (b) polluted dust events (PDEs), which involve short-range transported dust plumes originating from the arid western regions of the Indian subcontinent and traveling within the heavily polluted boundary layer of the Gangetic plains before reaching the Himalayas. Interestingly, compared to the background, the SSA and AAE decrease during PDEs but increase during MDEs. More importantly, we observe a twofold increase in black carbon concentrations and the aerosol absorption coefficient (relative to the background values) during the PDEs with negligible changes during MDEs. Consequently, the aerosol-induced snow albedo reduction (SAR) also doubles during MDEs and PDEs relative to background conditions. Thus, our findings provide robust observational evidence of substantial dust-induced snow and glacier melting over the Himalayas.
{"title":"In Situ Characterization of Dust Storms and Their Snow-Darkening Effect Over Himalayas","authors":"Amit Singh Chandel, Chandan Sarangi, Karl Rittger, Rakesh K. Hooda, Antti-Pekka Hyvärinen","doi":"10.1029/2024JD041874","DOIUrl":"https://doi.org/10.1029/2024JD041874","url":null,"abstract":"<p>In this study, we used satellite observations to identify 10 typical dust-loading events over the Indian Himalayas. Next, the aerosol microphysical and optical properties during these identified dust storms are characterized using cotemporal in situ measurements over Mukteshwar, a representative site in Indian Himalayas. Relative to the background values, the mass of coarse particles (size range between 2.5 and 10 μm) and the extinction coefficient were found to be enhanced by 400% (from 24 ± 15 to 98 ± 40 μg/m<sup>3</sup>) and 175% (from 89 ± 57 Mm<sup>−1</sup> to 156 ± 79 Mm<sup>−1</sup>), respectively, during these premonsoonal dust-loading events. Moreover, based on the air mass trajectory, these dust storms can be categorized into two categories: (a) mineral dust events (MDEs), which involve long-range transported dust plumes traversing through the lower troposphere to reach the Himalayas and (b) polluted dust events (PDEs), which involve short-range transported dust plumes originating from the arid western regions of the Indian subcontinent and traveling within the heavily polluted boundary layer of the Gangetic plains before reaching the Himalayas. Interestingly, compared to the background, the SSA and AAE decrease during PDEs but increase during MDEs. More importantly, we observe a twofold increase in black carbon concentrations and the aerosol absorption coefficient (relative to the background values) during the PDEs with negligible changes during MDEs. Consequently, the aerosol-induced snow albedo reduction (SAR) also doubles during MDEs and PDEs relative to background conditions. Thus, our findings provide robust observational evidence of substantial dust-induced snow and glacier melting over the Himalayas.</p>","PeriodicalId":15986,"journal":{"name":"Journal of Geophysical Research: Atmospheres","volume":"130 2","pages":""},"PeriodicalIF":3.8,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115353","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
<p>δ<sup>15</sup>N-NO<sub>3</sub><sup>−</sup> is widely used to trace the NO<sub>x</sub>/NO<sub>3</sub><sup>−</sup> emission sources without unique source tracers. However, there is still controversy regarding the <sup>15</sup>N fractionation effects during NO<sub>3</sub><sup>−</sup> formation, leading to uncertain source apportionment. To address this, this study introduces dual oxygen isotopes (∆<sup>17</sup>O and δ<sup>18</sup>O) to constrain the <sup>15</sup>N fractionation (∆<sup>15</sup>N-∆<sup>17</sup>O/∆<sup>15</sup>N-δ<sup>18</sup>O) of NO<sub>3</sub><sup>−</sup> formation and compare the impact of δ<sup>15</sup>N-NO<sub>x</sub> (∆<sup>17</sup>O) and δ<sup>15</sup>N-NO<sub>x</sub> (δ<sup>18</sup>O) on NO<sub>x</sub>/NO<sub>3</sub><sup>−</sup> source apportionment. Results found significant differences in ∆<sup>15</sup>N-∆<sup>17</sup>O (−3.7 ∼ +16.1‰) and ∆<sup>15</sup>N-δ<sup>18</sup>O (+8.5 ∼ +16.2‰) in haze, reflecting the ∆<sup>15</sup>N from three pathways (NO<sub>2</sub> + OH, NO<sub>3</sub> + HC, N<sub>2</sub>O<sub>5</sub> hydrolysis) and two pathways (NO<sub>2</sub> + OH and N<sub>2</sub>O<sub>5</sub> hydrolysis), respectively. The <sup>15</sup>N fractionation value differences obtained by dual oxygen isotopes increases with the increase of NO<sub>3</sub> + HC contribution (0.02–0.65). Additionally, different results of NO<sub>x</sub>/NO<sub>3</sub><sup>−</sup> sources apportionment were obtained by δ<sup>15</sup>N-NO<sub>x</sub>(<span></span><math>