Measurement report: Unexpected high volatile organic compounds emission from vehicles on the Tibetan Plateau

Vehicular emissions significantly impact atmospheric composition, yet their behavior in high-altitude environments remains poorly understood. Here, we present a comprehensive investigation of vehicular emissions conducted across ten tunnels on the Tibetan Plateau, encompassing nearly 3000 m of altitudinal variation. Our findings reveal that total emission factors of volatile organic compounds (VOCs) increase with altitude, driven primarily by a pronounced rise in evaporative emissions. Source apportionment indicates that evaporative emissions contribute 67 % of the total VOCs, markedly outweighing the 24 % from tailpipe exhaust, a pattern that diverges from those observed at lower altitudes, where tailpipe emissions dominate. We identify low atmospheric pressure as the key factor enhancing evaporative VOC release at high altitudes. This study offers critical new insights into the mechanisms of vehicular emissions in plateau regions, underscoring the necessity of accounting for specific environmental conditions and fuel evaporation when formulating emission control strategies for high-altitude settings.

1 Introduction

The Tibetan Plateau, with an average altitude of over 4000 m and a total area of approximately 2.5×10⁶ km², is known as the “Roof of the World” (Xu et al., 2014). The low pressure and cold climate have resulted in low levels of population and industrialization in the Tibetan Plateau, making it one of the most pristine regions in the world (Xue et al., 2013; Li et al., 2017). However, with the development of infrastructure and tourism, the vehicle ownership in Tibet has notably increased in recent years, leading to a rapid growth in vehicular emissions (Tang et al., 2022; Zhong et al., 2023). Yet, the characteristics of vehicular emissions on the Tibetan Plateau – both in terms of emission strengths and compositions, remain largely unknown.

Vehicular emissions comprise tailpipe exhausts and non-tailpipe emissions (Wu et al., 2023; Wang et al., 2022a; Hata et al., 2019; Jathar et al., 2014; Du et al., 2018). The former refers to gases emitted from engine systems due to incomplete combustion or unburned fuel (Zhang et al., 2024), while the latter mainly consists of evaporative emissions from fuels and fuel-related additives, released through processes such as running loss, hot soak, permeation, and diurnal breathing (Zhang et al., 2024; Pierson et al., 1999; Yue et al., 2017; Man et al., 2020; Liu et al., 2015; Harrison et al., 2021). In terms of volatile organic compounds (VOCs), the compositions of the two types of emissions differs significantly. With the upgrading engine and after-treatment technologies, tailpipe exhausts have been effectively reduced in recent decades (Cliff et al., 2023; Wang et al., 2022a). As a consequence, the relative contributions of non-tailpipe emissions are increasing, accounting for 32 % to 48 % of total vehicular VOC emissions (Liu et al., 2017; Wu et al., 2023; Chai et al., 2023; Huo et al., 2024; Sun et al., 2021).

The evaporation of fuels and solvents is an equilibrium process involving hydrocarbon molecules transitioning between the gas and liquid phases, governed primarily by temperature and pressure conditions (Huang et al., 2022). Elevated temperature and low pressures theoretically facilitate fuel evaporation (Huo et al., 2024; Wang et al., 2018). However, in plateau regions, low atmospheric pressure often coincides with cooler temperatures, resulting in competing influences that are not yet fully quantified. Computational models such as MOVES model estimate that tank vapor generation at 1700 m a.s.l. (above sea level) is approximately 1.4 times that at sea level, indicating potentially significant altitude-enhanced emissions (US EPA, 2024). They also underscore the critical need for empirical validation under real-world, high-altitude conditions, which remain severely limited.

Within the framework of the second scientific expedition and research program of the Tibetan plateau (STEP) (Yao et al., 2012; Ye et al., 2023), we conducted vehicular emission measurements in 10 tunnels across the Tibetan Plateau, spanning an altitude range of nearly 3000 m. This unique natural setting enabled us to investigate how vehicular emission characteristics respond to changing elevation, with a particular emphasis on the role of low atmospheric pressure. The study aims to enhance the current understanding of VOC emissions from vehicular in high-altitude regions, which remain poorly characterized in existing literature.

2 Materials and methods

2.1 Research site

Following the criteria of representative altitude, we specifically chose ten tunnels located between Lhasa and Nyingchi, two major cities in Tibet autonomous region, China (Fig. S1 in the Supplement). We prioritized selecting one-way tunnels, as well as the longest available tunnels. These tunnels cover a distance of approximately 700 km and span an altitude range of nearly 3000 m. By merging tunnels with similar altitudes (height difference less than 100 m), we categorized the ten tunnels into six groups: 2000, 3000, 3400, 3600, 4200, and 4750 m. Detailed information about the tunnels is summarized in Table S1 in the Supplement.

2.2 Mobile measurement system

A vehicle equipped with a CO₂/H₂O gas analyzer (850A, LiCOR, Inc., USA), an infrared absorption CO analyzer (48i-TL, Thermo, Inc., USA), and summa canisters (Entech Instrument, Inc., USA, 3.2 L) as offline VOCs sampler, served as the mobile measurement system for this study (Fig. S2). Experiments were conducted with the vehicle traveling along the right side of the road inside the tunnel at relatively constant speeds, ranging from 15 to 40 km h⁻¹, depending on the length of the tunnels. The sampling tubing was positioned on the left side of the mobile measurement vehicle to capture mid-tunnel conditions. The experimental process was documented using both driving recorder footage and manual record-keeping (Fig. S2).

The online instruments were pre-calibrated to minimize random errors. The VOC samples, collected over a 10 d sampling campaign across multiple tunnels, were analyzed within one week after transportation to our laboratory using a gas chromatography-flame ionization detection/mass spectrometry (GC-FID/MS) system (Sect. S1 in the Supplement), ensuring minimal pre-analysis storage time. In the one-way tunnels, the online data (i.e., CO₂ and CO) showed a noticeable piston effect (Fig. S3) (Chung and Chung, 2007), with concentrations gradually increasing towards the end of the tunnel. The air at the tunnel's tail end was assumed to represent a well-mixed plume from emissions of all vehicles in the tunnel (Hwang et al., 2023; Gillies et al., 2001). Therefore, in these tunnels, offline sampling was initiated in the rear section and lasted approximately 1 min to capture the accumulated air masses. Additionally, three tunnels in our study had bidirectional traffic, where the piston effect was less pronounced due to opposing flows. For these cases, sampling was conducted at the tunnel midpoint to ensure representative mixing of emissions from both directions. Background concentrations of VOCs were determined at the Yangbajing background site during the same field campaign of STEP (July–August 2022) (Tao et al., 2024).

Each tunnel underwent either parallel or bidirectional sampling, with a minimum of two tests conducted, resulted in a total of 46 valid samples. A total of 90 VOCs were quantitatively identified and classified into six categories: alkanes (29 species), alkenes (11 species), alkynes (1 species), aromatics (14 species), halocarbons (27 species), and oxygenated volatile organic compounds (OVOCs) (8 species) (Table S2). As illustrated in Figs. S4 and S5, the parallel samples from the same tunnels showed good consistency.

2.3 Emission factors and emission ratios

To compare the vehicular emission intensities under different conditions with a standardized scale, we calculated the fuel-based emission factors (EF) and emission ratios (ER) of VOC relative to CO, respectively.

$$\begin{array}{}\text{(1)}& {\displaystyle }{\mathrm{EF}}_i={\displaystyle \frac{\sum C_i\cdot C_{\mathrm{F}}\cdot \mathrm{1000}}{\left(\frac{\sum C_{{\mathrm{CO}}_{\mathrm{2}}}}{{\mathrm{MW}}_{{\mathrm{CO}}_{\mathrm{2}}}}+\frac{\sum C_{\mathrm{CO}}}{{\mathrm{MW}}_{\mathrm{CO}}}\right)\cdot {\mathrm{MW}}_{\mathrm{C}}}}\text{(2)}& {\displaystyle }{\mathrm{ER}}_i={\displaystyle \frac{\sum C_i}{\sum C_{\mathrm{CO}}}}\end{array}$$

In Eqs. (1) and (2), EF_i is the emission factor of VOC species i, µg kg${}_{\mathrm{fuel}}^{-\mathrm{1}}$, and ER_i is the emission ratio of VOC species i, ppb ppm⁻¹; C_i, $C_{{\mathrm{CO}}_{\mathrm{2}}}$ and C_CO are the concentration of VOC species i, CO₂ and CO respectively; C_F is the carbon mass fraction of the fuel, and a typical value of 0.86 used here (Wren et al., 2018). ${\mathrm{MW}}_{{\mathrm{CO}}_{\mathrm{2}}}$, MW_CO and MW_C are the molecular weights of CO₂, CO, and carbon atom respectively, g mol⁻¹.

2.4 Source profiles similarity analysis

An approach commonly used to assess the resemblance between mass spectrometry signals, known as vector similarity analysis (Humes et al., 2022; Zhu et al., 2021; Kostenidou et al., 2009; Ulbrich et al., 2009; Wang et al., 2024), is adopted here to evaluate the similarity among various source profiles. Specifically, relative molecular weight of each VOC species multiplied by its corresponding mass fraction is treated as a dimension in the vector. The cosine angle (θ) between the vector of two source profiles (x and y) is then calculated using Eq. (3) to determine their similarity. The θ angles are categorized into four groups: 0–15, 15–30, 30–50, and >50°, representing excellent consistency, good consistency, many similarities, and poor consistency, respectively (Wang et al., 2024).

$$\begin{array}{}\text{(3)}& \cos \mathit{\theta }={\displaystyle \frac{\mathit{x}\cdot \mathit{y}}{\left|\mathit{x}\right|\left|\mathit{y}\right|}}\end{array}$$

2.5 Positive matrix factorization and non-negative least squares analysis

The positive matrix factorization (PMF) model 5.0, developed by the United States Environmental Protection Agency (US EPA) (Paatero and Tapper, 1994), was used to identify sources of VOCs sampled from the tunnels. Three factors representing different sources of VOCs are extracted from the PMF model: evaporative emissions, tailpipe exhaust, and non-vehicular sources. A total of 66 measured VOC species were included in the PMF analysis, comprising high-concentration species and typical tracers of evaporative and combustion sources. As such, PMF factors can be identified with those well-known and long-lived emission tracers, such as methyl tertiary-butyl ether (MTBE) as a tracer for gasoline evaporation (Kirchstetter et al., 1999), alongside interpreting the sources of other species based on their covariance with those well-known tracers and their presence in PMF factors.

To corroborate the findings from the PMF analysis, a non-negative least squares (NNLS) multivariate linear regression was utilized to examine the relative contributions of evaporative and exhaust emissions (Wang et al., 2022b). Using characteristic indicators of evaporative and tailpipe emissions as independent variables, the optimal proportions of contributions from the three sources are determined to minimize the sum of squared residuals of the objective function. Detailed information can be found in Sect. S2 and Table S3. The data associated with this research can be accessed via the Zenodo repository (Huang, 2024).

3 Results and discussions

3.1 Vehicular VOC emissions increase in high-altitude tunnels

For cross-study comparison, a consistent subset of 59 VOC species, commonly detected in both our study and low-altitude tunnel studies (e.g., Ho et al., 2009), was selected to ensure comparability in EF and ER calculations. The total EF and ER of the detected VOCs generally exhibit higher values at increased altitudes (Figs. 1a and S6), with the averages across all plateau tunnels of 3.3 ± 3.1 µg kg${}_{\mathrm{fuel}}^{-\mathrm{1}}$ and 87 ± 92 ppb ppm⁻¹, respectively. When comparing the identical compositions, the average EF in plateau tunnels is 1.9 times higher than that in the Shing Mun tunnel in Hong Kong (50 m a.s.l.) (Ho et al., 2009). The determined ER is 3.9 and 1.9 times higher than those in the Shing Mun tunnel and the Chung-Liao tunnel in Taiwan (330 m a.s.l.) (Chiang et al., 2007), respectively. These results are significantly higher than the 1.4 times increase in fuel evaporation observed in Denver (1700 m a.s.l.) relative to sea level (US EPA, 2024).

Figure 1 Altitudinal trends of different VOC components. (a) Stacked plot illustrating variations of EF with altitude. The solid line in (a) represents a subset of 59 VOC species that overlapped with species reported in key low-altitude tunnel studies (Ho et al., 2009; Chiang et al., 2007). The gray dashed line on the left side represents results from other studies conducted in low-altitude tunnels. Specifically, data at 50 m from the Shing Mun tunnel in Hong Kong (Ho et al., 2009). (b) Altitudinal variations in source profiles (relative mass percentage contribution) observed in this study.

Alkanes are the predominant components of VOCs in plateau tunnels (Fig. 1b), with an average mass fraction of 65 ± 17 %, significantly higher than those reported in low-altitude studies, reporting the range of alkanes typically from 23 % to 57 % (Staehelin et al., 1998; Hwa et al., 2002; Zhang et al., 2018b; Chiang et al., 2007; Gentner et al., 2013; Sun et al., 2019; Feng et al., 2021; Song et al., 2020; Jin et al., 2021; Song et al., 2018). Meanwhile, the proportion of alkanes shows an increasing trend with altitude, growing from 60 % at 2000 m to 77 % at 4750 m. The tunnels at 3000 m deviates from this trend, potentially due to interference from other sources, which will be discussed later.

Butanes (n/i-butane) and pentanes (n/i-pentane) are identified as the major species in plateau tunnels, collectively contributing to 46 ± 20 % of the total VOCs mass and 67 ± 19 % of the total alkane emissions. This high proportion significantly exceeds those typically found in low-altitude tunnels, ranging from 3.0 % to 21 % (Staehelin et al., 1998; Hwa et al., 2002; Zhang et al., 2018b; Chiang et al., 2007; Gentner et al., 2013; Sun et al., 2019; Feng et al., 2021; Song et al., 2020; Jin et al., 2021; Song et al., 2018). Furthermore, their combined fraction exhibited an upward trend with rising altitude, consistent with the trend observed for alkanes. Clearly, the increased vehicular VOC emissions on the plateau are primarily attributable to the increase in alkanes, specifically butanes and pentanes.

3.2 Low atmospheric pressure significantly enhances evaporative emissions

Butanes and pentanes are primary constituents of fuel evaporation (Chin and Batterman, 2012; Wu et al., 2023; Zhang et al., 2021). As depicted in Fig. 2a, the EF of butanes and pentanes both exhibit distinct enhancement with higher altitude, consistent with their proportional trends. Notably, the average EF of butanes and pentanes in plateau tunnels are 4.5 to 21 times higher than those at 50 m altitude, respectively. Similar results can also be obtained when expressed in terms of the ER (Fig. S7a). To the best of our knowledge, this is the first report of enhanced evaporative emissions from vehicles in high-altitude regions.

Figure 2 Variation trends of EF for representative species with altitude. (a) Trends in EF for butanes and pentanes. (b) Trends in EF for ethylene and benzene. Data at 50 m is from the Shing Mun tunnel in Hong Kong (Ho et al., 2009).

As altitude increases, atmospheric pressure decreases, with pressure at 4750 m altitude being only about 55 % of sea level pressure. According to Dalton evaporation formula (Hage, 1975), the rate of liquid evaporation is inversely proportional to atmospheric pressure. In addition, based on the Clausius-Clapeyron equation (Koutsoyiannis, 2012), lower pressure decreases the boiling point of a solvent. Four-carbon and five-carbon compounds could be more sensitive to this effect, as their boiling points typically lie close to ambient temperatures under standard pressure. For example, as atmospheric pressure decreases from 101 to 56 kPa (corresponding to altitudes from sea level to 4750 m), the boiling point of n-pentane decreases from 35 to 19° (Sect. S3). This means that the n-pentane would reach its boiling point at 4750 m altitude at ambient temperature, which is 19–22° for daytime during our measurement campaign. Interestingly, our tests have confirmed these thermodynamic theories despite the discrepancies between fuel and ideal solutions. While lower temperatures in high-altitude regions are expected to reduce evaporation (Huang et al., 2022), our real-world tests indicate that the influence of low pressure completely outweighs the temperature effect on the Tibetan Plateau.

In contrast, the EF of certain tracers (i.e., ethylene, benzene) mainly from tailpipe exhausts exhibit a fluctuating distribution pattern at various altitudes (Figs. 2b and S7b), with some trend for even higher emissions at 50 m altitude compared to those on the plateau. Moreover, the $\mathrm{CO}/{\mathrm{CO}}_{\mathrm{2}}$ ratio, as an indicator of engine combustion efficiency (Vollmer et al., 2007; Ammoura et al., 2014; Hu et al., 2019), did not appear a discernible altitude-specific distribution across our dataset (Fig. S8), with average values ranging from 5.1 to 11 ppbv ppmv⁻¹. This range is comparable to values reported in low-altitude tunnel studies, spanning 4.4–15 ppbv ppmv⁻¹ across cities in Asia and Europe (Cui et al., 2018; Liu et al., 2022; Ammoura et al., 2014; Legreid et al., 2007). Although reduced oxygen at higher elevations may theoretically impair combustion efficiency and increase $\mathrm{CO}/{\mathrm{CO}}_{\mathrm{2}}$ ratios, other factors, such as vehicle type, engine operation conditions, and tunnel ventilation may obscure the effect of altitude alone.

3.3 Source profile of vehicular VOC emissions on the Tibetan Plateau

Although the vehicle types and quantities vary among the tunnels, the VOCs source profiles at six altitudes demonstrate remarkable consistency, except at 3000 m, where aromatics and OVOCs are relatively higher compared to other altitudes (Fig. S9). Quantitative similarity analysis also indicates that the θ angles between tunnels excluding 3000 m are less than 30°, indicating excellent or good consistency (Table S4). The tunnel at 3000 m is located in the urban areas of Nyngchi and is much shorter than the others (Table S1). As shown in discussion of source analysis below, these two factors both make the collected tunnel samples more susceptible to influenced by surrounding urban sources, hence exhibiting the poorest similarity.

Figure 3 Comparison of source profiles between this study and other studies. (a) Comparison of the tunnel measurement source profile in this study with the gasoline vapor source profile tested by SHED and the source profile of the low-altitude tunnel. (b) Comparison of source profiles for evaporative emissions derived from PMF, gasoline vapor emissions from SHED test results, and directly measured source profile in this study. (c) Comparison of source profiles for tailpipe exhaust derived from PMF, gasoline vehicle exhausts from chassis dynamometer test, and the low-altitude tunnel source profile. The background colors of red, yellow, blue, and purple represent alkanes, alkenes, aromatics, and OVOCs, respectively.

The average source profile of plateau tunnels significantly differs from those documented for low-altitude tunnels (Figs. 3a, S10) (Staehelin et al., 1998; Hwa et al., 2002; Zhang et al., 2018b; Chiang et al., 2007; Gentner et al., 2013; Zhang et al., 2018a; Sun et al., 2019; Feng et al., 2021; Song et al., 2020; Jin et al., 2021; Song et al., 2018). The former exhibit much higher contributions from evaporative-related species (such as butanes, pentanes, and MTBE), while the proportions of combustion-related VOC species (such as ethylene, acetylene, and aromatics) are substantially lower compared to the latter. Moreover, the source profiles at high altitudes align closely with the sealed housing for evaporative determination (SHED) tests (Harley et al., 2000; Na et al., 2004; Liu et al., 2008; Zhang et al., 2013; Wu et al., 2015; Man et al., 2020; Sun et al., 2021), whereas those at low altitudes correlate well with vehicle exhaust from chassis dynamometer test (Schauer et al., 2002; Na et al., 2004; Guo et al., 2011; Gao et al., 2012; Ou et al., 2014; Li et al., 2019; Wang et al., 2022a). This observation is further supported by the similarity analysis (Table S5). These results demonstrate that vehicular VOC emissions on the Tibetan Plateau are dominated by evaporative emissions, which significantly differ from low-altitude tunnels, where tailpipe exhaust prevail.

Table 1 The θ angles (°) among the source profiles of this study and those derived from SHED tests, gasoline vehicle exhaust from dynamometer tests, and low-altitude tunnel measurements. Direct measurement refers to the average VOC source profile based on all tunnel samples measured in this study.

^a Average profiles of gasoline vapors derived from SHED tests (Harley et al., 2000; Na et al., 2004; Liu et al., 2008; Zhang et al., 2013; Wu et al., 2015; Man et al., 2020; Sun et al., 2021). ^b Average profiles of gasoline vehicle exhaust derived from dynamometer tests (Schauer et al., 2002; Na et al., 2004; Guo et al., 2011; Gao et al., 2012; Ou et al., 2014; Li et al., 2019; Wang et al., 2022a). ^c Average profiles of vehicular emissions derived from low-altitude tunnel measurements (Staehelin et al., 1998; Hwa et al., 2002; Zhang et al., 2018b; Chiang et al., 2007; Gentner et al., 2013; Zhang et al., 2018a; Sun et al., 2019; Feng et al., 2021; Song et al., 2020; Jin et al., 2021; Song et al., 2018).

3.4 Source apportionment of VOCs in plateau tunnels

The factors derived from the PMF model are depicted in Fig. S11. One of these factors exhibits remarkable consistency with the composition of laboratory-measured gasoline vapor, with the angle between them being 15° (Fig. 3b, Table 1), thus designated as “evaporative emissions”. Specifically, the source profile of evaporative emissions comprises a majority of C₄–C₈ alkanes, C₄–C₆ alkenes (except for isoprene), and MTBE (Fig. S12, Table S6). In contrast, another factor contains a high proportion of C₂–C₃ alkanes, C₂–C₃ alkenes, C₉–C₁₂ alkanes, aromatics, as well as acetylene, a typical marker of combustion emissions. This agrees with the reported characteristics of tailpipe exhausts in the literature (Wang et al., 2020, 2022a), therefore, this factor is interpreted as representing “tailpipe exhausts”. The angle between the profiles of the PMF-resolved tailpipe exhaust and the chassis dynamometer-tested gasoline vehicle exhaust is 38° (Fig. 3c), showing a relatively poor similarity (Wang et al., 2024), which could be attributed to the influence of diesel vehicles, as well as potential influences from other sources. This discrepancy likely arises because the dynamometer data exclusively represent gasoline vehicle emissions, whereas the tunnel samples include contributions from diesel vehicles, which typically emit higher proportions of heavier alkanes and aromatics (Wang et al., 2022a; Schauer et al., 2002), differing from the VOC composition commonly observed in gasoline vehicle exhaust. The third factor contains substantial amounts of ethanol (59 %), halocarbons (5.4 %–65 %), and isoprene (36 %). Ethanol is unlikely to originate from fuel evaporation or combustion, as there is no ethanol-blended gasoline in Tibet. Additionally, halocarbon species are typically deemed insignificant in vehicular emissions (Kamiya et al., 2016). Isoprene, a tracer of plant emissions (Rasmussen, 1970), was only detected in the tunnels at 2000 m (Fig. S9), which aligns with the dense forests at 2000 m and sparse above 3000 m in Tibet (Fig. S13). It is thus considered to represent non-vehicular sources.

Figure 4 The relative contributions of different sources to VOCs measured in tunnels at various altitudes from PMF analysis. Gray, light red, and light blue bars represent evaporative emissions, tailpipe exhaust, and non-vehicular sources, respectively. The average results for PMF analysis and NNLS analysis are shown in pie charts as well.

As anticipated, evaporative emissions made the largest contributions, account for 67 % of total VOC emissions, whereas tailpipe exhausts and non-vehicular sources contributed to 24 % and 9.0 % (Fig. 4), respectively. For comparison, a tunnel study in Tianjin performed in 2020 shows that evaporative emissions and tailpipe emissions contribute 32 % and 68 % of VOCs, respectively (Sun et al., 2021). Source apportionment in Henan reported that evaporative and tailpipe emissions contributed 38 % and 62 % of vehicular VOC emissions, respectively (Liu et al., 2022). In Haikou, a tropical city, tunnel measurements found that high temperatures significantly enhanced evaporative emissions, leading to the highest reported contribution of evaporative emission (45 %) (Huo et al., 2024). In addition, emission inventory (EI) data for China in 2015 indicated that evaporative emissions and tailpipe exhaust accounted for 39 % and 61 % of vehicular emissions, respectively (Liu et al., 2017). The markedly higher proportion of evaporative emissions in the plateau tunnels reflects that evaporative emission is the most important source of vehicular emissions in high-altitude regions. While variations in vehicle fleet composition, fuel quality, and driving conditions across different cities and time periods can significantly influence VOC emission characteristics in low-altitude studies, our high-altitude tunnel measurements consistently demonstrate systematically elevated evaporative emissions.

The contribution of vehicular evaporative emissions exhibits an upward trend with ascending altitude, rising from 51 % at 2000 m to 72 % at 4750 m (Table S7). The peak at 3600 m (77 %) could be attributed to its relatively higher EF of evaporative species (Fig. 2). In contrast, tailpipe exhaust shows an overall decreasing trend, with contributions decreasing from 36 % to 16 %, although some altitudes deviate somewhat from this trend. It should be note that the decline in the proportion of tailpipe exhaust does not necessarily indicate a decrease in combustion emission intensity but rather be attributed to the overall increase in total vehicular emissions. The contribution of non-vehicular emissions reaches its maximum at 3000 m altitude, consistent with previous descriptions that tunnels at this altitude are more susceptible to impacts from other sources in urban regions.

As a complementary analysis, the results of the NNLS multivariate linear are consistent with those of the PMF analysis, showing the highest contribution from evaporative emissions (62 %), followed by tailpipe exhaust (35 %) and non-vehicular sources (3.0 %) (Fig. 4). Additionally, the representative species for evaporative emissions and tailpipe exhaust from both approaches also show overall consistency (Fig. S12, Table S6). The agreement between the two approaches and their alignment with literature results underscore the reliability of our data and analysis results.

3.5 Implications and perspectives

Over the past two decades (2000–2020), the number of civilian vehicles in Tibet has increased by a factor of 14, while road mileage has expanded fivefold and road freight volume has surged over 30 times, exerting substantial pressure on the region's fragile ecosystem (National Bureau of Statistics of China, 2021). This study demonstrates that low atmospheric pressure significantly amplifies vehicular VOC emissions on the Tibetan Plateau by enhancing evaporative emissions. However, current emission inventories predominantly rely on emission factors derived from low-altitude regions, resulting in a significant underestimation of VOC emissions in high-altitude and mountainous areas. For instance, Liu et al. (2017) estimated that evaporative emissions accounted for only 4 % of total vehicular VOC emissions in Tibet, a figure markedly lower than our experimental results. This discrepancy implies that the inventory of vehicular VOC emissions in Tibet may be underestimated by approximately 70 %. Given that high-altitude lands (above 2000 m a.s.l.) constitute around 32 % of China's total land area and about 11 % of the global land area (Central People's Government of the People's Republic of China, 2005), these findings could have far-reaching implications for regional and even global atmospheric environments and climate.

In recent decades, stringent emission standards, advances in emission control technologies, and improvements in fuel quality have significantly reduced vehicular emissions (Wu et al., 2019; Man et al., 2018). However, these measures primarily focus on tailpipe emissions and may not effectively mitigate VOC emissions from fuel evaporation (Huang et al., 2022). This deficiency, as our study demonstrates, is more pronounced in high-altitude environments. Thus, it is crucial to develop targeted strategies and technologies to mitigate vehicular evaporative emissions, particularly in such regions.

One promising approach for Tibet is the strategic promotion of electric vehicles (EVs), supported by China's mature EVs industry and Tibet's abundance renewable energy resources. The region's installed capacity of hydropower and solar energy exceeds 10×10⁶ kW (National Energy Administration, 2024), with ongoingdevelopment of large-scale solar and wind projects. However, local electricity demand remains low due to sparse population and limited industrial activity, resulting in surplus energy transmitted eastward via the “West–East Electricity Transmission” project (Xinhua News Agency, 2024). Promoting EVs could absorb this surplus, alleviating grid strain. Dispite this, the deployment of EVs in Tibet faces specific challenges, particularly due to the region's low-temperature environments. These conditions necessitate advancements in battery technology, such as the development of solid-state batteries with improved thermal resilience, as well as altitude-adaptive battery management systems. In summary, promotion of EVs in high-altitude regions like Tibet presents a win-win solution: it not only reduces vehicular emissions but also fosters a synergistic relationship with local renewable energy projects, yielding both environmental and economic benefits. This approach could serve as a model for integrating clean energy and transportation solutions in other high-altitude regions globally, contributing to broader efforts to mitigate climate change.

4 Conclusions

This study presents a detailed investigation of vehicular VOC emissions across the Tibetan Plateau, exploring their behavior in a high-altitude region. The findings reveal significant differences in emission characteristics compared to low-altitude regions. Total VOC emission factors and ratios increase significantly with altitude on the Tibetan Plateau, driven primarily by enhanced evaporative emissions, with alkanes, particularly butanes and pentanes, dominating at 67 % of emissions. Low atmospheric pressure at higher altitudes enhances evaporative emissions, with butanes and pentanes showing emission factors 4.5–21 times higher than at sea level, confirming pressure's dominant role over temperature. In contrast, tailpipe exhaust emissions, such as ethylene and benzene, show no clear altitude dependence. Source profile similarity analysis further elucidates differences among emission sources, revealing that VOC profiles in plateau tunnels closely resemble gasoline vapor from SHED tests, highlighting distinct high-altitude emission characteristics. Source apportionment analysis reveals that evaporative emissions contribute the most to VOCs, accounting for 67 %, followed by tailpipe exhaust (24 %) and non-vehicular sources (9 %). This is in stark contrast to low-altitude regions, where tailpipe emissions dominate. These results suggest that vehicular VOC emissions in high-altitude regions, such as the Tibetan Plateau, may be significantly underestimated. They underscore the critical need to account for unique environmental factors, such as low atmospheric pressure, when evaluating emissions in these areas. Furthermore, this study emphasizes the importance of developing altitude-specific emission inventories and leveraging the region's abundant renewable energy resources to reduce atmospheric impacts.