Effects of black carbon and boundary layer interaction on surface ozone in Nanjing, China

As an important solar radiation absorbing aerosol, the effect of black carbon (BC) on surface ozone, via reducing photolysis rate, has been widely discussed by “offline” model studies. However, BC–boundary layer (BL) interactions also influence surface ozone. Using the “online” model simulations and process analysis, we demonstrate the significant impact of BC–BL interaction on surface ozone in Nanjing. The absorbing effect of BC heats the air above the BL and suppresses and delays the development of the BL, which eventually leads to a change in surface ozone via a change in the contributions from chemical and physical processes (photochemistry, vertical mixing and advection). For chemical processes, the suppression of the BL leads to large amounts of ozone precursors being confined below the BL which has an increased effect on ozone chemical production and offsets the decrease caused by the reduction of the photolysis rate, thus enhancing ozone chemical formation from 10:00 to 12:00 LT. Furthermore, changes in physical processes, especially the vertical mixing process, show a more significant influence on surface ozone. The weakened turbulence, caused by the suppressed BL, entrains much less ozone aloft down to the surface. Finally, summing-up the changes in the processes, surface ozone reduces before noon and the maximum reduction reaches 16.4 ppb at 12:00 LT. In the afternoon, the changes in chemical process are small which inconspicuously influence surface ozone. However, change in the vertical mixing process still influences surface ozone significantly. Due to the delayed development of the BL, there are obvious ozone gradients around the top of BL. Therefore, high concentrations of ozone aloft can still be entrained down to the surface which offsets the reduction of surface ozone. Comparing the changes in the processes, the change in vertical mixing plays the most important role in impacting surface ozone. Our results highlight the great impacts BC–BL interactions have on surface ozone by influencing the ozone contribution from physical process. This suggests that more attention should be paid to the mechanism of aerosol–BL interactions when controlling ozone pollution. Published by Copernicus Publications on behalf of the European Geosciences Union. 7082 J. Gao et al.: Effects of black carbon


Introduction
Black carbon (BC) aerosol, also known as soot, is primarily formed by incomplete combustion of carbonaceous fuels, diesel fuels, and biomass (Bond et al., 2004). BC accounts for a small fraction, less than 15 %, of the total mass concentration of aerosol particles in atmosphere over urban areas (Yang et al., 2011). However, it is of great interest because of its significant influence on global radiation balance (Chameides and Bergin, 2002), both directly, by absorbing solar radiation (Liao and Seinfeld, 2005;Ramanathan and Carmichael, 2008), and indirectly, by affecting cloud formation (Lohmann and Feichter, 1997;Fan et al., 2015). Owing to impacts, such as the abovementioned, on radiation, BC plays an important role in global and regional climate change (Jacobson, 2001;Bond et al., 2013), weather (Qian et al, 2003;Saide et al., 2015), and the atmospheric environment (Li et al., 2005;Ding et al., 2016;Peng et al., 2016). In addition, the BC aging mechanism (Qiu et al., 2012), the performance of BC in heterogeneous reactions (Lei et al., 2004;Li et al., 2012), and its impact on human health (Atkinson et al., 2012) have been the focus of significant research in recent years.
Tropospheric ozone is a typical secondary air pollutant (Crutzen, 1973). It also has important environmental effects on the atmosphere (Monks et al., 2015), especially in the boundary layer (BL). The impact of aerosols, especially BC, on surface ozone has been attracting much attention from researchers. Dickerson et al. (1997) reported that BC decreases surface ozone concentration by reducing photolysis rates. Jacobson (1998) suggested that aerosols containing BC cores reduced photolysis rates and resulted in a decrease of ozone concentration by 5-8 % at ground level in Los Angeles. Castro et al. (2001) found a strong reduction in photolysis rate (10-30 %) due to BC-containing aerosols. They also reported that this photolysis rate reduction led to a decrease of surface ozone in Mexico City. Similar results have been found in other studies simulating the effects of BC on surface ozone in other places around the world (Li et al., 2005(Li et al., , 2011. In addition to reducing photolysis rates, the global warming effect of BC is significant, preceded only by CO 2 (Jacobson, 2002). Incident solar radiation is absorbed by BC in the atmosphere, leading to the air aloft being heated and the temperature being raised. Conversely, air at low levels is cooled and the temperature is decreased. Under this condition, atmospheric stability is increased and the development of the BL is suppressed during the daytime. Using an atmospheric model, Ding et al. (2016) demonstrated this effect and suggested that such BC-BL interactions would enhance the occurrences of haze pollution episodes. Tie et al. (2017) also discussed the relationships between aerosol (including BC) and BL. Moreover, they highlighted the serious impacts of this mechanism on air quality in winter. Owing to the close relationship between ozone and BL development during the daytime (Zhang and Rao, 1999;Zhu et al., 2015), BC-BL interactions may also influence surface ozone; however, relevant studies on this phenomenon are still lacking.
As one of the most developed regions in China, the Yangtze River Delta (YRD) has reported severe haze (Ding et al., 2013a) and ozone pollution (Tie et al., 2013) in recent decades. Nanjing is the capital of Jiangsu Province, which is an important industrial and economic center in the YRD region. Previous studies have reported that BC (Zhuang et al., 2014) and ozone (Ding et al., 2013b) concentrations are relatively high in Nanjing in October. In this study, measured data of BC and ozone, which were obtained in Nanjing in October 2015, were used to show variations in surface ozone under high BC concentration and low BC concentration conditions during daytime. Furthermore, the fully coupled "online" model system Weather Research and Forecasting with Chemistry (WRF-Chem) was applied to simulate the air pollutants (ozone, PM 2.5 , and BC) in the YRD in October 2015.
With the consideration of the aerosol-BL feedback mechanism in the model system, we demonstrate the BC-BL interaction mechanism affecting surface ozone in Nanjing.
2 Model setting and data description WRF-Chem, which is widely used to evaluate the impacts of aerosols on radiation (Zhang et al., 2010;Forkel et al, 2012), is a fully coupled online 3-D Eulerian chemical transport model considering both chemical and physical processes. We used version 3.4 in this study, and detailed descriptions of the meteorological and chemical aspects of WRF-Chem can be found in Skamarock et al. (2008) and Grell et al. (2005). Regarding simulation settings, two nested domains ( Fig. 1) were set up with horizontal resolutions of 36 and 12 km, and grids of 99 × 99 and 99 × 99 for the parent domain (D1) and nested domain (D2), respectively. The parent domain (D1), centered at (119.0 • E, 31.5 • N), covered most of China and the surrounding countries and ocean. The corresponding simulations provided meteorological and chemical boundary conditions for the nested domain (D2), which covered most of eastern China.
There were 38 vertical layers from the surface up to a pressure limit of 50 hPa, 12 levels of which were located below the lowest 2 km to fully describe the vertical structure of the BL. Carbon Bond Mechanism Z (CBM-Z), which includes 133 chemical reactions for 53 species and extends the framework to function for a longer time and at larger spatial scales than its predecessor; Carbon Bond IV, was used as the gas-phase chemical mechanism (Zaveri and Peters, 1999). The corresponding aerosol chemical mechanism was the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with eight bins (Zaveri et al., 2008), which is extremely efficient and does not compromise accuracy in aerosol model calculations. Other major model configuration options are listed in Table 1. The National Centers for Environmental Prediction (NCEP) final (FNL) Operational Global Analysis data files were used to provide the initial and boundary meteorological conditions for our simulations. The initial and boundary chemistry conditions were provided by the output of the Model of Ozone and Related Chemical Tracers (MOZART-4; Emmons et al., 2010). Both anthropogenic and natural emissions were inputted into the model system. Anthropogenic emissions were derived from the Multi-resolution Emission Inventory for China (MEIC) database (http://www. meicmodel.org/, last access: 8 May 2018). MEIC contains both gaseous and aerosol species, including SO 2 , NO x , NH 3 , CO, VOCs (volatile organic compounds), BC, OC (organic carbon), PM 10 , and PM 2.5 . Biogenic emissions were calculated using the Model of Emission of Gas and Aerosols from Nature (MEGAN; Guenther et al., 2006).
Two parallel experiments were conducted to investigate our subject: (1) simulation with aerosol feedback considering both the direct and indirect radiation effects from all chemical species (Exp_WF), in which aerosol optical properties were calculated at each time step and then coupled with the radiative transfer model for both shortwave and long wave radiation (Iacono et al, 2008); (2) simulation with aerosol feedback, excluding BC (Exp_WFexBC), in which only the optical property of BC was subtracted when calculating the shortwave radiation and optical properties of all other aerosols were retained and calculated in the same manner as in experiment EXP_WF . The two experiments started at 00:00 LT on 1 October 2015 and ended at 00:00 LT on 26 October 2015. In order to reduce the influence of initial conditions, the first nine days were designated as the spin-up period.
We collected observational data on both meteorology (temperature, wind direction, and wind speed) and air pollutants (Ozone, NO 2 , and PM 2.5 ) in five cities (Hefei, Maanshan, Nanjing, Zhejiang, and Wuxi, Fig. 1) in October 2015 to evaluate the model performance. Temperature, wind direction, and wind speed, with a temporal resolution of 3 h, were obtained from Meteorological Information Comprehensive Analysis and Process System (MICAPS). These data were measured by the national surface observation network operated by the China Meteorological Administra- These air pollutants and three other pollutants (SO 2 , PM 10 , and CO) were measured by the national air quality monitoring network operated by the Ministry of Environmental Protection of China. The China National Environmental Monitoring Center is responsible for ensuring data quality. More information on the air pollution measurements is available in Wang et al. (2014). In addition, measurements of shortwave irradiance and BC concentrations in the northern suburb of Nanjing (near the ozone monitoring site) were taken. Shortwave irradiance was measured using a pyranometer (MS-802F of EKO instrument, Japan). The measurement accuracy is 1 W m −2 and the temporal resolution is 1 min. More information is available in the instrument manual (https://eko-eu.com/files/ PyranometerManual20160926V11.pdf, last access: 10 December 2017). Hourly concentrations of BC were measured by using an Aethalometer (model AE-33 of Magee Scientific, USA). The sampling time interval was set to 1 min and the inlet flow was set to 5 L min −1 . More information is available in Drinovec et al. (2015).

Model evaluation
To evaluate model performance, the measured and simulated variables from 10 to 26 October 2015 were compared ( Fig. 2a, b). In addition, shortwave radiation and BC concentration data observed in Nanjing were used to evaluate the corresponding model predictions.

Meteorology evaluation
Regarding meteorological factors (Fig. 2a), which highly impact transport, deposition, and transformation in the atmosphere (Li et al., 2008), good performance of meteorological conditions will guarantee accuracy in the air pollutant simulations. In this study, the observed temperature, wind direction, and wind speed were used for evaluating the meteorological model. In addition, the shortwave irradiance measured in Nanjing was also used to evaluate the simulation of shortwaves. Model performance statistics, including index of agreement (IOA), mean bias (MB), and the root mean square error (RMSE), are shown in Table 2, along with benchmarks derived from Emery et al. (2001). The IOA of wind direction (WD) was based on the calculation suggested by Kwok et al. (2010), which considers the nature of WD. The other abovementioned performance statistics were calculated following the approach of Lu et al. (1997). WD showed similar patterns between the simulation and observations, but the MB values of WD in Zhenjiang, Hefei, and Wuxi were beyond the benchmark. However, the high values of IOA showed good agreement, indicating that the simulation successfully captured the wind direction during the period. Wind speed (WS) showed a slight overestimation in Hefei and Wuxi with positive MBs of 0.72 and 0.87. Furthermore, the WS statistics are acceptable, especially IOAs and RMSEs, which are within the scope of the benchmarks. The negative MBs of temperature at 2 m above the surface (T2) indicate that the model predictions slightly underestimated T2. High values of IOAs in these cities indicate acceptable agreements between measurements and simulations. The shortwave radiation (SW) time series showed similar patterns between the simulation and measurements in Nanjing. The high IOAs and low biases (MB and RMSE) suggest satisfactory model performance on SW, which is very important for the determination of photochemistry and the BC feedback mechanism.

Air pollutant evaluation
Measured hourly concentrations of O 3 , NO 2 , and PM 2.5 in the five cities were collected to evaluate the model performance on air pollutants. In addition, hourly BC concentrations measured in Nanjing were also used for evaluating the BC simulation. Statistical metrics on air pollutants, which include IOA, mean normalized bias (MNB), and mean fractional bias (MFB), are shown in Table 3, along with benchmarks derived from EPA (2005,2007).
For ozone, MNBs in Nanjing and Wuxi were slightly beyond the benchmarks. However, the time series shows similar patterns between the simulation and measurements (Fig. 2b), reflected by the high values of IOA. As an important ozone precursor, a good model performance on NO 2 is necessary. More than 0.7 of the IOA values reflect the good agreement between measurements and simulations. The values of MNBs and MFBs are close to those of other studies in this region in October . Owing to its complex chemical composition, accurate simulation of PM 2.5 was difficult. The IOA values of PM 2.5 were lower than those of O 3 ; however, MFBs were within the benchmarks, indicating that the model predictions of PM 2.5 are acceptable. For the BC time series, similar patterns were found between the simula- (WD is wind direction; WS is wind speed; T2 is temperature at 2 m above the surface; and SW is shortwave). (b) Same as Fig. 2a but for chemical variables (O 3 is ozone; NO 2 is nitrogen dioxide; PM 2.5 is particulate matter with diameter less than or equal to 2.5 µm; and BC is black carbon). In general, based on the comparisons between measurements and simulations of both meteorology and air pollutants, the similar time series patterns and acceptable statistical results demonstrate the good simulation capability of the model system in capturing the meteorological and chemical features. In this case, the meteorological and chemical features in eastern China from 10 to 26 October 2015 were well reproduced, which is favorable for our analysis and discussions.

Observational and model results of BC and ozone
Analyzing relevant observation data, ozone and BC concentrations on sunny days were selected. Subsequently, ozone concentrations were distributed into two sets based on the daytime mean BC concentrations being higher (black) and lower (red) than the monthly mean concentration (∼ 3 µg m −3 ). When the BC concentration was higher (Fig. 3a), the increase in ozone concentrations occurred later than when BC was lower during the 10:00 to 14:00 period. At 12:00 LT, the difference between the two patterns of ozone concentrations reached its maximum value. Because of the limited measurements, the formation of changes in ozone and the relationship between the BC-BL interaction and surface ozone could not be validated. In this case, numerical simulations using the online coupled chemistry transport model WRF-Chem provided an effective method to analyze and discuss this subject.
With the application of WRF-Chem, Ding et al. (2016) reported the suppression of the BL development induced by the warming effects of BC. By analyzing our model outputs, the maximum changes in the mean BL height ( BLH_MAX) over Nanjing were calculated. BLH_MAX is defined as the maximum difference in the hourly mean boundary layer height (BLH) over Nanjing between Exp_WFexBC and Exp_WF during morning (Eq. 1): · · ·· · ·, BLH t=12:00 In Eq. (1), BLH is the mean boundary layer height over Nanjing. The simulated mean BC columns (from the surface to 2 km) over Nanjing, during the occurrence of BLH_MAX, were also calculated. The relationship between the two variables is shown in Fig. 3b. In a similar fashion to Ding's study, the positive values suggest that BL development is suppressed by the warming effect of BC. This effect is more significant at higher BC concentrations. For example, the BC column was more than 9 mg m −2 on 17 October and the BLH_MAX was more than 400 m. Owing to the significant impact of BC on the BL in this case, we took the model result of 17 October as an example to study the relationship between surface ozone, BC, and BL development. From 10:00 to 14:00 LT on 17 October, the average distribution (Fig. 3c) of BC concentrations was inhomogeneous due to the inhomogeneous distribution of BC emissions (Qin and Xie, 2012) and wind fields at ground level. The northern areas were controlled by southerly winds, whereas northeast winds blew over the southern areas. In the central areas (most parts of Jiangsu and Anhui), winds were weak, which is favorable for the accumulation of air pollutants. BC was therefore mainly concentrated in Nanjing and the surrounding areas. Conversely, surface ozone displayed low concentrations over this region, suggesting that BC and ozone at ground level exhibit opposite distributions over Nanjing.

The BC-BL interaction and its effects on photolysis rate and ozone precursors
Because of the solar radiation absorption of BC in the atmosphere, the incident shortwave radiation is attenuated. Consequently, other meteorological elements are affected. Figure 4a presents the vertical profiles of BC concentration (black solid line) and changes in related meteorological elements (Exp_WF−Exp_WFexBC) induced by BC at 10:00, 12:00, and 14:00 LT. Once emitted from surface (Qin and Xie, 2012), BC concentration in the BL is normally higher than that above the boundary layer where the so-called residual layer is positioned in the morning. However, the lower concentrations of BC aloft could also significantly influence the meteorological elements. The vertical profiles of the atmospheric attenuation of shortwave radiation (SW; blue solid line), which is defined as the gradual vertical loss in intensity of the incident solar radiation induced by BC, showed a small reduction above 1.2 km. Shortwave radiation also increased with as BC increased downwards. A greater reduction in shortwave radiation could be observed at a height of 1 km than at the lower adjacent level. This suggests that, when BC concentrations are similar in adjacent layers, solar radiation absorption of BC is more efficient at higher altitudes than that at lower altitudes (Ferrero et al., 2014;Ding et al., 2016). It should also be noted that when BC concentrations are sufficiently high, large amounts of BC absorb more shortwave radiation; hence, the attenuation of shortwave radiation displayed a second peak below the height of 600 m in this study. Furthermore, owing to the absorption of solar radiation induced by BC, the heating rate of shortwave increased correspondingly. The significant shortwave absorption in the residual layer led to a greater increase in heating rate (magenta solid line), which caused a consistent heating effect on air aloft (Fig. 4a 10:00 and 12:00 LT). In this case, it was favorable to form the temperature inversion which showed that the equivalent potential temperature (EPT; grey solid line) in Exp_WF was lower at the top of the BL (black dashed line) than that in the residual layer. The temperature inversion increased the atmospheric stability and suppressed and delayed the development of the BL (black dashed line) at 10:00 and 12:00 LT. Conversely, the EPT in Exp_WFexBC (orange solid line) did not exhibit this behavior and the BL (red dashed line) developed rapidly, reaching 1 km in altitude by 12:00 LT. At 14:00 LT, the BL in Exp_WFexBC was almost fully developed, whereas the BL in Exp_WF was still developing and rising close to the BL in Exp_WFexBC. As one of the most important ozone precursors, the photolysis of NO 2 directly contributes to the chemical formation of tropospheric ozone. Due to the attenuation of incident solar radiation induced by BC, the photolysis rate of NO 2 was reduced in the daytime ( J [NO 2 ]; brown solid line in Fig. 4b), which is consistent with the results from previous studies (Dickerson et al., 1997;Li et al., 2005). NO 2 in the atmosphere is primarily formed through chemical production or is directly emitted from surface sources (e.g., industry, transportation, and soil) and is usually confined in BL. As a result of the BL suppression caused by BC, more NO 2 is confined below the BL, as shown in Fig. 4b ( NO 2 ; green solid line). Similarly to NO 2 , other ozone precursors (e.g., NO, VOCs) also showed the same variations in vertical distribution (Supplement Fig. S1). The chemical production of tropospheric ozone is affected by both photolysis rate and the concentrations of precursors (Crutzen, 1973;Tie et al., 2009). When only considering the photolysis rate, the reduction of J [NO 2 ] will weaken the photochemistry and re-duce ozone concentration. When considering the mechanism of the BC-BL interactions at the same time, however, the concentrations of ozone precursors will change which also influences the ozone chemical production. The results may consequently be different from previous studies which only consider the effects of photolysis rate reduction. In addition, the suppressed BL could also affect physical processes in the BL, which will therefore influence surface ozone. Quantitative changes in ozone contributions from chemical and physical processes will be discussed in the following section.

Changes in surface ozone resulting from BC-BL interaction
Since variations in ozone concentration are directly caused by physical and chemical processes (Zhu et al., 2015), the fact that BC-BL interactions affect the ozone concentration can be discussed using process analysis. The following processes were considered in this work: (1) advection (ADV) caused by transport, which is highly related to wind and ozone concentration gradients from upwind areas to downwind areas; (2) vertical mixing (VMIX) caused by atmospheric turbulence and vertical gradients of ozone concentra-tion, which is closely related to the development of, and variations in, the BL (Zhang and Rao, 1999;Gao et al., 2017); and (3) chemistry (CHEM), which is the result of chemical calculation including ozone chemical production and chemical reduction. Convection process, i.e., the ozone contribution caused by convection movements, was negligible in this case, and is not mentioned further in this paper. Complete details on the process analysis of WRF-Chem are available in Zhang et al. (2014), Gao et al. (2016), and in the WRF-Chem user guide. As shown in Fig. 5a, the simulated surface ozone concentrations in Exp_WF and Exp_WFexBC were similar, below 20 ppb before 10:00 LT, and increased slowly. From 10:00 LT, the surface ozone concentrations in both cases started to increase rapidly. However, owing to the impacts of BC, ozone in Exp_WF increased more slowly than that in Exp_WFexBC which displayed ozone reduction. At 12:00 LT, the ozone reduction reached a maximum value of 16.4 ppb. From 12:00 to 14:00 LT, the increasing trend of ozone in Exp_WFexBC slowed which reduced the difference in ozone between the two cases. After 14:00 LT, there was no significant differences between the two ozone patterns. By analyzing the process contributions of ozone in Exp_WFexBC (Fig. 5b), VMIX contributed the most whereas ADV contributed the least to surface ozone from 10:00 to 14:00 LT. The result of the three processes, also denoted as the net contribution of ozone (NET), shows the hourly variation of the ozone concentration. In Fig. 5b, the significant positive contributions of VMIX and the consistence in the variations of VMIX and NET suggested that vertical mixing processes were mainly responsible for the increase in ozone in Exp_WFexBC. Because of BC, ozone contributions from chemical and physical processes changed at the surface; this information (Exp_WF−Exp_WFexBC) is presented in Fig. 5c. From 10:00 to 12:00 LT, changes in the chemical contribution (CHEM_DIF) increased at a rate of 3.1 ppb h −1 . However, changes in VMIX contribution (VIMX_DIF) decreased more significantly with a rate of −8.2 ppb h −1 . NET_DIF in this study indicated that the hourly change of ozone concentration is caused by the impacts of BC. Thus, the negative values of NET_DIF from 10:00 to 12:00 LT suggested that surface ozone reduced continuously. In the afternoon (13:00 to 14:00 LT), the CHEM_DIF was near zero which seemed to have little influence on the change of surface ozone. However, VMIX in Exp_WF was enhanced (VMIX_DIF is positive) which led to a NET_DIF change from negative to positive. Further, the positive NET_DIF indicated that the reduction in surface ozone was offset which was because the increasing trend of surface ozone in Exp_WFexBC slowed down during this period.
Among all the changes in the contributions from processes, VMIX_DIF accounted for more than 50 % of NET_DIF, indicating that change in the ozone contribution from vertical mixing was the key factor impacting surface ozone in this study. In addition, J [NO 2 ] was reduced by BC but changes in CHEM at surface increased before noon. A reduction in the photolysis rate should reduce the chemical production of ozone; however, changes in ozone precursors will also influence the chemical production of ozone. As Fig. S2 shows, due to the suppression of the BL, the concentrations of both NO x and VOCs increased at surface. The ratios of VOCs / NO x also increased which suggested that VOCs increased more significantly than NO x did. In addition, the small change in the ratio of HCHO / NO y indicated that ozone still formed under VOC-limited conditions. In this case, the more significant increase of VOCs was favorable for ozone chemical production. Thus, it could be concluded that the enhancement of the contributions from chemistry was highly related to the increase of ozone precursors.
The development of the BL directly affects vertical mixing process and the distribution of ozone in the BL (Zhang and Rao, 1999). Herein, the distributions of VMIX contributions with and without the effects of BC are shown in Fig. 6a, b. The BLH with (black line) and without (red line) the impacts of BC are also presented in Fig. 6. In general, ozone concentration increases with height before sunrise and in the early morning over polluted regions (Wang et al., 2015). When the BL developed in the morning, ozone was exchanged by turbulence in the vertical direction resulting in the entrainment of the abundant ozone aloft down to the surface and the increase of surface ozone concentrations. Therefore, VMIXs in both cases showed positive values in lower layer and negative values in upper layer. Differences in VMIX (VMIX_DIF) between Exp_WF and Exp_WFexBC showed the changes in ozone contribution from vertical mixing caused by BC and the vertical distribution of VMIX_DIF is shown in Fig. 6c. With the effects of BC, the suppressed BL weakened turbulence which led to reduced vertical exchange of ozone from 10:00 to 12:00 LT. Thus, the VMIX_DIF showed negative values in the lower layer and positive values in the upper layer. After 12:00 LT, the BL in Exp_WFexBC was fully developed and VMIX became weak. By contrast, the BL in Exp_WF was still developing during this period. Under these conditions, VMIX was not weakened or even larger than that in Exp_WFexBC. Thus, VMIX_DIF showed positive values in lower layers and negative values in upper layers.
Similarly to the VMIX_DIF, changes in ozone contribution from chemistry (CHEM_DIF) are presented in Fig. 6d. The chemical production of ozone is related to both photolysis rate and ozone precursors. As discussed above, the BC-BL interaction leads to reduced photolysis rates but also to more accumulation of ozone precursors at surface. The reduction of photolysis rate will surely decrease photochemical rate; however, the increase of ozone precursors at the surface may offset the influence of the reduction in photolysis rate and enhance the chemical production of ozone. Thus, enhanced chemical contribution could be seen near the surface (with a positive CHEM_DIF of 3 ppb h −1 ) from 10:00 to 12:00 LT. In the residual layer, contrary to the situation near  the surface, less ozone precursors being transported up owing to the suppression of BL ( Fig. 4b and S1) and the reduction of photolysis rate (Fig. 4b) led to a decreased chemical contribution (with a negative CHEM_DIF between −9.4 and −2.1 ppb h −1 ) from 11:00 to 14:00 LT. Owing to the effects of BC, the contribution of ADV decreased slightly within the BL before noon (Fig. 6e). In Fig. 6f, NET_DIF, the result of VMIX_DIF, CHEM_DIF, and ADV_DIF, represents the hourly change of ozone concentrations caused by the impacts of BC. Comparing Fig. 6c-f, the similar distributions and the significant ratio of VMIX_DIF to NET_DIF suggest that the change in the ozone contributions from vertical mixing induced by BC had the greatest influence on surface ozone.
Due to BC-BL interactions, VMIX in Exp_WF is weakened before noon but enhanced in the afternoon (Fig. 7a). The contribution of VMIX depends on the vertical gradients (the difference of ozone concentrations between every two adjacent vertical layers) of ozone and the turbulent exchange coefficient. In Fig. 7b, vertical gradients of ozone (solid lines) showed similar profiles in both experiments dur- ing the early stage of BL development before noon. In the afternoon, the ozone gradient (red solid line in Fig. 7c) decreased significantly without the effects of BC in the mature stage of BL development, whereas the ozone gradient remained large around the top of the BL owing to the effects of BC (black line in Fig. 7c). The effects of BC can lead to a significantly smaller turbulence exchange coefficient (black line with square in Fig. 7b), lowering the entrainment of ozone from higher altitudes to the surface, which indicated that VMIX in Exp_WF was smaller before noon. In the afternoon, the faster developing BL without the effects of BC formed a uniformly mixed ozone profile, with very small vertical gradients. As a result, the contribution of VMIX in Exp_WFexBC was smaller, although the turbulence exchange coefficient was still large. In contrast, with the impacts of BC, larger ozone gradients around the top of the BL resulted in stronger entrainment of ozone down to the surface which led to the ozone contribution of VMIX in Exp_WF becoming larger at surface.

Conclusions
In this study, measured data of BC and ozone in Nanjing, an important industrial and economic center in the YRD region of China, showed that ozone concentration increased slowly from 10:00 to 14:00 LT when the BC concentration was relevantly high. With the application of WRF-Chem model, we studied the impacts of BC-BL interactions on surface ozone in Nanjing. Acceptable agreement was achieved between observed and simulated results for both meteorological and chemical variables, suggesting that the WRF-Chem model has the capability to accurately reproduce meteorological factors and air pollutants in this study.
The model results showed that, when high concentrations of BC were confined over Nanjing, the incident solar ra-diation was absorbed. The absorbed radiation heated the air above the BL, suppressed and delayed the development of the BL. Because of the BL suppression, it caused the changes in the ozone contributions from chemical and physical processes (photochemistry, vertical mixing and advection). In particular, BL suppression led to more ozone precursors (e.g., NO x and VOCs) being confined below the BL. Although the photolysis rate was reduced as a result of the solar radiation absorption induced by BC, more ozone precursors enhanced the contributions from chemistry before noon. However, as a more significant factor, the suppressed BL weakened the turbulence and entrained less ozone aloft down to the surface which reduced the contributions from vertical mixing process. As a result, surface ozone decreased before noon, with the maximum reduction reaching 16.4 ppb at 12:00 LT. In the afternoon, the changes in the chemical processes were small which inconspicuously influenced surface ozone. However, the vertical mixing process still significantly influenced the surface ozone. Due to the delay in the BL development, the vertical ozone structure showed obvious vertical gradients around the top of the BL, and the high concentration of ozone aloft could be entrained down to the surface to offset the surface ozone reduction before noon.
Comparing the changes in the contributions from the chemical and physical processes induced by the effects of BC, the change in vertical mixing process which was caused by BL suppression played a more important role in influencing surface ozone. Our results highlight the great impacts of BC-BL interactions on surface ozone via its influence on ozone contribution from physical processes, which suggests that more attention should be paid to the mechanism of aerosol-BL interactions when controlling ozone pollution in China.
Data availability. All observations and model outputs used in this study are available. Readers can access the data directly or by contacting Bin Zhu via binzhu@nuist.edu.cn.