Characteristics and chemical reactivity of biogenic volatile organic compounds from dominant forest species in the Jing-Jin-Ji area,China

Ying Lin,Xiaoxiu Lun*,Wei Tang,Zhongzhi Zhang,Xiaoxi Jing,Chong Fan and Qiang Wang*

Abstract

Keywords: Biogenic volatile organic compounds (BVOCs), Isoprene, Monoterpenes,Jing-Jin-Ji area,Spatiotemporal characteristics, Chemical reactivity

Background

Biogenic volatile organic compounds (BVOCs) are low boiling point compounds commonly synthesized by secondary metabolic pathways in plants. Many vascular plants can discharge BVOCs into the atmosphere (Loreto and Schnitzler 2010). Forest is one of the primary sources that emit BVOCs, which occupies about 70% of the total BVOCs amounts from vegetation. It was estimated that the annual emission of BVOCs in the world was about 106Gg C·year-1(Guenther et al. 2012), accounting for more than 90% of the total non-methane volatile organic compounds (NMVOCs) emission on the ground and far exceeding the anthropogenic compounds(Guenther et al. 1995). Isoprene (the simplest 5-carbon isoprenoid, C5H8, 2-methyl 1,3-butadiene) represents the highest emission component (Atkinson and Arey 2003)with an approximately 50% of the total annual global emissions of BVOCs in around 412–601 Tg C·year-1(Guenther et al. 2012). Monoterpenes are the 10-carbon isoprenoids that account for about 15% (32–157 Tg C·year-1) of global BVOCs emission (Guenther et al.2012). Both isoprene and monoterpenes are synthesized by the MEP pathway (Loreto and Schnitzler 2010).BVOCs are usually formed constitutively or after stress induction. Those components can improve plant tolerance towards abiotic stressors such as high temperature,oxidative stress and biotic stressors (e.g. competing plants and microorganisms) (Loreto and Schnitzler 2010;Filella et al. 2013).

BVOCs play an essential role in tropospheric atmospheric chemical reactions (Sartelet et al. 2012; Kulmala et al. 2013). BVOCs are the main precursors to form tropospheric ozone and atmospheric aerosols, promoting the formation of secondary pollutants such as peroxyacetyl nitrate (PANs), secondary organic aerosols (SOA),particulate matter (PM), aldehydes and ketones (Claeys et al. 2004; Laothawornkitkul et al. 2009). In particular,the formation of ozone occurs when the isoprene is dissociated and react with NOx(Fehsenfeld et al. 1992)while the formation medium of secondary organic aerosols (Claeys et al. 2004) is contributed by monoterpenes and sesquiterpenes via cloud condensation nuclei, hence affecting the local or global climate.

The Jing-Jin-Ji area locating on the North China Plain represents the core of north China and the most developed city cluster in China. With the rapid development of the economy and the acceleration of urbanization, the issue of air pollution requires immediate attention. The Jing-Jin-Ji area has been plagued by severe photochemical pollution and haze for many years (Tang et al. 2009;Han et al. 2013). It is known that the accumulation of particulate matter, primarily fine particulate matter PM2.5is the main factor that causes the haze (Li et al.2013; Zhai et al. 2016; Hsu et al. 2017; Zhao et al. 2013)showed that the hourly average concentration of PM2.5in Beijing reached 318 μg·m-3on hazy days. Several cities in Hebei province suffered more severe air pollution(Wang et al. 2012; Li et al. 2016a) used the PM2.5monitoring data of 161 cities to analyze the PM2.5pollution in mainland China. The results showed that the Jing-Jin-Ji and its surrounding areas were heavily polluted and ranked at one of China’s current four smoggy regions.Zhao et al. (2020) demonstrated that the average annual PM2.5concentration in the Jing-Jin-Ji area decreased by 8.66 μg·m-3from 2014 to 2018, indicating the population that exposed to high PM2.5concentration was decreasing. However, the average annual PM2.5concentration value was still far from the national standard limit (35 μg·m-3). Also, the ozone is another important pollutant that plagues the urban ambient air quality after PM2.5(Meng et al. 2017; Chen et al. 2013) reported that the ozone mixing ratio in the Jing-Jin-Ji area was very high, thus causing strong photochemical reactions from May to September. Wang et al. (2017) pointed out that the ozone concentration had exceeded the standard by 100–200% in the Jing-Jin-Ji area. During 2014–2018,the ozone concentration in the Jing-Jin-Ji area showed an upward trend with an average positive annual level of 4.90 μg·m-3as reported by Zhao et al. (2020). Some studies indicated that the contribution of BVOCs such as monoterpenes and sesquiterpenes to SOA formation was substantial (Steinbrecher et al. 2009; Aksoyoglu et al. 2011; Ghirardo et al. 2016) found that the contribution of BVOCs released by vegetation to SOA generation in Beijing was increased by two-fold within 2005–2010. Carlo et al. (2004) reported that the isoprene is highly reactive with hydroxyl radical (·OH) than most anthropogenic volatile organic compounds (AVOCs).Geng et al. (2011) suggested that BVOCs can contribute to the surface ozone concentrations. The above literatures show that the role of BVOCs in ozone formation cannot be ignored in the Jing-Jin-Ji area especially during summer (Xie et al. 2008; Ran et al. 2011), which would have led to severe impacts on human health, economic development, ecological environment, and climate change (Pierre et al. 2017). Considering the particular geographical location, diverse vegetation composition, and the enormous influence of BVOCs on ozone and secondary aerosols, it is crucial to clarify the BVOCs emission emitted from dominant forest species in the Jing-Jin-Ji area.

At present, the estimation of BVOCs emission from plants have been conducted through various studies using different methods, including models, land cover and meteorological data. Different results of the Jing-Jin-Ji area were obtained from different models and input parameters at home and abroad (Klinger et al. 2002;Song et al. 2012; Li et al. 2016b). Most of the previous estimations were conducted by grouping the forest species into several plant functional types (PFT) and then using the global average emission rates and biomass density of each vegetation type. However, such estimation was unsatisfactory because the differences between plant species and regions were ignored (Xia and Xiao 2019). Besides, most of the studies employed the recommended values of emission factors in the literature directly or selected the measured values of adjacent areas without the measured local primary data; hence the final results may not be representative (Song et al. 2012;Zhang et al. 2018). Different forest species or regions with different environments will cause different emission components and release rates of BVOCs (Owen et al.2002). The difference between the forest species could be the primary factor to determine the BVOCs emission,and different geographical locations own disparate environmental conditions (temperature, PAR) that will also distinguish the emissions eventually. Therefore, the precise emission rates and fluxes of BVOCs from specific forest species in different regions should be obtained before the evaluation of the spatiotemporal effects and interactions of BVOCs on the atmospheric environment.Until present, there has been limited evidence that supports the systematic comparison of the BVOCs measured emission rates of dominant forest species in the Jing-Jin-Ji area. According to the order of volume from high to low, we have selected Platycladus orientalis,

Quercus variabilis, Betula platyphylla, Populus tomentosa, Pinus tabuliformis, Robinia pseudoacacia, Ulmus pumila, Salix babylonica and Larix gmelinii as the dominant forest species in the Jing-Jin-Ji area. Based on the field vegetation investigation and measured forest species emission rates, this study used the algorithm in G93 model (Guenther et al. 1993) to estimate the BVOCs emission from dominant forest species in the Jing-Jin-Ji area to establish a more localized emission inventory and eliminate its uncertainty further, analyzed their spatiotemporal distribution characteristics and chemical reactivity to provide a theoretical basis for future air pollution prevention and control measures.

Materials and methods

Site description

Jing-Jin-Ji area is located in the North China Plain(Fig. 1), including Beijing, Tianjin, and Hebei (36°05′–42°40′ N, 113°27′–119°50′ E), with a total area of about 21.72 ha. This area is surrounded by Bohai Sea in the east, Tai-hang Mountains in the West and Yanshan Mountains in the north, 735 km from north to south,and 576 km from east to west, with various terrains and warm temperate continental monsoon climate. The annual average temperature is ranged from 0 to 13oC with an average yearly precipitation of 300–800 mm. There are various vegetation types in this area including forest,shrub, grassland, and so on. The sampling sites were selected in the forest parks with abundant plant resources and the collection of the air samples were performed on bright days.

Fig.1 The site locations of the Jing-Jin-Ji area

The sampling sites (solid circles) identified as follows:(1) Saihanba National Forest Park; (2) Heilongshan National Forest Park; (3) Labagou Primeval Forest Park; (4)Baicaowa National Forest Park; (5) Wuling Mountain Scenic Spot; (6) Yunmengshan National Forest Park; (7)Huangyangshan Forest Park; (8) Tianjin Jiulongshan National Forest Park; (9) Xishan National Forest Park; (10)Xiaolongmen National Forest Park; (11) Baihuashan National Nature Reserve; (12) Shangfangshan National Forest Park;(13) Wuyuezhai Scenic Spot.

Estimation model of BVOCs emission

Guenther series models are widely used to estimate BVOCs emission. Guenther combined the latest experimental data to deduce the G93 algorithm in 1993(Guenther et al.1993).The model was used to normalize the emission rate of BVOCs under various environmental conditions, including T=303 K and PAR=1000 μmol·m-2·s-1. The BVOCs emitted by dominant forest species were classified as isoprene, monoterpenes, and other VOCs (alcohols, aldehydes, ketones, esters, organic acids, low carbon alkanes and alkenes) in this paper. According to the model of light and temperature effect proposed by G93, the BVOCs emission classified by forest species were estimated respectively. Specific formulas are as follows:

Isoprene:

where, R is a gas constant (8.314 J·K-1·mol-1); Ts is the leaf temperature of the standard state (303 K). CT1(95,000 J·mol-1), CT2(230,000 J·mol-1) and TM(314 K) are empirical constants.

The monoterpenes and other VOCs emission rate in algorithm G93 can be calculated using Eqs.4 and 5:

where, M is the emission rate of monoterpenes and other VOCs (in C, μg·g-1·h-1) under a certain temperature T (K); MSis the standard emission rate of monoterpenes and other VOCs (in C, μg·g-1·h-1) under standard condition (Ts=303 K); CTMis the correction factor for the temperature of monoterpenes and other VOCs; β (0.09 K-1) is empirical constant.

Experimental results showed that the emission rate of isoprene is mainly controlled by leaf temperature and PAR. However, the main factor affecting monoterpenes and other VOCs emission by plants is temperature(Guenther et al. 1993). Therefore, the emission estimation method of isoprene is:

where, EISOPis the emission of isoprene (in C, μg·h-1);EMONOand EOVOCare the emissions of monoterpenes and other VOCs (in C, μg·h-1), respectively; and B is the leaf biomass (dry weight, g) of each species. The required meteorological data (temperature and PAR of every month in the Jing-Jin-Ji area in 2017) can be retrieved from MODIS satellite product data published on the NASA website (https://modis.ornl.gov/data.html).

Determination of model parameters

Emission rates

Air samples were collected using a dynamic headspace method (Jing et al. 2020). The sampling system consists of a transparent Teflon film sampling bag (10 L), an atmospheric sampler (LaoDong QC-1 S, Beijing Municipal Institute of Labor Protection, China), two drying towers(filled with activated carbon particles that have been pre-dried for more than 5 h and allochroic silica gel, respectively), an ozone removal column (Cleanert KI:1.4 g/2.5 mL,Agela Technologies,China) and an adsorption tube (tube type: stainless steel tube, Camsco company, USA, filled with Carbograph 2 (60/80),Carbograph 1 (40/60) and Carbosieve SIII (60/80) adsorbents). All parts were connected with polytetrafluoroethylene (PTFE) tubes (Fig. 2). Before sampling, the adsorption tubes were purged and activated at 270 °C for 2h under high-purity N2(purity: 99.99%), then stored under cold storage at 4 °C. When sampling, the branches with healthy foliage were selected and enclosed in the sampling bag. First, the air in the sampling bag was extracted, and then the air that has passed through the two drying towers and the ozone adsorption column was pumped in. Finally, the adsorption tube was connected to form a closed loop. The air samples were collected at a flow rate of 150 mL·min-1for 1 h. The temperature and relative humidity in the sampling bag were measured by a hand meteorometer instrument(NK4500, Kestrel, USA). The PAR in the sampling environment was measured using a light quantum meter(3415FQF, Spectrum, USA). Table 1 shows the sampling dates and averaged environmental conditions during the sampling. After collection, the adsorption tube was stored in the refrigerator at 4°C and analyzed within one week.

Fig.2 Flow chart of dynamic headspace sampling

The air samples were desorbed on PE TurboMatrix(650ATD-Clarus600) and analyzed by thermal desorption gas chromatography-mass spectrometry (GC-MS, Agilent 6890, USA). Most of the organics in the adsorption tube were released after 5 min in a thermal resolver at 260 °C,and then the substances were adsorbed into the cold trap(–25°C).The substances condensed in the cold trap were quickly heated from 260 to 300 °C at a heating rate of 40 °C·s-1then transferred to the chromatograph for further separation and analysis. The column used was DB-5MS (column height of 30 m,inner diameter of 0.25 mm,pore diameter of 0.25 μm) and the carrier gas used was helium (1.0 mL·min-1). The heating process was divided into three stages at 40°C,160 and 270°C that maintained at 2, 2 and 3 min, respectively, with a heating rate of 4 °C·min-1. The mass spectrometer adopts an electron bombardment ionization mode with energy of 70 eV and a scanning atomic mass range of 30–500 amu. The compounds were retrieved from the database of National Institute of Standards and Technology (NIST) based on their retention times and specific charge and quantified using the standard external method. The standard gases used were as follows: Photochemical Assessment Monitoring Stations (PAMS) (Spectra/Linde: 57) and n-hexane, isoprene, α-pinene, β-pinene, α-phellandrene, 3-carene, βmyrcene,α-terpinene,limonene,γ-terpinene,and ocimene made by the National Institute of Metrology, China. The TCT-GC-MS instrument has been calibrated using standard gases. The recovery rate of standard addition was between 90 and 110% and coefficients of variance were less than 5%.

The emission rates of dominant forest species in the Jing-Jin-Ji area were respectively calculated according to the following Equation:

where, ER is the chemical substance emission rate of each forest species (in C, μg·g-1·h-1); m is the mass ofthe chemical substance in adsorption tube (in C, μg); t is the sampling time (h); M is the leaf biomass of forest species in the sapmling bag (dry weight, g).

Table 1 Sampling dates and averaged environmental conditions during the samplings in 2017

Leaf biomass calculation

The leaf biomass data in the Jing-Jin-Ji area were obtained using the method of volume and biomass conversion. The statistical method of subdividing species enables them to correspond to more appropriate emission factors and biomass, by considering the difference of the same species in different regions. Based on the volumes obtained from the national forest resource inventory, the leaf biomass of the dominant forest species was calculated as follows:

where B is leaf biomass of forest species (dry weight, g);V is tree volume (m3); DTis the basic density of tree trunk (the ratio of absolute dry wood mass to raw wood volume); PTand PLare the proportion of stem and leaf in the total biomass of tree layer respectively. The data of PT, PLand DTwere retrieved from literature (Wang et al. 2001). In this study, the proportion of trunk density, leaf and stem biomass to the total biomass of dominant forest species in the Jing-Jin-Ji area were presented in Table A1. Tree volumes data were mainly based on national forest resources survey data, checked and supplemented by forestry network and field survey.

Methods of chemical activity evaluation

Different BVOCs components have different chemical compositions and physical properties, thus the atmospheric chemical reaction capacity is different (Goldan et al. 2004). The chemical reactivity of different components and their ability to generate ozone can provide reference for the control measures of BVOCs. This study adopted two methods namely the maximum incremental reactivity (MIR) method and the ·OH reaction rate(LOH) method, to comprehensively analyze the chemical activity of the BVOCs components of dominant forest species in the Jing-Jin-Ji area.

Ozone formation potential (OFP) was used to evaluate the potential release of BVOCs into the atmosphere under optimal reaction conditions for ozone generation and measure the reactivity of different BVOCs components (Carter 1991). The calculation of OFP is shown in Eq. 10:

where OFPiis the ozone generation potential of each component (μg·m-3); MIRiis the maximum incremental response factor of each component ((g O3)/(g VOCs))(Carter 2008); Ciis the mass concentration of each component (μg·m-3).

The chemical reaction of the troposphere during the daytime is mainly OH radicals (·OH). The volatile organic compounds will first react with ·OH, and then react with O2and NOxunder light condition to generate new free radicals to initiate the chain reaction. The first reaction is a key step that determines the rate of atmospheric photochemical reaction chain, so the consumption rate of ·OH can be used to evaluate the photochemical activity of BVOCs. The calculation is shown in Eq. 11:

where LiOHis the consumption rate of each component to atmospheric ·OH (s-1); KiOHis the reaction rate constant between each component and atmospheric ·OH(cm3·molecule-1·s-1) (https://kinetics.nist.gov/kinetics/index.jsp); [BVOCs]iis the atmospheric molecular concentration of every component (molecule·cm-3).

Results and discussion

Emission budgets and compositions of BVOCs

According to Eqs. 6 and 7, based on the measured standardized emission factors of each forest species, leaf biomass data obtained from the forest resource data and field survey, and the annual meteorological data through the NASA website,the annual BVOCs emission of every component of each forest species in the Jing-Jin-Ji area was estimated. BVOCs emitted from dominant forest species were divided into isoprene, monoterpenes (including αpinene, β-pinene, β-myrcene, limonene, 3-carene and so on), and other VOCs (OVOCs) (Table A2). The normalized emission rates of BVOCs in dominant forest species in the Jing-Jin-Ji area are tabulated in Table 2. The total annual emission of BVOCs from dominant forest species in the Jing-Jin-Ji area was estimated to be 70.8 Gg C·year-1, including 40.5% (28.7 Gg C·year-1) of isoprene,36.0% (25.5 Gg C·year-1) of monoterpenes and 23.4%(16.6 Gg C·year-1)of other VOCs,respectively.

Different forest species can produce various components of BVOCs. Figure 3 shows each component and proportion of BVOCs emitted by nine dominant forest species in the Jing-Jin-Ji area. Isoprene was emitted as the main component from broadleaf trees such as Betula platyphylla, Quercus variabilis, and Salix babylonica.Aside from isoprene, some of the broadleaf trees (i.e.Populus tomentosa, Robinia pseudoacacia, Ulmus pumila) and most coniferous trees like Platycladus orientalis, Pinus tabuliformis, and Larix gmelinii released monoterpenes such as α-pinene, β-pinene, β-myrcene and limonene. As shown in Fig. A1, the isoprene emission was mainly detected from Quercus variabilis andPopulus tomentosa that recording at 41.2 % and 31.1%,respectively. On the other hand, Platycladus orientalis and Pinus tabuliformis are the significant emitters of monoterpenes, providing 56.1 % and 36.0%, respectively.Therefore, it was revealed that these four forest species are the main contributors to BVOCs emission in the Jing-Jin-Ji area: 19.3 Gg C·year-1(27.2%) from Platycladus orientalis, 14.0 Gg C·year-1(19.8%) from Quercus variabilis, 13.4 Gg C·year-1(18.9%) from Populus tomentosa and 12.2 Gg C·year-1(17.2%) from of Pinus tabuliformis (Table A2).

Table 2 Normalized emission rates of biogenic volatile organic compounds (BVOCs)in dominant forest species in the Jing-Jin-Ji area (μg·g-1·h-1)

Monthly and seasonal variations

As shown in Eqs. 6 and 7, by multiplying the measured standardized emission factors of each forest species,monthly leaf biomass data obtained from the forest resource data and field survey, and the monthly meteorological data through the NASA website in 2017, the monthly and seasonal BVOCs emission of every component of each forest species in the Jing-Jin-Ji area were estimated. Table A3 and Fig. 4 indicate that the BVOCs emission and composition in the Jing-Jin-Ji area demonstrate significant variations for both monthly and seasonal. A distinct unimodal change with the month for the emissions of BVOCs was observed.

Fig.3 The BVOCs emission(106 g C·year-1) (a) and proportion of each component (%)(b) from dominant forest species in the Jing-Jin-Ji area(BP:Betula platyphylla; QV:Quercus variabilis; UP:Ulmus pumila; PoT:Populus tomentosa; RP:Robinia pseudoacacia; PO:Platycladus orientalis;PiT:Pinus tabuliformis; SB:Salix babylonica;LG: Larix gmelinii)

Fig.4 The monthly BVOCs emission from dominant forest species in the Jing-Jin-Ji area(106 g C)

The total BVOCs emission was peaked in July with a total amount of 15.8 Gg C·year-1, while recorded a minimum value in January of 0.4 Gg C·year-1with two orders of magnitude difference. In January, the emission of isoprene from vegetation hit the lowest point at 1.9×106g C·year-1due to low temperature, less sunlight,and the limited leaves of biomass. In contrast, the emission of isoprene achieved the highest point at 7.4 Gg C·year-1in July when the air temperature was high with presence of strong sunlight, and abundant leaves of biomass. Nevertheless, the emission of isoprene fluctuated in other months between the maximum and minimum values obtained. The same trend was obtained from the emission of monoterpenes and other VOCs in which the highest emission was recorded in July while the lowest emission was detected in January.

Overall, the total BVOCs emission was peaked in summer (June, July, and August) and bottomed out in winter(December,January and February). In spring,the BVOCs emission showed an upward trend in which the isoprene was increased the most. Despite the emissions of both monoterpenes and other VOCs also increased,the values recorded were comparatively lower than the isoprene that could be due to the more significant influence of sunlight and temperature on isoprene. Under these conditions in spring, the vegetation was still in the growing stage where the leaves were immature with regular enzyme activity. Therefore, the emission of BVOCs in spring only accounted for 17.8% in the whole year.Comparatively, the average daily temperature and sunshine time were significantly increased in summer than that in spring. Under these conditions, the vegetation leaves reached the mature stage and grown with maximum leaf area; consequently, the highest enzyme activity was achieved. As a result, high emission of BVOCs were observed in summer with the emission of 43.8 Gg C·year-1that accounting for 61.9% of the total annual emissions. In autumn, the temperature difference was noticeable. The leaves were transformed from the mature leaves to the decaying leaves. Therefore, the emission of BVOCs showed a decreasing trending in September, followed by a sharp decrease in October, and finally constant in November. The total emission of BVOCs recorded in autumn at only 18.5%. Due to low temperatures and limited solar irradiation in winter, the emission of BVOCs reached the lowest value of the year at only 1.8%.

Spatial distribution

The leaf biomass data of each dominant forest species from 13 cities in the Jing-Jin-Ji area were obtained based on the forest resources data and field survey, while the annual meteorological data of each city were retrieved from the NASA website. Then, the standardized emission factors,leaf biomass data and meteorological factors were calculated according to Eqs. 6 and 7 to obtain the annual BVOCs emission data of each forest species in 13 cities of the Jing-Jin-Ji area in 2017. The BVOCs emission of dominant forest species in 13 cities (Beijing,Tianjin, Baoding, Cangzhou, Chengde, Handan, Hengshui, Langfang, Qinhuangdao, Shijiazhuang, Tangshan,Xingtai, Zhangjiakou) throughout 2017 were calculated to study the spatial distribution of BVOCs emission in the Jing-Jin-Ji area. As shown in Table A4 and Fig. 5, the BVOCs emission fluxes and the compositions in the Jing-Jin-Ji area demonstrate an apparent spatial distribution. Given that Chengde and Beijing have high coverage of vegetation, and the presence of dominant species(Betula platyphylla, Quercus variabilis, Populus tomentosa and Pinus tabulaeformis) in Wuling Mountain Reserve and Saihanba Forest Farm in Chengde showed higher BVOCs emission rates. Therefore, Chengde contributed the highest BVOCs emission of dominant forest species in the Jing-Jin-Ji area with 20.4 Gg C·year-1(28.8%), followed by Beijing with a total discharge of 17.6 Gg C·year-1(24.9%). The remaining proportion of BVOCs emission was detected from Baoding, Tangshan,Hengshui and Zhangjiakou. Furthermore, the detected BVOCs emission of Cangzhou (0.4 Gg C·year-1, 0.6%)and Langfang (0.5 Gg C·year-1, 0.8%) were less than the others due to their smaller city area and lower vegetation coverage. In terms of compositions (Fig. A2), Chengde presents the largest isoprene emission with 26.8% (7.7 Gg C·year-1) of the total isoprene emissions released by dominant forest species in the Jing-Jin-Ji area. This could probably due to the extensive vegetation coverage of deciduous trees with high isoprene emissions such as Quercus variabilis. Asides from Chengde, Beijing also shows a high emissions of monoterpenes (9.4 Gg C·year-1) that accounting for 36.9% of the total monoterpenes emissions. This could be attributed to the emission of α-pinene and β-myrcene from coniferous such as Platycladus orientalis and Pinus tabuliformis,which could also explain the highest emissions detected for Chengde in spring, summer, and autumn, but lower than Beijing particularly in winter. Chengde possesses the most substantial difference in emissions in winter and summer because the dominant forest species here are mostly deciduous trees. The low temperature and less sunshine in winter would have caused the leaf biomass of deciduous trees to emit lower BVOCs; thus the emission was reduced significantly (Fig. 6). In general,the distribution of BVOCs emission fluxes was highly consistent with the distribution of vegetation.

Overall, the BVOCs emission estimated in this study were much lower than those estimated by Klinger et al.(2002) (Table 3). It is due to the difference in the scope of the study objects where this study focused on the specific dominant forest species while Klinger et al. (2002)covered all species in this area, including grasslands,shrublands, forests,and peatlands.

Chemical activity evaluation

The contribution of BVOCs to atmospheric chemical reaction depends on the level of its emissions and closely related to their chemical activity. Since the olefins with double bonds are more active compounds in addition to the main BVOCs components released by dominant forest species in the Jing-Jin-Ji area are isoprene and monoterpenes, the chemical activity of isoprene and monoterpenes by the maximum incremental reactivity(MIR) method and the ·OH reaction rate (LOH) method were analyzed comprehensively. Figure 7 shows the OFP and LiOHvalues and activity contribution rates of every dominant forest species. Overall, the activity contribution rate of isoprene and monoterpenes of each species calculated using two methods are basically consistent.Among them, the OFP values of Robinia pseudoacacia,

Fig.5 The BVOCs emission(106 g C·year-1) (a)and proportion of each component (%)(b) in each city in the Jing-Jin-Ji area(BJ:Beijing;TJ:Tianjin;BD:Baoding;CZ:Cangzhou;CD:Chengde;HD:Handan;HS:Hengshui;LF: Langfang; QHD:Qinhuangdao; SJZ:Shijiazhuang;TS:Tangshan;XT:Xingtai;ZJK:Zhangjiakou)

Fig.6 The seasonal BVOCs emission distribution of dominant forest species in the Jing-Jin-Ji area(106 g C)

Populus tomentosa and Quercus variabilis obtained by the MIR method were higher, followed by Betula platyphylla and Pinus tabuliformis while the lowest OFP value was obtained for Ulmus pumila. The forest species with higher LiOHvalue obtained by the LOHmethod were Robinia pseudoacacia, Populus tomentosa, Platycladus orientalis, Pinus tabuliformis and Quercus variabilis while Larix gmelinii had the lowest LiOHvalue. The difference between the two methods was due to their different principles in which the LOHmethod reflects the reactivity by calculating the ability of BVOCs and OH radical to produce RO2without considering the influence of other subsequent reactions. Although the MIRmethod considers a series of responses of BVOCs, the lack of some MIR coefficients also affect the results.Based on the two approaches, Robinia pseudoacacia, Populus tomentosa,Quercus variabilis and Pinus tabulaeformis are the dominant forest species that contribute higher BVOCs reaction activity in the Jing-Jin-Ji area. Among them, Robinia pseudoacacia,

Table 3 Comparison of estimated BVOCs emission in different regions

Populus tomentosa and Quercus variabilis have higher reactivity due to their higher isoprene activity contribution rate.The OFP and LOHvalue of Robinia pseudoacacia, Populus tomentosa, Quercus variabilis were 156.76 μg·m-3and 138.36×105, 87.19 μg·m-3and 78.86×105, 87.19 μg·m-3and 78.86×105,respectively.In addition,Pinus tabulaeformis was detected with a higher activity contribution rate of monoterpenes, especially β-myrcene, limonene, and αpinene, hence the OFP and LOHwere detected at 43.63 μg·m-3and 54.67×105,respectively.

Uncertainty in BVOCs emission estimates

Since the emissions of BVOCs by vegetation are influenced by multiple factors, the uncertainty needs to be taken into consideration for the estimation. The primary sources of uncertainty in BVOCs emission estimation include emission factors, leaf biomass, vegetation distribution, model algorithm and meteorological parameters.Estimation of total BVOCs emission released by vegetation on the ground in the Jing-Jin-Ji area was studied using different algorithms and data sources. Different models, emission factors or vegetation type data would produce considerably different emissions (approximately two-fold) (Carlton and Baker 2011; Hogrefe et al. 2011;Chen et al. 2018). Given the above, the BVOCs emission estimated in this study are based on the field vegetation investigation and measured forest species emission rates as an effort to eliminate the uncertainties further. However, the emission rates obtained through experiments are mostly related to plant physiological conditions.Hence, the differences in the values could be due to uncertainty or uncontrollable stressors, environmental conditions, or genetically related metabolic processes(Loreto et al. 2009; Monson et al. 2013). For the sampling method, the dynamic headspace method used in this study has more air circulation than the static system, thus the real existing environment of plants can be maintained in good condition. However, the disturbance to plants during the process and the unrealistic high concentration of H2O caused by transpiration may interfere with the leaf stomata; consequently the emissions could be impacted (Ortega et al. 2008). Moreover, the rapid reactivity of BVOCs would also affect the deviation of flux and the analysis performed by the GC-MS instrument may also possess potential measurement errors. As for the emission factors in this study, the influence of solar radiation on monoterpene emissions is not considered in the calculation. The emission mechanism of different monoterpenes varies greatly in which some of the monoterpenes are reported to be positively correlated with temperature (Guenther et al. 2012). The harvesting of leaf biomass was calculated via biological parameters such as forest vegetation accumulation and trunk density. It should be noted that the manual measurement errors would also affect the accuracy of the calculation. Furthermore, the meteorological data used in this study were retrieved from the MODIS satellite published on the NASA website. The acquisition, processing, analysis and conversion of satellite data would further introduce various types and degrees of uncertainty. For example, the positioning error of the ground station, the positioning error of the satellite and the tropospheric time delay will cause the accuracy and stability of the satellite calibration and the delay error in time transmission. In addition, in the inversion process of meteorological satellite remote sensing, atmospheric state deviation, model error, model parameter error, etc.,will also affect the inversion accuracy. In detail, there will be a systematic deviation between the reference temperature profile and the actual atmospheric profile during the seasonal transition, and solar radiation will cause large diurnal changes in the temperature of the troposphere bottom, which will cause inversion errors,and the errors caused are more complicated. Therefore,the estimated BVOCs emission using the satellite detection data may be higher or lower than the actual observed BVOCs emission. By comparing multiple sets of monthly measured data and satellite data, the difference between them is less than 5%. Overall, it is suggested that future work should be emphasized on the reduction and prevention of those uncertainties to obtain more accurate data.

Fig.7 The reaction activity and activity contribution rate of BVOCs of dominant tree species in the Jing-Jin-Ji area calculated by MIR method(a)and LOH method(b)

Conclusions

The total annual emission of BVOCs released from dominant forest species in the Jing-Jin-Ji area was estimated to be 70.8 Gg C·year-1. Isoprene, monoterpenes,and other VOCs contributed 40.5%, 36.0%, and 23.4%,respectively. As for monthly and seasonal variations, the emission of BVOCs was peaked in summer (June, July,and August) and bottomed out in winter (December,January, and February). This explains the summer represents the season with the most severe ozone pollution.In terms of spatial distribution, high BVOCs emission was mainly distributed in the Taihang Mountains and the Yanshan Mountains. Chengde contributed the most followed by Beijing, inferring that these two cities should have more attention to mitigate BVOCs emission to atmospheric ozone and particulate pollution. In terms of forest species, Platycladus orientalis, Quercus variabilis,Populus tomentosa, Pinus tabulaeformis and Robinia pseudoacacia were the primary contributors of BVOCs emission and atmospheric reactivity. In conclusion, it was recommended that the use of these forest species for greening configuration should be avoided since they negatively impact the atmosphere due to high emission or reactivity of BVOCs.

Abbreviations

BVOCs: Biogenic volatile organic compounds;AVOCs: Anthropogenic volatile organic compounds; G93: Guenther model developed in 1993;NMVOCs: Non-methane volatile organic compounds; PANs: Peroxyacetyl nitrate; SOA: Secondary organic aerosols; PM: Particulate matter;OH: Hydroxyl radical; PFT: Plant functional types; PAR:Photosynthetically active radiation; MIR:Maximum incremental reactivity; LOH: OH reaction rate;OFP: Ozone formation potential; OVOCs: Other VOCs; TVOCs: Total VOCs;NIST: National Institute of Standards and Technology; PAMS: Photochemical Assessment Monitoring Stations

Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s40663-021-00322-y.

Additional file 1: Table A1.The parameters of dominant forest species in the Jing-Jin-Ji area.Table A2.The annual emissions of BVOCs from dominant forest species in the Jing-Jin-Ji area (106g C·year–1).Table A3.The monthly emissions of BVOCs from dominant forest species in the Jing-Jin-Ji area (106g C).Table A4.The BVOCs emission of 13 cities in the Jing-Jin-Ji area (106g C·year–1).Figure A1.Contributions of dominant forest species to the components of BVOCs in the Jing-Jin-Ji area.Figure A2.The contribution of BVOCs emitted from the dominant forest species in each city in the Jing-Jin-Ji area.

Acknowledgements

Not applicable.

Authors’contributions

Xiaoxiu Lun, Qiang Wang and Ying Lin conceived the study; Ying Lin participated in study design, field measurements,data processing and writing the manuscript;Xiaoxi Jing and Chong Fan took part in the field measurements; Wei Tang and Zhongzhi Zhang took part in the procedure of data processing. All authors contributed critically to successive drafts and gave final approval for publication.

Funding

This work was supported by the grants from National Natural Science Foundation of China (No.42077454), National Research Program for Key Issues in Air Pollution Control (DQGG202126),National Natural Science Foundation of China (No.41605077).

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable requests.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1College of Environmental Science and Engineering, Beijing Forestry University, 100083 Beijing, China.2Institute of Atmospheric Environment,Chinese Research Academy of Environmental Sciences, 100012 Beijing, China.

Received: 22 October 2020 Accepted: 3 June 2021