• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Analysis of the Effect of Optical Properties of Black Carbon on Ozone in an Urban Environment at the Yangtze River Delta, China

    2021-06-22 07:29:24JunlinANHuanLVMinXUEZefengZHANGBoHUJunxiuWANGandBinZHU
    Advances in Atmospheric Sciences 2021年7期

    Junlin AN, Huan LV, Min XUE, Zefeng ZHANG, Bo HU, Junxiu WANG, and Bin ZHU

    1Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD), Nanjing University of Information Science and Technology, Nanjing 210044, China

    2Chinese Academy of Meteorological Sciences, Beijing 100081, China

    3Institute of Atmospheric Physics Chinese Academy of Sciences, Beijing 100029, China

    4Inner Mongolia Meteorological Bureau, Hohhot 010051, China

    ABSTRACT Black carbon (BC) reduces the photolysis coefficient by absorbing solar radiation, thereby affecting the concentration of ozone (O3) near the ground. The influence of BC on O3 has thus received much attention. In this study, Mie scattering and the tropospheric Ultraviolet and Visible radiation model are used to analyze the effect of BC optical properties on radiation. Combined with data of O3 precursors in Nanjing in 2014, an EKMA curve is drawn, and the variations in O3 concentration are further investigated using a zero-dimensional box mechanism model (NCAR MM). When O3 precursors are unchanged, radiation and O3 show a highly similar tendency in response to changing BC optical properties (R=0.997).With the increase of modal radius, the attenuation of fresh BC to radiation and O3 first trends upward before decreasing. In the mixing process, the attenuation of BC to radiation and O3 presents an upward tendency with the increase of relative humidity but decreases rapidly before increasing slowly with increasing thickness of coating. In addition, mass concentration is another major factor. When the BC to PM2.5 ratio increases to 5% in Nanjing, the radiation decreases by approximately 0.13%–3.71% while O3 decreases by approximately 8.13%–13.11%. The radiative effect of BC not only reduces O3 concentration but also changes the EKMA curve. Compared with the NOx control area, radiation has a significant influence on the VOCs control area. When aerosol optical depth (AOD) increases by 17.15%, the NOx to VOCs ratio decreases by 8.27%, and part of the original NOx control area is transferred to the VOCs control area.

    Key words: black carbon, ozone, radiation, optical properties, EKMA curve

    1. Introduction

    Tropospheric ozone (O) is a typical secondary pollutant in the atmosphere and poses a serious threat to human health, vegetation, and the atmospheric environment(Crutzen et al., 1973; Feng et al., 2015; Goodman et al.,2015). Urban air quality environmental monitoring data show that as the Oconcentration near the ground has increased year by year in China, Ohas gradually become the main pollutant affecting air quality in China (Li et al.,2019). Meteorological factors and synoptic patterns have great impacts on the variation of O(Dong et al., 2020).Also, the major contribution to Oconcentrations near the ground is well recognized to come from aerosol particles,which is clear from this study.

    Black carbon (BC) accounts for a small fraction of aerosols, but it exerts a large radiative effect on the energy balance. BC is produced by incomplete combustion of carbonaceous fuels (Bond et al., 2004). Strict emissions control and reduction measures implemented in China have made both aerosol (Zhang et al., 2019a) and BC mass concentrations (Zhang et al., 2019b) decrease sharply. By absorbing solar radiation, BC can not only change the global and regional radiative balance but also reduce the photolysis coefficient, thereby reducing the concentration of Onear the ground, which makes the influence of BC on Oalways of concern (Jacobson, 2001; Ramanathan et al., 2008; Yang et al., 2016).

    Chen et al. (2019) found a negative correlation between BC and Othrough observations in the Nanjing area. Dickerson et al. (1997) proposed that BC could reduce the surface Oconcentration by reducing the photolysis rate. Jacobson(1998) and Castro et al. (2001) found that aerosols containing BC significantly reduced the photolysis coefficient by absorbing solar radiation, and the near-surface Oconcentration decreased by 5%–8%. Li et al. (2005) simulated a period of heavy air pollution in Houston through a regional chemical transport model and found that BC reduced the photolysis coefficient in the boundary layer by 10%–30% and Oconcentration by 5%–20%. Gao et al. (2018) found that absorption of radiation by BC at the top of the boundary layer inhibits the development of the boundary layer in the daytime, and a large number of Oprecursors are limited to the boundary layer, thus promoting the photochemical generation of O. This would counteract the effect of BC lowering the photolysis coefficient to inhibit Oformation. However,the experimental results from studies of factors such as BC scale distribution and aging process are still highly uncertain (Huang et al., 2015). The optical properties of BC and its radiative effect are not fully understood, and the mechanisms of the influence of BC on Oneeds to be further studied.

    The optical properties of BC, which are closely related to its particle size distribution and mixing state, are the basis for simulating its radiative effect (Bond et al., 2013) and directly affect the simulated results for O(Saathoff et al.,2003; Chung et al., 2005; Ma et al., 2011). Based on a number of observations of BC particles, both BC number concentration and mass concentration are characterized by a logarithmic normal distribution (LND) (McMeeking et al., 2010;Huang et al., 2012; Wang et al., 2014). The peak particle size of approximately 120–210 nm (Kondo et al., 2006;Gong et al., 2015) is related to the environment and the emission source (Park and Lee, 2015). After being discharged into the atmosphere, BC tends to absorb other aerosols and mix with them (Oshima et al., 2009). Mixing mode, shell chemical composition, thickness, and environmental humidity are important factors that affect the mixing state of BC.Jacobson (2001) found that BC could significantly enhance radiation forcing when it was mixed with other groups, and the rate of increase was related to the hygroscopic growth effect of the aerosol. Zhang et al. (2008) found that, compared with fresh particles, the absorption efficiency of internal mixed particles was nearly two times higher and the scattering efficiency was approximately ten times higher in a high-humidity environment. Liu et al. (2015) proposed that the enhanced optical absorption effect of BC depends largely on the thickness of the shell layer, which is related to the emission source. Therefore, it is necessary to conduct sensitivity tests on BC optical parameters when studying the influence of BC on O.

    In this paper, to analyze the sensitivity of radiation and Oto BC optical properties, the mass concentration ratio of BC in particulate matters (PM) and the spectral distribution and the mixing state of BC were changed. As detailed in the following sections, Mie scattering theory and the program of coated spheres according to Bohren and Huffman(BHCOAT) scattering program were used to calculate optical parameters of BC, the tropospheric Ultraviolet and Visible radiation (TUV) model was used to calculate the value of radiation and photolysis coefficients, and a zerodimensional box mechanism model (NCAR MM) was used to calculate Oconcentrations under different optical parameters of BC. Finally, the variation in the Empirical Kinetics Modeling Approach (EKMA) curves after radiation attenuation was analyzed by using the observed data for nitrogen oxides (NO) and volatile organic compounds (VOCs) in Nanjing.

    2. Experimental methods and scheme design

    2.1. Experimental methods and data

    In this study, a box model with no source emission or advection and diffusion processes and only photochemical processes is used for simulation and calculation. The model is composed of three modules: aerosol optical parameter, radiative transfer, and the chemical mechanism related to phosgene (Fig. 1). The optical properties of BC are mainly calculated by Mie scattering theory. The radiation parameters and Oconcentration are respectively solved by the TUV model and the NCAR MM model. The optical parameters of fresh BC and PMare solved by the Lorenz-Mie Scattering code developed by Michael Mishchenko. Mie scattering assumes that the particles are spherical. Through the input of complex refractive index and spectral distribution, the following optical parameters of the aerosol particle swarm are obtained: single scattering albedo (

    ω

    ), asymmetric factor(g), extinction efficiency (Q), and volume extinction coefficient (k). Refractive index and density of PM(Dubovik et al., 2002; Pitz et al., 2003) and those of BC (Hess et al.,1998; Gong et al., 2015) are shown in Table 1. A logarithmic normal distribution (LND) is selected, and the analytical expression is shown in Eq. (1):

    Fig. 1. Structure diagram of box model.

    Table 1. Aerosol parameter settings.

    where ris the modal radius, σ is the width of the spectral distribution, and N is the particle number concentration.When considering the mixing process and the influence of relative humidity (RH), the k-EC-Mie model is adopted for the optical properties of BC (Zhang et al., 2017). The k-EC-Mie model is based on the core-shell model and assumes that 1)the physical and chemical properties of particles with the same scale are the same; 2) under dry conditions, the particle consists only of BC (core), located in the center of the particle, and an envelope (shell); and 3) regardless of the hygroscopicity of BC, the coating is evenly mixed with water before the hygroscopic growth process. Finally, the layered spherical Mie scattering method (BHCOAT scattering program) proposed by Bohren and Huffman (2008) is used for calculations. The core-shell volume ratio of each particle is set as constant when calculating the optical properties of particle swarm.

    The TUV model is a radiation transmission model developed and continuously improved by the National Center for Atmospheric Research (NCAR) and embedded in the NCAR MM model (Madronich and Flocke, 1999). It can calculate various radiation quantities such as solar incident radiation intensity and photolysis rate coefficients of 73 photochemical reactions with wavelengths in the range of 121–1000 nm (Stamnes et al., 1988). The TUV model is widely used in the NCAR MM model, WRF-Chem model,and other regional and global atmospheric chemistry models (Shao et al., 2016; Ryu et al., 2018). Palancar et al.(2013) present a comparison between measured actinic fluxes and those simulated with the TUV model and found that the simulated result is larger and the difference between the two is about 10%–40%. In addition to input parameters such as w and g, aerosol optical depth (AOD) with a wavelength of 0.55 μm near the ground should also be calculated using Eqs. (2) and (3) (Hess et al., 1998):

    The NCAR MM model is used to calculate the chemical mechanism related to phosgene. The model includes nearly 5000 gas phase chemical reactions between more than 2000 substances. Assuming no other emission sources,no dilution, and no recombination process, it can calculate the time-dependent chemical evolution of air with known initial components (Madronich and Calvert, 1990). Users define input parameters such as initial gas component and volume fraction and then calculate the Oconcentration by combining with photolysis coefficient output from the TUV model. Lingaswamy et al. (2017) used the TUV model to calculate the photolysis rate coefficients [

    J

    (O),

    J

    (NO)], used the NCAR MM model to simulate diurnal variation of CO,and reported their results.

    The volume fraction of gas and meteorological data required by the NCAR MM model are provided as hourly means from March 2014 to February 2015. The observation site is located in the Key Laboratory of Aerosol and Cloud Precipitation in Nanjing University of Information Technology (32°12′N, 118°42′E; altitude 62 m). O, NO,NO, and CO were observed using an air pollution environmental monitoring analyzer produced by Thermo Fisher Scientific in the United States, including an ultraviolet luminescence Oanalyzer (Model 49i), infrared absorption CO analyzer (Model 48i), and chemiluminescence NO-NO-NOanalyzer (Model 42i). Detailed instrument parameters and calibration methods can be found in Hansen et al. (1984). VOCs were observed using a GC5000 automatic online gas chromatography-flame ionization detector produced by AMA Company in Germany, and detailed instrument parameters and calibration methods can be found in An et al. (2014).

    2.2. Experimental scheme design

    The sensitivity experiment consists of a control experiment and four variant schemes. In the experiments, PMand BC are fully independent, and the sum of their mass concentration is always 75 μg m. The mixing mode between the two is external mixing. The external mixing method is actually independent scattering between particles. The aerosol optical parameters of PMand BC need to be calculated independently first, and then the optical properties of the mixed aerosol can be calculated according to the weight of the mass concentration.

    In the control test, PMis regarded independently without the influence of BC. Then, the net radiation irradiance (579.85 W m) and the Oconcentration (46.76 ppb)near the ground can be calculated, where the concentration of PMis determined by the Chinese ambient air PMstandard (75 μg m). In sensitivity tests, mass concentration ratio of BC to PM, the spectral distribution, and mixing state are respectively changed (Table 2). Based on the observed data in Nanjing from 2015 to 2018, Tan et al.(2020) found that the annual mean ratio of BC to PMwas about 4.92%, and the monthly mean was about 1.5% to 10.1%. Therefore, in this study, the ratio of BC to PMis set at 3%–9% when discussing the influence of the mass concentration of BC on O, and the BC to PMratio is set at 5% when discussing the influence of the spectrum distribution and mixing process of BC on O. REO (%) is defined as the relative change between a sensitivity test and control test [ Eq.(4)]:

    If REO is negative, then BC has an attenuation effect on net radiation irradiance and O. The greater the absolute value of REOand REO, the greater the attenuation of BC to radiation and O. The above Oconcentration refers to the average of a continuous eight hours of an Oconcentration at the largest daily values, and the net irradiance is the average value in the corresponding period.

    2.3. Model validation

    To verify the applicability of the above PMspectral distribution in the Nanjing area and the accuracy of the NCAR MM model, the observed concentration of Onear the ground is compared with the simulation results. In this study, the measured data during sunny weather with a wind speed less than 4 m sin 2014 are selected as the input to the NCAR MM model, resulting in six days being selected.By using these criteria, the simulation error caused by the transmission effect, wet deposition, and the influence of clouds will be significantly reduced (Shao et al., 2016).According to the above PMspectral distribution and the daily mean value of PMobservation data, the daily ground-level AOD can be calculated. The calculated AOD and other observed data are input into the NCAR MM model for simulation. The 6-hour average value of volume fraction of Oprecursors is used as the limit of each simulation. According to the six selected samples, six simulations are performed. The simulation results show the diurnal variation of Oand are highly correlated with the observations(R = 0.87); all reported correlation coefficients in this paper pass the significance test (P < 0.05). As shown in Fig. 2, compared with the observation data, the simulation results for tendencies in Ovariation and concentration perform well during the daytime (0600–1800 LST, LST = UTC + 8), but simulated Oconcentration is higher at noon, which may be related to real-time cloud changes (Zhang et al., 2015). In general, the observed tendency is consistent with the simulated tendency. The slope of the fitting curve is 1.085, indicating that the simulated values are relatively close to the actual values. This shows that the NCAR MM model can accurately simulate the variation tendencies and concentrations of O,and the spectral distribution setting of PMis reasonable.

    Table 2. Schematic setting for the sensitivity analysis.

    3. Results

    3.1. The influence of BC on radiation

    Fig. 2. Comparison of observed ozone and modeled ozone in the northern suburbs of Nanjing for the selected six days in 2014.

    Fig. 3. Effects of BC mass concentration on (a) net radiation and on (b) concentration of O3. REONRI and REOO3 represent the relative changes of net radiation and O3 concentration between sensitivity test and control test, respectively. “BC”,“BC-H2SO4”, and “BC-H2SO4” stand for fresh BC particles,old BC that fresh BC particles have coated with sulfate under dry conditions, and old BC that fresh BC particles have coated with sulfate under an environmental RH of 50%, respectively.

    Table 3. Optical properties of particles in three states.

    To explore the effect of fresh BC spectral distribution on radiation, rand width (

    σ

    ) spectral distribution are changed. The larger the ror

    σ

    is, the larger the particle size of the particle population. The net radiation attenuation of fresh BC is approximately 1.75%–3.71%. As shown in Fig. 4b, REOfirst decreases and then increases with increasing r. When the value of

    σ

    is 1.8, 2.0, and 2.2,the ris approximately 0.041 μm, 0.026 μm, and 0.017 μm, respectively (i.e., a shift towards small particles when REOis the smallest). The corresponding effective radius(reff) is 0.097 μm, 0.088 μm, and 0.079 μm, respectively(Table 4). That is, when the reff of particle swarm is approximately 0.1 μm, BC has the strongest attenuation of radiation, which is close to the results of Zhang et al. (2017). It is found that the ratio of σreffπ to wavelength is approximately 1 in the three cases, and the above optical parameter is calculated when wavelength is 0.55 μm. The attenuation of BC is strongest when the ratio of σ reffπ is close to wavelength.

    Fig. 4. Effects of spectral distribution on (a) changes in optical parameters when σ is 2 and λ is 0.55 μm and on (b) changes in radiation and concentration of O3.

    Table 4. The values of optical parameters with maximum k ′ext .

    Two enveloping substances are considered next,namely, sulfuric acid (HSO) and organic carbon (OC).Table 5 shows the relevant parameters of the two enveloping substances, in which k is the moisture absorption parameter. The real and imaginary parts of the complex refraction index of OC are both greater than that of HSO, which makes the

    k

    of BC enclosing OC larger. Therefore, the attenuation of BC to radiation is stronger when the wrapper is OC. The hysteretic parameter of HSOis 1.13, which is much higher than OC. In other words, when the coating is HSO, RH has a more significant influence on BC (Fig. 5).Taking the wrapper of BC being HSOas an example,Fig. 7 describes the change of particle density and other parameters with changes in RH under different particle sizes.With the increase of RH, the density and refractive index of the shell show a monotonically decreasing tendency (Fig. 7b).However, due to the increase of shell hygroscopicity and particle size, the particle volume and extinction coefficient show a monotonically increasing tendency, which is consistent with other studies (Zhang et al., 2017). When the thickness of coating is small (0.012 μm < r< 0.015 μm), the contour lines are dense and nearly straight. In this case, the hygroscopic growth effect of the shell is weak, kand mass are almost unchanged (Fig. 7a), and the effect of RH on

    k

    is very small. When RH increases from 0 to 95%, the attenuation value of net radiation only changes by 0.36%, and the thickness of coating is the main influencing factor. When the thickness of coating is close to the low center value(0.02 μm < r< 0.025 μm), the contour thinning and curvature increase, the influence of RH on BC begins to increase, and the attenuation value of BC to radiation increases with the increase in RH. When the thickness of coating is large (r> 0.025 μm), the hygroscopic growth effect of the coating is more obvious, and mass and kincrease significantly (Fig. 7a). In this case, BC reduces the net radiation by 1.25%, with this attenuation effect mainly being caused by the change in RH.

    Fig. 5. Effects of mixing state on (a) changes in net radiation when the coating is H2SO4, (b) changes in O3 when the coating is H2SO4, (c) changes in net radiation when the coating is OC, and (d) changes in O3 when the coating is OC.The horizontal coordinate is rmod,N in a dry state, which is proportional to the thickness of the coating, while the colored contours represent the values of REONRI and REOO3.

    Fig. 6. Changes in particle mass and kext when RH is 20%(top) and the coating is H2SO4 (bottom).

    3.2. The influence of BC on ozone

    When using the NCAR MM model to calculate the change of O, the hourly concentration value of Oprecursors is taken as the limiting condition of the model and as input to the model. In the model simulation, the data of Oprecursors from the days with no precipitation and wind speed less than 4 m sin the northern suburbs of Nanjing in 2014 are used. Figure 8 describes the diurnal change process of NO, VOCs, and CO in the northern suburbs of Nanjing for the selected days. The diurnal variations of the three factors are consistent: the concentration reaches the highest value at 0800 LST and the lowest at 1400 LST.

    Table 5. Parameter settings for wrappers.

    In addition to the influence of Oprecursors, radiation is the most critical factor for O(Sillman, 1995). In the calculations, when conditions for Oprecursors are consistent,the net irradiance is significantly correlated with O(R=0.997), so radiation and Oin all schemes show similar variation tendencies (Figs. 3–5) (the impact of Oconcentration changing with BC radiative effect is described in section 3.1). The higher the mass concentration of BC, the greater the attenuation to Owill be (Fig. 3b). With the increase of r, the decay of fresh BC to Ofirst increases and then decreases, and Oconcentration decreases by approximately 10.13%–13.11%. When BC undergoes a mixing process, Oconcentration is reduced by approximately 8.13%–11.94%. The attenuation value of Ois consistent with that of radiation with respect to changing with the thickness and relative humidity of the coating.From what has been discussed above, when the BC to PMratio is 5% in Nanjing, the change in the attenuation of BC to radiation is approximately 0.98%–1.71%, and the change in attenuation of BC to Ocan reach approximately 1.63%–2.57% due to different spectral distributions of BC;changes in the attenuation of BC to radiation and Oare approximately 0.94%–1.97% and 1.54%–3.06%, respect-ively, due to different mixing processes of BC.

    Fig. 7. The change process of optical parameters with RH when the coating is H2SO4.

    Fig. 8. The diurnal variation process of NOx, VOCs, and CO in Nanjing in 2014.

    After the maximum point,

    k

    shows a decreasing tendency and the irradiance increases with it. At this time, w and g, which show a slight decreasing tendency for the mixed aerosol, have the effect of reducing the irradiance, and that causes the drastic fluctuations of REOand REO.To further analyze the influence of radiation on the pathway of Ogeneration after the attenuation effect of BC on radiation, the NCAR MM model is used to simulate the concentration of O(

    φ

    ) in the northern suburbs of Nanjing in 2014, and an EKMA curve for the Nanjing area is established (Fig. 9). The volume fraction of Oprecursors, such as NO, and the proportion of alkanes, olefins, and aromatic hydrocarbons in VOCs, etc., are all based on average values of their observational data from 0700 to 0900 LST(Table 6). Other meteorological elements, such as temperature, are based on hourly mean values of the annual observation data in 2014. In NOcontrol areas, OH mainly reacts with VOCs, which have little influence on O. Ois more sensitive to NO, so reducing NOis beneficial for hindering the generation of O. At VOCs > 200 ppbc and NO<15 ppb (i.e., when the VOCs to NOratio is large), the radiation has little impact on O, and an approximately 0.95 ppb reduction in NOachieves the

    φ

    in the control experiment. However, when the volume fractions of NOand VOCs are close to the ridge line, the effect of radiation is greater. When

    φ

    remains the same, the difference in NOconcentration corresponding to the two experiments is as high as 6.00 ppb. In the VOCs control area, the reaction of OH and NOis dominant, and Ois sensitive to VOCs. In this case, reducing NOwill promote the formation of O.In this region, radiation has a significant effect on

    φ

    . At VOCs < 150 ppbc, NOneeds to increase by approximately 2.14 ppb for

    φ

    to be consistent with its value when the radiation dose decreases; at VOCs > 250 ppbc, NOneeds to increase by approximately 9.47 ppb. In addition, radiation can also cause the ridge line of the EKMA curve to shift.When the AOD increases from 0.88 to 1.08, the slope of ridge line (the NOto VOCs ratio) decreases from 0.090 to 0.083. The original part of the region located in the NOcontrol area is transferred to the VOCs control area under the influence of radiation.

    Fig. 9. Influence of radiation on diurnal maximum volume fraction curve of O3. The control test (AOD=0.88, PM2.5) is shown as the red curve while the sensitivity test (AOD=1.08,9%BC+91%PM2.5) is shown as the black curve. The dotted line is the ridge line of the EKMA curve, which is formed by fitting the turning point of each O3 contour line and the boundary between the VOCs control area and the NOx control area.

    Table 6. Initial values of ozone precursors in NCAR MM.

    4. Summary

    BC can reduce the photolysis coefficient by absorbing solar radiation, thus reducing the near-surface Oconcentration. In this paper, the attenuation degree of BC to radiation is discussed in response to changing the optical properties of BC, including the mass concentration ratio of BC to PM, the spectrum distribution of BC, the thickness of coating, and RH in the internal mixing process. Combined with observed data for Oprecursors in Nanjing in 2014, an EKMA curve was calculated to examine the change in Oconcentration near the ground and pathways for Ogeneration after radiation attenuation caused by BC.

    For BC, the attenuation generated by radiation is closely related to its mass, and the greater the mass, the greater the influence on radiation. When the mass concentration of the BC to PMratio is 5% in Nanjing, the attenuation of BC particles to net irradiance is approximately 0.1%–3.70%, and the range of this attenuation is approximately 0.98%–1.71% and 0.94%–1.97%, respectively, due to different spectral distribution and mixing processes of BC.With the increase of r, the radiation attenuation of fresh BC tends to increase first and then decrease. With the increase of

    σ

    , the radiation attenuation of fresh BC tends to shift to the direction of small particles. The attenuation of BC is strongest when the product of spectral width, effective radius, and π is similar to the wavelength. After BC goes through a mixing process, the attenuation value of radiation presents an upward tendency with an increase in RH,but it first decreases rapidly and then increases slowly with the increasing thickness of coating. When r< 0.015μm, the effect of RH is very small and the thickness of coating is the main factor affecting the attenuation capacity of BC. When r> 0.025 μm, the moisture absorption of coating increases significantly, with RH as the main influencing factor.It is found that when Oprecursors are unchanged, net irradiance is significantly correlated with O(R=0.997), so in the process of changing with changes in BC optical parameters, radiation and Oshow highly similar change tendencies. When the ratio of BC to PMincreases to 5% in Nanjing, Ocan be reduced by approximately 8.13%–13.11%,and the range of this attenuation is approximately 1.63%–2.57% and 1.54%–3.06%, respectively, due to different spectral distribution and mixing processes of BC. The attenuation effect of BC on radiation not only significantly reduces Oconcentration but also changes the EKMA curve. In the VOCs control area, the effect of radiation on

    φ

    is relatively significant. When radiation changes while

    φ

    stays the same, the difference in NOconcentration is approximately 2.14–9.47 ppb. In the NOcontrol area, only when the volume fraction of NOand VOCs is close to the ridge line, the effect of radiation is obvious, and the difference in NOconcentration is approximately 6.00 ppb, at the most, when

    φ

    stays the same. In addition, when the AOD increases from 0.88 to 1.08, the slope of the ridge line (the NOto VOCs ratio) decreases from 0.090 to 0.083 (i.e.,when the AOD changes by 17.15%, the NOto VOCs ratio will decrease by 8.27%), and part of the original NOcontrol area will be transferred to the VOCs control area.

    Acknowledgements

    . This work was supported by grants from the National Key Research and Development Program of China (Grant No. 2017YFC0210003), the National Natural Science Foundation of China (Grant No. 42075177), and the Qing Lan Project.

    亚洲丝袜综合中文字幕| 久久国产亚洲av麻豆专区| 22中文网久久字幕| 精品久久久噜噜| 国产精品久久久久久精品古装| 日日啪夜夜撸| 欧美日本视频| 狠狠精品人妻久久久久久综合| 色哟哟·www| 亚洲精品视频女| 三级经典国产精品| 国产午夜精品一二区理论片| 国产一区有黄有色的免费视频| 国产精品99久久99久久久不卡 | 国产精品一区二区性色av| 久久鲁丝午夜福利片| 赤兔流量卡办理| 观看免费一级毛片| av在线app专区| 黑人猛操日本美女一级片| 国产免费福利视频在线观看| 亚洲精华国产精华液的使用体验| 国产片特级美女逼逼视频| 九九久久精品国产亚洲av麻豆| 日韩,欧美,国产一区二区三区| av线在线观看网站| 亚洲精品国产av成人精品| 国产一区二区在线观看日韩| 亚洲欧美精品专区久久| 久久亚洲国产成人精品v| av在线老鸭窝| 中文在线观看免费www的网站| 国产无遮挡羞羞视频在线观看| 少妇人妻精品综合一区二区| 国产亚洲精品久久久com| 欧美日韩国产mv在线观看视频 | 国产又色又爽无遮挡免| 亚洲精品456在线播放app| 一级毛片我不卡| 热99国产精品久久久久久7| 欧美日韩视频高清一区二区三区二| 丰满乱子伦码专区| 久久青草综合色| 最近手机中文字幕大全| 狂野欧美激情性bbbbbb| 纯流量卡能插随身wifi吗| 插阴视频在线观看视频| 女人十人毛片免费观看3o分钟| 亚洲精品自拍成人| 亚洲欧美日韩东京热| 亚洲三级黄色毛片| 一级av片app| 热99国产精品久久久久久7| 亚洲欧美日韩卡通动漫| 久久久精品94久久精品| 亚洲怡红院男人天堂| 日韩av在线免费看完整版不卡| 久久久久久久久久人人人人人人| 国产一区二区三区综合在线观看 | 国产黄色视频一区二区在线观看| 校园人妻丝袜中文字幕| 街头女战士在线观看网站| 国产白丝娇喘喷水9色精品| 国产一区二区三区综合在线观看 | 观看免费一级毛片| 直男gayav资源| 欧美成人a在线观看| 有码 亚洲区| 搡女人真爽免费视频火全软件| 亚洲真实伦在线观看| 亚洲av二区三区四区| 欧美变态另类bdsm刘玥| 80岁老熟妇乱子伦牲交| 国产伦精品一区二区三区四那| 亚洲第一av免费看| 一个人免费看片子| 男人添女人高潮全过程视频| 直男gayav资源| 精品亚洲成a人片在线观看 | 婷婷色综合大香蕉| 亚洲精品一二三| 国产乱来视频区| 亚洲婷婷狠狠爱综合网| 亚洲精品色激情综合| 午夜福利网站1000一区二区三区| 国产伦精品一区二区三区四那| 国产免费视频播放在线视频| 一本色道久久久久久精品综合| 男女啪啪激烈高潮av片| 九草在线视频观看| 欧美激情国产日韩精品一区| 国产伦理片在线播放av一区| 人妻一区二区av| 99国产精品免费福利视频| 亚洲久久久国产精品| 精品少妇黑人巨大在线播放| 国产在线视频一区二区| 女人久久www免费人成看片| 黄色日韩在线| 亚洲av中文字字幕乱码综合| 欧美变态另类bdsm刘玥| 女的被弄到高潮叫床怎么办| 久久久久精品性色| 国产高清有码在线观看视频| 日韩大片免费观看网站| 成人国产麻豆网| 亚洲成人中文字幕在线播放| 一区二区三区精品91| 在线观看av片永久免费下载| 亚洲精品aⅴ在线观看| 午夜福利网站1000一区二区三区| 寂寞人妻少妇视频99o| 亚洲色图av天堂| 欧美激情国产日韩精品一区| 国产白丝娇喘喷水9色精品| 国产白丝娇喘喷水9色精品| 日韩国内少妇激情av| 久久久久国产网址| av天堂中文字幕网| av.在线天堂| 夜夜看夜夜爽夜夜摸| 日韩欧美 国产精品| 成人黄色视频免费在线看| 国产白丝娇喘喷水9色精品| 高清视频免费观看一区二区| 老熟女久久久| 又爽又黄a免费视频| a级一级毛片免费在线观看| 国产成人午夜福利电影在线观看| 内射极品少妇av片p| 小蜜桃在线观看免费完整版高清| 国产欧美日韩精品一区二区| 欧美老熟妇乱子伦牲交| 色5月婷婷丁香| 色综合色国产| 亚洲电影在线观看av| 国内少妇人妻偷人精品xxx网站| 亚洲精品,欧美精品| 少妇被粗大猛烈的视频| 草草在线视频免费看| 亚州av有码| 久久久亚洲精品成人影院| 高清黄色对白视频在线免费看 | 中文天堂在线官网| 18禁裸乳无遮挡动漫免费视频| 欧美日韩视频高清一区二区三区二| 久久精品夜色国产| 18禁动态无遮挡网站| 成人午夜精彩视频在线观看| 亚洲精品视频女| 欧美少妇被猛烈插入视频| 成年人午夜在线观看视频| 伊人久久国产一区二区| 国产人妻一区二区三区在| 亚洲欧美中文字幕日韩二区| 欧美一区二区亚洲| 中文精品一卡2卡3卡4更新| 啦啦啦在线观看免费高清www| 日本与韩国留学比较| 亚洲欧美精品自产自拍| freevideosex欧美| av天堂中文字幕网| 亚洲电影在线观看av| 中文资源天堂在线| 大香蕉久久网| 欧美日韩一区二区视频在线观看视频在线| 两个人的视频大全免费| 中文欧美无线码| 99热网站在线观看| 国产亚洲午夜精品一区二区久久| 久久影院123| 日韩免费高清中文字幕av| 欧美丝袜亚洲另类| 亚洲精品第二区| 韩国av在线不卡| 国产日韩欧美亚洲二区| 在线观看一区二区三区激情| tube8黄色片| 小蜜桃在线观看免费完整版高清| 欧美+日韩+精品| 人妻制服诱惑在线中文字幕| 黄色配什么色好看| 免费大片黄手机在线观看| 成人高潮视频无遮挡免费网站| 亚洲av二区三区四区| 国产高清不卡午夜福利| 免费观看av网站的网址| 午夜福利在线观看免费完整高清在| 大话2 男鬼变身卡| 街头女战士在线观看网站| 国内揄拍国产精品人妻在线| 十分钟在线观看高清视频www | 成人综合一区亚洲| 精品亚洲成a人片在线观看 | 我的女老师完整版在线观看| 亚州av有码| 天堂8中文在线网| 国产黄色视频一区二区在线观看| 国产久久久一区二区三区| 蜜桃久久精品国产亚洲av| av视频免费观看在线观看| 国产成人精品久久久久久| 国产男女内射视频| 观看美女的网站| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av天美| 亚洲第一av免费看| 午夜福利高清视频| 亚洲精品国产成人久久av| 国产亚洲精品久久久com| 亚洲精品成人av观看孕妇| 亚洲精品乱码久久久v下载方式| 亚洲国产精品国产精品| 欧美一级a爱片免费观看看| kizo精华| 制服丝袜香蕉在线| 午夜福利在线在线| 色婷婷av一区二区三区视频| 少妇的逼好多水| 人人妻人人添人人爽欧美一区卜 | 人妻系列 视频| 亚洲丝袜综合中文字幕| 麻豆成人av视频| 人人妻人人澡人人爽人人夜夜| 永久免费av网站大全| 一区在线观看完整版| 各种免费的搞黄视频| 国产精品一区二区在线不卡| 最近中文字幕高清免费大全6| 国产免费一区二区三区四区乱码| 91精品国产九色| 人妻一区二区av| 亚洲激情五月婷婷啪啪| 亚州av有码| 久久人人爽av亚洲精品天堂 | 在线观看av片永久免费下载| 伦理电影免费视频| 在线免费十八禁| 一个人免费看片子| 日本免费在线观看一区| 免费av中文字幕在线| 99热这里只有是精品在线观看| 人人妻人人爽人人添夜夜欢视频 | 欧美日韩一区二区视频在线观看视频在线| 十分钟在线观看高清视频www | 亚洲精品久久久久久婷婷小说| 3wmmmm亚洲av在线观看| 日韩av免费高清视频| 亚洲色图综合在线观看| 免费观看的影片在线观看| 成年免费大片在线观看| 新久久久久国产一级毛片| 91精品一卡2卡3卡4卡| 精品酒店卫生间| 亚洲欧美一区二区三区黑人 | 在线免费观看不下载黄p国产| 国产精品蜜桃在线观看| 精品一区二区三卡| 天堂中文最新版在线下载| 亚洲欧洲国产日韩| 国产精品爽爽va在线观看网站| 六月丁香七月| a 毛片基地| 97在线人人人人妻| 国产伦理片在线播放av一区| 美女主播在线视频| 久久国产亚洲av麻豆专区| 在线 av 中文字幕| 国产精品国产av在线观看| 婷婷色麻豆天堂久久| 男女国产视频网站| 下体分泌物呈黄色| 亚洲第一av免费看| 你懂的网址亚洲精品在线观看| 国内揄拍国产精品人妻在线| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的 | 3wmmmm亚洲av在线观看| 亚洲精品色激情综合| 国产精品久久久久久av不卡| 亚洲人与动物交配视频| 青春草视频在线免费观看| 亚洲国产高清在线一区二区三| 国产一区二区在线观看日韩| 国产综合精华液| 99热这里只有是精品在线观看| 亚洲无线观看免费| 大话2 男鬼变身卡| 日韩大片免费观看网站| 久久6这里有精品| 国产精品久久久久久久久免| 看非洲黑人一级黄片| 免费久久久久久久精品成人欧美视频 | 精品少妇黑人巨大在线播放| 中文字幕av成人在线电影| 十分钟在线观看高清视频www | 观看av在线不卡| 国产精品女同一区二区软件| 麻豆精品久久久久久蜜桃| 国产伦精品一区二区三区四那| 日韩欧美 国产精品| 亚洲av福利一区| 亚洲av免费高清在线观看| 在线观看一区二区三区激情| 免费观看在线日韩| 18禁裸乳无遮挡动漫免费视频| 成人无遮挡网站| 纯流量卡能插随身wifi吗| 亚洲欧美日韩东京热| 国产精品无大码| 91久久精品国产一区二区成人| 亚洲va在线va天堂va国产| 国产中年淑女户外野战色| 青春草国产在线视频| 精品人妻一区二区三区麻豆| 免费少妇av软件| 在线精品无人区一区二区三 | 一本色道久久久久久精品综合| 多毛熟女@视频| av在线播放精品| 久久久久网色| 韩国高清视频一区二区三区| 美女中出高潮动态图| 亚洲av二区三区四区| 国内揄拍国产精品人妻在线| 成人18禁高潮啪啪吃奶动态图 | 国产精品免费大片| 高清欧美精品videossex| 狠狠精品人妻久久久久久综合| 少妇的逼水好多| 交换朋友夫妻互换小说| 亚洲国产欧美在线一区| 一本久久精品| 久久久久久久亚洲中文字幕| 美女脱内裤让男人舔精品视频| 国产精品福利在线免费观看| 国精品久久久久久国模美| 蜜桃久久精品国产亚洲av| 亚洲,欧美,日韩| 两个人的视频大全免费| 成人国产av品久久久| 激情 狠狠 欧美| 丰满少妇做爰视频| 国产欧美日韩一区二区三区在线 | 精品少妇久久久久久888优播| 性色av一级| 99re6热这里在线精品视频| 看非洲黑人一级黄片| 狂野欧美激情性xxxx在线观看| av女优亚洲男人天堂| 精品一区二区三卡| 亚洲精品久久久久久婷婷小说| 成人漫画全彩无遮挡| 国产无遮挡羞羞视频在线观看| 亚洲欧美日韩东京热| 少妇人妻 视频| 18禁裸乳无遮挡免费网站照片| 夜夜爽夜夜爽视频| av网站免费在线观看视频| 国产亚洲最大av| 免费看av在线观看网站| www.av在线官网国产| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 精品一区二区三卡| 天天躁夜夜躁狠狠久久av| 国产亚洲最大av| 韩国av在线不卡| 日韩成人伦理影院| 国产精品无大码| av卡一久久| 国产一区有黄有色的免费视频| 国产免费一级a男人的天堂| 夫妻午夜视频| 国产片特级美女逼逼视频| 老女人水多毛片| 国产高潮美女av| 国产精品精品国产色婷婷| 如何舔出高潮| 亚洲经典国产精华液单| 我的女老师完整版在线观看| 女性被躁到高潮视频| 亚洲精华国产精华液的使用体验| 亚洲欧美日韩无卡精品| 国产免费视频播放在线视频| 丝袜喷水一区| 久久久亚洲精品成人影院| 久久久久久久精品精品| av在线app专区| .国产精品久久| 日韩中字成人| 国产 精品1| 欧美一级a爱片免费观看看| 亚洲欧美精品自产自拍| 99热网站在线观看| 美女福利国产在线 | 夜夜爽夜夜爽视频| 色网站视频免费| 亚洲欧洲国产日韩| 国产精品久久久久久av不卡| 国产精品久久久久久久久免| 久久鲁丝午夜福利片| 九九爱精品视频在线观看| 九草在线视频观看| 99热网站在线观看| 人妻少妇偷人精品九色| 亚洲在久久综合| 一级毛片 在线播放| 欧美xxxx性猛交bbbb| 精品酒店卫生间| 黑人高潮一二区| 少妇人妻精品综合一区二区| 国产探花极品一区二区| 亚洲国产精品专区欧美| 欧美精品亚洲一区二区| 亚洲最大成人中文| 丰满人妻一区二区三区视频av| 国产中年淑女户外野战色| 伦理电影大哥的女人| 一级片'在线观看视频| 欧美成人精品欧美一级黄| 国产永久视频网站| 一级二级三级毛片免费看| 美女中出高潮动态图| 中文欧美无线码| 一区二区三区乱码不卡18| 日韩视频在线欧美| 国产高潮美女av| 国产欧美另类精品又又久久亚洲欧美| 人妻少妇偷人精品九色| 伊人久久精品亚洲午夜| 丰满人妻一区二区三区视频av| a 毛片基地| 伊人久久精品亚洲午夜| 久久久久久久久大av| 熟女人妻精品中文字幕| 另类亚洲欧美激情| 亚洲欧美成人综合另类久久久| 一本一本综合久久| 欧美+日韩+精品| 最近最新中文字幕免费大全7| 免费黄色在线免费观看| 亚洲色图av天堂| 麻豆国产97在线/欧美| 香蕉精品网在线| 欧美xxⅹ黑人| 在线观看国产h片| 高清毛片免费看| 男女边摸边吃奶| 91精品国产国语对白视频| 久热这里只有精品99| 久久久久精品久久久久真实原创| 久久久久精品性色| 黄色一级大片看看| 久久人人爽人人爽人人片va| 97精品久久久久久久久久精品| 人人妻人人澡人人爽人人夜夜| av一本久久久久| 国产乱人偷精品视频| 汤姆久久久久久久影院中文字幕| 欧美日韩国产mv在线观看视频 | 日韩,欧美,国产一区二区三区| 建设人人有责人人尽责人人享有的 | 久久精品人妻少妇| av国产久精品久网站免费入址| 国产精品99久久久久久久久| 看非洲黑人一级黄片| 一区二区三区免费毛片| 精品久久久噜噜| 亚洲一区二区三区欧美精品| 91久久精品国产一区二区成人| 极品少妇高潮喷水抽搐| 亚洲av中文av极速乱| 精品少妇久久久久久888优播| 国产久久久一区二区三区| 免费观看av网站的网址| 夫妻午夜视频| 欧美日韩亚洲高清精品| 最近中文字幕2019免费版| 国产精品精品国产色婷婷| 在线亚洲精品国产二区图片欧美 | 精品一区二区三区视频在线| 黄色一级大片看看| 欧美少妇被猛烈插入视频| 精品亚洲成a人片在线观看 | 亚洲aⅴ乱码一区二区在线播放| 午夜免费观看性视频| 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 一级毛片黄色毛片免费观看视频| 校园人妻丝袜中文字幕| 精品视频人人做人人爽| 色婷婷av一区二区三区视频| 一级av片app| 天天躁夜夜躁狠狠久久av| 午夜福利影视在线免费观看| kizo精华| 大片电影免费在线观看免费| 久久久久久人妻| 亚洲精品久久久久久婷婷小说| videossex国产| 国产精品国产av在线观看| 精品久久久久久电影网| 男人狂女人下面高潮的视频| 午夜福利高清视频| 午夜老司机福利剧场| 大香蕉久久网| 18禁在线播放成人免费| 国产成人a区在线观看| 精品国产乱码久久久久久小说| 精品人妻熟女av久视频| 成人亚洲欧美一区二区av| 偷拍熟女少妇极品色| 免费少妇av软件| 亚洲美女视频黄频| 亚洲成人一二三区av| 亚洲激情五月婷婷啪啪| 毛片女人毛片| 国产高清不卡午夜福利| 一区二区三区四区激情视频| 日韩在线高清观看一区二区三区| 大片电影免费在线观看免费| 老师上课跳d突然被开到最大视频| 亚洲美女黄色视频免费看| 三级国产精品片| 国内少妇人妻偷人精品xxx网站| 午夜福利影视在线免费观看| 国产欧美亚洲国产| 天堂8中文在线网| 熟女av电影| 丰满乱子伦码专区| 18禁动态无遮挡网站| 中文欧美无线码| 下体分泌物呈黄色| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| 91狼人影院| 18禁在线播放成人免费| 国产淫语在线视频| 一本—道久久a久久精品蜜桃钙片| 狠狠精品人妻久久久久久综合| 国产真实伦视频高清在线观看| 国产黄色视频一区二区在线观看| 男女免费视频国产| 日本一二三区视频观看| 久久久久久久久久人人人人人人| 国产av国产精品国产| 伊人久久精品亚洲午夜| 校园人妻丝袜中文字幕| 中文资源天堂在线| 又粗又硬又长又爽又黄的视频| av在线app专区| 国产一区二区在线观看日韩| 国产中年淑女户外野战色| 日韩av不卡免费在线播放| 国产老妇伦熟女老妇高清| 天堂中文最新版在线下载| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 国产精品久久久久久av不卡| 久久久久国产网址| 免费大片18禁| 日韩成人伦理影院| 亚洲精品国产色婷婷电影| 久久精品熟女亚洲av麻豆精品| 国产亚洲91精品色在线| 亚洲三级黄色毛片| 国产成人精品婷婷| 国产色婷婷99| 午夜激情久久久久久久| 国产免费又黄又爽又色| 国产精品人妻久久久影院| 在线观看人妻少妇| 免费不卡的大黄色大毛片视频在线观看| 蜜桃久久精品国产亚洲av| 少妇丰满av| 欧美日韩亚洲高清精品| 看非洲黑人一级黄片| 亚洲国产成人一精品久久久| 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| 天天躁夜夜躁狠狠久久av| 老司机影院成人| 久久热精品热| 麻豆成人午夜福利视频| 日本-黄色视频高清免费观看| 国产亚洲精品久久久com| 黑丝袜美女国产一区| 91精品伊人久久大香线蕉| 黄色日韩在线| 91精品国产国语对白视频| 国产一区二区三区综合在线观看 | 搡女人真爽免费视频火全软件| 香蕉精品网在线| 亚洲美女视频黄频| 成年av动漫网址| 在线观看一区二区三区| 女的被弄到高潮叫床怎么办| 婷婷色综合大香蕉| 如何舔出高潮| 日韩不卡一区二区三区视频在线| 欧美高清性xxxxhd video| 亚洲伊人久久精品综合| 老师上课跳d突然被开到最大视频| 久久久久视频综合| 精品视频人人做人人爽| 精品一区在线观看国产| 中文资源天堂在线| 日韩成人av中文字幕在线观看| 国产在线视频一区二区| 91精品国产国语对白视频| 久久国产精品大桥未久av | 国产精品熟女久久久久浪| 男女边摸边吃奶| 麻豆成人午夜福利视频|