Impact of Climate-Driven Land-Use Change on O 3 and PM Pollution by Driving BVOC Emissions in China in 2050

: This study predicted three future land-use type scenarios in 2050 (including the Shared Socioeconomic Pathway SSP126, SSP585, and carbon scenario) based on the Land-Use Harmonization (LUH2) project and the future evolution of land-use types considering China’s carbon neutrality background. The contribution of land-use changes to terrestrial natural source biogenic volatile organic compounds (BVOCs), as well as O 3 and PM concentrations, were determined. Under the SSP126 pathway, meteorological changes would increase BVOC emissions in China by 1.0 TgC in 2050, compared with 2015, while land-use changes would increase them by 1.5–7.1 TgC. The impact of land-use changes on O 3 and PM concentrations would be less than 3.6% in 2050 and greater in summer. Regional differences must be considered when calculating future environmental background concentrations of pollutants. Due to more afforestation measures under the SSP126 scenario, the impact of land-use change on pollutants was more obvious under the SSP126 pathway than under the SSP585 pathway. Under the carbon scenario, the increase in PM concentration caused by land-use changes would pose a risk to air quality compliance; thus, it is necessary to consider reducing or offsetting this potential risk through anthropogenic emission control measures.


Introduction
Air pollution has received widespread attention in recent years as a result of its significant impact on the ecological environment and human health [1,2]. Although air quality in China has improved significantly [3], PM 2.5 concentrations in most cities still exceed the World Health Organization (WHO) standards (an annual average value of less than 5 µg m −3 ). Furthermore, the concentration of O 3 has been increasing in recent years. As a result, China is still confronted with the challenge of coordinating O 3 and PM 2.5 control [4].
A significant contribution to the generation mechanism of O 3 and PM 2.5 results from the transformation of pollutants discharged from natural sources [5][6][7][8]. Since emissions from natural sources are heavily influenced by meteorological conditions (such as temperature, wind speed, etc.) [9], and emissions from terrestrial natural sources are also influenced by land-use types [10][11][12][13], future climate and land-use types are important factors in quantifying changes in emissions from natural sources. The greening trend in China is causing an increase in BVOC and soil-source NO x emissions [11]. Jiang et al. [14] found that a 62% increase in urban land-use area from the 2000s to 2050s exerted more influence on O 3 than climate change. Feng and Li [15] and Liu et al. [16] also proved that concentrations of O 3 and PM 2.5 are related to their urban or suburban location based on the land-use regression model. However, land-use types will change as a result of future climate change, affecting the emission of terrestrial natural source emissions [10][11][12]; however, research on this topic is relatively limited. The study of future land-use changes is critical for more accurate quantification of future terrestrial natural source emissions and their pollutant contributions. Future changes in climate and land-use types will drive emissions from natural sources and further affect the background concentration of the atmospheric environment, which requires systematic quantification.
The method for calculating land-use type evolution in future years is based on base annual characteristics multiplied by the variation coefficient [17]. Predictions of future land-use scenarios vary greatly between models [18]. Many studies have been conducted in an attempt to forecast future land-use types [17,19,20]. Verburg et al. [19] developed two scenarios for assessing the uncertainty in agricultural land allocation in Europe between 2000 and 2030. The Land-Use Harmonization (LUH2) project [17] published the evolution of land-use types from the Shared Socioeconomic Pathway (SSP) to 2100, with a global resolution of 0.25 • × 0.25 • . The SSP126 scenario was created using the IMAGE 3.0 integrated assessment model. The SSP585 scenario's comprehensive assessment modeling framework consists of a regionalized investment and development model (REMIND) and a model of agricultural production and its environmental impact (MAgPIE) [17].
The Intergovernmental Panel on Climate Change released a report in 2018 calling on countries to take action to limit global warming to 1.5 K. At the 75th United Nations General Assembly in 2020 [21], the Chinese government proposed, "China will increase its nationally determined contribution, adopt more powerful policies and measures, strive to peak carbon dioxide emissions by 2030, and strive to achieve carbon dioxide emissions neutralization by 2060". However, research on the impact of land-use changes on atmospheric pollution at present is insufficient, particularly due to a lack of scenarios for carbon-neutral policies. To compensate for the shortcomings, the goal of this study is to forecast the evolution trend of land-use types in China in 2050, as well as quantify the impact of land-use changes on terrestrial natural source emissions as well as on pollution contributions.

Model Application
We used the Weather Research and Forecasting Model (WRF, version 4.2, http://www. mmm.ucar.edu/wrf/users/document.html, accessed on 11 June 2022) to simulate meteorological fields and the Meteorology-Chemistry Interface Processor (MCIP, version 5.1) [22] to process meteorological data in the format required by the Model of Emissions of Gases and Aerosols from Nature (MEGAN) and the Community Multiscale Air Quality Model (CMAQ). MEGAN (version 2.10) [23] was used to calculate biogenic volatile organic compound (BVOC) emissions, and the CMAQ model (version 5.3.2) [24] was configured to simulate O 3 and PM concentrations for various scenarios using the AERO6 aerosol module and the CB6 gas-phase chemical mechanism. Compared with the newer version of MEGAN3+, MEGAN2.10 does not take into account the effects of drought and non-drought period [25]. Based on limited observations from a study [26], results based on v2 and v3 were closer to observations with v2 overestimating isoprene concentration in almost all sites and v3 overestimating in urban areas. Further, when comparing total column tropospheric formaldehyde concentration, both v2 and v3 showed similar performance with respect to satellite retrievals, specifically in April [26].
Liu et al. [27] provided the input meteorological field data for the two SSP scenarios (SSP126 and SSP585), which are ensembles of five climate models: the Beijing Climate Center Climate System Model version 2 (BCC-CSM2-MR), a Climate Model developed at the  (CNRM-CM6-1), an Earth System Model  with Climate feedbacks and processes (EC-Earth3), a global Climate Model developed at  Institut Pierre-Simon Laplace (IPSL-CM6A-LR, LR stands for low resolution), and the Meteorological Research Institute Earth System Model version 2 (MRI-ESM2-0). Our 2015 upper air and surface meteorological and observation data analysis was based on datasets from the National Centers for Environmental Prediction (NCEP, datasets ds083.2, ds351.0, and ds461.0). To approximate the reanalysis data, the grid-nudging four-dimensional data assimilation method was used. The year-specific leaf area index (LAI) data were derived from a moderate-resolution imaging spectroradiometer (MODIS) LAI product (MOD15A2) [28]. The anthropogenic emissions inventory for 2015 in use for China is the Air Benefit and Cost and Attainment Assessment System-Emission Inventory (ABaCAS-EI) [29], which includes pollutants such as SO 2 , NO x , PM 10 , PM 2.5 , black carbon (BC), organic carbon (OC), VOCs, and NH 3 , as detailed by Zhao et al. [30] and Ding et al. [3], as well as CO, downscaled based on the ECLIPSE v5a [31]. In our simulations for 2050, we also used the anthropogenic emissions for 2015, which may be overestimated. In fact, the absolute contribution of natural source emissions to pollutant concentrations may decrease in 2050, as anthropogenic emissions are reduced, but the proportion will increase, compared with the results of this study. The boundary conditions for the WRF simulation were derived from the downscaling of the global model, and those for the CMAQ simulation were generated using the default configurations of CMAQ (clean atmosphere).

Preprocessing of LUH2 Data
In this study, SSP126 and SSP585 scenarios were chosen as climate change pathways to study the effect of future land-use change on atmospheric pollution. Land-use data based on the global LUH2 project [17] at 0.25 • × 0.25 • resolution were interpolated to a 27 km × 27 km resolution grid in order to perform regional simulations in China. The time scale, both in years, was consistent with the global data.
In this study, 12 different land-use types were considered-namely, primary forest, primary non-forest land, secondary forest land, secondary non-forest vegetation land, urban, C3 annual crops, C4 annual crops, C3 perennial crops, C4 perennial crops, C3 nitrogen-fixing crops, managed pastures, and other pastures-denoted as primf, primn, secdf, secdn, urban, c3ann, c4ann, c3per, c4per, c3nfx, pastr, and range, respectively. The Chinese government released a plan in 2016 to increase forest areas to 1.2 times those of the 2015 level by 2050. The proportion of forests in China may increase further in the future in order to achieve "carbon neutrality", which is an important way to increase carbon sinks. Even the SSP126 scenario with the lowest CO 2 emissions fails to meet the "carbon neutrality" goal and does not meet the "1.8 billion mu red line" for arable land. As a result, in order to achieve the "dual-carbon goal", this study increased the forest area in 2050 to 1.5 times that of 2015, controlled the decreasing trend of cultivated land to ensure the "red line of 1.8 billion mu" and developed a land-use "carbon-neutral" scenario, denoted as "carbon", to study the potential contribution to pollutants when vegetated land is increased as much as possible.
LAI is linearly related to vegetation cover, according to satellite inversion data. As a result, in this study, the LAI was calculated in 2050 using vegetation cover data from 2015 and 2050, as well as the LAI obtained through satellite inversion in 2015. The emission factors (EF) in 2050 were estimated based on the vegetation cover information in 2050 and the global EF in Guenther et al. [23]. The LAI and fraction of photosynthetically active radiation (FPAR) of vegetation in 2015 were fitted using the Beer-Lambert law; thus, the FPAR in 2050 could be calculated [32].
Since LUH2 and MEGAN use different vegetation classification methods, we referred to the method of Liao et al. [20] and summarized forest land and pasture into forest and grassland, respectively. LAI was used to fit 7 vegetation types, including 5 crops, and the correlation coefficient R 2 was greater than 0.7. This research employed two processing methods: ratio and grid. The former is based on the full grid data's average fitting result as the national average level, and the calculation for 2050 is in the form of the full grid multiplied by the same proportion, while the latter is based on the full grid data. The fitting coefficients were calculated one by one for each grid. The former method (ratio) is useful for capturing feature regions or feature grids, while the latter method (grid) is useful for reflecting the overall change trend between 2015 and 2050. In this study, we chose the grid processing method for the SSP scenario and the ratio processing method for the carbon scenario based on the principle of maximum data utilization and more reliable results (SSP grid). Table 1 shows the fitting coefficients for each variable in the carbon scenario using the ratio method (carbon ratio). The following equation can represent the relationship between FPAR and LAI: Since the coefficient k is related to the solar altitude, in this study, the simulated area was divided into four parts for fitting, from low to high latitude, and the overall correlation coefficient R 2 is around 0.95. Table S1 displays the fitting coefficient k and correlation coefficient R 2 for each month in the Supplementary Material. The fitting effect is strongest in low latitude regions and during the summer, which is related to the degree of vegetation coverage in the time and space where the input data is located. Typically, at low latitudes, or during the summer, the values of LAI and FPAR are higher, and data are more abundant. The LAI, EF, and FPAR calculated in 2050 based on the land-use type and vegetation cover information will be used as input data for MEGAN and CMAQ simulations to further quantify their impacts on terrestrial natural source emissions in the future.

Configurations of Parameters
The modeling domain covered East Asia, northern Southeast Asia, and northeastern India, with a 27 km × 27 km grid resolution. A Lambert projection with a domain origin of 34 • N, 110 • E, and two true latitudes of 25 • N and 40 • N was used [33]. The meteorological field simulation employed 35 vertical layers, while the chemical field employed 14 layers (from ground to 100 hPa). As in previous studies, a five-day simulation spin-up was performed to minimize the effects of initial conditions [3,34,35]. The concentrations in the ground layer were analyzed. As shown in Figure 1, the region-specific analysis covered five typical regions: the Jin-Jin-Ji region (JJJ), the Yangtze River Delta (YRD), the Pearl River Delta (PRD), the Sichuan Basin (SCH), and central China (HUZ) [36]. Data from eastern China including five typical regions (ECH) and the entire Mainland China (CHA) were used for specific calculations to analyze overall trends. The main parameter settings of the WRF-CMAQ model were the same as those of Liu et al. [36].  Based on the downscaled LUH2 data and the evolution of land-use types considering China's "carbon neutrality" background, we quantified the emission terrestrial natural sources under two SSP meteorological pathways, as well as co tions to O3 and PM (hereafter, PM refers to PM2.5 and PM10 in this study) concent The inventory of anthropogenic source emissions was kept consistent with the ba and the meteorological field used the 2015 baseline scenario and the future scen SSP126 and SSP585. In the scenario named "carbon", we considered China's "carb trality" background, and its meteorological field followed the SSP126 scenario, wh closer to low-carbon evolution. Table 2 lists the specific simulation cases. The dif in each group represents the impact of land-use change in each scenario in which bon group is compared with the SSP126 group. For example, the difference b ssp126 and ssp126g represents the impact of land-use change under the SSP126 change pathway.

Model Performance
The meteorological predictions were compared with data from the National C Based on the downscaled LUH2 data and the evolution of land-use types in 2050 considering China's "carbon neutrality" background, we quantified the emissions from terrestrial natural sources under two SSP meteorological pathways, as well as contributions to O 3 and PM (hereafter, PM refers to PM 2.5 and PM 10 in this study) concentrations. The inventory of anthropogenic source emissions was kept consistent with the base year, and the meteorological field used the 2015 baseline scenario and the future scenarios of SSP126 and SSP585. In the scenario named "carbon", we considered China's "carbon neutrality" background, and its meteorological field followed the SSP126 scenario, which was closer to low-carbon evolution. Table 2 lists the specific simulation cases. The difference in each group represents the impact of land-use change in each scenario in which the carbon group is compared with the SSP126 group. For example, the difference between ssp126 and ssp126g represents the impact of land-use change under the SSP126 climate change pathway.

Model Performance
The meteorological predictions were compared with data from the National Climatic Data Center (NCDC, https://www.ncdc.noaa.gov/data-access/land-based-station-data/, accessed on 11 June 2022). Surface observations from the China National Environmental Monitoring Center (http://beijingair.sinaapp.com/, accessed on 11 June 2022) were used to validate air quality predictions. We compared the simulated temperature and humidity levels at 2 m (T2 and Q2), as well as wind speed and wind direction levels at 10 m (WS10 and WD10), with observations from 2015. The statistical indicators used include mean observation (Mean OBS), mean simulation (Mean SIM), bias, gross error (GE), root mean square error (RMSE), index of agreement (IOA), and normalized mean bias (NMB). The benchmark values were derived based on Emery et al. [37]. We also compared the simulated maximum daily average of 8 h O 3 (MDA 8 h O 3 ) and PM 2.5 from CMAQ with observation data. The error level of the WRF and CMAQ models was considered to be within acceptable limits for this study, as described in Liu et al. [36].  T benchmark values were derived based on Emery et al. [37]. We also compared the sim lated maximum daily average of 8 h O3 (MDA 8 h O3) and PM2.5 from CMAQ with obs vation data. The error level of the WRF and CMAQ models was considered to be with acceptable limits for this study, as described in Liu et al. [36]. Figure 2 depicts the evolution of China from 2015 to 2100 under the SSP126 a SSP585 scenarios. In general, the forecast trends of forest land use under the two scenar were consistent, i.e., primary forest land decreased year after year, while secondary for land increased year after year, but the SSP126 scenario was more effective in protecti the forest and other vegetation regions. Furthermore, when compared with the SSP5 scenario, the proportion of pastures and range land decreased significantly under t SSP126 scenario. scenarios (primary forest, secondary forest land, managed pastures, other pastures, C3 ann crops, C3 nitrogen-fixing crops, C3 perennial crops, C4 annual crops, C4 perennial crops, urb primary non-forest land, and secondary non-forest vegetation land, denoted as primf, secdf, pa range, c3ann, c3nfx, c3per, c4ann, c4per, urban, primn, and secdn, respectively).

Land-Use Changes of China from 2015 to 2050 under SSP Scenarios
The differences in the evolution of land-use types under the SSP126 and SSP585 s narios were determined by the constraints used. With the SSP126 scenario, we aimed forecast future scenarios based on the green growth paradigm, comprising moderate po ulation growth, economic growth, technological progress, and increased agricultu productivity. Furthermore, this scenario depicted a world concerned with limiting bio versity loss and reducing the demand for animal products. This scenario considered re vant mitigation policies in order to limit warming to no more than 2 °C. Important polic in terms of global land-use change parameters include increasing the use of bioener combined with carbon capture and storage to reduce deforestation and restore degrad forests. In this scenario, the combination of socioeconomic trends and climate policies sulted in a significant reduction in total agricultural land. Simultaneously, large are were being dedicated to bioenergy production, and the forest area expanded. For exa ple, Figure 3 shows the changes in secdf (secondary forest land) distribution from 2015 2050. Increased secdf in southern China may increase BVOC emissions in correspondi regions. scenarios (primary forest, secondary forest land, managed pastures, other pastures, C3 annual crops, C3 nitrogen-fixing crops, C3 perennial crops, C4 annual crops, C4 perennial crops, urban, primary non-forest land, and secondary non-forest vegetation land, denoted as primf, secdf, pastr, range, c3ann, c3nfx, c3per, c4ann, c4per, urban, primn, and secdn, respectively).
The differences in the evolution of land-use types under the SSP126 and SSP585 scenarios were determined by the constraints used. With the SSP126 scenario, we aimed to forecast future scenarios based on the green growth paradigm, comprising moderate population growth, economic growth, technological progress, and increased agricultural productivity. Furthermore, this scenario depicted a world concerned with limiting biodiversity loss and reducing the demand for animal products. This scenario considered relevant mitigation policies in order to limit warming to no more than 2 • C. Important policies in terms of global land-use change parameters include increasing the use of bioenergy combined with carbon capture and storage to reduce deforestation and restore degraded forests. In this scenario, the combination of socioeconomic trends and climate policies resulted in a significant reduction in total agricultural land. Simultaneously, large areas were being dedicated to bioenergy production, and the forest area expanded. For example, Figure 3 shows the changes in secdf (secondary forest land) distribution from 2015 to 2050. Increased secdf in southern China may increase BVOC emissions in corresponding regions.  The SSP585 was distinguished by rapid and resource-intensive development that c responds to material-intensive consumption patterns, as well as technological advanc including agricultural productivity. In the SSP585 scenario, which represented the high level of emissions, we considered very high levels of fossil fuel use, doubling global fo demand, and tripling GHG emissions by 2050. Under the SSP585 scenario, there would a strong expansion of global cropland into pastures and woodlands, with a global incre of about 300 Mha (20%), a trend that also applied to the downscaling results in China.

Changes in BVOC Emissions Driven by Land-Use Changes
Land-use changes directly affect emissions from terrestrial natural sources. Increas vegetation coverage, for example, results in a significant increase in BVOC emissions. examined the changes in terrestrial natural source emissions driven by land-use chan under SSP scenarios from 2015 to 2050. The findings revealed that, after changing lan use type, BVOC emissions changed significantly, whereas changes in soil-source ni oxide (NO) and dust emissions were still controlled by meteorology. This is because land-use change data used in this study primarily concern forests, pastures, crops, a other vegetation-related types but do not include the prediction of refined soil types, m ing BVOC the most sensitive to such changes. Figure 4 depicts the monthly variations soil-source NO and BVOC emissions in China under various scenarios. The soil-sou NO emissions from the five simulation sets showed only three magnitudes that cor sponded to their respective input meteorological data, indicating that there were no s nificant differences between simulations using the same set of meteorological data a different land-use data. BVOC emissions, on the other hand, produced different resu For example, the magnitudes of ssp126 and ssp126g simulation results using the same The SSP585 was distinguished by rapid and resource-intensive development that corresponds to material-intensive consumption patterns, as well as technological advances, including agricultural productivity. In the SSP585 scenario, which represented the highest level of emissions, we considered very high levels of fossil fuel use, doubling global food demand, and tripling GHG emissions by 2050. Under the SSP585 scenario, there would be a strong expansion of global cropland into pastures and woodlands, with a global increase of about 300 Mha (20%), a trend that also applied to the downscaling results in China.

Changes in BVOC Emissions Driven by Land-Use Changes
Land-use changes directly affect emissions from terrestrial natural sources. Increased vegetation coverage, for example, results in a significant increase in BVOC emissions. We examined the changes in terrestrial natural source emissions driven by land-use change under SSP scenarios from 2015 to 2050. The findings revealed that, after changing land-use type, BVOC emissions changed significantly, whereas changes in soil-source nitric oxide (NO) and dust emissions were still controlled by meteorology. This is because the landuse change data used in this study primarily concern forests, pastures, crops, and other vegetation-related types but do not include the prediction of refined soil types, making BVOC the most sensitive to such changes. Figure 4 depicts the monthly variations in soil-source NO and BVOC emissions in China under various scenarios. The soil-source NO emissions from the five simulation sets showed only three magnitudes that corresponded to their respective input meteorological data, indicating that there were no significant differences between simulations using the same set of meteorological data and different land-use data. BVOC emissions, on the other hand, produced different results. For example, the magnitudes of ssp126 and ssp126g simulation results using the same set of meteorological data but different land data were not the same. In general, the simulated dust emission fluxes of different land-use scenarios under the SSP126 and SSP585 pathways differed by less than 1%, and the annual total soil-source NO emissions were 15.3 Mmol and 17.0 Mmol, respectively, and the changing ratio affected by land-use was 1.3% of the absolute value of emission, which is negligible compared with the proportion of changes caused by meteorological differences. Therefore, we focused on analyzing the response of BVOC in this study.   (18.4 TgC) in 2050, compared with 2015, which was compara to the SSP126 pathway. This is because the increase in forest coverage under the SSP land-use change scenario was greater than that under the SSP585 scenario; therefore, EF corresponding to the BVOC was greater. In comparison, under the SSP585 scena temperature and humidity increased more rather than SSP126, leading to more BVOC crease due to meteorological changes than under the SSP126 scenario. However, BVOC emission of ssp585g in summer was about 1.5% higher than that of ssp126g, du the more significant climate change in summer under the SSP585 scenario. Overall, BVOC emissions of each scenario can be divided into four levels based on their tempo variation as in Figure 4b: (1) 2015, (2) SSP126 without land-use changes, (3) SSP126 w land changes and SSP585 without land-use changes, and (4) SSP585 with landchanges; the last two scenarios were close in magnitude. This study concluded, that un the SSP126 scenario, the impact of land-use changes on BVOC emissions combined w the impact of meteorological changes could equal or even exceed the impact of only teorological changes under the SSP585 scenario. However, the impact of land-use chan on BVOC emissions under the SSP585 scenario was not obvious. This can be attributed the impact of meteorological changes, which might have already pushed BVOC emissi This is because the increase in forest coverage under the SSP126 land-use change scenario was greater than that under the SSP585 scenario; therefore, the EF corresponding to the BVOC was greater. In comparison, under the SSP585 scenario, temperature and humidity increased more rather than SSP126, leading to more BVOC increase due to meteorological changes than under the SSP126 scenario. However, the BVOC emission of ssp585g in summer was about 1.5% higher than that of ssp126g, due to the more significant climate change in summer under the SSP585 scenario. Overall, the BVOC emissions of each scenario can be divided into four levels based on their temporal variation as in Figure 4b: (1) 2015, (2) SSP126 without land-use changes, (3) SSP126 with land changes and SSP585 without land-use changes, and (4) SSP585 with land-use changes; the last two scenarios were close in magnitude. This study concluded, that under the SSP126 scenario, the impact of land-use changes on BVOC emissions combined with the impact of meteorological changes could equal or even exceed the impact of only meteorological changes under the SSP585 scenario. However, the impact of land-use changes on BVOC emissions under the SSP585 scenario was not obvious. This can be attributed to the impact of meteorological changes, which might have already pushed BVOC emissions to a high level such that the impact of land-use changes was comparatively negligible. In addition, the lower increase in forest land-use type in SSP585, compared with in SSP126 land-use scenario may have also resulted in a lower impact of land-use changes on BVOC emissions in the SSP585 scenario.
Regional differences in the spatial distribution of changes in BVOC emissions were also observed. Figure 5 depicts the impact of land-use changes on BVOC emission rates under the SSP126 and SSP585 pathways, as calculated from ssp126g/ssp585g minus ssp126/ssp585, respectively. Since the land-use changes for each grid were independent, some singularities of abrupt emissions were included, but overall excessive increase and excessive decrease canceled each other out, resulting in average emissions in CHA that were within a reasonable range. Under the SSP126 pathway, land-use changes would reduce BVOC emissions in northern and central China while increasing emissions in southern and northeast China. Such changes are directly related to the evolution of land-use types; for example, the increase in BVOC emissions in the south was linked to the increase in land-use type secdf (secondary forest land) (Figure 3) under the SSP126 scenario. Land-use changes under the SSP585 pathway would reduce BVOC emissions in the south while increasing emissions in the southwest. The BVOC emissions in typical regions were extracted from the impact of land-use changes under two SSP pathways, as listed in Table 3.   Overall, the SSP126 land-use scenario had a greater impact on BVOC emissions creased by 9.0%) than the SSP585 land-use scenario (increased by 1.8%), which matc the evolution pattern of downscaled LUH2 for the two scenarios indicating a higher crease in forest cover under the SSP126 pathway. The structural adjustment of landtypes was more visible in the SSP126 scenario, such as a significant increase in forest la Since the vegetation-related land-use type structure in the SSP585 scenario changed than SSP126, its impact on BVOC emissions was less severe. The downscaling meth could capture the characteristics of land-use type changes in different regions under  Overall, the SSP126 land-use scenario had a greater impact on BVOC emissions (increased by 9.0%) than the SSP585 land-use scenario (increased by 1.8%), which matches the evolution pattern of downscaled LUH2 for the two scenarios indicating a higher increase in forest cover under the SSP126 pathway. The structural adjustment of land-use types was more visible in the SSP126 scenario, such as a significant increase in forest land. Since the vegetation-related land-use type structure in the SSP585 scenario changed less than SSP126, its impact on BVOC emissions was less severe. The downscaling method could capture the characteristics of land-use type changes in different regions under the SSP126 scenario, i.e., forest land increased in southern and northeast China.

Changes in O 3 and PM Concentrations in Response to Land-Use Changes
Land-use changes further affect the concentration of secondary pollutants by causing terrestrial natural source emissions. Since this study was focused on the impact of land-use change by driving terrestrial natural source emissions, the "effects of land-use change on pollutant concentrations" mentioned in the following text are shorthand for the impact of this mechanism. Under both SSP scenarios, land-use changes decreased O 3 concentration by about 0.05 ppbv, increased O 3 MAX concentration by 0.07 ppbv, and decreased O 3 MAX concentration by 0.12 ppbv, compared with no land-use change scenarios, on annual average of China, while changes (increased and decreased grids) under SSP126 scenario were more significant. Land-use changes reduced PM 2.5 and PM 10 concentrations by about 0.05 µg m −3 , especially in northwestern China, where PM 10 concentrations decreased by 0.24 µg m −3 . Figure 6 depicts the annual average concentration changes in O 3 MAX and PM caused by land-use change under the SSP126 and SSP585 scenarios. The impact of land-use changes was calculated by comparing the simulated concentrations produced using the same meteorological data but different land-use data. The variation trend and degree of O 3 and PM concentrations differed due to the varying degree of land-use types in different regions, as listed in Table S2. When compared with the SSP585 scenario, land-use change under the SSP126 scenario had more obvious impacts on O 3 and PM, with a decrease in the north and an increase in the south. This suggests that it was easier to control pollutant concentrations under the SSP126 pathway by optimizing land-use types. Furthermore, the effects of land-use changes on O 3 and PM concentrations were largely consistent with their effects on BVOC emissions, demonstrating that causing BVOC emissions was the primary mechanism by which land-use changes affected O 3 and PM concentrations. Although the overall reduction was still much smaller than the impact of meteorological changes (for example, the impact of meteorological changes on O 3 , O 3 MAX, PM 2.5 , and PM 10 under the SSP126 scenario was about 19.6, 2.2, 14.1, and 8.7 times the impact of land-use changes, respectively), it still demonstrates the feasibility of reducing the contribution of terrestrial natural sources by adjusting land-use type structure.
The impact of land-use changes on O 3 and PM concentrations was significantly weaker when compared with the impact of meteorological changes. In terms of O 3 , the annual average change range of each typical area in the two SSP scenarios was roughly 0.4%. In terms of PM, the average annual change in each typical area under the two SSP scenarios was roughly 1.2%. It should be noted that since this study employed relatively high 2015 anthropogenic emissions for 2050 simulations, the impacts of land-use changes may be much higher, considering the significant future anthropogenic emission reductions. Since terrestrial natural source emissions are concentrated in the summer, when land-use type changes, the impact on O 3 and PM concentrations are more pronounced in summer due to higher changes in terrestrial natural source emissions. We further extracted the simulation results in typical regions during summer for analysis, as listed in Table S3. Land-use change had a greater impact on the O 3 and PM concentrations in summer, and the maximum of summer average change range of the SSP scenario in typical regions expanded to 0.7% and 1.6%, respectively.
The impact of the land-use changes in China on O 3 and O 3 MAX concentrations in summer was approximately twice the annual average level, while the impact on PM 2.5 and PM 10 was approximately 1.5 times on average of two SSP scenarios. The specific results of each scenario are listed in Table 4.
BVOC emissions was the primary mechanism by which land-use changes affected O3 and PM concentrations. Although the overall reduction was still much smaller than the impact of meteorological changes (for example, the impact of meteorological changes on O3, O3MAX, PM2.5, and PM10 under the SSP126 scenario was about 19.6, 2.2, 14.1, and 8.7 times the impact of land-use changes, respectively), it still demonstrates the feasibility of reducing the contribution of terrestrial natural sources by adjusting land-use type structure.

Projections under Carbon-Neutral Scenario of China
The land-use change under the carbon scenario was similar to that of SSP126 in terms of vegetation changes. Assuming that forest and cultivated land area conform to national policy, fossil energy would be reduced, biomass ratio would be increased, and the trend of urban development would increase and then decrease. Based on the Global Change Assessment Model (GCAM), we adjusted the land-use type to achieve a peak in carbon emissions in 2030 and carbon neutrality in 2060, as shown in Figure 7. The classification method for land-use types under this scenario differed from that under the SSP scenarios, including biomass (land used for planting energy crops including sugar and oil crops), crops, deserts, forests, grassland, other arable lands, pasture, shrubs, tundra, and urban, which are expressed as biomass, crops, desert, forest, grass, other arable, pasture, shrubs, tundra, and urban in turn, from which the vegetation-related land-use change proportions were extracted to correspond to the SSP scenarios. As shown in Figure 8, BVOC emissions in all regions would increase under the carbon scenario when compared with the ssp126 case without accounting for land-use change. Table 5 shows the changes in total BVOC emissions caused by land-use changes in typical regions. Meteorological changes under the SSP126 pathway would increase BVOC emissions in China by 1.0 TgC in 2050, compared with 2015, whereas land-use changes under the carbon scenario would increase them by 7.1 TgC, which is 6.8 times the meteorological impact. On the one hand, under the carbon scenario used in this study, a large number of carbon sinks were added by increasing vegetation coverage in order to achieve the goals of "carbon peaking" and "carbon neutrality," and the vegetation growth estimation was relatively aggressive. On the other hand, if vegetation coverage is greatly increased in the future, BVOC emissions will continue to rise and cannot be ignored.  Since for the carbon scenario in this study, we used the same meteorological data as those of the SSP126 scenario, the simulated soil-source NO emissions (total 16.1 Mmol) were the same as the results of scenarios using the SSP126 meteorological data but were still 5.3% higher than those of the scenarios that did not consider land-use changes. The annual BVOC emission under the carbon scenario was 24.0 TgC, exceeding the results of all other scenarios, and was 1.4 times higher than that of the SSP126 scenario without considering land-use changes. As a result, land-use changes revealed a higher potential to increase BVOC emissions under the SSP126 pathway, compared with the SSP585 pathway. The monthly variation trends of BVOC and soil-source NO emissions under the carbon scenario were consistent with other scenarios, and the specific emissions are listed in Table S4.
As shown in Figure 8, BVOC emissions in all regions would increase under the carbon scenario when compared with the ssp126 case without accounting for land-use change. Table 5 shows the changes in total BVOC emissions caused by land-use changes in typical regions. Meteorological changes under the SSP126 pathway would increase BVOC emissions in China by 1.0 TgC in 2050, compared with 2015, whereas land-use changes under the carbon scenario would increase them by 7.1 TgC, which is 6.8 times the meteorological impact. On the one hand, under the carbon scenario used in this study, a large number of carbon sinks were added by increasing vegetation coverage in order to achieve the goals of "carbon peaking" and "carbon neutrality", and the vegetation growth estimation was relatively aggressive. On the other hand, if vegetation coverage is greatly increased in the future, BVOC emissions will continue to rise and cannot be ignored. meteorological impact. On the one hand, under the carbon scenario used in large number of carbon sinks were added by increasing vegetation coverag achieve the goals of "carbon peaking" and "carbon neutrality," and the veget estimation was relatively aggressive. On the other hand, if vegetation cover increased in the future, BVOC emissions will continue to rise and cannot be  Under the carbon scenario, land-use changes increased O 3 and O 3 MAX concentrations by 0.06 ppbv and 0.70 ppbv, respectively, with an average annual change of less than 1.6%, much less than the impact of meteorological changes. The land-use changes resulted in increased PM 2.5 and PM 10 concentrations by 0.16 µg m −3 and 0.28 µg m −3 , respectively, with an average annual change range of less than 3.6%. Figure 9 depicts the impact of land-use change on pollutant concentrations under the carbon scenario, and Table 6 lists specific results in typical regions. Under the carbon scenario, land-use change showed the potential to increase pollutant concentrations. The land-use change had a greater impact on O 3 and PM concentrations during the summer, with maximum ratios of about 3.2% and 6.2%, respectively. The ratios of the impact of land-use changes on the concentration of O 3 , O 3 MAX, PM 2.5 , and PM 10 concentrations in summer and the annual average level were 3.77, 2.68, 2.16, and 1.61, respectively. The PM 2.5 concentration in PRD and HUZ was significantly increased by more than 0.5 µg m −3 in summer. To ensure air quality (the WHO annual PM 2.5 limit of 5 µg m −3 ), potential risks must be mitigated or offset through anthropogenic emission control. of O3, O3MAX, PM2.5, and PM10 concentrations in summer and the annual average lev were 3.77, 2.68, 2.16, and 1.61, respectively. The PM2.5 concentration in PRD and HUZ w significantly increased by more than 0.5 μg m −3 in summer. To ensure air quality (the WH annual PM2.5 limit of 5 μg m −3 ), potential risks must be mitigated or offset through anth pogenic emission control.

Summary
This study predicted three future land-use type scenarios in 2050 (SSP126, SSP585, and carbon scenario) based on downscaled LUH2 data and carbon-neutral policy in China. We quantified the contribution of land-use changes to terrestrial natural source BVOC emissions, as well as O 3 and PM concentrations.
Land-use changes directly affect emissions from terrestrial natural sources. The variation ranges of wind-blown dust and soil-source NO emissions were less than 1.3%, far less than the impact of meteorological changes. BVOC emission changes are sensitive to land-use type. After adjusting the land-use type, the rate of increase in BVOC emissions under the SSP126 pathway could reach the level caused only by meteorological changes under the SSP585 pathway. Under the SSP126 pathway, meteorological changes would increase BVOC emissions in China by 1.0 TgC in 2050 compared with 2015, while land-use changes would increase them by 1.5-7.1 TgC, up to 6.8 times the meteorological impact. If vegetation coverage is significantly increased in the future, the resulting increase in BVOC emissions cannot be ignored. The SSP126 land-use scenario had a greater impact on BVOC emissions than the SSP585 land-use scenario, which was consistent with the evolution characteristics obtained after downscaling the two LUH2 scenarios.
The impact of land-use changes on O 3 and PM concentrations (by causing terrestrial natural source emissions) was essentially the same trend as their impact on BVOC emissions, which was greater in summer (about twice the annual average level), demonstrating that the driving force behind BVOC emissions is the impact mechanism of land-use changes on O 3 and PM concentrations. The land-use change caused a significant decrease in PM 10 concentration in northwestern China under the SSP126 scenario, demonstrating the feasibility of reducing the contribution of terrestrial natural sources by adjusting the structure of land-use types. Overall, land-use changes had a weaker impact on O 3 and PM concentrations than meteorological changes. The average annual change in each scenario was within a 3.6% limit and larger in summer (within 6.2%). Regional differences must be considered when calculating future pollutant background concentrations in the environment. The impact of land-use change on pollutants was more obvious under the SSP126 pathway than under the SSP585 pathway, indicating that more attention should be paid to optimizing land-use types under the SSP126 pathway. Under the carbon scenario, the increase in PM concentration caused by terrestrial natural sources as a result of landuse change would pose a risk to air quality compliance; thus, it is necessary to consider reducing or offsetting this potential risk through anthropogenic emission control measures.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/atmos13071086/s1, Table S1: Fitting results of the relationship between FPAR and LAI; Table S2: Impact of land−use changes on annual average O 3 (unit: ppbv) and PM (unit: µg m −3 ) concentra-tions in typical regions under SSP scenarios; Table S3