3 Luc J.J. Wismans * 4 5 University of Twente 6
Faculty of Engineering Technology 7
Centre for Transport Studies. 8
P.O. Box 217, 7500 AE, Enschede, the Netherlands 9
T. +31.570.666840, F. +31.570.666888, E. l.j.j.wismans@ctw.utwente.nl 10
11
Robert M.M. van den Brink 12
13
Goudappel Coffeng 14
Research and Development department 15
P.O. Box 161, 7400 AD Deventer, The Netherlands 16 T. +31.570.666839, F. +31.570.666888, E. rvdbrink@goudappel.nl 17 18 Luuk J.N. Brederode 19 20 Goudappel Coffeng 21
Department Transport innovation & modelling 22
P.O. Box 161, 7400 AD Deventer, The Netherlands 23 T. +31.570.666801, F. +31.570.666888, E. lbrederode@goudappel.nl 24 25 Kobus J. Zantema 26 27 Goudappel Coffeng 28
Department Transport innovation & modelling 29
P.O. Box 161, 7400 AD Deventer, The Netherlands 30
T. +31.570.666868, F. +31.570.666888, E. kzantema@goudappel.nl 31
32
Eric. C. van Berkum 33
34
University of Twente 35
Faculty of Engineering Technology 36
Centre for Transport Studies. 37
P.O. Box 217, 7500 AE, Enschede, the Netherlands 38
T. +31.53.4894886, F. +31.53.4894040, E. e.c.vanberkum@ctw.utwente.nl 39
40
Paper prepared for the 92th Annual Meeting of the Transportation Research Board, 2013. 41 * corresponding author 42 Word count: 43 44 Text 5,503 45 Tables 1 x 250 = 250 46 Figures 8 x 250 = 2.000 47 Total 7,753 48
ABSTRACT 1
For the estimation of air quality and noise nuisance often the output of static traffic assignment models is 2
used. However, static traffic model have several limitations related to the realism of dealing with 3
saturated traffic conditions, which are important for the estimation of emissions. Static traffic assignment 4
models predict congestion on the wrong locations and can not deal with traffic flow phenomena like 5
spillback. As a result also air quality or noise bottlenecks can be predicted on wrong locations, which may 6
lead to erroneous policy decisions concerning new infrastructural measures or implementing mitigating 7
measures on wrong locations. Dynamic traffic assignment models can take these traffic flow phenomena 8
into account and are therefore more suitable to estimate emissions and can in the long term replace the use 9
of static traffic assignment models for larger networks. To provide insights in the differences between 10
estimates based on static and dynamic traffic assignment models in realistic cases, a highway case and 11
urban case are used applying the ARTEMIS model (air emissions) and AR-INTERIM-CM model (noise 12
emissions). The comparison for both cases shows that large differences are found in the estimates for 13
noise (i.e. > 3 dB(A)) and also for the local estimates of air emissions (e.g. up to 45% for NOx emissions).
14
When an infrastructural measure is taken the comparison shows that not only the location and extent of 15
effect, but also the direction of effect can be different (i.e. increase versus decrease and vice versa). 16
17 18
INTRODUCTION 1
2
Sustainability has become a key aspect in the field of traffic and transport and considering the 3
externalities of traffic is therefore of increasing interest when policy decisions are made. Most countries 4
have regulations regarding air quality and noise, but also in earlier planning stages these aspects are 5
important in the appraisal of alternatives (e.g. using cost-benefit analysis). These aspects, but also energy 6
use and CO2 emissions, are becoming more and more important. To assess the effects of measures
7
transport models can be used. Traffic assignment is the step in transport modeling in which trips are 8
assigned to the network by confronting demand with supply resulting in route choice, loads and traffic 9
conditions. Different types of assignment models can be used for the assessment of measures, and can be 10
classified into static and dynamic models. Static traffic assignment (STA) models describe the interaction 11
between travel demand and infrastructure supply assuming that the traffic demand and supply are time-12
independent, hence constant during the considered time period (stationary). The basic output of these 13
models are link loads (amount of traffic using individual roads) and average link travel times or speeds. 14
STA models are generally used at the strategic level in order to carry out long-term studies into effects of 15
(mainly mobility-) measures. Dynamic traffic assignment (DTA) models are typically flow propagation 16
models over time that calculate the resulting traffic conditions, taking changes in supply and demand over 17
time into account. 18
In most cases the output of STA models is used to assess the effects on emissions [1]. The output 19
used is related to traffic volume, traffic composition (i.e. passenger cars and trucks) and the level of 20
service (i.e. average traffic conditions). The level of service is mostly based on the volume/capacity (V/C) 21
ratio, direct or indirect via the used travel time functions within a STA or using a post-processing 22
procedure (e.g. estimation of modal activity). However, the limitations of STA particular for over-23
saturated traffic conditions and the importance of taking traffic dynamics (i.e. realistic flow characteristics 24
as well as temporal information) into account when assessing emissions are widely recognized. Many 25
traffic flow phenomena cannot be taken into account which DTA typically can which are important for 26
the estimates of emissions [2]. In general, when using STA, only the average traffic condition per road 27
can be determined, queues appear in the bottleneck instead of upstream the bottleneck, and spillback is 28
typically neglected. Spillback is the phenomenon that a queue on a downstream link affects the possible 29
output volume of the upstream link or links connected to it. Different studies have shown that there is a 30
proven relation between traffic dynamics and externalities. High speeds, significant speed differences 31
between vehicles, and speed variation (accelerating, braking) have for instance a negative effect on air 32
and noise emissions [3,4,5,6]. Using the outcome of DTA, show high potential to improve the emission 33
estimates significantly. Although microscopic models can capture most traffic dynamics, there are several 34
issues in using these models related to the realism of vehicle operations on this microscopic level and the 35
needed computation times. Therefore, especially macroscopic or mesoscopic models may in the long term 36
replace the use of STA models for larger networks. These models can capture many traffic flow 37
phenomena which are important for the emission estimates and can be used on larger networks than 38
microscopic models [7]. 39
Although it has been acknowledged that DTA models provide more realistic flow characteristics 40
and the estimates of emissions improves when DTA models are used, there is limited knowledge on the 41
possible errors made in the estimation of emissions for realistic cases when STA models are used 42
especially when the effects of measures are assessed. The contribution of this paper is therefore to provide 43
these insights by comparing STA and DTA in estimating air or noise emissions for two realistic cases (i.e. 44
a highway network as well as an urban network). The estimates based on these models are also compared 45
with using measurements (i.e. volumes and speeds) as an input for the emission calculations. To be able 46
to perform this comparison both model types are connected with a suitable emission model for air and a 47
suitable emission model for noise which can be used in conjunction with DTA as well as STA assuring 48
consistency in emission modeling. The second section of this paper provides background on emission 49
modeling and current practice. The importance of capturing traffic dynamics for the emission estimates is 1
described in the third section. The fourth section describes the used framework in this research to 2
compare STA and DTA for the estimation of emissions. The comparison using a highway and an urban 3
case is described in the fifth section. This paper ends with the conclusions and discussion. 4 5 6 EMISSION MODELING 7 8 Air emissions 9
To determine air emissions of substance E on road links (a) the next general formula is used: 10
,
a a
E =αK (1)
where α is the emission factor and K amount of vehicle kilometers. 11
This method is generally used in conjunction with traffic models. The factors that are determinative are 12
averaged or used to differentiate the emission factors and vehicle kilometers. There are many different 13
emission models that can be connected to traffic models to estimate emissions. The types and main 14
characteristics that can be distinguished are presented in Table 1 (adapted from [7]). Aggregated emission 15
functions use single emission factors representing a particular vehicle type and a general driving type. 16
Vehicle operation is therefore only taken into account at a very rudimentary level. More sophisticated 17
aggregated emission functions are based on traffic conditions in which cycle average emission rates are 18
correlated with various driving cycle parameters and are referenced to specific traffic situations. These 19
traffic situations relate to certain conditions (e.g. level of service). Average speed models are based upon 20
the principle that average emission factors vary according to the average speed during a trip. In general a 21
continuous average-speed emission function is fitted for several vehicles over a range of driving cycles. In 22
principle, the input is the trip-based average speed, although in practice it is also common to be used for 23
local speed at discrete locations. Within regression models, each driving cycle used is characterized by a 24
(large) number of descriptive parameters (e.g. average speed, relative positive acceleration and number of 25
stops per km) and their derivatives. A regression model is fitted to the average emission values over the 26
various driving cycles, resulting in the determination of the best predictors for emissions. In modal 27
models, emission factors are allocated to the specific modes of vehicle operation encountered during a 28
trip. Different types of modal models exist, and the terminology used can be confusing. A possible simple 29
type defines vehicle operation in terms of modes like ‘idle’, ‘acceleration’, ‘deceleration’ and ‘cruise’. 30
Instantaneous speed based models relate fuel consumption and/or emissions to vehicle speed and 31
acceleration during a driving cycle, typically at one-second intervals. Instantaneous power based models 32
use a description of the engine power requirement in which the most complex models represent physical 33
and chemical phenomena that generate emissions. Methods part of these types are load based using 34
engine maps or methods using a surrogate of engine power by using e.g. the product of speed and 35
acceleration instead of acceleration, and relate fuel consumption and/or emissions to vehicle speed and 36
this product. 37
In current practice often STA models are used to assess traffic and transport policies on 38
emissions. Until recently many countries used an average speed based emission model (e.g. COPERT in 39
Europe and MOBILE in the United States). The past few years the modal approach (e.g. MOVES) and 40
traffic situation based emission models (e.g. ARTEMIS) are emerging [1,8]. The calculated V/C-ratios’ 41
are mainly used to estimate the average traffic conditions. These ratios’ can for instance be used to 42
determine the fractions of traffic situations or operational modes of vehicles. Other often used approaches 43
are calculating link speeds based on the travel time functions possibly in combination with post-44
processing methods to improve speed estimates or to derive other congestion indicators (e.g. delay or 45
queue length). 46
47 48
1
Table 1. Overview model types 2
Model type Function type
Input level Road
Characteristics Vehicle characteristics Driving characteristics Traffic characteristics Example Emissions Aggregated emission function
Discrete Trip or road section
Road type Vehicle class Flow NAEI, CAR,
EMFAC, VFEM
Traffic situation
Discrete Road
section
Road type Vehicle class Flow, Level of
service
HBEFA, ARTEMIS
Average speed Continuous Trip - Vehicle class Flow,
(Adjusted) Average speed
COPERT, MOBILE, TEE,
QUARTET
Regression Continuous Trip Road design Vehicle type Driving cycle - VERSIT+
Driving mode/modal
Discrete Road
section
- Vehicle class Distribution of driving modes - UROPOL, MEASURE, MOVES Instantaneous/ speed based Discrete/ continuous
Continuous - Vehicle class Driving-cycle
(speed/accelera tion)
- DGV, MODEM
Instantaneous/ power based
Continuous Continuous - Vehicle type,
engine load Driving-cycle, used gear - CMEM, EMIT, PHEM, VeTESS Noise Instantaneous/ speed based Continuous Road section Surface type, wetness, gradients.
Vehicle class Acceleration/ deceleration, (used gear)
Flow, speed Harmonoise, CNOSSOS, NORDIC, MITHRA, RMVII, NMPB, ASJ Traffic situation /speed based Continuous Road section Surface type, gradients
Vehicle class Traffic situation Flow, average speed AR-INTERIM-CM TNM 3 Emissions of noise 4
To determine sound emissions, generally functional relationships are used based on statistical analysis per 5
vehicle type and possibly in 1/3-octave frequency bands. In most source models, also called emission 6
laws, rolling noise and aerodynamic noise are combined, but can also be separated functions. Often, the 7
function used is fitted for reference conditions, and for other situations deviating correction factors (γ) are 8
used. These correction factors in most models include influence of road surface, weather and driving 9
conditions (acceleration/deceleration). The general formula for the sound power level (Lw) of a single 10
vehicle is: 11
(
, ref)
corrections,Lw= +α βf v v +γ (2)
where f is either a logarithmic function of the vehicle speed in case of rolling and aerodynamic noise, and 12
a linear function of vehicle speed for propulsion noise. The parameter α is the noise production of a 13
vehicle at the reference speed at a specific distance from the road centre. If the sound power level is 14
calculated per 1/3-octave frequency band, the parameters ,α β and γ are also per 1/3-octave frequency 15
band. In general two types of emission functions can be distinguished, instantaneous or traffic situation 16
based. The instantaneous based assumes smooth and homogeneous flow and possibly uses a correction 17
factor for acceleration. Traffic situation based distinguishes traffic conditions for which different 18
functions are used. In the latter the influence of accelerations is incorporated in the emission functions. 19
The types and main characteristics that can be distinguished are presented in Table 1 (adapted from [7]). 20
Also for noise the current practice is to use the outcome of STA models to estimate emissions. In 1
contrast to the air emissions most countries do no use an estimate of travel speed based on the STA, but 2
use the posted speed limit or operating speed. The new European CNOSSOS model [9] also recommends 3
to use the posted speed limits when estimating noise emissions using STA. Although this is a 4
conservative assumption for highways, for urban road networks this can result in underestimating noise 5
emissions. The correction for driving conditions is often connected with facilities (e.g. intersections, toll 6
plaza’s), which means that for noise the predicted traffic conditions based on the STA model is rarely 7
used. 8
9 10
IMPORTANCE OF TRAFFIC DYNAMICS 11
12
Traffic dynamics and especially saturated traffic conditions are, next to traffic volume and traffic 13
composition, important explanatory variables for the air and noise emissions [7]. In current practice this 14
has also resulted in taking traffic dynamics to a certain level into account when STA are used. However, 15
STA models can only determine the average traffic conditions and STA models have several limitation in 16
capturing traffic flow phenomena and therefore possibly predicts incorrect traffic conditions. Figure 1 17
presents the relation between traffic condition and externalities for different road types based on the 18
ARTEMIS emission model [10]. Four different traffic conditions are distinguished and 10 road types. The 19
importance of traffic dynamics is obvious, based on the differences in emission factors dependent on 20
traffic condition. Emissions increase when traffic conditions get saturated and can be higher at high 21
speeds. Similar relations can be found for other substances like particulate matter (PM10) and carbon
22
dioxide (CO2). For trucks, not presented in this paper, the emissions are especially higher for stop & go
23
traffic conditions. 24
25
Emission factors dependent on traffic condition and road type
Passenger cars, NOx, 2008, Dutch fleet
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 80 100 120 50 60 70 80 30 50 70
Highway Rural road Urban road
Road type E m is s io n f a c to r N Ox ( g r/ k m ) free flow heavy saturated stop & go 26
Figure 1. NOx emission factors dependent on traffic condition 27
28
The relevance of capturing realistic traffic dynamics has been subject of research for many years now. In 29
[11] for example it was found that an error of 5 km/h in the estimate of speed could result in a 42% 30
difference in CO emission predictions. Smit et al. [12] found that considering speed distribution instead of 31
average single mean speed can improve emission predictions since this would be closer to reality and can 32
lead to up to 24% higher emissions of CO, HC, NOx, PM10 and CO2. In [13] it is shown that the ability to
33
take modal activity into account quadruples the estimated emission benefits of signal coordination. In the 34
past years there is an increase in research in which the output of DTA models is used to estimate the air 35
emissions. Often this research is based on using microscopic simulation models [14,15,16,17,18]. Less 1
research has been done in connecting macroscopic or mesoscopic DTA models with emission models. In 2
[19] for example, a framework has been developed for dynamic estimation of traffic emissions using 3
REALITY, a road activity emission calculation model, together with LADTA, an analytical dynamic 4
traffic assignment model. A traffic situation emission model VERSIT+macro, in conjunction with the 5
macroscopic DTA model INDY was used in [12]. The VERSIT+macro emission model consists of 6
composite emission factors for discrete traffic situations. The mesoscopic DTA model Dynasmart-P and 7
the average speed based emission model EMFAC was used in [20] to analyze trip-based versus link-based 8
traffic data for emissions estimation. 9
10
Noise emission AR-INTERIM-CM
Passenger cars (1 vehicle, reference-distance 1 m)
38.0 40.0 42.0 44.0 46.0 48.0 50.0 52.0 54.0 0 20 40 60 80 100 120 140 Speed (km/h) S o u n d p o w e r le v e l (d B (A )) Free flow Saturated 11
Figure 2. Noise emission dependent on traffic condition 12
13
The relation between speed and noise is presented in Figure 2 and based on the AR-INTERIM-CM model 14
[21]. The figure illustrates that emissions increase in saturated conditions, because of an increase of 15
vehicles accelerating and because at low speeds the propulsion noise becomes the dominant source for 16
noise. For trucks a similar relation exist, however in this case propulsion noise stays dominant at higher 17
speeds and as a result the emissions for trucks in saturated traffic conditions can be higher than in free 18
flow conditions on highways. 19
In [22] and [23] it is concluded that there are several benefits in considering temporal information 20
on dynamics of traffic in order to improve road noise emissions estimates, because of identifying local 21
peaks and variations. The relevance of traffic dynamics for noise emission estimates has been subject of 22
research by [6] and [24]. Both compared noise assessments using static and dynamic traffic models in 23
conjunction with a noise emission model for a urban network (single intersection or corridor). Both 24
concluded that STA is not refined enough to guarantee precise noise estimates and that incorporating 25
more dynamics could improve noise estimates. The French ‘Guide du Bruit’ (also used for the AR-26
INTERIM-CM model) is used in [6] and [24] the FHWA noise emission database, in which the sound 27
power level depends on vehicle speed and throttle conditions (cruising or full). Earlier [25] concluded that 28
a macroscopic dynamic traffic model is sufficient for noise assessments in urban traffic conditions but can 29
be improved by using microscopic traffic models. Research in which DTA models are used in conjunction 30
with noise emission models are however less available. Examples are [23,26,27]. 31
In general, literature agrees that traffic dynamics are important for emission estimates and that the 32
estimates improve when DTA models are used instead of STA models. Most efforts made use 33
microscopic DTA models in conjunction with emission models. However, as argued by [7], [13] and [19] 1
these type of models are less suitable for large networks because of calculation times and can also result 2
in apparent accuracy when individual vehicle operations are used. As a result still STA is mainly used to 3
estimate emission. However, there is less research available on the possible error of using STA especially 4
when assessing the effects of measures taking into account the limitations of STA for realistic cases. 5
Using STA models can result in predicting congestion on wrong locations namely downstream the 6
bottleneck instead of upstream, STA models do not model spillback and do not take into account that 7
bottlenecks will influence downstream traffic demand. As a result the errors in the emission estimates 8
based on STA can be possibly large especially in oversaturated traffic conditions. 9 10 11 MODELING FRAMEWORK 12 13
Based on the extensive review [7] emission models are selected which are suitable to use in conjunction 14
with DTA models, but which can also be used with STA. The ARTEMIS model for air emissions and 15
AR-INTERIM-CM model for noise emissions were selected. Both models are traffic situation based 16
model which can take traffic dynamics into account. For the comparison an emission module has been 17
developed based on these models and connected with the OmniTRANS transport planning software. The 18
Streamline DTA model [28] is part of this software, which is a second order multi-class macroscopic 19
DTA model with physical queuing and spillback. The used STA as well as DTA uses junction modeling 20
to take delay at intersection into account. 21
The estimation of air emissions is based on the ARTEMIS model which distinguishes 4 traffic 22
conditions to calculate hot engine emissions (cold start emissions can be calculated but are not considered 23
in this comparison, assuming these to be equal in STA and DTA). Based on Dutch figures on the 24
distribution of vehicle types for different base and forecast years weighted emission factors are derived 25
per vehicle class, for 10 road types and these 4 traffic conditions. Figure 1 shows the resulting weighted 26
NOx emission factors for passenger cars (2008). The calculation of emissions based on the outcome of the
27
DTA model is done by determining the traffic condition based on calculated speed and fundamental 28
diagram (input of DTA model) per link for each time interval. The traffic condition determines the which 29
emission factors should be used and based on the volumes and composition for that same time interval the 30
emissions are calculated. The total emission of a link is the summation over all time intervals: 31
(
)
( ) ( ), , s s free C a am ad md am am am a t m d E =∑∑∑
q t δ E v t v v l (3) with: ( ) amq t : Volume vehicle class m for link a at time t (veh) ( )
am
v t : Average speed of vehicle class m on link a at time t (km/h)
free am
v : Free flow speed of vehicle class m on link a (km/h)
C am
v : Speed at capacity of vehicle class m on link a (km/h) ( )
s md
E ⋅ : Emission factor of substance s for vehicle class m, depending on speed (grams/(veh*km))
a
l : Length of link a (km)
ad
δ : Road type indicator, equals 1 if link a is of road type d, and 0 otherwise 32
For the calculation based on STA the average traffic condition for the period considered (e.g. morning 33
peak hour) is based on the V/C-ratio. Note, that in this case the outflow capacity, based on junction 34
modeling, is used to determine the V/C-ratios on urban links at intersections if the outflow capacity is 35
lower than the link capacity. The emission factors are interpolated (see Figure 3), to avoid large 36
(unrealistic) changes in emissions in the case of small changes in the V/C-ratio. The interpolation method 37
is also used in the Dutch prescribed standard method to estimate emissions. The boundaries for the traffic 38
condition bins are based on Dutch default values used in current practice. The total emission of a link is: 39
(
, ,)
s s outflow a am ad md a a a a m d E =∑∑
Q δ E Q C C l (4) with: amQ : Volume vehicle class m for link a (veh)
a
Q : Total volume all vehicle classes for link a (pcu)
a
C : Capacity for link a (pcu)
outflow a
C : Outflow capacity for link a (pcu) ( )
s md
E ⋅ : Emission factor of substance s for vehicle class m, depending on V/C (grams/(veh*km))
a
l : Length of link a (km)
ad
δ : Road type indicator, equals 1 if link a is of road type d, and 0 otherwise 1 0.3 E m is s io n fa c to r Stop & Go Saturated Heavy Free flow V/C-ratio 0.0 0.6 0.7 0.8 0.9 1.0 2
Figure 3. Interpolation emission factors when applying STA 3
4
For noise two emission functions are used and derived from the AR-INTERIM-CM model: one for 5
saturated and one for non-saturated conditions. In Figure 2 the noise emission functions for passenger cars 6
are shown. The AR-INTERIM-CM is based on the French ‘Guide du Bruit’ and like the ARTEMIS 7
model developed in a project of the European Union. The emission functions are road type independent, 8
however correction factors can be used to correct for road surfaces. There are different emission functions 9
per vehicle class. The calculation of emissions based on the outcome of the DTA model is done by 10
determining the traffic condition based on calculated speed and fundamental diagram (input of DTA 11
model) per link for each time interval. The traffic condition determines which emission function should 12
be used and based on the traffic volumes and traffic composition (energetically summed) for that same 13
time interval the noise emissions are calculated. The average sound power level (i.e. noise emission) of a 14
link is calculated by energetically averaging the emissions over time: 15
( )
(
)
(
)
ref(
)
( ) 10 ( ) ( ) ( ) log 10 log ( ) 10 10 log am am am m am m am am m L t t m a v t L t v t v t q t v t L T α β = + + ∆ = ∑∑
(5)with:
( ) am
q t : Volume vehicle class m for link a at time t (veh/hour)
( ) am
v t : Average speed of vehicle class m on link a at time t (km/h) ( )
am
L t : Sound power level of vehicle class m on link a at time t (dB(A)) a
L : Average sound power level for link a (dB(A))
,
m m
α β : Parameters dependent on vehicle category and traffic condition for noise calculations ref
m
v : Reference speed dependent on vehicle category (km/h)
t
∆ : Time interval output data DTA model (h) T : Total time period (h)
1
For the calculation based on STA there are two options. The first option uses the posted speed limit (in 2
our case equal to the free flow speed) which is often used in current practice. The second option uses the 3
calculated speed based on the travel time functions. In the latter case interpolation is used by energetically 4
weighing the two emission functions (default between VC-ratios of 0.8 and 1.0). The average sound 5
power level (i.e. noise emission) of a link is calculated by: 6
Based on posted speed
( )
( )
ref(
)
10 log 10log 10 log 10 am free free free am am m am m am am m L a m v L v v Q T v L α β = + + = ∑
(6)Based on calculated speed
( )
( )
(
)
( ) ref ( ) 10 10 10 log 10 log 10log 10 , if , 10 log 10 1 10 , if , , if am freeflow saturated a a am am m am m am am m L a m freeflow a a a a a L E a a a a a saturated a v L v v Q T v L Q L C a Q Q b C b C Q L a C b b a b a L α β = + + = < − − = − + − − < < ∑
i i . a a Q b C > (7) with: amQ : Volume vehicle class m for link a (veh)
a
Q : Total volume all vehicle classes for link a (pcu)
a
C : Capacity for link a (pcu)
am
v : Average speed of vehicle class m on link a based on travel time function (km/h)
am
L : Sound power level of vehicle class m on link a (dB(A)) a
L : Average sound power level for link a (dB(A))
,
m m
α β : Parameters dependent on vehicle category and traffic condition for noise calculations ref
m
v : Reference speed dependent on vehicle category (km/h)
,
CASES: HIGHWAY AND URBAN NETWORK 1
2
To determine the consequences of using STA versus DTA to estimate emissions the presented methods 3
are applied on a highway corridor, namely the highly congested A12 highway between Gouda and The 4
Hague and a moderate congested urban corridor the Nassaukade in Amsterdam (see Figure 4) with 5 5
traffic controlled intersections. Network and demand were derived from existing calibrated models of the 6
regions (i.e. morning peak of Spitsmijden model respectively evening peak of tactical model Amsterdam). 7
For the sake of clarity corridors are selected to present the results of the comparison to assure the same 8
route choice behavior in both models. This means that in the cases we focus on the difference as a result 9
of the differences between static and dynamic loading and therefore the differences in capturing the traffic 10
flow phenomena described earlier. Total demand in STA and DTA are set to be equal, however the DTA 11
model uses a distribution over time. For the urban case of Amsterdam the demand is slightly changed, 12
which means that the outcomes are hypothetical. The green splits of the signal plans of the 5 intersections 13
are fixed and the same for both models, which means that the capacity given to each direction is exactly 14
the same in the STA as well as the DTA model. Because the differences between STA, measurements and 15
DTA in emissions comprises the difference in loads and traffic dynamics, the total air emissions were also 16
corrected for the differences in total loads per link. This means that for the corrected effects the weighted 17
average emission factors are calculated. 18
19
4a: Highway network 20
4b: Urban network
1
Figure 4. Used cases 2
3
Results highway network 4
For the highway network also measured flow and speeds (using detector data) are available for all 5
highway links of the morning peak. Note that these measurements were not used to calibrate the model. 6
These measurements were also used to estimate emissions using the ARTEMIS and AR-INTERIM-CM 7
model. In the remainder of this paper we will refer to these results as measured emissions. Figure 5 shows 8
the results of the total emissions of all highway links (corrected and not corrected) as well as the 9
differences per link (only corrected, direction of travel from right, link 15 near Gouda, to the left, link 1 10
near The Hague). In both figures the measured emissions are the reference case (index=100), which 11
means that if the value is closer to 100 the estimated emission based on the model is closer to the 12
measurements. 13
The results show that the average speed and link load as well the emission estimates based on the 14
DTA model are closer to the measurements than those based on the STA model. The total emission of 15
NOx estimated by the STA is for example 23% higher than the measured emissions (corrected 11%
16
higher) and for PM10 13% higher (corrected 2%). The estimates based on the DTA model are 8%
17
(corrected 4%) higher for NOx and for PM10 5% (corrected 1%). For noise a comparison is made in which
18
the noise emissions are calculated based on the posted speed limit as well as the calculated (and 19
measured) speed. The average noise emissions based on the posted speed limit is 5% (in dB(A)) higher 20
compared to measured emissions. Note that, noise is measured on a logarithmic scale. This means that the 21
average noise emission is 4 dB(A) lower in absolute values, which equals a reduction of traffic volume of 22
more than 50%. When the calculated speed is used, the DTA model performs slightly better than the STA 23
model (0.3% versus 0.5%). 24
Differences model versus measured
Total of highway links, total peak period
85 90 95 100 105 110 115 120 125
Loads Speed CO2 NOx PM10 dB(A)
Indicator In d e x ( M e a s u re d = 1 0 0 )
Static average emission factor Dynamic average emission factor
Static total Dynamic total
Static sound calculated speed
1
Static versus dynamic
Average weighted NOx emission factor, individual links, total peak period
60 80 100 120 140 160 180 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Linknumber (numbered West - East)
In d e x ( m e a s u re d = 1 0 0 ) Static Dynamic 2
Figure 5. Results assignments on emissions A12, Gouda –The Hague 3
4
For the individual links the differences are larger and it is even possible that the STA estimates are lower 5
than measured emissions while the DTA estimates are higher and vice versa. At almost every link the 6
DTA model performs better. There is one link where the DTA model performs significantly worse 7
because there DTA model predicts a totally different traffic condition than the measurements indicate. 8
The differences between the DTA and measurements is mainly caused by the differences in demand and 9
resulting traffic conditions. If the DTA model would have been calibrated using the measurements the 10
differences would be much smaller. The comparison between the estimates based on DTA and STA 11
therefore provide an indication of the possible error which is being made when using a STA. This 12
comparison shows that the emissions of NOx are between 30% lower and 42% higher (corrected 38%
13
lower and 28% higher) when a STA is used. For PM10 the differences are smaller but still 15% lower or
14
15% higher (corrected 10% lower and 7% higher). For CO2 the differences are locally even larger than for
15
NOx, however for CO2 the local emissions are not relevant. For noise the sound emissions can be 8%
16
lower based on the DTA model compared to STA using the posted speed limit. Note that higher is not 17
possible on highways. If the calculated speeds based on the travel time functions are used for the STA it is 18
possible that the noise emissions can be 5% higher or 6% lower compared to the DTA model (in absolute 19
values these are enormous differences of over 4 db(A)). These differences found in the estimates of 20
emissions for peak hours are large and can have a significant influence even when translated to 24 hour 1
values for policy decisions (e.g. within cost-benefit analysis or comparison with limit values). 2
3
Effect measure
Total network, total peak hour
80 85 90 95 100 105 110
CO CO2 NOx PM10 Fuel HC dB
maxspeed dB calcspeed Indicator In d e x ( re fe re n c e = 1 0 0 ) Static Dynamic 4
Figure 6. Effect measures on total emissions based on STA and DTA 5
6
The differences between DTA and STA are also of interest when the effects of a measure is assessed. We 7
assume that an additional lane is realized between Zoetermeer and Nootdorp and demand stays fixed. In 8
this case the STA estimates a decrease of emissions and also only for those links where the capacity is 9
increased. Within the DTA model also other links are affected by this measure, because the interaction 10
between links is considered. The DTA shows that for some links emissions increase (especially 11
downstream) and for some links the emissions decrease (especially upstream). On the total network the 12
DTA model estimates an increase of PM10 emissions while the STA model estimates a decrease of PM10
13
emissions. This indicates that not only the location at which the STA predicts an effect is unreliable, but 14
also the extent and even the direction of effect. 15
16
Results urban network 17
For the urban case only the comparison between STA and DTA are presented. Figure 7 shows the total 18
emissions and average emission factors. The comparison shows that for the urban case the emission 19
estimates using the STA model are lower than the DTA model. One of the reasons is that in the urban 20
case the junctions are determinative for the traffic conditions. The possible congestion problems as a 21
result of these junctions are taken into account by translating them to the first upstream link. However, the 22
STA model can not model spillback, which means that in saturated conditions and relative short link 23
length the STA model can underestimate the congestion problems. On the other hand, STA models do not 24
take into account that traffic is metered by upstream bottlenecks, which means that traffic demand can be 25
overestimated in STA models at specific locations. Because in our urban case the total traffic demand in 26
this corridor network is chosen to be the same for STA and DTA and the congestion problems are 27
moderate (i.e. only two links with V/C ratio’s above 1, namely 1.07 and 1.03) the emission estimates of 28
the STA model are lower. As a result the differences between STA and DTA in total emissions and 29
average emission factors are similar. As depicted in Figure 7 the NOx emissions are 9% lower based on
1
the STA and the PM10 emissions 1% lower. The weighted average sound emissions are in this case equal
2
when based on the calculated speed instead of the posted speed limit. 3
The local differences can be larger. The total NOx emissions can be 10% higher or 23% lower for
4
STA compared with DTA (corrected for load differences 5% higher or 24% lower. For PM10 the total
5
emissions can be 11% higher or 5% lower (corrected 3% higher or 6% lower). For sound the emissions 6
based on the calculated speed versus the posted speed limit (both for STA) can be 2% higher 7
(approximately 1 dB(A) higher) or 3% lower (approximately 2 dB(A) lower). The difference between 8
STA based on the posted speed limit and DTA can be 2% higher or 4% lower. This shows that in the 9
urban case the emissions based on the calculated speed can be higher than based on the posted speed 10
limit. 11
12
Static versus dynamic
Total network, total peak period
80 85 90 95 100 105 110
CO CO2 NOx PM10 Fuel HC dB
maxspeed dB calcspeed Indicator In d e x ( D y n a m ic = 1 0 0 ) Average emissionfactor Total emissions 13
Figure 7. Results assignments urban case Amsterdam 14
15
When introducing a measure, namely realizing an additional lane for the direction with the highest V/C-16
ratio (i.e. traveling from north to south, first intersection (V/C-ratio 1.07)), similar results are found as in 17
the highway case. This means that the STA model only predicts a lower V/C-ratio for the link with an 18
additional lane, while for the DTA model this measure influences the traffic conditions on many links. In 19
the STA case the emissions slightly decrease on network level (less than half percent for all substances), 20
while in the DTA model the emissions slightly increase (approximately half percent dependent on 21
substance) because of worsened traffic conditions downstream. These changes are also small because, the 22
green split of all intersections remains the same and congestion problems are moderate. For the DTA 23
model traffic conditions can only change downstream of the second intersection due to changes in 24
composition of traffic (i.e. distribution of directions). Although, the effects on network level are small, in 25
this case the STA model predicts a reduction for all substances, while the DTA predicts an increase. 26
Again the difference in the estimated effect by STA or DTA on a local level can be large. In 27
Figure 8 the estimated effect for links on corridor from north to south which are numbered from 1 to 11. 28
As shown, the STA model only predicts an effect on a single link (i.e. link 5). The STA model predicts a 29
reduction in total emissions on this link, while the DTA model predicts an increase in total emissions 30
which is mainly due to an increase of volume. This is also the main reason of the increase in emissions of 31
link 2, 3 and 4, while for link 6, 7, 8 and 9 (i.e. downstream) the increase is also due to deteriorating 32
traffic conditions. This means that also the urban case shows that location, extent and direction of effect 33
estimated by STA can be unreliable, given the realism of estimating traffic conditions using both types of 1
assignment models. 2
3
Effect measure on total NOx emissions individual links
75 80 85 90 95 100 105 110 115 1 2 3 4 5 6 7 8 9
Linknumber (numbered north - south)
In d e x ( re fe re n c e = 1 0 0 ) NOx Static NOx Dynamic 4
Figure 8. Effect measures on total emissions based on STA and DTA 5
6 7
DISCUSSION AND CONCLUSIONS 8
9
In current practice mainly STA models are used for the estimation of emissions. In many countries there 10
are regulations or protocols on specific models to use and also often the prescribed use of these STA 11
models. However the unreliability of this model type at link level can be large especially in saturated 12
traffic conditions. Eventually this can result in erroneous policy decisions or implementing mitigating 13
measures on wrong locations or not implementing them on locations where they should. DTA models can 14
take the traffic flow phenomena into account which are important explanatory variables for emissions. 15
The current generation of DTA models and the increasing computation power (e.g. using parallel 16
computing) has improved the scalability of this model type significantly and are therefore capable in 17
assessing large networks. However, computation times may remain an issue. Developments in quasi-18
dynamic assignments like STAQ [2] are therefore of interest because they further improve the scalability 19
of DTA while improving the realism of STA. Comparison of emission estimates with these type of 20
models is therefore of interest. 21
For noise emissions often the operational (free) speed or posted speed limits are used. For 22
highways this means that the estimates are often to conservative (i.e. can be more than 4 dB(A) lower in 23
reality). However for urban roads this can result in an underestimate of emissions, depending of the 24
source model used and way of application. The urban case showed that noise emissions can be 1 dB(A) 25
higher. For emissions the estimated traffic conditions based on STA can result in large over- and 26
underestimates compared to DTA (e.g. for NOx up to 42% higher and 30% lower for the highway case
27
and 10% higher and 23% lower for the urban case). In general, the differences will be larger when the 28
congestion problems are larger, because the limitations of STA in capturing realistic traffic flow 29
phenomena is mainly associated with these traffic conditions. The possible errors made may be even 30
larger than found in the presented cases, because the STA model does not take into account that upstream 31
bottlenecks not part of the cases presented in this research influence local demand. In addition, the 32
behavioral responses, not considered in this comparison, in both model types can be different as well 1
because of differences in the estimated effects of the measures taken on travel costs. As a result the errors 2
in the emission estimates based on STA can be even larger. In addition, intersections which are more 3
difficult to take into account in emission estimates based on STA models, are determinative for traffic 4
handling on urban networks. When a measure is taken to alleviate congestion problems it is shown in 5
both cases that not only the location at which the STA predicts an effect is unreliable, but also the extent 6
and even the direction of effect. 7
The peak hours are only a part of the day and for noise the night is for example more important 8
for the nuisance than the peak hours. However, given the found differences between STA and DTA, these 9
can still be significant even for the consequences on a daily average. Improvements in the incorporation 10
of traffic dynamics in the current practice of emission estimates are therefore important to assess the 11
livability and health effects of traffic measures. 12
13 14
REFERENCES 15
1. Smit, R (2006). An examination of congestion in road traffic emission models and their application 16
to urban road networks. Doctoral dissertation, School of Environmental Planning. Faculty of 17
Environmental Sciences. Griffith University, Queensland, Australia. 18
2. Bliemer, M, L. Brederode, L. Wismans and E.S Smits (2012). Quasi-dynamic traffic assignment: 19
static traffic assignment with queueing and spillback. In proceedings 91th Annual Meeting of the 20
Transportation Research Board, January 22-26, Washington, D.C., USA. 21
3. Rakha, H, and Ahn, K. (2003). Integration Modelling Framework for Estimating Mobile Source 22
Emissions. Journal of transportation engineering, 130, pp. 183-193. 23
4. Beek, W. Derriks, H., Wilbers, P., Morsink, P., Wismans, L. and Van Beek, P. (2007). The effects 24
of speed measures on air pollution and traffic safety. In: Proceedings European Transport 25
Conference 2007, 17-19 October 2007 (Noordwijkerhout, The Netherlands, Leeuwenhorst 26
Conference centre). 27
5. Barth, M. and Boriboonsomsin, K. (2008). Real-World CO2 Impacts of Traffic Congestion. In:
28
Proceedings 87th Annual Meeting of Transportation Research Board, 13-18 January 2008 29
(Washington, DC, Transportation Research Board). 30
6. Can, A., Leclercq, L., Lelong, J. and Defrance, J. (2009). Accounting for traffic dynamics improves 31
noise assessment: Experimental evidence. Applied Acoustics, Vol. 70, pp. 821–829. 32
7. Wismans, L. J. J., E.C. van Berkum and M. C. J. Bliemer (2011). Modelling Externalities using 33
Dynamic Traffic Assignment Models: a review. Transport Reviews, Volume 31, Issue 4, July 2011, 34
pages 521-545 (ISSN 0144-1647). DOI:10.1080/01441647.2010.544856. 35
8. Nesami, K.S., L. Chu, M.G. McNally and R. Jayakrishnan (2007). Estimation of vehicular emissions 36
by capturing traffic variations. Atmospheric Environment, 41, pp. 2996-3008. 37
9. CNOSSOS (2010). Draft JRC Reference Report on Common NOise ASSessment MethOdS in EU 38
(CNOSSOS-EU). To be used by the EU Member States for strategic noise mapping after adoption as 39
specified in the Directive 2002/49/EC. Version 3, November 2010 40
10. INFRAS (2007). Artemis: Assessment and reliability of transport emission models and inventory 41
systems. ROAD EMISSION MODEL – Model Description. Workpackage 1100. Deliverable No. 13. 42
11. TRB, 1997. Improving transportation data for mobile-source emissions estimates. In: Chatterjee, A., 43
Miller, T., Philpot, J., Wholley, T., Guensler, R., Hartgen, D., Margiotta, R., Stopher, P. (Eds.), 44
Transportation Research Board. National Research Council, National Academy Press, Washington, 45
DC, NCHRP Report No. 394. 46
12. Smit, R., M. Poelman and J. Schrijver (2008). Improved road traffic emission inventories by adding 47
mean speed distributions. Atmospheric Environment, 42, pp. 916–926. 48
13. TRB, 2005. Predicting air quality effects of traffic flow improvements. Final report and user’s guide. 49
In: Dowling, R., Ireson, R., Skarbadonis, A., Gillen, D. and Stopher, P., Transportation Research 50
Board. National Cooperative Highway Research Program, Washington, DC., NHCRP Report No. 1
535. 2
14. Ahn, K. and Rakha, H. (2008). The effects of route choice decisions on vehicle energy consumption 3
and emissions. Transportation Research Part D, Vol. 13, pp. 151-167. 4
15. Akcelik, R. and Besley, M. (2003). Operating cost, fuel consumption, and emission models in 5
aaSidra en aaMotion. Proceedings 25th conference of Australian Institutes of Transport Research, 6
Adelaide. 7
16. Kun, C. and Lei, Y. (2007). Microscopic Traffic-Emission Simulation and Case Study for 8
Evaluation of Traffic Control Strategies. J Transpn Sys Eng & IT, Vol. 7(1), 93−100. 9
17. Int Panis, L., Broekx, S., Liu, R. (2006). Modelling instantaneous traffic emission and the influence 10
of traffic speed limits. Science of the Total Environment, Vol. 371, pp. 270–285. 11
18. Mensink, C. and Cosemans, G. (2008). From traffic flow simulations to pollutant concentrations in 12
street canyons and backyards. Environmental Modelling & Software Vol. 23, pp 288 – 295. 13
19. Aguiléra, V. and M. Lebacque (2010). Dynamic estimation of traffic emissions in metropolitan road 14
networks. In proceeding 12th WCTR, July 11-15, Lisbon, Portugal. 15
20. Bai, S., Chiu, Y.C. and Niemeier, D.A. (2007). A comparative analysis of using trip-based versus 16
link-based traffic data for regional mobile source emissions estimation. Atmospheric Environment, 17
Vol. 41, pp. 7512–7523. 18
21. AR-INTERIM-CM (2003). Adaptation and revision of the interim noise computation methods for 19
the purpose of strategic noise mapping. WP 3.1.1 Road traffic noise – description of the calculation 20
method and WP 3.1.2 Road traffic noise – Noise emission: databases. Wolfel et al. 21
22. Leclercq, L. (1999). Use of dynamic traffic modelling to improve estimate of noise emissions. A 22
feasibility study. Recherche transports securite, Vol. 65, pp. 60–77. 23
23. De Coensel., B., De Muer, T., Yperman, I. and Botteldooren, D. (2005). The influence of traffic 24
flow dynamics on urban soundscapes. Applied Acoustics, Vol. 66, pp. 175 – 194. 25
24. Chevallier, E., Can, A., Nadji, M. and Leclercq, L. (2009). Improving noise assessment at 26
intersections by modeling traffic dynamics. Transportation Research Part D, Vol. 14, pp. 100–110. 27
25. Can, A., Leclercq, L. and Lelong, J. (2008). Dynamic estimation of urban traffic noise: Influence of 28
traffic and noise source representation. Applied Acoustics, Vol. 69, pp 858 - 867. 29
26. King, E.A. and Rice H.J. (2009). The development of a practical framework for strategic noise 30
mapping. Applied Acoustics, doi:10.1016/j.apacoust.2009.01.005. 31
27. Bhaskar, A., Chung, E. and Kuwhara, M. (2007). Development and implementation of the areawide 32
Dynamic Road traffic NoisE (DRONE) simulator. Transportation Research Part D, Vol. 12, pp. 33
371 – 378. 34
28. Raadsen, M.P.H., H.E. Mein, M.P. Schilpzand and F. Brandt (2010). Implementation of a single 35
dynamic traffic assignment model on mixed urban and highway transport networks including 36
junction modeling. In proceedings DTA symposium, Takayama, Japan, 2010. 37