Simulation and measurement of air quality in the traffic congestion area

The traffic congestion in the Hsuehshan tunnel and at the Toucheng interchange has led to traffic-related air pollution with increasing concern. To ensure the authenticity of our simulation, the concentration of the last 150m in Hsuehshan tunnel was simulated using the computational fluid dynamics fluid model. The air quality at the Toucheng interchange along a 2 km length highway was simulated using the California Line Source Dispersion Model. The differences in air quality between rush hours and normal traffic conditions were also investigated. An unmanned aerial vehicle (UAV) with installed PM2.5 sensors was developed to obtain the three-dimensional distribution of pollutants. On different roads, during the weekend, the concentrations of pollutants such as SOx, CO, NO, and PM2.5 were observed to be in the range of 0.003–0.008, 7.5–15, 1.5–2.5 ppm, and 40–80 μgm , respectively. On weekdays, the vehicle speed and the natural wind were 60 kmh 1 and 2.0 m s , respectively. On weekdays, the SOx, CO, NO, and PM2.5 concentrations were found to be in the range of 0.002–0.003, 3–9, 0.7–1.8 ppm, and 35–50 μgm , respectively. The UAV was used to verify that the PM2.5 concentrations of vertical changes at heights of 9.0, 7.0, 5.0, and 3.0m were 45–48, 30–35, 25–30, and 50–52 μgm , respectively. In addition, the predicted PM2.5 concentrations were 40–45, 25–30, 45–48, and 45–50 μgm 3 on weekdays. These results provide a reference model for environmental impact assessments of long tunnels and traffic jam-prone areas. These models and data are useful for transportation planners in the context of creating traffic management plans.


Introduction
Indoor and outdoor air quality is a hot research topic and a major determinant of human health. More attention has been given to indoor air quality in recent years, as it has more significant effects on respiratory disease and cardiovascular health than outdoor air pollution [1]. According to previous studies, almost 90% of our time is spent in enclosed areas (e.g., homes, schools, offices, transport, and meeting places) [2][3][4][5][6]. Many people are exposed to high concentrations of traffic contaminants when they drive in heavy traffic and spend time at places near roads that have a high amount of traffic, especially if the location is downwind of the road [7,8]. Air pollutants (e.g., volatile organic compounds, nitrogen oxides, and particulate matter) are complex and dynamic, and vehicle emissions have become the dominant source of air pollutants [9,10]. Campagnolo et al. [11] demonstrated that vehicle emission reductions have a strong impact on the effective control of human in-cabin exposure and improve the air quality in traffic environments. The increasing severity and duration of traffic congestion have the potential to increase pollutant emissions and degrade air quality, particularly near large roadways. Therefore, a study on air quality caused by vehicle emissions is necessary to help understand the effect of air pollution as a health risk to drivers.
There are several methods for modeling the dispersion of different pollutants emitted from a roadway reported in the literature [12,13]. The Gauss dispersion California Line Sources Dispersion Model (CALINE4) treats emissions as a continuous line source, either without any adjustment or using a simple enhancement [14]. Dhyani and Sharma [15] pointed out that the CALINE4 requires relatively lower levels of expertise and comparatively fewer input data than other vehicular dispersion models. The computational fluid dynamics fluid model (CFD) is popularly used in environmental engineering [16,17]. For example, Wang et al. [18] recommended the optimal deflect angle for tunnel ventilation using CFD. Qin et al. [19] recommended the numerical simulation of airborne HCHO pollution in vehicle cabins using CFD. Through CFD simulation, it can be an accurate prediction of air pollution flow properties and species transport [20][21][22]. The estimation of the airflow patterns and flows in tunnels and other complex structures was mathematically modeled by CFD [23]. To accurately predict the species transport and flow properties, governing flow equations were solved to simulate atmospheric turbulence [24]. In addition, unmanned aerial vehicles (UAVs) are commonly used in the military, agricultural surveillance, and transportation fields [25]. The UAV has flexibility and mobility to be used at metropolises or suburban areas. UAVs are low-cost technology and have application potential in various fields, including collecting detailed military, traffic volumes, and meteorological data [26,27]. A UAV with sensors was a new application method to measure the threedimensional spatial distribution of air pollutants [28].
The rapid development of tourism in Yilan County, Taiwan, has led to traffic-related air pollution with increasing concern. The 12.9 km Hsuehshan Tunnel is located between Taipei and Yilan. The Hsuehshan Tunnel provides convenient travel but simultaneously causes air pollution due to poor dispersion conditions compared to an ordinary road [29]. The motivation of this study was to use the CALINE4 and CFD to investigate the air quality in the traffic congestion area and provide a reference method that can be applied to estimate the health risks to drivers due to exposure as well as obtain useful data for stakeholders. The first objective was to simulate the concentrations of the last 150 m in the Hsuehshan tunnel using CFD. The second objective was to predict the air quality at the Toucheng interchange over a 2 km length of the highway using CALINE4 and obtain the threedimensional distribution of pollutants using a UAV for a validated air quality forecast. Many countries have fewer than three monitoring stations/million inhabitants, and their location is often restricted. If these models can predict air quality, they would help reduce the enormous workload involved in on-site measurement. The results of this study could provide stakeholders with information on the air quality inside long tunnels or near tunnel exits/inlets and help them draw up useful traffic management plans.

Study region and data collection
In this study, the air quality was assessed along a highway at different points, i.e., the last 150 m in the Hsuehshan tunnel and over a 2 km length of the highway at the Toucheng interchange (as shown Fig. 1). Near the Toucheng interchange, there are many recreational landmarks (e.g., hot springs and dolphin watching) and tourist attractions with more than 1 hundred hotels, and many fresh seafood restaurants. Nevertheless, the level of bus ridership in Yilan County is lower than that in other metropoles in Taiwan. Therefore, traffic congestion in the Hsuehshan tunnel and the Toucheng interchange has long been a characteristically obvious reference.
The traffic volume data were sourced from the databank of the Ministry of Transportation and Communications (MOTC) and Taiwan Area National Expressway Engineering Bureau (TANEEB). Basic information on the Hsuehshan tunnel and Toucheng interchange was obtained from the databank of TANEEB. Due to safe driving considerations, vehicle information at 5-min intervals can be collected using cameras taking images of vehicle discs/hubcaps in the Hsuehshan tunnel. The information contains the total number of heavy-duty vehicles (e.g., bus) and light-duty vehicles (e.g., sedan and van). According to the data in the databank of MOTC and our previous research, the number of vehicles in 2019 was two times higher than that in 2009 [29]. The emission rate of pollutants was sourced from the Taiwan Emission Data System (TEDS) of Environmental Protection Administration (EPA). The weather information was obtained from operational numerical forecasting models. The meteorological data were summarized from the Toucheng station of the Central Weather Bureau. The coordinates and height of the monitoring site are 121.81 and 24.85 and 6.0 m, respectively.

Model design and establishment
A modified CALINE4 was used to predict local traffic emissions within 2 km for four traffic-related pollutants (CO, NO, SO 2 , and PM 2.5 ) based on the Gaussian diffusion model. The data that the simulation software needs to enter are road type pollutants, emission rate, receptors, and meteorology (e.g., wind direction, wind speed, atmospheric stability, and mixing height) [15], as shown in Table 1. The grid area was established within 200 × 200 m along with the line source. The blue points are receptor points, and the red star is monitoring site, as shown in Fig. 1. All of the receptor points are extended from the center. In addition, each interval distance is 200 m to compare the pollutant concentration at each receptor point. The simulated results with CALINE4 were also compared with monitored data during holidays and weekdays in March, 2019.
In this paper, the CO, NO, SO 2 , and PM 2.5 pollutant concentrations were simulated at the tunnel exit during the traffic jam period and normal period, and the relationship between air quality and traffic volume was established. ANSYS Fluent 7.0 software has well-validated physical and chemical modeling capabilities, accurate results across the widest range of CFD, and multiphysics applications. ANSYS Fluent 7.0 is the numerical simulation software for Windows Version. The step is the geometric description, mesh, selecting master equation and boundary conditions, and numerical solution [19]. In addition, ANSYS Fluent software has very good performance optimization capabilities and can predict flow, turbulence, heat transfer, and reactions for industrial applications. Therefore, this software was used in this study for the air quality simulations. A numerical analysis was used to assess the air quality of tunnels in different conditions to further improve the associated ventilation equipment, pollution prevention, and control strategies. In this study, we used the k-ε turbulence model as [20,22,23]:  where K is the thermal diffusivity (m 2 s − 1 ), Xi is Cartesian coordinates, V is the kinematic viscosity (m 2 s − 1 ) of air, β is the volumetric coefficient of expansion (K − 1 ), ε is the dissipation rate, Ui is instantaneous velocity (m s − 1 ) component i, and θ is the instantaneous temperature difference.
In addition, an in-tunnel simulation was performed using CFD, including (i) a pollutant concentration simulation during holiday and nonholiday periods; (ii) nitrogen oxides, carbon monoxide, and particulate matter concentration prediction near the exit of Toucheng tunnel; and (iii) a flow field condition estimation in the tunnel. The geometry was created using the real tunnel specifications and vehicle sizes, which were applied as the differential and integral terms in the fluid dynamics' fundamental equations, including time/space variables and physical variables. CFD calculation processes include the analysis of balance and fluid dynamics equations, which are the momentum balance, mass balance, and energy balance equations. The governing equations are listed as follows [20,22,23]: Momentum Equation Conservation Equation Conservation of Energy

Conservation of Contaminant
where C is the instantaneous concentration of passive contaminant (kg m − 3 ), D is the molecular diffusion coefficient (m 2 s − 1 ), g is the gravitational acceleration (m s − 2 ), P is the instantaneous static pressure difference (N m − 2 ), ρ is the density (kg m − 3 ) of air, H is the volume heat source generation rate (kW m − 3 ), S is the flow rate of contaminant generation source (m 3 s − 1 ), and k is the kinetic energy (kg m 2 s − 2 ). The temporal and spatial variables and the physical variables had to be discretized to replace these integral or differential terms with discrete algebraic forms. Discrete spatial variables corresponding to the solution domain were divided into a series of lattices, called the unit body or control body (mesh, cell, control volume). The grid corresponding to the lattice boundary is called the grid. In addition, the intersection of the grid is called the grid point. Similar to algebraic methods, differential methods are also used to generate grids, as used by Thompson et al. [30]. For differential methods, discrete physical variables are often defined at grid points. The differential operation on a grid point can be approximated as the algebraic relationship between the physical point on the grid point, the adjacent grid points, and the grid point coordinates. The numerical method, in this case, is called the finite difference method. In this study, the mesh type was a tetrahedron, and the total grid number was approximately 2.9 million. Bhautmage and Gokhale [24] pointed out that modeling the shape of the automobile is complex. Therefore, simplified vehicle body shapes of average dimensions were built. As shown in Fig. 2, the vehicle objects were categorized as a bus, a sedan, a van, and a small truck, which have variations in size, architecture, and pollution discharge rate. The aspect ratio of the vehicle models was 4.5 × 2 × 2 m for the van, 4.5 × 1.5 × 1.5 m for the sedan, 8 × 2.5 × 3 m for the bus, and 5.5 × 1.5 × 2.5 m for the truck. The exhaust pipe diameter for a bus and a van was set to 6 cm. In addition, the exhaust diameter of the sedan and small truck was 4 cm. The exhaust pipe was positioned behind the sedan's right bottom corner (which is 0.2 m above the road surface), van (which is 0.4 m high), and truck (which is 0.3 m high) in the realtime data for model validation. 0.5 m high at the bottom in the middle of the bus. An exhaust pipe in the vehicle was defined as an emission source. The effect of the shape on velocity estimation was neglected due to the minor variation in the physical properties between the small truck and van for further analysis.

Measurement and data analysis
In the CFD simulation and calculation, the monitored concentration data were sourced from the EPA. In the CALINE4, the blue points are receptor points, and the red star is the monitoring site (5 m away from the border of the road), as shown in Fig. 1b. This study measured the CO, NO, SO 2, and PM 2.5 concentrations near the Toucheng intersection. The sampling site address is No. 6, Sec. 1, Qingyun Rd., Toucheng Township, Yilan County 261, Taiwan. The data were collected every hour from 18th March to 28th April in 2019. PM 2.5 sampling was performed using the atmospheric particulate matter analysis high volume sampling method. Atmospheric particulates were collected on quartz filter paper after sampling. A UAV with ambient pollutant sensors is emerging as a new means to produce three-dimensional observations of ambient air pollutants, as shown in Fig. 3. The objective of this study was to validate air quality forecasting by UAVs with PM 2.5 sensors. PM 2.5 was measured only since the bias and variation of the other sensors are too large. The sensors are two types, one type is 1.45 kg, and the other type is 100 g. The sensors were made in two companies. All of the flying missions are all along the highway. The UAV can remain at the same point for 20 s and the time interval of measuring is a second. In addition, to ensure accurate measurement values, the sensor was used to collect sampling site data with EPA pollutant station data (in Yilan Fushing Junior High School) to calibrate the sensor. To collect threedimensional air pollutant concentration data, the experiment was designed around a 2 × 2 km area. In addition, the points have the same coordinates as the simulation, and data were collected at heights of 3, 5, 7, and 9 m to monitor the vertical profile of air pollutants. Normal mean square error (NMSE) is the standard deviation of the residuals (prediction errors) and a measure of how to spread out these residuals [14]. It is traditionally used in climatology, forecasting, and regression analyses to verify experimental results. The fractional bias (FB) is a measure of the correlation between averages of predicted and observed values. The FB is based on the averages of the predicted and observed concentrations (macrostatistics) rather than a datum-by-datum comparison (microstatistics) [14]. In this study, the statistical performance of the CALINE4 was validated using NMSE and FB.

Vehicle volume and vehicles emission rate
The road between Shihting Interchange and Toucheng Interchange prohibits the carriage of dangerous goods in both directions, including long, wide, high, and overweight vehicles. Small cars and buses can drive on the road between Shihting Interchange and Toucheng Interchange. Small cars include small passenger cars, small trucks, and small special passenger vehicles. Vehicle emission rates change at different speeds. Information about emission rates was created from TEDS. The particulate pollutant ratio sourced from TEDS for Toucheng is shown in Fig. 4a. Linear pollutants accounted for approximately 48.4%. Accounting for the largest number, the diesel truck represents 27.0%. The contributions of restrained dust and diesel small truck emissions were 12.3 and 4.9%, respectively. The rest of the pollution sources from the area were commercial (4.2%) and agricultural burning (26.3%). The data were sourced from TEDS; the emission rates of buses, vans, trucks, and sedans for all kinds of pollutants are shown in Fig.  3. It also shows that vehicle emissions have a strong impact on the dispersion of air pollutants, such as CO and NOx, near Toucheng.
The Hsuehshan tunnel is one of the longest tunnels in the world. The road connects Taipei through New Taipei to Yilan County and shortens the journey from 2 h to just half an hour. Most frequently, this trip is made by car or by public transportation (though there are only a few services) through Hsuehshan Tunnel and Highway No. 5. As a result, the traffic through the tunnel is severely congested during the traffic jam period, increasing the travel time from 30 min up to 2 h. The alternative method is a winding shoreline highway, or an equally dangerous mountain pass, both of which take much longer and are more dangerous because of their cliff-side locations. The data from the databank and monitoring system of MOTC regarding the traffic volume are shown in Fig. 4b and c. The average traffic volume was approximately 33,000 passenger car units (PCUs) d − 1 on the weekend and 1800-2000 PCU h − 1 during rush hour. It also indicated that the maximum was more than 35,000 PCU d − 1 on the weekend. In addition, the tunnel standard of pollutants is shown in Table 2.
In this study, the emission factor with different vehicle speeds was used in Fluent software from TEDS. The CO, NO, SO 2 , and PM 2.5 pollutant emission rates changed when the vehicle speed changed. For example, the CO emission rate decreased with increasing traffic speed, especially from the exhaust of vans and buses, as shown in Fig. 4; with the increase of mean vehicle speed, vehicle emissions of a comprehensive factor gradually decreased. In this study, the vehicle speeds ranged from 20 to 100 km h − 1 . In addition, the SO 2 emission rate was lower than before because of lower gasoline sulfur content. The vehicle with the lowest emission rate was the sedan [31,32]. The SO 2 emission rate increased with increasing mean vehicle speed when the vehicle speed was below 50 km h − 1 . In contrast, the emission rate decreased with increasing mean vehicle speed when the vehicle speed was above 60 km h − 1 since diesel engines are used in van. In addition, diesel engines are also used in buses and trucks in similar   Both heavy-duty and light-duty types of vehicles emit large amounts of PM 2.5 , CO, and NO and affect the tunnel's air quality. The pollutant emission factor of heavyduty vehicles is larger than that of light-duty vehicles. However, the vehicle number of heavy-duty vehicles is smaller than that of light-duty vehicles. The total emission rate of light-duty vehicles is smaller than that of heavy-duty vehicles. According to our previous study and the databank of MOTC, traffic volume during nonrush hours on weekends can be approximately 60% higher than that during non-rush hours on weekdays [29]. Therefore, traffic volume on the weekend could be regarded as rush hour, and on the weekday, it could be regarded as non-rush hours. On weekdays, the pollutant concentrations are mainly distributed in the tunnel outlet. In addition, the SO 2 , CO, NO, and PM 2.5 concentrations ranged from 0.001 to 0.002, 2 to 6, 0.36 to 1.08 ppm, and 20 to 55 μg m − 3 , respectively, at ambient temperatures of 23 to 25°C, as shown in Fig. 6. From the simulation results, it was found that the concentration of various pollutants in the tunnel increases with a higher number of vehicles. The concentration of various pollutants becomes higher in the tunnel when the number of vehicles entering the tunnel increases. Thus, it is recommended to control the flow of vehicles within 4000 h − 1 to control the CO, NO, and PM 2.5 concentrations to meet the air quality standards.

Concentrations simulated at the Toucheng intersection
In this section, the comparison with modeled pollutants concentration in different height level is discussed. In this work, the air quality at the Toucheng interchange along a 2 km length of the highway was simulated by CALINE4. To analyze the PM 2.5 pollutants along the highway at different heights, both height and concentration data were used, as shown in Fig. 7. To evaluate the vertical concentration of particle matter along the highway, the highest concentration of PM 2.5 was seen at 9 m at a level of 74 μg m − 3 , located at the intersection of tunnel exit, Toucheng intersection exit, and province freeway No. 2. The levels at 7 m located at the intersection of tunnel exit, Toucheng intersection exit, and province freeway No. 2 were 50-58 μg m − 3 .
The results show that the pollutant concentrations decreased with increasing distance from the highway and tunnel. The concentration of PM 2.5 at 5 m was almost the same as that at 7 m. However, the pollutants were stuck in the middle since the concentration spread from the center. In addition, the highest concentration was 54-58 μg m − 3 . At 3 m, the situation was different because the concentration increased again. The highest concentration, 58 μg m − 3 , was found at Provincial Road No. 9. The highest PM 2.5 concentrations are located at 9 m since the height of the highway is near 9 m, and the PM 2.5 concentrations are highest near the highway. The PM 2.5 located at 3 m was slightly higher than that at 7-5 m since the data located at 3 m were affected by the vehicle's emissions on Provincial Road No. 9 with a height of 0 m. The dispersion concentration became lower as the distance increased from where the traffic jam occurred. At 3 m, the concentration becomes higher than at 7 m because human activity and large trucks raise particulate matter levels. In addition, the concentration of CO at 9 m was approximately 7 ppm, and the concentrations at 7.0, 5.0, 3.0, and 1.5 m were approximately 5.8, 5.8, 8.5, and 10.0 ppm, respectively, as shown in Fig. 8.

UAV monitoring
It is difficult to measure the vertical pollutant concentration near a high level without UAVs. Therefore, it is good to verify the prediction data along the vertical direction with UAV data. To obtain the characteristics of the three-dimensional space distribution, a UAV with a sensor hanging system was used to monitor pollutants' levels. The vertical air quality distribution was monitored with UAV carrying sensors. Using the UAV, the air quality monitoring instruments are taken up and down vertically to obtain spatial distribution data. The rotor of UVA affects the measured values with the UAV carrying sensors. As shown in Fig. 9a, the relative correlation of the hanging length of the sensor goes to 0.85 0.73, 0.32, and 0.09 at 1.5, 1.0, 0.5, and 0.3 m, respectively. In addition, the wind speed increases as the height increases. Therefore, positioning the sensor too high is dangerous. Therefore, the best hanging height was approximately 1.5 m. As shown in Fig. 9b, the UAV PM 2.5 monitoring was compared with the EPA station in Yilan Fushing Junior High School (the height is approximately 14 m), and the data for each hour were very close between them. However, UAV monitoring was also affected by weather conditions. Figure 10 presents the concentration of PM 2.5 with a 3D distribution. Generally, the concentration of PM 2.5 decreased as the height increased, and the results were similar to those of another study [33]. Due to the minimal ground emission source, no difference in all flights was observed in the horizontal distribution. On weekday mornings, the PM 2.5 concentrations was approximately 40-45 μg m − 3 at 1.5 m. The human activity caused the emissions of pollutants to accumulate [34,35]. The vehicle volume on the highway during the week days was lower than that on weekends, and the concentration decreased as the height increased. There is an incremental temperature in the air, while there is a rise in the solar energy and radiation levels. Until the early afternoon, there was a gradual increase in the concentration of the

Model sensitivity analysis
The inputs for the CALINE4 package are roadway geometry, meteorological parameters, such as temperature and humidity, emission rate, traffic volume (in vehicles per hour), wind speed, wind direction and its standard deviation. In CALI NE4, the wind direction standard deviation (S) was calculated according to the Eq. (6) and it ranged from 0 to 20 0 .
where X i is the wind direction, X is the average of wind direction, and n is the monitored number of wind direction. The concentration rapidly dropped as the wind speed increased, especially at 0.5-1 m s − 1 , and the CO concentration changed from 10 to 50 ppm. The wind direction standard deviation is also important when the S is more than 20 0 . The CO concentration tended to be moderate, and the standard deviation of the wind direction was 5 0 at a CO concentration of approximately 120 ppb. The standard deviation of wind direction was more than 20 0 when the concentration of CO remained at 50 ppb. At a high mixing layer height (1000 m), the concentration of CO was 1.4-2.5 ppb. At a lower mixing layer height (500 m), the maximum CO concentration was approximately 2.5-4.0 ppb. The contaminants were uniformly mixed by turbulence in the mixed layer since the concentration decreased as the mixing height increased as the contaminants were transported and diluted. In contrast, the low mixing layer height promotes the formation of atmospheric heavy pollution processes. Atmospheric stability has obvious effect on the simulation results. Therefore, the correlation coefficient of the meteorological conditions affects the simulated CO concentration, the mixed layer height, the wind speed, and wind direction, as shown in Fig. 11. The traffic volume data are from TANEEB. They offer real-time traffic information from a camera. The daily traffic volume in hours was collected by Electronic Toll Collection. From the simulation results, it was found that the concentration of various pollutants in the tunnel increased with the increase in the number of vehicles. The concentration of various pollutants was higher in the tunnel when the number of vehicles entering the tunnel was greater. The correlation coefficient between volume and concentration in the tunnel increased with the increasing number of vehicles entering the tunnel, as shown in Fig. 12a. The correlation coefficients between CO, PM 2.5 , NO, and SO 2 concentrations and the number of vehicles are 0.99, 0.94, 0.89, and 0.89, respectively.

Validation against the measurement data
The comparison of the modeled concentrations with the monitored concentrations is a process termed verification [18,20]. Model verification investigates the discrepancies between modeled and measured concentrations, which can arise as a result of the presence of inaccuracies or uncertainties in model input data and/or modeling and monitoring data assumptions [36,37]. The following are examples of potential causes of such discrepancies: estimates of background pollutant concentrations, meteorological data uncertainties, traffic data uncertainties, model input parameters, and overall limitations of the dispersion model [38]. In the CFD simulation and calculation, it was necessary to check the correlation coefficient. The monitored concentration data were sourced from the EPA. The sampling positions were located in one safety passing bay and the inlet/outlet of the tunnel. In the CALINE4 model, we compared the measured and predicted concentrations near the Toucheng intersection. In addition, it is difficult to measure the vertical pollutant concentration near a high level without UAVs. Therefore, a UAV with a sensor hanging system was used to monitor the levels of pollutants, and it was good to verify the prediction data along the vertical direction. Figure 12b indicates that the correlation coefficient of CFD was 0.9374. This is very close to the simulation data; a reason for this might be that indoor concentrations from the inside tunnel are more controllable.
In this study, we compared the measured and predicted concentrations, for example, the concentrations of NO, CO, and PM 2.5 that were predicted near the Toucheng intersection. On a weekday, the highest NO concentrations were measured at 8 am, as shown in Fig. 13.   [15]. This could be because vehicle heat fluxes increase as there is an increment in traffic volumes, providing higher vertical dispersion and dilution. PM 2.5 concentrations at the measurement locations during weekdays and weekends are shown in Fig. 13. The highest concentration of PM 2.5 is in the morning. It has two peaks in Weekday. The predicted data exceed the measurement at the Toucheng intersection.
In addition to the above comparative analysis, this model was tested using NMSE and FB to understand the statistical performance. The analysis results of the performance of the model data are shown in Tables 3 and 4. The simulation results demonstrate the satisfactory performance of the CALI NE4 and show that the model could be used for exposure assessments in future environmental impact assessments. Hence, according to statistical analysis of simulation model performance for air quality prediction, the values (e.g., FB, and NMSE) are well within the acceptable limits. It could be proposed that the CALINE4 model for the prediction of vehicular emissions in traffic congestion areas is valid.

Conclusions
This study validated the methodology using NMSE, FB, and correlation coefficient to understand the statistical performance. The CALINE4 model and CFD are suitable to estimate traffic-related pollutant concentrations, and a flexible UAV with a sensor hanging system can be chosen to monitor the levels of pollutants. The results show that the correlation coefficient and the other   parameters are in a reasonable range, despite the traffic in this area being complicated. In the Hsuehshan tunnel, the vehicle speed is 40 km h − 1 when the fan is opened on the weekend, and SO 2 , CO, NO, and PM 2. The results provide information on vehicular exhaust emissions in a long tunnel and along a 2 km length of the highway, including information on air quality together with information on the three-dimensional spatial distribution of PM 2.5 concentrations. These data are useful for transportation planners working towards a traffic management plan. Meanwhile, this study provides a reference method to be used for environmental impact assessments for long tunnels and areas prone to traffic jams. For future research, CFD could simulate volatile organic compounds in indoor environments to evaluate the relationships between volatile organic compound use and health. Compared with the traditional measurement pollutant concentration method, UAVs with sensor hanging systems have many advantages in measuring vertical pollutants, thus they address the limitations of future studies.