1. Introduction
The ramp acts as a connection between highways and local roads. Significant differences in operating conditions exist between the roadway facilities that ramps connect, in terms of speed, elevation, and roadway configuration. Hence, vehicles accelerate, decelerate, and navigate up and down on these ramps. The Motor Vehicle Emission Simulator (MOVES), the regulatory emission model invented by the U.S. Environmental Protection Agency [
1], defines ramp activity as vehicle operations that occur between the unrestricted road (local road) and the restricted road (highway/freeway). MOVES assumes that ramp activity comprises about 8% of nationwide freeway operating hours for all source types [
1], based upon the ramp fraction value originally developed by Systems Applications International, Inc. for the MOBILE model [
2]. The large fraction of ramp mileage and the high-power operations that occur on ramps lead to a significant contribution of ramps to the transportation emissions inventory. Ramps are potential emissions hotspots if significant localized emissions occur at these locations. In Environmental Protection Agency (EPA) project-level Particulate Matter (PM) conformity hotspot analysis guidance [
3], ramps are presented as an introductory example for MOVES and AERMOD (American Meteorological Society/Environmental Protection Agency Regulatory Model) analysis.
Generally, light-duty vehicles on-ramp emit about two to three times more than vehicles on highways (excluding ramps) and local cycles with an average speed greater than 30 mph (48.28 km/h) [
4]. Despite the potential impact of ramps on emissions and air quality, relevant literature is quite limited. Liu et al. [
5] compared energy and emissions of the passengers using aircraft, bus, and automobile in the U.S. To collect the operations data for the cases of Detroit on-ramp, Thornton et al. [
6] used a laser rangefinder, presented speed–distance plots of the on-ramp activity, and estimated the emissions of ramp metering schemes via the MEASURE emission rate model. This research was a first step in understanding ramp operations and emissions, although the extension of these results to the general description of on-ramp activities is limited.
Recently, the influence of ramp metering on emission reduction has been considered in research. Bae et al. [
7] adopted ramp metering to evaluate the changes in traffic flow and CO
2 emission variations, which suggested that ramp metering makes a great contribution to the reduction of CO
2 emissions. Du et al. [
8,
9] found that the isolated and integrated ramp metering approaches can improve mobility and reduce travel delay. However, relevant studies are limited and emissions vary with vehicle type. The emissions of light-duty vehicles have been studied a lot [
10,
11,
12,
13], but research about the emission of heavy-duty vehicles is relatively less.
The research presented in this paper assesses the observed on-road operations and modeled emissions of transit bus activity on-ramps. The road grade has been considered in the model in view of its great influence [
14,
15]. Analyses include on-ramps and off-ramps, in both diagonal and loop ramp configurations [
16]. The analyses employ the Global Positioning System (GPS) data, including second-by-second speed and position data, collected from a fleet of the Detroit transit system. Researchers extracted ramp-associated operations from the data using known ramp locations and Geographic Information System (GIS) tools, and then linked ramp operations data with applicable MOVES emission rates for the analyses.
This paper presents the data preparation in
Section 2 and analytical methods in
Section 3, respectively.
Section 4 provides analytical results for on-road operations and emissions including presentation and discussion of: speed–acceleration plots, ramp emission rates based on the average speed, hotspot analysis, and road grade effects. Research results are useful for understanding transit bus operations and emission characteristics on various types of ramps, the potential impacts of ramp grade on emissions, and the ramp hotspot analysis. Finally,
Section 5 concludes the research.
2. Data Preparation
To evaluate the emissions and fuel consumption of transit fleets, transit operations data were collected by GPS from 12 buses operating in the Detroit Metropolitan Area. The data included 381 bus-days between October 28th, 2014 and October 27th, 2015. At the time of data collection, the Detroit transit fleet included more than 500 transit buses (seating capacity of 33–42). About half of the fleet were compressed natural gas (CNG) buses, and the other half were diesel buses.
2.1. Data Preprocessing
Before emission analysis, collected GPS data need to go through the Kalman filter and spline process [
17] to ensure their validity, accuracy, and continuity. GPS data quality varies with the GPS signal condition, characterized by the number of satellites and positional dilution of precision (PDOP) values. The modified discrete Kalman filter algorithm effectively controls GPS random errors. The Kalman gain matrix is generated using GPS quality criteria, including the quantity of satellites and the PDOP values. If the number of satellites is smaller than four and the PDOP is greater than eight, the quality of the speed is determined to be poor; therefore, the error variance in the matrix is high. The modified discrete Kalman filter algorithm was applied to the Detroit transit operation data to minimize the random errors from GPS loggers.
Another factor affecting the quality of GPS data is missing data. Missing segments exist in almost all of the GPS trip files due to obstruction and signal interference. If the missing data segments are of short duration, missing data can be reasonably interpolated from the speed profile. The research team used a cubic spline algorithm [
18] to infill speed data when missing data segments were three seconds or less. Six good speeds points surrounding the missing data were required. Spline data consist of the three closest seconds before the start of the missing segment, all within 10 seconds before the missing segment, and the three closest seconds after the end of the missing segment, all within 10 seconds after the missing segment. If the good speeds could not be found within the 20 seconds surrounding the candidate segment, this segment was identified as missing and the splining process was abandoned. Missing segments longer than three seconds terminated the current drive trace (spitting the trace into two separate traces at that breakpoint).
After applying the Kalman filter and spline process, cycles less than 60 seconds in duration were discarded. The remaining 4717 vehicle-hours of transit bus operation were used to extract ramp operation data.
2.2. ArcGIS-Based Ramp Data Extraction
In this research, ArcGIS [
19] analysis routines identify vehicle operations on-ramps. The shapefiles of the most recent Detroit road network were downloaded from Topologically Integrated Geographic Encoding and Referencing (TIGER) [
4,
20], and imported with the bus GPS data into ArcGIS. The roadway shapefiles indicate ramp locations, since the ramp acts as a transition between arterials and highways. To identify the road class shift, a 22 m buffer (3.66 m width/lane × 6 lanes per direction) was created for each restricted-access highway centerline in each direction. To identify the physical rap boundaries, satellite images were overlaid on the ArcGIS base maps. Transit bus GPS data and roadway buffers identified the portions of the vehicle traces that were running on restricted highways, on the ramps, and on the local roads. The following boundary identification rules applied:
Start of Freeway on-ramp
- ○
For diagonal ramp interchanges
- ○
For loop ramp interchanges
End of Freeway on-ramp
- ○
For diagonal ramp interchanges
- ○
For loop ramp interchanges
Start of Freeway off-ramp
- ○
For diagonal ramp interchanges
- ○
For loop ramp interchanges
End of Freeway off-ramp
- ○
For diagonal ramp interchanges
- ○
For loop ramp interchanges
Located at the free-flow connection with highway or local road
- ○
For parallel-type ramp terminals, the ramp boundary was taken as the midpoint of the lane changing area (see
Figure 1a))
- ○
For taper-type ramp terminals, the end of solid traffic line was taken as the ramp boundary (see
Figure 1b))
Ramp free-flow connection with local road
- ○
If no extension lanes applied, the edge of the access was taken as the boundary (see
Figure 1c))
- ○
If extension lane applied, the midpoint of the lane changing area was taken as the boundary (see
Figure 1d))
A lot of manual labor was involved in establishing the ramp boundaries to ensure that GPS data were properly assigned to each ramp. Through this process, we also flagged points with ramp types. A total of 161 drive traces on-ramps were collected and analyzed (
Table 1).
Figure 2 shows the ramps layout covered by the transit GPS data and the transit bus in Metro Detroit.
3. Methodology
3.1. Ramp Cycle Determination or Ramp Effect Area
After entering the highway from an on-ramp, the vehicle speed is normally lower than the design speed on the highway and the speed of traffic flow, forcing drivers to accelerate continuously and merge into the mainline traffic. In the same way, drivers tend to slow down when they leave the highway off-ramps. As such, highway operating speeds immediately after exiting on-ramps or before entering an off-ramp are likely to be affected by the presence of the ramps and the traffic volumes entering/exiting the freeway. For the analysis of emissions, including a portion of the highway as part of the ramp operations cycle is necessary, since high-speed and high-power operations may be associated with ramp effects. The ramp-associated driving traces were determined as follows:
On-ramp cycle: speed trace within ramp plus the extension after the vehicle enters highway, i.e., post-ramp effect, in which vehicle accelerates continuously to the highway design speed;
Off-ramp cycle: speed trace within off-ramp plus the extension of trace before the vehicles enter ramp (pre-ramp effect), in which vehicle decelerates beforehand to the ramp design speed;
Figure 3 presents two examples of speed traces assigned to the local road, ramp, and highway.
The duration of the extended cycle is dependent on factors such as an individual’s driving behavior, ramp geometry, and traffic conditions. Based on the speed and acceleration characteristics, we created a method to determine the duration of the ramp cycle extension. A bracket was used to comprise five-consecutive second points on highways at the cycle right before or after the ramp boundary, and the speed/acceleration values of these five were judged. If the five points satisfied the criteria in
Table 2, these points could be treated as part of the ramp cycle extension, and the bracket was moved forward by a one-second step at the post-ramp-associated trace (or backward for pre-ramp associated trace) and the same judgment was made. This process was carried out in a loop until points in bracket broke the criteria, where the midpoint (third point) will be regarded as the boundary of the ramp cycle.
Figure 4 below presents the diagram of the bracket used to implement the algorithm and the movement along the speed trace. According to the ArcGIS process, the bold dashed lines are determined as the physical boundary of the ramp. Additionally, the light dashed line represents the boundary of the ramp operating trace. The algorithm is applicable to data obtained in this study. For other types of data sources, such as light-duty vehicles, the standard needs to be modified.
After leaving a diagonal on-ramp, it took an average of 15.6 s of post-ramp extension (standard deviation of 13.3 s) for bus drivers to reach the design speed of the highway. By contrast, it took longer for a post-ramp extension (an average value of 24.5 s and a standard deviation of 14.1 s) to approach the design speed after leaving a loop on-ramp. The reason is drivers tried to keep a lower speed when driving through a loop ramp, and a longer post-acceleration time was required. The pre-ramp extension of a diagonal off-ramp cycle was just 2.9 s on average, and the extension of a loop off-ramp was 6.5 s.
3.2. STP Calculation and Operating Mode Bin Generation
In this research, emissions were determined through emission inventory from MOVES 2014. For each second, scaled tractive power (STP) and an operating mode bin as defined in MOVES needed to be generated as a bridge between bus activity and emission inventory [
1]. The second-by-second STP formulation is expressed as:
where
A: the tire rolling resistance coefficient (kilowatt second)/(meter)
B: the mechanical rotating friction coefficient (kilowatt second2)/(meter2)
C: the aerodynamic drag coefficient (kilowatt second3)/(meter3)
M: the source type in metric tons (metric tons)
: scaling factor for heavy-duty vehicles
v: the speed of the vehicle (m/s)
a: the vehicle acceleration (m/s2)
sin θ: the road grade (percent)
g: the acceleration rate due to gravity (9.8 m/s2)
Detroit transit bus models correspond to transit bus (source use type ID = 42) in the MOVES model, so A, B, C, M, and
were determined as 1.03968, 0, 0.00358702, 16.556, and 17.1, respectively from the MOVES database
sourceusetypephysics. The second-by-second operating mode bin was generated by STP, speed, and acceleration value. The criteria of each operating mode bin were summarized in [
20], as well as in
Table 3. Most of the results in this research were obtained by assuming a flat gradient (sin θ = 0), considering the terrain in Detroit is relatively flat. These results can be regarded as the benchmark and the ramp grade effects will be specifically investigated in the last section as a sensitivity analysis.
3.3. Emission Rates Assignment
Emissions results were based on MOVES emission rates with real-world speed trace and GPS speed data. We obtained second-by-second emissions for each ramp cycle as the fundamental data set for the statistical analysis. As the sample size (the number of ramp cycles) becomes larger, which in this case was 161, it would take much time in input preparation and running MOVES. In addition, a lot of replicated calculations would be needed in the sensitivity analysis of the ramp gradient effect, and running MOVES is neither convenient nor economical for multi-scenarios with only one parameter (grade) changed.
Instead of running ramp cycle in county-level MOVES, we ran MOVES 2014 in project-level to get the emission rates for bus applied with national default age distribution [
21], and obtain the emission rates in g/h for each operating mode bin of the fleet. Here is the scenario we set to get the emission rates from MOVES:
Month and Year: July 2019
Date and Time: weekday, 5:00–6:00 p.m.
Region: Detroit
Meteorology (default value):
Temperature: 85.6 °F (27.8 °C)
Humidity: 50.3%
Fuel: default value from MOVES 2014
Inspection/Maintenance (I/M) Strategy: default 2019 I/M strategy from MOVES
Source Type: transit bus (source type ID = 42)
Age Distribution: 2019 national default values
23 Links: each link was applied with 100% fraction of one operating mode bin, and the time was scaled to 1-hour operation using the link length and link average speed
The emission rate was assigned to each second from the matrix table of bus emission rate, based on the operating mode bin of each second. For example, the CO emission rate of a bus per second with operating mode bin ID as 15 (9 kW/metric ton ≤ STP < 12 kW/metric ton, 1 mph ≤ speed < 25 mph, or 1.61 km/h ≤ speed < 40.23 km/h) was assigned as the CO emission factor 447.698 g/h for operating mode bin 15 from look-up table, i.e., emission rate as 184.84 (g/h) ÷ 3600 (second/hour) = 0.124 g/s. This process was more efficient than directly running MOVES since it saved time in duplicated calculations. Moreover, it has been validated that the aggregated emission results from second-by-second emission rates are the same as running the MOVES model. We accomplished the emission assignment through Java and its connection with MySQL [
22].
Five criteria pollutant emissions (HC, CO, NOX, PM10, and PM2.5), greenhouse gases (GHGs), and fuel consumption rate were estimated. GHGs is quantified as CO2e, or carbon dioxide equivalent, a measure to compare the emissions from three greenhouse gases based on their global warming potential from vehicle sources using CO2 as the reference over 100 years. For example, the global warming potential for methane over 100 years is 21. This means that emissions of one million metric tons of methane is equivalent to emissions of 21 million metric tons of carbon dioxide. Based on MOVES calculation, PM10 is linearly related to PM2.5, and fuel and CO2e are linearly related to CO2. It is known that the majority of HC emissions come from a cold start and evaporations instead of from the running mode. The result analysis in this research mainly focused on CO, NOX, PM2.5, and CO2.
From the second-by-second data, we aggregated operations and emissions by ramp cycle, by ramp functions, and by roadway segments. Furthermore, the speed trace and aligned emission rates can help to understand the emission factor of the ramp cycle in different segments, which is useful for the hotspot analysis.
4. Results and Discussion
4.1. Speed–Acceleration Plots
Figure 5 shows the speed–acceleration results from on-ramp cycles and off-ramp cycles, where points from within ramps and extension areas are flagged with different colors.
For on-ramp cycles, 87% of points were classified as in acceleration mode, 35.2% of on-ramp cycle points were contributed from extension areas that defined where vehicles continue to accelerate right after entering highways to reach the desired highway speed. Points in extension areas (red points) were in high-speed level (≥ 64.37 km/h or 40 mph) and acceleration mode, in which vehicles tend to generate high emissions, and this will be reflected in the aggregated emission results and hotspot analysis.
For cycles in off-ramp, 78.8% of points were in deceleration mode (see
Figure 5). Only 6.6% of off-ramp cycle points resulted from pre-ramp extension (yellow points) defined as where drivers decelerate right before entering ramps to adjust from the highway speed to the desired ramp speeds. Emissions from off-ramp cycles are expected to be relatively low compared with cycles from other road types or ramps.
4.2. Emission Rates Based on the Average Speed
Emission rates in g/km/vehicle by ramp cycle were calculated by summing up the second-by-second emissions and dividing by distance. In
Figure 6, the scatterplots indicate the average speed and emission rates (g/km) on the ramp cycle for CO, PM
2.5, NO
X, and CO
2, flagged with different colors for on-ramp and off-ramp, and with Locally Estimated Scatterplot Smoothing (LOESS) applied [
23], or local polynomial regression, to represent the overall trends of speed–emission rates. The emission rate of the same bus fleet under MOVES highway and local cycles was also calculated and shown below as black and blue curves. From the comparison, at the same speed level higher than 64.37 km/h or 40 mph, since a lot of acceleration modes involved in on-ramps, emission rates for the on-ramp cycle were approximately twice higher than emission rates from MOVES highway cycles. Moreover, on-ramp emissions were diverse under the same speed level, indicating unknown factors that need to be examined further in detailed on-ramp emission modeling, in which factor might be individual behaviors, ramp configuration, traffic condition, etc. At the low-speed level (< 64.37 km/h or 40 mph), as the average speed decreased, on-ramp emission rates decreased to approaching the highway emission rates, and this was due to the lower fraction of aggressive accelerations involved. Off-ramp emission rates were almost half of the highway emissions since a lot of deceleration modes were involved.
4.3. Hotspot Analysis
Since transit buses produce a lot of emissions when operating on the on-ramp cycle, second-by-second speed (km/h), second-by-second emissions (g/s), and emission rates (g/km) by sub-cycle segment (i.e., within ramp, post-extension, and highway) of on-ramp cycle were presented to explore the significance of a potential hotspot.
Figure 7 provides an example of hotspot analysis by illustrating the second-by-second speed and MOVES-based CO and PM
2.5 emission rates for a transit bus driving through a diagonal on-ramp and a loop on-ramp. It should be noted that acceleration is the main operation mode that contributes to on-road emissions, which mainly happened on the on-ramps. There were not much emissions coming from off-ramps, since the operation modes during off-ramps are mainly braking and deceleration. Therefore, we only present the on-ramp results. Based on the MOVES estimation, a lot of emissions were produced within a diagonal on-ramp. By contrast, for a loop on-ramp, lots of emissions were generated in the post-ramp extension area. This is because drivers generally keep a comparatively stable and low speed as they go through the loop ramp, and then are compelled to accelerate to the highway design speed after entering the highway, resulting in highway emissions in the extension area. This indicated the importance of involving the extension area into the ramp operation hotspot analysis. The results shown in
Figure 7 are clear to help us identify the hotspots of emissions.
Within each on-ramp cycle, CO and PM
2.5 emission rates were measured for each sub-cycle segment, i.e., within the ramp, by post-ramp extension, and connected highway, based on real-world operations. The emission rates of all 96 on-ramp cycles were aggregated into the average values as shown in
Figure 8 below, and the standard deviation was flagged on top of each bar chart.
In general, for diagonal on-ramps, emission rates of transit buses reached a peak within the on-ramp area. While for loop on-ramps, emission rates of transit buses in the post-ramp extension area reached the peak. Emission rates of transit buses ramp sub-cycle within both ramp area and post-ramp extension (in the highway) were higher than those of vehicles in highway operations. This indicated that for the emission analysis of on-ramps, besides on-ramp segments, emission rate in a certain extension area on connected highways could also be high, particularly for the loop on-ramp. Therefore, we should include post-extension areas as well as on-ramps while conducting the hotspot analysis.
4.4. Sensitivity Analysis of the Ramp Grade
All previous results were based on the assumption of flat terrain with road grade set at zero, it needs to be noticed that road gradient is indeed a factor that could significantly affect emission results and is especially likely to contribute to emissions on-ramp, since the connected local road and highway tend to be at different elevations. However, it is not easy to obtain robust altitudinal data or road slope data. Therefore, this research conducted a sensitivity analysis. By doing this, road grades in the ramps were set from −6% to +6%, with 1% interval, which covered most cases, since ramp grade is required to be within ±5%, and preferred within ±4% [
16]. Emission rates in g/km/vehicle of each ramp cycle in each grade scenario were calculated based on the same method. Grades of post-ramp sub-cycle were still set at zero since the connected highway and local road were assumed to be relatively flat.
Emissions from “flat” ramp cycle (grade = 0%) were regarded as reference values. Ramp grade effect could be reflected by comparing emission rate changes of all scenarios with the “flat” scenario. The emission change is expressed as:
where
Figure 9 below presents the average emission changes in the on-ramp cycle and off-ramp cycle under different ramp grade values. A significant contribution of emissions could be due to the effect of road grades: for example, when the road grade changed from 0% to +4%, the overall CO, NO
X, and CO
2 emission rate of transit driving through upward on-ramp could be increased by 35–45%, and the PM
2.5 rate could be increased by as much as 50%. The estimated emissions grow linearly with the grade, which might be due to the capped operating mode bin in MOVES that restricts further emissions estimation on extremely high engine load. Real-world data are required to better describe the impact of road grade.
5. Conclusions
This research presented operations and emission characteristics of transit buses driving through ramps. Second-by-second speed and position data were collected in 2015 from the Detroit transit system, the principal rapid transit system in Detroit. Ramp-associated operations were extracted from daily trips and assigned with emission rates from MOVES 2014. Some key findings are summarized below:
In ramp operation and emission analysis, it is significant to consider not only bus behaviors within the ramp, but also pre-ramp decelerations before buses entering the off-ramps from the highway, and post-ramp accelerations after vehicle enter the highway from the on-ramps, where vehicle operations are different from operations of the mainstream highway traffic.
In this research, based on the transit bus, on loop on-ramp operations, CO and PM2.5 rate of the bus during post-ramp operations are even greater than operations within ramps, with the presumption of zero road grade, which indicates it is necessary to treat the post-ramp area separately from the highway links in the hotspot analysis.
At a speed level of higher than 40 mph (64.37 km/h), emission rates of transit buses under on-ramp operations are about twice higher than emission rates from the MOVES highway cycles. At the low-speed level (< 40 mph or 64.37 km/h), as the average speed decreased, on-ramp emission rates reduced to approach the highway emission rates due to the lower fraction of aggressive accelerations involved. Off-ramp emission rates are almost half of the highway emissions. On-ramp emissions are varying, indicating unknown factors that need further research in detailed on-ramp emission modeling.
The sensitivity analysis of the ramp grade indicated that, for project-level analysis, ramp grade is an unavoidable factor. When ramp grade changed from 0% to +4%, the overall CO, NOX, and CO2 emission rates of transit buses driving through upward on-ramp could be increased by 35–45%, and the PM2.5 rate could be increased by as much as 50%.
This study focused on heavy-duty buses. Although the bus emissions are significantly higher than those of light-duty vehicles, the results show similar trends with MOVES-based study on the light-duty vehicle emissions [
4]. The on-ramp emission rates at on-ramps are 1.5–2 times larger than freeway emissions, depending on the average speed, which is consistent with the real-world based study conducted by [
10].
This research was useful for presenting driving characteristics and emissions of transit buses on different types of ramps. It is a further step for ramp classification, investigating the influence of ramp gradient, and the development of the ramp cycle for emission modeling. Some suggestions were provided in terms of implementing hotspot analysis for the on-ramps. While only transit data were applied in this study, the method can be extended to the heavy-duty trucks or the light-duty vehicles.