
Article 1 
Regulation (EC) No 692/2008 is amended as follows:

((1)) In Article 2, the following points 43 and 44 are added:
'
43.. ‘base emission strategy’ (hereinafter ‘BES’) means an emission strategy that is active throughout the speed and load operating range of the vehicle unless an auxiliary emission strategy is activated;
44.. ‘auxiliary emission strategy’ (hereinafter ‘AES’) means an emission strategy that becomes active and replaces or modifies a BES for a specific purpose and in response to a specific set of ambient or operating conditions and only remains operational as long as those conditions exist.'.
((2)) In Article 3(10) the third paragraph shall be replaced by the following text:
'Until three years after the dates specified in Article 10(4) and four years after the dates specified in Article 10(5) of Regulation (EC) No 715/2007 the following provisions shall apply:'.
((3)) Article 3(10)(a) shall be replaced by the following text:
'The requirements of point 2.1 of Annex IIIA shall not apply.'.
((4)) In Article 5, the following paragraphs 11 and 12 are inserted:
'
11.. The manufacturer shall also provide an extended documentation package with the following information:

((a)) information on the operation of all AES and BES, including a description of the parameters that are modified by any AES and the boundary conditions under which the AES operate, and indication of the AES or BES which are likely to be active under the conditions of the test procedures set out in this Regulation;
((b)) a description of the fuel system control logic, timing strategies and switch points during all modes of operation.
12.. The extended documentation package referred to in paragraph 11 shall remain strictly confidential. It may be kept by the approval authority, or, at the discretion of the approval authority, may be retained by the manufacturer. In the case the manufacturer retains the documentation package, that package shall be identified and dated by the approval authority once reviewed and approved. It shall be made available for inspection by the approval authority at the time of approval or at any time during the validity of the approval.'.
((5)) Appendix 6 to Annex I is amended as set out in Annex I to this Regulation.
((6)) Annex IIIA is amended as set out in Annex II to this Regulation.
Article 2 
This Regulation shall enter into force on the twentieth day following that of its publication in the Official Journal of the European Union.
This Regulation shall be binding in its entirety and directly applicable in all Member States.Done at Brussels, 20 April 2016.
For the Commission
The President
Jean-Claude JUNCKER
ANNEX I

In Appendix 6 to Annex I to Regulation (EC) No 692/2008, Table 1 is amended as follows:

((1)) rows ZD, ZE, ZF are replaced by the following:
'
ZD Euro 6c Euro 6-2 M, N1 class I PI, CI  1.9.2018 31.8.2019
ZE Euro 6c Euro 6-2 N1 class II PI, CI  1.9.2019 31.8.2020
ZF Euro 6c Euro 6-2 N1 class III, N2 PI, CI  1.9.2019 31.8.2020'
((2)) the following rows are inserted after row ZF:
'
ZG Euro 6d-TEMP Euro 6-2 M, N1 class I PI, CI 1.9.2017 1.9.2019 31.12.2020
ZH Euro 6d-TEMP Euro 6-2 N1 class II PI, CI 1.9.2018 1.9.2020 31.12.2021
ZI Euro 6d-TEMP Euro 6-2 N1 class III, N2 PI, CI 1.9.2018 1.9.2020 31.12.2021
ZJ Euro 6d Euro 6-2 M, N1 class I PI, CI 1.1.2020 1.1.2021 
ZK Euro 6d Euro 6-2 N1 class II PI, CI 1.1.2021 1.1.2022 
PLN Euro 6d Euro 6-2 N1 class III, N2 PI, CI 1.1.2021 1.1.2022 '
((3)) in the key to the table, the following paragraphs are inserted after the paragraph concerning the ‘Euro 6b’ emissions standard:
'“Euro 6c” emissions standardFull Euro 6 emission requirements but without quantitative RDE requirements, i.e. Euro 6b emission standard, final particle number standards for PI vehicles, use of E10 and B7 reference fuel (where applicable) assessed on regulatory lab test cycle and RDE testing for monitoring only (no NTE emission limits applied);“Euro 6d-TEMP” emissions standardFull Euro 6 emission requirements, i.e. Euro 6b emission standard, final particle number standards for PI vehicles, use of E10 and B7 reference fuel (where applicable) assessed on regulatory lab test cycle and RDE testing against temporary conformity factors;';
((4)) in the key to the table, the paragraph concerning the ‘Euro 6c’ emissions standard is replaced by the following:
'“Euro 6d” emissions standardFull Euro 6 emission requirements, i.e. Euro 6b emission standard, final particle number standards for PI vehicles, use of E10 and B7 reference fuel (where applicable) assessed on regulatory lab test cycle and RDE testing against final conformity factors;'.

ANNEX II

Annex IIIA to Regulation (EC) No 692/2008 is amended as follows:

((1)) point 2.1 is replaced by the following:
'
2.1 
Throughout the normal life of a vehicle type approved according to Regulation (EC) No 715/2007, its emissions determined in accordance with the requirements of this Annex and emitted at any possible RDE test performed in accordance with the requirements of this Annex, shall not be higher than the following not-to-exceed (NTE) values:

NTEpollutant = CFpollutant × TF(p1,…, pn) × EURO-6

where EURO-6 is the applicable Euro 6 emission limit laid down in Table 2 of Annex I to Regulation (EC) No 715/2007.
';
((2)) the following points 2.1.1, 2.1.2 and 2.1.3 are inserted:
'
2.1.1 
The conformity factor CFpollutant for the respective pollutant is specified as follows:


Pollutant Mass of oxides of nitrogen (NOx) Number of particles (PN) Mass of carbon monoxide (CO) Mass of total hydrocarbons (THC) Combined mass of total hydrocarbons and oxides of nitrogen (THC + NOx)
CFpollutant 1 + margin with margin = 0,5 to be determined — — —


“margin” is a parameter taking into account the additional measurement uncertainties introduced by the PEMS equipment, which are subject to an annual review and shall be revised as a result of the improved quality of the PEMS procedure or technical progress.

2.1.2 
By way of exception to the provisions of point 2.1.1, during a period of 5 years and 4 months following the dates specified in Article 10(4) and (5) of Regulation (EC) No 715/2007 and upon request of the manufacturer, the following temporary conformity factors may apply:


Pollutant Mass of oxides of nitrogen (NOx) Number of particles (PN) Mass of carbon monoxide (CO) Mass of total hydrocarbons (THC) Combined mass of total hydrocarbons and oxides of nitrogen (THC + NOx)
CFpollutant 2,1 to be determined — — —


The application of temporary conformity factors shall be recorded in the certificate of conformity of the vehicle.

2.1.3 
The transfer function TF(p1,…, pn) referred to in point 2.1 is set to 1 for the entire range of parameters pi (i = 1,…,n).

If the transfer function TF(p1,…, pn) is amended, this shall be done in a manner which is not detrimental to the environmental impact and the effectiveness of the RDE test procedures. In particular the following condition shall hold:

∫ TF (p1,…, pn) * Q (p1,…, pn) dp = ∫ Q (p1,…, pn) dp

Where:


— dp represents the integral over the entire space of the parameters pi (i = 1,…,n)
— Q(p1,…, pn), is the probability density of an event corresponding to the parameters pi (i = 1,…,n) in real driving.
';
((3)) the following point 3.1.0 is inserted:
'
3.1.0 

3.1.0.1 Mgas,d,t ≤ NTEpollutant and Mgas,d,u ≤ NTEpollutant with the definitions of point 2.1 of this Annex and points 6.1 and 6.3 of Appendix 5 and the setting gas = pollutant.
3.1.0.2 Mw,gas,d ≤ NTEpollutant and Mw,gas,d,U ≤ NTEpollutant with the definitions of point 2.1 of this Annex and point 3.9 of Appendix 6 and the setting gas = pollutant.
';
((4)) point 5.3 is deleted;
((5)) point 5.4 is replaced by the following:
'
5.4. 
The dynamic conditions encompass the effect of road grade, head wind and driving dynamics (accelerations, decelerations) and auxiliary systems upon energy consumption and emissions of the test vehicle. The verification of the normality of dynamic conditions shall be done after the test is completed, using the recorded PEMS data. This verification shall be conducted in two steps:


5.4.1 The overall excess or insufficiency of driving dynamics during the trip shall be checked using the methods described in Appendix 7a to this Annex.
5.4.2 If the trip results as valid following the verifications according to point 5.4.1, the methods for verifying the normality of the dynamic conditions and laid down in Appendices 5 and 6 to this Annex must be applied. Each method includes a reference for dynamic conditions, ranges around the reference and the minimum coverage requirements to achieve a valid test.
';
((6)) point 6.8 is replaced by the following:
'
6.8 The average speed (including stops) of the urban driving part of the trip should be between 15 and 40 km/h. Stop periods, defined as vehicle speed of less than 1 km/h, shall account for 6-30 % of the time duration of urban operation. Urban operation shall contain several stop periods of 10 s or longer. If a stop period lasts more than 180 s, the emission events during the 180 s following such an excessively long stop period shall be excluded from the evaluation.
';
((7)) in point 6.11, the following sentence is added:
'In addition, the proportional cumulative positive altitude gain shall be less than 1 200 m/100km) and be determined according to Appendix 7b.';
((8)) point 9.5 is replaced by the following:
'
9.5. If during a particular time interval the ambient conditions are extended in accordance with point 5.2, the emissions during this particular time interval, calculated according to Appendix 4, shall be divided by a value of 1,6 before being evaluated for compliance with the requirements of this Annex.
';
((9)) Appendix 1 is amended as follows:

((a)) in point 3.4.6, the following sentence is added:
'It is permitted to power any safety-related illumination of fixtures and installations of PEMS components outside of the vehicle's cabin by the vehicle's battery.';
((b)) in point 4.5, the following sentence is added:
'To minimise analyser drift, one should conduct the zero and span calibration of analysers at an ambient temperature that resembles, as closely as possible, the temperature experienced by the test equipment during the RDE trip.';
((10)) in Appendix 2, footnote 2 to Table 4 in point 8 is replaced by the following:
'
(2) This general requirement applies to the speed sensor only; if vehicle speed is used to determine parameters like acceleration, the product of speed and positive acceleration, or RPA, the speed signal shall have an accuracy of 0,1 % above 3 km/h and a sampling frequency of 1 Hz. This accuracy requirement can be met by using the signal of a wheel rotational speed sensor.
';
((11)) in Appendix 6 point 2 the following definition is deleted:
'aiActual acceleration in time step i, if not other defined in an equation:
ai=vi+1−vi3,6×ti+1−ti, [m/s2]';
((12)) in Appendix 6 point 2 the following definitions are inserted:
'm–gas,UWeighted emission value of an exhaust gas component ‘gas’ for the subsample of all seconds i with vi < 60 km/h, g/sMw,gas,d,UWeighted distance-specific emissions for the exhaust gas component ‘gas’ for the subsample of all seconds i with vi < 60 km/h, g/kmv–UWeighted vehicle speed in the wheel power class j, km/h';
((13)) in Appendix 6 point 3.1 the first paragraph is replaced by the following text:
'The actual wheel power Pr,i shall be the total power to overcome air resistance, rolling resistance, road gradients, longitudinal inertia of the vehicle and rotational inertia of the wheels.';
((14)) in Appendix 6 point 3.2 is replaced by the following text:
'
3.2 
The standard power frequencies are defined for urban driving and for the total trip (see paragraph 3.4) and a separate evaluation of the emissions shall be made for the total trip and for the urban part. The three second moving averages calculated according to paragraph 3.3 shall therefore be allocated later to urban and extra-urban driving conditions according to the velocity signal (vi) from the actual second i as outlined in Table 1-1.


 Urban Rural Motorway
vi [km/h] 0 to ≤ 60 > 60 to ≤ 90 > 90
'
((15)) in Appendix 6 point 3.9 is replaced by the following text:
'
3.9. 
The time-based weighted averages of the emissions in the test shall be converted into distance-based emissions once for the urban data set and once for the total data set as follows:


 For the total trip:Mw,gas,d=1000×m–gas×3600v–
 For the urban part of the trip:Mw,gas,d,U=1000×m–gas,U×3600v–U

Using these formulas, weighted averages shall be calculated for the following pollutants for the total trip and for the urban part of the trip:

Mw,NOx,dweighted NOx test result in [mg/km]Mw,NOx,d,Uweighted NOx test result in [mg/km]Mw,CO,dweighted CO test result in [mg/km]Mw,CO,d,Uweighted CO test result in [mg/km]
';
((16)) the following Appendices 7a and 7b are inserted:

Appendix 7a  1. 
This Appendix describes the calculation procedures to verify the overall trip dynamics, to determine the overall excess or absence of dynamics during urban, rural and motorway driving.
 2. 
RPArelative positive acceleration‘acceleration resolution ares’minimum acceleration > 0 measured in m/s2T4253Hcompound data smoother‘positive acceleration apos’acceleration [m/s2] greater than 0,1 m/s2

Index (i) refers to the time step

Index (j) refers to the time step of positive acceleration datasets

Index (k) refers to the category (t = total, u = urban, r = rural, m = motorway)

Δdifference>larger≥larger or equal%per cent<smaller≤smaller or equalaacceleration [m/s2]aiacceleration in time step i [m/s2]apospositive acceleration greater than 0,1 m/s2 [m/s2]apos,i,kpositive acceleration greater than 0,1 m/s2 in time step i considering the urban, rural and motorway shares [m/s2]aresacceleration resolution [m/s2]didistance covered in time step i [m]di,kdistance covered in time step i considering the urban, rural and motorway shares [m]Mknumber of samples for urban, rural and motorway shares with positive acceleration greater than 0,1 m/s2Nktotal number of samples for the urban, rural and motorway shares and the complete tripRPAkrelative positive acceleration for urban, rural and motorway shares [m/s2 or kWs/(kg × km)]tkduration of the urban, rural and motorway shares and the complete trip [s]vvehicle speed [km/h]viactual vehicle speed in time step i [km/h]vi,kactual vehicle speed in time step i considering the urban, rural and motorway shares [km/h](v · a)iactual vehicle speed per acceleration in time step i [m2/s3 or W/kg](v · apos)j,kactual vehicle speed per positive acceleration greater than 0,1 m/s2 in time step j considering the urban, rural and motorway shares [m2/s3 or W/kg].(v · apos)k_[95]95th percentile of the product of vehicle speed per positive acceleration greater than 0,1 m/s2 for urban, rural and motorway shares [m2/s3 or W/kg]v–kaverage vehicle speed for urban, rural and motorway shares [km/h]
 3.  3.1.  3.1.1. 
Dynamic parameters like acceleration, v · apos or RPA shall be determined with a speed signal of an accuracy of 0,1 % above 3 km/h and a sampling frequency of 1 Hz. This accuracy requirement is generally fulfilled by wheel (rotational) speed signals.

The speed trace shall be checked for faulty or implausible sections. The vehicle speed trace of such sections is characterised by steps, jumps, terraced speed traces or missing values. Short faulty sections shall be corrected, for example by data interpolation or benchmarking against a secondary speed signal. Alternatively, short trips containing faulty sections could be excluded from the subsequent data analysis. In a second step the acceleration values shall be ranked in ascending order, in order to determine the acceleration resolution ares = (minimum acceleration value > 0).

If ares ≤ 0,01 m/s2, the vehicle speed measurement is accurate enough.

If 0,01 < ares ≤ rmax m/s2, smoothing by using a T4253 Hanning filter.

If ares > rmax m/s2, the trip is invalid.

The T4253 Hanning filter performs the following calculations: The smoother starts with a running median of 4, which is centred by a running median of 2. It then re-smoothes these values by applying a running median of 5, a running median of 3, and Hanning (running weighted averages). Residuals are computed by subtracting the smoothed series from the original series. This whole process is then repeated on the computed residuals. Finally, the smoothed residuals are computed by subtracting the smoothed values obtained the first time through the process.

The correct speed trace builds the basis for further calculations and binning as described in paragraph 3.1.2.
 3.1.2. 
The following calculations shall be performed over the whole time-based speed trace (1 Hz resolution) from second 1 to second tt (last second).

The distance increment per data sample shall be calculated as follows:

di = vi/3,6, i = 1 toNt

Where:


 di is the distance covered in time step i [m]
 vi is the actual vehicle speed in time step i [km/h]
 Nt is the total number of samples

The acceleration shall be calculated as follows:

ai = (vi + 1 – vi – 1)/(2 · 3,6), i = 1 toNt

Where:

ai is the acceleration in time step i [m/s2]. For i = 1: vi – 1 = 0, for i = Nt: vi + 1 = 0.

The product of vehicle speed per acceleration shall be calculated as follows:

(v · a)i = vi · ai/3,6, i = 1 toNt

Where:

(v · a)i is the product of the actual vehicle speed per acceleration in time step i [m2/s3 or W/kg].
 3.1.3. 
After the calculation of ai and (v · a)i, the values vi, di, ai and (v · a)i shall be ranked in ascending order of the vehicle speed.

All datasets with vi ≤ 60 km/h belong to the ‘urban’ speed bin, all datasets with 60 km/h < vi ≤ 90 km/h belong to the ‘rural’ speed bin and all datasets with vi > 90 km/h belong to the ‘motorway’ speed bin.

The number of datasets with acceleration values ai > 0,1 m/s2 shall be bigger or equal to 150 in each speed bin.

For each speed bin the average vehicle speed v–k shall be calculated as follows:

v–k=∑ivi,k∕Nk, i = 1 toNk,k = u,r,m

Where:

Nk is the total number of samples of the urban, rural, and motorway shares.
 3.1.4. 
The 95th percentile of the v · apos values shall be calculated as follows:

The (v · a)i,k values in each speed bin shall be ranked in ascending order for all datasets with ai,k ≥ 0,1 m/s2 and the total number of these samples Mk shall be determined.

Percentile values are then assigned to the (v · apos)j,k values with ai,k ≥ 0,1 m/s2 as follows:

The lowest v · apos value gets the percentile 1/Mk, the second lowest 2/Mk, the third lowest 3/Mk and the highest value Mk/Mk =100 %.

(v · apos)k_[95] is the (v · apos)j,k value, with j/Mk = 95 %. If j/Mk = 95 % cannot be met, (v · apos)k_[95] shall be calculated by linear interpolation between consecutive samples j and j + 1 with j/Mk < 95 % and (j + 1)/Mk > 95 %.

The relative positive acceleration per speed bin shall be calculated as follows:

RPAk = Σj(Δt · (v · apos)j,k)/Σidi,k, j = 1 toMk,i = 1 toNk,k = u,r,m

Where:

RPAk is the relative positive acceleration for urban, rural and motorway shares in [m/s2 or kWs/(kg*km)]

Δttime difference equal to 1 secondMkthe sample number for urban, rural and motorway shares with positive accelerationNkthe total sample number for urban, rural and motorway shares.
 4.  4.1.1. 
If v–k≤74,6 km∕h

and
v×aposk_95>0,136×v–k+14,44
is fulfilled, the trip is invalid.

If v–k>74,6 km∕h and v×aposk_95>0,0742×v–k+18,966 is fulfilled, the trip is invalid.
 4.1.2. 
If v–k≤94,05 km∕h and RPAk<−0,0016×v–k+0,1755 is fulfilled, the trip is invalid.

If v–k>94,05 km∕h and RPAk < 0,025 is fulfilled, the trip is invalid.

Appendix 7b  1. 
This Appendix describes the procedure to determine the cumulative elevation gain of an RDE trip.
 2. 
d(0)distance at the start of a trip [m]dcumulative distance travelled at the discrete way point under consideration [m]d0cumulative distance travelled until the measurement directly before the respective way point d [m]d1cumulative distance travelled until the measurement directly after the respective way point d [m]dareference way point at d(0) [m]decumulative distance travelled until the last discrete way point [m]diinstantaneous distance [m]dtottotal test distance [m]h(0)vehicle altitude after the screening and principle verification of data quality at the start of a trip [m above sea level]h(t)vehicle altitude after the screening and principle verification of data quality at point t [m above sea level]h(d)vehicle altitude at the way point d [m above sea level]h(t-1)vehicle altitude after the screening and principle verification of data quality at point t-1 [m above sea level]hcorr(0)corrected altitude directly before the respective way point d [m above sea level]hcorr(1)corrected altitude directly after the respective way point d [m above sea level]hcorr(t)corrected instantaneous vehicle altitude at data point t [m above sea level]hcorr(t-1)corrected instantaneous vehicle altitude at data point t-1 [m above sea level]hGPS,iinstantaneous vehicle altitude measured with GPS [m above sea level]hGPS(t)vehicle altitude measured with GPS at data point t [m above sea level]hint(d)interpolated altitude at the discrete way point under consideration d [m above sea level]hint,sm,1(d)smoothed interpolated altitude, after the first smoothing run at the discrete way point under consideration d [m above sea level]hmap(t)vehicle altitude based on topographic map at data point t [m above sea level]Hzhertzkm/hkilometre per hourmmetreroadgrade,1(d)smoothed road grade at the discrete way point under consideration d after the first smoothing run [m/m]roadgrade,2(d)smoothed road grade at the discrete way point under consideration d after the second smoothing run [m/m]sintrigonometric sine functionttime passed since test start [s]t0time passed at the measurement directly located before the respective way point d [s]viinstantaneous vehicle speed [km/h]v(t)vehicle speed of data point t [km/h].
 3. 
The cumulative positive elevation gain of an RDE trip shall be determined based on three parameters: the instantaneous vehicle altitude hGPS,i [m above sea level] as measured with the GPS, the instantaneous vehicle speed vi [km/h] recorded at a frequency of 1 Hz and the corresponding time t [s] that has passed since test start.
 4.  4.1. 
The cumulative positive elevation gain of an RDE trip shall be calculated as a three-step procedure, consisting of: (i) the screening and principle verification of data quality; (ii) the correction of instantaneous vehicle altitude data; and (iii) the calculation of the cumulative positive elevation gain.
 4.2. 
The instantaneous vehicle speed data shall be checked for completeness. Correcting for missing data is permitted if gaps remain within the requirements specified in Point 7 of Appendix 4; else, the test results shall be voided. The instantaneous altitude data shall be checked for completeness. Data gaps shall be completed by data interpolation. The correctness of interpolated data shall be verified by a topographic map. It is recommended to correct interpolated data if the following condition applies:

|hGPS(t) – hmap(t)| > 40 m

The altitude correction shall be applied so that:

h(t) = hmap(t)

where:

h(t)vehicle altitude after the screening and principle verification of data quality at data point t [m above sea level]hGPS(t)vehicle altitude measured with GPS at data point t [m above sea level]hmap(t)vehicle altitude based on topographic map at data point t [m above sea level].
 4.3. 
The altitude h(0) at the start of a trip at d(0) shall be obtained by GPS and verified for correctness with information from a topographic map. The deviation shall not be larger than 40 m. Any instantaneous altitude data h(t) shall be corrected if the following condition applies:

|h(t) – h(t – 1)| > (v(t)/3,6 * sin45°)

The altitude correction shall be applied so that:

hcorr(t) = hcorr(t-1)

where:

h(t)vehicle altitude after the screening and principle verification of data quality at data point t [m above sea level]h(t-1)vehicle altitude after the screening and principle verification of data quality at data point t-1 [m above sea level]v(t)vehicle speed of data point t [km/h]hcorr(t)corrected instantaneous vehicle altitude at data point t [m above sea level]hcorr(t-1)corrected instantaneous vehicle altitude at data point t-1 [m above sea level].

Upon the completion of the correction procedure, a valid set of altitude data is established. This data set shall be used for the final calculation of the cumulative positive elevation gain as described in point 4.4.
 4.4.  4.4.1. 
The total distance dtot [m] covered by a trip shall be determined as sum of the instantaneous distances di. The instantaneous distance di shall be determined as:
di=vi3,6
Where:

diinstantaneous distance [m]viinstantaneous vehicle speed [km/h]

The cumulative elevation gain shall be calculated from data of a constant spatial resolution of 1 m starting with the first measurement at the start of a trip d(0). The discrete data points at a resolution of 1 m are referred to as way points, characterised by a specific distance value d (e.g., 0, 1, 2, 3 m…) and their corresponding altitude h(d) [m above sea level].

The altitude of each discrete way point d shall be calculated through interpolation of the instantaneous altitude hcorr(t) as:
hintd=hcorr0+hcorr1−hcorr0d1−d0×d−d0
Where:

hint(d)interpolated altitude at the discrete way point under consideration d [m above sea level]hcorr(0)corrected altitude directly before the respective way point d [m above sea level]hcorr(1)corrected altitude directly after the respective way point d [m above sea level]dcumulative distance travelled until the discrete way point under consideration d [m]d0cumulative distance travelled until the measurement located directly before the respective way point d [m]d1cumulative distance travelled until the measurement located directly after the respective way point d [m].
 4.4.2. 
The altitude data obtained for each discrete way point shall be smoothed by applying a two-step procedure; da and de denote the first and last data point respectively (Figure 1). The first smoothing run shall be applied as follows:

roadgrade,1d=hintd+200 m−hintdad+200 mfor d ≤ 200 m

roadgrade,1d=hintd+200 m−hintd−200 md+200 m−d−200 mfor 200 m < d < (de – 200 m)

roadgrade,1d=hintde−hintd−200 mde−d−200 mfor d ≥ (de – 200 m)

hint,sm,1(d) = hint,sm,1(d – 1 m) + roadgrade,1(d), d = da + 1 to de

hint,sm,1(da) = hint(da) + roadgrade,1(da)

Where:

roadgrade,1(d)smoothed road grade at the discrete way point under consideration after the first smoothing run [m/m]hint(d)interpolated altitude at the discrete way point under consideration d [m above sea level]hint,sm,1(d)smoothed interpolated altitude, after the first smoothing run at the discrete way point under consideration d [m above sea level]dcumulative distance travelled at the discrete way point under consideration [m]dareference way point at a distance of zero metres [m]decumulative distance travelled until the last discrete way point [m].

The second smoothing run shall be applied as follows:

roadgrade,2d=hint,sm,1d+200 m−hint,sm,1dad+200 mfor d ≤ 200 m

roadgrade,2d=hint,sm,1d+200 m−hint,sm,1d−200 md+200 m−d−200 mfor 200 m < d < (de – 200 m)

roadgrade,2d=hint,sm,1de−hint,sm,1d−200 mde−d−200 mfor d ≥ (de – 200 m)

Where:

roadgrade,2(d)smoothed road grade at the discrete way point under consideration after the second smoothing run [m/m]hint,sm,1(d)smoothed interpolated altitude, after the first smoothing run at the discrete way point under consideration d [m above sea level]dcumulative distance travelled at the discrete way point under consideration [m]dareference way point at a distance of zero metres [m]decumulative distance travelled until the last discrete way point [m].

Figure 1 4.4.3. 
The positive cumulative elevation gain of a trip shall be calculated by integrating all positive interpolated and smoothed road grades, i.e. roadgrade,2(d). The result should be normalised by the total test distance dtot and expressed in meters of cumulative elevation gain per 100 kilometres of distance.
 5. 
Tables 1 and 2 show the steps performed in order to calculate the positive elevation gain on the basis of data recorded during an on-road test performed with PEMS. For the sake of brevity an extract of 800 m and 160 s is presented here.
 5.1. 
The screening and principle verification of data quality consists of two steps. First, the completeness of vehicle speed data is checked. No data gaps related to vehicle speed are detected in the present data sample (see Table 1). Second, the altitude data are checked for completeness; in the data sample, altitude data related to seconds 2 and 3 are missing. The gaps are filled by interpolating the GPS signal. In addition, the GPS altitude is verified by a topographic map; this verification includes the altitude h(0) at the start of the trip. Altitude data related to seconds 112-114 are corrected on the basis of the topographic map to satisfy the following condition:

hGPS(t) – hmap(t) < – 40 m

As result of the applied data verification, the data in the fifth column h(t) are obtained.
 5.2. 
As a next step, the altitude data h(t) of seconds 1 to 4, 111 to 112 and 159 to 160 are corrected assuming the altitude values of seconds 0, 110 and 158 respectively since the following condition applies:

|h(t) – h(t – 1)| > (v(t)/3,6 * sin45°)

As result of the applied data correction, the data in the sixth column hcorr(t) are obtained. The effect of the applied verification and correction steps on the altitude data is depicted in Figure 2.
 5.3.  5.3.1. 
The instantaneous distance di is calculated by dividing the instantaneous vehicle speed measured in km/h by 3,6 (Column 7 in Table 1). Recalculating the altitude data to obtain a uniform spatial resolution of 1 m yields the discrete way points d (Column 1 in Table 2) and their corresponding altitude values hint(d) (Column 7 in Table 2). The altitude of each discrete way point d is calculated through interpolation of the measured instantaneous altitude hcorr as:
hint0=120,3+120,3−120,30,1−0,0×0−0=120,3000hint520=132,5+132,6−132,5523,6−519,9×520−519,9=132,5027 5.3.2. 
In Table 2, the first and last discrete way points are: da = 0 m and de = 799 m, respectively. The altitude data of each discrete way point is smoothed by applying a two-step procedure. The first smoothing run consists of:
roadgrade,10=hint200 m−hint00+200 m=120,9682−120,3000200=0,0033
chosen to demonstrate the smoothing for d ≤ 200 m
roadgrade,1320=hint520−hint120520−120=132,5027−121,9808400=0,0288
chosen to demonstrate the smoothing for 200 m < d < (599 m)
roadgrade,1720=hint799−hint520799−520=121,2000−132,5027279=−0,0405
chosen to demonstrate the smoothing for d ≥ (599 m)

The smoothed and interpolated altitude is calculated as:

hint,sm,1(0) = hint(0) + roadgrade,1(0) = 120,3 + 0,0033 ≈ 120,3033 m

hint,sm,1(799) = hint,sm,1(798) + roadgrade,1(799) = 121,2550 – 0,0220 = 121,2330 m

Second smoothing run:
roadgrade,20=hint,sm,1200−hint,sm,10200=119,9618−120,3033200=−0,0017
chosen to demonstrate the smoothing for d ≤ 200 m
roadgrade,2320=hint,sm,1520−hint,sm,1120520−120=123,6809−120,1843400=0,0087
chosen to demonstrate the smoothing for 200 m < d < (599 m)
roadgrade,2720=hint,sm,1799−hint,sm,1520799−520=121,2330−123,6809279=−0,0088
chosen to demonstrate the smoothing for d ≥ (599 m)
 5.3.3. 
The positive cumulative elevation gain of a trip is calculated by integrating all positive interpolated and smoothed road grades, i.e. roadgrade,2(d). For the presented example the total covered distance was dtot = 139,7 km and all positive interpolated and smoothed road grades were of 516 m. Therefore a positive cumulative elevation gain of 516 × 100/139,7 = 370 m/100 km was achieved.


Time t [s] v(t)[km/h] hGPS(t)[m] hmap(t)[m] h(t)[m] hcorr(t)[m] di[m] Cum. d[m]
0 0,00 122,7 129,0 122,7 122,7 0,0 0,0
1 0,00 122,8 129,0 122,8 122,7 0,0 0,0
2 0,00 - 129,1 123,6 122,7 0,0 0,0
3 0,00 - 129,2 124,3 122,7 0,0 0,0
4 0,00 125,1 129,0 125,1 122,7 0,0 0,0
… … … … … … … …
18 0,00 120,2 129,4 120,2 120,2 0,0 0,0
19 0,32 120,2 129,4 120,2 120,2 0,1 0,1
… … … … … … … …
37 24,31 120,9 132,7 120,9 120,9 6,8 117,9
38 28,18 121,2 133,0 121,2 121,2 7,8 125,7
… … … … … … … …
46 13,52 121,4 131,9 121,4 121,4 3,8 193,4
47 38,48 120,7 131,5 120,7 120,7 10,7 204,1
… … … … … … … …
56 42,67 119,8 125,2 119,8 119,8 11,9 308,4
57 41,70 119,7 124,8 119,7 119,7 11,6 320,0
… … … … … … … …
110 10,95 125,2 132,2 125,2 125,2 3,0 509,0
111 11,75 100,8 132,3 100,8 125,2 3,3 512,2
112 13,52 0,0 132,4 132,4 125,2 3,8 516,0
113 14,01 0,0 132,5 132,5 132,5 3,9 519,9
114 13,36 24,30 132,6 132,6 132,6 3,7 523,6
… … … … … … … 
149 39,93 123,6 129,6 123,6 123,6 11,1 719,2
150 39,61 123,4 129,5 123,4 123,4 11,0 730,2
… … … … … … … 
157 14,81 121,3 126,1 121,3 121,3 4,1 792,1
158 14,19 121,2 126,2 121,2 121,2 3,9 796,1
159 10,00 128,5 126,1 128,5 121,2 2,8 798,8
160 4,10 130,6 126,0 130,6 121,2 1,2 800,0
- denotes data gaps


d[m] t0[s] d0[m] d1[m] h0[m] h1[m] hint(d)[m] roadgrade,1(d)[m/m] hint,sm,1(d)[m] roadgrade,2(d)[m/m]
0 18 0,0 0,1 120,3 120,4 120,3 0,0035 120,3 – 0,0015
… … … … … … … … … …
120 37 117,9 125,7 120,9 121,2 121,0 – 0,0019 120,2 0,0035
… … … … … … … … … …
200 46 193,4 204,1 121,4 120,7 121,0 – 0,0040 120,0 0,0051
… … … … … … … … … …
320 56 308,4 320,0 119,8 119,7 119,7 0,0288 121,4 0,0088
… … … … … … … … … …
520 113 519,9 523,6 132,5 132,6 132,5 0,0097 123,7 0,0037
… … … … … … … … … …
720 149 719,2 730,2 123,6 123,4 123,6 – 0,0405 122,9 – 0,0086
… … … … … … … … … …
798 158 796,1 798,8 121,2 121,2 121,2 – 0,0219 121,3 – 0,0151
799 159 798,8 800,0 121,2 121,2 121,2 – 0,0220 121,3 – 0,0152

Figure 2
Figure 3

d[m] t0[s] d0[m] d1[m] h0[m] h1[m] hint(d)[m] roadgrade,1(d)[m/m] hint,sm,1(d)[m] roadgrade,2(d)[m/m]
0 18 0,0 0,1 120,3 120,4 120,3 0,0035 120,3 – 0,0015
… … … … … … … … … …
120 37 117,9 125,7 120,9 121,2 121,0 – 0,0019 120,2 0,0035
… … … … … … … … … …
200 46 193,4 204,1 121,4 120,7 121,0 – 0,0040 120,0 0,0051
… … … … … … … … … …
320 56 308,4 320,0 119,8 119,7 119,7 0,0288 121,4 0,0088
… … … … … … … … … …
520 113 519,9 523,6 132,5 132,6 132,5 0,0097 123,7 0,0037
… … … … … … … … … …
720 149 719,2 730,2 123,6 123,4 123,6 – 0,0405 122,9 – 0,0086
… … … … … … … … … …
798 158 796,1 798,8 121,2 121,2 121,2 – 0,0219 121,3 – 0,0151
799 159 798,8 800,0 121,2 121,2 121,2 – 0,0220 121,3 – 0,0152

