Next Article in Journal
Pannonian Basin Nocturnal Boundary Layer and Fog Formation: Role of Topography
Next Article in Special Issue
Theoretical Chemistry and the Calculation of the Atmospheric State
Previous Article in Journal
Revisiting Forest Effects on Winter Air Temperature and Wind Speed—New Open Data and Transfer Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

O3 Concentration and Its Relation with BVOC Emissions in a Subtropical Plantation

LAGEO, Institute of Atmospheric Physics, Chinese Academy of Sciences (IAP, CAS), Beijing 100029, China
Atmosphere 2021, 12(6), 711; https://doi.org/10.3390/atmos12060711
Submission received: 26 March 2021 / Revised: 25 May 2021 / Accepted: 27 May 2021 / Published: 31 May 2021
(This article belongs to the Special Issue Theoretical Chemistry of Atmospheric Processes)

Abstract

:
An empirical model of O3 is developed using the measurements of emissions of biogenic volatile organic compounds (BVOCs), O3 concentration, global solar radiation, photosynthetically active radiation (PAR) and meteorological variables in a subtropical Pinus plantation, China, during 2013–2016. In view of the different structures of isoprene and monoterpenes, two empirical models of O3 concentration are developed, considering PAR absorption and scattering due to gases, liquids and particles (GLPs), as well as PAR attenuation caused by O3 and BVOCs. The estimated O3 is in agreement with the observations, and validation of the O3 empirical model is conducted. O3 concentrations are more sensitive to changes in PAR and water vapor than S/Q (horizontal diffuse to global solar radiation) and BVOC emissions. O3 is positive to changes in isoprene emission at low light and high GLPs, or negative at high light and low GLPs; O3 is negative to changes in monoterpene emissions. O3 are positive with the changes of PAR, water vapor and S/Q. It is suggested to control human-induced high BVOC emissions, regulate plant cutting, and reduce NOx and SO2 emissions more strictly than ever before. There are inverted U-shape interactions between O3 and its driving factors, and S/Q controls their turning points.

1. Introduction

Ozone is one of the important greenhouse gases influencing climate change on global scale [1], and a key constituent of pollutants produced through chemical and photochemical reactions (CPRs) in the atmosphere. Emissions of biogenic volatile organic compounds (BVOCs) dominate the global total volatile organic compounds entering the atmosphere [2,3]. BVOCs, sensitive precursors of O3, are oxidized to produce new gases, liquids and particles (GLPs) in the atmosphere, e.g., blue haze in the forest [4], O3 and secondary organic aerosol (SOA) [5,6,7,8,9,10]. In addition, volatile organic compounds (VOCs, including anthropogenic VOCs (AVOCs and BVOCs) and nitrogen oxides (NOx) are the main precursors of ozone [5]. SOA plays a significant role in cloud condensation nuclei formation and cloud formation [11,12]. OH radicals react with almost all atmospheric trace gases and are rapidly recycled. VOCs, OH radicals, NOx, SO2 and other GLPs react through CPRs and produce new GLPs (O3, PM2.5, SOA, HCHO, etc.). During CPRs, solar ultraviolet and visible energy (UV, VIS) are absorbed and utilized by GLPs through OH radicals and H2O [13], and OH radicals are important reactants and bridges connecting almost all GLPs in the atmosphere.
There are some limitations and uncertainties associated with OH radicals, O3, BVOCs and SOA. For example, OH measurement and model simulation [14,15,16]; the missing and unexplained OH sink/reactivity [16,17,18,19]; O3 simulation [20,21,22,23]; BVOC emission measurements and simulations [24,25]; unknown mechanisms of SOA formation from BVOC oxidation [26,27,28], missing UV energy [13]. Thus, more challenges are still in realizing the atmosphere, biosphere and their interactions [29,30,31,32], and better understanding the photochemistry of O3–BVOCs. In view of the large uncertainties in BVOC measurements and model simulations [2,24,25], accurate simulations of the total amount of O3 and SOA still have large uncertainties, though the air quality models, such as weather research forecasting coupled with Community Multi-scale Air Quality (WRF-CMAQ) and chemistry model (WRF-Chem), or other models of similar structures and working principles have been used to simulate O3 and SOA and much progress have been achieved. Are there any important driving factors in controlling the changes in the relations in O3–BVOCs and O3–aerosol? More investigations on interactions and mechanisms in O3–BVOCs, O3–aerosol (e.g., SOA) and O3–PAR are necessary, especially under natural atmospheric conditions. It is an urgent need to accurately evaluate the total amount of O3 and SOA contributed from BVOC oxidation. Finally, the interactions and mechanisms in O3 with its driving factors (BVOCs, PAR, aerosols, etc.) under realistic atmospheric conditions can be used as references in future model studies and pollution control of O3.
BVOC emissions, O3, fine particle formations and solar radiation change synchronously, i.e., all kinds of GLPs and solar radiation are multiple processes that interact. The sun provides UV and visible energy and triggers the changes of GLPs (gases, H2O, particles, clouds, haze, etc.) [13]. In view of the challenge in calculating all atmospheric constituent amounts, especially their total PAR absorption and indirect utilization, energy is an important source for all atmospheric constituents and controls their changes, so an energy method is selected. The progress of measurements and an empirical model of BVOC emissions make it possible to develop an empirical model of O3 concentration and study the photochemistry of O3 and BVOCs. The main objective of this study is to investigate the interactions between O3 and its driving factors, e.g., BVOCs, absorbing substances, scattering substances and PAR. Based on the measurements of O3, BVOC emissions and solar radiation, two empirical models of O3 concentration for considering the individual roles of isoprene or monoterpenes are proposed and evaluated. The sensitivity of the O3 changes with its driving factors, and the relationships between the calculated and observed O3 and its driving factors are investigated. Several methods to reduce O3 and fine particle pollution are suggested.

2. Instrumentation and Methods

2.1. Site Description

O3 concentrations, BVOC emissions and solar radiation were measured at the Qianyanzhou subtropical Pinus plantation, Taihe County, Jiangxi province, China (26°44′48″ N, 115°04′13″ E, 110.8 m) from 22 May 2013 to 4 January 2016. In the center of the flux tower, the forest coverage is about 70% within 100 Km2, and the mean forest canopy height is 18 m. Local dominant trees are Pinus massoniana, Pinus elliottii, Cunninghamia lanceolata. The shrubs are Loropetalum chinense, Adinandra millettii, and Lyonia compta. The average slope is 2.8°–13.5°. The annual average temperature is 17.9 °C [33]. The annual global solar irradiance is 4578 MJ m−2, and the annual photosynthetically active radiation (PAR) is 7997 mol m−2 in 2013 [34]. A 45-m tower was built in this Pinus plantation for flux, solar radiation and meteorological measurements.

2.2. Methods and Instruments

O3 was measured using an ozone monitor (Model 205, 2B Technologies Inc., Boulder, CO, USA) installed at the same height (23 m) as the relaxed eddy accumulation (REA) and sonic anemometer. The solar radiation measurement system consists of 3 spectral sensors at 270–3200, 400–3200 and 700–3200 nm, respectively, and a recorder with an accuracy of 5% (model TBQ-4-1, 322 Institute of Jinzhou, Jinzhou, China). Global horizontal radiation was observed by the sensor at 270–3200 nm (referred to as Q). Direct normal solar radiation (referred to as D) was measured with a radiometer (270–3200 nm, Model TBS-2, China). Diffuse horizontal radiation (referred to as S) was derived from Q-D × cos(Z), where Z is the solar zenith angle (degree). UV and visible radiation were derived from the difference over 400–3200 nm, 700–3200 nm, from 270–3200 nm, respectively. The solar radiation sensors were placed at the top of a building (10 m above the ground), located at the Qianyanzhou Experimental Station of Red Soil and Hilly Land, Chinese Academy of Sciences, 800 m away from the flux tower. Temperature and relative humidity were also measured in this Pinus plantation using a HOBO weather station (Model H21, Onset Company, Menlo Park, CA, USA) during this campaign. The NOx concentrations were measured at the bottom of the tower using a NO/NO2/NOx analyzer (Model EC9841, Ecotech Company, Knoxfield, Australia) from 12 June–16 October, 2014.
BVOC emission fluxes were measured using an REA system. It consists of a three-dimensional sonic anemometer, a data logger, and a data acquisition and control unit that were used to collect air samples onto absorbents (Tenax GR and Carbograph 5TD) in stainless steel cartridges [35,36]. The REA system was located at a platform (23 m above ground level). Air samples were collected from 23 July 2014 using a new gradient measurement system, including the same three-dimensional sonic anemometer and a data logger located at the same position during the REA and gradient measurements, and some new pumps (HL-2 sampler, Beijing municipal institute of labor protection, Beijing, China) for collection of air samples at 20 m and 28 m. Air samples were collected in 30 min using the REA and gradient methods. More information about the system of the REA and gradient, cartridge usage, and solar radiation is reported in [34].
The emission fluxes of a given BVOC species (Fi) from the REA technique is Fi = bσw (Cup − Cdown), where σw is the standard deviation of the vertical wind velocity, b is an empirical coefficient, and Cup and Cdown are the concentrations (μg m−3) of the BVOC species in the up- and down-stream cartridges, respectively. The Fi measured using the gradient method is: Fi = Kdiff × (ΔC/Δz), where Kdiff is the eddy diffusion coefficient; Kdiff = (k) × (u*) × (z − d) for neutral atmospheric stability, k = 0.4 (von Karmons constant); u* is the friction velocity (m/s); z is the measurement height (geometric value, m); d is the displacement height and assumed to be 2/3 canopy height (m). The mean canopy height (18 m) was used for calculation [34].
Air samples collected using the REA and the gradient method were shipped to a lab at the National Center for Atmospheric Research (NCAR) in Boulder, U.S.A. and in Beijing at the Institute of Atmospheric Physics, Chinese Academy of Sciences (IAP, CAS) for analysis. The procedures for sample analyses have been reported in [34,37,38,39]. The procedures for sample analyses by gas chromatographs equipped with flame ionization detector at the laboratory (IAP, CAS) in Beijing, China, were similar to those described by Greenberg et al. [36]. A summary of all sampling periods as well as O3 measurements are shown in Table 1.

2.3. Empirical Model of O3 Concentration

Based on the empirical model of BVOC emissions for Qianyanzhou Pinus plantation, China [34], and further application of the principle of PAR energy balance, isoprene and monoterpene roles are considered in the development of the O3 empirical model, respectively. Four processes associated with PAR transfer above the canopy are considered:
(1) PAR attenuation due to BVOCs. Isoprene and monoterpene terms represent PAR attenuation by isoprene and monoterpenes, which are assumed to obey the exponential law and described as e a 1 I S O m and e a 2 M T m , respectively. a1 and a2 are the absorption coefficients of isoprene and monoterpenes (presumed to be united, mg−1 m2), ISO and MT are the total isoprene and monoterpene emissions in the sampling period, ISO = t × EFI × 0.1, MT = t × EFM, EFI and EFM are isoprene and monoterpene emission fluxes (mg m−2 h−1), 0.1 is a normalizing coefficient for isoprene, t (0.5 h) is the sampling period, and m is the optical air mass in the center of the averaging window.
(2) PAR absorption and consumption (photochemical term) due to other GLPs (e.g., NOx, SO2) when they are taking part in CPRs except O3 and isoprene or O3 and monoterpenes, with emphasis on BVOCs through OH radicals and H2O. For example, when isoprene is displayed in the empirical model (Equations (1) and (3)) but monoterpenes are not, the PAR utilization due to monoterpenes is considered in this term and connected through the CPRs with OH radicals. It is similar to monoterpenes when they are considered in the empirical model (Equations (2) and (4)). This term also includes chemical energy converted from PAR absorption by the absorbers. The fraction of global solar radiation absorbed by water vapor in 0.70–2.845 μm is calculated as ΔS′ = 0.172 (mW)0.303, W (water vapor content in the whole atmospheric column) = 0.021E × 30, and E is the mean water vapor pressure (hPa) at the ground during the sampling period. Assuming only water vapor is considered and the atmosphere is plane-parallel, the solar radiation on the horizontal surface is I0 cos(Z) − ΔS′ = I0 e kWm cos(Z), i.e., e kWm = 1 − ΔS′/Io, where Io = 1.94 cal min−1 cm−2 (1367 W m−2, solar constant), and k is the mean absorption coefficient of water vapor in 0.70–2.845 μm.
In the visible region (400–700 nm), the important OH radicals in the troposphere are produced from several ways: (a) NO2 + hv→NO2*, NO2* + H2O→HONO + OH, which is 50% of that assumed by the traditional O(1D) + H2O reaction [40]. The reaction rate reported by Li et al. [40] impacts O3 formation in high NOx emissions significantly [41]. (b) The photodissociation of CH3OOH [42], (c) the photolysis of O3∙H2O clusters [43], (d) the photochemistry of SO2 “on water” [44] and (e) NO2 and H2O upon photoinitiated by 410 nm light [45]. The “missing OH sinks” in the forest [17,27,46] and related mechanisms and other OH sources should be investigated continuously [47] to improve our deep understanding of the OH photochemistry in the visible region. In view of the above mechanisms, water vapor is used to represent OH radicals, NO2*, H2O, SO2 and other chemical constituents.
During the CPRs, gases (e.g., BVOCs, O3 and NOx), liquids (e.g., H2O) and particles (mainly BVOC oxidation products, i.e., SOA) react with each other, especially OH radicals, change in concentrations and gas, liquid and particle phases. Considering the concentrations of each GLPs cannot be measured, and many mechanisms associated with O3, BVOCs, SOA formation are unknown at present, the direct PAR absorption and indirect PAR utilization due to GLPs are calculated by the photochemical term using an energy method because water vapor plays an important reactant/bridge in CPRs, especially in the processes of OH radical formation. Thus, water vapor is used as a surrogate to represent all absorbing GLPs, and the photochemical term expresses the entire energy use during the sampling period, associating with all atmospheric constitutes, e.g., NOx, SO2 and other VOCs, when they are not explicitly described in the empirical Equations (1)–(4), apart from the isoprene or monoterpene terms. Whether the photochemical term be representative of the total energy use caused by absorbing GLPs is explained in Section 4.1 when the O3 empirical model considers NOx roles and the UV empirical model considers 2 terms and 3 terms.
(3) O3 absorption in the region (400–700 nm) (O3 term) is described as e 0.1 k 3 O 3 t m , where k3 is the O3 absorption coefficient in 400–700 nm and k3 = 0.053 × 105 Pa−1 cm−1; 0.1 is a normalized coefficient. O3 term can be used to study the interactions between O3 and BVOCs in the visible region and is an application from a previous study, the interactions between O3 and NOx in the UV region [48].
(4) The scattering of GLPs (scattering term) is described as e S/Q. Scattering factor S/Q describes the relative amount of GLPs in the atmosphere, including gas molecules, aerosols, clouds, haze and rain. In more detail, when air mass moves from other regions to the study region, the total scattering from the air mass and other GLPs (e.g., NOx, SO2, VOCs, O3 and aerosols) in the atmospheric column during the sampling period (0.5 h) were described using this term objectively.
Isoprene (C5H8) and monoterpene ((C5H8)2) have different carbon structures and behaviors [49] and are two dominant BVOCs emitted from plants, contributing about 50% and 15% of the global BVOCs, respectively [50]. They also show different contributions in the formation of O3 and SOA [51,52,53,54]. Therefore, it is necessary to study their individual interactions with O3. It should be pointed out that BVOCs, O3, photochemical and scattering terms express the whole atmospheric column amount and their total energy in the whole atmospheric column during the sampling period. In specific, isoprene and monoterpene terms express their individual PAR extinction by the law of exponential attention. Meanwhile, the photochemical and scattering terms describe the total absorbing and scattering energy associated with entire column GLPs, including air mass moves to the study region in the sampling period. The energy distribution in each term is objectively quantified by analyzing the observational data under realistic atmospheric conditions. PAR is an energy source that drives the absorbing and scattering processes, as well as BVOC emissions, no matter how chemical constituents change in the atmosphere. The PAR energy balance at a horizontal plane during the sampling period is described by Equations (1) and (2) for considering the roles of isoprene or monoterpenes, respectively:
PAR = A 1 e a 1 I S O m × cos ( Z ) + A 2 e kWm   ×   cos ( Z ) + A 3 e 0.1 k 3 O 3 t m ×   cos ( Z ) +   A 4 e S / Q + A 0
PAR = B 1 e a 2 M T m × cos ( Z ) + B 2 e kWm   ×   cos ( Z ) + B 3 e 0.1 k 3 O 3 t m ×   cos ( Z ) +   B 4 e S / Q + B 0
The O3 empirical models considering the roles of isoprene or monoterpenes are:
e 0.1 k 3 O 3 t m × cos ( Z ) = A 1 PAR + A 2 e a 1 I S O m   ×   cos ( Z ) + A 3   e kWm   ×   cos ( Z ) + A 4   e S / Q + A 0
e 0.1 k 3 O 3 t m × cos ( Z ) = B 1 PAR + B 2 e a 2 M T m   ×   cos ( Z ) + B 3   e kWm   ×   cos ( Z ) + B 4   e S / Q + B 0
where cos(Z) = 1/m. In order to obtain the optimal PAR energy interactions and distributions in the whole atmosphere between O3 and its driving factors that represent most realistic situations and reduce the influences of errors in the measurements of O3, emission flux and solar radiation, strict criteria for observational data were used in the analysis: (1) BVOC emission fluxes (i.e., isoprene and monoterpenes) less than twice the standard deviation of measured values, respectively, and (2) solar zenith angle < 55°, (3) S/Q < 0.5 (i.e., relative clean atmospheric conditions and low cloudiness).
The development of the empirical model of O3 was based on the measurement and empirical model of BVOC emissions (EMBE) at this subtropical forest [34]. According to the above data selection criteria, 18 and 8 samples were obtained and used in the EMBE models for considering isoprene and monoterpenes, respectively. It means these data are most representative of the ideal interactions between BVOCs and their driving factors under good light and atmospheric conditions (see Section 3.2.1). It also provides a reliable dataset to further study the interactions between O3 and BVOCs. Later, 12 and 11 samples were determined by combining these BVOCs and O3 data for considering the roles of isoprene and monoterpenes, respectively. It still aims to find and determine the basic law in the energy interaction between O3 and its driving factors by using more strict criteria data, similar to that in the EMBE model development [34]. The empirical model of O3 is evaluated in Section 3.1 and 4.3 and modified in Section 4.3. All terms (except PAR) were normalized individually, so as to express their energy equally and objectively [48].
Water vapor as a function of air temperature and relative humidity was calculated using an equation described by Lowe [55], i.e., water vapor can represent air temperature and relative humidity to some extent. In addition, water vapor can represent a dynamic balance between the atmosphere, plants and soil, and the changes of water and water vapor when they take part in CPRs. Therefore, water vapor has more important roles than the air temperature and relative humidity and is used in the empirical model as a proxy of air temperature and relative humidity.
During the sampling period (0.5 h), when air mass associated with the wind speed and direction transports to the study region, the concentrations of BVOCs, O3 and other GLPs (NOx, SO2, etc.), as well as PAR change simultaneously. Their individual roles are described in different terms (i.e., O3, isoprene or monoterpenes, photochemical, scattering terms and PAR) and considered in model development and later simulation of ozone. It should be noted that all terms in Equations (1)–(4) describe the total column amount of all parameters (e.g., O3, BVOCs, PAR, water vapor and S/Q) in the sampling period. Their representative area can be several decade km2, as the large area representatives of the solar radiation. In a further application of the O3 and BVOC empirical models [34], they describe the O3 and BVOC emissions for the entire forest region from this half-hour to the next half-hour in the daytime, and these empirical models describe the air mass transportation and associated chemical and photochemical processes continuously when the air mass is passing through the study region.

3. Main Results

3.1. O3 Estimations and Its Validations

The observed O3 concentrations, BVOC emission fluxes measured by the REA technique, solar radiation and meteorological parameters were analyzed to develop the O3 empirical models and determine the coefficients and constants (Equations (3) and (4), samples, n = 12 and 11 for considering isoprene or monoterpenes). The statistical results calculated using Equations (3) and (4), i.e., coefficients and constants, coefficient of determination (R2), average and maximum of the absolute relative bias, δavg and δmax (%) (δ = |ycal − yobs| × 100/yobs, ycal and yobs are calculated and observed O3 concentrations), normalized mean square error (NMSE = ( y c a l y o b s ) 2 ¯ /( y c a l ¯ × y o b s ¯ ) ) [56] and standard deviations of the calculated and observed O3 concentrations (σcal and σobs), are shown in Table 2.
Generally, the simulated O3 concentrations agreed with the observation (Table 2, Figure 1 and Figure 2) for considering the roles of isoprene and monoterpenes. The mean calculated and observed O3 concentrations were 46.85 and 46.88 ppb with a relative bias of 0.07%, and the root mean square error (RMSE) was 3.94 ppb for situation A, and the corresponding values were 38.74 and 38.76 ppb, 0.06% and 4.37 ppb for situation B.
To evaluate the empirical model performance and reduce the simulation error, the air mass used in BVOC and O3 terms in the early morning and late afternoon was doubled to meet a little larger optical air mass in view of the current atmospheric conditions, according to previous simulations of BVOC emissions in some forests. In view of the REA measurements of BVOC emissions are extensively carried out in typical forests in China, e.g., a temperate forest, subtropical Pinus and bamboo forests [34], the observed dataset obtained using the REA technique were used for validation, as well as model development. Two methods were used to evaluate the O3 empirical model. Firstly, the mean half-hourly variation was calculated by averaging half-hourly O3 at each time of the day for 22 May 2013–11 Sep 2013 (BVOC emission measurements using the REA technique, n = 104) and is given in Figure 3. The observed O3 concentrations within 2σ of the average were selected. The simulated O3 concentrations overestimated the observations by about 30%, and their averages were 53.8 and 41.2 ppb, respectively. It is reasonable considering the uncertainties in BVOC measurements of about 25% [24] as well as O3 and solar radiation measurements. Secondly, the calculated and observed O3 concentrations in single measurements (e.g., in half-hour) in 2013 were also compared, the mean calculated and observed O3 concentrations were 57.7 and 41.8 ppb (n = 104), respectively, i.e., the simulated O3 overestimated the observed O3 by 36%. The NMSE and RMSE values were 0.1534 and 22.9 ppb. The increased relative bias (i.e., 30% and 36% compared to 0.07% and 0.06%) was caused by using more observed emission fluxes with larger uncertainties in measurements and O3 simulations under all sky conditions (low PAR, high GLP loads, such as clouds and aerosols, see Section 3.2.1, whereas the data used in the development of the empirical model were measured under optimal atmospheric conditions (i.e., high PAR, clean atmosphere and low S/Q) and had smaller uncertainties in observations and simulations. In addition, a large simulation error was caused by the changes in S/Q (GLP loads, including aerosols and clouds) than isoprene and monoterpenes (see Section 3.2.1 and Section 3.2.2), especially for high GLP load conditions.
Similarly, validation for the O3 empirical model considering monoterpene roles was conducted using the observational data in 2013. (1) The mean half-hourly variation was estimated and given in Figure 4 (n = 104). The simulated O3 overestimated the observed O3 (43.8 vs. 41.0 ppb) by 7%. (2) The calculated and observed O3 mean concentrations were 45.3 and 41.6 ppb (n = 104) for single measurements during 2013, respectively, i.e., the simulated O3 overestimated O3 by 9%. The NMSE and RMSE values were 0.1063 and 18.1 ppb, respectively. These uncertainties were reduced compared to that considering the isoprene role, which is mainly due to the monoterpene emissions are the summation of main monoterpene compositions and the random error reduction in the measurements. The comparisons of RMSE values for the O3 empirical model and chemical models are shown in Section 4.3.
Based on the above validations, the simulated O3 concentrations were in reasonable agreement with those observed, and the empirical models (Equations (3) and (4)) can be used to estimate O3. The representative relationships for O3 and its driving factors were also quantified.

3.2. Sensitivity Study of O3 to Its Affecting Factors

To investigate the responses of O3 to its factors under different situations thoroughly, the BVOC emission data (i.e., half-hour measurement) used in the development of the BVOC empirical model (described by REA 12), measured by all REA (described by REA) and all gradient techniques (described by GRA) were selected for the sensitivity tests. These three situations reflect different light and atmospheric conditions (Section 3.2.1).

3.2.1. O3 Responses with Its Affecting Factors Considering the Roles of Isoprene

The sensitivity test was performed using a half-hourly observed dataset under realistic atmospheric conditions, i.e., how simulated O3 responses with one driving factor (isoprene, or PAR, E, S/Q) using Equation (3) when the others keep at their original values.
For the above three situations, the average changes in O3 concentrations (%) caused by changes (%) of each factor while keeping all others at their original levels are shown in Figure 5. When extreme changing rates of O3 were removed, sample points were about 12, 111 and 175 for REA 12, REA and GRA, respectively, when considering isoprene role; when considering monoterpene roles, the sample points were 11 (expressed as REA 11), 110 and 150, respectively.
The mean responses of O3 concentration with the changes of each driving factors (isoprene emission, PAR, E and S/Q) by 20% are reported in Table 3. (1) For REA 12, REA and GRA measuring periods, O3 concentrations were more sensitive to E than PAR, S/Q and finally isoprene. E (as well as temperature) and PAR are dominant driving factors and make larger changes in O3 than other factors, reflecting that O3 is strongly produced/destroyed through chemical and photochemical processes triggered by light and temperature, together with water and water vapor supply associated with the OH production. (2) It is well known that O3 formation is contributed by BVOC oxidation at most sites in China [22,23,57,58,59,60,61,62]. O3 responds (a) negatively with the changes of isoprene (REA12) in good light and atmospheric conditions, i.e., high PAR, high temperature and low S/Q (low GLP loads), and (b) positively with the changes of isoprene (REA, GRA) in bad light and atmospheric conditions, i.e., low PAR, low temperature and high S/Q (Figure 5, Table 3 and Table 4). Similar modeled and observed results, i.e., positive and negative responses, are also reported in [63,64]. Therefore, the positive or negative response of O3 to the isoprene is dependent on light and atmospheric conditions. More and long-time measurements in different forests are necessary to better understand the complex O3–BVOCs interactions. O3 responses to the change of isoprene were larger for GRA than REA and finally REA12, corresponding well to the low, medium and high states of GLPs–PAR, i.e., high GLPs and low PAR (S/Q = 0.75, PAR = 1.26 mol m−2), medium GLPs and medium PAR (S/Q = 0.59, PAR = 1.88 mol m−2), low GLPs and high PAR (S/Q = 0.35, PAR = 2.75 mol m−2). These three situations represent three states: high GLP loads with low air temperature (T) and high relative humidity (RH); medium atmospheric conditions with medium T and RH; clean atmospheric conditions with high T and low RH. Therefore, O3 showed different responses to isoprene, which depended on the state of atmospheric substances and PAR energy.
It also reveals that O3 is destroyed (−0.9 ppb) by O3 photolysis and reactions with OH radicals and other GLPs with the increase of isoprene when PAR is high at 2.75 mol m−2 (1527.8 μmol m−2 s−1). The chamber experiment shows a similar result; the OH radical initiated photooxidation of isoprene produces methyl vinyl ketone, methacrolein and formaldehyde [65]. O3 depleted a little (−0.3 ppb) with the increase of isoprene emission when PAR is at medium level, 1.26 mol m−2 (700 μmol m−2 s−1). Large O3 is produced (7.8 ppb) with the increase of isoprene emission when PAR is low at 1.26 mol m−2 (700 μmol m−2 s−1). It is clear that PAR = 1.80 mol m−2 (1000 μmol m−2 s−1) or its corresponding state of GLPs–PAR is a controlling/turning point for the positive/negative response of O3 to isoprene (Table 3 and Table 4).
The increase of O3 is larger than its decrease when isoprene changes at the same rate (comparing the situation REA 12 to GRA, Table 3), and it is mainly caused by the initial isoprene emissions, 0.39 and -0.61 mg m−2 h−1. Thus, it is most effective to control isoprene emission at low O3 (39.1 ppb) than at high O3 (46.9 ppb) and control human-induced BVOC emissions, i.e., plant-cutting and biomass burning. For example, plant-cutting (branches and grasses) is suggested to be carried out after 16:00 in big cities, as injured leaves and the grass cutting enhance BVOC emissions dramatically. This mechanism is a reference in future O3 pollution control.
It is interesting that the ratios of REA to GRA for the O3 response to each factor were similar to the ratios of REA to GRA for the average of each factor (Table 3 and Table 4), e.g., the above two ratios were −0.1 and −0.2 for isoprene, 1.4 and 1.5 for PAR, 1.3 and 1.4 for E, 0.6 and 0.8 for S/Q, meaning that O3 responses to its driving factor strongly depend on the mean values of the corresponding factor, and the REA-measured BVOC emission fluxes and gradient techniques have similar features for O3–isoprene photochemistry. Similarly, these features were also found for O3–monoterpenes photochemistry (Table 4 and Table 5).
(3) The mean response of O3 concentration to the changes of PAR was positive, indicating O3 formation is driven by solar energy. (4) The mean response of O3 to the changes of water vapor was positive, revealing that O3 formation is beneficial from sufficient and increased water supply, and more OH radicals are produced from H2O. 5) The mean response of O3 concentration to the changes of S/Q was positive, reflecting fine particles are produced with O3 through CPRs. The O3 responses to all affecting factors were linear, except isoprene.

3.2.2. O3 Responses with Its Affecting Factors Considering the Roles of Monoterpenes

Under realistic atmospheric conditions, the O3 responses with each affecting factor (monoterpenes or PAR, E and S/Q) were calculated using Equation (4). The calculating results are given in Figure 6 and Table 5.
For REA 11, REA and GRA measuring campaigns, O3 concentrations were more sensitive to the changes of PAR than E, S/Q and finally, monoterpenes. PAR and E (as well as temperature) are still dominant driving factors and result in larger changes of O3 than other factors. In general, O3 responses with the changes of monoterpenes were negative for three situations, linear for REA11 and REA and non-linear for GRA. A similar modeled result is reported by Nishimura et al. [52]. O3 responses with the changes of monoterpenes were larger for REA11 than REA and smaller for GRA, corresponding well to their clean levels of the atmospheric conditions, which demonstrates that the cleaner of the atmosphere, the more production of O3 contributed from high monoterpene emissions and their oxidation.
For monoterpenes, the mean responses of O3 concentration to the changes of PAR, water vapor and S/Q factors were also positive, indicating that they have similar interactions with O3 as isoprene does with O3. It should be emphasized that the CPRs are homogeneous and heterogeneous, and S/Q can represent the fine particle formation when S/Q and cloud amounts are low.
When monoterpenes increased by 20%, O3 responses (in percent and ppb) with the changes of monoterpenes were larger (REA11 and REA) in good light and clean atmospheric conditions (Figure 6 and Table 5) and smaller (GRA) in low light and polluted atmospheric conditions (Figure 6 and Table 5). The O3 responses with other driving factors (except monoterpenes) were linear and positive.
In general, the responses of O3 to its driving factors are similar when considering the roles of isoprene and monoterpenes and also similar to that in a subtropical bamboo forest in China [66].
Comparing the response of O3 to the scattering substances (expressed by S/Q), the changing rate of O3 was larger for monoterpenes than isoprene. This similar feature is also found in a subtropical bamboo forest [66]. Both indicate that monoterpenes are more productive for the formation of fine particles and O3 than isoprene under the natural atmosphere. Similar results are observed using the chamber experiments, i.e., SOA yields from the oxidation of monoterpenes are larger than isoprene [51,53]. The response of O3 to monoterpenes was larger than isoprene for REA measurements (Figure 5 and Figure 6), and it is in agreement with the studies that monoterpenes have higher O3 formation potential than isoprene [54,67].
The linear fitting line between O3 and its affecting factors for the O3 sensitivity test is described by DO3 = C1 × Dfactori + C0 and shown in Table 6 and Table 7. DO3 and Dfactori denote the difference of O3 and factori, respectively. C1 and C0 are the coefficient and constant, respectively. A comparison of the two situations in model development considering the roles of isoprene (REA12) and monoterpenes (REA11) provides the following conclusions: (1) the negative response of O3 to the changes of isoprene was larger than monoterpenes, indicating more O3 is produced/destroyed from isoprene than monoterpenes through CPRs under good light and atmospheric conditions. This characteristic still existed for the GRA situation, with smaller O3 changes under bad light and atmospheric conditions, except the positive response of O3. In contrast, the negative response of O3 to the changes of monoterpenes was larger than isoprene for REA measurement, and the corresponding mean values were PAR= 1.88 mol m−2 (1044.4 μmol m−2 s−1), T = 32.6 °C, RH = 63% and S/Q = 0.59. These light and atmospheric conditions are best suitable for O3 photochemical production and close to the turning point (S/Q = 0.55, Section 4.3), the best interactions between O3 and its influencing factors. It is obvious that the larger sensitivity of O3 to the changes of isoprene or monoterpenes depends on the light and atmospheric conditions (Table 3, Table 4 and Table 5, Figure 5 and Figure 6). (2) O3 changes with the changes of other driving factors were positive. (3) PAR and E are the primary important factors to make O3 changes.
The O3 response rates (C1/C0) when considering the roles of isoprene or monoterpenes, respectively, are given in Table 8 and Table 9. It is more interesting that for three situations, O3 response rates to isoprene were larger than monoterpenes, and O3 response rates to other corresponding factors (PAR or E and S/Q) were the same for the roles of isoprene or monoterpenes, respectively. Therefore, O3 response rates per unit of O3 to its key influencing factors (PAR or E and S/Q) are the same for all kinds of atmospheric conditions. This characteristic may be a useful reference to evaluate O3 formation from BVOC oxidation.
O3 concentrations were most sensitive to the changes in PAR and E (representative of air temperature and relative humidity) than S/Q and isoprene/monoterpenes, implying that total energy (i.e., PAR, sensible and latent heat) is the dominant driving factor in O3 photochemistry (also shown in Section 4.4), then the scattering factor, and finally the BVOCs. BVOCs are recognized as important O3 precursors, but their energy role in PAR utilization is the smallest. The stronger correlations between O3 and its driving factors (PAR, T, RH and E) also show that energy plays a key role in controlling O3 formation and destruction (Section 4.4). Similarly, model results also show that changes in atmospheric conditions dominate the interannual variations of O3 and SOA, and the interannual variations in BVOCs alone lead to small differences (2%–5%) in the calculated O3 and SOA in the summer [68].

4. Discussion

4.1. O3 Empirical Model

Based on the energy balance between PAR and its transfer processes associated with GLPs in the whole atmosphere, especially considering the PAR attenuations of O3, isoprene or monoterpenes, their energy interactions and distributions were determined in good light and atmospheric conditions. O3 empirical models considering isoprene and monoterpene roles give reasonable estimates of O3, and their advantages are that the interactions between O3 and its driving factors, PAR, E, isoprene/monoterpenes and atmospheric GLPs (especially fine particles at low S/Q), can be studied extensively.
It is found that the observed PAR and UV have positive correlations with the photochemical term [13], but water does not absorb visible and UV radiation according to previous knowledge (discarding the debate of water absorption in the UV region); thus, the photochemical term represents absorption and indirect use by GLPs through OH radicals during CPRs. More detailed explanations are reported in the UV region [13]. In short, OH radicals are produced in many ways in the visible region (Section 2.3), e.g., H2O plays an energy use/transfer bridge role in OH radical formation through NO2* + H2O→OH, implying that H2O utilizes the visible energy from NO2*, which is similar to that in other situations (b–e, Section 2.3). Similarly, H2O utilizes UV energy from O3 in OH radical formation because of a strong positive correlation between UV and the photochemical term [13].
Under S/Q < 0.8 conditions, strong positive correlations were also found between observed monthly UV, VIS, PAR and water vapor pressure at this forest (n = 14, corresponding to 2779 hourly values), their correlation coefficients were 0.953, 0.927 and 0.897, respectively, at the confidence level of 0.001. These results reveal the point of view that UV and VIS utilization is caused by all GLPs through photochemical reactions with OH radicals and H2O, and UV plays more important roles than VIS because of its higher frequency and higher energy.
It is known that OH radicals react with almost all atmospheric GLPs, especially VOCs [69], and visible energy is absorbed and consumed by GLPs in a single GLP phase and gas–particle conversions during CPRs. Many GLPs (O3, NO2, glyoxal, CH3CO radical, NO3 radical, OClO, CHOCHO, biacetyl, butenedial, BC and other aerosols) are visible radiation absorber [13,70,71]. Others without visible radiation absorption react with OH radicals, and these absorbers, thus, consume visible radiation indirectly. The most important thing is that OH radical recycles quickly. This part energy is expressed by the photochemical term and determined by analyzing observational data and using the multiple-fitting. The photochemical term is an application from a previous study in the UV region [13] to the visible region.
It is necessary to discuss the meaning of the photochemical term displayed in Equations (1)–(4), direct absorption and indirect utilization due to all GLPs in the whole atmospheric column, except O3 and isoprene or O3 and monoterpenes. When a UV empirical model is used to calculate UV at the surface and consider (1) two terms, i.e., photochemical and scattering terms, which are described similar to this study; (2) three terms, photochemical, O3 and scattering terms (as O3 is one important absorber in the UV region); the relative errors of monthly UV are 4.1% for two terms and 3.7% for three terms; the annual mean UV loss (UV at the top of the atmosphere and UV at the ground) caused by photochemical terms using a two-term equation, and O3 and photochemical terms using a three-term equation are 0.57 and 0.54 MJ m−2, respectively. These similar results for using two-term and three-term equations indicate that photochemical terms can express the role of O3 when O3 is not explicitly displayed in the two-term equation [13]. This feature was assumed to be existed in the visible region, i.e., photochemical terms described PAR utilization due to all GLPs in the whole atmospheric column except O3 and isoprene in Equation (1), or O3 and monoterpenes in Equation (2), respectively. Therefore, Equations (3) and (4) were used to study the interactions between O3 and isoprene and O3 and monoterpenes, respectively.
One important issue should be mentioned that there is no evident correlation between the absorbing term and scattering term, and their correlation coefficients were 0.368 for using the data of O3 and isoprene in the development of the O3 model (n = 12), and 0.362 for S/G < 0.8 and 0.059 for S/G ≥ 0.8 conditions using observed monthly averages during 2013–2016. Therefore, the photochemical and scattering terms are independent and can be used to describe the PAR absorbing and scattering roles related to the absorbing and scattering GLPs separately.
NOx (NO and NO2) are also important precursors of O3, their roles in O3 photochemistry in a subtropical evergreen broad-leaved forest (Dinghushan, Guangdong province, China) are studied using a similar method and empirical model as this study [48], but the difference is that it is in the UV region and without the consideration of BVOCs. The relative bias and NMSE of hourly ozone concentration are 6.82% and 0.01 for clear sky conditions (n = 113) and 11.30% and 0.02 for all sky conditions [48]. This can be as a reference for the representative of the photochemical term, the total energy absorption and use caused by all of the absorbing GLPs through OH radicals and H2O in CPRs. In this broad-leaved forest, O3 is more sensitive to its precursors (NO2, NO) than the other factors (UV, E, S/Q), more sensitive to NO2 than NO; the responses of O3 to the changes of driving factors (NO, NO2, UV, E, S/Q) are higher in summer than autumn and higher in clear skies than cloudy skies. If NOx is displayed in the O3–BVOCs empirical model, it would be progressive for understanding O3 photochemistry in the future when NOx is available. If more specific variables and their roles are needed to be simulated and investigated, e.g., BVOC emissions [34], O3 and BVOCs in this study, O3 and NOx [48], these variables can be picked out from the photochemical term and expressed explicitly, and let the roles of the other variables (not described explicitly) described in the photochemical term. When AVOCs are available in the future, their roles to O3 photochemistry can be expressed additionally and studied as a further application of this empirical model. The more variable variables are expressed, the deeper the understanding of O3–BVOCs–aerosols interactions achieved.
BVOC emissions vary with PAR and temperature. BVOCs react with OH radicals and AVOCs and produce new GLPs (e.g., cloud condensation nuclei, SOA, contributing to cloud formation). BVOCs play critical bridge roles to connect the atmospheric substances in gases, liquids and particles to the GLPs and solar energy. More than 30000 BVOCs are released from vegetation [72] and interact with other GLPs and PAR. The interactions between atmospheric substances (O3, NOx, VOCs, H2O, particles, etc.) and light are in many dimensions/directions and non-linear. Numerous atmospheric constituents change in three phases all the time, and many mechanisms associated with CPRs are still not clear, e.g., SOA formation from BVOCs, OH reactivity, but the one important point is that UV and VIS radiation is a very important energy source triggering the CPRs. No matter how the chemical compositions and reactions change, the energy associated with their main processes is a basis to drive their changes in the whole atmosphere. Therefore, the empirical energy method was selected to study the complicated interactions in O3–BVOCs–H2O–other GLPs–PAR. PAR and the energy interactions (Equations (1)–(4)) controlled the changes and interactions of O3, BVOCs and other GLPs (e.g., NOx, AVOCs, and stratospheric O3 through OH radicals). It is an advantage using the energy method because we pay attention to only their energy roles and do not need high or low concentrations of BVOCs and O3, the specific mechanisms discussed in the introduction, how the BVOCs and O3 correlated, etc. The actual roles with the changes of BVOCs and O3 under realistic atmospheric conditions are expressed by their related terms (Equations (1)–(4)). It saves much time in the calculations of the concentrations, chemical and photochemical reaction rates, and PAR utilization of GLPs. A similar O3 empirical model is developed, and similar positive and negative interactions are also found between O3 and BVOCs in a subtropical bamboo forest as in this study [66]. The empirical model of O3 is a specific model for one site currently. It is a beginning using energy method to deal with the complicated O3 and BVOCs photochemistry, and more studies need to be conducted in other forests. Vegetation under low or high accumulated O3 can lead to increased or decreased isoprene emission [20]. This feature in the vegetation is somewhat similar to that in the atmosphere (Section 3.2.1, Table 3 and Table 4). Therefore, they should be investigated together for the interaction mechanisms between O3 and isoprene and other BVOCs in the atmosphere and plants.

4.2. Implications of Sensitivity Analysis

Equations (3) and (4) represent multi-directional interactions in the natural atmosphere. The positive responses of O3 to the change of S/Q for considering isoprene and monoterpenes indicate that the harmful fine particulate matter (e.g., PM2.5, SOA) and O3 can be reduced simultaneously, especially at good light and clean atmospheric conditions. O3 responses to isoprene and monoterpenes are positive in most environmental conditions, indicating that it is critical to control or regulate high O3 formation from BVOC oxidation. Understanding how to regulate high BVOC emissions in cities is essential. It was found that wounded leaves and grass cutting enhanced BVOC emissions dramatically [73,74], and biomass burning results in high BVOC emissions and O3 in this subtropical plantation [34], so it is practical to reduce artificially enhanced BVOC emissions, e.g., changing plant-cutting to after 16:00, reducing biomass burning (e.g., straw) or planting trees and grass with no or low BVOC emissions for large cities so as to reduce high O3 and fine particle formation.
It is realized that PM2.5 and O3 are produced by VOCs reacting with OH radicals and other GLPs triggered by UV and VIS [74]:
AVOCs + BVOCs + OH+ NO2 + SO2 +UV + VIS…→new GLPs (O3, PM…)
China’s forests or tree and grass planting areas have been expanding fast. For example, the planting area in Beijing was 26 × 107 m2 in 2001 and increased by 6 × 107 m2 compared to 1995. Forests in Hebei province cover 32% of its land territory in 2017 compared to 26% in 2011. The forest area would be 1.14 × 1010 m2 and the coverage above 35% in Beijing–Tianjin–Hebei region in 2020. The increasing trees and grasses will enhance BVOC emissions and consume more NOx and SO2, then produce a large number of air pollutants (O3, SOA, etc.). Therefore, stricter control of NOx and SO2 emissions is highly recommended for polluted regions in China, such as the Beijing–Tianjin–Hebei region. In addition, AVOC emissions should also be controlled as suggested [22,74].
Satellite results indicate a greening pattern in 2000–2017 worldwide that is strikingly prominent in China, and China contributes 25% of the global net increase in leaf area [75]. To reach carbon peaking by 2030 and carbon neutrality by 2060, more trees and grasses will need to be planted in China. Considering the current situation and future changes, more stringent measures are much needed in emission control of human-induced BVOCs, AVOCs, NOx and SO2. High techniques to meet higher needs in emission control of AVOCs, NOx and SO2 in the near future are very necessary.
A similar O3 empirical model considering isoprene’s role is developed in a bamboo forest as does in this study [66]. The sensitivity shows that O3 responds a) negatively with the changes of isoprene at PAR (1.83 mol m−2) and low temperature (T = 13 °C) and (b) positively with the changes of isoprene at PAR (1.97 mol m−2) and low temperature (21 °C). It is deduced that O3 responds to isoprene positively at suitable air temperatures (21–32 °C), implying that both increased O3 and isoprene can be achieved by using the potential energy from outside of the O3–isoprene system (i.e., PAR), and negatively at high air temperature (>32 °C) and low air temperature (<13 °C), implying that an energy transport from O3 to isoprene (in O3–isoprene system) through an O3 decrease at states of very high and very low atmospheric system energy (solar global irradiance at the ground together with CpT; Cp is specific heat at constant pressure) when no more extra energy can be used from the outside of the O3–isoprene system. It means that air temperature is an important factor in the photochemistry of O3–isoprene. There were no outputs of simulated O3 considering isoprene or monoterpenes roles when PAR decreased by 100%, indicating that O3 decreased to zero and no more photochemical O3 production during the night (Figure 5 and Figure 6).

4.3. Improved Empirical Model of O3 Concentration and Validation

To improve the simulation of O3, the air mass used in the O3 term was multiplied by two(S/Q) for hourly S/Q > 0.7 (2 times S/Q in the development of the O3 model), considering the increase of the optical path caused by the GLP scattering in high GLP loads. Then, O3 concentrations were calculated using these modified empirical models for considering the roles of isoprene and monoterpenes, respectively. The new validation results are given briefly using the same data in Section 3.1. For considering isoprene’s role, the calculated O3 still overestimated the observation, but the relative bias decreased to 23% (51.2 vs. 41.8 ppb), RMSE was 18.9 ppb for single measurements (n = 104) and 15% (47.5 vs. 41.2 ppb) for the half-hourly variation. For considering monoterpene roles (n = 104), the calculated O3 underestimated the observed by 5% (39.8 vs. 41.6 ppb), RMSE was 16.0 ppb for single measurements and 8% (37.9 vs. 41.0 ppb) for the half-hourly variation. Comparing to the previous validation (Section 3.1), the modified empirical models improved the simulations of O3, especially for considering the isoprene situation.
It is necessary to understand chemical model performance over China and East China, i.e., their comparisons between the simulated and observed O3 concentrations. In short, the RMSE values for hourly O3 simulation are 21.6 and 41.0 ppb using the Models-3 Community Multi-scale Air Quality (CMAQ) [76], in the range of 9.4–20.1 ppb using nested air quality prediction modeling system (NAQPMS) [77], and 12.2–62.8 μg m−3 using the Regional Atmospheric Modeling System CMAQ [78]. RMSE values range from 9.9 to 28.1 ppb for daily O3 simulation using the eight regional Eulerian chemical transport models (CTMs) [79], and 10.0 to 32.7 ppb for annual O3 simulation using 14 state-of-the-art chemical transport models (CTMs) [80]. Ye et al. report that the RMSE value is 37 μg m−3 for O3 simulation in 7 days using the weather research forecasting coupled with chemistry model (WRF-Chem) [81]. In general, the performance of the empirical model of O3 is in agreement with these chemical models, though there are differences in time and space scales. With the more reliable BVOC emission fluxes and emission factors available in representative forests in China, model simulations of O3 and SOA would be improved in the future [34,68,82].
During the validation of the O3 empirical model, the observed ambient O3 concentrations, i.e., mixing ratios, were used for considering the individual role of isoprene or monoterpenes. For example, when considering the isoprene role, the empirical model expressed that O3 varies with isoprene explicitly and other chemical constituents (i.e., monoterpenes, NOx, SO2, etc.) non-explicitly in the photochemical term through the OH radicals and H2O (Section 4.1). It is similar to the situation when considering monoterpene roles.
To further investigate the performance of the energy method and the O3 empirical model, the emission fluxes of isoprene and monoterpenes were considered together as BVOC term to develop a new O3 empirical model (similar to Equation (4)) using the same dataset for considering the monoterpenes (n = 11). The corresponding coefficients and constant, as shown in Table 2, were 0.234, −0.274, 2.080, 1.555 and −1.089. R2 = 0.993, the mean and maximum of the relative bias were 7.43% and 17.87%, respectively. NMSE = 0.007, RMSE values were 3.32 ppb and 8.57%. It can be seen that similar simulations of O3 were also obtained in comparison with considering the roles of isoprene or monoterpenes, respectively (Section 3.1, Table 2). However, the specific roles that isoprene and monoterpenes play were mixed and changed in this new expression. To accurately understand their actual roles, it is better to describe the specific roles of isoprene and monoterpenes explicitly.

4.4. Relationships between Ozone and Its Influencing Factors

To understand the relationships between O3 and its influencing factors in natural atmospheric conditions, O3 concentrations during the REA measurements were calculated using the modified empirical models for considering the roles of isoprene (Section 4.3). The atmospheric GLPs (S/Q values) were divided into small groups with the interval of 0.1 between 0 and 1, and the corresponding sample points were 0, 9, 31, 40, 36, 37, 28, 17, and 76 for S/Q and other variables. All other parameters (PAR, air temperature, water vapor pressure, S/Q) were also divided into small groups with the same interval as the S/Q values.
The relationships between the calculated and observed O3 and each of the factors (PAR, water vapor (E), S/Q, air temperature (T), relative humidity (RH)) were non-linear (Figure 7). Generally, the calculated and observed O3 showed similar features and interactions with each of the influencing factors non-linearly, exhibiting that O3 was produced with the increases of PAR, T, RH, E and GLP loads when these factors are at low levels, and then O3 decreased with the increases of PAR, T, RH, E and GLPs after the driving factors reach their peaks. At which, the mean values were about 2.0 mol m−2 (=1110 μmol m−2 s−1) for PAR, 30 °C for T and 62% for RH, and 26 hPa for E, corresponding to a turning point S/Q = 0.5. It is evident that O3 and fine particles are produced simultaneously at low GLP loads (S/Q < 0.5, clean atmosphere), and O3 formation is inhibited at high GLP loads (S/Q > 0.5, polluted atmosphere). It should be emphasized that the production and destruction of O3 and fine particles depend on light and atmospheric conditions. A similar relationship (i.e., Figure 7) between measured O3 and relative humidity and O3 peaks at RH 50%–60% are also found in the Beijing–Tianjin–Hebei region, China, in 2014–2017 [83]. The features shown in Figure 7 are popular or not are needed to be investigated in other regions for better understanding the interactions of O3 and its influencing factors.
The relationships between PAR and S/Q were non-linear negatively (Figure 7), reflecting much PAR is attenuated with the increase of GLP loads, including the formation of aerosols, clouds and other GLPs.
The relationships between the estimated BVOC emissions (isoprene + monoterpenes) using the emission model of BVOC emissions [34] and S/Q was inverted U-shape and described by BVOCs = −13.924 × (S/Q)2 + 14.295 × (S/Q) + 0.3275 (R2 = 0.8813), and the turning point was also at S/Q = 0.55.
The relationships between hourly NOx concentrations (n = 1367) and S/Q were analyzed using the data measured on 12 June–16 October 2014. NO, NO2 and NOx decreased with the increases of S/Q at low S/Q (<0.75) and increased with the increase of S/Q at high S/Q (>0.75) (Figure 8). It reflects that the NOx is also important O3 precursors when GLPs are low, but the turning points of NOx (S/Q = 0.75) lag a little compared with other factors, implying an important mechanism that more NOx still participate in the CPRs to produce O3 even at low BVOC emissions, PAR, T and E, after S/Q > 0.55. It is the reason that the stricter reduction of NOx emissions is adopted. After the point of 0.75, no NOx reacted with BVOCs and other GLPs to produce O3, and NOx accumulated in the atmosphere, associating with the decrease of PAR, T, RH and E.
O3, BVOCs, water vapor and NOx vary with S/Q non-linearly, and O3, BVOCs and GLPs interact through the light. The GLP amounts are also important factors in controlling the interactions of O3–BVOCs–H2O–GLPs–radiation. The Sun provides the UV and visible radiation to trigger the GLPs taking part in homogeneous and heterogeneous reactions, and UV radiation play more important roles than visible radiation because of its high frequency; thus, much UV radiation is absorbed and utilized by atmospheric GLPs than visible radiation [13]. Water and water vapor are important constituents and sources of OH radicals in the UV and visible regions. Several studies report that organic aerosol (OA) contributes more than 50% of the total mass of fine particulate matter (PM) during haze events in China, including North China [84,85,86,87]; thus, the contributions from VOCs (especially BVOCs) to formation of O3 and fine particles will be more evident and important, in view of the growth of plants in China in the current and future [74,75]. The different S/Q levels represent the approximate equilibrium states of the interactions between the total atmospheric substances and solar radiation energy. The state at S/Q around 0.5 reflects the highest energy state for the atmospheric GLPs, i.e., highest PAR, sensible heat and latent heat, and an optimal interaction between GLPs–light.
To further understand the GLPs–solar radiation system, the changes of GLPs and solar radiation were analyzed when S/Q increased from 0.2 to 0.6: the BVOC emissions and O3 increased (2.6 mg m−2 h−1 and 6.0 ppb), along with the decreases of NO, NO2 and NOx (2.0, 0.6 and 2.6 ppb), water vapor (3.3 hPa), temperature (9.9 °C), humidity (18.2%), PAR and global solar radiation (632.7 μmol m−2 s−1 and 239.2 W m−2). It is obvious that BVOCs, NOx and water vapor contributed to the O3 and fine particle photochemical formation, alone with much PAR consumption. During new GLP production, the air temperature dropped, which is associated with the decreases of global solar radiation at the surface and the increase in GLPs (T = 25.09 × (S/Q)2 − 48.232 × (S/Q) + 46.907, R2 = 0.7733). When S/Q increased from 0.6 to 1.0, almost all variables decreased, including BVOC emissions (by 2.7 mg m−2 h−1), water vapor (2.6 hPa), temperature (5.4 °C), PAR and global solar radiation (740.2 μmol m−2 s−1, 427.3 W m−2), O3 (13.0 ppb), NO, NO2 and NOx (0.0, 0.5, 0.5 ppb); however, humidity increased (18.7%). It demonstrates that an increase in GLPs accumulated, which is associated with the increase in NOx, destruction or low production of O3 and low emissions of BVOCs, a great loss of solar radiation in the atmosphere and the drop of air temperature, implying a mechanism that the accumulation of air pollutants results in a more stable atmosphere. Therefore, BVOCs–other GLPs–solar radiation, i.e., GLPs–light, should be studied as a whole system.
Lee et al. [88] report that UK surface NO2 levels dropped by 42% during the COVID-19 lockdown, but O3 increased compared to previous years, which is attributed to the increased isoprene, UV and temperature. These observed facts provide evidence that the increase of BVOC emissions and UV (together with PAR) result in the production of O3 from BVOC oxidation. These results are in good agreement with the above analyses (Section 3.2 and 4.3). It should be emphasized that the interactions between O3 and its driving factors are very complicated, and the control strategies of O3 should consider the actual states of S/Q (GLPs)–solar radiation, especially in good light and atmospheric conditions and S/Q at low levels. Apart from the driving factors discussed above, other factors, e.g., UV, UV-A, climate change, CO2, warming, bidirectional exchange of BVOCs [89,90,91,92,93,94], influence BVOC emissions positively or negatively and are suggested to be studied together in the future.
Under UV radiation: the ozone reacts with H2O to produce OH radical, which reacts with most GLPs in the atmosphere, e.g., BVOCs and AVOCs, NOx and SO2, and then, they produce SOA [13,48,95], and the references therein and benefit cloud formation. Further, it influences the solar UV radiation balance in the atmosphere and on the ground. Similarly, under visible radiation, OH radicals produced from excited NO2* react with H2O and other OH sources [40,41,42,43,44], implying that H2O utilizes/transfers visible energy from NO2* and other GLPs (Section 2.3). In more detail, numerous GLPs absorb visible energy, e.g., glyoxal, CH3CO Radical, NO3 radical, OClO, CHOCHO, biacetyl, butenedial, NOCl and black carbon [96] and the references therein, and all of these absorbers react with OH radicals and H2O and transfer absorbed visible energy to other GLPs. Later, absorbers and non-absorbers take part in CPRs, exchange/consume visible energy and contribute to the formation of aerosols and clouds. Visible and UV light play key roles in the formation of SOA and clouds and then air motion but with differences because of their different electromagnetic frequency [97]. Figure 9 shows the mechanism of OH radical production through H2O and BVOC oxidation in CPRs, aerosols and cloud formation under UV and visible light, as well as the interactions in BVOC emissions and anthropogenic emissions of NOx, SO2, aerosol formation, solar UV and visible light. These interactions are in multiple directions and occur simultaneously in the natural atmosphere. Given the numerous limitations in measurements and simulations, together with large uncertainties and unknown mechanisms in OH [98], SOA production from BVOC’s reaction with O3, NO3 and OH + NOx in aqueous photochemistry [99] and from isoprene oxidation [100], missing OH sinks and unmeasured VOCs [46,96], the uncertainties in the application of the laboratory results (e.g., chamber) under controlled conditions to the natural atmosphere and thousands of BVOC compounds and their heterogeneous CPRs, it is impossible and a challenge to measure and simulate each chemical constituent and chemical and photochemical feature (Section 1) [28,72,73,96,97,98,99,100]. Therefore, it is an optional and practical method to study so complicated interactions and systems by grasping the key absorbing and scattering energy processes in the atmosphere. The rationality of this point of view is well proven by clear evidence that PAR energy drives the changes of all chemical compositions as well as all GLPs (e.g., O3, water vapor, S/Q, NOx, Figure 7 and Figure 8) and some results from O3 empirical models were in agreement with the studies of the chamber and chemical models.
In terrestrial vegetation regions, the enhancements of tropospheric column HCHO concentration is consistent with BVOC responses [100], and increased HCHO vertical column density (VCD), surface O3 and aerosol optical density (AOD) are contributed by BVOC oxidation [74], revealing that BVOCs play critical roles in the formation of HCHO, O3 and aerosols in a realistic atmosphere. According to the above discussion, it is better to combine the chemical models focusing on the specific compositions and empirical models focusing on energy use and distribution to study the mechanisms of O3–BVOCs–aerosols–solar radiation thoroughly.
In a short summary, the O3 empirical model show similar performance to chemical models, and several similar results, e.g., O3 and SOA formation, O3 response to its driving factors that obtained in the laboratory and using different chemical models. There is also agreement in the interactions between calculated and observed O3 and its driving factors. All of the above results indicate that the empirical model of O3 can be used to study the photochemical mechanisms of O3 and BVOCs. The energy method has a unique advantage, but the empirical model based on only the correlation between pure numbers may do not have. Visible and UV radiation provides an important energy source to the GLPs in the photochemical reactions in the whole atmosphere and is a critical connection to grasp and understand O3–BVOCs–aerosols–radiation interactions.

5. Conclusions

Based on the principle of PAR energy balance, the empirical models of O3 concentration for considering the roles of isoprene or monoterpenes were developed for a subtropical coniferous forest. The calculated O3 concentrations were in agreement with those observed with relative biases of 7.4% and 9.1% for considering isoprene and monoterpenes, respectively. Reasonable validation results of the O3 empirical models were obtained. O3 concentrations were more sensitive to the changes in PAR and E than S/Q and finally isoprene or monoterpenes. It implies that O3 is produced through CPRs. Most responses of O3 with its affecting factors were positively linear, except monoterpenes. The responses of O3 to isoprene were negative at high light and low atmospheric GLP loads or positive at low light and high atmospheric GLP loads. The air temperature also played a key role in the negative or positive responses. Moreover, the responses of O3 to monoterpenes were negative. Monoterpenes are easier to oxidize to fine particles than isoprene. The positive response of O3 to the change of S/Q indicates that O3 and fine particles are produced synchronously, and artificially enhanced BVOC emissions should be controlled.
In natural environmental conditions, the relationships between O3 and its driving factors were non-linear. O3 increased with the increase of PAR, temperature, humidity, water vapor and GLPs simultaneously during low GLP loads (S/Q < 0.55). S/Q at 0.55 is a turning point between positive and negative interactions for most variables (BVOCs, O3, T, RH, E, NOx and PAR). It is the highest atmospheric substances-energy state, i.e., the sufficient and optimal interaction in GLPs (O3–BVOCs–H2O)–light, representing the best conditions for O3 photochemistry and controlling the directions of the GLP changes. After this point, O3 decreased with the increase of PAR, temperature, humidity, water vapor and GLPs.
It is suggested to regulate human-induced BVOC emissions, and adopt a more stringent reduction standard in AVOCs, NOx and SO2 emissions, so as to reduce O3 and fine particle formation. The energy method deserves to be studied for a better understanding of the O3–BVOCs–H2O–GLPs–radiation system and interactions between the atmospheric substances and solar radiation.

Funding

This research was supported by the National Natural Science Foundation of China (grant no. 41275137), Dragon 4 and 5 projects (ID 32771 and 59013) and the European Union (EU) 7 framework programme MarcoPolo (grant no. 606953).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Data used in this research measured at Qianyanzhou Station are published on a Big Earth Data Platform for Three Poles, http://poles.tpdc.ac.cn/zh-hans/ (accessed on 13 March 2021). The Relaxed Eddy Accumulation system was provided by the National Center for Atmospheric Research, which is sponsored by the US National Science Foundation. The author thanks all the people for their great assistance, including Alex Guenther at the University of California, Andrew Turnipseed at 2B Technologies, Inc. Boulder, CO 80301, USA, James Greenberg and Tiffany Duhl at the National Center for Atmospheric Research, Boulder, CO 80307, USA; H.M., Wang, F.T., Yang, Q.K., Li, G.Z., Liu, L. Huang, Y.G., Wang, S.Y., Yin, J.D., Zou, J.Z., Zhang, Y.F., Huang, G.L. Zhu at Taihe County, Jiangxi province, and X.W. Wan and Y.M. Wu from the Institute of Atmospheric Physics, CAS. The author thanks the Qianyanzhou Experimental Station of Red Soil and Hilly Land (CAS) for providing meteorological and solar radiation data from January to May 2013. The author also thanks all reviewers for their constructive comments and beneficial suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Tignor, M.; Allen, S.K.; Boschung, J.; Midgley, P.M. (Eds.) Intergovernmental Panel on Climate Change (IPCC). Climate change 2013: The physical science basis. In Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar]
  2. Guenther, A.; Hewitt, C.N.; Erickson, D.; Fall, R.; Geron, C.; Graedel, T.; Harley, P.; Klinger, L.; Lerdau, M.; McKay, W.A.; et al. A global model of natural volatile organic compound emissions. J. Geophys. Res. Atmos. 1995, 100, 8873–8892. [Google Scholar] [CrossRef]
  3. Lathière, J.; Hauglustaine, D.; Friend, A.D.; De Noblet-Ducoudré, N.; Viovy, N.; Folberth, G.A. Impact of climate variability and land use changes on global biogenic volatile organic compound emissions. Atmos. Chem. Phys. Discuss. 2006, 6, 2129–2146. [Google Scholar] [CrossRef] [Green Version]
  4. Went, F.W. Blue Hazes in the Atmosphere. Nat. Cell Biol. 1960, 187, 641–643. [Google Scholar] [CrossRef]
  5. Atkinson, R. Atmospheric chemistry of VOCs and NOx. Atmos. Environ. 2000, 34, 2063–2101. [Google Scholar] [CrossRef]
  6. Claeys, M.; Graham, B.; Vas, G.; Wang, W.; Vermeylen, R.; Pashynska, V.; Cafmeyer, J.; Guyon, P.; Andreae, M.O.; Artaxo, P.; et al. Formation of Secondary Organic Aerosols Through Photooxidation of Isoprene. Science 2004, 303, 1173–1176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kanakidou, M.; Seinfeld, J.H.; Pandis, S.N.; Barnes, I.; Dentener, F.J.; Facchini, M.C.; Dingenen, R.V.; Ervens, B.; Nenes, A.; Nielsen, C.J.; et al. Organic aerosol and global climatemodelling: A review. Atmos. Chem. Phys. 2005, 5, 1053–1123. [Google Scholar] [CrossRef] [Green Version]
  8. Böge, O.; Miao, Y.; Plewka, A.; Herrmann, H. Formation of secondary organic particle phase compounds from isoprene gas-phase oxidation products: An aerosol chamber and field study. Atmos. Environ. 2006, 40, 2501–2509. [Google Scholar] [CrossRef]
  9. Wałaszek, K.; Kryza, M.; Werner, M. The role of precursor emissions on ground level ozone concentration during summer season in Poland. J. Atmos. Chem. 2017, 75, 181–204. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, K.; Yang, X.; Chen, D.; Gu, S.; Lu, Y.; Jiang, Q.; Wang, K.; Ou, Y.; Qian, Y.; Shao, P.; et al. Estimation of biogenic VOC emissions and their corresponding impact on ozone and secondary organic aerosol formation in China. Atmos. Res. 2020, 231, 104656. [Google Scholar] [CrossRef]
  11. Riipinen, I.; Yli-Juuti, T.; Pierce, J.R.; Petäjä, T.; Worsnop, D.R.; Kulmala, M.; Donahue, N.M. The contribution of organics to atmospheric nanoparticle growth. Nat. Geosci. 2012, 5, 453–458. [Google Scholar] [CrossRef]
  12. Wright, T.P.; Hader, J.D.; Mcmeeking, G.R.; Petters, M.D. High Relative Humidity as a Trigger for Widespread Release of Ice Nuclei. Aerosol Sci. Technol. 2014, 48, i. [Google Scholar] [CrossRef] [Green Version]
  13. Bai, J. UV extinction in the atmosphere and its spatial variation in North China. Atmos. Environ. 2017, 154, 318–330. [Google Scholar] [CrossRef]
  14. Thompson, A.M. Measuring and Modeling the Tropospheric Hydroxyl Radical (OH). J. Atmos. Sci. 1995, 52, 3315–3327. [Google Scholar] [CrossRef] [Green Version]
  15. Kanaya, Y.; Cao, R.; Akimoto, H.; Fukuda, M.; Komazaki, Y.; Yokouchi, Y.; Koike, M.; Tanimoto, H.; Takegawa, N.; Kondo, Y. Urban photochemistry in central Tokyo: 1. Observed and modeled OH and HO2 radical concentrations during the winter and summer of 2004. J. Geophys. Res. 2007, 112, D21312. [Google Scholar] [CrossRef]
  16. Praplan, A.P.; Pfannerstill, E.Y.; Williams, J.; Hellén, H. OH reactivity of the urban air in Helsinki, Finland, during winter. Atmos. Environ. 2017, 169, 150–161. [Google Scholar] [CrossRef]
  17. Di Carlo, P.; Brune, W.H.; Martinez, M.; Harder, H.; Lesher, R.; Ren, X.R.; Thornberry, T.; Carroll, M.A.; Young, V.; Shepson, P.B.; et al. Missing OH reactivity in a forest: Evidence for unknown biogenic VOCs. Science 2004, 304, 722–725. [Google Scholar] [CrossRef]
  18. Williamsa, J.; Keßela, S.U.; Nölschera, A.C.; Yang, Y.; Lee, Y.; Yáñez-Serrano, A.M.; Wolff, S.; Kesselmeier, J.; Klüpfel, T.; Lelieveld, J.; et al. Opposite OH reactivity and ozone cycles in the Amazon rainforest and megacity. Atmos. Environ. 2016, 125, 112–118. [Google Scholar] [CrossRef] [Green Version]
  19. Zannoni, N.; Gros, V.; Lanza, M.; Sarda, R.; Bonsang, B.; Kalogridis, C.; Preunkert, S.; Legrand, M.; Jambert, C.; Boissard, C.; et al. OH reactivity and concentrations of biogenic volatile organic compounds in a Mediterranean forest of downy oak trees. Atmos. Chem. Phys. Discuss. 2016, 16, 1619–1636. [Google Scholar] [CrossRef] [Green Version]
  20. Biesenthal, T.; Wu, Q.; Shepson, P.; Wiebe, H.; Anlauf, K.; Mackay, G. A study of relationships between isoprene, its oxidation products, and ozone, in the Lower Fraser Valley, BC. Atmos. Environ. 1997, 31, 2049–2058. [Google Scholar] [CrossRef]
  21. Tao, Z.; Williams, A.; Huang, H.-C.; Caughey, M.; Liang, X.-Z. Sensitivity of Surface Ozone Simulation to Cumulus Parameterization. J. Appl. Meteorol. Clim. 2008, 47, 1456–1466. [Google Scholar] [CrossRef]
  22. Li, N.; He, Q.; Greenberg, J.; Guenther, A.; Li, J.; Cao, J.; Wang, J.; Liao, H.; Wang, Q.; Zhang, Q. Impacts of biogenic and anthropogenic emissions on summertime ozone formation in the Guanzhong Basin, China. Atmos. Chem. Phys. Discuss. 2018, 18, 7489–7507. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, Y.; Li, L.; An, J.; Huang, L.; Yan, R.; Huang, C.; Wang, H.; Wang, Q.; Wang, M.; Zhang, W. Estimation of biogenic VOC emissions and its impact on ozone formation over the Yangtze River Delta region, China. Atmos. Environ. 2018, 186, 113–128. [Google Scholar] [CrossRef]
  24. Guenther, A.; Baugh, W.; Davis, K.; Hampton, G.; Harly, P.; Klinger, L.; Vierling, L.; Zimmerman, P. Isoprene fluxes measured by enclosure, relaxed eddy accumulation, surface layer gradient, mixed layer gradient, and mixed layer mass balance tech-niques. J. Geophys. Res. 1996, 101, 18555–18567. [Google Scholar] [CrossRef] [Green Version]
  25. Guenther, A. Biological and Chemical Diversity of Biogenic Volatile Organic Emissions into the Atmosphere. ISRN Atmos. Sci. 2013, 2013, 786290. [Google Scholar] [CrossRef]
  26. Turpin, B.J.; Saxena, P.; Andrews, E. Measuring and simulating particulate organics in the atmosphere: Problems and prospects. Atmos. Environ. 2000, 34, 2983–3013. [Google Scholar] [CrossRef]
  27. Goldstein, A.H.; Galbally, I.E. Known and Unexplored Organic Constituents in the Earth’s Atmosphere. Environ. Sci. Technol. 2007, 41, 1514–1521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Carlton, A.G.; Wiedinmyer, C.; Kroll, J.H. A review of Secondary Organic Aerosol (SOA) formation from isoprene. Atmos. Chem. Phys. Discuss. 2009, 9, 4987–5005. [Google Scholar] [CrossRef] [Green Version]
  29. Hu, X.M.; Fuentes, J.D.; Toohey, D.; Wang, D. Chemical processing within and above a loblolly pine forest in North Carolina, USA. J. Atmos. Chem. 2015, 72, 235–259. [Google Scholar] [CrossRef]
  30. Sadiq, M.; Tai, A.P.K.; Lombardozzi, D.; Martin, M.V. Effects of ozone–vegetation coupling on surface ozone air quality via biogeochemical and meteorological feedbacks. Atmos. Chem. Phys. 2017, 17, 3055–3066. [Google Scholar] [CrossRef] [Green Version]
  31. Grote, R.; Sharma, M.; Ghirardo, A.; Schnitzler, J.-P. A New Modeling Approach for Estimating Abiotic and Biotic Stress-Induced de novo Emissions of Biogenic Volatile Organic Compounds from Plants. Front. For. Glob. Chang. 2019, 2, 1–13. [Google Scholar] [CrossRef]
  32. Gong, C.; Lei, Y.; Ma, Y.; Yue, X.; Liao, H. Ozone–vegetation feedback through dry deposition and isoprene emissions in a global chemistry–carbon–climate model. Atmos. Chem. Phys. Discuss. 2020, 20, 3841–3857. [Google Scholar] [CrossRef] [Green Version]
  33. Liu, Y.; Yu, G.; Wen, X.; Wang, Y.; Song, X.; Li, J.; Sun, X.; Yang, F.; Chen, Y.; Liu, Q. Seasonal dynamics of CO2 fluxes from subtropical plantation coniferous ecosystem. Sci. China Ser. D Earth Sci. 2006, 49, 99–109. [Google Scholar] [CrossRef]
  34. Bai, J.; Guenther, A.; Turnipseed, A.; Duhl, T.; Greenberg, J. Seasonal and interannual variations in whole-ecosystem BVOC emissions from a subtropical plantation in China. Atmos. Environ. 2017, 161, 176–190. [Google Scholar] [CrossRef]
  35. Baker, B.; Guenther, A.; Greenberg, J.; Goldstein, A.; Fall, R. Canopy fluxes of 2-methyl-3-buten-2-ol over a ponderosa pine forest by relaxed eddy accumulation: Field data and model comparison. J. Geophys. Res. 1999, 104, 26107–26114. [Google Scholar] [CrossRef] [Green Version]
  36. Greenberg, J.P.; Guenther, A.; Harley, P.; Otter, L.; Veenendaal, E.M.; Hewitt, C.N.; James, A.E.; Owen, S.M. Eddy flux and leaf-level measurements of biogenic VOC emissions from mopane woodland of Botswana. J. Geophys. Res. 2003, 108, 8466. [Google Scholar] [CrossRef]
  37. Duhl, T.R.; Gochis, D.; Guenther, A.B.; Ferrenberg, S.; Pendall, E. Emissions of BVOC from lodgepole pine in response to mountain pine beetle attack in high and low mortality forest stands. Biogeosciences 2013, 10, 483–499. [Google Scholar] [CrossRef] [Green Version]
  38. Greenberg, J.P.; Guenther, A.; Zimmerman, P.R.; Baugh, W.; Geron, C.; Davis, K.; Helmig, D.; Klinger, L.F. Tethered balloon measurements of biogenic VOCs in the atmospheric boundary layer. Atmos. Environ. 1999, 33, 855–867. [Google Scholar] [CrossRef] [Green Version]
  39. Greenberg, J.P.; Guenther, A.B.; Madronich, S.; Baugh, W.; Ginoux, P.; Druilhet, A.; Delmas, R.; Delon, C. Biogenic volatile organic compound emissions in central Africa during the Experiment for the Regional Sources and Sinks of Oxidants (EX-PRESSO) biomass burning season. J. Geophys. Res. 1999, 104, 30659–30671. [Google Scholar] [CrossRef]
  40. Li, S.; Matthews, J.; Sinha, A. Atmospheric Hydroxyl Radical Production from Electronically Excited NO2 and H2O. Science 2008, 319, 1657–1660. [Google Scholar] [CrossRef]
  41. Ensberg, J.J.; Carreras-Sospedra, M.; Dabdub, D. Impacts of electronically photo-excited NO2 on air pollution in the South Coast Air Basin of California. Atmos. Chem. Phys. 2010, 10, 1171–1181. [Google Scholar] [CrossRef] [Green Version]
  42. Matthews, J.; Sinha, A.; Francisco, J.S. The importance of weak absorption features in promoting tropospheric radical pro-duction. Proc. Natl. Acad. Sci. USA 2005, 102, 7449–7452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Buckley, P.T.; Birks, J.W. Evaluation of visible-light photolysis of ozone-water cluster molecules as a source of atmospheric hydroxyl radical and hydrogen peroxide. Atmos. Environ. 1995, 29, 2409–2415. [Google Scholar] [CrossRef]
  44. Martins-Costa, M.T.; Anglada, J.M.; Francisco, J.S.; Ruiz-López, M.F. Photochemistry of SO2 at the Air-Water Interface: A Source of OH and HOSO Radicals. J. Am. Chem. Soc. 2018, 140, 12341–12344. [Google Scholar] [CrossRef] [PubMed]
  45. Fang, Q.; Han, J.; Jiang, J.; Chen, X.; Fang, W. The Conical Intersection Dominates the Generation of Tropospheric Hydroxyl Radicals from NO2 and H2O. J. Phys. Chem. A 2010, 114, 4601–4608. [Google Scholar] [CrossRef] [PubMed]
  46. Kim, S.; Karl, T.; Helmig, D.; Daly, R.; Rasmussen, R.; Guenther, A. Measurement of atmospheric sesquiterpenes by proton transfer reaction-mass spectrometry (PTR-MS). Atmos. Meas. Tech. 2009, 2, 99–112. [Google Scholar] [CrossRef] [Green Version]
  47. Bauerle, S.; Moortgat, G.K. Absorption cross-sections of HOCH2OOH vapor between 205 and 360 nm at 298 K. Chem. Phys. Lett. 1999, 309, 43–48. [Google Scholar] [CrossRef]
  48. Bai, J.H. Study on surface O3 chemistry and photochemistry by UV energy conservation. Atmos. Pollut. Res. 2010, 1, 118–127. [Google Scholar]
  49. Feng, Z.; Yuan, X.; Fares, S.; Loreto, F.; Li, P.; Hoshika, Y.; Paoletti, E. Isoprene is more affected by climate drivers than monoterpenes: A meta-analytic review on plant isoprenoid emissions. Plant Cell Environ. 2019, 42, 1939–1949. [Google Scholar] [CrossRef]
  50. Guenther, A.B.; Jiang, X.; Heald, C.L.; Sakulyanontvittaya, T.; Duhl, T.; Emmons, L.K.; Wang, X. The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): An extended and updated framework for modeling biogenic emissions. Geosci. Model Dev. 2012, 5, 1471–1492. [Google Scholar] [CrossRef] [Green Version]
  51. Griffin, R.J.; Cocker, D.R., III; Flagan, R.C.; Seinfeld, J.H. Organic aerosol formation from the oxidation of biogenic hydro-carbons. J. Geophys. Res. 1999, 104, 3555–3567. [Google Scholar] [CrossRef] [Green Version]
  52. Nishimura, H.; Shimadera, H.; Kondo, A.; Bao, H.; Shrestha, K.L.; Inoue, Y. Evaluation of light dependence of monoterpene emission and its effect on surface ozone concentration. Atmos. Environ. 2015, 104, 143–153. [Google Scholar] [CrossRef]
  53. Kroll, J.H.; Ng, N.L.; Murphy, S.M.; Flagan, R.C.; Seinfeld, J.H. Secondary organic aerosol formation from isoprene photooxidation under high-NOx conditions. Geophys. Res. Lett. 2005, 32, L18808. [Google Scholar] [CrossRef] [Green Version]
  54. Gomez, M.C.; Durana, N.; García, J.A.; de Blas, M.; de Camara, E.S.; García-Ruiz, E.; Gangoiti, G.; Torre-Pascual, E.; Iza, J. Long-term measurement of biogenic volatile organic compounds in a rural background area contribution to ozone formation. Agric. For. Meteorol. 2020, 287, 107979. [Google Scholar] [CrossRef]
  55. Lowe, P.R. An Approximating Polynomial for the Computation of Saturation Vapor Pressure. J. Appl. Meteorol. 1977, 16, 100–103. [Google Scholar] [CrossRef] [Green Version]
  56. Chang, J.C.; Hanna, S.R. Air quality model performance evaluation. Theor. Appl. Clim. 2004, 87, 167–196. [Google Scholar] [CrossRef]
  57. Tie, X.; Li, G.; Ying, Z.; Guenther, A.; Madronich, S. Biogenic emissions of isoprenoids and NO in China and comparison to anthropogenic emissions. Sci. Total Environ. 2006, 371, 238–251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Wang, Q.; Han, Z.; Wang, T.; Zhang, R. Impacts of biogenic emissions of VOC and NOx on tropospheric ozone during summertime in eastern China. Sci. Total Environ. 2008, 395, 41–49. [Google Scholar] [CrossRef]
  59. Tang, X.; Wang, Z.; Zhu, J.; Gbaguidi, A.E.; Wu, Q.; Li, J.; Zhu, T. Sensitivity of ozone to precursor emissions in urban Beijing with a Monte Carlo scheme. Atmos. Environ. 2010, 44, 3833–3842. [Google Scholar] [CrossRef]
  60. Geng, F.; Tie, X.; Guenther, A.; Li, G.; Cao, J.; Harley, P. Effect of isoprene emissions from major forests on ozone formation in the city of Shanghai, China. Atmos. Chem. Phys. Discuss. 2011, 11, 10449–10459. [Google Scholar] [CrossRef] [Green Version]
  61. Lee, K.-Y.; Kwak, K.-H.; Ryu, Y.-H.; Lee, S.-H.; Baik, J.-J. Impacts of biogenic isoprene emission on ozone air quality in the Seoul metropolitan area. Atmos. Environ. 2014, 96, 209–219. [Google Scholar] [CrossRef]
  62. Mo, Z.; Shao, M.; Wang, W.; Liu, Y.; Wang, M.; Lu, S. Evaluation of biogenic isoprene emissions and their contribution to ozone formation by ground-based measurements in Beijing, China. Sci. Total Environ. 2018, 627, 1485–1494. [Google Scholar] [CrossRef]
  63. Tagaris, E.; Sotiropoulou, R.; Gounaris, N.; Andronopoulos, S.; Vlachogiannis, D. Impact of biogenic emissions on ozone and fine particles over Europe: Comparing effects of temperature increase and a potential anthropogenic NOx emissions abatement strategy. Atmos. Environ. 2014, 98, 214–223. [Google Scholar] [CrossRef]
  64. Castell, N.; Stein, A.F.; Salvador, R.; Mantilla, E.; Millán, M. The impact of biogenic VOC emissions on photochemical ozone formation during a high ozone pollution episode in the Iberian Peninsula in the 2003 summer season. Atmos. Environ. 2016, 145, 326–337. [Google Scholar] [CrossRef]
  65. Miyoshi, A.S.; Hatakeyama, S.; Washida, N. OH radical-initiated photooxidation of isoprene: An estimate of global CO pro-duction. J. Geophys. Res. 1994, 99, 18779–18787. [Google Scholar] [CrossRef]
  66. Bai, J.H. Ozone and Its Affecting Factors in the Subtropical Bamboo Forest. Adv. Geosci. 2017, 07, 536–551. (In Chinese) [Google Scholar] [CrossRef]
  67. Xie, X.; Shao, M.; Liu, Y.; Lu, S.H.; Chang, C.C.; Chen, Z.M. Estimate of initial isoprene contribution to ozone formation po-tential in Beijing. Atmos. Environ. 2008, 42, 6000–6010. [Google Scholar] [CrossRef]
  68. Yu, F.; Liao, H. Simulation of the interannual variations of biogenic emissions of volatile organic compounds in China Impacts on tropospheric ozone and secondary organic aerosol. Atmos. Environ. 2012, 59, 170–185. [Google Scholar]
  69. Ruppert, L.; Becker, K.H. A product study of the OH radical-initiated oxidation of isoprene: Formation of C5-unsaturated diols. Atmos. Environ. 2000, 34, 1529–1542. [Google Scholar] [CrossRef]
  70. Horowitz, A.; Meller, R.; Moortgat, G.K. The UV–VIS absorption cross sections of the α-dicarbonyl compounds: Pyruvic acid, biacetyl and glyoxal. J. Photochem. Photobiol. A Chem. 2001, 146, 19–27. [Google Scholar] [CrossRef]
  71. Moosmüller, H.; Chakrabarty, R.; Arnott, W. Aerosol light absorption and its measurement: A review. J. Quant. Spectrosc. Radiat. Transf. 2009, 110, 844–878. [Google Scholar] [CrossRef]
  72. Harley, P.C. The Roles of Stomatal Conductance and Compound Volatility in Controlling the Emission of Volatile Organic Compounds from Leaves. In Biology, Controls and Models of Tree Volatile Organic Compound Emissions; Niinemets, Ü., Monson, R.K., Eds.; Tree Physiology, 5; Springer: Dordrecht, The Netherlands; New York, NY, USA, 2013; pp. 181–208. [Google Scholar]
  73. Fuentes, J.D.; Lerdau, M.; Atkinson, R.; Baldocchi, D.; Bottenheim, J.W.; Ciccioli, P.; Lamb, B.; Geron, C.; Gu, L.; Guenther, A.; et al. Biogenic Hydrocarbons in the Atmospheric Boundary Layer: A Review. Bull. Am. Meteorol. Soc. 2000, 81, 1537–1575. [Google Scholar] [CrossRef] [Green Version]
  74. Bai, J.; de Leeuw, G.; van der A, R.; De Smedt, I.; Theys, N.; Van Roozendael, M.; Sogacheva, L.; Chai, W. Variations and photochemical transformations of atmospheric constituents in North China. Atmos. Environ. 2018, 189, 213–226. [Google Scholar] [CrossRef]
  75. Chen, J.; Avise, J.; Guenther, A.; Wiedinmyer, C.; Salathe, E.; Jackson, R.B.; Lamb, B. Future land use and land cover influences on regional biogenic emissions and air quality in the United States. Atmos. Environ. 2009, 43, 5771–5780. [Google Scholar] [CrossRef] [Green Version]
  76. Zhang, M.; Uno, I.; Zhang, R.; Han, Z.; Wang, Z.; Pu, Y. Evaluation of the Models-3 Community Multi-scale Air Quality (CMAQ) modeling system with observations obtained during the TRACE-P experiment: Comparison of ozone and its related species. Atmos. Environ. 2006, 40, 4874–4882. [Google Scholar] [CrossRef]
  77. Yang, W.Y.; Chen, H.S.; Wang, W.D.; Wu, J.B.; Li, J.; Wang, Z.F.; Zheng, J.Y.; Chen, D.H. Modeling study of ozone source appor-tionment over the Pearl River Delta in 2015. Environ. Pollut. 2019, 253, 393–402. [Google Scholar] [CrossRef]
  78. Liu, H.; Zhang, M.; Han, X.; Li, J.; Chen, L. Episode analysis of regional contributions to tropospheric ozone in Beijing using a regional air quality model. Atmos. Environ. 2019, 199, 299–312. [Google Scholar] [CrossRef]
  79. Han, Z.; Sakurai, T.; Ueda, H.; Carmichael, G.; Streets, D.; Hayami, H.; Wang, Z.; Holloway, T.; Engardt, M.; Hozumi, Y.; et al. MICS-Asia II: Model intercomparison and evaluation of ozone and relevant species. Atmos. Environ. 2008, 42, 3491–3509. [Google Scholar] [CrossRef]
  80. Li, J.; Nagashima, T.; Kong, L.; Ge, B.; Yamaji, K.; Fu, J.S.; Wang, X.; Fan, Q.; Itahashi, S.; Lee, H.-J.; et al. Model evaluation and intercomparison of surface-level ozone and relevant species in East Asia in the context of MICS-Asia Phase III—Part 1: Overview. Atmos. Chem. Phys. Discuss. 2019, 19, 12993–13015. [Google Scholar] [CrossRef] [Green Version]
  81. Ye, L.; Wang, X.; Fan, S.; Chen, W.; Chang, M.; Zhou, S.; Wu, Z.; Fan, Q. Photochemical indicators of ozone sensitivity: Application in the Pearl River Delta, China. Front. Environ. Sci. Eng. 2016, 10, 15. [Google Scholar] [CrossRef]
  82. Wang, X.; Situ, S.; Chen, W.; Zheng, J.; Guenther, A.; Fan, Q.; Chang, M. Numerical model to quantify biogenic volatile organic compound emissions: The Pearl River Delta region as a case study. J. Environ. Sci. 2016, 46, 72–82. [Google Scholar] [CrossRef]
  83. Wang, M.; Zheng, Y.F.; Liu, Y.J.; Li, Q.P.; Ding, Y.H. Characteristics of ozone and its relationship with meteorological factors in Beijing-Tianjin-Hebei Region. Chin. Environ. Sci. 2019, 39, 2689–2698. [Google Scholar]
  84. Huang, R.-J.; Zhang, Y.; Bozzetti, C.; Ho, K.-F.; Cao, J.-J.; Han, Y.; Daellenbach, K.R.; Slowik, J.G.; Platt, S.M.; Canonaco, F.; et al. High secondary aerosol contribution to particulate pollution during haze events in China. Nat. Cell Biol. 2014, 514, 218–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Sun, Y.; Jiang, Q.; Xu, Y.; Ma, Y.; Zhang, Y.; Liu, X.; Li, W.; Wang, F.; Li, J.; Wang, P.; et al. Aerosol characterization over the North China Plain: Haze life cycle and biomass burning impacts in summer. J. Geophys. Res. Atmos. 2016, 121, 2508–2521. [Google Scholar] [CrossRef] [Green Version]
  86. An, Z.; Huang, R.-J.; Zhang, R.; Tie, X.; Li, G.; Cao, J.; Zhou, W.; Shi, Z.; Han, Y.; Gu, Z.; et al. Severe haze in northern China: A synergy of anthropogenic emissions and atmospheric processes. Proc. Natl. Acad. Sci. USA 2019, 116, 8657–8666. [Google Scholar] [CrossRef] [Green Version]
  87. Xu, W.; Xie, C.; Karnezi, E.; Zhang, Q.; Wang, J.; Pandis, S.N.; Ge, X.; Zhang, J.; An, J.; Wang, Q.; et al. Summertime aerosol volatility measurements in Beijing, China. Atmos. Chem. Phys. 2019, 19, 10205–10216. [Google Scholar] [CrossRef] [Green Version]
  88. Lee, J.D.; Drysdale, W.S.; Finch, D.P.; Wilde, S.E.; Palmer, P.I. UK surface NO2 levels dropped by 42% during the COVID-19 lockdown: Impact on surface O3. Atmos. Chem. Phys. 2020, 20, 15743–15759. [Google Scholar] [CrossRef]
  89. Pirjola, L. Effects of the Increased UV Radiation and Biogenic VOC Emissions on Ultrafine Sulphate Aerosol Formation. J. Aerosol Sci. 1999, 30, 355–367. [Google Scholar] [CrossRef]
  90. Pallozzi, E.; Fortunati, A.; Marino, G.; Loreto, F.; Agati, G.; Centritto, M. BVOC emission from Populus x canadensis saplings in response to acute UV-A radiation. Physiol. Plant. 2013, 148, 51–61. [Google Scholar] [CrossRef]
  91. Niinemets, Ü.; Fares, S.; Harley, P.; Jardine, K. Bidirectional exchange of biogenic volatiles with vegetation: Emission sources, reactions, breakdown and deposition. Plant Cell Environ. 2014, 37, 1790–1809. [Google Scholar] [CrossRef]
  92. Fraser, W.T.; Blei, E.; Fry, S.C.; Newman, M.F.; Reay, D.S.; Smith, K.A.; McLeod, A.R. Emission of methane, carbon monoxide, carbon dioxide and short-chain hydrocarbons from vegetation foliage under ultraviolet irradiation. Plant Cell Environ. 2015, 38, 980–989. [Google Scholar] [CrossRef] [Green Version]
  93. Guidolotti, G.; Rey, A.; Medori, M.; Calfapietra, C. Isoprenoids emission in Stipa tenacissima L.: Photosynthetic control and the effect of UV light. Environ. Pollut. 2016, 208, 336–344. [Google Scholar] [CrossRef]
  94. Holopainen, J.K.; Virjamo, V.; Ghimire, R.P.; Blande, J.D.; Julkunen-Tiitto, R.; Kivimäenpää, M. Climate Change Effects on Secondary Compounds of Forest Trees in the Northern Hemisphere. Front. Plant Sci. 2018, 9, 1445. [Google Scholar] [CrossRef] [Green Version]
  95. Fiore, A. No equatorial divide for a cleansing radical. Nat. Cell Biol. 2014, 513, 176–178. [Google Scholar] [CrossRef]
  96. Bai, J. Photosynthetically active radiation loss in the atmosphere in North China. Atmos. Pollut. Res. 2013, 4, 411–419. [Google Scholar] [CrossRef] [Green Version]
  97. Bai, J.; Zong, X. Global Solar Radiation Transfer and Its Loss in the Atmosphere. Appl. Sci. 2021, 11, 2651. [Google Scholar] [CrossRef]
  98. Patra, P.K.; Krol, M.C.; Montzka, S.A.; Arnold, T.; Atlas, E.L.; Lintner, B.R.; Stephens, B.B.; Xiang, B.; Elkins, J.W.; Fraser, P.J.; et al. Observational evidence for interhemispheric hydroxyl-radical parity. Nature 2014, 513, 219. [Google Scholar] [CrossRef]
  99. Romonosky, D.E.; Li, Y.; Shiraiwa, M.; Laskin, A.; Laskin, J.; Nizkorodov, S.A. Aqueous Photochemistry of Secondary Organic Aerosol of α-Pinene and α-Humulene Oxidized with Ozone, Hydroxyl Radical, and Nitrate Radical. J. Phys. Chem. A 2017, 121, 1298–1309. [Google Scholar] [CrossRef]
  100. Naimark, J.G.; Fiore, A.M.; Jin, X.; Wang, Y.; Klovenski, E.; Braneon, C. Evaluating Drought Responses of Surface Ozone Precursor Proxies: Variations with Land Cover Type, Precipitation, and Temperature. Geophys. Res. Lett. 2021, 48, e2020GL091520. [Google Scholar] [CrossRef]
Figure 1. (a) The simulated and observed O3 concentrations (O3 sim and O3 obs), with error bars indicating standard deviation of observed O3 (n = 12) (hereinafter, the same) in a subtropical Pinus plantation in China in 2013, and the relative biases (δ, %) (left, situation A, considering isoprene role). (b) Scatter plot of O3 concentrations calculated vs. observed (right, situation A).
Figure 1. (a) The simulated and observed O3 concentrations (O3 sim and O3 obs), with error bars indicating standard deviation of observed O3 (n = 12) (hereinafter, the same) in a subtropical Pinus plantation in China in 2013, and the relative biases (δ, %) (left, situation A, considering isoprene role). (b) Scatter plot of O3 concentrations calculated vs. observed (right, situation A).
Atmosphere 12 00711 g001
Figure 2. (a) Same as Figure 1a, but for considering monoterpene roles (situation B); (b) Same as Figure 1b, but for considering monoterpene roles (situation B).
Figure 2. (a) Same as Figure 1a, but for considering monoterpene roles (situation B); (b) Same as Figure 1b, but for considering monoterpene roles (situation B).
Atmosphere 12 00711 g002
Figure 3. The simulated and observed mean half-hourly O3 concentrations with error bars indicating 2 standard deviations of observed O3 (n = 104) in 2013 and the relative biases (δ, %) for situation A.
Figure 3. The simulated and observed mean half-hourly O3 concentrations with error bars indicating 2 standard deviations of observed O3 (n = 104) in 2013 and the relative biases (δ, %) for situation A.
Atmosphere 12 00711 g003
Figure 4. Same as Figure 3 but for situation B.
Figure 4. Same as Figure 3 but for situation B.
Atmosphere 12 00711 g004
Figure 5. O3 changing rates (%) with the change of one factor (isoprene, or PAR, E, S/Q for situations of (ad), respectively) and other factors kept at their original levels under realistic atmospheric conditions, REA and GRA denote the emission data measured by REA and gradient technique, respectively. REA12 denotes the isoprene emission data used in model development (O3 empirical model considering isoprene role).
Figure 5. O3 changing rates (%) with the change of one factor (isoprene, or PAR, E, S/Q for situations of (ad), respectively) and other factors kept at their original levels under realistic atmospheric conditions, REA and GRA denote the emission data measured by REA and gradient technique, respectively. REA12 denotes the isoprene emission data used in model development (O3 empirical model considering isoprene role).
Atmosphere 12 00711 g005
Figure 6. O3 changing rates (%) with the change of one factor (monoterpenes, or PAR, E, S/Q for situations of (ad), respectively) and other factors kept at their original levels under realistic atmospheric conditions. REA and GRA denote the emission data were measured by REA and gradient technique, respectively. REA11 denotes the monoterpene emission data used in model development (O3 empirical model considering monoterpene roles).
Figure 6. O3 changing rates (%) with the change of one factor (monoterpenes, or PAR, E, S/Q for situations of (ad), respectively) and other factors kept at their original levels under realistic atmospheric conditions. REA and GRA denote the emission data were measured by REA and gradient technique, respectively. REA11 denotes the monoterpene emission data used in model development (O3 empirical model considering monoterpene roles).
Atmosphere 12 00711 g006
Figure 7. The relationships between O3 concentration and its influencing factor, as well as between the observed PAR and S/Q, dark dot and circle denote the calculated and observed O3, respectively. The regression equations were obtained for a polynomial fit.
Figure 7. The relationships between O3 concentration and its influencing factor, as well as between the observed PAR and S/Q, dark dot and circle denote the calculated and observed O3, respectively. The regression equations were obtained for a polynomial fit.
Atmosphere 12 00711 g007
Figure 8. The relationships between NO, NO2, NOx and S/Q; the regression parameters obtained from a polynomial fit are for NO, NO2 and NOx from top to bottom, respectively.
Figure 8. The relationships between NO, NO2, NOx and S/Q; the regression parameters obtained from a polynomial fit are for NO, NO2 and NOx from top to bottom, respectively.
Atmosphere 12 00711 g008
Figure 9. Under solar UV and visible light, OH radical production and chemical and photochemical reactions between OH, BVOCs, SOA and O3 in gas, liquid and particle phases in the atmosphere. SOA contributes to clouds formation, and then clouds attenuate solar UV, visible and near-infrared radiation and influence the regional radiation balance and air motion. The emissions of BVOCs and their oxidation, along with the above associated potential effects, are interacted with and controlled by UV and visible energy.
Figure 9. Under solar UV and visible light, OH radical production and chemical and photochemical reactions between OH, BVOCs, SOA and O3 in gas, liquid and particle phases in the atmosphere. SOA contributes to clouds formation, and then clouds attenuate solar UV, visible and near-infrared radiation and influence the regional radiation balance and air motion. The emissions of BVOCs and their oxidation, along with the above associated potential effects, are interacted with and controlled by UV and visible energy.
Atmosphere 12 00711 g009
Table 1. Observation periods for O3 and BVOC emission flux measurements. The numbers in parentheses are the number of emission flux measurements.
Table 1. Observation periods for O3 and BVOC emission flux measurements. The numbers in parentheses are the number of emission flux measurements.
YearObservational Periods
201322 May–28 May (26)29 June–6 July (29)6 Aug.–13 Aug. (36)7 Sep.–11 Sep. (30)
201418 Jan.–19 Jan. (16)23 July–27 July (9)
201514 Jan.–19 Jan. (39)22 Apr.–30 Apr. (30)6 June–16 June (43)23 Aug.–4 Sep. (30)2 Nov.–7 Nov. (36)31 Dec. 2015–4 Jan. 2016 (37)
Table 2. The coefficients and constants, coefficient of determination (R2), average and maximum of the absolute relative bias (δavg, δmax (%)), normalized mean square error (NMSE) and standard deviations of calculated and observed fluxes (σcal and σobs) for the O3 empirical model considering isoprene and monoterpene roles (situation A and B), respectively.
Table 2. The coefficients and constants, coefficient of determination (R2), average and maximum of the absolute relative bias (δavg, δmax (%)), normalized mean square error (NMSE) and standard deviations of calculated and observed fluxes (σcal and σobs) for the O3 empirical model considering isoprene and monoterpene roles (situation A and B), respectively.
SituationA1A2A3A4A0R2δavgδmaxNMSEσcalσobs
A−0.038−0.5621.9760.226−0.1640.9967.416.50.0074.645.97
B1B2B3B4B0
B−0.049−0.0261.3990.219−0.1860.9939.128.00.0134.655.71
Table 3. O3 changing rates for REA12, REA and GRA (in % and ppb) caused by the changes of one factor at 20%, while other factors kept at their originally simulated levels under realistic atmospheric conditions.
Table 3. O3 changing rates for REA12, REA and GRA (in % and ppb) caused by the changes of one factor at 20%, while other factors kept at their originally simulated levels under realistic atmospheric conditions.
SituationISOPARES/Q
%ppb%ppb%ppb%ppb
REA 12−2.0−0.920.09.220.49.510.14.7
REA−0.1−0.312.87.017.09.712.07.6
GRA2.17.89.15.113.58.519.47.3
Table 4. Averages of the parameters measured in the establishment of the BVOC empirical model, by all REA and all gradient techniques.
Table 4. Averages of the parameters measured in the establishment of the BVOC empirical model, by all REA and all gradient techniques.
MeasurementPARTRHES/QO3
REA 122.7537.244.127.90.3546.9
REA1.8832.663.230.10.5940.5
GRA1.2621.876.622.10.7539.1
Table 5. Same as Table 3 but considering monoterpene (MT) roles.
Table 5. Same as Table 3 but considering monoterpene (MT) roles.
SituationMTPARES/Q
%ppb%ppb%ppb%ppb
REA 11−1.6−0.631.211.917.56.711.14.3
REA7.73.521.18.715.26.513.66.1
GRA11.018.16.315.15.519.66.7
Table 6. The coefficients and R2 of the fitting lines in responses of O3 with its driving factors considering isoprene’s role.
Table 6. The coefficients and R2 of the fitting lines in responses of O3 with its driving factors considering isoprene’s role.
Situation ISO PAR E S/Q
C1C0R2C1C0R2C1C0R2C1C0R2
REA12−0.9898.9700.99910.54−92.530.96610.06−101.20.9518.486−85.820.939
REA−0.0220.3430.2796.719−59.090.9978.689−85.90.9694.954−47.930.959
GRA0.876−6.4320.9534.686−41.480.9986.747−67.720.9544.648−44.080.986
Table 7. The coefficients and R2 of the fitting lines in responses of O3 with its driving factors considering monoterpene roles.
Table 7. The coefficients and R2 of the fitting lines in responses of O3 with its driving factors considering monoterpene roles.
Situation MT PAR E S/Q
C1C0R2C1C0R2C1C0R2C1C0R2
REA11−0.1572.5070.36716.79−146.40.9948.671−87.560.9425.143−48.790.985
REA−0.29511.300.93311.26−98.570.9957.480−75.930.9435.907−57.850.965
GRA−0.6907.2410.9059.515−83.740.9977.480−75.930.9438.563−86.000.946
Table 8. The O3 response rates with its driving factors considering isoprene’s role.
Table 8. The O3 response rates with its driving factors considering isoprene’s role.
SituationISOPARES/Q
REA12−0.110−0.114−0.099−0.099
REA−0.064−0.114−0.101−0.103
GRA−0.136−0.113−0.100−0.105
Table 9. The O3 response rates with its driving factors considering monoterpene roles.
Table 9. The O3 response rates with its driving factors considering monoterpene roles.
SituationMTPARES/Q
REA12−0.063−0.115−0.099−0.105
REA−0.026−0.114−0.099−0.102
GRA−0.095−0.114−0.099−0.100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bai, J. O3 Concentration and Its Relation with BVOC Emissions in a Subtropical Plantation. Atmosphere 2021, 12, 711. https://doi.org/10.3390/atmos12060711

AMA Style

Bai J. O3 Concentration and Its Relation with BVOC Emissions in a Subtropical Plantation. Atmosphere. 2021; 12(6):711. https://doi.org/10.3390/atmos12060711

Chicago/Turabian Style

Bai, Jianhui. 2021. "O3 Concentration and Its Relation with BVOC Emissions in a Subtropical Plantation" Atmosphere 12, no. 6: 711. https://doi.org/10.3390/atmos12060711

APA Style

Bai, J. (2021). O3 Concentration and Its Relation with BVOC Emissions in a Subtropical Plantation. Atmosphere, 12(6), 711. https://doi.org/10.3390/atmos12060711

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop