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., H
2O plays an energy use/transfer bridge role in OH radical formation through NO
2* + H
2O→OH, implying that H
2O utilizes the visible energy from NO
2*, which is similar to that in other situations (b–e,
Section 2.3). Similarly, H
2O utilizes UV energy from O
3 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 (O
3, NO
2, glyoxal, CH
3CO radical, NO
3 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 O
3 and isoprene or O
3 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, O
3 and scattering terms (as O
3 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 O
3 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 O
3 when O
3 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 O
3 and isoprene in Equation (1), or O
3 and monoterpenes in Equation (2), respectively. Therefore, Equations (3) and (4) were used to study the interactions between O
3 and isoprene and O
3 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 NO
2) are also important precursors of O
3, their roles in O
3 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 H
2O in CPRs. In this broad-leaved forest, O
3 is more sensitive to its precursors (NO
2, NO) than the other factors (UV, E, S/Q), more sensitive to NO
2 than NO; the responses of O
3 to the changes of driving factors (NO, NO
2, UV, E, S/Q) are higher in summer than autumn and higher in clear skies than cloudy skies. If NOx is displayed in the O
3–BVOCs empirical model, it would be progressive for understanding O
3 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], O
3 and BVOCs in this study, O
3 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 O
3 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 O
3–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 (O
3, NOx, VOCs, H
2O, 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 O
3–BVOCs–H
2O–other GLPs–PAR. PAR and the energy interactions (Equations (1)–(4)) controlled the changes and interactions of O
3, BVOCs and other GLPs (e.g., NOx, AVOCs, and stratospheric O
3 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 O
3, the specific mechanisms discussed in the introduction, how the BVOCs and O
3 correlated, etc. The actual roles with the changes of BVOCs and O
3 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 O
3 empirical model is developed, and similar positive and negative interactions are also found between O
3 and BVOCs in a subtropical bamboo forest as in this study [
66]. The empirical model of O
3 is a specific model for one site currently. It is a beginning using energy method to deal with the complicated O
3 and BVOCs photochemistry, and more studies need to be conducted in other forests. Vegetation under low or high accumulated O
3 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 O
3 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 O
3 to the change of S/Q for considering isoprene and monoterpenes indicate that the harmful fine particulate matter (e.g., PM
2.5, SOA) and O
3 can be reduced simultaneously, especially at good light and clean atmospheric conditions. O
3 responses to isoprene and monoterpenes are positive in most environmental conditions, indicating that it is critical to control or regulate high O
3 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 O
3 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 O
3 and fine particle formation.
It is realized that PM
2.5 and O
3 are produced by VOCs reacting with OH radicals and other GLPs triggered by UV and VIS [
74]:
China’s forests or tree and grass planting areas have been expanding fast. For example, the planting area in Beijing was 26 × 10
7 m
2 in 2001 and increased by 6 × 10
7 m
2 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 × 10
10 m
2 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 SO
2, then produce a large number of air pollutants (O
3, SOA, etc.). Therefore, stricter control of NOx and SO
2 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 SO
2. High techniques to meet higher needs in emission control of AVOCs, NOx and SO
2 in the near future are very necessary.
A similar O
3 empirical model considering isoprene’s role is developed in a bamboo forest as does in this study [
66]. The sensitivity shows that O
3 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 O
3 responds to isoprene positively at suitable air temperatures (21–32 °C), implying that both increased O
3 and isoprene can be achieved by using the potential energy from outside of the O
3–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 O
3 to isoprene (in O
3–isoprene system) through an O
3 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 O
3–isoprene system. It means that air temperature is an important factor in the photochemistry of O
3–isoprene. There were no outputs of simulated O
3 considering isoprene or monoterpenes roles when PAR decreased by 100%, indicating that O
3 decreased to zero and no more photochemical O
3 production during the night (
Figure 5 and
Figure 6).
4.3. Improved Empirical Model of O3 Concentration and Validation
To improve the simulation of O
3, the air mass used in the O
3 term was multiplied by two(S/Q) for hourly S/Q > 0.7 (2 times S/Q in the development of the O
3 model), considering the increase of the optical path caused by the GLP scattering in high GLP loads. Then, O
3 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 O
3 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 O
3 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 O
3, 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 O
3 concentrations. In short, the RMSE values for hourly O
3 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 O
3 simulation using the eight regional Eulerian chemical transport models (CTMs) [
79], and 10.0 to 32.7 ppb for annual O
3 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 O
3 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 O
3 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 O
3 and SOA would be improved in the future [
34,
68,
82].
During the validation of the O
3 empirical model, the observed ambient O
3 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 O
3 varies with isoprene explicitly and other chemical constituents (i.e., monoterpenes, NOx, SO
2, etc.) non-explicitly in the photochemical term through the OH radicals and H
2O (
Section 4.1). It is similar to the situation when considering monoterpene roles.
To further investigate the performance of the energy method and the O
3 empirical model, the emission fluxes of isoprene and monoterpenes were considered together as BVOC term to develop a new O
3 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. R
2 = 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 O
3 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 O
3 and its influencing factors in natural atmospheric conditions, O
3 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 O
3 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 O
3 showed similar features and interactions with each of the influencing factors non-linearly, exhibiting that O
3 was produced with the increases of PAR, T, RH, E and GLP loads when these factors are at low levels, and then O
3 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 O
3 and fine particles are produced simultaneously at low GLP loads (S/Q < 0.5, clean atmosphere), and O
3 formation is inhibited at high GLP loads (S/Q > 0.5, polluted atmosphere). It should be emphasized that the production and destruction of O
3 and fine particles depend on light and atmospheric conditions. A similar relationship (i.e.,
Figure 7) between measured O
3 and relative humidity and O
3 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 O
3 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 (R
2 = 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, NO
2 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 O
3 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 O
3 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 O
3, and NOx accumulated in the atmosphere, associating with the decrease of PAR, T, RH and E.
O
3, BVOCs, water vapor and NOx vary with S/Q non-linearly, and O
3, BVOCs and GLPs interact through the light. The GLP amounts are also important factors in controlling the interactions of O
3–BVOCs–H
2O–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 O
3 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 NO
2 levels dropped by 42% during the COVID-19 lockdown, but O
3 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 O
3 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 O
3 and its driving factors are very complicated, and the control strategies of O
3 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, CO
2, 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 H
2O to produce OH radical, which reacts with most GLPs in the atmosphere, e.g., BVOCs and AVOCs, NOx and SO
2, 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 NO
2* react with H
2O and other OH sources [
40,
41,
42,
43,
44], implying that H
2O utilizes/transfers visible energy from NO
2* and other GLPs (
Section 2.3). In more detail, numerous GLPs absorb visible energy, e.g., glyoxal, CH
3CO Radical, NO
3 radical, OClO, CHOCHO, biacetyl, butenedial, NOCl and black carbon [
96] and the references therein, and all of these absorbers react with OH radicals and H
2O 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 H
2O 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, SO
2, 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 O
3, NO
3 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., O
3, water vapor, S/Q, NOx,
Figure 7 and
Figure 8) and some results from O
3 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 O
3 and aerosol optical density (AOD) are contributed by BVOC oxidation [
74], revealing that BVOCs play critical roles in the formation of HCHO, O
3 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 O
3–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.