Search

Air quality over the North China Plain (NCP) and rainfall over southern China exhibit robust and correlated interannual variations during winter from 1979 to 2018. A correlation coefficient of 0.7 between the PM2.5 concentration over the NCP and rainfall in southern China indicates that poor air quality and rainfall are linked (referred to as the south rainfall–north pollution pattern), which is likely due to the strong eastern Pacific (EP) El Niño. In the mature phase of the strong EP El Niño, the Northeast Asian anomalous anticyclone (NAAA) is strengthened by two Rossby wave trains modulated by warmer sea surface temperatures. One wave train is directly related to the strong EP El Niño and originates in the tropical eastern Pacific, propagating eastward to East Asia along the middle and high latitudes, which weakens the Ural ridge and strengthens the NAAA. The other wave train is derived from the northern Indian Ocean and propagates into East Asia along the great circle route, which strengthens an anomalous cyclone in southern China and the NAAA. This results in an anomalous horizontal trough of 850 hPa being captured in eastern China at ∼30°N. The southwesterly flow along the southern margin of the trough transports abundant moisture to southern China, which leads to heavy rainfall in southern China combined with an anomalous ascending motion related to the anomalous cyclone. The southeasterly flow to the north of the trough weakens the East Asia winter monsoon due to a weakened Ural ridge, which then leads to a high PM2.5 concentrations over the NCP with the help of an anomalous descending motion related to the NAAA.

The cold air outbreak (CAO) is the most important way to reduce air pollution during the winter over northern and eastern China. However, a rapid reappearance of air pollution is usually observed during its decay phase. Is there any relationship between the reappearance of air pollution and the properties of CAO? To address this issue, we investigated the possible connection between air pollution reappearance and CAO by quantifying the properties of the residual cold air mass after CAO. Based on the analyses of recent winters (2014–2022), we found that the rapid reappearance of air pollution in the CAO decay phase has an occurrence frequency of 63 %, and the air quality in more than 50 % of CAOs worsens compared to that before CAO. The reappearance of air pollution tends to occur in the residual cold air mass with a weak horizontal flux during the first 2 d after CAO. By categorizing the CAOs into groups of rapid and slow air pollution reappearance, we found that the residual cold air mass with a moderate depth of 150–180 hPa, a large negative heat content, and small slopes of isentropes is favorable for the rapid reappearance of air pollution. Among these factors, the cold air mass depth is highly consistent with the mixing layer height, below which most air pollutants are found; the negative heat content and slope of isentropes in the cold air mass jointly determine the intensity of low-level vertical stability. The rapid reappearance of air pollution is also attributed to the maintenance of the residual cold air mass and the above conditions, which are mainly regulated by the dynamic transport process rather than diabatic cooling or heating. Furthermore, analysis of the large-scale circulation of CAOs in their initial stage shows that the anticyclonic (cyclonic) pattern in northern Siberia (northeastern Asia) can be recognized as a precursor for the rapid (slow) reappearance of air pollution after the CAO.

New particle formation (NPF) and subsequent particle growth are important sources of condensation nuclei (CN) and cloud condensation nuclei (CCN). While a number of observations have shown positive contributions of NPF to CCN at low supersaturation, negative NPF contributions were often simulated. Using the observations in a typical coastal city of Qingdao, we thoroughly evaluate the simulated number concentrations of CN and CCN using a NPF-explicit parameterization embedded in WRF-Chem model. In terms of CN, the initial simulation shows large biases of particle number concentrations at 10–40 nm (CN10–40) and 40–100 nm (CN40–100). By adjusting the process of gas-particle partitioning, including mass accommodation coefficient of sulfuric acid, the phase changes of primary organic aerosol emissions and the condensational amount of nitric acid, the concomitant improvement of the particle growth process yields a substantial reduction of overestimates of CN10–40 and CN40–100. Regarding CCN, SOA formed from the oxidation of semi-volatile and intermediate volatility organic vapors (SI-SOA) yield is an important contributor. In the original WRF-Chem model with 20 size bins setting, the yield of SI-SOA is too high without considering the differences in oxidation rates of the precursors. Lowering the SI-SOA yield results in much improved simulations of the observed CCN concentrations. On the basis of the bias-corrected model, we find substantial positive contributions of NPF to CCN at low supersaturation (~0.2 %) in Qingdao and over the broad areas of China, primarily due to the competing effects of increasing particle hygroscopicity surpassing that of particle size decrease. This study highlights the potentially much larger NPF contributions to CCN on a regional and even global basis.

The canonical view of the northeast Asian anomalous anticyclone (NAAA) is a crucial factor for determining poor air quality (i.e., higher particulate matter, PM2.5 concentrations) in the North China Plain (NCP) on the interannual timescale. However, there is considerable intraseasonal variability in the NAAA in early winter (November–January), and the corresponding mechanism of its impacts on PM2.5 pollution in the NCP is not well understood. Here, we find that the intraseasonal NAAA usually establishes quickly on day 3 prior to its peak day with a duration of 8 d, and its evolution is closely tied to the Rossby wave from upstream (i.e., the North Atlantic). Moreover, we find that the NAAA with a westward tilt might be mainly related to the wavenumbers 3–4. Further results reveal that against this background, the probability of regional PM2.5 pollution for at least 3 d in the NCP is as high as 69 % (80 % at least 2 d) in the Nov–Jan (NDJ) period 2000–2021. In particular, air quality in the NCP tends to deteriorate on day 2 prior to the peak day and reaches a peak on the next day with a life cycle of 4 d. In the course of PM2.5 pollution, a shallower atmospheric boundary layer and stronger surface southerly wind anomaly associated with the NAAA in the NCP appear 1 d earlier than poor air quality, which provides dynamic and thermal conditions for the accumulation of pollutants and finally occurrence of the PM2.5 pollution on the following day. Furthermore, we show that the stagnant air leading to poor air quality is determined by the special structure of temperature in the vertical direction of the NAAA, while weak ventilation conditions might be related to a rapid build-up of the NAAA. The present results quantify the impact of the NAAA on PM2.5 pollution in the NCP on the intraseasonal timescale.

The autumn Intertropical Convergence Zone (ITCZ) over the South China Sea (SCS) is typically held south of 10°N by prevailing northeasterly and weakening southwesterly winds. However, the ITCZ can move north, resulting in heavy rainfall in the northern SCS (NSCS). We investigate the mechanisms that drove the northward movement of the ITCZ and led to heavy non-tropical-cyclone rainfall over the NSCS in autumn of 2010. The results show that the rapid northward movement of the ITCZ on 1 and 2 October was caused by the joint influence of the equatorial easterlies (EE), southwesterly winds, and the easterly jet (EJ) in the NSCS. A high pressure center on the east side of Australia, strengthened by the quasi-biweekly oscillation and strong Walker circulation, was responsible for the EE to intensify and reach the SCS. The EE finally turned southeast and together with enhanced southwesterly winds associated with an anticyclone, pushed the ITCZ north. Meanwhile, the continental high moved east, which reduced the area of the EJ in the NSCS and made room for the ITCZ. Further regression analysis showed that the reduced area of the EJ and increased strength of the EE contributed significantly to the northward movement of the ITCZ. The enhancement of the EE preceded the northward movement of the ITCZ by six hours and pushed the ITCZ continually north. As the ITCZ approached 12°N, it not only transported warm moist air but also strengthened the dynamic field by transporting the positive vorticity horizontally and vertically which further contributed to the heavy rainfall.

Source apportionments of urban aerosols identified with positive matrix factorization (PMF) are sensitive to input variables. So far, total elements were frequently included as effective factors in PMF-based source apportionment. We investigated the advances to involve soluble parts of elements in the source apportionment with four data sets of PM

China has been experiencing severe ozone pollution problems in recent years. While a number of studies have focused on the ozone-pollution-prone regions such as the North China Plain, Yangtze River Delta, and Pearl River Delta regions, few studies have investigated the mechanisms modulating the interannual variability of ozone concentrations in Shandong Province, where a large population is located and is often subject to ozone pollution. By utilizing both the reanalysis dataset and regional numerical model (WRF-CMAQ), we delve into the potential governing mechanisms of ozone pollution in Shandong Province—especially over the major port city of Qingdao—during summer 2014–2019. During this period, ozone pollution in Qingdao exceeded the tier II standard of the Chinese National Ambient Air Quality (GB 3095-2012) for 75 days. From the perspective of meteorology, the high-pressure ridge over Baikal Lake and to its northeast, which leads to a relatively low humidity and sufficient sunlight, is the most critical weather system inducing high-ozone events in Qingdao. In terms of emissions, biogenic emissions contribute to ozone enhancement close to 10 ppb in the west and north of Shandong Province. Numerical experiments show that the local impact of biogenic emissions on ozone production in Shandong Province is relatively small, whereas biogenic emissions on the southern flank of Shandong Province enhance ozone production and further transport northeastward, resulting in an increase in ozone concentrations over Shandong Province. For the port city of Qingdao, ship emissions increase ozone concentrations when sea breezes (easterlies) prevail over Qingdao, with the 95th percentile reaching 8.7 ppb. The findings in this study have important implications for future ozone pollution in Shandong Province, as well as the northern and coastal areas in China.

Both low clouds and elevated ducts are common phenomena in the oceanic atmosphere. Low clouds affect elevated ducts by changing the structure of atmospheric temperature and humidity. However, due to the limitation of met-ocean measurements, research on them is still scattered. This paper presents the distribution of elevated ducts and clouds over the central Western Pacific Ocean (WPO) based on Global Position System (GPS) sounding data and Himawari-8 satellite products from November 2015 to January 2016. Results show that the frequency of elevated ducts detected by ship-based GPS soundings was as high as 77% over the central WPO. The height and frequency of elevated ducts are closely related to the low clouds. If there are no clouds, the occurrence probability and mean base height of the elevated ducts are 14% and 730 m, respectively. By comparison, the occurrence probability and mean base height increase up to 24% and 1471 m, respectively, in the presence of cumulus (Cu) clouds, and 22% and 1511 m, respectively, in the presence of stratocumulus (Sc) clouds. Elevated ducts occur near the cloud top. The analysis of geopotential height and wind fields from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset (ERA-interim) shows that the study area is covered by a strong and stable subtropical high, and slowly sinking dry air masses inside the subtropical high are above the moist boundary-layer air mass. The appearance and evolution of low clouds will adjust the temperature and humidity structure of the lower troposphere. If there are no clouds, the marine boundary layer (MBL) is the classic mixed boundary layer. Humidity gradient and subsidence inversion are formed atop the mixed layer. When low clouds are present, long wave radiation and entrainment atop clouds form a strong temperature inversion and humidity gradient, which strengthen elevated ducts. However, when Sc clouds are decoupled, a weaker temperature inversion and humidity gradient may occur between the surface mixed layer and subcloud layer, leading to a weak elevated duct atop the mixed layer.

The impact of sea surface temperature anomalies (SSTAs) over the North Pacific (32°–45° N, 140° E–150° W) and North Atlantic (north: 52°–68° N, 60°–20° W; south: 0°–30° N, 100°–40° W) on the summer surface air temperature (SAT) over the Mongolian Plateau (MP) is studied using NCEP/NCAR and Climatic Research Unit (CRU) reanalysis data. The results show that the circumglobal teleconnection (CGT) wave train related to the SSTAs in the North Atlantic propagates along with the westerly jet, resulting in a geopotential height anomaly at 200 hPa that favors an SAT anomalies (SATAs) on the MP. In addition, an anticyclonic (cyclonic) circulation anomaly with the characteristics of an atmospheric ultralong wave at 200 hPa over the Eurasia–Pacific region associated with the SSTAs in the North Pacific is also responsible for the positive (negative) SATAs on the MP. The results further reveal that the SSTAs in the North Pacific and the northern North Atlantic act synergistically on the SATAs on the MP. If there are strongly positive (negative) SSTAs in the North Pacific and the northern North Atlantic at the same time, the probability of a strongly positive (negative) SATAs on the MP is 67% (100%) for the research period of this paper. However, if the strongly positive (negative) SSTAs only appear in the North Pacific or the northern North Atlantic, the probability of a strongly positive (negative) SATAs on the MP is less than 34% (31%). The synergism of the SSTAs in the North Pacific and the northern North Atlantic enhances Rossby wave energy in Eurasia, with more Rossby wave energy spreading to the MP that in turn causes the circulation anomaly on the MP. Meanwhile, baroclinic conditions at 850 hPa north of the North Pacific, the northern and southern North Atlantic favor a northward shift of the westerly jet that guides the Rossby wave downstream so that the circulation anomaly that produces the SATAs on the MP is located in the mid–high latitudes.

The long retention of dust air masses in polluted areas, especially in winter, may efficiently change the physicochemical properties of aerosols, causing additional health and ecological effects. A large-scale haze-to-dust weather event occurred in the North China Plain (NCP) region during the autumn-to-winter transition period in 2018, affecting the coastal city Qingdao several times between Nov. 27th and Dec. 1st. To study the evolution of the pollution process, we analyzed the chemical characteristics of PM

The cold air outbreak (CAO) is the most important way to reduce air pollution during the winter over North China. However, a rapid reappearance of air pollution is usually observed during its decay phase. Is there any relationship between the reappearance of air pollution and the properties of CAO? To address this issue, we investigated the possible connection between air pollution reappearance and CAO by quantifying the properties of the residual cold airmass after CAO. Based on the analyses in recent winters (2014–2019), we found that the rapid reappearance of air pollution in the CAO decay phase has an occurrence frequency of 74 %, and the air quality in more than 50 % of CAOs worsens than that before CAO. The reappearance of air pollution tends to occur in the residual cold airmass with weak horizontal flux during the first two days after CAO. By categorizing the CAOs into groups of rapid and slow air pollution reappearance, we found that the residual cold airmass with moderate depth of 150–180 hPa, large negative heat content and small slopes of isentropes is favorable for the rapid reappearance of air pollution. Among these factors, the cold airmass depth is highly consistent with the mixing layer height, below which contains most air pollutants; the negative heat content and slope of isentropes in cold airmass jointly determine the intensity of low-level vertical stability. The rapid reappearance of air pollution is also attributed to the maintenance of residual cold airmass and the above conditions, which is mainly regulated by the dynamic transport process rather than diabatic cooling or heating. Furthermore, analysis of the large-scale circulation of CAOs in their initial stage shows that the anticyclonic (cyclonic) pattern in northern Siberia (northeastern Asia) can be recognized as a precursor for the rapid (slow) reappearance of air pollution after the CAO.

The canonical view of the northeast Asian anomalous anticyclone (NAAA) is a crucial factor for determining poor air quality in the North China Plain (NCP) on the interannual timescale. However, there is considerable intraseasonal variability in the NAAA in early winter (November to January), and the corresponding mechanism of its impacts on air pollution in the NCP is not well understood. Here, we find that the intraseasonal NAAA usually establishes quickly on day −3 with a life span of eight days, and its evolution is closely tied to the Rossby wave from upstream (i.e., the North Atlantic). Moreover, we find that the NAAA with a westward tilt might be mainly related to the wavenumbers 3−4. Further results reveal that under this background, the probability of regional air pollution for at least three days in the NCP is as high as 69 % (80 % at least two days) in NDJ period 2000−2021. In particular, air quality in the NCP tends to deteriorate on day 2 prior to the peak day of the NAAA and reaches a peak on day −1 with a life cycle of four days. In the course of air pollution, a shallower atmospheric boundary layer and stronger surface southerly wind anomaly associated with the NAAA in the NCP appear 1 day earlier than poor air quality, which provides dynamic and thermal conditions for the accumulation of pollutants and finally occurrence of the air pollution on the following day. Furthermore, we show that the stagnant air leading to poor air quality is determined by the special structure of temperature in the vertical direction of the NAAA, while weak ventilation conditions might be related to a rapid buildup of the NAAA. The present results quantify the impact of the NAAA on air pollution in the NCP on the intraseasonal timescale.

Cold-front systems provide scavenging mechanisms for air pollution in the North China Plain (NCP), but the transport of pollutants with cold fronts aggravates air quality downstream. The impact of cold fronts on PM

Previous studies on ozone pollution primarily focus on the characterization of ozone on a large regional scale, yet much less attention has been paid towards the contrast between urban and surrounding suburban-rural areas. As anthropogenic emissions are projected to decrease in the coming decades, the evolutions of ozone concentrations over urban and suburban-rural areas are compared using the Community Multiscale Air Quality (CMAQ) model coupled with the Weather Research and Forecasting (WRF). The top 25 city clusters are firstly identified across China based on the amount of NO

New particle formation (NPF) induces a sharp increase in ultrafine particle number concentrations and potentially acts as an important source of cloud condensation nuclei (CCN). As the densely populated area of China, the Yangtze River Delta (YRD) region shows a high frequency of observed NPF events at the ground level, especially in spring. Although recent observational studies suggested a possible connection between NPF at the higher altitudes and ground level, the role played by vertical mixing, particularly in the planetary boundary layer (PBL) is not fully understood. Here we integrate measurements in Nanjing on 15-20 April 2018, and the NPF-explicit Weather Research and Forecast coupled with chemistry (WRF-Chem) model simulations to better understand the governing mechanisms of the NPF and CCN. Our results indicate that newly formed particles at the boundary layer top could be transported downward by vertical mixing as the PBL develops. A numerical sensitivity simulation created by eliminating aerosol vertical mixing suppresses both the downward transport of particles formed at a higher altitude and the dilution of particles at the ground level. The resulting higher Fuchs surface area at the ground level, together with the lack of downward transport, yields a sharp weakening of NPF strength and delayed start of NPF therein. The aerosol vertical mixing, therefore, leads to a more than double increase of surface CN

Based on hourly high-density precipitation data in Guangdong Province, China, 134 warm-sector heavy rainfall (WSHR) events were selected from 2016 to 2018. The synoptic weather patterns of these WSHR events were objectively classified using T-mode principal component analysis. Six WSHR weather patterns were identified, as follows: Type 1-southwest (T1-SW), Type 2-southeast (T2-SE), Type 3-coastal jets Ⅰ (T3-CJI), Type 4-coastal jets Ⅱ (T4-CJ II), Type 5-western low vortex (T5-WL), and Type 6-high-pressure (T6-HP). Three high-occurrence WSHR centers were finally extracted: the areas of Yangjiang and Shanwei, and the urban agglomeration of Guangdong–Hong Kong–Macao Greater Bay Area (GBA). Compared with the other five patterns, T6-HP is a newly identified WSHR weather pattern, which is related to a local/small-scale weather system in the context of anomalous northward movement of the western Pacific subtropical high. Relative to events in the other five patterns, the area of this T6-HP type of WSHR event is smaller, which can only be captured by high-density observations. In addition, the occurrence locations of six large-scale extreme precipitation events were closely associated with the urban agglomerations in GBA, implying that urbanization plays an important role in extreme magnitudes of large-scale WSHR events and their occurrence centers.

In this study, a methodology for satellite retrieval of cloud condensation nuclei (CCN) in shallow marine boundary layer clouds is presented and validated. This methodology is based on retrieving cloud base drop concentration (Nd) and updrafts (Wb), which are used for calculating supersaturation (S). The Nd is the activated CCN concentration in clouds at a given S. The accuracy of the satellite retrieval is validated against the surface-measured CCN of a cruise campaign over the heavily polluted northwest Pacific Ocean. Clouds which are coupled with the sea surface have good agreement between satellite retrieved Nd and surface-measured CCN after performing corrections for temperature and adiabatic fraction. This study broadens the applicability of the methodology from aerosol-limited to contaminated regions. The validation shows ±30% accuracy in retrieving CCN of both clean and polluted regions. The results further demonstrate the strong dependence of marine shallow cloud Nd on CCN number concentrations and updraft, which allows us to further apply this methodology to quantify the relationships between CCN and cloud microphysical properties and reduce the uncertainty of radiation forcing caused by aerosol cloud interaction (ACI).

Attributing sources of air pollution events by deploying an efficient observational network is an important and interesting problem in air quality control and forecast studies, but it is very challenging. In order to estimate the sensitivities of pollution events to emission sources, a comprehensive framework is built based on a horizontal 2-dimensional transport model and its adjoint in solving this problem. In an analysis of an idealized air pollution event of PM2.5 over the region of North China, an objective function is defined to optimally estimate the initial concentrations and emission sources through a series of minimization procedures. Results by means of the 4-dimensional variational approach show that, with the optimal initial conditions and emission sources, the model can successfully forecast the pollution event in a few days. The optimal observing network based on sensitivity analysis takes only one third of the cost but greatly retains predictability skill compared to the full-grid observing system, while nearly no predictability skill is detectable if the same number of observational sites is randomly deployed. We evaluate air pollution predictability in the point of focusing on to what degree the root mean square errors between the modeled concentration and the targeted air pollution are limited by the optimal observational network. Results show that air pollution predictability in association with the optimal observational network is limited in the time scales about 6 days. With the high efficiency and in an economic fashion, such a sensitivity-based optimal observing system holds promise for accurately predicting an air pollution event in the targeted area once the adjoint and variational procedure of a realistic atmosphere model including transport and chemical processes is performed.

Trend turning (or trend change) is a type of structural change that is common in climate data, and methods for detecting it in time series with multiple turning‐points need to be developed. A recently developed method for this, the running slope difference (RSD) t‐test, examines trend differences in sub‐series of the sample time series to identify the trend turning‐points. In this paper, we use Monte Carlo simulation to evaluate this method's detection ability. Evaluation results show the method to be an effective tool for detecting trend turning time series and identify three major advantages of the RSD t‐test: ability to detect multiple turning‐points, capacity to detect all three types of trend turning, and great performance of reducing false alarm rate.

Global mean surface air temperature (GMT) rose roughly 0.85 °C from 1880 to 2012 (IPCC 2013), attributing mainly to an increase in atmospheric greenhouse gases. For different decadal timescale periods in the past 100 years, the warming rate of different periods may significantly different. For example, IPCC AR1 (1990) point out that GMT between 1910-1940 and 1975-1990 are significantly warming, meanwhile GMT stay nearly constant between 1940 and early 1970. The phenomenon of two nearby periods showing significantly different trends is knowing as trend turning, this phenomenon is common in climate time series and crucial when climate change is investigated. However, the available detection methods for climate trend turnings are relatively few, especially for the methods which have the ability of detecting multiple trend turnings. We propose a new methodology named as the running slope difference (RSD) t-test to detect multiple trend turnings. This method employs a t-distributed statistic of slope difference to test the sub-series trends difference of the time series, thereby identify the turning-points. We compare the RSD t-test method with some other existing trend turning detection methods with an idealized time series case and several climate time series cases. And we also report the Monte Carlo simulation used to evaluate this method’s detection ability. Results show that the RSD t-test method is an effective tool for detecting trend turning in time series, and this method has three major advantages: ability to detect multiple turning-points, capacity to detect all three types of trend turning, and great performance of avoiding false alarm.

In this study, the spatial pattern and temporal evolution of PM2.5 over North China Plain (NCP) and Northeast China (NEC) during 2014-2018 was investigated. The annual mean PM2.5 shows clear decreasing trends over time, but the seasonal mean PM2.5 as well as the seasonal total duration and frequency of haze days shows large inter-annual fluctuation. Based on the atmospheric stagnation index (ASI), this study examined the correlation between ASI and haze events over NCP and NEC. Detailed analysis indicates that location dependency exists of ASI in the capability of capturing the haze events, and the ability is limited in NCP. Therefore, we first propose two alternative methods in defining the ASI to either account for the lag effect or enlarge the threshold value of wind speed at 500 hPa. The new methods can improve the ability of ASI to explain the haze events over NEC, though marginal improvement was achieved in NCP. Furthermore, this study constructed the equation based on the boundary layer height and wind speed at 10-meter, apparently improving the ability in haze capture rate (HCR), a ratio of haze days during the stagnation to the total haze days. Based on a multi-model ensemble analyses under Representative Concentration Pathway (RCP) 8.5, we found that by the end of this century, climate change may lead to increases in both the duration and frequency of wintertime stagnation events over NCP. In contrast, the models predict a decrease in stagnant events and the total duration of stagnation in winter over NEC.

Nitrous acid (HONO) is an important reservoir of the hydroxyl radical (OH) and thus plays a central role in tropospheric chemistry. Exhaust from engines has long been known as a major primary source of HONO, yet most previous studies focused on vehicle emissions on land. In comparison, ship emissions of HONO have been rarely characterized, and their impacts on the tropospheric oxidation chemistry have not been quantified. In this study, we conducted cruise measurements of HONO and related species over the East China Sea. Contrasting air masses from pristine marine background air to highly polluted ship plumes were encountered. The emission ratio of ΔHONO/ΔNOx (0.51 ± 0.18%) was derived from a large number of fresh ship plumes. Using the in-situ measured emission ratio, a global ship emission inventory of HONO was developed based on the international shipping emissions of NOx in the Community Emission Data System inventory. The global shipping voyage emits approximately 63.9 ± 22.2 Gg yr−1 of HONO to the atmosphere. GEOS-Chem modelling with the addition of ship-emitted HONO showed that HONO concentrations could increase up to 40–100% over the navigation areas, leading to about 5–15% increases of primary OH production in the early-morning time. This study elucidates the potentially considerable effects of ship HONO emissions on the marine atmospheric chemistry, and calls for further studies to better characterize the ship emissions of HONO and other reactive species, which should be taken into account by global and regional models.