Abstract

Real-time performance measures are important for agencies to maintain their roadways during the winter season. However, infrastructure can be expensive to deploy and maintain and may be sparse in rural areas, while speed data alone may not provide enough fidelity in borderline conditions. This study looks at high-frequency vehicle data from the in-vehicle bus to detect changes in the vehicle and driver behavior during changing winter road conditions. The data is reported to the cloud via cellular communication and is viewable in real-time using a map-based web dashboard. Three winter weather events are assessed using in-vehicle data collected from the February–March 2018 period. Using a two-sample Kolmogorov-Smirnov test, we found high significance in the reductions of applied brake pressures and rates of braking in winter versus fair weather conditions before vehicle intervention is necessary. The paper concludes that pairwise comparison of driver change in brake pressure may be a valuable data source indicative of deteriorating road conditions before more severe indicators such as traction control, antilock brake, and/or hazard indicators are activated.

Introduction

Winter road conditions are hazardous for drivers to navigate and are a challenge for agencies to manage. According to the US Federal Highway Administration’s (FHWA’s) Road Weather Management Program (RWMP), out of approximately 5,748,000 vehicle crashes each year, about 22% of these crashes are weather related (FHWA 2018). From 2005 to 2014, an average of 445,303 persons were injured and 5,897 persons were killed annually in weather-related crashes. Of these weather-related crashes, 17% of occurred during snow or sleet, 13% occurred on icy pavement, and 14% occurred on snowy or slushy pavement. The estimated annual cost of weather-related crashes ranges from $22 to $51 billion (Pisano et al. 2008). Additionally, it is estimated that weather-related delays cost the trucking industry about $8 to $9 billion annually (FHWA 2017a).
In 2015, the Road Weather Management Performance Measures report (FHWA 2015) evaluated a set of 27 performance measures aimed to assess the status of the transportation system during inclement winter weather conditions. The objectives identified in an updated report (FHWA 2017b) included “advancement of the state of the art for mobile sensing and integrating vehicle data into road weather applications,” and “improve integration of weather-related decision-support technologies into traffic operations and maintenance procedures.”

Vehicle Data

Before the advent of onboard vehicle communication electronics in the late 1980s (McMillan 1984; Cox 1984), it was difficult to obtain any type of vehicle telematics data without dedicated instrumentation. However, modern vehicles have come equipped with localized networks that facilitate communication between the various onboard components (Kiencke et al. 1986), and report vehicle dynamics information such as wheel speed, engine controls, drivetrain, and brake pressures at subsecond frequencies. The data are typically used for controlling and coordinating vehicle features such as traction control and antilock braking systems [Gerstenmeier and Weber 1987; Antilock Braking System, “Teldix,” UK Patent No. 1,245,620 (1971)]. However, if the same data were to be packaged and delivered via a system that can integrate with data-driven software interfaces to generate weather-related performance measures using modern-day cellular networks, it can provide tremendous value to stakeholders and agencies. A previous study has shown that using segment-based connected vehicle data, speed drops can be used to detect inclement weather conditions and recommend winter maintenance strategies (Hainen et al. 2012). Expanding on this concept, because modern production vehicles can be used to detect loss of friction conditions using signals from integrated onboard components (Zhao et al. 2014; Casselgren 2018), and if the data can be harvested and calibrated using high-fidelity data logging instruments and crowdsourcing, may have potential to achieve the objectives outlined by the 2015 RWMP (FHWA 2015).

Methodology

This study developed methods to extract on board vehicle data to assess whether current vehicle technologies can provide informative road condition data during winter storms.

Data System

The experimental setup (Fig. 1) was comprised of a series production vehicle (2017 Audi A4 Allroad) and a laptop. Each test was performed with one occupant (the driver), and no additional loading. The brake pad and rotors were in good condition [under 16,000 km (10,000 mi)] and tire pressures were checked before each test. The vehicle comes with a build-in wireless internet access point. A laptop was placed in the trunk with an inexpensive controller area network universal serial bus (CANUSB) device, which is a device that feeds the controller area network (CAN) signals forwarded from the various communication buses in the vehicle to the laptop, and is natively supported in Linux as a socket (Hartkopp 2012).
Fig. 1. System architecture for collecting enhanced probe data.
The laptop was used to decode and record the following CAN signals: antilock brake (ABS) activation, traction-control intervention, wheel ticks, hazard lights, windshield wiper state, and brake pressure. Depending on the bus on which the device is emitting, the frequency of the signal can be anywhere between 100 and 200 milliseconds (ms). Once a notable status change of the identified signal was detected, the signal was uploaded with a time stamp, together with a global positioning system (GPS) point where the change occurred, to a secure server via the onboard internet access point. The message was a hypertext transfer protocol (HTTP) post (Fielding 1996) and encoded as a JavaScript object notation (JSON) message (Bray 2014) sent through the cloud. Through a web interface, the server then processed the message and stored the signal into an structured query language (SQL) database. The laptop also had the ability to store the selected in-vehicle signals as comma-separated variable (CSV) files instead of uploading to the cloud.
A web dashboard was developed to verify the uploaded vehicle signals in the database. Fig. 2(a) shows a dashboard interface on a Google Maps graphical layer. The windshield wiper (Callout i), hazard light activation (Callouts ii and iii), and traction control (Callout iv) signals from the test vehicle were plotted as they were received while traveling eastbound on I-70 during a winter storm on February 5, 2018 in central Missouri. According to local reports, heavy snowfall and icing on the roadways caused over 650 crashes [Fig. 2(b)] and multiple fatalities. The interface also has the ability to look at historic events.
Fig. 2. Inclement weather event on I-70 in Missouri, February 5, 2018: (a) web dashboard showing real-time sensor data (Map data: © 2018 Google Imagery; © 2018 TerraMetrics); and (b) image from the storm. (Image by authors.)

User Road Condition Rating

A dashboard camera captured images from the forward view of the vehicle at 5-s intervals and the road condition was then assessed based on a qualitative rating scale by the driver. The ratings are as follows: a rating of 1 is dry or slightly damp where the driver would expect good friction; a rating of 2 is wet and the driver expects slightly slippery conditions; and pa rating of 3 is where there is snow or slush and where the driver would expect some sliding with a hard maneuver. Because the driver was not able to record the observation during the drives, this was performed as an after-action assessment. User ratings were recorded using software that displays an image from the experiments at random and records the user rating input. Over 9,000 images from the experiments were assessed.

Mobile Road Weather Information Sensor

Friction data was collected using a mobile road weather information sensor (MARWIS) system. The sensor was mounted on the rear of a chase vehicle about 1 m (3 ft) above the pavement surface. The coefficient of friction (μ) measured from the road surface condition was averaged and recorded at 10-s intervals. The friction estimates and based on lasers have relatively low fidelity, but these devices are commonly used by highway agencies as mobile probes during winter storms. The geographic location of the friction data was also saved using the associated GPS reported from the MARWIS device.

Wheel Ticks

Modern vehicles are equipped with wheel speed sensors (Hartley 1974) that report the position of each wheel as an integer, such as in the range of 0 to 1,000, to the various onboard control units and devices via the in-vehicle bus [T. Svantesson, “Detection of short term irregularities in a road surface,” US Patent No. 14,913,684 (2015)]. As the vehicle moves, the integer position, or wheel tick, increments at a rate equal to the angular velocity of each wheel. The number is reset to 0 once it reaches the maximum within its allocated range (wmax). Depending on the vehicle’s equipment, the range of the integer and the rotation angle per unit can be different. The data collection system used in this study recorded the wheel tick value of each wheel at a rate of 100 ms. To estimate the angular velocity of one wheel, the difference between two consecutive wheel tick values wt and wt+1 are used
ω={wt+1wt0.1  sif  wt+1wtwt+1+(wmaxwt)0.1  sotherwise
(1)
where ω = angular velocity in wheel ticks per second. The angular velocities of all four wheels of the vehicle were measured using this method and recorded at each time interval to check for large discrepancies.
Because the wheels on a vehicle travel at different velocities while turning due to the difference of the inner and outer track radii, the discrepancies of the wheels across the front and rear axles were ignored. Only the differences along each side of the vehicle were assessed. Fig. 3(a) shows the recorded angular velocities of each wheel over a 3.4 s time interval during a deceleration event. The front driver (fd), front passenger (fp), rear driver (rd), and rear passenger (rp) angular velocities are plotted on the graph. Between 02:49.9 and 02:50.2, a noticeable divergence in wheel speed is evident for the front passenger wheel. This wheel started to travel slower than the other three wheels, and at Callout i in Fig. 3(a) this was most pronounced. The angular velocity difference (Δω) along each side of the vehicle is computed as follows:
Δωdriverside=|ωfdωrd|
(2)
Δωpassengerside=|ωfpωrp|
(3)
Fig. 3. Vehicle slip detection using wheel ticks: (a) measured angular velocity of each wheel; (b) differential angular velocity on passenger and driver side wheels; and (c) location of deceleration and ABS trigger. (Image by authors.)
Fig. 3(b) plots the differences along each side of the vehicle. The greatest value was recorded at 02:50.3 with an absolute difference of eight wheel ticks between the front passenger and rear passenger wheels [Fig. 3(b), Callout i]. This was associated with a braking event before a stop sign that triggered an ABS intervention by the vehicle [Fig. 3(c), Callout i].

Braking

The amount of braking applied was measured using static brake line pressures sampled in bar units from the vehicle bus at 200-ms intervals. Fig. 4(a) shows a profile of the pressure during a test procedure where the brake was applied with heavy force over a 10-s interval with the vehicle parked. Callout i indicates the moment when the brake was first applied. The force hit a maximum value of 140 bars at Callout ii and was released at Callout iii. Typical braking pressures from the winter driving conditions tested in this study were recorded to be under 25 bars.
Fig. 4. Brake pressures application test: (a) hard press on the brake pedal; and (b) change in brake pressure over time.
In addition to static pressure, the change in pressure (Δbar/s) can also be calculated over two consecutive data samples pt and pt+1 to measure the rate at which the driver applies or releases braking force
Δbar/s=pt+1pt0.2  s
(4)
Fig. 4(b) plots the change in braking pressure in the same braking test. Two rapid changes in brake application were detected during the experiment as indicated by Callout i about 7–10 s into the experiment. The change in brake pressure reached 53 and 87  bars/s in each of the peaks, respectively. The braking force was held constant between 10 and 17 s, which reflects a zero value for the rate of change (Callout ii). At Callout iii, the brake was rapidly released and registered a rate of 301  bars/s. A benefit of using change in brake pressure over static brake pressure is that it is able to characterize the speed at which the brake pedal is applied and released.

Data Collection

Test Circuit Road Conditions

The test vehicle was driven in a predetermined loop of 11.3 km in the area around West Lafayette, IN. The experiment was repeated on four different days with different conditions in 2018:
February 11 (icy),
March 24 (dry to snowy),
July 4 (dry and sunny), and
July 9 (dry and sunny).
On February 11, the test started and ended in icy conditions. On March 24, the test circuit was dry at the beginning of the day and became snow covered and slippery as the day progressed. In the days with winter conditions, a separate vehicle with a MARWIS pavement friction sensor was also included. Two drivers were assigned to a pair of winter and dry condition days; one driver for February 11 and July 9, and the other driver for March 24 and July 4. The drivers were aware and consented to the data collection, and were told to drive normally as they would under each weather condition. The experiments collected measurements from multiple in-vehicle bus signals, and were repeated along the same road links in the test loop per day. This was to essentially emulate fleet data as time progressed. For example, using vehicle count data from traffic signal controller event logs, the runs were approximately 2% of the total volume for March 24.

Map Segmentation

The vehicle provided a constant stream of wheel ticks from which the differential angular velocity from the same side of the vehicle was calculated. Each data point is also associated with a GPS coordinate identifying the vehicle’s position. In order to visualize the segment of roadway affected with the highest slip differentials, the measured data was attributed to a road link. This process would also help hide precise locations for future use cases where privacy is a concern. Moreover, GPS coordinates are not always perfectly aligned with the road. In order to match measurements of repeated passes to a road, links from a HERE Global B.V navigation map were leveraged for this study. Road links consist of geometrical lines defined by several latitude and longitude points. Each road link has a unique identifier with starting and ending points along the roadway. The measurement of the vehicle was mapped to the nearest road link identifier using a perpendicular distance (dashed line in Fig. 5). If multiple slip measurements fall onto the same road link, they are averaged. The 48-road links used for this study are shown in Table 1.
Fig. 5. Assigning road link from vehicle path.
Table 1. Segment details for test loop
Link IDLength (km)Limit (km/h)End of segmentMovementLink IDLength (km)Limit (km/h)End of segmentMovement
10.1980No controlThough250.0432No controlThrough
20.2080No controlRight turn260.0432No controlThrough
30.2064No controlThrough270.2732YieldRoundabout
40.1964No controlThrough280.0532No controlThrough
50.3664Stop signThrough290.1056No controlThrough
60.2164Stop signRight turn300.0456No controlThrough
70.2464No controlThrough310.0556No controlThrough
80.2064No controlThrough320.0656No controlThrough
90.1664No controlThrough330.0956No controlThrough
100.1564No controlThrough340.4064SignalRight turn
110.0864No controlThrough352.5688SignalLeft turn
120.6280No controlThrough360.1264No controlThrough
130.1664No controlThrough370.0180SignalThrough
140.2664No controlThrough380.1256No controlThrough
150.1256No controlThrough390.0280SignalThrough
160.1556No controlThrough401.2180SignalLeft turn
170.2756SignalThrough410.5280SignalThrough
180.0956No controlThrough420.0180No controlThrough
190.0756No controlThrough430.4256SignalRight turn
200.0856SignalRight turn440.0332RoundaboutRoundabout
210.0556No controlThrough450.0132RoundaboutRoundabout
220.1656No controlThrough460.0432No controlRoundabout
230.3080No controlThrough470.3080No controlThrough
240.2480No controlThrough480.0480No controlThrough

Friction Data

Road surface friction data was collected by a chase vehicle equipped with a MARWIS sensor. Using the map segmentation scheme described in Table 1, the pavement friction measurements were averaged into 15-min bins and assigned to the closest road link from the vehicle at that moment in time. In general, the friction estimates provide reasonable overall segment characterization of friction ranges, but do not have the fidelity to identify localized very short segments with low friction.
An table of friction values for the experiment conducted on March 24 is shown in Table 2. A dash ("—") indicates locations and time bins where no data was collected. On that day, at around 10:15 the coefficient of friction (μ) on the test loop begins to drop below 0.8 due to accumulating snowfall for most links. The links that had the lowest detected friction values were 11 and 34. By 14:00, μ recovers to above 0.7 for most links in the loop.
Table 2. Average friction values reported by MARWIS for each segment on test loop, March 24, 2018
Link ID8:459:009:159:309:4510:0010:1510:3010:4511:1511:4512:0012:1512:3012:4513:0013:1513:3013:4514:0014:15
10.820.820.820.820.820.760.730.730.730.720.750.740.750.750.760.760.76
20.820.820.820.820.820.760.750.730.690.720.700.730.710.740.750.740.76
30.820.820.820.800.560.750.650.680.690.730.700.730.770.77
40.820.820.820.800.700.760.670.650.650.740.710.730.760.770.77
50.820.820.810.790.680.750.590.640.670.650.660.660.770.780.78
60.820.820.820.790.710.710.650.630.690.610.620.660.760.750.78
70.820.820.820.800.760.760.740.630.690.750.690.720.760.760.71
80.820.820.820.810.800.780.760.66 0.690.630.740.750.730.760.77
90.820.820.820.820.780.730.720.660.580.690.610.680.760.720.780.79
100.820.820.820.820.820.780.740.760.660.580.740.700.670.770.770.760.77
110.820.820.810.790.620.660.670.640.640.460.480.580.680.750.76
120.820.820.820.820.820.820.800.750.750.750.750.770.760.750.750.750.760.77
130.820.820.820.820.820.820.780.750.740.740.730.760.770.730.750.760.77
140.820.820.820.820.820.820.800.750.740.740.740.740.760.720.740.760.760.77
150.820.820.820.820.820.790.760.790.770.800.810.810.790.740.790.820.80
160.820.820.820.820.820.800.760.780.790.790.800.800.790.780.810.810.81
170.820.820.820.820.810.800.760.760.760.760.770.780.790.810.760.740.770.780.82
180.820.820.760.770.770.770.79
190.780.77
200.820.810.820.820.820.810.760.730.740.740.750.760.770.770.760.790.780.81
210.810.820.820.820.810.770.730.730.780.760.770.770.760.790.780.81
220.820.820.820.820.820.790.780.770.800.790.780.780.770.780.780.77
230.820.820.820.820.820.790.750.750.730.740.770.760.740.730.730.730.74
240.820.820.820.820.820.780.740.760.750.730.760.760.750.760.750.750.77
250.820.820.820.820.820.800.790.760.770.780.770.790.690.740.790.79
260.820.820.820.820.800.760.760.780.770.780.770.790.770.740.790.790.81
270.820.820.820.820.820.800.760.770.750.760.740.770.760.740.770.780.78
280.820.820.820.820.810.790.770.780.780.770.780.780.760.700.780.790.79
290.820.820.820.820.820.800.740.760.710.760.710.720.790.750.740.740.770.770.81
300.820.820.820.820.820.790.770.790.770.760.810.780.740.740.800.810.81
310.820.820.820.820.810.790.780.810.800.810.800.800.780.810.810.82
320.820.820.820.820.810.770.750.770.770.780.800.770.790.790.800.780.81
330.820.820.820.820.810.770.750.810.770.790.800.790.790.790.800.810.82
340.820.820.820.790.660.590.600.640.620.530.620.720.670.750.770.77
350.820.820.820.820.820.820.790.760.710.760.750.740.750.760.730.740.780.780.79
360.820.820.820.820.820.820.780.770.710.750.720.730.780.770.740.670.770.720.77
370.820.820.820.820.820.820.780.760.710.760.720.730.730.750.730.670.760.720.77
380.820.820.820.820.810.780.760.760.770.780.720.710.730.750.710.740.760.780.77
390.820.820.820.820.810.790.740.770.770.780.770.770.790.79
400.820.820.820.820.820.820.800.770.710.770.770.720.760.770.770.760.760.780.790.79
410.820.820.820.820.820.820.800.780.760.770.770.780.740.770.760.770.770.780.78
420.820.820.820.820.820.810.790.760.770.780.770.790.760.770.770.780.75
430.820.820.820.820.820.820.790.760.750.770.750.770.770.770.750.770.770.78
440.820.820.820.820.820.800.760.760.78 0.770.770.770.770.770.780.78
450.820.820.820.820.820.800.760.760.780.780.770.770.770.770.770.780.78
460.820.820.820.820.820.790.780.760.770.750.780.770.780.770.770.780.78
470.820.820.820.820.820.790.740.760.750.730.750.770.750.750.750.760.79
480.820.820.820.820.790.750.760.750.740.770.760.740.730.730.730.74

Wheel Slip

In contrast to MARWIS data, wheel slip can provide very high-fidelity information on location where low friction levels result in vehicle wheels sliding. In this study, wheel slippage is calculated using the wheel tick method, and is defined by the maximum absolute differential angular velocity of the front and rear wheels on either side of the vehicle. Fig. 6(a) shows a time series plot of the max wheel tick differential during the experiment conducted on February 11. The data is also colorized in grayscale using μ values collected from the MARWIS device (darker colors indicating higher friction values). Callout i indicates the wheel slip event from the example in Fig. 3(b) with the maximum differential value of 8. This was also where the ABS intervened on an approach to a stop sign, and corresponded to a μ value of 0.6 to 0.7.
Fig. 6. Implementation of wheel tick on February 11, 2018: (a) differential angular velocity with MARWIS friction scale; and (b) image of localized snow/ice patch. (Image by authors.)
Callouts ii and iii indicate locations where there was a high μ segment friction value reported by the MARWIS unit, but had some slip among the wheels. Inspecting these locations from photos found isolated patches of snow or ice [Fig. 6(b), Callout ii] that the vehicle likely slipped on but were not identified from the friction data. These discrepancies were due to the time-resolution and difference in sensing technology between the MARWIS lasers and wheel rotation data. High wheel rotation differentials were only detected when the vehicle was braking or accelerating hard in slippery areas, such as during the approach to or departing from an intersection. However, under unfavorable road conditions a driver may be less inclined to perform aggressive maneuvers like they would in dry conditions to intentionally prevent wheel slip. It is impractical to assume that as conditions deteriorate, drivers would proportionally increase the vehicle slip, and therefore this is an inconclusive metric for actual road conditions. As shown by the data from February 11, this may be how the driver behaved in this case. A related but more promising method uses wheel vibrations [T. Svantesson, “Detection of short term irregularities in a road surface,” US Patent No. 14,913,684 (2015)] which we may investigate in future research.
The maximum absolute differential angular velocity of one wheel can be converted to a slip percent by dividing it by the average angular velocities of the other wheels not experiencing slippage. The resulting value can be colorized based on a gradient scale and plotted on a map based on the segmentation scheme defined in Table 1. The map shown in Fig. 7 represents the slip percentages during the drive on March 24 at 13:31. Each link is assigned a color based on the last slip percentage value recorded on that link. The direction of vehicle travel is shown by the white arrows.
Fig. 7. Map of test loop with slip ratio colorization. (Map data © 2018 Google.)

Braking

Static brake pressure (bar) and change in brake pressure (Δbar/s) were recorded in each experiment to determine if a driver changed his braking behavior based on actual road conditions. Fig. 8 shows the user road condition rating [Fig. 8(a)], friction [Fig. 8(b)], speed [Fig. 8(c)], static brake pressure [Fig. 8(d)], and change in brake pressure [Fig. 8(e)] from the experiment conducted on March 24 for Link 34 of the test loop. As indicated in Table 2, the link experienced substantial drops in friction over the middle of the day. The data for all 18 passes at this location were combined together, and each sequential data point was assigned a time index. Photos from the dashboard camera corresponding to notable time indices are shown in Fig. 9. The condition at the beginning of the test was dry with good friction as indicated by Fig. 9(a).
Fig. 8. Eighteen laps at Segment 34 during deteriorating conditions.
Fig. 9. Changing conditions at Segment 34: (a) Index 29; (b) Index 121; (c) Index 191; (d) Index 263; (e) Index 467; and (f) Index 64. (Images by authors.)
The callouts on Fig. 8 document the findings from 18 passes at Segment 34. In Fig. 8(a), Callout i indicates when the driver first noticed the road becoming wet [Fig. 9(b)]. However, the wet road did not impact the friction readings (Callout ii) but was reflected in the driver’s slight reduction in speed of about 1.6 to 4.8  km/h [Fig. 8(c), Callout iii]. The brake pressure applied during the fourth lap [Fig. 8(d), Callout iv] was not any less than the previous three laps. An important note is that the vehicle arrived at the intersection for the third lap indicated by Fig. 8(d), Callout v, where the vehicle did not come to a complete stop before making the turn. However, comparing the rate of change in brake pressure of the previous three laps with the fourth lap, there is a slight increase in the rate of application (positive pressure) of the brakes [Fig. 8(e), Callout vii].
By the sixth lap, the driver likely observed the link had started becoming snow-covered [Fig. 8(a), Callout viii]. The measured friction also started to drop [Fig. 8(b), Callout ix]. The top speed on the link also dropped about 13  km/h from top speeds recorded in the first three laps [Fig. 8(c), Callout x]. During this lap, the brake pressure applied remains high at its peak (Callout xi); however, the rate of change of brake pressure peaked at 43  bars/s [Fig. 8(e), Callout xii], which was more than double of the first three laps. Fig. 9(c) shows the road condition during this period.
By the eighth lap, the friction had dropped to about 0.6 and the roadway was substantially covered with snow or slush [Fig. 9(d)]. The brake pressures start dropping below 15 bars at their peak [Fig. 8(d), Callout xiii] and the rate of brake pressure change also decreased in both the positive and negative directions [Fig. 8(e), Callout xiv].
By the 13th lap, the driver noticed the conditions improving slightly (Callout xv), and the friction was sustained above 0.7 briefly [Fig. 8(b), Callout xvi]. The road condition became more wet than snowy [Fig. 9(e)] and the driver started applying a greater amount of pressure [Fig. 9(d), Callout xvii]. The rate of change in brake pressure increased in both the positive and negative directions [Fig. 8(e), Callout xviii].
By the 17th lap, the road condition on the link had been rated consistently wet for a few laps by the driver [Fig. 8(a), Callout xix] and was reflected in the dashboard camera [Fig. 9(f)]. The friction value also approached 0.8 [Fig. 8(b), Callout xx], close to values reported on dry pavement. While the static brake pressures applied remained consistently higher during the high friction laps compared to the low friction laps, the rate of change of brake pressure by qualitative observation was not substantially different between laps of high and low friction [Fig. 8(e), Callout xxi]. The laps with the greatest variation in rate of change in brake pressure appeared to be during the shoulder periods where the conditions start deteriorating substantially or when the conditions start improving; however, it is difficult to identify clear trends from the time series representation of the data. This behavior is explored further using statistical methods in the next section.

Analysis

Braking Behavior by Speed

The static brake pressure and change of brake pressure data collected from icy conditions on February 11 and dry conditions from July 9 were compared using cumulative frequency distribution (CFD) diagrams. Data from the same driver and on the same circuit was collected for each day. Fig. 10 shows the CFDs for static brake pressure [Figs. 10(a–e)] and braking rate [Figs. 10(f–j)] with the data grouped into five speed categories. Braking data below 1.6  km/h were excluded from the analysis to filter out any arbitrary braking application, such as at a traffic signal when the vehicle was stationary.
Fig. 10. Distribution of brake pressures and braking rates in icy and dry conditions at various speeds: (a) 1.614  km/h; (b) 1630  km/h; (c) 3246  km/h; (d) 4862  km/h; (e) >64  km/h; (f) 1.614  km/h; (g) 1630  km/h; (h) 3246  km/h; (i) 4862  km/h; and (j) >64  km/h.
Figs. 10(a–e) shows that for all of the speed categories, the pressure applied was greater during the dry condition than the icy condition. At 1.614  km/h, the pressure was 3 bars greater at the median during the dry condition, and this doubles to 6 bars at 1630  km/h, or about 60%. At 3246  km/h, the difference was 5 bars at the median, but at 4863  km/h and above 64  km/h, this difference decreases to about 2 bars. These figures suggest that at higher speeds and lower speeds, there tend to be smaller differences in braking application than when the vehicle was traveling at speeds between the two ranges. The speed range at which the greatest median difference in brake pressure applied between the dry and icy condition was from 16 to 30  km/h.
Figs. 10(f–j) shows the change of brake pressure during icy conditions was more conservative overall than during dry conditions. In the range of 1.614  km/h, the change of brake pressure applied was about 3  bars/s faster at the 80th percentile than in icy conditions, and the change of brake pressure released was also about 3  bars/s faster at the 20th percentile. In the range of 1630  km/h, the change of brake pressure applied was similar to 1.614  km/h at 3  bars/s at the 80th percentile, but the change of brake pressure released had almost no difference at the 20th percentile. In the range of 3246  km/h, the difference in the change of brake pressure between the dry and icy conditions increased to 4.5  bars/s at the 80th percentile. However, the change of brake pressure release in dry conditions continued to approach the rates at icy conditions at these higher speeds. At speeds from 48 to 62  km/h, the change of brake pressure applied continued to diverge from the dry conditions at the higher pressures; at the 80th percentile the change of brake pressure in icy conditions was only about 1.5  bars/s, or 25% of the change of brake pressure in dry conditions. Braking pressures during the icy conditions did not reach over 1.5  bars/s until well above 95th percentile of braking samples. At speeds over 64  km/h, the braking application rate in the dry conditions were slightly lower than at lower speeds. The data from these results suggest that the change of brake pressure applied at speeds between 32 and 62  km/h can be used as an indicator to differentiate road conditions based on user behavior.
Using a two-sample Kolmogorov-Smirnov (KS) test, there was found to be high significance (P-value<0.001) disproving the null hypothesis between the braking samples of the two dates at all speed categories.

Braking Behavior by Friction Level

To see if there was any correlation between the braking rate and the different friction levels experienced on the roadway, the samples were joined spatially and temporally and plotted against each other in Fig. 11. The graph in Fig. 11(a) shows data from the ice storm on February 11 while Fig. 11(b) shows data from the deteriorating winter storm condition on March 24. Because one driver was collecting data on February 11 and another for March 24 on different routes, there will be differences in braking behavior. Normalizing different drivers’ behavior across different routes, under different vehicle and equipment setups are subjects for many future research topics.
Fig. 11. Change in brake pressure over nondeteriorating and deteriorating winter events: (a) icy conditions with no deterioration; and (b) deterioration from dry to snowy conditions.
On February 11, the friction levels experienced were over a much wider range, from a μ value of 0.82 down to 0.11. The experiment also started in icy conditions with no reduction or improvement of friction conditions through the test period other than on different road segments on different parts of the loop. The range of braking rates were mostly within 50 to 50  bars/s and were pretty consistent throughout the range of friction values. There was a slight increase in the change of brake pressure application and release below μ of 0.2.
On March 24, the experiment started under dry conditions with high friction values (μ>0.8) and within 2 h had deteriorated into snowy conditions with some areas experiencing μ<0.6. In Fig. 11(b), the graphs show the range of friction values for that day were sampled between 0.82μ and 0.38μ. The scatter plot of brake pressure rates follows a tapering conical shape—as the friction values decrease, the extremities of the change of brake pressure also appeared to be curtailed. This may be due to the driver adjusting his braking behavior as the conditions deteriorated throughout the day. Also, for the friction values in the range of 0.65μ0.70μ, there appears to be a noticeable decrease of the vertical bulge in the data, which may hint to the limit where the driver makes their adjustment to changing road conditions over the test period.

Significance of Braking Variance under Different Friction Levels

To verify whether the trends shown in Fig. 11(b) were significant, the Brown-Forsythe test (Brown and Forsythe 1974) for the equality of group variances was used to compare the rates of braking in each friction group. The test determines whether the variances of each group were equal based on each sample’s absolute difference from the group median. The change of brake pressure data from March 24 was divided into four groups using the floor value of the corresponding friction sample to the nearest 0.1μ. Friction values below 0.4μ were not included due to few samples from this day.
Table 3 shows the P-values of the test output for the 10 combinations of the four friction groups. At 95% confidence level, the variance of braking rates was found to be statistically significant between most friction groups except for (0.5μ,0.6μ), and (0.7μ,0.8μ). This suggests that in a scenario where road conditions deteriorated, the rate of change in brake pressure applied is different for different road conditions. However, the data does not suggest whether different drivers will respond accordingly in deteriorating road conditions.
Table 3. P-values from Brown-Forsythe test comparing variance of braking rates at different friction conditions
μ0.50.60.70.8
0.40.0060.036<0.001<0.001
0.50.1770.0010.002
0.6<0.001<0.001
0.70.676

Results and Discussion

The high-frequency in-vehicle data collected for this study under inclement weather conditions suggested that extreme but rare events, such as large differentials of wheel angular velocity, traction control, and ABS intervention, were useful indicators for specific areas of the roadway where there were friction concerns. Because of their frequency of occurrence and pragmatic driver behavior to adjust to such conditions to prevent them, these events were not sufficient to approximate road conditions over a wide area. Moreover, at the local arterials level, because of stochastic variations of controlled stops such as from traffic signals and turning movements, speed may not be the best indicator of deteriorating conditions over a short time period.
The high-frequency braking data examined in this study over two winter storm events suggested that the change of brake pressure may be a useful early indicator of deteriorating road conditions based on the behavioral response of the driver to adjust to the perceived conditions. However, during conditions that were already slippery, the driver may be more careful from the beginning of the drive and not adjust the change of brake pressure too much. Therefore, it would be difficult to assess the conditions from the braking data alone. Fig. 12 shows a time-series plot of the sampled MARWIS friction values and the variance of the change of brake pressure from the March 24 experiment. Only braking data at speeds between 32 and 62  km/h were used, because from the analysis this range was found to have the greatest difference between inclement and noninclement conditions. The data samples were aggregated into 5-min values, taking the minimum friction value sampled and the average variance of the change of brake pressure during each of the 5 min. Callout i indicates the period before any substantial deterioration in conditions were detected; there was a slow ramp-up to around 40  bar2/s2, which then hovered around this rate for 30 min. Callout ii shows the first instance where there was a relatively sharp jump in the variance and a subsequent steep drop within a 5-min period, compared with the previous 60 min. At this time, the minimum friction values obtained from the MARWIS were still excellent at 0.82. This indicates that the driver was perceiving the road conditions changing either due to falling snow or pavement condition and adjusting their braking behavior.
Fig. 12. Braking variance at 3262  km/h and change in friction conditions.
Callout iii shows another jump in the variance to a maximum of 110  bar2/s2, but this time the corresponding minimum friction value also dropped to about 0.75μ, and then down to 0.46μ in the next 5 min. As the minimum friction values fluctuated during this time, the braking variance tapered off and did not show any sharp drops or increases past the initial two spikes at Callout i and ii. Later in the day (Callout iv), as the vehicle traveled between segments of roadway with higher and lower levels of friction, the braking variance showed fewer dramatic spikes compared to the morning, but still had some peaks between 40 and 60  bar2/s2. Later in the afternoon, as the friction levels recovered, the variance also started to peak less, with minimums below 20  bar2/s2. These observations suggest that the change of brake pressure can be used as a surrogate for the observed road condition at the leading edge of a winter storm and during the recovery period as conditions improve. However, these may not be reliable indicators during ongoing events where drivers are already accustomed to the winter conditions. Furthermore, the response behavior may be different for different drivers, so that a longer analysis period would be needed to characterize each driver’s driving style, and identify the thresholds at which nontypical driving conditions start occurring.

Conclusion and Future Work

The experiments conducted in this study demonstrated the application of high-frequency vehicle CAN data for observing winter weather road conditions. A system was deployed to collect wheel tick and brake pressure data at 100 to 200 ms time resolution. A web application was developed to display windshield wipers, hazard lights, traction control, and ABS activations for both real-time and historic events during an ice storm. Additionally, two separate winter events with independent road friction validation were used as test cases to demonstrate the applicability of using CAN data to monitor changing road conditions.
Key findings in this study suggest the following:
1.
A driver may reduce their applied braking pressure in deteriorating road conditions up to 60% at the median intensity;
2.
Braking pressures applied during wintry conditions were most different compared to dry conditions at speeds between 32 and 62  km/h where the heaviest braking was performed;
3.
The variance of the change of brake pressure was found to be significantly different during deteriorating road conditions;
4.
Wheel slip data alone may not account for adjustments to driving behavior that would mask actual slippery road conditions;
5.
Extreme vehicle intervention events such as traction control and ABS were typically rare, even at locations where very low friction values were measured; and
6.
Speed data alone may not be sufficient to characterize changing road conditions on arterials.
Finding 3 is particularly interesting. The driver anticipated deteriorating conditions at least 15 minutes before the road was actually slippery. This has potential for predicting deteriorating conditions ahead of time. Furthermore, it shows that images labeled by humans do not represent the ground truth friction. There is tremendous opportunity in the area of winter road research to leverage the in-vehicle data collection, processing and analysis methods presented in this study over a large fleet and a larger driver pool, and for different length journeys. For instance, different vehicles with front- or rear-wheel drivetrains and type of tires equipped may reveal slip conditions even more effectively. Having data over different equipment and demographic groups would validate or disprove using the data elements currently employed in this study and estimate the sample size needed to correctly identify changing road conditions. For example, some drivers may prefer to brake at the same rate even as conditions deteriorate.
Further tests may also suggest incorporating new data elements such as gas pedal usage, lateral acceleration, and air temperature sensors for determining road conditions. Moreover, applying these and similar methods in a real-time system with an informative interface, given enough protections to encrypt and secure anonymized user data, can be leveraged by agencies to operate and maintain their roads during winter events without costly infrastructure. These technologies can be further synergized with existing networks of automatic vehicle location (AVL) fleets such as snowplows.

Notation

The following symbols are used in this paper:
p
static brake pressure (bars);
w
wheel tick count integer;
μ
coefficient of friction; and
ω
angular velocity (ω/s).

Data Availability Statement

Some or all data, models, or code generated or used during the study are proprietary or confidential in nature and may only be provided with restrictions. Items include the following:
CAN data—wheel tick, brake pressure, speed, and location data can be provided with restrictions;
MARWIS data—no restrictions; and
Dashcam image data—no restrictions.

Acknowledgments

This work was supported by the Indiana Department of Transportation. A test vehicle for data collection, software, and hardware components were provided by the Volkswagen Electronics Research Lab for this study. HERE Global BV provided the API reading the geometry of road links for this research. Google Maps was used for the mapping interface in accordance with guidelines detailed at: https://www.google.com/permissions/geoguidelines/.

Disclaimer

The contents of this paper reflect the views of the authors, who are responsible for the facts and the accuracy of the data presented in this paper, and do not necessarily reflect the official views or policies of the sponsoring organizations. These contents do not constitute a standard, specification, or regulation.

References

Bray, E. T. 2014. “The Javascript object notation (JSON) data interchange format.” Accessed April 7, 2020. https://www.rfc-editor.org/info/rfc7159.
Brown, M. B., and A. B. Forsythe. 1974. “Robust tests for the equality of variances.” J. Am. Stat. Assoc. 69 (346): 364–367. https://doi.org/10.1080/01621459.1974.10482955.
Casselgren, J. 2018. RSI fas II—Road status information. Dynamisk väglagsinformation, AS II—Slutrapport, 154, Borlänge, Sweden: Trafikverket.
Cox, R. W. 1984. “Local area network technology applied to automotive electronic communications.” IEEE Trans. Ind. Electron. Control Instrum. 32 (4): 327–333. https://doi.org/10.1109/TIE.1985.350105.
FHWA (Federal Highway Administration). 2015. Road weather management performance measures survey, analysis, and report. Washington, DC: FHWA.
FHWA (Federal Highway Administration). 2017a. Regional assessment of weather impacts on freight. Washington, DC: FHWA.
FHWA (Federal Highway Administration). 2017b. Road weather management performance measures update. Washington, DC: FHWA.
FHWA (Federal Highway Administration). 2018. “How do weather events impact roads?” Accessed April 7, 2020. ops.fhwa.dot.gov/Weather/q1_roadimpact.htm.
Fielding, R. T. 1996. “Hypertext transfer protocol—HTTP/1.0.” Accessed April 7, 2020. https://www.rfc-editor.org/info/rfc1945.
Gerstenmeier, J., and W. Weber. 1987. “Traction control (ASR): An extension of the anti-lock braking system (ABS).” In Vol. 280 of Proc., 6th Int. Conf. on Automotive Electronics, 29–35. London: Institution of Electrical Engineers.
Hainen, A. M., S. M. Remias, T. M. Brennan, C. M. Day, and D. M. Bullock. 2012. “Probe vehicle data for characterizing road conditions associated with inclement weather to improve road maintenance decisions.” In Proc., IEEE Intelligent Vehicles Symp., 730. New York: IEEE.
Hartkopp, O. 2012. “The CAN networking subsystem of the Linux kernel.” In Proc., 13th Int. CAN Conf. (iCC). Nuremberg, Germany: CAN in Automation.
Hartley, J. 1974. “Anti-lock braking.” In Vol. 13 of Automotive design engineering. London: Design Engineering Publications.
Kiencke, U., S. Dais, and M. Litschel. 1986. “Automotive serial controller area network.” In Proc., SAE Technical Papers, International Congress and Exposition. Warrendale, PA: SAE International.
McMillin, R. J. 1984. “Communications, control, and display.” In Proc., Int. Congress on Transportation Electronics, 22–24. New York: IEEE.
Pisano, P., L. Goodwin, and M. Rossetti. 2008. “Highway crashes in adverse road adverse weather conditions.” In Proc., 88th Annual Meeting American Meteorological Society. Boston: American Meteorological Society.
Zhao, J., J. Zhang, and B. Zhu. 2014. “Development and verification of the tire/road friction estimation algorithm for antilock braking system.” In Mathematical problems in engineering. New York: Hindawi Publishing.

Information & Authors

Information

Published In

Go to Journal of Transportation Engineering, Part A: Systems
Journal of Transportation Engineering, Part A: Systems
Volume 146Issue 8August 2020

History

Received: Nov 30, 2018
Accepted: Jan 6, 2020
Published online: May 29, 2020
Published in print: Aug 1, 2020
Discussion open until: Oct 29, 2020

Authors

Affiliations

Senior Software Engineer, Purdue Univ., 207 S Martin Jischke Dr., West Lafayette, IN 47906 (corresponding author). ORCID: https://orcid.org/0000-0001-5230-1006. Email: [email protected]
Senior Embedded Software Engineer, Volkswagen Group of America, 500 Clipper Dr., Belmont, CA 94002. ORCID: https://orcid.org/0000-0001-5695-3490. Email: [email protected]
Transportation Research Engineer, Purdue Univ., 207 S Martin Jischke Dr., West Lafayette, IN 47906. ORCID: https://orcid.org/0000-0003-0670-6234. Email: [email protected]
Nitin Navali [email protected]
Research Assistant, 238B Transportation Engineering Division, Dept. of Civil Engineering, Indian Institute of Technology Madras, Chennai 600036, Tamil Nadu. Email: [email protected]
Design Engineer, United Consulting, 8440 Allison Pt. Blvd. #200, Indianapolis, IN 46250. ORCID: https://orcid.org/0000-0002-3869-6539. Email: [email protected]
Brandon L. Hardin [email protected]
Project Manager, E&B Paving, Inc., 6211 Gheens Mills Rd., Jeffersonville, IN 47130. Email: [email protected]
Darcy M. Bullock, M.ASCE [email protected]
Professor, Dept. of Civil Engineering, Purdue Univ., 207 S Martin Jischke Dr., West Lafayette, IN 47906. Email: [email protected]

Metrics & Citations

Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

View Options

Media

Figures

Other

Tables

Share

Share

Copy the content Link

Share with email

Email a colleague

Share