Impact of Land Cover and Leaf Area Index on BVOC Emissions over the Korean Peninsula

: Biogenic volatile organic compound (BVOCs) emissions are the largest VOC emission source globally, and are precursors to ozone and secondary organic aerosols, both of which are strong, short-lived climate pollutants. BVOC emissions are usually estimated using the Model of Emissions of Gases and Aerosols from Nature (MEGAN), which requires Plant Functional Types (PFTs) and Leaf Area Indexes (LAIs) as inputs. Herein, the e ﬀ ects of reﬁned input data on regional BVOC emission estimates are analyzed. For LAIs, lower resolution MODerate-resolution Imaging Spectroradiometer (MODIS), and higher spatio-temporal resolution Spatial and Temporal Adaptive Reﬂectance Fusion Model (STARFM) LAI were generated. For PFTs, local land cover maps were developed, in addition to MODIS PFT. In South Korea, annual emissions of isoprene and monoterpenes in 2015 were estimated as 384 and 160 Gg / year, respectively, using STARFM LAI and Local PFT (Case 4). For North Korea, 340 Gg / year isoprene and 72 Gg / year monoterpenes emissions were estimated using STARFM LAI and MODIS PFT. These estimates were 14–110% higher than when using MODIS LAI and MODIS PFT (Case 1). Inter-comparison with satellite-based inverse isoprene emission estimates from GlobEmission shows 32% (North Korea) to 34% (South Korea) overestimation in bottom-up data. Our new vegetation inputs improve MEGAN performance and resulting BVOC emission estimations. Performance of Weather Research and Forecasting (WRF) meteorological modeling requires improvement, especially for solar radiation, to avoid overestimation of isoprene emissions.


Introduction
Climate change theory suggests that ecosystems can become potentially large volatile organic compound (VOC) sources [1]. Biogenic VOCs (BVOCs) emitted from plants are involved in ozone generation through photochemical reactions, and are precursors for the formation of secondary organic aerosols [2][3][4]. VOC emissions from vegetation are comparable to those emitted from anthropogenic sources, with relative emission values varying across regions and seasons [5]. The Korean Peninsula is located in the far east of the Eurasian continent where many industrial countries are located, and anthropogenic emissions are regionally very strong. The combination of these anthropogenic VOCs and natural BVOCs acts as a precursor to ozone and secondary organic aerosol formation.
where F i is the emissions of chemical species i (µg/m 2 h), γ i is the emissions activity of chemical species i, ε i , j is the standard emission factor for chemical species i in vegetation type j, and χ j is the percentage of vegetation type j occupying a grid square. Several vegetation environmental conditions are considered through the genetic factor expression reaction (2), when estimating the emission response, including light reaction by electron transfer, temperature response by in vivo enzyme activity, in vivo metabolite change, enzyme activity, and CO 2 levels. The main consideration is response to light, γ P,i , which indicates the change of emission activity of chemical species i with light intensity, and is represented by the light-dependent fraction (LDP) and photosynthetic photon flux density (PPFD). Variable γ T,i indicates the change in the activity of chemical species i depending on temperature, expressed as its light-independent fraction (LIF); γ A,i is the discharge activity according to the age of the leaf-which is classified by MEGAN v2.1 into either New, Grow, Mature, or Old. The variable γ SM,i is the emission activity value according to the soil moisture content of the vegetation area, while γ C,i is the vegetation emission activity value according to atmospheric CO 2 concentration; γ SM,i and γ C,i are 1 for all chemical species apart from isoprene. The canopy environment coefficient (C CE ) and LAI are also used for calculating γ values: γ i = C CE LAIγ P,i γ T,i γ A,i γ SM,i γ C,i .
In addition, the PFT fraction and EF are used to calculate standard compound emission factors, expressed as grid emission factors, under standard conditions (temperature: 303 K; photosynthetic Atmosphere 2020, 11, 806 3 of 19 photon density: 1000 µmol/s) [8]. Emission factors provided by MEGAN v2.1 [8] using various observations and experiments have been reported previously.

Plant Functional Type
PFT is a scheme used to classify vegetation and enables emissions estimation from vegetation species on a global basis. It is used in research land use and vegetation responses to climate change, and to classify vegetation according to climatological, phylogenic, and biological characteristics. In general, PFT can be obtained from satellite data, vegetation lists, land cover maps, and ecosystem maps, but is mainly composed of satellite data [13]. In MEGAN v2.1, the PFT scheme involves 16 vegetation types, with each type designated using information on vegetation type and climate [8].

Leaf Area Index
Remote sensing, which has been widely used in recent vegetation and forest observational studies, provides information on a periodic basis, and is suitable for studying large areas at low cost [14]. Here, satellite-based remote sensing data are useful as the scope includes the entire Korean Peninsula-including North Korea, which cannot be easily accessed on the ground. Vegetation Index is a radiant value that represents the relative distribution and activity of green plants on the surface. LAI is a parameter related to plant physiology activity, and can be calculated using the normal distribution vegetation index (NDVI)-which is itself calculated using reflectance values from the red visible (Red) and near infrared (NIR) regions, for green plants [15,16]: NDVI is the most commonly used vegetation index; values range from −1 to 1, with those between 0.1 and 0.6 being most common. The equation for calculating LAI using NDVI is shown as Equation (4) [17]: LAI = 6.7537 × NDVI + 0.8384 (4)

BVOC Emissions Modeling Framework and Meteorological Input
The MEGAN modeling domain for BVOC emission estimation covered the Korean Peninsula ( Figure 1). The modeling year was 2015, and meteorological data were produced by the Weather Research and Forecasting (WRF) model [18]. MEGAN v2.1 Emission Factor (EF), provided by the MEGAN development group, was used as the emission factor dataset, and detailed modeling framework information is listed in Tables 1 and 2. In this study, four emission estimation cases were used to examine the effects of PFT and LAI on BVOCs estimation (

BVOC Emissions Modeling Framework and Meteorological Input
The MEGAN modeling domain for BVOC emission estimation covered the Korean Peninsula ( Figure 1). The modeling year was 2015, and meteorological data were produced by the Weather Research and Forecasting (WRF) model [18]. MEGAN v2.1 Emission Factor (EF), provided by the MEGAN development group, was used as the emission factor dataset, and detailed modeling framework information is listed in Tables 1 and 2.

Plant Functional Type
PFT data used in this study were MODIS MCD12Q1 products [8]. MCD12Q1 has a resolution of 500 m and provides 5 types of land cover, with Type 5, which provides 14 PFT types, being used in this study. Mapping of data with different classification schemes is required to apply the type to MEGAN v2.1, and so applying the MODIS Land Cover information into the MEGAN PFT format enabled mapping of 7 of the 16 vegetation categories in MEGAN [8] Figure 2. Natural grasslands GS_C3_WARM 14 420 Artificial grasslands  The most abundant vegetation type in South Korea was identified as Broadleaf forests using both MODIS PFT and local PFT data, accounting for 53%, and 41%, respectively. In MODIS PFT, other major categories include Farmland (25%) and Needleleaf forest (6%); with local PFT, other major categories were Needleleaf forest (35%) and Farmland (20%). Compared with MODIS PFT, the largest change in local PFT data was an increase of 30% in Needleleaf forest, followed by a 13% difference for Broadleaf forest, 9% for non-vegetation, and there 6% less was categorized as Farmland ( Figure  3). The most abundant vegetation type in South Korea was identified as Broadleaf forests using both MODIS PFT and local PFT data, accounting for 53%, and 41%, respectively. In MODIS PFT, other major categories include Farmland (25%) and Needleleaf forest (6%); with local PFT, other major categories were Needleleaf forest (35%) and Farmland (20%). Compared with MODIS PFT, the largest change in local PFT data was an increase of 30% in Needleleaf forest, followed by a 13% difference for Broadleaf forest, 9% for non-vegetation, and there 6% less was categorized as Farmland ( Figure 3).  In this study, we did not produce a local PFT for North Korea, as closed access makes it very difficult to collect ground-survey data for land cover, terrain, weather, etc. However, the MODIS PFT data itself is supportive in understanding North Korean biogenic emissions. In this study, we did not produce a local PFT for North Korea, as closed access makes it very difficult to collect ground-survey data for land cover, terrain, weather, etc. However, the MODIS PFT data itself is supportive in understanding North Korean biogenic emissions.

Leaf Area Index
Most BVOC studies use LAI created using the MODIS MCD15A2H product [15]. The MODIS LAI is composed of 46 data sets, created at 8-day intervals, each with a spatial resolution of 500 m. However, in order to produce LAI with higher spatial resolution than MODIS LAI, Landsat satellite data, with a spatial resolution of 30 m, were used. A total of 23 images were needed to cover the entire Korean Peninsula area, meaning that 276 Landsat satellite images were needed annually, to achieve monthly peninsula coverage. Weather interference meant that not all 276 images were usable, and so a simulated satellite data reconstruction method, known as STARFM (Spatial and Temporal Adaptive Reflectance Fusion Model), was used.
The STARFM method is a satellite data fusion technique, developed by Gao [22] to predict Landsat satellite data for a desired day, by fusing high spatial resolution Landsat satellite data with high temporal resolution MODIS sensor data. From the Landsat satellite, we used 11 band images with a resolution of 30 m for the whole Earth, while the MODIS product used in STARFM was MCD43A4, which provides 7 band images with 50 m resolution. Among these, NDVI was created using the relationship between the Red and NIR bands, and LAI was obtained through the relationship between NDVI and LAI [17].
MEGAN v2.1 LAI v was calculated by dividing LAI (the leaf area index by grid using satellite data) by the vegetation share, F v (5). LAI v should have a value between 0 and 6 [6]: Our MODIS LAI used the MODIS MCD15A2H product [15], and our STARFM LAI used the LAI created using the STARFM method. As discussed, LAI is based on the reflectance of vegetation observed from the satellite, and LAIs for August and December in 2015 are shown in Figure 4. It can be seen that, overall, STARFM LAIs were higher than those from MODIS, which show fewer dark shaded areas (i.e., areas near zero).
For STARFM LAI, distinct edge effects were detected after creating a peninsular-wide mosaic. To produce the STARFM LAI for the Korean Peninsula, 23 pieces of LANDSAT data were required-although simultaneous imagery could not always be used owing to weather conditions (e.g., cloud interference). There were differences between some LANDSAT images owing to image timings, and these showed as a distinct edge. However, these differences should disadvantage emissions estimation much less than the benefit gained by using this more refined, higher resolution data.
Total LAIs, divided into South Korea and North Korea, are summarized in Table 5. Since LAI represents leaf activity, it shows seasonal change, with lower values in winter.
Atmosphere 2020, 11, 806 8 of 19 Atmosphere 2020, 11, x FOR PEER REVIEW 8 of 20 these showed as a distinct edge. However, these differences should disadvantage emissions estimation much less than the benefit gained by using this more refined, higher resolution data.

Meteorological Input Data
Meteorological data were prepared by matching spatial ranges with the biogenic emissions model. The meteorological data used in this study were written every hour in 2015, the target period, using the WRF v3.7.1 model. The Meteorology-Chemistry Interface Processor (MCIP) was then used to convert this meteorological model output into a MEGAN-ready format. (Meteorology model configuration: Table 6). We conducted verification by estimating the correlation between WRF meteorological model data used in MEGAN modeling and daily measurements (averaged from hourly measurements) from ground meteorological sites. Since there are no ground measurements available for North Korea, the selection of meteorological measurement site was limited to South Korea. There are 590 ground measurement sites (45 of which are manned) in South Korea. Seoul is one of the largest megacities in the world and the capital city of South Korea; it represents developed land in Korea. Chuncheon and Wonju are located in the Gangwon province, the province closest to North Korea, and represent a forest rich area. Chupungnyeong is located in the central hill of South Korea, and represents a forest rich area of the southern peninsula.
Correlations for temperature and radiation between model and weather station measurements are shown in Figures 5 and 6. Data distributions were calculated from daily averages during 2015.
the world and the capital city of South Korea; it represents developed land in Korea. Chuncheon and Wonju are located in the Gangwon province, the province closest to North Korea, and represent a forest rich area. Chupungnyeong is located in the central hill of South Korea, and represents a forest rich area of the southern peninsula.
Correlations for temperature and radiation between model and weather station measurements are shown in Figures 5 and 6. Data distributions were calculated from daily averages during 2015.  For temperature ( Figure 5), slopes were mostly >0.9; the WRF model results showed a tendency for slight underestimation, resulting in a negative bias. Correlation coefficient R showed high correlation (≥0.9). For solar radiation (Figure 6), slopes were 0.9~1.1; the WRF model had a tendency to overestimate, resulting in a positive bias. Correlation coefficient R showed high correlation (0.9). These results confirm that the WRF model produces reliable meteorological information for MEGAN modeling.  Table 7, the R of temperature data was very high (0.98-0.99) but the model tended to underestimate temperature slightly (MBs of −3.47 K~−2.32 K, except for Chupungnyeong at 0.59 K). In contrast, Solar radiation (  For temperature ( Figure 5), slopes were mostly >0.9; the WRF model results showed a tendency for slight underestimation, resulting in a negative bias. Correlation coefficient R showed high correlation (≥0.9). For solar radiation (Figure 6), slopes were 0.9~1.1; the WRF model had a tendency to overestimate, resulting in a positive bias. Correlation coefficient R showed high correlation (0.9). These results confirm that the WRF model produces reliable meteorological information for MEGAN modeling. Tables 7 and 8 Table 7, the R of temperature data was very high (0.98-0.99) but the model tended to underestimate temperature slightly (MBs of −3.47 K~−2.32 K, except for Chupungnyeong at 0.59 K). In contrast, Solar radiation (Table 8)     For BVOC emissions estimation, large seasonal variations would be expected in each scenario. To investigate the causes of seasonal differences, we analyzed LAI input data by comparing monthly average MODIS LAI and STARFM LAI data ( Table 5) with 2-m-above-ground temperature (TEMP2) data (Figure 7). Both MODIS LAI and STARFM LAI showed low values in winter and high values in summer, with STARFM LAI values higher than those produced by MODIS LAI, particularly in winter, where values were 2~5 times higher (  For BVOC emissions estimation, large seasonal variations would be expected in each scenario. To investigate the causes of seasonal differences, we analyzed LAI input data by comparing monthly average MODIS LAI and STARFM LAI data ( Table 5) with 2-m-above-ground temperature (TEMP2) data ( Figure 7). Both MODIS LAI and STARFM LAI showed low values in winter and high values in summer, with STARFM LAI values higher than those produced by MODIS LAI, particularly in winter, where values were 2~5 times higher (Table 5). For STARFM LAI, unlike MODIS LAI, LAI values were > 1 when TEMP2 was below freezing (273 K). According to [23], the NDVI can be positive, even below 273 K (0 °C). Furthermore, at 293 K (20 °C), the NDVI value can be as low as 0.4. Even at the same temperature, the NDVI value should be supported by precipitation data. The fact that NDVI is positive means that the LAI will also be According to [23], the NDVI can be positive, even below 273 K (0 • C). Furthermore, at 293 K (20 • C), the NDVI value can be as low as 0.4. Even at the same temperature, the NDVI value should be supported by precipitation data. The fact that NDVI is positive means that the LAI will also be positive, according to Equation (4). Since the average temperature in the Korean Peninsula is over 20 • C in summer, and rainfall is frequent (~83 days per year in 2015), the relationship between temperature, precipitation, and NDVI shows that STARFM LAI is a properly estimated input. In Siberia, vegetation is completely dominated by evergreen needleleaf spruce, fir, and pine; in below-freezing conditions, needleleaf forest that can tolerate cold, and evergreen forest that photosynthesizes regardless of season are required for NDVI and LAI to exhibit positive values [23].
From Figure 3, the composition of evergreen needleleaf forest increased from 5.9% in MODIS PFT to 32.7% in local PFT. With the local PFT, the ratio of evergreen forest to deciduous forest was three-fourth, which is more than four times that of MODIS PFT. From Figure 7, STARFM shows higher LAI values than MODIS during winter months, which also supports the high composition of evergreen forest in South Korea

BVOC Emissions Estimation
Using the four modelling approaches, monthly vegetation emission estimates were calculated for isoprene and monoterpenes emissions (Tables 9 and 10). Overall, monthly emissions were low in winter and high in summer, as expected, and there was also an overall difference in emission amounts. South Korea exhibited its highest emission levels in August, and North Korea showed its highest emissions in July; mean temperatures and the monsoon are potential causes of these differences. In South Korea, four cases of model emissions (Case 1, Case 2, Case 3, and Case 4) were estimated, based on a combination of input data (Tables 1 and 2). Isoprene emissions were lowest in January (0.36, 0.30, 0.56, and 0.62 Gg/month for the respective cases). In August, the highest emission month, values were 118.83, 95.98, 123.13, and 109.56 Gg/month for the respective cases. For monoterpenes, the month with the lowest emissions was classified according to the type of LAI. When using MODIS LAI (Case 1, Case 2), January was the lowest month (0.42 and 0.66 Gg/month, respectively). When STARFM LAI (Case 3, Case 4) was used, emissions were 0.85, and 1.64 Gg/month, respectively. in February.
Monoterpene emissions from PFT showed the opposite trend to that of isoprene emissions. The lowest isoprene emission figure was 328.49 Gg/year (Case 2), and the highest was 435.48 Gg/year (Case 3). The lowest monoterpene emission reading was 79.71 Gg/year (Case 1) and the highest was 160.33 Gg/year (Case 4).
The emission estimate amounts differed between cases and species. When the LAI changed from MODIS LAI to STARFM LAI, emissions of both isoprene and monoterpenes increased (by 50% for isoprene and 70% for monoterpenes). When the PFT changed from MODIS PFT to Local PFT (South Korea), isoprene emissions decreased by 20%, while monoterpenes emissions increased by 60%.
When using the same LAI, the use of Local PFT in cases 2 and 4 reduced isoprene emissions and increased monoterpenes emissions because the composition of the PFT had changed (Figure 3). The portion of broadleaf trees decreased by 10% in the local PFT compared with the MODIS PFT, and the needleleaf trees increased by~6 times. When we examined the Emission Factor (EF) for each PFT used in MEGANv2.1 [8], isoprene EF decreases by 95% when changing from broadleaf to needleleaf trees; in monoterpenes, it increases by 50% for the same change. For this reason, when the Local PFT was used, isoprene emissions were reduced, while monoterpene emissions increased.
North Korea produced two emission cases (Case 1, Case 4) to changes in LAI (Table 10). The lowest monthly isoprene emission estimates were for January (0.08 Gg/month in Case 1, and 0.18 Gg/month in Case 4), and the highest were in July (84.49 and 92.40 Gg/month, for cases 1 and 4, respectively). Total annual emission estimates for North Korea (Case 1 and Case 3) were 297.63 and 340. 16 Gg/year, respectively, with Case 4 estimates 14% higher than those of Case 1 values.
For monoterpenes, January produced the smallest emission estimate in North Korea, with Case 1 yielding an estimate of 0.17 Gg/month, and Case 4 yielding 0.46 Gg/month. The largest monthly monoterpene emission estimates were in July, at 13.97 Gg/month and 16.26 Gg/month. Comparing Case 1 (59.43 Gg/year) to Case 4 (71.98 Gg/year). Case 4 values were 21% larger than Case 1 monoterpene emissions in North Korea. In most cases, North Korea showed the highest emissions in July (Table 10). Unlike the rest of the Korean Peninsula, temperatures in North Korea were higher in July than in August, causing North Korea's BVOC emissions to peak then, rather than in August.
Case 4 BVOC emission estimates were −8-101% higher than Case 1 estimates, depending on whether isoprene or monoterpene data were reviewed. However, when compared on a monthly basis, the differences were larger, with cold season estimates (January-March and November-December) more than twice as high; between May and August, differences were approximately 10%. In the cold season, local PFT and STARFM LAI estimates were relatively high compared with those of MODIS PFT and LAI, as described in Section 3.1, reflecting a relative increase in Evergreen needle leaf forests in the local PFT.

Comparison with Inverse Estimates of BVOC Emissions
We used inverse emissions estimates from GlobEmission [24] and observational data for 2015 to compare vegetation emissions from the Korean Peninsula. The isoprene emissions, provided by GlobEmission, were estimated using HCHO (formaldehyde) column concentrations, measured with an OMI (Ozone Monitoring Instrument) sensor. As of June 2020, GlobEmission provides monthly Global Biogenic Isoprene emissions for 2005-2014 (but not 2015). To compare isoprene emissions over the Korean Peninsula, the results of this study were compared with GlobEmission estimations for 2014. To make this comparison, we performed WRF modeling for 2014.
South Korea's monthly isoprene emissions in 2014 and 2015 were compared using Case 4 ( Figure 8). Emissions in 2014 were 349.88 Gg/yr, 10% less than in 2015. Monthly variations showed greater variation than annual estimates. Usually, August has the highest temperature and humidity, which induces high vegetation emissions. In August 2014, isoprene emissions were 37% less than those of 2015 ( Figure 8). In 2014, two large typhoons affected South Korea, causing heavy rain in most parts of the country. Therefore, the average temperature was lower and the average precipitation was higher than the 30year average (year 1981~2010). Furthermore, more days of precipitation were recorded [25].
Low daylight hours (74% of a normal year), were reflected in the lower isoprene emissions in 2014. As such, emissions in August 2014 and 2015 showed a marked difference; however, Figure 9 shows that monthly emissions in 2014 were well estimated by the bottom-up modeling, since emissions showed a similar trend to GlobEmission. 8). Emissions in 2014 were 349.88 Gg/yr, 10% less than in 2015. Monthly variations showed greater variation than annual estimates. Usually, August has the highest temperature and humidity, which induces high vegetation emissions. In August 2014, isoprene emissions were 37% less than those of 2015 ( Figure 8). In 2014, two large typhoons affected South Korea, causing heavy rain in most parts of the country. Therefore, the average temperature was lower and the average precipitation was higher than the 30year average (year 1981~2010). Furthermore, more days of precipitation were recorded [25]. Low daylight hours (74% of a normal year), were reflected in the lower isoprene emissions in 2014. As such, emissions in August 2014 and 2015 showed a marked difference; however, Figure 9 shows that monthly emissions in 2014 were well estimated by the bottom-up modeling, since emissions showed a similar trend to GlobEmission.   Monthly isoprene emissions in South Korea (a) and North Korea (b) are inter-compared with GlobEmissions in Figure 10. In general, bottom-up emissions showed similar monthly variations to GlobEmission emissions in 2014, with the highest emissions in July. In South Korea, the LAI and PFT for Case 1 and Case 4 were different, causing a difference in monthly emissions between these cases. Monthly emissions from May to July, when vegetation is active, showed a noticeable difference between cases, with estimations close to those of GlobEmissions in Case 4. Although emission amounts for Case 4 were larger than those of GlobEmission values (134%), the difference was smaller, considering that solar radiation shows a positive bias of ~15-20 (Table 8). However, for North Korea, the difference between Case 1 and Case 4 was only observed for LAI, with no significant difference in emissions (7%) or patterns.  Monthly isoprene emissions in South Korea (a) and North Korea (b) are inter-compared with GlobEmissions in Figure 10. In general, bottom-up emissions showed similar monthly variations to GlobEmission emissions in 2014, with the highest emissions in July. In South Korea, the LAI and PFT for Case 1 and Case 4 were different, causing a difference in monthly emissions between these cases. Monthly emissions from May to July, when vegetation is active, showed a noticeable difference between cases, with estimations close to those of GlobEmissions in Case 4. Although emission amounts for Case 4 were larger than those of GlobEmission values (134%), the difference was smaller, considering that solar radiation shows a positive bias of~15-20 (Table 8). However, for North Korea, the difference between Case 1 and Case 4 was only observed for LAI, with no significant difference in emissions (7%) or patterns. To compare emission estimation performances, bottom-up emissions for Case 1 and Case 4 were plotted with those of GlobEmission ( Figure 11). For Case 1, the slope was 1.59, and for Case 4, the slope was 1.09, with bottom-up emission results showing a tendency towards slight overestimation (i.e., positive bias). The correlation coefficient R showed high correlation (0.90). Even though both bottom-up emission cases overestimated isoprene emissions in June and July of 2014 for South Korea (Table 11), the degree of agreement was much better for Case 4, which is consistent with the previous analysis. To compare emission estimation performances, bottom-up emissions for Case 1 and Case 4 were plotted with those of GlobEmission ( Figure 11). For Case 1, the slope was 1.59, and for Case 4, the slope was 1.09, with bottom-up emission results showing a tendency towards slight overestimation (i.e., positive bias). The correlation coefficient R showed high correlation (0.90). Even though both bottom-up emission cases overestimated isoprene emissions in June and July of 2014 for South Korea (Table 11), the degree of agreement was much better for Case 4, which is consistent with the previous analysis.

Conclusions
The purpose of this study was to understand changes in BVOC emissions estimates using different input data. The biogenic emissions model used was MEGAN v2.1, and the target VOC species were isoprene and monoterpenes. Data related to vegetation were: 1) MODIS PFT and LAI; 2) local PFT, which is a fusion of Korean land cover and forest mapping, and 3) LAI using STARFM, a satellite data convergence algorithm. An important input dataset, meteorology, was obtained by taking a full year of WRF modeling and converting it using the MCIP. BVOC emission input data series were prepared by converting each of the above datasets into a form suitable for MEGAN v2. We attempted a comparison with biogenic isoprene emission data provided by GlobEmission, but could not make a direct comparison, as GlobEmission only provides emissions information up to 2014. For direct inter-comparison, we ran the WRF model for 2014 and performed emissions evaluation using inverse emissions estimation from GlobEmissions. For Case 1, emission estimates

Conclusions
The purpose of this study was to understand changes in BVOC emissions estimates using different input data. The biogenic emissions model used was MEGAN v2.1, and the target VOC species were isoprene and monoterpenes. Data related to vegetation were: (1) MODIS PFT and LAI; (2) local PFT, which is a fusion of Korean land cover and forest mapping, and (3) LAI using STARFM, a satellite data convergence algorithm. An important input dataset, meteorology, was obtained by taking a full year of WRF modeling and converting it using the MCIP. BVOC emission input data series were prepared by converting each of the above datasets into a form suitable for MEGAN v2.  16 Gg/year and 71.98 Gg/year, respectively, were estimated for North Korea. Emissions estimates between Case 1 and the other cases varied by −21-101% for isoprene and monoterpenes. Changes in LAI increased emission estimates for both isoprene and monoterpenes, while changes in PFT led decreased isoprene emissions and increased monoterpenes emissions estimates.
We attempted a comparison with biogenic isoprene emission data provided by GlobEmission, but could not make a direct comparison, as GlobEmission only provides emissions information up to 2014. For direct inter-comparison, we ran the WRF model for 2014 and performed emissions evaluation using inverse emissions estimation from GlobEmissions. For Case 1, emission estimates for 2014 were 373.85 Gg/year for South Korea and 332.82 Gg/year for North Korea, which is 143% and 125% of the GlobEmission isoprene emissions, respectively. For Case 4, bottom-up emissions were 349.88 Gg/year for South Korea and 351.83 Gg/year for North Korea, which is 134% of GlobEmissions for South Korea and 132% for North Korea. Correlation analysis and monthly inter-comparison indicate that our bottom-up emissions represent reasonable spatial and temporal variations using Case 4, but generally overestimated the amount of BVOCs.
We have tested BVOCs emissions estimation over the Korean Peninsula using multiple PFTs and LAIs. The results indicate that the enhancement of BVOC model input data using local PFT and refined LAI can effectively improve BVOC emissions estimation accuracy. However, the accuracy of the meteorological model requires improvement; currently, it is a major reason for uncertainties owing to overestimation of solar radiation.
In future study, the following improvements will be made. First, BVOCs emission factors require improvement by using more locally developed ones. Second, LAI estimation formulae need to be improved to be more specific to each PFTs. Third, soil moisture needs to be improved, especially for monoterpenes, since our present estimation does not include the effect of drought on terpenoid emissions. Fourth, understanding the performance of monoterpenes emissions estimation by direct evaluation using measurement data.