Exploring E ﬀ ective Chemical Indicators for Petrochemical Emissions with Network Measurements Coupled with Model Simulations

: A large petrochemical complex, dubbed Petro–complex, situated in a rather rural region of Taiwan, was used as a test bed to detect emissions from the Petro–complex to its surroundings. Hourly observations of speciated non–methane hydrocarbons (NMHCs) by the photochemical assessment monitoring stations (PAMSs), as well as the total amounts of NMHCs, SO 2 , and NO x provided by the air quality stations (AQSs), were utilized to ﬁnd useful petro–emission indication methods. The analytical aspect of NMHCs either as a speciated form or as total amounts was demonstrated through ﬁeld comparison to illustrate data quality. Using ethyne to o ﬀ set tra ﬃ c inﬂuence, the ratios of ethene to ethyne (acetylene) (E / A) and propene to ethyne (P / A) were proven to be e ﬀ ective indicators of petro–emissions owing to pronounced emissions of ethene and propene, revealed as tall spikes in PAMS measurements. SO 2 and NO x were also explored as petro–emission indicators mainly for stack sources. By coordinating with three–dimensional modeling, SO 2 from petro–emissions could be distinguished from other prominent sources, such as coal–ﬁred power plants. An attempt was also made to use SO 2 and NO x as indicators of broader petro–emissions with stringent criteria to minimize tra ﬃ c interference and increase speciﬁcity. Similar ﬁndings were observed with the three indicators, that is, volatile organic compounds (VOCs) ratios, SO 2 and NO x , to identify the southwest area of the Petro–region as the most a ﬀ ected area, as represented by Taisi station (F2). The percent a ﬀ ected time of a year at F2 was 10%–14%, owing to the dominant wind ﬁeld of northeast monsoonal (NEM) in the region, as compared with other sites in the east and north of 1–5%. Using VOC ratios as petro–emission indicators is more e ﬀ ective than using other gases, owing to the compositional advantage to minimize tra ﬃ c interference.


Introduction
The petrochemical industry contributes significantly to the global economy. However, air pollutants from petrochemical facilities, such as volatile organic compounds (VOCs), particulate matter, SO 2 , NO x , heavy metals, polycyclic aromatic hydrocarbons (PAHs), and so on, have been frequently reported, some of which have been recognized as environmentally hazardous [1][2][3][4][5]. VOCs of stationary sources are mainly emitted from flaring, feedstock storage and distribution, solvent/reagent usage, as well as other chemical processes [6][7][8]. For a region with a major petrochemical facility, fugitive emissions, flaring, and stack emissions could be the major sources of VOCs dominating the ambient levels [9][10][11]. The AQS network began monitoring hourly mixing ratios (or concentrations) of NOx, NO, SO2, CO, PM10, O3, as well as the total level of NMHCs (called AQS-TNMHC hereafter), from February 2013 to the present. The working theory of AQS-TNMHC is based on a subtractive method. The TNMHC concentration in an air sample aliquot is obtained from subtracting methane concentration from the total hydrocarbons (THCs). No separation of NMHC is needed. Usually, one measurement of TNMHC can be completed within minutes. Thus, minutely data can be averaged to become hourly data. The monitoring instruments at the nine AQS conformed to the specifications of the Taiwan Environmental Protection Administration (Taiwan EPA, see Supporting Material Table S1). Hourly measurements of 54 VOCs at the nine PAMS also began in October 2013. The working theory for PAMS measurements is based on gas chromatography (GC) to separate individual NMHCs of C2 to C12 with GC capillary columns. In our PAMS measurements, the separation needs an hour to complete, and thus the PAMS data are in hourly resolution. Both AQS-TNMHC and PAMS measurements involve the same detection principle, that is, flame ionization detection (FID). As FID basically works as a carbon counter, less complication can be anticipated when inter-compared between the two datasets.
Quality assurance and quality control (QC/QA) of PAMS observations were performed on a permonth basis. With the exclusion of routine maintenance, rare malfunctions, and QA injections of standard gas, data completeness was higher than 85% on average for all stations throughout the year [33].
The data length in this study is from October 2013 to September 2015. Hourly measurements for The AQS network began monitoring hourly mixing ratios (or concentrations) of NO x , NO, SO 2 , CO, PM 10 , O 3 , as well as the total level of NMHCs (called AQS-TNMHC hereafter), from February 2013 to the present. The working theory of AQS-TNMHC is based on a subtractive method. The TNMHC concentration in an air sample aliquot is obtained from subtracting methane concentration from the total hydrocarbons (THCs). No separation of NMHC is needed. Usually, one measurement of TNMHC can be completed within minutes. Thus, minutely data can be averaged to become hourly data. The monitoring instruments at the nine AQS conformed to the specifications of the Taiwan  Environmental Protection Administration (Taiwan EPA, see Supplementary Materials Table S1). Hourly measurements of 54 VOCs at the nine PAMS also began in October 2013. The working theory for PAMS measurements is based on gas chromatography (GC) to separate individual NMHCs of C 2 to C 12 with GC capillary columns. In our PAMS measurements, the separation needs an hour to complete, and thus the PAMS data are in hourly resolution. Both AQS-TNMHC and PAMS measurements involve the same detection principle, that is, flame ionization detection (FID). As FID basically works as a carbon counter, less complication can be anticipated when inter-compared between the two datasets.
Quality assurance and quality control (QC/QA) of PAMS observations were performed on a per-month basis. With the exclusion of routine maintenance, rare malfunctions, and QA injections Atmosphere 2020, 11, 439 4 of 15 of standard gas, data completeness was higher than 85% on average for all stations throughout the year [33].
The data length in this study is from October 2013 to September 2015. Hourly measurements for wind field, SO 2 , NO x , and AQS-TNMHC mixing ratios were obtained from the AQS network, whereas the 54 individual hourly NMHC mixing ratios were obtained from the PAMS network. Summation of the 54 NMHCs, dubbed PAMS-TNMHC, can be compared with AQS-TNMHC for mutual validation. When calculating PAMS-TNMHC, the true mixing ratios of all target compounds in parts per billion by volume (ppbv) are multiplied by their own carbon numbers to become methane-based mixing ratios (ppbC), which are then summed together to become the value of TNMHC. For AQS-TNMHC, the FID response of TNMHC as a single response is calibrated by a methane standard to be the value of AQS-TNMHC. As a result, both types of measurements are methane based and thus can be directly compared.

Wind Characteristics in the Petro-Region
In the Petro-region of this study, the wind field is highly seasonal and greatly affecting the transport routes of air pollutants. In order to systematically define air pollution characteristics from the Petro-complex, three major wind flow patterns were identified according to variations in daily wind direction (WD), that is, northeast monsoonal (NEM) flow, southwest monsoonal (SWM) flow, and local circulation/diffusion (LCD) [37]. These three types of wind fields are very characteristic throughout the island of Taiwan, but vary in strength and duration from north to south. In the central region, where the study domain was located, the three wind patterns occurred 79% of the time for the two year period from October 2013 to September 2015, during which the NEM flows accounted for 274 days or 39% of the time, occurring mostly from late fall to early spring (October to February); whereas the SWM flows accounted for 70 days or 10%, occurring mostly in summer; and the LCD accounted for 209 days or 30% of the time, occurring mostly in warm seasons during the period. For the remaining days, the wind types were too vague (such as typhoon or ambiguous wind fields) to be clearly categorized, and hence were not included in the assessment. The methods of identifying petro-emissions versus non-petro emissions in the following discussion are tested under the three characteristic wind fields.

Model Setup
The Fifth-Generation Penn State/NCAR Mesoscale Model (MM5) and the Taiwan air quality model (TAQM) were used for meteorological and air quality simulations. The horizontal grid sizes for the nested domains are 81 × 81 km 2 (D1), 27 × 27 km 2 (D2), 9 × 9 km 2 (D3), and 3 × 3 km 2 (D4), as shown in Figure S1. There are a total of 15 non-uniform sigma levels in the vertical direction from the surface to approximately 13.5 km above the surface, with the lowest level having a thickness of approximately 20 m. The simulation setup has been described in earlier studies [40,41].

Validation of NMHC Observations
Despite separation by substantial distances in the Petro-region, occasionally, very coherent variations could be observed between PAMS sites under similar source influences and prevailing wind ( Figure 2). We use this phenomenon to help validate the data quality of PAMS in addition to the routine quality assurance (QA) procedure. The high synchronicity in TNMHC between the three PAMS sites separated by 7-12 km from each other suggests the required quality in the PAMS's calibration and measurements. Calibration and maintenance under one single protocol for all PAMS sites played a central role to ensure consistent data quality across all PAMS sites. Note that the synchronicity may not occur or disappear when these sites experienced different emission or meteorological influences. However, when the synchronicity does occur, it becomes very beneficial in terms of inspecting internal consistency and the overall data quality. For the two-year dataset, a number of relatively short periods of synchronicity occurred mostly under NEM when the wind field was sufficiently monotonous. In general, PAMS-TNMHC measurements are highly correlated with AQS-TNMHC measurements, despite two different instruments [42]. We use site F1 as an illustration, which is located within the campus of Petro-complex, and is thus closest to the emissions. The two datasets, hourly mixing ratios of PAMS-TNMHC and AQS-TNMHC at site F1, were highly coherent with an R 2 = 0.79 (see Figure 3a,b). Note that, despite the high synchronicity of the two measurements, the PAMS-TNMHC values were consistently lower than the AQS-NMHC values by at least 20%, as can be observed by the much lower baseline of PAMS than of AQSs (Figure 3a), as well as the smaller than 1:1 ratio in the regression of the correlation plot in Figure 3b. The lower abundance of PAMS-TNMHC than of AQS-TNMHC can be rationalized in that only 54 NMHCs are targeted, whereas AQS-TNHC is responsive to all organics. For instance, oxygenated VOCs (denoted as OVOCs), such as alcohols, ketones, ethers, esters, and acids, produced and released in large amounts in the petrochemical industry, can be accounted for by AQS-TNMHC, but not by PAMS-TNMHC, despite a lower response than hydrocarbons. As a result, when it comes to reporting ambient VOC levels, AQS-TNMHC is a much better proxy than PAMS-TNMHC, because all organic compounds can be detected by FID, despite generally higher sensitivity for NMHC than for OVOCs [43,44]. Even though there is a shortfall of PAMS-TNMHC compared with AQS-TNMHC, the high coherence of the two, as shown in Figure 3a, suggests the 54 target NMHCs measured by PAMS played a dominant role in ambient VOC variabilities in the Petro-region. In general, PAMS-TNMHC measurements are highly correlated with AQS-TNMHC measurements, despite two different instruments [42]. We use site F1 as an illustration, which is located within the campus of Petro-complex, and is thus closest to the emissions. The two datasets, hourly mixing ratios of PAMS-TNMHC and AQS-TNMHC at site F1, were highly coherent with an R 2 = 0.79 (see Figure 3a,b). Note that, despite the high synchronicity of the two measurements, the PAMS-TNMHC values were consistently lower than the AQS-NMHC values by at least 20%, as can be observed by the much lower baseline of PAMS than of AQSs (Figure 3a), as well as the smaller than 1:1 ratio in the regression of the correlation plot in Figure 3b. The lower abundance of PAMS-TNMHC than of AQS-TNMHC can be rationalized in that only 54 NMHCs are targeted, whereas AQS-TNHC is responsive to all organics. For instance, oxygenated VOCs (denoted as OVOCs), such as alcohols, ketones, ethers, esters, and acids, produced and released in large amounts in the petrochemical industry, can be accounted for by AQS-TNMHC, but not by PAMS-TNMHC, despite a lower response than hydrocarbons. As a result, when it comes to reporting ambient VOC levels, AQS-TNMHC is a much better proxy than PAMS-TNMHC, because all organic compounds can be detected by FID, despite generally higher sensitivity for NMHC than for OVOCs [43,44]. Even though there is a shortfall of PAMS-TNMHC compared with AQS-TNMHC, the high coherence of the two, as shown in Figure 3a, suggests the 54 target NMHCs measured by PAMS played a dominant role in ambient VOC variabilities in the Petro-region.

Ratios of P/A and E/A as Petro-Emission Indicators
It was found that ethene and propene were the two most pronounced compounds observed by PAMS, possibly owing to the extremely high volatility and the large quantities produced as feedstock of downstream plastic manufacture. Air masses coming from the Petro-complex are likely to be laden with these two very gaseous species. Although these two species do have a traffic source, the amounts were far less than the typical levels observed by the downwind PAMS sites. Concentrations of ethene and propene at F2 appeared as distinct spikes, as illustrated in Figure 4. Comparing F2 with F4, which are at the opposite directions of the Petro-complex on the northeast axis (see Figure 1), more spikes were encountered at F2 than at F4 owing to the dominant wind field of NEM in the region (upper two panels in Figure 4). These spike were often brief in time, for example, within a few hours, to suggest that the emissions coming off the Petro-complex were very direction sensitive. It appears that

Ratios of P/A and E/A as Petro-Emission Indicators
It was found that ethene and propene were the two most pronounced compounds observed by PAMS, possibly owing to the extremely high volatility and the large quantities produced as feedstock of downstream plastic manufacture. Air masses coming from the Petro-complex are likely to be laden with these two very gaseous species. Although these two species do have a traffic source, the amounts were far less than the typical levels observed by the downwind PAMS sites. Concentrations of ethene and propene at F2 appeared as distinct spikes, as illustrated in Figure 4. Comparing F2 with F4, which are at the opposite directions of the Petro-complex on the northeast axis (see Figure 1), more spikes were encountered at F2 than at F4 owing to the dominant wind field of NEM in the region (upper two panels in Figure 4). These spike were often brief in time, for example, within a few hours, to suggest that the emissions coming off the Petro-complex were very direction sensitive. It appears that when the wind direction was in alignment with F2, ethene and propene as high as several hundred ppbv could be detected. The mixing ratios fell off rapidly and when the wind direction shifted even slightly to miss the monitoring site. As the area surrounding F2 is extremely rural, mixing ratios of ethene and propene were very close to the detection limits of the PAMS instrument. As a result, the sharp contrast between the sharp spikes and the low baseline was exploited as a form of indication of petro-emissions.
Atmosphere 2020, 11, x FOR PEER REVIEW 7 of 15 when the wind direction was in alignment with F2, ethene and propene as high as several hundred ppbv could be detected. The mixing ratios fell off rapidly and when the wind direction shifted even slightly to miss the monitoring site. As the area surrounding F2 is extremely rural, mixing ratios of ethene and propene were very close to the detection limits of the PAMS instrument. As a result, the sharp contrast between the sharp spikes and the low baseline was exploited as a form of indication of petro-emissions. To exclude the traffic influence, a reference compound, ethyne (acetylene), which has no petrosources and is mostly from car exhaust, is used as the common denominator in the ratios of ethene/ethyne (E/A) and propene/ethyne (P/A) to cancel off traffic contribution, as well as variability owing to air mixing. These three species can be obtained simultaneously from the PAMS network To exclude the traffic influence, a reference compound, ethyne (acetylene), which has no petro-sources and is mostly from car exhaust, is used as the common denominator in the ratios of ethene/ethyne (E/A) and propene/ethyne (P/A) to cancel off traffic contribution, as well as variability owing to air mixing. These three species can be obtained simultaneously from the PAMS network with hourly resolution. Thus, E/A and P/A were used as effective indicators of petro-emissions in the region, and the thresholds of E/A > 3 and P/A > 1.5 were set to differentiate air masses affected by the petro-emissions from those that are not affected [33]. The usefulness of the two ratio indicators is illustrated for the three characteristic seasonal wind fields as follows.
Under the NEM flow, which prevailed mostly in fall and winter, wind direction in the region was mainly between 315 and 45 degrees; hence, Site F2, in the southwest of the Petro-complex, would be more affected by the Petro-complex's emissions than F4~F9, which are spread over from the north to the east of the Petro-complex and further inland. Under the SWM flow, which mostly prevailed in summer, wind direction was mostly between 135 and 225 degrees; hence, petro-emissions were more likely to affect site F4, situated at the north of the Petro-complex (see Figure 1 for site F4 location). Under the LCD flow, the affected downwind sites varied considerably owing to the rapidly changing wind direction, and the influences from the Petro-complex would be more towards the inland area than the coastal. Some typical cases are demonstrated in Figure 5. Under the NEM flow (Figure 5a), the simultaneous increase in PAMS-TNMHC mixing ratios as well as P/A and E/A ratios at site F2 indicated that this area was more affected by the petro-emissions during the evening hours. Although a concentration peak also appeared at F9, the lower-than-threshold in P/A and E/A values suggests that F9 was not affected by petro-emissions, but more likely by traffic sources. Under the SWM flow (Figure 5b), the wind field is nearly the opposite of NEM, and high PAMS-TNMHC mixing ratios appeared at night for all sites, but only at noon for site F4; however, the E/A and P/A ratios suggest that only site F4 was affected by petro-emissions at noon. Thus, F2 and F4 can be either the upwind or downwind site depending upon the seasonality and the prevailing wind field.
Atmosphere 2020, 11, x FOR PEER REVIEW 8 of 15 with hourly resolution. Thus, E/A and P/A were used as effective indicators of petro-emissions in the region, and the thresholds of E/A > 3 and P/A > 1.5 were set to differentiate air masses affected by the petro-emissions from those that are not affected [33]. The usefulness of the two ratio indicators is illustrated for the three characteristic seasonal wind fields as follows. Under the NEM flow, which prevailed mostly in fall and winter, wind direction in the region was mainly between 315 and 45 degrees; hence, Site F2, in the southwest of the Petro-complex, would be more affected by the Petro-complex's emissions than F4~F9, which are spread over from the north to the east of the Petro-complex and further inland. Under the SWM flow, which mostly prevailed in summer, wind direction was mostly between 135 and 225 degrees; hence, petro-emissions were more likely to affect site F4, situated at the north of the Petro-complex (see Figure 1 for site F4 location). Under the LCD flow, the affected downwind sites varied considerably owing to the rapidly changing wind direction, and the influences from the Petro-complex would be more towards the inland area than the coastal. Some typical cases are demonstrated in Figure 5. Under the NEM flow (Figure 5a), the simultaneous increase in PAMS-TNMHC mixing ratios as well as P/A and E/A ratios at site F2 indicated that this area was more affected by the petro-emissions during the evening hours. Although a concentration peak also appeared at F9, the lower-than-threshold in P/A and E/A values suggests that F9 was not affected by petro-emissions, but more likely by traffic sources. Under the SWM flow (Figure 5b), the wind field is nearly the opposite of NEM, and high PAMS-TNMHC mixing ratios appeared at night for all sites, but only at noon for site F4; however, the E/A and P/A ratios suggest that only site F4 was affected by petro-emissions at noon. Thus, F2 and F4 can be either the upwind or downwind site depending upon the seasonality and the prevailing wind field.
Under the LCD flow (Figure 5c), the PAMS-TNMHC mixing ratios varied dramatically within a time frame of only a few hours owing to the rapidly changing wind directions. The VOC-based ratios indicated that site F2 and F4 were highly influenced by the petro-emissions, whereas F9 was less affected mostly around noontime.  ratios indicated that site F2 and F4 were highly influenced by the petro-emissions, whereas F9 was less affected mostly around noontime.

SO 2 and NO x with Modeling to Indicate Stack Plumes
In addition to fugitive emissions and flaring, which can be detected by PAMS-derived indicators of E/A and P/A, flue gas released as plumes from stacks of the Petro-complex containing both VOCs as well as SO 2 and NO x can be put to the test for the effectiveness as indicators. The idea is that SO 2 and NO x in the plumes should be much higher in concentrations than outside the plumes. Although stack plumes from the petro-industry contain relatively high levels of SO 2 and NO x , their specificity is generally low owing to the lack of compositional features. For instance, differentiation of plumes from the petro-industry from the plumes coming off the fire power plants can be difficult, as both types of sources contain high levels of the two pollutants. Unless plume detection by these two gases can be supported by other means, the specificity could be in question. As a supporting tool, air-quality modeling was used to enhance the petro-plume indication with SO 2 and NO x . For demonstration, three characteristic wind fields of NEM, SWM, and LCD are exemplified as follows.
Despite the rather monotonous wind characteristics with the NEM flow, the source-receptor relationship can still be ambiguous at times. Figures 6 and 7 illustrate an example of the wind direction affecting plume source identification. A large coal burning fire power plant (5824 mWatts), the second largest worldwide, located 60 km north of the Petro-region (as show in Figure 1), releases SO 2 and NO x from stacks of 250 m height. At a very specific angle of the NEM flow, the smoke plumes from the power plant may merge with the plume tracks from the Petro-complex and become indistinguishable from the petro-plumes, and thus are often perceived as local petro-emissions. However, this misconception can be easily divulged by the combination of on-site measurements and modeling. Thus, hourly data of SO 2 and NO x provided by AQS were examined for the plume influences on the Petro-region. Under the prevailed northeasterly wind, site F2, a downwind site of both the power plant and the Petro-complex, can receive the merged plumes or single plume form the Petro-complex. When the wind slightly shifted in direction, the merged plumes then diverged, and only F2, at the southwest of the Petro-complex, was affected by the Petro-complex's stack emissions, as revealed by the elevated SO 2 level (black box in Figures 6 and 7a). Moreover, site F4 is at the northeast of the Petro-region, and thus could not possibly be influenced by the Petro-complex under the prevailed northeasterly wind. Hence, the stack plumes elevated the SO 2 and NO x concentrations at F4 at around 09:00, as indicated by the yellow box in Figure 6. Site F2, however, can receive the merged plumes. Figure 7b provides a complete pictorial view of the plume tracks from both sources under the particular wind flow. The slightly elevated SO 2 and NO X levels at other sites to the east were mostly contributed by the power plant plumes. The divergence of the two plume tracks became even more obvious in the evening hours as the wind shifted from northeasterly to northerly. F9, the most inland and eastward site from the Petro-complex, was entirely affected by the power plant emissions and not by petro-plumes to show the highest levels of SO 2 and NO X among all sites, as indicated by the red box in Figure 6 and the simulation in Figure 7c. As a result, using SO 2 and NO x as plume indicators to investigate the source-receptor relationship can be obscure at times if not aided by other means such as modeling to provide a more realistic view (See Figure S2 for the comparison of simulations with observations on 10 October 2013).    Figure 8 for the observed concentration peaks (black box in Figure   310 8a) accompanied by simulation results (Figure 8b). Again, measurements and modeling combined Under the SWM flow, clean marine air was brought from the south to the Petro-region. All stations except F4 would show comparable and relatively low levels of pollutants owing to the influence from the clean marine air and greater mixing heights in warm seasons. However, site F4 would occasionally show elevated levels of SO 2 and NO x when the petro-plumes arrives at the downwind site F4, as shown in Figure 8 for the observed concentration peaks (black box in Figure 8a) accompanied by simulation results (Figure 8b). Again, measurements and modeling combined are able to help investigate the source-receptor relationship with clarity (See Figure S3 for the comparison of simulations with observations on 19 June 2014).
Note that, although Figure 6, Figure 7, and Figure 8 only show brief periods of data, they are intended to illustrate typical features of the source-receptor relationship under the very repetitive and seasonal wind patterns. With the tools of modeling and densely deployed monitoring facilities, these features can be adequately characterized. Under the SWM flow, clean marine air was brought from the south to the Petro-region. All stations except F4 would show comparable and relatively low levels of pollutants owing to the influence from the clean marine air and greater mixing heights in warm seasons. However, site F4 would occasionally show elevated levels of SO2 and NOx when the petro-plumes arrives at the downwind site F4, as shown in Figure 8 for the observed concentration peaks (black box in Figure  8a) accompanied by simulation results (Figure 8b). Again, measurements and modeling combined are able to help investigate the source-receptor relationship with clarity (See Figure S3 for the comparison of simulations with observations on 19 June 2014).
Note that, although Figures 6,7, and 8 only show brief periods of data, they are intended to illustrate typical features of the source-receptor relationship under the very repetitive and seasonal wind patterns. With the tools of modeling and densely deployed monitoring facilities, these features can be adequately characterized.

Test of SO 2 and NO X as Petro-Emission Indicators
The PAMS derived ratios of E/A > 3 and P/A > 1.5 have proven to be effective petro-emission indicators. However, PAMS is usually not a typical establishment for air quality assessment worldwide.
For general air quality indication, SO 2 and NO x are more readily available as standardized methods adopted by environmental protection agencies. When assisted with model simulations, SO 2 and NO x can be used as plume tracers, as illustrated earlier. However, when traffic source comes into play in more populated areas, the conditions of using SO 2 and NO 2 as petro-indicators have to be restricted by more stringent boundaries to be indicative enough of petro-emissions.
The attempt of using SO 2 and NO x as petro-emission indicators was bound by the following three criteria. First, the Petro-complex has to be upwind of the given monitoring sites in the Petro-region under a given prevailing wind direction. Second, the SO 2 and NO x data used for analysis are required to be significantly greater than the baseline values, that is, twice the standard deviations or higher. Note that the baseline values for the two gases were determined by the lowest 25% of the dataset. Third, SO 2 and NO x should be free from prominent sources other than the Petro-complex, such as fire power plants, steel mills, and cement plants. For instance, if SO 2 and NO x at both F2 and F4 exceeded the baseline criteria, then the downwind site (F2) was considered as polluted by other sources than the petro-emissions, and thus the data were excluded.
The percent time affected by Petro-complex emissions during 2014 and 2015 for the areas represented by F2, F4, and F9, representing south, north, and east, were calculated with the two years of hourly monitoring data based on the three different indicators of VOC ratios, SO 2 , and NO x . The percent time affected at F2 was 14.0%, 10.7%, and 9.7% for VOC ratios, SO 2 , and NO x , respectively; 1.8%, 3.4%, and 4.0% at F4, respectively; and 1.4%, 3.1%, and 5.1% at F9, respectively. Because the wind characteristics in the region are very seasonal and repetitive year after year, we believe the findings based on two years of hourly measurements are sufficiently representative of the region. It is evident that all three indicators were able to identify F2 as the most affected site by the Petro-complex, as it was in alignment with the Petro-complex under NEM, which is the dominant wind field in colder seasons with lower mixing heights and rather monotonous wind directions. In contrast, the other sites to the east (as represented by F9) or north (as represented by F4) were remarkably less affected by the Petro-complex. Although all three indicators resulted in similar findings of identifying F2 as the most affected site, the use of SO 2 and NO x was subjected to many restrictions when to minimize traffic interference. Furthermore, the discrepancy between the three indicators for a given site could result from the incomplete removal of traffic interference despite strict reconditions set for SO 2 and NO x . In this regard, VOC-based indicators can be more useful than other gases, as they can be extremely sensitive and better exempted from traffic interference owing to the compositional advantage.

Conclusions
A large Petro-complex in a rather rural area was exploited as a test ground for investigating petro-emission indicators. The extensively deployed PAMS and AQS sites in the region providing hourly data of speciated NMHCs, as well as SO 2 and NO x , were utilized for indication of petro-emissions. Data quality of PAMS played a pivotal role in the data utilization, and thus was validated by comparison among PAMS sites and between AQS-TNMHC and PAMS-TNMHC to show internal consistency. The VOC ratios of E/A and P/A were found very useful as petro-emission indicators owing to the prominent emissions of ethene and propene from refinery. SO 2 can be used as an indicator for stack emissions if assisted with model simulations to decouple contribution from other prominent SO 2 sources such as distant fire power plants. When there is no PAMS available, SO 2 and NO x can be used as emission indicators of the Petro-complex, provided that they are bound by stringent criteria to minimize traffic interference. Similar findings were observed with the three indicators of VOC ratios, SO 2 , and NO x to identify F2 as the most affected site owing to the dominant wind field of NEM in the region. The VOC-based indicators are more effective in tracing petro-emissions than other gases such as SO 2 or NO x owing to the compositional advantage with the VOC data to minimize traffic interference.