Evaporative Emissions from Onroad
Vehicles in MOVES5
£% United States
Environmental Protect
Agency
-------
Evaporative Emissions from Onroad
This technical report does not necessarily represent final EPA decisions
or positions. It is intended to present technical analysis of issues using
data that are currently available. The purpose in the release of such
reports is to facilitate the exchange of technical information and to
inform the public of technical developments.
Vehicles in MOVES5
Assessment and Standards Division
Office of Transportation and Air Quality
U.S. Environmental Protection Agency
NOTICE
4>EPA
United States
Environmental Protection
Agency
EPA-420-R-24-014
November 2024
-------
Table of Contents
1 Introduction and Background 1
2 Test Programs and Data Collection 6
3 Design, Analysis and MOVES Inputs 8
3.1 Fuel Tank Temperature Generator 8
3.1.1 Fuel Tank Temperature for Hot and Cold Soaks 10
3.1.2 Fuel Tank Temperature while Running 12
3.2 Permeation 14
3.2.1 Temperature Adjustment 17
3.2.2 Fuel Adjustment 18
3.3 Tank Vapor Venting 19
3.3.1 Altitude 21
3.3.2 Cold Soak 21
3.3.3 Hot Soak 37
3.3.4 Running Loss 45
3.4 Inspection/Maintenance (I/M) Program Effects 48
3.4.1 I/M Factor (Relative Program Effectiveness) 49
3.4.2 Leak Prevalence 52
3.5 Liquid Leaks 53
3.6 Refueling 54
3.6.1 Vapor Loss 57
3.6.2 Spillage 61
Appendix A Notes on Selected Evaporative Emission Data 62
Appendix B Tank Fuel Generator 64
Appendix C List of Acronyms 68
Appendix D Glossary 69
-------
References
-------
1 Introduction and Background
EPA's Office of Transportation and Air Quality (OTAQ) has developed the Motor Vehicle
Emission Simulator (MOVES). The MOVES model estimates emissions for mobile sources
covering a broad range of pollutants and allows multiple scale analysis. MOVES currently
estimates emissions from cars, trucks and motorcycles.
Evaporative processes can account for a significant portion of gaseous hydrocarbon emissions
from gasoline vehicles. Volatile hydrocarbons evaporate from the fuel system while a vehicle
is refueling, parked or driving. MOVES does not include estimates for emissions from non-
fuel sources such as window washer fluid, paint, plastics, and rubber. Evaporative processes
differ from exhaust emissions because they don't directly involve combustion, which is the
main process driving exhaust emissions. For this reason, evaporative emissions require a
different modeling approach. In the previous MOBILE models and in certification test
procedures, evaporative emissions were quantified by the test procedures used to measure
them:
• Running Loss - Vapor lost during vehicle operation.
• Hot Soak - Vapor lost after turning off a vehicle.
• Diurnal Cold Soak - Vapor lost while parked at ambient temperature.
• Refueling Loss - Vapor lost and spillage occurring during refueling.
For MOVES, we instead model the underlying physical processes involved in evaporation of
fuels. This "modal" approach characterizes the emissions by different emissions generation
processes, each having its own engineering design characteristics and failure rates. This way,
certain physical processes can be isolated, for example, ethanol has a unique effect on
permeation, which occurs in all the above modes. The approach used in MOVES categorizes
evaporative emissions based on the evaporative mechanism, using the following processes:
• Permeation - The migration of hydrocarbons through materials in the fuel system.
• Tank Vapor Venting (TVV) - Vapor generated in fuel system lost to the
atmosphere, when not contained by evaporative emissions control system.
• Liquid Leaks - Liquid fuel leaking from the fuel system, ultimately evaporating.
• Refueling Emissions - Spillage and vapor displacement as a result of refueling.
Figure 1 illustrates the evaporative emission processes. Permeation occurs continuously
through the tank walls, hoses, and seals. It is affected by fuel tank temperature and fuel
properties.
1
-------
Figure 1: Illustration of Evaporative Processes
Vapor^/enting
Vapor Leaks
Vapor
"\
0
Fuel
J
Permeation
*
Liquid Leaks
These processes occur in each operating mode (Running Loss, Hot Soak, Cold Soak) used in
the MOVES model. Each emission process can be modeled over a user-defined mix of
operating modes, shown in Table 1. This makes for more accurate modeling of scenarios that
do not replicate test procedures. The processed values for the evaporative emission processes
used by MOVES are shown in Table 2.
Table 1: MOVES Evaporative Process Operating Modes
opModelD
Operating mode description
150
Hot Soaking
151
Cold Soaking
300
Engine Operation
Table 2: MOVES Evaporative Emission Processes
processID
Emission process description
11
Evap permeation
12
Evap vapor venting losses
13
Evap liquid leaks
18
Refueling displacement vapor losses
19
Refueling fuel spillage
2
-------
Evaporative emissions are a function of many variables. In MOVES, these variables
include:
• Ambient temperature
• Fuel Tank temperature
• Model year group (as a surrogate for technology and certification standard)
• Vehicle age
• Vehicle class
- Passenger Vehicle
- Motorcycle
- Short/Long-haul Trucks
• Fuel Properties
- Ethanol content
- Reid Vapor Pressure (RVP)a
• Failure Modes
• Presence of inspection and maintenance (I/M) programs
Both ambient temperature and engine operation cause increases in fuel tank temperature. An
increase in fuel tank temperature generates more vapor in the tank. Activated charcoal
canisters are a control technology commonly used to adsorb the generated vapor. During
engine operation, the canister is purged periodically, and the captured vapor is diverted to the
engine and burned as fuel. The emission certification standards for a vehicle (associated with
model year and vehicle class) influence the capacity of the canister system. When the
generated vapor exceeds the capacity of the canister, the vapor is vented to the atmosphere.
This can occur when a fuel undergoes a large ambient temperature increase, or if a fuel with
higher volatility is used, or when a vehicle canister collects vapor for many days without
purging. In calculating vapor venting, MOVES accounts for co-mingling ethanol and non-
ethanol gasoline, and for RVP weathering of in-use fuel. Details on these Tank Fuel Generator
calculations are provided in Appendix B.
Fuel systems can develop liquid and vapor leaks that circumvent the vehicle emissions control
system. Some inspection and maintenance (I/M) programs explicitly intend to identify
vehicles in need of evaporative system repairs. Some states also implement Stage II programs
at gas stations to capture the vapors released during refueling. These programs capture
refueling vapor with technology installed at the pump rather than internal to the vehicle.
The model year groups for evaporative emissions are shown in Table 3. They reflect
evaporative emission standards and related technological improvements. Early controls
included the introduction of activated charcoal canisters for controlling fuel vapor emissions.
Later controls included fuel tanks and hoses built with more advanced materials less prone to
a The MOVES fuel supply table provides the characteristics of gasoline sold in each county and month. For vapor venting
calculations, the MOVES Tank Fuel Generator uses the fuel supply information to account for the effects of "comingling" ethanol
with non-ethanol gasoline and for the "weathering" effect on RVP for in-use fuel. See appendix for details.
3
-------
permeation. Also, reduction of fittings and connections became an important consideration for
vapor mitigation.
Table 3: Model Year Groups for Evaporative Emissions in MOVES
Model year group
Evaporative emissions standard or technology level
1971-1977
Pre-control
1978-1995
Early control
1996
80% early control, 20% enhanced evap
1997
60% early control, 40% enhanced evap
1998
10% early control, 90% enhanced evap
1999-2003
100% Enhanced evap
2004-2015
Tier 2, LEV II
2016-2017
40% Tier 3
2018-2019
60% Tier 3
2020-2021
80% Tier 3
2022 +
Tier 3
This report documents the evaporative emission rates measured in terms of total hydrocarbons
(THC). Total hydrocarbon gases are defined as the measurement of gaseous hydrocarbons by a
flame ionization detector (FID). Evaporative emissions also contain oxygenated hydrocarbons
such as alcohols and aldehydes.
MOVES estimates organic gas aggregate species (e.g., Volatile Organic Compounds, Total
Organic Gases) from the THC emissions as documented in the speciation report.1 MOVES
estimates specific hydrocarbon species as fractions of VOC and TOG emissions. Eight
important mobile source air toxics (MSATs), including benzene and ethanol, are calculated
from evaporative VOC emissions as documented in the air toxics report.2 Evaporative
emissions are not directly affected by the combustion process, and does not estimate any
emissions from combustion products. Table 4 contains a list of the evaporative pollutants
calculated by MOVES. MOVES calculates additional chemical mechanism species from
evaporative emissions used for air quality modeling as documented in the speciation report.1
The data used for this evaporative analysis was collected on light-duty gasoline vehicles but
were also applied to heavy-duty gasoline vehicles since heavy-duty gasoline data was not
available at time of analysis.
4
-------
Table 4: MOVES Evaporative Pollutants
pollutantID
pollutantName
NEIPollutantCod
shortName
1
Total FID Hydrocarbons
HC
THC
20
Benzene
71432
Benzene
21
Ethanol
ETOH
40
2,2,4-Trimethylpentane
540841
2,2,4-Trimethylpentane
41
Ethyl Benzene
218019
Ethyl Benzene
42
Hexane
206440
Hexane
45
Toluene
85018
Toluene
46
Xylene
123386
Xylene
79
Non-Methane Hydrocarbons
NMHC
NMHC
80
Non-Methane Organic Gases
NMOG
NMOG
86
Total Organic Gases
TOG
TOG
87
Volatile Organic Compounds
VOC
VOC
185
Naphthalene gas
91203
Naphthalene Gas
For diesel vehicles, it is assumed that there are no evaporative emission losses except for
refueling spillage. Due to the low vapor pressure of diesel fuel, diesel evaporative losses are
considered negligible.
At the time of this analysis, there was no relevant evaporative emissions data for compressed
natural gas (CNG) vehicles. CNG fuel systems and refueling procedures are significantly
different from those of liquid petroleum-based fuels. For the current release of MOVES, all
evaporative emission rates for CNG vehicles are set at zero.
We significantly updated the evaporative emission calculations and rates in MOVES2014
based on updated emissions data, failure rates, and vehicle activity in MOVES2014. Because
of the significant updates, the MOVES2014 version of this report was subject to peer review
under EPA's peer review guidelines. More information about this peer review, including peer
reviewer comments and EPA response is available on the web.3
Evaporative emission inputs for MOVES2014 were also reviewed by the Coordinating
Research Council.4 Based on our evaluation, most of the issues pointed out in the CRC report
are expected to have very little impact on the magnitude of the evaporative emissions
computed by MOVES. However, we continue to look for opportunities to improve how
MOVES estimates evaporative emissions.
Updates for MOVES3 and MOVES4 changed the calculation of refueling emissions as
explained in Section 3.6 below. For MOVES5, we updated information on Stage II vapor
recovery programs, explained in the same section.
As explained in the MOVES Population and Activity report5, the activity associated with
evaporative emissions remains the same as in MOVES2014.
5
-------
2 Test Programs and Data Collection
The modeling of evaporative emissions in MOVES is based on data from a large number of
studies (Table 5). Over a decade of research greatly modernized evaporative emissions
modeling. New test procedures provided modal emissions data that greatly advanced the state
of the science. For example, the CRC E-77 test programs6'7'8'9 measured permeation emissions
separately from vapor emissions. Implanted leak testing from these studies along with further
field research provided the first large database regarding the prevalence and severity of
evaporative leaks and other mal- functions. The studies applied an innovative sampling design
which preferentially recruited high emissions vehicles with the aid of infrared ultraviolet
remote sensing devices. The field studies used a portable test cell (PSHED) to measure in-use
hot soak emissions on a large number of vehicles. Findings from these studies were introduced
in MOVES2014 with the explicit modeling of vapor leaks.10'11
Appendix A has a more detailed summary of these test programs.
6
-------
Table 5: Evaporative Emission Research Programs
Program
# of Vehicles
CRC E-9
Measurement of Diurnal Emissions from In-Use
Vehicles12
151
CRCE-35
Measurement of Running Loss Emissions in In-Use
Vehicles13
150
CRCE-41
Evaporative Emissions from Late-Model In-Use
Vehicles14'15
50
CRCE-65
Fuel Permeation from Automotive Systems16
10
CRC E-65-3
Fuel Permeation from Automotive Systems: EO, E6,
E10, and E8517 1717[25]
10
CRCE-77
Vehicle Evaporative Emission Mechanisms: A Pilot
Study6
8
CRC E-77-2
Enhanced Evaporative Emission Vehicles7
8
CRC E-77-2b
Aging Enhanced Evaporative Emission Vehicles8
16
CRC E-77-2c
Aging Enhanced Evaporative Emission Vehicles with
E20 Fuel9
16
High Evap field studies1011
Thousands
Fourteen Day Diurnal study18
5
Running Loss Testing with Implanted Leaks19
5
API Leakage Study20
Not Avail.
API Gas Cap Study21
Not Avail.
EPA Compliance Testing22
Thousands
7
-------
3 Design, Analysis and MOVES Inputs
This chapter provides detailed information on how evaporative emissions in MOVES are
calculated, and the analysis used to determine appropriate default inputs for the model.
As emission standards have tightened, fuel system materials and connections have become
more efficient at containing fuel vapors. Purge systems and canister technologies have also
advanced, resulting in less vented emissions. Fuel tank temperature is an important
consideration in modeling permeation and vapor emissions. However, liquid leaks occur
regardless and are not dependent on temperature.
3.1 Fuel Tank Temperature Generator
Fuel tank temperature is closely correlated with permeation and vapor venting as observed in
the CRC E-77 pilot testing program 6. This program tested ten vehicles in model years 1992
through 2007. The results showed that fuel temperature strongly influences evaporative
emissions in all testing regimes. Fuel tank temperature is dependent on the daily ambient
temperature profile and vehicle operation patterns. Modern vehicles (enhanced-evap, 1996 &
later) do not recirculate fuel from the engine to the fuel tank and therefore have a lower
temperature rise than older vehicles during operation. In Figure 2, the permeation emissions
are plotted over a 3-day California diurnal tests with (65-105°F) as the low temperature range and
85-120°F as the high temperature range. Both the effects of temperature and fuel volatility can
be observed.
8
-------
Figure 2: Permeation Temperature and RVP effects
Test Time hours
The MOVES "Fuel Tank Temperature Generator" calculates fuel temperature (also referred to
as fuel tank temperature) for a given ambient temperature profile and vehicle trip schedule
based on the vehicle type and model year. Different equations are used depending on the
operating mode of the vehicle: running, hot soak, or cold soak. Fuel tanks are warmer during
ninning operation than the ambient temperature. The routing of hot exhaust, vehicle speed,
and airflow can all affect tank temperature. Immediately after the engine is turned off, the
vehicle is in a hot-soak condition, and the fuel tank begins to cool to ambient temperature. In
cold soak mode, the vehicle has reached ambient temperature.
Input parameters for the fuel tank temperature generator are:
o Hourly ambient temperature profile (zoneMonthHour table)
o Key on and key off times (sampleVehicleTrip table)
o Day and hour of first KeyON (hourDay table)
o Vehicle Type (Light-duty vehicle, Light-duty truck, Heavy-duty gas truck and
Motorcycles)
o Pre-enhanced or enhanced evaporative emissions control system
The MOVES algorithm iterates through a set of "typical" vehicle trips based on information
from instrumented vehicles. This data is stored in two tables in the MOVES default database.
9
-------
SampleVehicleDay lists a sample population of vehicles, each with an identifier (vehID), an
indication of vehicle type (sourceTypelD), and an indication (daylD) of whether the vehicle is
part of the weekend or weekday vehicle population. To represent vehicles that may sit unused
for multiple days, some vehicles in this table do not have any trips. The second table,
SampleVehicleTrip, lists the trips in a day made by each of the vehicles in the
SampleVehicleDay table. It records the vehID, daylD, a trip number (tripID), the hour of the
trip(hourlD), the trip number of the prior trip (priorTripID), and the times at which the engine
was turned on and off for the trip. The keyOnTime and keyOffTime are recorded in minutes
since midnight of the day of the trip. For more information on the activity data used to
determine the time of keyOn and keyOff events, see the MOVES technical report on vehicle
populations and activity5 and supporting contractor reports.23' 24
Coefficients for the Fuel Tank Temperature Generator are recorded in MOVES in the
tankTemperatureGroups and tankTemperatureRise tables.
3.1.1 Fuel Tank Temperature for Hot and Cold Soaks
Equation 3-1 is used to model tank temperature as a function of ambient temperature.
TTank is the fuel tank temperature, Tair is the ambient temperature, and k is a constant
proportionality factor (k = 1.4 hr ). The value of k was established from EPA compliance
data on 77 vehicles that underwent a 2-day diurnal test and had a 1-hour hot soak (See
Appendix A). No distinction was made between hot and cold soak for this derivation. We
assumed that during any soak, the only factor driving change in the fuel tank temperature was
the difference between the tank temperature and the ambient temperature.
This equation only applies during parked conditions, which include the following time
intervals:
• From the start of the day (midnight) until the first trip (keyOnTime)
• From a keyOffTime until the next keyOnTime
• From the final keyOffTime time until the end of the day
Mathematical steps:
1. At time to = 0 or keyOffTime (start of soak), TTank = TThis value will either be the
ambient temperature at the start of the day, or the fuel tank temperature at the end of a
trip.
Equation 3-1
10
-------
2. Then, for all t >0 and keyOffTime, the next tank temperature is calculated by integrating
numerically13 over the function for temperature change, using Equation 3-2.
(TTank)n+l — Txank + k(Tair Txank)At Equation 3-2
where:
TTank = Tank temperature
Tair = Ambient air temperature
t = Time
-l
k = Temperature constant (1.4 hr )
Figure 3 demonstrates the Euler approximation for calculating the tank temperature based on
ambient temperature.
60 1
6 10 14 18 22 2 6
Hourof Day
Figure 3: Example Day Modeled with Euler Method
—•— Calculated Tank
Temperature
Ambient Temperature
b Numerical integration is used to perform this step using the Euler method, one of the simplest methods of integration. The
smaller the time step At, the more accurate the solution. MOVES uses a At of 15 minutes, which is accurate enough for our
modeling purposes without causing tremendous strain on computing resources.
11
-------
3.1.2 Fuel Tank Temperature while Running
Vehicle trips are short compared to the length of the day. Therefore, we assume a linear
temperature increase during a trip to improve model performance with minimal compromise to
accuracy.
In this algorithm, to determine AT tank, the increase in tank temperature for an arbitrary trip
length with arbitrary starting tank temperature, we must first find ATtank95, the average
increase in tank temperature during a standard 4300 second, 95°F running loss test. The
algorithm models the increase in fuel tank temperature using the tank temperature at
KeyOnTime, the amount of running time and the vehicle type and technology. Newer
technologies reduce the heat transferred to the fuel tank. The MOVES ATia/ik95
temperatures are as follows:
• If the vehicle is pre-enhanced (model year pre-1996), vehicle type affects ATtank95 '- 13
LDV ATtank95 = 35°F
LDT ATtank95 = 29°F
• If the vehicle is evap-enhanced (model year 1996+):
ATtank95 =24°F
These values are used to calculate the change in temperature for a 4300 second test
(ATtank4300) for arbitrary starting fuel tank temperatures using Equation 3-3
^l ]'ank43()() = 0.352(95 7']auk,KeyON ) + kTTcmk95 Equation 3-3
The parameters in Equation 3-3 are derived from regression analyses of light-duty vehicles
driving the running loss drive cycle with varied starting temperatures 25. The lower the initial
tank temperature, the larger the increase over a given drive cycle.
The average ratio of fuel temperature increase to initial fuel temperature is -0.352. This gives
us the increase in tank temperature, so we can create a linear function that models fuel tank
temperature for each trip.
TTank ~ 4300/36OO ^ _ ^Tank,KeyON Equation 3-4
12
-------
Where:
TTank = Tank temperature
t = Time
tkeyON = Time of engine start
The 4300/3600 in the Equation 3-4 denominator converts seconds to hours (4300 seconds in
the running loss certification test), maintaining temporal consistency in the algorithm. The
resultant tank temperatures for an example temperature cycle are illustrated in Figure 4.
Running operation is shown as a red line, and hot soak operation is shown as a blue line
Figure 4: Modeled Vehicle Tank Temperature During a Day of Operation
120
110
GT 100
&D
01
3 90
— Ambient Temperature
Hot Soak Operation
Running Operation
70
60-
-1
9 14
Hour of Day
19
24
Assumptions:
• The first trip is assumed to start halfway into the hour stated in the first trip's
HourDaylD.
• The effect of a change in ambient temperature during a trip is negligible compared to
the temperature change caused by operation.
• The KeyOn tank temperature is known from calculation of tank temperature from the
previous soak.
13
-------
3.2 Permeation
Permeation emissions are specific hydrocarbon compounds that escape through micro-pores in
pipes, fittings, fuel tanks, and other vehicle components (typically made of plastic or rubber).
They differ from leaks in that they occur on the molecular level and do not represent a
mechanical/material failure in a specific location. In MOVES, base permeation rates are
estimated and then adjusted to the modelled tank temperatures and fuel properties.
The base rates in MOVES represent rates for gasoline with no ethanol. As detailed in Section
3.2.2, while keeping other fuel properties constant, the presence of ethanol increases
permeation emissions approximately two-fold.
For model years before 1999, permeation base rates are developed using the mg/hour emission
data from the last six hours of a 72-96-72T diurnal test (also known as cold soak/resting loss)
The diurnal tests were measured on the federal cycle (72F-96T) for the CRC E-9 and E-41
programs 12'14'15. Together, these two programs tested a total of 151 vehicles with model
years ranging from 1971 to 1997. The final six hours of the diurnal are an appropriate
surrogate for a permeation-only test because the emission rate, ambient temperature, and fuel
temperature are relatively stable or constant. Permeation should be the only evaporative
process occurring. The rates were developed for distinct model year and age groups. Model
years 1996-1998 are represented individually to reflect the phase-in of Enhanced Evaporative
Emissions Standards (20 percent Enhanced Evap in 1996, 40 percent in 1997, 90 percent in
1998).
For the 1999 - 2003 model years (full implementation of the Enhanced Evaporative Emissions
Standards) and the 2004 - 2015 model years (Tier 2 Evaporative Emissions Standards) the E-
77-2 Static test data was analyzed using the 86° F tests for non-leaking vehicles, corrected to
72°F. We used the Static test for these vehicles because the last six hours of Diurnal test was a
surrogate to represent permeation rates in the older vehicle data, but the static test is intended
to measure permeation and only permeation. The values in MOVES are based on 54 tests in
the 1999-2003 model year group taken from EPA's Compliance Testing Program. Later data
from the E-77 programs served to validate the Tier 2 permeation base rates already used in
MOVES. The E-77-2 data points expanded the range of the age groups, but the data was not
sufficient to differentiate the estimates for the age and model year groupings, thus, we kept the
0.0102 rate for all ages of both 1999-2003 and 2004+. Table 6 summarizes the data, analysis
and resulting emission rates for model years 1999-and-later for the number of tests within each
age group to understand how these decisions were made. There were three vehicles in the 0-3
age group for model year 2004+ that had an average of 0.003 but ten vehicles in the next age
group had essentially the same average emissions rate of 0.01. Assuming there is deterioration
for these vehicles, averaging all ages into one group made sense to characterize the fleet
without data to support finer age groups.
The Tier 3 evap standards apply starting in model year 2017, and phase in over model years
2016-2022, with early allowances. The Tier 3 permeation standard reflects a 40 percent
reduction from the previous standard and the introduction of 10 percent ethanol to the
certification fuel discussed earlier.
14
-------
Permeation base rates are presented in Table 7. These rates are recorded in the MOVES
emissionRateByAge table.
15
-------
Table 6: Permeation Analysis for 1999 - 2003 and 2004 and later Model Years
Previous Data
(Compliance data)
Decision for
MOVES2014
Model Year
Age
# of tests
MOVES2010
# of E-77 tests
E-77 avg rate
weighted avg of ALL
Enhanced Evap
1999-2003
0-3
52
0.0102
0.0102
0.0102
4-5
0
0.0102
0.0102
0.0186
0.0140
0.0102
0.0102
0.0102
6-7
2
0.0102
7
0.021
8-9
0
0.0102
2
0.014
10-14
0
0.0102
15-19
0
0.0102
20+
0
0.0102
Model Year
Age
# of tests
MOVES2010
# of E-77 tests
E-77 avg rate
weighted avg of ALL
Tier 2 / LEV II
2004-2015
0-3
0
0.0102
3
0.003
0.003
0.01
0
0
0
0
0
0.0102
4-5
0
0.0102
10
0.01
6-7
0
0.0102
8-9
0
0.0102
10-14
0
0.0102
15-19
0
0.0102
20+
0
0.0102
Model Year
Age
# of tests
MOVES2010
weighted avg of ALL
40% of Tier 3
2016-2017
0-3
0
0.0102
0
0
0
0
0
0
0
0.007
4-5
0
0.0102
6-7
0
0.0102
8-9
0
0.0102
10-14
0
0.0102
15-19
0
0.0102
20+
0
0.0102
Model Year
Age
# of tests
MOVES2010
weighted avg of ALL
60% of Tier 3
2018-2019
0-3
0
0.0102
0
0
0
0
0
0
0
0.006
4-5
0
0.0102
6-7
0
0.0102
8-9
0
0.0102
10-14
0
0.0102
15-19
0
0.0102
20+
0
0.0102
Model Year
Age
# of tests
MOVES2010
weighted avg of ALL
80% of Tier 3
2020 - 2021
0-3
0
0.0102
0
0
0
0
0
0
0
0.004
4-5
0
0.0102
6-7
0
0.0102
8-9
0
0.0102
10-14
0
0.0102
15-19
0
0.0102
20+
0
0.0102
Model Year
Age
# of tests
MOVES2010
# of E-77 tests
E-77 avg rate
weighted avg of ALL
100% of Tier 3
2022 and later
0-3
0
0.0102
1
0.005
0.005
0.002
0
0
0
0
0
0.003
4-5
0
0.0102
1
0.002
6-7
0
0.0102
8-9
0
0.0102
10-14
0
0.0102
15-19
0
0.0102
20+
0
0.0102
16
-------
Table 7: MOVES Base Permeation Rates at 72°F
First Model Year
Last Model Year
Age Group
Base Rate
1960
1970
20+
0.311
1971
1977
10-14
0.192
15-19
0.229
20+
0.311
1978
1995
0-5
0.0554
6-9
0.0913
10-14
0.124
15-19
0.148
20+
0.201
1996
1996
0-5
0.0464
6-9
0.0751
8-9
0.0751
15-19
0.12
20+
0.163
1997
1997
0-5
0.0373
6-9
0.0589
10-14
0.0784
15-19
0.0929
20+
0.125
1998
1998
0-5
0.0147
6-9
0.0183
10-14
0.0216
15-19
0.024
20+
0.0293
1999
2015
all ages
0.0102
2016
2017
all ages
0.0072
2018
2019
all ages
0.0056
2020
2021
all ages
0.0041
2022
2060
all ages
0.0026
3.2.1 Temperature Adjustment
The E-65 permeation study found that permeation rates, on average, double for every 18°F
increase in temperature. 16 This study tested 10 vehicle fuel systems (the vehicle body was cut
away from the fuel system, which remained intact on a frame) at 85°F and 105°F. The vehicles
ranged in model year from 1978-2001. In MOVES, the base permeation rates are calculated at
72°F, the same temperature as the certification test. Equation 3-5 is derived from this study and
used to adjust the base permeation rate to account for the tank temperature described in
17
-------
Section 3.1.
P = p* e00385(^TTank Tbase'y Equation 3-5
adj base
Where:
Pbase Base Permeation Rate
Tiank = Tank Temperature
Tbase = Base Temperature for a given cycle (e.g., 72°F for a federal diurnal test)
3.2.2 Fuel Adjustment
Ethanol affects evaporative emissions from gasoline vehicles due to the increased permeation
of specific hydrocarbon compounds through tanks and hoses. Modeling permeation emissions
separately from vapor venting emissions allows us to apply these effects only on permeation
where the complex chemistry of ethanol-gasoline blends increases permeation through the
tanks and hoses of the fuel system.
Permeation fuel effects were developed from the CRC E-6516 and E-65-317 programs, which
measured evaporative emissions from ten fuel systems that were removed from the vehicles
and filled withEO, E5.7, and E10 fuels. This method assures that the emissions measured are
purely from permeation (assuming the systems were not leaking). Additional data was
provided from the CRC E-77-27 and E-77-2b8 programs, which measured evaporative
emissions from sixteen intact vehicles. For this analysis, vehicles certified to enhanced-
evaporative and Tier 2 standards were analyzed separately from vehicles certified to earlier
standards. Enhanced evaporative standards were phased in from 1996-1999 and imposed a 2.0
gram standard over a 24-hour diurnal test. Standards previously in effect applied a 2.0 gram
standard to a 1-hour simulated diurnal.
The ethanol effect is estimated with a mixed model shown in Table 8. The evaporative
certification level, ethanol content, and RVP were modeled as fixed effects and the vehicle was
modeled as a random effect. The natural logarithm of the emission rates over the 65-105-65°F
diurnal cycle provided a normally distributed dataset to the model. The dataset was not large
enough to distinguish the three ethanol levels within each evaporative certification bin.
Therefore, E5.7 and E10 test results were combined into a single bin of ethanol-containing
fuel which had a significant effect compared to E0 fuel.
18
-------
Table 8: Mixed Model for Ethanol Effects on Permeation
Mixed Model for ETOH Effects on Permeation
logHC-1 ~ Standard Group + ETOHYN * Standard Group
Standard
Value
Std. Error
DF
t-value
p-value
3.0279465
0.1887947
59
16.038305
0
Tier 2
-0.761386
0.3095233
15
-2.459867
0.0265
Zero Evap
-1.050358
0.6297706
15
-2.938146
0.0102
The percent difference between the ethanol rate and the EO rate is used in MOVES as the fuel
adjustment. Due to the enhanced-evaporative certification standards phase in from 1996-1999
(20/40/90/100 percent), the two fuel adjustments must also be phased in for those model years.
Table 9 lists the fuel adjustments used for E5 through E85 for the model year groupings used
in MOVES. These values are recorded in the MOVES hcPermeationCoeff table as
multiplicative factors.
Table 9: Ethanol effect for Permeation Emissions
Model Years
Percent increase due to
ethanol (5-85%)
1995 andearlier
65.9
1996
75.5
1997-2000
107.3
2001 and later
113.8
3.3 Tank Vapor Venting
Vapor generated in the tank can escape to the atmosphere during a process labeled "Evap
Vapor Venting" or "Tank Vapor Venting" (TVV). Hydrocarbons emitted by this process
originate from a variety of sources. As tank temperature rises and vapor is generated within the
tank, the vapors are forced out of the tank from increased pressure. Fully sealed gas tanks are
rare as they must be constructed with metal to prevent bloating. Using metal as a tank
material can be expensive, heavy, and difficult to shape for tightly packed modern vehicles.
Instead, most vehicles are equipped with an activated charcoal canister to adsorb the vapors as
they are generated. Later, the vapors are consumed as they purge to the engine (through the
intake manifold) during vehicle operation. The canister is open (or vented) to the atmosphere
to prevent pressure from building within the fuel system. Consequently, if the engine is not
operated for several days, fuel vapors can diffuse through the charcoal or even freely pass
through a completely saturated canister. Tampering, mal-maintenance, and system failure can
result in excess evaporative emissions. Inspection and maintenance (I/M) programs can also
19
-------
influence how leaks and other problems are controlled over the life of a vehicle. MOVES
models vapor venting for vehicles using gasoline and E85 fuels.
In MOVES, vapor venting default values are recorded in several tables. The
emissionratebyage table stores the vapor venting rates (process 12) for hot soak (opmode 150)
and running (opmode 300). Cold soak (opmode 151) emission equations and parameters are
recorded in cumTvvCoeffs as explained below.
Integral to the understanding of Tank Vapor Venting (TVV) is the calculation of Tank Vapor
Generated (TVG). Tank vapor generated depends on the rise in fuel tank temperature (F),
ethanol content (vol. percent), Reid vapor pressure (RVP, psi) and altitude. Calculations in
MOVES use the Wade-Reddy equation for vapor generation (Equation 3-6).26
TYQ — AeB*RVP (eCTx — eCTl ) Equation 3-6
Where:
T1 = Initial temperature
Ta = Temperature at timex
In Equation 3-6 coefficients A, B, C vary by altitude and fuel ethanol content. These
coefficients, stored in the database table TankVaporGenCoeffs, are shown in Table 10.
Table 10: TVG Constants for Equation 3-6
E0 Gasoline
E10 Gasoline
Constant
Sea Level
Denver alt.
Sea Level
Denver alt.
A
0.00817
0.00518
0.00875
0.00665
B
0.2357
0.2649
0.2056
0.2228
C
0.0409
0.0461
0.0430
0.0474
The vapor venting emission process occurs during all three operation modes: running, hot soak
and cold soak. While running, vapors are generated as the fuel system is warming and active.
During hot soak, vapor generation is caused by latent heat transfer due to fuel recirculation and
other convective processes. Cold soak vapor generation is concurrent with ambient
temperature increases. MOVES modeling of fuel system warming is detailed in Section 3.1.
20
-------
3.3.1 Altitude
Evaporative vapor generation is affected by the lower ambient pressure at high altitudes.
MOVES accounts for this effect during the calculation of tank vapor generated. This process
relies on the values in Table 10 for high altitude (Denver, CO) and a low altitude (Sea Level).
MOVES applies linear interpolation/extrapolation based on the barometric pressure in each
county.
3.3.2 Cold Soak
Cold soak vapor emissions occur while a vehicle is not operating, and the engine and fuel
system have cooled to ambient temperature. Emissions occurring under these conditions are
also referred to as diurnal emissions.
Emission quantification started with the measurement of emissions based solely on standard
regulatory test cycles. As the emissions levels over the test cycles became more controlled,
concern grew about "off-cycle emissions"- real-world emissions that occur outside of the test
procedure constraints on ambient temperatures, fuel RVP and soak time.
3.3.2.1 Multi-day Diurnals
As a vehicle sits through multiple diurnal cycles, the carbon canister accumulates vapor every
day. It can only adsorb vapor until it reaches its capacity; then it begins to vent to the
atmosphere. A canister with degraded/damaged carbon may have reduced capacity, and
eventually every canister will vent to the atmosphere once it reaches saturation. During cooling
hours, a canister backpurges to the fuel tank and regains some capacity. Then, during the
subsequent warming period the canister is re-filled with vapor and any vapor generated
beyond capacity will escape to the atmosphere.
21
-------
Figure 5: Multiday Vapor Accumulation in Charcoal Canister
Vapor
generated
1
Vapor
after
backpurge
Day 2
breakthrough
Empty warming cooling
Canister Day 1
warming cooling
Day 2
Day 3
breakthrough
warming cooling
Day 3
Figure 5 illustrates the dynamic behavior of vapor within a charcoal canister over three days of
continuous cold soaking. During the first day, vapor accumulates within the canister but does
not exceed the canister capacity. During the cooling period of day L we observe backpurge
when some of the fuel vapors that were previously adsorbed to the charcoal flow back into the
cooling tank. The fresh air is drawn in through the canister vent while the vapor condenses in
the tank during the cooling portion of the cycle. During warming on day 2, we see generated
fuel vapors that exceed the canister capacity (though some canisters may be constructed to
hold more than 2 days of vapor). These emissions are lost to the atmosphere, and only what
remains in the canister can be backpurged during the subsequent cooling cycle. In day 3, more
vapor is generated and consequently lost to the atmosphere. Any additional days without
engine purge during normal driving (i.e. inactivity) will exhibit the same behavior as day 3. It
should be mentioned that plug-in hybrid electric vehicles that are mainly driven on short
(electric only) trips, may also exhibit similar breakthrough over time. However, modeling of
these vehicles is currently beyond the scope of MOVES since the penetration rates of these
technologies are low, and we are not aware of any multi-day diurnal data collected on PHEVs.
3.3.2.2 Cold Soak Emissions for the Average Fleet
Modeling a fleet of vehicles involves a diverse population of canisters with differing capacities.
A given amount of vapor will be fully contained by some vehicles but exceed the canister
capacity in others. Figure 6 illustrates the approach for calculating the tank vapor vented
(TVV) as a function of the tank vapor generated (TVG).
22
-------
Figure 6: Vapor Vented Curve
Several factors accommodate this modeling approach. The following variables, explained in
more detail below, are included in the MOVES default database in the 'cumTvvCoeffs'
table:
• BackPurge Factor
• Average Canister Capacity
• Tank Size
• Tank Fill Fraction
• Leak Fraction (Prevalence)
• Leak Fraction (Prevalence)IM
• TVV Equation
• Leak Equation
3.3.2.3 B ackPurge F actor
The backpurge factor is the percent of hydrocarbon vapor that is desorbed from a vehicle's
canister during cooling hours when pressure decreases within the tank, drawing ambient air in
through the canister vent. In the real-world, this process occurs nightly as temperatures cool
and restore some canister capacity. In the Multiday Diurnal Study 18, test vehicles soaked for
14 consecutive 72F-96F diurnals (the Federal Test Procedure temperature cycle). During this
time, the vehicle canister mass was measured continuously. During the cooling period, the
measured mass of the vehicle canisters decreased. This cyclical effect can be observed in
Figure 7.
An average value of 23.8 percent backpurge was developed from these results and is used in
the MOVES model. For example, a vehicle canister with 100 grams of hydrocarbons will
backpurge 23.8 grams and begin the next day with 76.2 grams. A more complex model for
backpurge was considered (similar to tank vapor generation), but it would require significant
computational resources and potentially slow model performance considerably. As diurnal
23
-------
temperatures are relatively symmetrical, detailed modeling of tank vapor generation has
already provided a high level of precision, justifying a simpler model here.
Figure 7: Vehicle Canister Mass, 14-day Diurnal Test
Eh
d
15€
C/5 C
If)
an
^ 3
<3 D
!-h
C/5 3
10g
0 Oh
-5 s
43 H
;>
5€
Legend
lORVPfuel
9 RVP fuel
SHED Temperature
5000 10000
Test Time (Minutes)
15000
20000
3.3.2.4 A verage Canister Capacity
The canister capacity reflects how much vapor generated in the tank can be contained by the
canister before breaking through. Note that while average canister capacity is stored in the
MOVES database, it is not actually used in the MOVES model. Instead, it was used in the
DELTA pre-processor 27
To calculate a sales-weighted average canister size, we used sales data28 and EPA evaporation
certification data.22 Certification data includes the evaporative family code which contains the
Butane Working Capacity (BWC) of the canister; it is found in digits 7, 8 and 9 for enhanced
evap vehicles, and in digits 5, 6 and 7 for pre-enhanced vehicles. The BWC represents the
ability of a canister to capture butane vapor, rather than gasoline vapor, so it must be adjusted
by a factor of 0.9229.
Evaporative control was introduced in 1971, so canisters are not modelled for pre-1971
vehicles. For model years beyond 2010, the 2010 average canister capacity was used. The
calculated average canister capacities for cars and trucks combined are listed in Table 11.
A peak in average canister size at model year 2005 corresponds to greater sales of cars with
larger fuel tanks. Motorcycles are modeled without canisters, because they were not being
24
-------
used at the time of the analysis. We hope to review and revise in future versions of the model.
Heavy-duty gasoline vehicles are modeled with the same canister capacity as light-duty
vehicles.
Table 11: Average Canister Capacity by Model Year
Model Year Group
Average Canister Capacity (grams)
1960-1970
0
1971-1977
64.7
1978-1995
72.8
1996
78.7
1997
83
1998
115.4
1999-2003
122.9
2004
145
2005
150.7
2006
145.3
2007
142.9
2008
138.6
2009
136.2
2010+
137.5
3.3.2.5 Tank Size
The average tank size for a given model year is an important facet of the vapor generation
calculation because a larger tank will have more space in which vapor can accumulate. Both
sales data28 and tank size information 30 were required to calculate a sales-weighted average
tank size for model years 1990-2010. For this analysis, car and truck sales, and tank sizes were
combined. For vehicles with multiple styles (i.e. different cab sizes on pick-up trucks) with
different tank sizes, the average available tank size was used as sales information was
unavailable by style. Data sources only span from 1990-2010, so past and future values were
projected. Vehicles in the 1990-2010 range have tanks with an average capacity of 1.25 times
greater than a calculated 300-mile range, so this ratio was applied using fuel economy data
going back to model year 1975. 31. Pre-1975 vehicles use the 1975 fuel tank size. For future
vehicles, tank size is assumed to stay constant from 2010 on. It is also possible that
manufacturers will maintain range with a shrinking fuel tank. In future versions of MOVES,
we will reexamine this assumption. The calculated sales-weighted tank sizes are in Table 12.
25
-------
Table 12: Sales-Weighted Average Fuel Tank Size
Model Year Group
Tank Size (gal)
1960-1970
28
1971-1977
27.3
1978-1995
18.6
1996-1997
19.1
1998
19.5
1999-2003
19.9
2004
20.5
2005
20.3
2006
20
2007
19.7
2008
19
2009-2030
19.1
HD Vehicles
38
Motorcycles
3
3.3.2.6 Tank Fill Fraction
The tank fill fraction is an important input used in calculating tank vapor generation. The more
vapor space above the liquid fuel, the more capacity there is for vapors to accumulate. The
average tank fill fraction used in the model is 40 percent fill. This is a typical fill level for
certification procedures and many of the test programs from which our data originates. It is
also a figure supported by existing research on tank filling behavior by consumers 32.
3.3.2.7 Vapor Leak Prevalence
In order to accurately quantify emissions from leaking vehicles, one must not only estimate
emission rates from leaks of various sizes, but also prevalence of leaks in the fleet. This
corresponds to an emissions rate and its corresponding activity. Our estimates of leak
prevalence are informed by the analysis of a field study which took place at the Ken Caryl IM
Station in Denver, CO during the summer of 2009 u. In this study, a remote sensing device
(RSD) was used to recruit high emitting vehicles which were then tested in a Portable Sealed
Housing for Evaporative Detection (PSHED). The vehicle's hydrocarbon emissions were
measured over 15 minutes during hot-soak conditions, and vehicles were inspected to identify
the cause/source of the leaks when possible. The set of hot-soak measurement from individual
vehicles, with inverse-probability sampling weights and solicitation response weights applied
to all vehicles, allows the prevalence of leaks in the fleet to be estimated.
We have defined a vapor leaker as any vehicle that would fail the enhanced evaporative
standard of 2 grams. The standard sums the emissions from the worst day of a 3-day diurnal
test and the hot soak. To develop a surrogate standard for a 15-minute hot soak test, we used
knowledge of certification testing to attribute 0.4 grams (g) of the 2 g standard to the hot soak
26
-------
portion, and 76 percent of 0.4 g to the first 15 minutes of the hour-long hot soak test. This
approach suggests that 0.3 g can be taken as a surrogate standard for a 15-minute hot soak.
Figure 8: Prevalence of Vapor Leaks Above a Given Threshold in the 2009 Ken Caryl
Fleet
Vapor Leak
Size (g/
15min)
Table 13 (plotted in Figure 8) displays leak prevalence at various emission thresholds for what
constitutes a "leak". Observing the difference between any two points determines how many
vehicles fall into a particular range. Looking at Table 13, in model year group 1981-1995, 2.6
percent of vehicles are leaking at more than 20 grams and 4.2 percent of vehicles are leaking at
more than 10 grams. Subtracting these two values yields that 1.6 percent of vehicles in the
model year group have a leak between 10 and 20 grams.
The data only contain prevalence rates for PSHED measurements as low as 1.0g/15min.
Failure rates are extrapolated to 0.3g/15min. Using aggregate data from the Ken Caryl station,
it is found that 0.3g/15min PSHED measurements are 50 percent more prevalent than
1.0g/15min PSHED measurements.
27
-------
Table 13: Prevalence of Leaks Above a Given Threshold (g/15min)
Leak Threshold (g/15 min)
Denver
100
50
20
10
5
2
1
0.3
Sea Level (MOVES)
70.9
35.5
14.1
7
3.6
1.4
.7
.2
Prevalence by MY Group
1961 - 1970
0
0
0.53
0.53
0.68
0.68
1
1
1971 - 1980
0
0
0
0.3
0.85
1
1
1
1981 - 1995
0.004
0.004
0.026
0.042
0.083
0.22
0.26
0.39
1996- 2003
0
0
0
0.02
0.021
0.029
0.033
0.064
2004- 2010
0
0
0
0
0
0
0
0
Because the data used to estimate leak prevalence was collected in Denver, Colorado at an
altitude of 5,280 feet above sea level, measurements must be adjusted to sea level. At sea
level, the amount of vapor generated will be less due to higher atmospheric pressure. To
determine the appropriate correction factor, we performed the Wade-Reddy calculation and
found that under identical conditions, the higher altitude will generate 41 percent more vapor.
Colorado is a strategic location to perform a leak quantification program because a given
vapor leak will produce higher levels of emissions at a higher altitude, therefore making it
easier to detect. Each of the leak magnitude bins have been corrected for altitude by this factor.
For example, the prevalence of leaks at lg-2g levels in Denver will be the same prevalence of
leaks at ,71g-1.42g levels at sea level.
Because this was a cross-sectional study, we must populate the model for many model year and
age group combinations that were not measured. A set of linear regressions was used to model
vapor leak prevalence for ages and model years where data is not available. We divided model
year groups in years when new technologies or standards were introduced. Modeling was
based on the assumption that newer cars will have lower leak prevalence than older cars due to
the advancing technology and use of more durable materials. Therefore, data from the 1996-
2003 model year group was used as a surrogate for new vehicles in the 1971-1980 and 1981-
1995 model year groups. However, because vapor leaks also occur due to tampering and mal-
maintenance, deterioration is not the only factor involved in occurrence of vapor leaks. The
regressions from the older model year show more rapid vehicle deterioration rates than newer
model years.
Figure 9 shows the vapor leak prevalence as the percent of the vehicle fleet with a leak larger
than 0.3g/15min. For model years 1996 and later, the estimate for leak prevalence at ages 0-3
was developed with I/M data from five states. The analysis revealed that 1-2 percent of
vehicles consistently arrived at I/M stations with an evap Diagnostic Trouble Code (DTC) set.
The vast majority of the DTCs set specifically indicated a vapor leak detected. The green
diamonds in the 1971-1980 and 1981-1995 model year groups are an assumption made based on
the 1996-2003 data to describe these vehicles' leak rates when they were new. The slope of the
28
-------
2004-2010 prevalence rates was developed by applying the 5-10 year-old 1996-2003 data point
to the 10-15 year old 2004-2010 point.
Figure 9: Non-IM Vapor Leak Prevalence, Extrapolated from data
Tier 3 and LEY III Leak Prevalence
To model the leak prevalence rates of LEY III and Tier 3 vehicles, the effectiveness of
improved OBD systems and the efficacy of vehicle leak testing were quantified. In the
above mentioned field study performed in Colorado, it was found that 70 percent of
evaporative leaks were due to deterioration of the evaporative system (e.g. corroded fuel lines,
filler neck, cracked hoses etc.) that could be improved with new design and material
considerations. The remaining 30 percent of evaporative leaks were beyond manufacturer
control, (e.g. Improper maintenance, tampering, missing gas caps, etc). See Table 14.
OBD effectiveness and OBD readiness are also important factors in the detection and repair of
leaks after they occur. OBD effectiveness refers to the ability of diagnostic systems to identify
leaks within the fuel system and alert the driver by illuminating a warning light. OBD
readiness refers to the time during which vehicle diagnostics are actively assessing the
integrity of the vehicle fuel system.
Our reference case assumed 40 percent OBD effectiveness and 95 percent OBD readiness.
These numbers were based on an assessment of vehicles with OBD-detectable leaks and
whether or not the leak was identified by the vehicle and the driver alerted via a check engine
light. 33
29
-------
We estimated the implementation of LEV III would immediately reduce the 70 percent of
deterioration-caused leaks by 33 percent simply due to the lower emissions standard.
Longitudinally, we see reductions in leak prevalence associated with lower emissions
standards. We also estimated that due to improved vehicle diagnostic systems, 80 percent of
detectable leaks will be discovered and reported by the vehicle. In addition, we assumed that
with the increased rigor of requirements the readiness will increase to 99 percent.
We estimated that the implementation of Tier 3 would immediately reduce the 70 percent of
deterioration-caused leaks by 66 percent due to the additional benefit of the Tier 3 leak
standard. As in LEV III estimates, we also estimated that 80 percent of detectable leaks will
be discovered and reported by Tier 3 vehicles, as well as an increase of 99 percent readiness.
These estimates result in an overall reduction of leak frequency of 26 percent for the LEV III
program and 49 percent for the Tier 3 program
Table 14: Summary of Tier 3 and LEV Assumptions
Base Inputs
# of Leaks >0.020"
100
% Mai-Repair
30%
% Durability
70%
Tier 2 Case
OBD Ready %
95%
OBD Effectiveness
40%
LEV III Control Case
% of "durability" leaks prevented
33%
OBD Ready %
99%
OBD Effectiveness
80%
Tier 3 Control Case
% of "durability" leaks prevented
66%
OBD Ready %
99%
OBD Effectiveness
80%
30
-------
Figure 10: LEV III, Tier 3 Leak Prevalence Estimates
u
JZ
01
>
u-
o
c
o
• hi
£
(0
100
TO
-------
Figure 11: SHED Leak Emissions for one Severity Bin
14
Vapor generated in the tank (TVG) is calculated using the Wade-Reddy equation (Equation
3-6), thus requiring fuel RVP, fuel ethanol content, and temperature data. Two datasets
containing this information were used in developing leak emission rates. The E-77 suite of
programs 6"7"8-9 measured high-emitting vehicles, with known fuel properties and artificially
implanted leaks, on the California (65°F-105°F) diurnal cycle. In another effort, the Colorado
Department of Public Health and Environment (CDPHE) carried out a repair effectiveness
program during the summer of 2010 in collaboration with the Regional Air Quality Council
(RAQC). This program 34 measured 16 vehicles with identified leaks. A 6-hour test was
performed with a temperature increase of 72°F- 96°F. This effort was less resource-intensive
than the full diurnal procedure and still provides the necessary information to calculate TVG.
The SHED measurements of Tank Vapor Vented (TVV) and calculated TVG form the basis
for a linear regression of TVV vs. TVG for each vehicle. The resulting slope represents the
mass of vapor vented per mass of vapor generated. The average of the regressions becomes the
leak rate for that severity bin. This approach can be observed in Figure 11. This test procedure
could not distinguish permeation and leak vapor emissions. However, permeation for these
vehicles is assumed to be negligible during the 6-hour test given the severity of the leak
emissions. In the E-77 program, TVV emissions were collected in a canister external to the
SHED. The external canister was connected to the vent on the vehicle canister. No permeation
was included in the measurement.
Because the emissions measured were highly variable, spanning several orders of magnitude,
the emissions data for leaking vehicles was binned by magnitude. Accordingly, both emission
rates and prevalence were calculated within these bins. As the leak prevalence estimates were
measured at high altitude in Denver, it is essential to develop adjustments to apply the binning
process at lower altitudes, such as sea level. Application of the Wade-Reddy equation
(Equation 3-6) suggests that an E10 fuel in Denver generates 1.41 times as much vapor as at
32
-------
sea level. For example, a vapor leak at 0.3g/15min in Denver would have an equivalent rate of
0.21g/15min at sea level. The bins used to categorize leak severity as well as the average leak
emission rate for that bin are listed in Table 15.
Table 15: Leak Emission Rates by Bin
Denver bins
(g/15min)
Sea Level bins
(g/15min)
Average Ratio of Grams vented /
Grams generated
0.3 - 2
0.2 - 1.4
0.12
2 - 5
1.4- 3.6
0.27
5 - 10
3.6- 7.1
0.65
>10
>7.1
1.33
Each data point was binned by its hot soak measurement from the E-77 programs or PSHED.
(Portable SHED) measurement from the Denver program. The PSHED tests are 15-minute hot
soak measurements.
Figure 12 illustrates the leak emission rates for each leak severity bin. The average emission
rate for vehicles with 15-min hot soak measurements greater than lOg exceeds 1. It is possible
to measure more fuel vapor in the SHED than is calculated with Equation 3-6 It is known that
the equation is less reliable at higher temperatures. Also, complicated factors such as fuel
sloshing and tank geometry can influence vapor generation beyond the estimation capabilities
of the Wade-Reddy equation.
Figure 12: Leak Emission Rates by Leak Severity Bin
5 10 15 20 25 30 35
Tank VaporGenerated
33
-------
An average emission coefficient for an "average severity" leak (Table 16) was computed based
on the prevalence information in Section 3.3.2.7. These values are recorded in the leakEquation
field in the cumTVVCoeffs table of the MOVES default database.
Table 16 Leak Emission Coefficients
Regulatory
Class
Model Years
Age Group
Leak Coefficient
All Except
1960-1970
20+
0.952
Motorcycles
1971-1977
10-14
0.782
15-19
0.79
20+
0.796
1978-1995
0-3
0.524
4-5
0.408
6-7
0.388
8-9
0.376
10-14
0.365
15-19
0.357
20+
0.351
1996 and Later
All
0.524
Motorcycles
All
All
0.814
Then the estimated vapor venting from leaking for each of soaking day is a multiplicative
function of the first day vapor generation (TVGDaily) as shown in Equation 3-7.
LeakEmissions = LeakCoefficient * TVGDaily Equation
TVGDaily is the same for each soakDaylD, and is also the same as Xn[soakDayID=l] (see
below), such that TVGDaily = tvgSumlH.
3.3.2.9 Tank Vapor Venting Equation
For normally operating non-leaking vehicles, tank vapor vented (TVV) from the canister was
calculated. First MOVES estimates the quantity of vapor generated using Equation 3-6 in
34
-------
g/gal-headspace, based on the maximum temperature for the day, interpolated for the
appropriate altitude and gasoline ethanol level, and allocated to hour of the day.
The model then estimates the vapor vented for consecutive days of parked (soak) activity. The
algorithm accounts for average canister capacity (ACC) and backpurge factor. Daily backpurge
removes fuel vapors from the canister, increasing capacity to store vapor generated during
successive days. Vapor generated above the ACC is lost to the atmosphere, therefore backpurge
only applies to what remains in the canister. For each hour of continuous parking time, the
model calculates the vapor generated on soak day "n", represented as Xn in Equation 3-8 below.
Xn[soakDayID=l] = tvgSumlH
Equation
3-8
Xn[soakDayID=2] = ((1-
backPurgeFactor)*least(tvgSumI24,averageCanisterCapacity))
+tvgSumlH
Xn[soakDayID>2] = ((l-backPurgeFactor)*least(Xn[soakDayID-l] +
tvgSumH24,averageCanisterCapacity))+tvgSumlH
Xn is calculated for soak days 1 through 5. Beyond five days, the algorithm assumes that
breakthrough has occurred and that behavior over additional days has stabilized. The vapor
emissions are fleet averages by model year group. Vapor venting emissions rises on successive
days as more vehicles exceed their canister capacities and begin venting fuel vapors.
The vapor generated values are then summed for each hour. The vapor leak and the cold soak
emissions are totaled using the following equations for vehicle groups with and without relevant
Inspection and Maintenance (I/M) programs.
tankVaporVented = (l-leakFraction)*(TW equation) + leakFr action* (Leak
equation). q"aq10n
tankVaporVentedIM = (l-leakFractionIM)*(TW equation) +
leakFractionIM*(Leak equation).
Where:
leakFraction is the appropriate leak prevalence as described in Section 3.3.2.7
TW equation is a function of Xn as listed in the leakEquation field in the
cumTWCoeffs table of the MOVES default database.
Leak equation is the result of Equation 3-8, above.
leadFractionIMis the leak prevalence described in Section 3.4.2.
For more detail on the vapor venting calculations, see the MOVES code at
database\MultidayTankVaporVentingCalculator.sql and documentation of the DELTA pre-
processor.27
35
-------
3.3.2.10 Activity
In order to properly account for off-cycle emissions, MOVES must account for the different
emissions rates of short (several hours) and long (multiple day) soaks. For any modeled day,
there is a sub-population of vehicles exhibiting 1st, 2nd, 3rd, nth day diurnal emissions. The
fractional allocations for 1st, 2nd, 3rd, and nth day diurnals are calculated from the
sampleVehicleTrip and sampleVehicleDays tables in MOVES. SampleVehicleTrip assigns
numbers of first starts during each hour of the day. For the fraction of vehicles having soaked
since at least midnight, the first engine start ends the cold soak episode. SampleVehicleDay
contains the population of vehicles for each sourceTypelD.5 Combining information for both
tables, it is simple to calculate the fraction of vehicles having soaked since midnight at any
given hour. For example, at 1:00AM, some fraction of vehicles less than 100 percent have not
yet started. The fraction continuously decreases throughout the day as more and more vehicles
start. At 12:00AM, the fraction only represents vehicles that were not driven.
Once the fraction of vehicles soaking at a given hour has been calculated, it must be estimated
how many prior days each has been soaking. We classify vehicles as 1st day, 2nd day, 3rd day, 4th
day, or 5+ days. We assume that after the 5th day, vehicles will exhibit repeat emissions since
the evaporative canister will either have broken through or be in conditions that will never
cause breakthrough. An activity study performed by Georgia Technological University 35
suggest that 16 percent of vehicles drive less than 3,000 miles per year. The MOVES inputs
are based on the conservative estimate that 50 percent of these low-mileage vehicles, or 8
percent of all vehicles, have been soaking for more than 5 days on any given day.
The sampleVehicleSoakingDayBasis table establishes the fraction of vehicles soaking for 5+
days. It contains five values, one for each soak day. The value for SoakDaylD 1 is the
percentage of vehicles soaking at the final hour of day 1. The product of SoakDayID=l and
SoakDayID=2 is the percent of vehicles soaking at the final hour of day 2. The product of all
five values is the percent of vehicles soaking for five days or longer.
Figure 13 presents the fraction of soaking vehicles throughout the day. The majority of
vehicles were driven the previous day and are on their first day soaking. The fractions of
vehicles on 2nd through 4th day soaking are developed from the remainder of 1st day soaking
vehicles at hour 24. The fraction of vehicles soaking for 5 days or longer is 8 percent at hour
24. This method models bimodal vehicle usage, with most vehicles being driven almost daily
and the remaining vehicles being driven more intermittently.
36
-------
Figure 13. Passenger Car Soak Distribution on a Weekday
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour of Day
3.3.3 Hot Soak
Hot-soak vapor emissions begin immediately after a car ceases operation and continue until
the fuel tank reaches ambient temperature. In MOVES, the process of calculating hot-soak
vapor emissions is simpler than that for cold soak. Base rates exist for each model year and
age group and are expressed in units of grams per hour. They represent emissions at sea
level with RVP assumed at 9.0 psi. In developing the rates, leak and non-leak rates are
weighted together to form the base rate, similar to cold soak.
3.3.3.1 Hot Soak Data
Hot soak data comes from several programs with diverse testing procedures, vehicle model
years and technology, fuel parameters, and altitude. These programs include three summer
programs in Colorado and the E-77-2 programs in Arizona. See Table 17.
37
-------
Table 17: Hot Soak Evaporative Test Programs
Program
Location
Hot Soak Length
Fuel RVP
Altitude (ft)
No. Obs.
High Evap
Lipan IM station, CO
15 min
Fuel Supply
5130
100
High Evap
Ken Caryl IM station, CO
15 min
Fuel Supply
5130
175
High Evap
Denver IM station, CO
15 min, 1 hour
Fuel Supply
5130
100
E-77-2
Mesa, AZ
1 hour
7, 9, 10
1243
100
As explained below, the data collected in these programs was adjusted to the MOVES baseline
of a one-hour rate on 9.0 RVP fuel at sea level.
In addition, some tests were removed from our analysis. The vehicles in Colorado that
participated in the studies were recruited in-situ and therefore were subject to a wide range of
leak mechanisms. It was observed that some vehicles emitting more than 50 grams in 15
minutes in the PSHED had liquid leaks present. All vehicles with a calculated 15-minute
measurement greater than 50g/15min were removed from vapor leak analysis.
Furthermore, vehicles in the E-77 program were tested multiple times with different fuels,
whereas each vehicle in the Colorado population was tested once. In order to not over-
represent the E-77 vehicles in our sample, one measurement from each vehicle was selected
with preference given to the measurements on 9 RVP, E10 fuels (where available).
3.3.3.2 Test Duration Conversion
Every datum required a 15-minute mass and a one-hour mass because base rates in the
MOVES input table must be expressed in grams per hour; however, our method for
distinguishing leaks from non-leaks uses the 15-minute rate. Furthermore, if a measurement is
designated as from a leaking vehicle, the 15-minute measurement is used to project its rate of
occurrence in the fleet.
Because engines and fuel systems do not cool at a uniform rate, existing data was used to
develop this test-length conversion factor. In the E-77 suite, the cumulative time series data
for hot-soak tests on a minute-by-minute scale was readily available, enabling estimation of
vapor emissions over 15 minutes. Each set of vehicle data also contained a permeation rate.
The permeation rate was subtracted from the 15-minute hot soak measurement. The result is
the assumed vapor emissions during 15-minutes of hot soak. Similarly, hourly permeation
was subtracted from the 1-hour hot soak measurement. After compiling the 15-minute and 1-
hour values, the fraction of emissions occurring in the first 15 minutes can be calculated.
All of the Denver testing programs provided similar vehicle measurements to augment the E-
77 dataset. A subset of the vehicles were transported to a lab to receive a Hot Soak test.
Readings were taken at both 15 and 60 minutes.
Figure 14 illustrates the evaporative emissions occurring during a Hot Soak test. Vapor
emitted by permeation is assumed to accumulate at a linear rate while vapor emissions
attributed to the hot soak accumulate rapidly following engine shutoff but more slowly as the
engine cools.
38
-------
Figure 14: Cumulative Hot Soak and Permeation Emissions
O Hot Soak
Time (m)
Using the combined data from E-77 and Denver testing, we estimated the average fraction of
emissions in the first fifteen minutes following engine shutoff. At first, it was thought that this
fraction would vaiy among groups of vehicles certified to different evaporative standards.
However, analysis of test results by certification groups did not seem to yield notably different
results. Instead a single fraction developed from all available data was applied fleet-wide. It
was estimated that 54 percent of emissions from a one-hour hot soak occur in the first 15
minutes. Conversely, emissions from a 15-minute hot soak must be multiplied by 1.85 to
estimate a full hour's emissions.
3.3.3.3 Correction for RVP and Altitude
MOVES base emission rates are intended to represent emissions on 9.0 RVP gasoline at sea
level so, the hot soak test data must also be corrected to account for the RVP and altitude of
each test.
Emissions in the available datasets were measured at varying levels of RVP. Some programs
recorded RVP, while other data has no explicit RVP information. Our first step is to estimate
the RVP for all measurements that do not contain this information.
The majority of the data with unknown RVP was gathered in the summer months in locations
with available fuel survey data. The mean RVP for June through August 2010 in Denver was
8.40 RVP (standard deviation 0.20 RVP), and this value was assumed for all vehicles tested
from May through September. For non-summer months, RVP information was collected with a
39
-------
small subset of the vehicle measurements. In the case of a non-summer measurement without
RVP information, the mean of all non-summer months is assumed. The mean RVP for non-
summer vehicles is 10.67 (standard deviation 1.75 RVP). The testing at the Lipan station was
all performed in the summer, so the RVP of the Lipan dataset is assumed to be 8.4.
Associating an RVP value with every measurement enables calculation of corrections for
altitude. All vehicles were tested either in Colorado (Elev. = 5,130 ft) or Mesa, AZ (Elev. =
1,243 ft). Both locations are far enough above sea-level that it would be erroneous to assume
their emissions are representative of sea-level emissions. Our approach is to use the E10
values in Table 10 to solve the Wade-Reddy equation (Equation 3-6) for RVP and thus
calculate the equivalent RVP at sea level that would generate the same amount of emissions as
measured at higher altitude, as shown in Equation 3-10
RVPSeaLevel - (BsgaLgvj * ln (ASeaLevel * (
-------
TVG
measRVP
= AgQdiQygiG^Sealjevel*^^>meas (Q^SecLLevel*^! — Q^SeaLevel*^o^
TVGmoves = AhigheB^*9-°(ec^^ - ec^*T°)
/T~\f G cisR VP \
HotSoakMOVES HotSoakMeasured * [ I
1 v umoves '
Equation
3-12
Thus, for each measurement we have an estimated emission rate for both 15 minutes and 60
minutes, at sea level, with 9 RVP fuel. As a quality check, the results of our 15-minute
emissions to 60-minute conversion and the results for data at both durations are plotted in
Figure 17. As expected, the estimated hourly emissions (red circles) from the 15-minute
measurements closely match the measurements (blue triangles) where data at both test lengths
were available.
Figure 17: Hot Soak Measurement Test Length
10
0
2 101
=5
(/)
03
CD
| 10°
10
10"
A •
Legend
~ 1 hour estimated
A Data available for both lengths
10"
10" 10"
15 minute Measurement (g)
—r~
101
Quality assurance checks were also performed on the emissions values before and after
calculating their equivalences at Sea Level and 9.0 psi fuel. As expected, the tests measured
with higher RVP fuels at high altitude were reduced by wider margins under the influence of
the two corrections.
41
-------
Figure 18: Hot Soak Measurement Normalization to 9.0 RVP
RVP
1 1 1 1 1 1 1 1
0 5 10 15 20 25 30 35
Normalized Emissions (g)
3.3.3.4 Extrapolation to Missing Model Years and Ages
After normalizing the complete dataset, the results were incorporated into the MOVES
database. In the MOVES emission rates tables, emission rates must exist for all model year
and age group combinations. As with most cross-sectional datasets, this required additional
modeling. For example, there is no data for 20-year-old, model year 2010 vehicles, or brand
new 1980 vehicles. To address this problem, we extrapolated the emission rate values. Table
18 describes the data.
Table 18: Hot Soak Measurements by Model Year and Age
Age Group
0-3
4-5
6-7
8-9
10-14
15-19
20+
Leak?
N
Y
N
Y
N
Y
N
Y
N
Y
N
Y
N
Y
Total
1961-1970
5
5
re
1971-1980
8
8
>¦ =
1981-1995
6
15
46
55
8
39
169
« £
^ c
1996-2003
1
26
6
36
6
53
30
158
o v-
s
2004-2010
12
3
26
2
5
48
Total
12
3
27
2
31
6
36
6
59
45
46
55
8
52
388
In ranges where no data was collected, leak and non-leak measurements are extrapolated from
similar MY/age groups. In MY/age groups where very small amounts of data were collected,
the measurements are combined with similar MY/age groups. Figure 18 illustrates how we
populated model year and age group emission rates where there was no data.
42
-------
Figure 19: Measurement Averaging
ro
m
>
Q.
3
13
O
T3
O
(J
S
1961-1970
2004-2010
0-3
Age G roups Age Groups
4-5 6-7 8-9 10-14 15-19 20-99 0-3 4-5 6-7 8-9 10-14 15-19 20-99
(a) Non-Leak
(b) Leak
• A darker shaded cell represents a bin where data is present.
• An enclosed area represents one rate. The rate is calculated by averaging all enclosed
data.
For example, one non-leak rate exists for model years 1996-2003, ages 0-7. The rate is
calculated by averaging available data, which only exists at age 6-7. For every model year and
age group, there is a leaking rate and non-leaking rate. The two rates, weighted by leak
prevalence, form the average hourly hot soak emission rate for a given bin. Figure 20
demonstrates how leak rates and non-leak rates are combined to form a final weighted rate for
a given model year and age combination.
43
-------
Figure 20: Calculate Weighted Evaporative Emissions
> lOg / Qhr > 5g / Qhr >2g/Qhr >0.3g/Qhr <0.3g/Qhr
n
Leak Measurements
mean
mean
I
mean
I
mean
X X X X
Prevalence Prevalence Prevalence Prevalence
Leak Prevalence
Non-leak
measurem
ents
* weighted
mean non-
5 leak rate
leak rate
Non-Leak Prevalence
ij
MOVES
-> Emission
Rate
For every model year and age group combination, the calculation outlined in Figure 20 is
performed. Figure 21 shows the Hot Soak rates. The inclusion of leaking vehicle resulted in
higher emissions, particularly for older model years where leaks are more prevalent.
44
-------
Figure 21: Hot Soak Emission Base Rates (9.0 RVP at Sea Level)
20
15'
E 10
LU
1961-1970
1971-1980
1981-1995
1996-2003
2004-2010
MYG
10 15 20 0
10 15 20
0 5 10 15 20
Vehicle Age
10 15 20 0
10 15 20
1961-1970
1971-1980
1981-1995
1996-2003
2004-2010
3.3.4 Running Loss
3.3.4.1 Pre-Tier 2 Emission Rates
Running Loss emissions consist of vapor venting during vehicle operation. Data used to
develop running loss emission rates for Pre-Tier 2 vehicles is from CRC E-3513 and CRC E-
41.14-15 These two programs tested 200 vehicles with model years ranging from 1971-1997.
For each vehicle, fuel tank temperature was calculated at the end of the running loss test using
the fuel tank temperature algorithm (See Section 3.1). The running loss test performed in E-35
consisted of a UDDS, a two minute idle, a New York City Cycle (NYCC), a two minute idle, a
second NYCC, a two minute idle, a UDDS, and a final two minute idle.
The data was filtered/reduced such that each test meets the following requirements:
1. Non-liquid-leakers (emissions <137.2 g/hourc)
• As received vehicles (no retests)
• Fuel system pressure test result must be pass, fail, or blank
The average tank temperature was calculated by assuming a linear increase in temperature.
c Converted from 7.0 g/mile used in MOBILE6 31
45
-------
Thus, the average is calculated by averaging the start temperature of the test and the final
temperature. The average temperature is used to estimate the permeation rate using default
permeation rates and the permeation temperature adjustment.
Gram/hour rates were calculated by dividing total emissions by the duration of the running loss
test (4300 seconds). Permeation was subtracted for each hour to segregate tank vapor venting
(TVV) emissions. After analysis of the TVV data, running loss TVV rates were distinguished
by model year only. Table 19 shows the results of the analysis.
An I/M effect was not observable from this data, so the MOVES running loss TVV rates for
I/M and non-I/M are the same.
Table 19: Pre-Tier 2 Running Loss Emission Rates by Model Year and Age
Model year group
TW mean [g/hr]
Pre-1971
12.59
1971-1977
12.59
1978-1995
11.6
1996-2003
0.72
3.3.4.2 Tier 2 & Later Emission Rates
Running loss emission rates for Tier 2 and later vehicles were developed from a 2014 study of
five Tier 2 vehicles.25 In this study, vehicles were tested at two fuel RVP levels (7.51 psi and
10.33 psi) with and without implanted vapor leaks. Vapor leaks were installed at either the
canister or top of fuel tank, and with either 0.020" or 0.040" diameters, for a total of 4
possible leak configurations. The canister and fuel tank locations were chosen due to their high
rate of occurrence in the fleet.34
MOVES running loss emission rates are expressed in grams per hour and with a fuel vapor
pressure of 9 psi. Results from this testing are expressed in grams per test (4300 seconds) and
at two fuel vapor pressures (7.51 and 10.33). Therefore, the reported results must be
normalized to MOVES dimensions.
As in the development of Pre-Tier 2 emission rates, gram/hour rates were calculated by
dividing total emissions by the duration of the running loss test (4300 seconds). The
measurements were then adjusted to a 9-RVP equivalent emissions measurement using the
equations and coefficients described in Section 3.3.4.3.
Because our determination of a given vapor leak's rate of occurrence among all vapor leaks is
based on its hot soak emissions, each running loss test was immediately followed by a standard
one-hour hot soak procedure. Using the same process as in Section 3.3.3.2, the one-hour hot
soak results were multiplied by .54 to estimate the emissions at the 15-minute point. With this
result, each measurement was binned and the weighted average leak emissions rate
determined.
46
-------
Using the average non-leak value, the weighted average leak value, and the leak prevalences
from Section 3.3.2.7, an average emissions rate is calculated. Tier 2 and later running loss
emission rates are the first running loss rates in MOVES to account for vapor leak emissions
Tier 3 rates were estimated using ratios to Tier 2 standards. The resulting running rates are
shown in Figure 22.
Figure 22: Tier 2 & Tier 3 Running Loss Rates
Vehicle Age
3.3.4.3 Running Loss Fuel & Temperature Effects
Running Losses are affected by both temperature and fuel Reid Vapor Pressure (RVP). The
adjustments used in MOVES are taken from MOBILE6 and are applied to all model years and
source types. MOBILE6 was run for a series of temperatures and RVP levels for passenger
cars. A linear model was fit to the MOBILE6 results. The mean base emission rates for
running losses in MOVES are recorded in the 'EmissionRciteByAge ' table. Running loss
rates were assumed to be measured at 9 RVP and 95°F. The results from MOBILE6 were
normalized to the MOVES emission rates as multiplicative adjustments to the mean base rates.
For example, a multiplicative adjustment of 1 would be applied to a 9 RVP fuel at 95°F.
The running loss adjustments:
• Are multiplicative adjustments.
• Apply to all gasoline source types and model years.
• Are the same at temperatures at or below 40 F.
47
-------
• Are applied as a function of both RVP and ambient temperature.
• Use the 7 RVP coefficients for RVP values below 7 psi.
• Use the 10 RVP coefficients for RVP values above 10 psi.
• Are not adjusted for RVP at temperatures of 40 F or below.
AdjustedRunningLoss
= Running Loss * Adjustment (Temperature, RVP )
Equation
3-13
The adjustment coefficients are recorded in the evapRvpTemperatureAdjustment table in the
default MOVES database. The RVP adjustment range is dynamic; if new sets of coefficients
for RVP values greater than 10 or less than 7 are added to the table, MOVES will use those
values and set new minimum and maximum RVP values. Figure 23 shows the calculated
adjustment values.
Figure 23: Running Loss Temperature and RVP Effect
1.2
1.1
0.9
0.8
0.7
0.6
40 50 60 70 80 90
Ambient Temperature (degrees Fahrenheit)
100
3.4 Inspection/Maintenance (l/M) Program Effects
Inspection and Maintenance program efforts vary in their procedures for testing evaporative
emissions. Some locations use a fill pipe pressure check and gas cap check, others use just a
48
-------
scan of the onboard diagnostics (OBD), and others will use all three approaches. These types
of tests are not expected to reduce permeation or liquid leaks and do not guarantee the
detection of a vapor leak within a vehicle. MOVES assumes tank vapor venting is the only
evaporative process where I/M benefits are realized.
The I/M effects for evaporative emissions are calculated as a function of three types of
information:
• I/M program coverage, which describes program coverage by location and year (see
the MOVES Emissions Adjustments report36),
• I/M factor which describes the relative effectiveness of different I/M program designs
for different vehicles (see Section 3.4.1)
• Vapor venting emission rates with a reference I/M program, in comparison to rates
with no program. For the cold soak operating mode (opModelD =151) These are
calculated based on differences in leak prevalence as explained in Section 3.4.2. For
hot soak (opModeID=150) and running loss (opModeID=300), they are calculated
based on rates in the EmissionRateByAge table.
3.4.1 I/M Factor (Relative Program Effectiveness)
In MOVES, the I/M factor describes the relative effectiveness of an I/M program; a higher I/M
factor indicates a more effective I/M program. This data is stored in the MOVES IMFactor
table for gasoline and E85 passenger cars, passenger trucks and light commercial vehicles. The
I/M factors vary with inspection frequency, type of test (testStandardsID), model year group
and vehicle age.
Data from four I/M programs were used in the development of MOVES I/M factors for
evaporative emissions. The Phoenix, AZ program contained the most extensive data, so we
have used it to represent a reference condition relative to which other programs can be
assessed. Data from programs in Tucson, AZ, Colorado, and North Carolina were used to
determine the effectiveness of other I/M program designs. See Table 20.
Table 20: Description of I/M Programs 37
Location
Gas Cap
Test
OBD
Pressure
test
Frequency
Network
Years
Colorado
Y
Advisory
N
Biennial
Hybrid
Annual
2003-2006
N.Carolina
N
Y
N
Annual
Annual
2002-2006
Phoenix
Y
Y
Y
Biennial
Centralized
2002-2006
Tucson
Y
Y
N
Annual
Centralized
2002-2006
In order to develop I/M factors, failure data was used from I/M programs. These failure
frequencies were only used to estimate the relative effectiveness of differing evaporative I/M
49
-------
programs. They were not used to model the actual prevalence of evaporative leaks. For
information on the modeling of vapor leak prevalence please see Section 3.3.2.7. For
information on liquid leaks, see Section 3.5.
The Phoenix evaporative I/M program performed gas-cap tests on all vehicles, OBD scans on
OBD-equipped vehicles, and fill-pipe pressure tests on pre-OBD vehicles. The OBD codes
used to assign evaporative failures are listed in Table 21for all vehicle makes and additionally
P1456 and P1457 for Honda and Acura vehicles. Vehicles with one or more of these faults
were flagged as failing vehicles, analogous to pre-OBD vehicles that failed the pressure test.
Very few vehicles failed both the gas cap test and the pressure/OBD test. Therefore, the total
number of failures is the sum of gas cap and pressure/OBD failures.
Table 21: OBD Evaporative Emission Trouble Codes
OBD Code
Description
P0440
Evaporative Emission Control System Malfunction
P0442
Evaporative Emission Control System Leak Detected (small leak)
P0445
Evaporative Emission Control System Purge Control Valve Circuit Shorted
P0446
Evaporative Emission Control System Vent Control Circuit Malfunction
P0447
Evaporative Emission Control System Vent Control Circuit Open
P1456
EVAP Emission Control System Leak Detected (Fuel Tank System)
P1457
EVAP Emission Control System Leak Detected (Control Canister System)
To determine the reference program effectiveness, we estimated failure frequencies with I/M
and without I/M ("non-I/M"). The I/M failure frequencies were developed from the Phoenix
data using initial and final results for a vehicle in a given I/M cycle. For passing vehicles, the
initial and final tests are the same. The initial and final failure frequencies were averaged to
develop an I/M failure frequency for each model year and age group. Using the initial failure
frequencies alone would neglect the required repairs occurring on most failing vehicles, while
using only final failure frequencies would neglect the prior existence of failing vehicles.
To develop non-I/M failure frequencies, the sample was restricted to vehicles that were
inspected, but were registered in states that do not have any I/M programs.
The Tucson data was used to determine the effect of I/M program frequency (annual vs.
biennial). For OBD-equipped vehicles, Tucson performs gas-cap and OBD tests annually,
while Phoenix performs them biennially. Therefore, we were able to develop an the
effectiveness ratio for Annual/Biennial programs by analyzing the Tucson data.
The North Carolina data was used to estimate the effectiveness of using the OBD scan as the
sole test in a program. In North Carolina, expansion of I/M program boundaries led to many
vehicles being tested for the first time. These vehicles were effectively non-I/M until their first
test. Vehicles were flagged as non-I/M tests if they were tested before the official start of the
I/M program or were registered in a new I/M county. Failure frequencies of the non-I/M
vehicles were compared to vehicles tested in I/M areas. The I/M effectiveness of an OBD only
I/M program is estimated to be a 63 percent reduction in failures or a non-I/M to I/M failure
50
-------
ratio of 1.6. This ratio was then applied to Phoenix OBD and pressure test failure frequencies
to determine non-I/M failure frequencies.
The Colorado data was used to determine the effectiveness of gas cap tests. In Colorado, the
I/M data is primarily from the Denver and Boulder metropolitan areas. However, some
residents are new to this area, having moved from non-I/M counties and states. These vehicles
were effectively non-I/M until their first test. Vehicles were flagged as non-I/M if they were
registered in a state without an I/M program, or in a non-I/M county within Colorado.
Colorado OBD data was not used, because OBD in Colorado is only advisory and does not
pass or fail a vehicle. The failure rates of the non-I/M vehicles were compared to those in the
I/M fleet. The effectiveness of a gas-cap-only I/M program is estimated to be a 45 percent
reduction in failures or a non-I/M to I/M failure ratio of 1.2. This ratio was then applied to gas
cap failure frequencies to determine non-I/M failure frequencies.
The I/M factor in MOVES adjusts emission rates depending on the characteristics of a given
county's I/M program. Our reference program, Phoenix, has an IM factor of 1. Non-I/M areas
have an IM factor of 0. The failure frequencies from the other counties are used to calculate
I/M factors for the diverse types of evaporative I/M procedures. The I/M factor is assumed to
have a linear relationship with failure frequency. Figure 24 illustrates how the I/M factor varies
with different I/M programs. Different programs fall on the line as determined by the analysis
above, based on specific evaporative tests performed. For the vehicles in Figure 24, Tucson's
OBD and gas cap tests are annual, compared to Phoenix's biennial requirement, which gives
Tucson a lower failure frequency, thus a higher I/M factor. Colorado's frequency is biennial,
but their OBD test is non-enforcing. As a result, their data shows a higher failure frequency,
resulting in a lower I/M factor.
51
-------
Figure 24: I/M Factor, MY 1999-2003, Age 4-5
l/M Factor
3.4.2 Leak Prevalence
To estimate the impact of an I/M program, the appropriate I/M factor is applied to the
estimated vapor leak prevalence rates. However, because the leak prevalence rates were
developed from a test program in the Denver, CO area, and because, as explained above, the
Denver program is not the MOVES reference program, the Denver vapor leak prevalence rates,
developed in Section 3.3.2.7 must be adjusted for use in MOVES.
As illustrated in Figure 24, the with-I/M failure frequency in Denver is about 30 percent less
than non-IM (I/M factor = 0) and 30 percent higher than Phoenix (I/M factor = 1) so the leak
prevalence rates developed from Denver data were adjusted accordingly before being added to
the MOVES database. This adjustment reflects the analysis described in the previous section
and can be observed in Figure 25. During a MOVES run for the Denver area, the Denver I/M
factor will be applied to the adjusted leak prevalence rates and emissions will be modeled with
the same prevalence rates originally estimated for Denver.
52
-------
Figure 25: Adjusting Denver Leak Prevalence Data
Denver Leak Frequency
Denver Non-Failure
Frequency
Phoenix (I/M)
noa-I/M
3.5 Liquid Leaks
Liquid leaks include any non-vapor fuel escaping the fuel system. The average liquid leaking
rate is determined using the leaking vehicles excluded from the analysis above. Because the
testing methods used did not distinguish the different evaporative emission processes,
permeation and tank vapor venting are estimated using the calculation methods described in
Section 3.2 and Section 3.3 and subtracted from the total measurement. The remaining
emissions after permeation and vapor venting are subtracted are assumed to be caused by
liquid leaks. Due to limitations in the data quality and quantity, the measurements are
averaged across all vehicles by the three different modes and shown in Table 22.
Table 22: Liquid Leak Emission Rates (g/hr)
Operating Mode
Liquid leak rate
Cold Soak
9.85
Hot Soak
19.0
Operating
178
53
-------
The liquid leak emission rates must be multiplied by the percentage of leakers in the fleet to
get an average liquid leaking emission rate. The studies by BAR 20 and API 21 provided
this data. The estimates of liquid leak prevalence are shown in Table 23 It is assumed that
most leaks do not occur until vehicles are 15 years or older.
Table 23: Percentage of Liquid Leaks by Age
Age group
Percentage of leakers in fleet
0-9
0.09%
10-14
0.25 %
15-19
0.77%
20+
2.38%
Table 24 contains the fleet-weighted liquid leak rate. There is insufficient data to conclude
that these rates change with model year or are affected by I/M programs.
Table 24: Weighted Liquid Leak Emissions(g/hr)
Age group
Cold soak
Hot soak
Operating
0-9
0.009
0.017
0.158
10-14
0.025
0.048
0.450
15-19
0.075
0.145
1.360
20+
0.235
0.452
4.230
As with vapor leaks, we expect a reduction in liquid leak prevalence due to improved system
design and integrity under the Tier 3 regulations. However, liquid leaks in advanced evaporative
systems are primarily caused by tampering and mal-maintenance. Therefore, we estimate Tier
3 will prevent half as many liquid leaks as vapor leaks, shown in Table 25.
Table 25: Weighted Tier 3 Liquid Leak Emissions (g/hr)
Age group
Cold soak
Hot soak
Operating
0-9
0.007
0.013
0.123
10-14
0.019
0.037
0.342
15-19
0.058
0.113
1.054
20+
0.180
0.348
3.258
3.6 Refueling
Refueling emissions are the fuel vapors when liquid fuel is added to the tank. In MOVES these
are reported as Refueling Displacement Vapor Loss (processid =18) and Refueling Spillage
54
-------
Loss (processid = 19). Refueling emissions are estimated from the total volume of fuel
dispensed (gallons). This volume is estimated from the average daily distance travelled (VMT)
and the estimated fuel consumption. Both the spillage and the vapor displacement associated
with refueling events are in terms of grams per gallon of fuel dispensed. Diesel vehicles are
assumed to have negligible vapor displacement, but fuel spillage is calculated for diesel
vehicles.
Two emission control strategies exist to limit fuel lost during refueling for gasoline and E85
vehicles. First, there are programs designed to capture refueling vapors at the pump. These are
often referred to as Stage II vapor control programs and vary by location. Second, vehicles
manufactured since 1998 have onboard refueling vapor recovery (ORVR) systems that store
refueling vapors in the vehicle's evaporative emission canister. Stage II programs are assumed
to reduce emissions from vehicles not equipped with ORVR and additionally, any refueling
emissions that are not captured by the ORVR systems. MOVES does not account for any
interaction between ORVR systems and gasoline dispensing stations equipped with Stage II
equipment. The technology adjustment for both refueling vapor loss and spillage loss is the
same in all locations.
The implementations of Stage II systems vary from area to area and affect the fuel vapor
displacement and the amount of spillage. MOVES uses two factors to adjust the refueling losses
and account for this variation.
1. The refueling vapor program adjustment is a value between zero and one indicating the
percent reduction of total potential vapor losses by state or local programs (such as Stage
II recovery programs).
2. The refueling spill program adjustment is a value between zero and one indicating the
percent reduction of refueling spillage losses by state or local programs (such as Stage II
recovery programs).
These program adjustments in MOVES are applied by county and calendar year. Each county
has a unique value for vapor and spillage program adjustments. The MOVES default database
contains information about all of the existing Stage II programs by county based on the
parameters used for the 2005 National Emission Inventory (NEI), with updates accounting for
states that have decommissioned Stage II programs through early 2024. The program
adjustment values for each county and calendar year are stored in the default MOVES
'CountyYear' table. Table 26 lists the first and last year of Stage II coverage by state as
recorded in MOVES.
55
-------
Table 26: Summary of Stage II Programs by State
State*
Begin Year"
End Year***
Arizona
2011
2017
California
1999
N/A
Connecticut
1999
2014
Delaware
1999
2020
District of Columbia
1999
2021
Florida
1999
2010
Georgia
1999
2015
Illinois
1999
2016
Indiana
1999
2016
Kentucky
1999
2017
Louisiana
1999
2021
Maine
1999
2012
Maryland
1999
2018
Massachusetts
1999
2016
Missouri
1999
2014
Nevada
1999
N/A
New Hampshire
1999
2014
New lersey
1999
2019
New York
1999
2022
Ohio
1999
2016
Oregon
1999
N/A
Pennsylvania
1999
2021
Rhode Island
1999
2016
Tennessee
1999
2018
Texas
1999
2017
Vermont
1999
2014
Virginia
1999
2016
Washington
1999
2022
Wisconsin
1999
2016
Stage II programs are applied by county in MOVES. This table provides a summary of
which states contain at least one county with a Stage II program.
MOVES cannot model calendar years between 1991-1999. As such, this table is not
intended to indicate when Stage II programs began, but simply the calendar years in
which MOVES will model a Stage II program.
*" In MOVES, Stage II programs that have not ended are modeled as continuing for all
future years (i.e. through 2060).
MOVES uses separate factors to address the on-board refueling vapor recovery (ORVR)
systems on gasoline and E85 vehicles. The effects of ORVR technology is phased in beginning
in model year 1998, as shown in Table 27. There is no ORVR for motorcycles. In MOVES,
ORVR reduces emissions from both vapor and spillage.
56
-------
Table 27: Phase-In of Onboard Refueling Vapor Recovery
Model
Passenger
Light
Light
Heavy Duty
Heavy Duty
Heavy Duty
Year
Cars
Trucks
Trucks
Trucks
Trucks
Trucks
<6,000 lbs
6,000-
8,501-
10,001-
>14,000 lbs
GVWR
8,500 lbs
10,000 lbs
14,000 lbs
GVWR
GVWR
GVWR
GVWR
1998
40%
0%
0%
0%
0%
0%
1999
80%
0%
0%
0%
0%
0%
2000
100%
0%
0%
0%
0%
0%
2001
100%
40%
0%
0%
0%
0%
2002
100%
80%
0%
0%
0%
0%
2003
100%
100%
0%
0%
0%
0%
2004
100%
100%
40%
40%
0%
0%
2005
100%
100%
80%
80%
0%
0%
2006-
100%
100%
100%
100%
0%
0%
2017
2018-
100%
100%
100%
100%
100%
0%
2026
2027+
100%
100%
100%
100%
100%a
100%a
aNote that for MY2027+, MOVES models a 100% ORVR requirement for all gasoline and E85 vehicles other
than motorcycles. However, the requirement for engine-certified Class 3 trucks (10,001-14,000 lbs GVWR]
does not actually begin until model year 2030. This modeling simplification underestimates refueling
emissions for a small number of MY 2007-2029 regulatory class 42 vehicles.
3.6.1 Vapor Loss
As noted above, the Refueling Displacement Vapor Loss process (processID 18) estimates the
displaced grams of fuel vapor per gallon of liquid fuel consumed. For an uncontrolled system,
the vapor displaced by refueling with gasoline or E85 is a function of the fuel's Reid Vapor
Pressure (RVP), the dispensed gasoline temperature, and the temperature difference between
the vehicle's fuel tank and the dispensed gasoline, as shown in Equation 3-14:38
DVR = exp (-1.2798 - 0.0049dT + 0.02037^ + 0.1315WP) Equation
3-14
Where:
DVR
df
Tdf
RVP
= base displaced vapor rate (g/gal)
= difference in temperature between the vehicle's fuel tank and the dispensed fuel
(°F) (see Equation 3-16)
= temperature of the dispensed fuel (°F) (see Equation 3-15)
= Reid Vapor Pressure of the dispensed fuel (psi)
In MOVES, the coefficients of Equation 3-14 are stored in the RefuelingFactors table of the
default database as vaporTermA, vaporTermB, vaporTermC, and vaporTermD, respectively.
57
-------
The temperature of the dispensed fuel (TDF) depends on the temperature of the fueling
station's storage tank temperature, which is modeled as a function of the ambient temperature,
as shown in Equation 3-15:39
Tdf — 20.3 + 0.81 X Tamb
Equation
3-15
Where:
Tdf = temperature of the dispensed fuel (°F)
Tamb = ambient temperature (°F), bounded between 45°F and 90°Fd
In MOVES, the lower bound of 45°F and upper bound of 90°F are stored in the
RefuelingFactors table of the default database as vaporLowTLimit and vaporHighTLimit,
respectively. The coefficients of Equation 3-15 are not stored in the default database; instead,
they are hard coded directly in the MOVES algorithm. The ambient temperature (Tamb) is
stored in the ZoneMonthHour table of the default database, which varies by county, month,
and hour of day.
The difference in temperature between the vehicle's fuel tank and the dispensed fuel (TDF) is
calculated as shown in Equation 3-16:40
dT = 0.418 X Tdf - 16.6 Equation
3-16
Where:
dT = difference in temperature between the vehicle's fuel tank and the dispensed fuel
(°F), up to a maximum difference of 20°F
Tdf = temperature of the dispensed fuel (°F)
In MOVES, the coefficients of Equation 3-16 are stored in the RefuelingFactors table of the
default database as vaporTermE and vaporTermF, respectively. The upper bound of dT is
stored as tankTDiffLimit. This upper bound is implemented because in the Amoco study40 that
this equation relies on, the difference in temperature was never greater than 20 degrees.
Base displaced vapor values calculated for a warm summer day with an ambient temperature
of 90°F and 7.8 psi RVP fuel are about 4.7 grams per gallon. With unseasonable weather,
extreme values could range between about 2 grams per gallon with low temperatures (45°F)
and low summer RVP values (6.5 psi), and 9.2 grams per gallon with high temperatures (90°F)
and high winter RVP (13 psi).
We limit the base displaced vapor rate to a minimum value of 1.8 grams per gallon, which is
stored in the RefuelingFactors table of the default database as minimumRefuelingVaporLoss.
This is equivalent to the value calculated with a low ambient temperature (less than 45°F) and
d The lack of data on dispensed fuel temperatures for ambient temperatures above 90 degrees Fahrenheit is a
limitation in our modeling that may be important in areas with sustained high temperatures.
58
-------
an RVP of 5.7.
After calculating the base displaced vapor rate for uncontrolled systems, we calculate an
adjusted displaced vapor rate, which accounts for controls on displaced vapors from refueling.
This step is presented as two separate calculations for clarity, as shown in Equation 3-17 and
Equation 3-18.
Equation 3-17 calculates the adjusted displaced vapor rate for uncontrolled vehicles, which
accounts for Stage II controls and the fraction of vehicles without functioning ORVR controls:
ADVRy = DVR x (1 - S2C) X (1 - RTA)
3-17
Where:
ADVRu = adjusted displaced vapor rate for uncontrolled vehicles (g/gal)
DVR = base displaced vapor rate (g/gal)
S2C = Stage II controls described as refueling vapor program adjustment above,
dimensionless
RTA = refueling technology adjustment, dimensionless
For the above equation, DVR is calculated as shown in Equation 3-14. S2C represents the
presence and effectiveness of Stage II controls. A hypothetical value of 1 represents full
penetration and effectiveness of Stage II controls, which would result in an ADVRy of 0 g/gal.
Conversely, a value of 0 for S2C represents a lack of Stage II controls, which would result in
no adjustment to the base displaced vapor rate due to Stage II.
RTA is stored in the RefuelingControlTechnology table of the default database, and is
determined by regulatory class, source type, fuel type, vehicle age, and model year. This value
incorporates both the ORVR technology adoption rates from Table 27 and information on real-
world ORVR function derived from the "High Evaporated Emission Investigation Field
Study."41 This EPA study investigated the percentage of vehicles which had observable
refueling emissions plumes. A model year-dependent regression was developed to determine
the percent increase in observed plumes (and thus the decrease in ORVR effectiveness) as a
proxy of age. The decrease was calculated as -0.68 percent per year6 which was applied for
ages 0 through 19f and used for all model years, regulatory classes, and source types with
ORVR implemented. The RTA value is the same in all locations.
Since RTA essentially represents the penetration of functioning ORVR technology, 1 — RTA
represents the fraction of vehicles that do not have functioning ORVR and is thus used to
adjust the base displaced vapor rate to calculate ADVRy.
Equation 3-18 calculates the adjusted displaced vapor rate for the fraction of ORVR-controlled
e The MOVES value is slightly different than that listed in the report (-0.56%/ year) because the MOVES value
included an extra year and some vehicles determined to be light-duty after additional VIN decoding work.
f For ages 19 and older we set a minimum value of 87 percent since this was the oldest age with fully phased-in
ORVR in the field study.
59
-------
vehicles:
ADVRC = CRR X (1 - S2C) X RTA
Equation
3-18
Where:
ADVRC = adjusted displaced vapor rate for ORVR-controlled vehicles (g/gal)
CRR = controlled refueling rate (g/gal)
S2C = Stage II controls described as refueling vapor program adjustment above,
dimensionless
RTA = refueling technology adjustment, dimensionless
In the above equation, CRR represents the refueling vapor losses in a controlled ORVR
system. This value is based on measurements from the In-Use Verification Program (IUVP).
The mean HC emissions of the functional ORVR systems (defined as those with emission
below the certification level of 0.2 g/gal) as measured in the 2009-2020 light-duty IUVP data
was 0.0361 g/gal.42 This value was used for all vehicles with controlled ORVR systems. It is
stored in the RefuelingControlTechnology table of the default database, but does not vary by
age, model year, regulatory class, or source type.
Equation 3-18 has the same Stage II control adjustment for CRR as Equation 3-17 has for
DVR, as Stage II programs are assumed to reduce any refueling emissions not captured by
ORVR systems.
As described above, RTA represents the penetration of functioning ORVR technology and is
therefore applied directly to adjust the controlled refueling rate to calculate ADVRC.
The total adjusted displaced vapor rate is then calculated by combining Equation 3-17 and
Equation 3-18, as shown in Equation 3-19:
Where:
ADVR = total adjusted displaced vapor rate (g/gal)
ADVRu = adjusted displaced vapor rate for uncontrolled vehicles (g/gal)
ADVRC = adjusted displaced vapor rate for ORVR-controlled vehicles (g/gal)
The result of Equation 3-19 is an adjusted vapor rate of grams of HC emissions per gallon of
fuel consumed. While MOVES does not directly calculate fuel consumption, it does calculate
total energy consumption. To calculate refueling displacement vapor loss, MOVES
implements Equation 3-20:
ADVR = ADVRy + ADVRC
Equation
3-19
TEC
Equation
RDVL = ADVR X
Energy/Content x FuelDensity
3-20
60
-------
Where:
RDVL = refueling displacement vapor loss (g)
ADVR= total adjusted displaced vapor rate (g/gal)
TEC = total energy consumption (kJ)
EC = energy content of the fuel (kJ/g)
FD = fuel density (g/gal)
The total energy consumption (TEC) is calculated as described in the MOVES GHG and
Energy Consumption technical report.43 The fuel energy content (EC) and density (FD) are
stored in the FuelSubType and FuelType tables (respectively) of the default database and are
also documented in the same technical report.
3.6.2 Spillage
AP-42 Volume I Section 5.2.2.344 lists the spillage as 0.7 lb/1000 gallons, which is
0.31g/gallon of dispensed fuel. MOVES uses this value for uncontrolled spillage from
gasoline, E85, and diesel fueled-vehicles.
MOVES adjusts spillage to account for Stage II programs with refuelingSpillProgram
Adjust, a value between zero and one as listed in the county Year table.
MOVES also adjusts spillage to account for ORVR. In MOVES4 and later, the refueling tech
adjustment for refueling spillage losses reflect the fuel type and regulatory class
population fractions® from MOVES3.0, coupled with an assumption that gasoline vehicle
spillage emissions are reduced using a 50% technology adjustment factor when ORVR
technology is fully implemented. They are listed in the MOVES
'SourceTypeTechAdjustment table. The technology adjustment values take into account the
technology adoption rates from Table 27 and the population fraction of the vehicles by
fuel type and regulatory class (GVWR for heavy-duty vehicles) within each source type
as documented in the MOVES3 Population and Activity Report.45 Note that the factors
for heavy-duty trucks 8,500 lbs and higher are reflected in the values for single-unit short-
haul and long-haul vehicles (source types 52 and 53). For example, LHD2b3 vehicles
(GVWR < 14,000 lbs) compose 49% of the population of gasoline source type 52 and 53
vehicles in model year 2014 and later in MOVES3.h
Refueling Spillage Loss = UncontrolledSpillage* (1- refuelingSpillProgramAdjust) * Equation
(1-refuelingTechAdjustment) 3-21
g The population weightings in SourceTypeTechAdjustment are consistent with MOVES3. We hope to update the
spillage calculations in a future version of MOVES.
61
-------
Appendix A Notes on Selected Evaporative Emission Data
Parameters: Vehicle Numbers, Test No., Ambient Temperature, RYP, Model Year, Fuel
System, Purge, Pressure, Canister, Gram HC, Retest
E-41 CRC Late Model In-Use Evap. Emission Hot Soak Study (1998)
• 50 vehicles (30 passenger cars and 20 light duty trucks)
• Model years 1992 to 1997
• Average RYP: 6.5 psi
• Diurnal Temperature: 72 to96°F
• Fuel System: Port Fuel Injection, Throttle Body Injection
• Vehicle fuel tank drained and refilled to 40% of capacity with Federal Evaporative
Emission
Test Fuel
• Driving schedule will be a full LA-4-NYCC-NYCC-LA4 sequence, with two minute idle
periods
following the first LA-4, the second NYCC, and the final LA-4.
• Hydrocarbon readings will be taken continuously throughout the running loss test.
• Cumulative mass emissions will be reported at one-minute intervals.
• Ambient Temperature in running loss enclosure: 95 °F
E-9 CRC Real Time Diurnal Study (1996)
• 151 vehicles (51 vehicles MY 1971-1977, 50 vehicles MY 1980-1985, 50 vehicles MY 1986-
1991)
• Odometers range from 39,000 to 439,000 miles
• Fuel tank volume was 15% of the rated capacity
• RVP: 6.62 psi (average sum of 47 vehicles)
• Diurnal temperature: 72 to 96°F
• Fuel System: Port Fuel Injection, Carburetor, Throttle Body Injection
CRC E-35 Running Loss Study (1997)
• 150 vehicles (50 vehicles MY 1971-1977, 50 vehicles MY 1980-1985, 50 vehicles MY 1986-
1991)
62
-------
• Ambient Temperature in running loss enclosure: 95 °F
• RYP: 6.8 psi
• Fuel System: Port Fuel Injection, Carburetor, Throttle Body Injection
EPA Compliance Data
• 2-Day Test
• Length of the hot soak: 1 hour
• 77 vehicles
• RYP: average 8.81 psi
• Ambient Temperature:
• Federal Standard (72 to 96° F) Diurnal
• Cal. (65 to 105°F)Diurnal
Hot Soak: 81.67'F
• Fuel System: Port Fuel Injection
MSOD (Mobile Source Observation Database):
Hot Soak 1 hour hot soak evaporative test
FTP Federal test procedure (19.53 mph), also referred to as the UDDP schedule
NYCC New York City Cycle Test (7.04 mph)
BL1A 1 hour Breathing Loss Evap. Test Gas Cap left On BL1B 1 hour Breathing Loss
Evap. Test Canister as reed. ST01 Engine Start cycle test
4HD 4 hour Diurnal test
24RTD 24 Hour Real Time Diurnal 33RTD 33 Hour Real Time Diurnal 72RTD 72 Hour
Real Time Diurnal
3Rest 3 Hour Resting Loss Evap. Emission Test (follows 1 HR Hot Soak) CY6084 Real time
diurnal temperature pattern: range 60 to 84 F CY7296 Real time diurnal temperature
pattern: range 72 to 96 F CY8210 Real time diurnal temperature pattern: range 82 to
102 F
DIURBL Standard temperature rise for 1 hour diurnal or breathing loss evaporative emis-
sion test
F505 Bag 1 of federal test procedure (25.55 mph)
ASM Acceleration Simulation Mode Test Procedure
ATD Ambient Temperature diurnal evaporative Test, shed temp constant, vehicle begins 24
degree cooler
63
-------
Appendix B Tank Fuel Generator
In MOVES, we account for changes to the RVP of gasoline fuels due to weathering and comingling as a
preliminary step before calculating vapor venting emissions. These calculations are handled in the Tank
Fuel Generator (TFG).
TFG-la: Calculate Average Pump Gasoline and Ethanol Blend Type
Inputs:
marketShare from FuelSupply (county, fuelYear, monthGroup, fuelFormulationID)
ETOHvolume from FuelFormulation (fuelFormulationID)
RVP from FuelFormulation (fuelFormulationID)
fuelSubTypelD from FuelFormulation (fuelFormulationID)
fuelTypelD from FuelSubType (fuelsubtypelD)
Outputs:
averageRVP (county, fuelYear, monthGroup)
tankAverageETOHVolume (county, fuelYear, monthGroup)
Calculations:
averageRVP = For all FuelFormulations in county, fuel year & monthGroup where
fuelType = "gasoline" (ie fuelTypelD =1))
(Sum (RVP*marketshare)) / (Sum (marketshare))
tankAverageETOHVolume = For all Fuel Formulations in county, fuel year & monthgroup
where fuelType = "gasoline" (ie fuelTypelD = 1))
(Sum (ETOH Volume*marketshare)) / (Sum (marketshare))
TFG-lb: Calculate Ethanol Market Share and Ethanol BlendType
Inputs:
marketShare from FuelSupply (county, fuelYear, monthGroup, fuelFormulation)
fuelSubTypelD from FuelFormulation
fuelTypelD from FuelSubType
Outputs:
gasoholMarketShare (countyID, fuelYearlD, monthGroupID)
ethanolBlendType (county, fuelYear, monthGroup)
Calculation:
gasoholMarketShare: For all FuelFormulations in county, fuelyear & monthgroup where
ETOHVolume >= 4
gasoholMarketShare =Sum (marketShare)
lowETOHRVP: For all FuelFormulations in county, fuel year & monthgroup WHERE
64
-------
fuelType = "gasoline" (ie fuelTypelD = 1) and ETOHVolume <4
IF (sum (marketshare) = 0,
lowETOHRVP=AverageRVP
ELSE
lowETOHRVP=(Sum (RVP*marketshare)) / (Sum (marketshare))
highETOHRVP: For all FuelFormulations in county, fuel year & monthgroup WHERE
fuelType = "gasoline" (ie fuelTypelD = 1)) and ETOHVolume >= 4
IF gasoholMarketShare = 0,
highETOHRVP = AverageGasolineBlendRVP
ELSE
highETOHRVP =(Sum (RVP*marketshare)) / gasoholMarketShare
ethanolBlendType:
IF absolute value (highETOHRVP -lowETOHRVP) <= 0.2,
ethanolBlendType ="Match"
ELSE ethanolBlendType = "Splash"
TFG-lc: Calculate Commingled Tank Fuel RVP
Inputs:
gasoholMarketShare (countylD, fuelYearlD, monthGroupID) from TFG-lb
averageRVP (countyID, fuelYearlD, monthGroupID) from TFG-la
Commingling Lookup (stored in MOVES)
LookupMarketShare
Commingling
RVP Factor
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
1.000
1.016
1.028
1.035
1.039
1.040
1.038
1.034
1.027
1.018
1.000
0.8
0.9
1.0
Outputs:
commingledRVP (countyID, fuelYearlD, monthGroupID)
65
-------
Calculation:
comminglingFactor (countylD, fiielyearlD, monthgroupID) = lookup from table using
smallest value of "LookupMarketShare" that is greater than or equal to the
gasoholMarketShare.
commingledRVP = averageRVP * comminglingFactor
TFG-2: WeatheredRVP
TFG-2a: Calculate "EvapTemp" by ZonelD, MonthGroupID
Inputs:
temperature (zonelD, hourlD monthgroupID)
zonelD from masterloopcontext
Outputs:
zoneEvapTemp (zonelD, monthgroupID)
Calculation :
zoneMin(zoneID, monthgroupID) = MIN (temperature(zoneID, monthgroupID, hourlD))
zoneMax (zonelD, monthgroupID) = MAX(temperature(zoneID, monthgroupID, hourlD)
zoneEvapTemp =
IF zoneMax <40 or zoneMax-zoneMin <=0, (zoneMin+zoneMax)/2
ELSE zoneEvapTemp(zoneID, monthGroupID) =
-1.7474+1.029*zoneMin+ 0.99202* (zoneMax-zoneMin)-0.0025173*zoneMin*
(zoneMax-zoneMin)
TFG-2b: Calculate ratio of weathering loss for gasoline by Zone, Year & Month at actual ambient
temperatures relative to a diurnal swing of 72-96 F
Inputs:
zoneEvapTemp (zonelD, monthgroupID) from previous step
commingledRVP (countylD, fuelYearlD, monthgroupID) from TFG-lc
zone(countyID, zonelD)
Outputs:
ratioGasolineRVPLoss(zoneID, fuelYearlD, monthgroupID)
Calculation :
ratioGasolineRVPLoss =MAX (0, [-2.4908 + 0.026196 * zoneEvapTemp + 0.00076898 *
zoneEvapTemp * commingledRVP]/[-0.0860 + 0.070592 * commingledRVP ] )
TFG-2c: Calculate weathering loss for average fuel for standard temperatures
66
-------
Inputs:
ethanolBlendType (county, fuelYear, monthGroup) from TFG-la
gasoholMarketShare (countylD, fuelYearlD, monthGroupID) from TFG-lb
Outputs:
avgWeatheringConstant (countylD, fuelYearlD, monthGroupID)
Calculations:
IF ethanolBlendType = "Match", avgWeatheringConstant = 0.049 - 0.0034 *
gasoholMarketShare
ELSE avgWeatheringConstant = 0.049 - 0.0116 * gasoholMarketShare
TFG-2d: Calculate weathered RVP for county-average fuel adjusted for zone temperatures
Inputs:
ratioGasolineRVPLoss (zonelD, fuelYearlD, monthgroupID) from TFG-2b
avgWeatheringConstant (countylD, fuelYearlD, monthGroupID)
from previous step
commingledRVP (countylD, fuelYearlD, monthGroupID) from TFG-lc
zone(countyID, zonelD)
Outputs:
tankAverageGasolineRVP(zoneID, fuelYearlD, monthgroupID)
Calculation :
tankAverageGasolineRVP (zoneID) = commingledRVP(countylD) * (1 -
ratioGasolineRVPLoss (zonelD) * avgWeatheringConstant (countylD))
67
-------
Appendix C List of Acronyms
API
American Petroleum Institute
BWC
Butane Working Capacity
CNG
Compressed Natural Gas
CRC
Coordinating Research Council
DTC
Diagnostic Trouble Code
EO
Gasoline containing 0 percent ethanol by volume
E10
Gasoline containing 10 percent ethanol by volume
E-65
CRC fuel permeation from automotive systems study
E-77-2
evaporative emission/permeation test program
EPA
U.S. Environmental Protection Agency
ETOH
Ethanol
FTP
Federal Test Procedure
HC
Hydrocarbons
HD
Heavy-Duty
I/M
Inspection and Maintenance program
LDGV
Light-Duty Gasoline Vehicle
LEV
Low emission vehicle
LEV III
California Tier 3 light-duty emission standards of 2012
MOBILE
Original Highway Vehicle Emission Factor Model pre-2004
MOBILE6
Versioned Highway Vehicle Emission Factor Model
MOVES
Motor Vehicle Emission Simulator Model
MSAT
Mobile Source Air Toxics
MTBE
Methyl tertiary-butyl ether
NMHC
Non-Methane Hydrocarbons
NMOG
Non-methane organic gases
OBD
Onboard Diagnostics
ORVR
Onboard Refueling Vapor Recovery
OTAQ
Office of Transportation and Air Quality
PI
Petroleum Institute
RVP
Reid Vapor Pressure
PSHED
Portable Sealed Housing for Evaporative Determination
SHED
Sealed Housing for Evaporative Determination
THC
Total Hydrocarbon
Tier 2
Vehicle emissions certification standards phased in from 2002
Tier 3
Vehicle emissions certification standards phased in from 2017
TOG
Total Organic Gases
TVG
Tank Vapor Generated
TVV
Tank Vapor Venting
voc
Volatile Organic Compound
VSP
Vehicle specific power
68
-------
Appendix D Glossary
backpurge - as the temperature decreases a vacuum is created in the fuel system which pulls
the hydrocarbons from the charcoal canister into the fuel tank, creating more space in the
canister for hydrocarbons to adhere during the next heating period
breakthrough - when the vapor generated by the fuel system overwhelms the charcoal can-
ister and uncontrolled hydrocarbons are released into the atmosphere
canister - the device in an evaporative emission control system that captures and stores
evaporative emissions generated within the vehicle for later combustion by the engine; a
canister typically contains activated carbon as a storage medium
CRC - Coordinating Research Council, a consortium of auto and oil industry members which
sponsors common research programs
diurnal cold soak - Vapor lost while vehicles are parked at ambient temperature.
HC - hydrocarbon, an organic compound consisting entirely of hydrogen and carbon; a
com- bustible fuel source which can be either gaseous or liquid
hot soak - Vapor lost in the time period immediately after turning off a vehicle.
I/M - Inspection and Maintenance program run by States to find and correct emissions
problems for vehicles registered in the State
light-duty vehicle/LDGV - passenger cars
MOVES - MOtor Vehicle Emissions Simulator; official US EPA model for estimating
emis- sions from national fleet of onroad vehicles
MSAT - Mobile Source Air Toxic rule which regulates toxic mobile source emissions such
as benzene and ethanol
permeation - the migration of hydrocarbons through materials in the fuel system
OBD - Onboard Diagnostics, an electronic automotive system with the ability to continually
track the functionality of emissions control and other components, and alerts the driver and/or
vehicle inspector when a problem is found
ORVR - Onboard refueling vapor recovery system which is designed to capture fuel vapors
at time of refueling
PSHED - portable SHED for evaporative emissions field measurements
purge - evaporative emissions control system that creates a vacuum in the fuel system to pull
the hydrocarbons from the charcoal canister while the engine is running for combustion
refueling loss - Vapor lost and spillage occurring during refueling
69
-------
running loss - Vapor lost during vehicle operation.
RVP - Reid Vapor Pressure, a measure of volatility in the gasoline at 100 degrees
Farenheit, as determined by the test method ASTM-D-323
SHED - Sealed Housing for Evaporative emissions Determination; structure for
evaporative testing in alaboratory
Stage II - vapor control programs at refueling stations to recover fuel vapor losses from fuel
displacement at the refueling pump
tank vapor generated (TVG) - vapor generated in the fuel system as temperature rises
tank vapor vented (TVV) - vapor generated in fuel system lost to the atmosphere, when
not contained by evaporative emissions control systems
Tier 2 - vehicle emissions certification standards phased in from 2004 through 2007
Tier 3 - vehicle emissions certification standards will phase in from 2017 through 2025
70
-------
References
1 USEPA (2020). Air Toxic Emissions from Onroad Vehicles in MOVES3. EPA-420-R-20-022.
Office of Transportation and Air Quality. US Environmental Protection Agency. Ann
Arbor, Ml. November 2020. https://www.epa.gov/moves/moves-onroad-technical-
re ports
2 USEPA (2023). S peciation of Total Organic Gas and Particulate Matter Emissions from Onroad
Vehicles in M0VES4. EPA-420-R-23-006. Office of Transportation and Air Quality. US
Environmental Protection Agency. Ann Arbor, Ml. August 2023.
https://www.epa.gov/moves/moves-onroad-technical-reports
3 EPA Science Inventory, MOVES2014: Evaporative Emissions Report, Record Type: DOCUMENT
(INTERNAL REPORT) Product Published Date: 09/30/2014, Record Last Revised:
02/26/2016,0MB Category: Other, Record ID: 263651
https://cfpub.epa.gov/si/si public record report.cfm?dirEntryld=263651
4 J. Heiken, T. Carlson, D. McClement and J. Lyons. CRC Project No. E-116: Assessment of MOVES
Model Evaporative Emission Inputs, 2017. http://crcsite.wpengine.com/wp-
content/uploads/2019/05/CRC E-116 MOVES Final-Report 2017-06-14.pdf
5 USEPA (2024). Population and Activity of Onroad Vehicles in M0VES5. EPA-420-R-24-019 Office of
Transportation and Air Quality. US Environmental Protection Agency. Ann Arbor, Ml.
November 2024. https://www.epa.gov/moves/moves-onroad-technical-reports
6 Harold M. Haskew and Thomas F. Liberty. CRCE-77 Vehicle Evaporative Emission
Mechanisms: A Pilot Study. June 2008. https://crcao.org/wp-content/uploads/2019/05/E-
77-Pilot-Study-Final-Report-6.24.08.pdf
7 Harold M. Haskew and Thomas F. Liberty. CRC E-77-2 Enhanced Evaporative Emission
Vehicles. March 2010. https://crcao.org/wp-content/uploads/2019/05/E-77-
2 Final Report March 2010.pdf
8 Harold M. Haskew and Thomas F. Liberty. CRC E-77-2b Evaporative Emissions From In-use
Vehicles: Test Fleet Expansion. EPA-420-R-10-025. October 2010.
https://nepis.epa.gov/Exe/ZyPDF.cgi/P100XCJC.txt?Dockev=P100XCJC.txt
9 Harold M. Haskew and Thomas F. Liberty. CRC E-77-2c Study to Determine Evaporative
Emission Breakdown, Including Permeation Effects and Diurnal Emissions Using E20 Fuels
on Aging Enhanced Evaporative Emissions Certified Vehicles. December 2010.
https://crcao.org/wp-content/uploads/2019/05/E-77-2c-Final-Report-for-sure-l-28-ll.pdf
10 Timothy H. DeFries, J. Lindner, S. Kishan, and C. Palacios. Investigation of Techniques for
High Evaporative Emissions Vehicle Detection: Denver Summer 2008 Pilot Study at
71
-------
Lipan Street Station, 2009.
https://nepis.epa.gov/Exe/ZyPURL.cgi?Dockev=P100KPX8.TXT
11 Timothy H. DeFries. Estimated Summer Hot-Soak Distributions for Denver Ken Caryl IM
Station Fleet, 2013. https://nepis.epa.gov/Exe/ZyPURL.cgi?Dockev=P100LCFS.TXT
12 CRC E-9 Measurement of Diurnal Emissions from In-Use Vehicles. September 1998.
13 Harold M. Haskew and Associates Inc. CRC E-35 Measurement of Running Loss Emissions from
In-Use Vehicles. February 1998.
14 CRC E-41-1 Real World Evaporative Testing of Late-Model In-Use Vehicles. October 1999.
15 CRC E-41-2 Evaporative Emissions from Late-Model In-Use Vehicles. October 1999.
16 Harold M. Haskew, Thomas F. Liberty, and Dennis McClement. CRC E-65 Fuel
Permeation from Automotive Systems. September 2004. https://crcao.org/wp-
content/uploads/2019/05/E-65-Final-Report-9-2-04.pdf
17 Harold M. Haskew, Thomas F. Liberty, and Dennis McClement. CRC E-65-3 Fuel
Permeation from Automotive Systems: E0, E6, E10, E20, and E85. December 2006.
https://crcao.org/wp-content/uploads/2019/05/CRC-E-65-3-Final-Report.pdf
18 J Linder, M. Sabisch, G Glinsky, J Stewart, M. St. Denis, and J. Roeschen. Multi-Day Diurnal
Testing, Eastern Research Group and SGS Environmental Testing Corporation. EPA
Contract No. EP-C-06-0-80, Work Order No.5-11, EPA-420-R-14-006, March 2014
https://nepis.epa.gov/Exe/ZyPDF.cgi/P100LC95.PDF?Dockev=P100LC95.PDF
19 M. Sabisch, S. Kishan, Inc. Eastern Research Group, J. Stewart, G. Glinsky, and SGS
Environmental Testing Corporation. Running loss testing with implanted leaks, EPA-420-R-
14009, March 2014.
https://nepis.epa.gov/Exe/ZyPDF.cgi/P100KPPC.PDF?Dockev=P100KPPC.PDF
20 EPA. Transportation And Marketing Of Petroleum Liquids. AP-42, Fifth Edition Compilation of Air
Pollutant Emissions Factors, Volume 1: Stationary Point and Area Sources. July 2008.
https://www.epa.gov/sites/default/files/2020-
09/documents/5.2 transportation and marketing of petroleum liquids.pdf
21 E.M. Liston, American Petroleum Institute, and Stanford Research Institute. A Study of
Variables that Effect the Amount of Vapor Emitted During the Refueling of Automobiles.
1975.
22 EPA Certification Data, https://www.epa.gov/compliance-and-fuel-economy-data/annual-
72
-------
certification-data-vehicles-engines-and-equipment
23 Sierra Research. Development of Trip and Soak Activity Defaults for Passenger Cars and
Trucks in MOVES2006. March 2006..
24 Sierra Research. Development of Trip and Soak Activity Defaults for Passenger Cars and
Trucks in MOVES. June 2007.
25 T. Cam, K. Cullen, S. Baldus, and K. Sime. Running Loss Temperature Profiles. 1993.
https://www.sae.org/publications/technical-papers/content/930078/
26 S.R. Reddy and Society of Automotive Engineers. Prediction of Fuel Vapor Generation from
a Vehicle Fuel Tank as a Function of Fuel Reid Vapor Pressure and Temperature.
1989.
27 J. Brown. The DELTA Model: Improved Evaporative Emissions Modeling for EPA MOVES. 2011.
https://cfpub.epa.gov/si/si public file download.cfm?downloadlD=513755
28 Wards Automotive Group. Summary Sales Data.
29 H.R. Johnson and R.S. Williams. Performance of Activated Carbon in Evaporative Loss Control
Systems. SAE Technical Paper 902119, 1990, https://doi.org/10.4271/902119.
30 Edmund's. Vehicle Fuel Tank Data, http://www.edmunds.com/.
31 USEPA. MOBILE6 Emissions Model. https://nepis.epa.gov/Exe/ZvPURL.cgi?Dockev=9100W3E2.TXT.
32 P. Caffrey and P. Machiele. In-use Volatility Impact of Commingling Ethanol and Non-Ethanol
Fuels. SAE Technical Paper 940765, 1994, https://doi.org/10.4271/940765.
33 M. Weatherby, M. Sabisch, and S. Kishan. Analysis of Evaporative On-Board Diagnostic
(OBD) Readiness and DTCs Using l/M Data. EPA Contract No. EP-C-12-017, EPA-420-R-14-
007, March 2014. https://nepis.epa.gov/Exe/ZvPDF.cgi/P100LCCC.PDF?Dockev=P100LCCC.PDF
34 Jim Kemper. High Evaporative Emitter Repairs. December 2010.
35 R. Guensler, S. Yoon, H. Li, V. Elango. Atlanta Commute Vehicle Soak and Start Distributions
and Engine Starts per Day: Impact on Mobile Source Emission Rates. EPA Contract No.
EP-05C-000127, EPA-600-R-07-075, April 2007.
https://cfpub.epa.gov/si/si public record report.cfm?Lab=NRMRL&dirEntryld=167006
36 USEPA (2024). Emission Adjustments for Onroad Vehicles in MOVES5. EPA-420-R-24-013. Office
of Transportation and Air Quality. US Environmental Protection Agency. Ann Arbor, Ml.
73
-------
November 2024. https://www.epa.gov/moves/moves-onroad-technical-reports
37 Sierra Research. United States Motor Vehicle Inspection and Maintenance (l/M) Programs.
December2005.
38 P. Cingle and D. McClement. Study of Uncontrolled Automotive Refueling Emissions. Final Report,
January 1988. https://www.osti.gov/biblio/6977672
39 G. Schremp and N. Janusch, Fuel Delivery Temperature Study, California Energy Commission.
CEC-600-2009-002-CM F, 2009.
http://citeseerx.ist. psu.edu/viewdoc/down load?doi=10.1.1.146.4412&rep=repl&type=pdf
40 "A Study of Variables that Effect the Amount of Vapor Emitted During the Refueling of
Automobiles," API Report CEA-21, May 16, 1975.
41 USEPA (2022). High Evaporative Emission Investigation Field Study. EPA-420-R-22-003. Office of
Transportation and Air Quality. US Environmental Protection Agency. Ann Arbor, Ml.
November 2022. https://www.epa.gov/moves/mobile-source-emission-factors-research.
42 USEPA (2021) Manufacturer-Run In-Use Testing Program Data for Light-Duty Vehicles, Light-Duty Trucks.
(2021) https://www.epa.gov/compliance-and-fuel-economy-data/manufacturer-run-use-
testing-program-data-light-duty-vehicles
43 USEPA (2024). Greenhouse Gas and Energy Consumption Rates for Onroad Vehicles in M0VES5.
EPA-420-R-24-018. Office of Transportation and Air Quality. US Environmental Protection
Agency. Ann Arbor, Ml. November 2024.
44 USEPA (1995). AP-42, Compilation of Air Pollutant Emissions Factors, Volume 1: Stationary Point
and Area Sources, Fifth Edition, Chapter 5, Petroleum Industry.
https://gaftp.epa.gov/ap42/ch05/s02/final/c05s02 1995.pdf
45 USEPA (2020). Population and Activity of Onroad Vehicles in M0VES3. EPA-420-R-20-023. Office
of Transportation and Air Quality. US Environmental Protection Agency. Ann Arbor, Ml.
November 2020. https://www.epa.gov/moves/moves-onroad-technical-reports
74
------- |