<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.2 20190208//EN" "https://jats.nlm.nih.gov/publishing/1.2/JATS-journalpublishing1-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.2" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">3069</journal-id>
<journal-title-group>
<journal-title>Traffic Safety Research</journal-title>
</journal-title-group>
<issn pub-type="epub">n/a</issn>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">143435</article-id>
<article-id pub-id-type="doi">10.55329/uzno7733</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>A data-driven framework for assessing horizontal curve safety using friction and reliability analysis: a case study on an Indian highway</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<contrib-id contrib-id-type="orcid" authenticated="false">https://orcid.org/0009-0005-0286-3568</contrib-id>
<name>
<surname>Kumar</surname>
<given-names>R. Srinivasa</given-names>
</name>
<xref ref-type="aff" rid="author-aff-1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kumar</surname>
<given-names>Vemula Rajesh</given-names>
</name>
<xref ref-type="aff" rid="author-aff-1"><sup>1</sup></xref>
</contrib>
</contrib-group>
<aff id="author-aff-1">
<label>1</label>
<institution-wrap>
<institution content-type="edu">Osmania University, India</institution>
<institution-id institution-id-type="ror">https://ror.org/030sjb889</institution-id>
</institution-wrap>
</aff>
<author-notes>
<corresp id="author-note-1">Corresponding author: <email>rungoz@osmania.ac.in</email></corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="13-06-2026">
<day>13</day>
<month>06</month>
<year>2026</year>
</pub-date>
<volume>10</volume>
<fpage>e000139</fpage>
<lpage>e000139</lpage>
<history>
<date date-type="received" iso-8601-date="13-04-2026">
<day>13</day>
<month>04</month>
<year>2026</year>
</date>
<date date-type="accepted" iso-8601-date="02-06-2026">
<day>02</day>
<month>06</month>
<year>2026</year>
</date>
</history>
<permissions>
<license license-type="open-access">
<ali:license_ref xmlns:ali="http://www.niso.org/schemas/ali/1.0/">
http://creativecommons.org/licenses/by/4.0
</ali:license_ref>
<license-p>
This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0">Creative Commons Attribution License (4.0)</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
</license-p>
</license>
</permissions>
<abstract>
<p>Horizontal curves are frequently associated with elevated crash risk because of the combined effects of vehicle speed, roadway geometry, and tire–pavement friction. This study examines these interactions using field data collected from five horizontal curves on NH-340C in India. An empirical model for estimating the maximum available side friction was developed using measured vehicle speeds, superelevation, and pavement surface texture represented by Mean Texture Depth (MTD). Curve safety was evaluated using three complementary indicators: the Safety Index (SI), Change in Safety Index (ΔSI), and Dynamic Curve Safety Index (DCSI), which together describe local stability conditions and safety variations along the curve. The calibrated friction model shows a decrease in available friction with increasing speed and a moderate increase with pavement texture, producing values in the range of 0.16–0.28 that are consistent with international design guidelines. Model reliability was evaluated using uncertainty and sensitivity analysis techniques, including the delta method, bootstrap resampling, the First-Order Reliability Method, and global sensitivity analysis using Sobol and Morris approaches. The results indicate that pavement texture and operating speed are the dominant factors influencing friction variability, while superelevation has a comparatively smaller effect. In addition, a regression-based model was developed to estimate mid-curve operating speeds, which can support the determination of advisory speeds for curve safety management.</p>
</abstract>
<kwd-group>
<kwd>safety index</kwd>
<kwd>reliability analysis</kwd>
<kwd>global sensitivity</kwd>
<kwd>advisory speed prediction</kwd>
</kwd-group>
<funding-group>
<funding-statement>No agency was involved in funding this study.</funding-statement>
</funding-group>
</article-meta>
</front>
<body>
<sec sec-type="level-A" id="introduction">
  <title>1. Introduction</title>
  <p>Road safety is a global public health priority, as motorization and travel demand continue to rise faster than improvements in infrastructure and enforcement (<xref ref-type="bibr" rid="WHO2023">WHO, 2023</xref>). Horizontal curves are consistently identified as high-risk roadway elements because they require drivers to adjust speed, maintain steering precision, and adapt to varying tire-pavement friction. Failures in these tasks, particularly under low-friction conditions, are a major cause of loss-of-control crashes (<xref ref-type="bibr" rid="AASHTO2018">AASHTO, 2018</xref>; <xref ref-type="bibr" rid="Donnell2016">Donnell et al., 2016</xref>). Although curves represent only a small share of road length, they account for 20–30% of severe crashes, often intensified by excessive speed, inadequate superelevation, surface deterioration, and insufficient friction (<xref ref-type="bibr" rid="Shalkamy2021">Shalkamy et al., 2021</xref>).</p>
  <p>India continues to face a severe road safety challenge. In the year 2022, 461,312 crashes caused 168,491 deaths and 443,366 injuries, equivalent to 53 crashes and 19 deaths every hour, with year-on-year increases of 11.9%, 9.4%, and 15.3% (<xref ref-type="bibr" rid="MoRTH2022">MoRTH, 2022</xref>; <xref ref-type="bibr" rid="MoRTH2023">MoRTH, 2023</xref>; <xref ref-type="bibr" rid="MoRTH2024">MoRTH, 2024</xref>). Curve-related crashes alone accounted for 54,593 incidents and 20,573 deaths, while overturning crashes added 20,070 incidents and 9,829 deaths (<xref ref-type="bibr" rid="MoRTH2022">MoRTH, 2022</xref>). Field records from NH-340C further highlight this risk on five study curves experienced 22 crashes and 12 fatalities between 2020 and May 2025. The Krishnapuram Anjineya Swamy Statue curve recorded nine crashes and three deaths due to rash manoeuvres, while Venkatapuram, despite only three crashes, caused eight fatalities linked to over-speeding and animal crossings. These observations underscore the urgent need for localized safety diagnosis and targeted interventions.</p>
  <p>Key contributing factors in India include heterogeneous traffic ranging from two-wheelers to heavy trucks, rapid pavement deterioration due to monsoons and poor drainage, mismatches between design and operating speeds, and risky behaviours such as late braking and curve-overtaking (<xref ref-type="bibr" rid="Singh2017">Singh, 2017</xref>; <xref ref-type="bibr" rid="Dash2020">Dash et al., 2020</xref>; Shailaja et al., 2022; Rosas-López et al., 2021). Among curve types, simple circular curves pose higher risks than spirals because lateral forces depend heavily on speed, radius, superelevation, and pavement condition (<xref ref-type="bibr" rid="Chen2007">Chen et al., 2007</xref>; Sharf Aldeen et al., 2022).</p>
  <p>Safe negotiation of curves depends on balancing the demand and supply of side friction. The maximum side friction factor (<italic>F<sub>max</sub></italic>) reflects available supply and is shaped by pavement texture, materials, climate, and tire properties (<xref ref-type="bibr" rid="McLean1981">McLean, 1981</xref>; <xref ref-type="bibr" rid="AASHTO2018">AASHTO, 2018</xref>; <xref ref-type="bibr" rid="Gattis2005">Gattis et al., 2005</xref>; <xref ref-type="bibr" rid="Abohassan2022">Abohassan et al., 2022</xref>). Current Indian standards (<xref ref-type="bibr" rid="IRC2023">IRC:73-2023</xref>) assume a fixed <italic>F<sub>max</sub></italic> of 0.15 regardless of speed, texture, or traffic heterogeneity, risking both under- and over-design (<xref ref-type="bibr" rid="Nicholson1998">Nicholson, 1998</xref>; Morrall &amp; Talarico, 1994). Research has consistently shown friction to be speed-dependent. <xref ref-type="bibr" rid="Barnett1938">Barnett (1938)</xref> recommended reductions beyond 60 mph, <xref ref-type="bibr" rid="McLean1981">McLean (1981)</xref> linked failures to inadequate superelevation, and Mudry (1999) demonstrated that 85th percentile friction often exceeded AASHTO limits (<xref ref-type="bibr" rid="Shalkamy2021">Shalkamy et al., 2021</xref>). <xref ref-type="bibr" rid="Bonneson1999">Bonneson (1999</xref>, <xref ref-type="bibr" rid="Bonneson2000">2000)</xref> confirmed higher friction acceptance at lower speeds and suggested using the 95th percentile of speeds for design, while <xref ref-type="bibr" rid="Fitzpatrick2000">Fitzpatrick et al. (2000)</xref> showed that stability is most compromised at midpoints. These findings stand against static assumptions and call for dynamic speed–friction models in safety evaluation.</p>
  <p>Speed profiling has been widely used to study driver behavior, consistently showing deceleration before entry, minimum speed at midpoints, and acceleration at exits. These trends are influenced by curve radius, superelevation, sight distance, shoulder provision, and approach speed (<xref ref-type="bibr" rid="Bennett1994">Bennett, 1994</xref>; <xref ref-type="bibr" rid="Gong2008">Gong &amp; Stamatiadis, 2008</xref>; <xref ref-type="bibr" rid="Bonneson2009">Bonneson &amp; Pratt, 2009</xref>; Shallam &amp; Ahmed, 2016; <xref ref-type="bibr" rid="Llopis-Castello2018">Llopis-Castelló et al., 2018</xref>). Extended deceleration zones often cause discomfort and instability, particularly on circular compared with spiral curves (<xref ref-type="bibr" rid="Chen2007">Chen et al., 2007</xref>; <xref ref-type="bibr" rid="Montella2015">Montella &amp; Imbriani, 2015</xref>; <xref ref-type="bibr" rid="Himes2019">Himes et al., 2019</xref>; <xref ref-type="bibr" rid="Han2023">Han et al., 2023</xref>). Approach speed strongly governs mid- and exit speeds, which in turn influence vehicle stability and crash risk (<xref ref-type="bibr" rid="Himes2019">Himes et al., 2019</xref>). Previous models have attempted to predict operating speeds using geometric variables such as shoulder width, radius, superelevation, and sight distance (<xref ref-type="bibr" rid="Bennett1994">Bennett, 1994</xref>; <xref ref-type="bibr" rid="Gong2008">Gong &amp; Stamatiadis, 2008</xref>; <xref ref-type="bibr" rid="Bonneson2009">Bonneson &amp; Pratt, 2009</xref>; Shallam &amp; Ahmed, 2016; <xref ref-type="bibr" rid="Semeida2014">Semeida, 2014</xref>; <xref ref-type="bibr" rid="Llopis-Castello2018">Llopis-Castelló et al., 2018</xref>). Others stressed the role of approach speed in enabling timely driver alerts, since long deceleration zones reduce comfort and safety (<xref ref-type="bibr" rid="Montella2015">Montella &amp; Imbriani, 2015</xref>).</p>
  <p>With advances in data-driven approaches, machine learning has been increasingly applied in transportation safety. Applications include lane-change prediction (<xref ref-type="bibr" rid="Zhao2022">Zhao et al., 2022</xref>), traffic speed forecasting with SHAP explanations (<xref ref-type="bibr" rid="Lee2023">Lee, 2023</xref>), crash detection from speed analysis (<xref ref-type="bibr" rid="Parsa2020">Parsa et al., 2020</xref>), freeway crash prediction using Bayesian networks and XGBoost (<xref ref-type="bibr" rid="Yang2022">Yang et al., 2022</xref>), real-time cab speed estimation (Shen, 2022), travel time prediction (<xref ref-type="bibr" rid="Chen2021">Chen &amp; Fan, 2021</xref>), and high-speed urban traffic flow modelling (<xref ref-type="bibr" rid="Lu2023">Lu et al., 2023</xref>). However, limited research has focused on predicting reliable operating speeds at approach, midpoint, and exit sections of horizontal curves under mixed-traffic and pavement conditions. This gap highlights the need for models that integrate speed–friction dynamics with local roadway and behavioural factors.</p>
  <p>Advisory speed plaques remain a standard tool for curve safety management (Traffic, 2009). Research has linked advisory speeds to crash–radius relationships (Bonneson et al., 2007; <xref ref-type="bibr" rid="Fitzpatrick2000">Fitzpatrick et al., 2000</xref>; <xref ref-type="bibr" rid="Milstead2011">Milstead et al., 2011</xref>; <xref ref-type="bibr" rid="Montella2015">Montella &amp; Imbriani, 2015</xref>) and factors such as speed differentials (<xref ref-type="bibr" rid="FigueroaMedina2007">Figueroa Medina &amp; Tarko, 2007</xref>; <xref ref-type="bibr" rid="Pratt2015">Pratt et al., 2015</xref>), curve geometry (<xref ref-type="bibr" rid="Dixon2008">Dixon &amp; Rohani, 2008</xref>), side friction (<xref ref-type="bibr" rid="Dixon2008">Dixon &amp; Rohani, 2008</xref>), superelevation (<xref ref-type="bibr" rid="Bonneson2009">Bonneson et al., 2009</xref>; Pratt et al., 2019), and signage (<xref ref-type="bibr" rid="Hummer2010">Hummer et al., 2010</xref>). However, static advisory speeds are often ineffective because of design inconsistencies (<xref ref-type="bibr" rid="Nicholson1998">Nicholson, 1998</xref>), variable driver behaviour (<xref ref-type="bibr" rid="Fleiter2005">Fleiter &amp; Watson, 2005</xref>), attitudes (<xref ref-type="bibr" rid="Hassen2011">Hassen et al., 2011</xref>), and decision-making patterns (Warner &amp; Åberg, 2006). Recent advances favour adaptive systems that continuously monitor vehicles and deliver dynamic feedback. Customized techniques create adjustable alerts (<xref ref-type="bibr" rid="Davis2018">Davis et al., 2018</xref>), while automated systems provide real-time visual or audio warnings when drivers approach curves at unsafe speeds (<xref ref-type="bibr" rid="Davis2018">Davis et al., 2018</xref>; <xref ref-type="bibr" rid="Mahmud2023">Mahmud et al., 2023</xref>).</p>
  <p>Despite international research, little work has focused on developing field-based friction models or adaptive speed frameworks for Indian highways, where mixed traffic, heterogeneous driving behaviour, and pavement surface texture variability exacerbate safety risks. To address this gap, the present study pursues four objectives:</p>
  <list list-type="order">
    <list-item>
      <p>Develop an India-specific maximum side friction function dependent on speed, texture, and superelevation, calibrated against AASHTO envelopes.</p>
    </list-item>
    <list-item>
      <p>Address the need for field-based <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> modelling under Indian traffic and validating by uncertainty analysis.</p>
    </list-item>
    <list-item>
      <p>Evaluate curve safety using three indices: Safety Index (SI), Change in Safety Index (ΔSI), and Dynamic Curve Safety Index (DCSI).</p>
    </list-item>
    <list-item>
      <p>Predict mid-curve operating speeds through regression, with advisory speeds suggested for real-time alerts via adaptive signage.</p>
    </list-item>
  </list>
  <p>By integrating a calibrated <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> model, SI-based diagnostics, and detailed speed profiles, this study offers a data-driven framework for identifying hazardous curves and guiding India-specific design standards, safety audits, and real-time risk management strategies (<xref ref-type="bibr" rid="Montella2015">Montella et al., 2015</xref>; <xref ref-type="bibr" rid="Abohassan2022">Abohassan et al., 2022</xref>).</p>
</sec>
<sec sec-type="level-A" id="methodology">
  <title>2. Methodology</title>
  <p>This study develops a framework for evaluating horizontal curve safety under Indian roadway conditions through four steps: (i) field data collection, (ii) development of a locally calibrated maximum side friction model (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>), (iii) safety index evaluation, and (iv) mid-curve speed prediction for real-time warnings.</p>
  <p>Field surveys recorded curve geometry (radius, superelevation) and operating speeds (tangent, entry, middle, and exit) as well as pavement texture (Mean Texture Depth, MTD, via the sand patch method). These inputs were used to develop a nonlinear empirical model of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>, expressed as <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>a</mml:mi><mml:mo>*</mml:mo><mml:msup><mml:mrow><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow><mml:mi>b</mml:mi></mml:msup><mml:mo>*</mml:mo><mml:msup><mml:mrow><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>−</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msup><mml:mo>*</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:msup></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>, where parameters <italic>a, b, c,</italic> and <italic>d</italic> were calibrated against <xref ref-type="bibr" rid="AASHTO2018">AASHTO (2018)</xref> values using nonlinear regression to capture Indian conditions. The calibrated <italic>F<sub>max</sub></italic> values were compared with the side friction demand:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq1">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(1)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:msup><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mrow><mml:mi>g</mml:mi><mml:mi>.</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mi>e</mml:mi></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>The actual side friction demand values were quantified for safety assessment using three indices: the Safety Index (SI), Change in Safety Index (ΔSI), and Dynamic Curve Safety Index (DCSI). These indices enabled a systematic assessment of how safety varies from curve entry to exit. Subsequently, mid-curve operating speeds identified as critical for stability were predicted using a regression model, based on geometric and field data. This model generated advisory speed values derived from SI, supporting proactive warnings through adaptive roadside signs or in-vehicle systems.</p>
  <p>The analytical framework integrates deterministic modelling, uncertainty quantification, and safety diagnostics. The nonlinear regression establishes the friction model, while the delta method and bootstrap quantify parameter uncertainty. First-Order Reliability Method (FORM) translates these uncertainties into reliability-based safety probabilities, and global sensitivity analysis identifies dominant influencing variables.</p>
</sec>
<sec sec-type="level-A" id="study-sites-and-data-collection">
  <title>3. Study sites and data collection</title>
  <p>Five horizontal curves were selected on National Highway 340C (NH-340C) in Andhra Pradesh for detailed investigation, considering their accident history and geometric complexity. The study curves represent a broad range of radii, superelevation levels, pavement texture conditions, and crash histories commonly observed on rural two-lane highways in India. The selected sites include the Pamulapadu curve, Banumukula curve, Krishnapuram statue curve, Krishnapuram curve, and Venkatapuram curve. For each location, geometric characteristics such as chainage, radius, superelevation, roadway width, and deflection angle were surveyed to serve as input variables for the development of side friction and safety models. The radii of the selected curves ranged from 50 m to 180 m, while superelevation varied between 4.0% and 7.0%, and carriageway widths were consistently about 7.0–7.2 m. Although the curves were designed for a nominal speed of 80 km/h, field observations revealed notable deviations in actual operating speeds.</p>
  <p>All geometric, traffic, pavement texture, and speed data were collected separately for each individual curve location. Independent field surveys were conducted at each site to capture roadway-specific operating characteristics, pavement conditions, and traffic behavior under mixed-traffic conditions.</p>
  <p><xref ref-type="fig" rid="fig1">Figure 1</xref> presents the geographical location of the five selected study curves along NH-340C in Andhra Pradesh, India. The map illustrates the spatial distribution of the investigated sites between Pamulapadu and Venkatapuram.</p>
<fig id="fig1">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 1</label>
<caption>
<title>Geographical location and representative field photographs of the selected study curves along NH-340C in Andhra Pradesh, India</title>
</caption>
<graphic xlink:href="e000139-Figure1.png" />
</fig>
  <p>The field investigations comprised measurements of spot speeds at tangent (pre-start), entry, middle, and exit locations; derivation of acceleration–deceleration profiles; measurement of Mean Texture Depth (MTD) using the sand-patch method; and detailed surveying of curve geometry. The MTD values at curve midpoints ranged from 0.380 mm to 1.123 mm, reflecting pavement surface variability. <xref ref-type="table" rid="tab1">Table 1</xref> summarizes the geometric parameters and surface conditions of the study curves.</p>
<table-wrap id="tab1">
<object-id pub-id-type="publisher-id"></object-id>
<label>Table 1</label>
<caption>
<title>Geometric parameters of the five study curves on NH-340C</title>
</caption>
<table>
<thead>
        <tr>
          <th style="background-color:rgb(204,204,204)" align="left">Curve #</th>
          <th style="background-color:rgb(204,204,204)" align="center"><italic>#1</italic></th>
          <th style="background-color:rgb(204,204,204)" align="center"><italic>#2</italic></th>
          <th style="background-color:rgb(204,204,204)" align="center"><italic>#3</italic></th>
          <th style="background-color:rgb(204,204,204)" align="center"><italic>#4</italic></th>
          <th style="background-color:rgb(204,204,204)" align="center"><italic>#5</italic></th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <td align="left">Location</td>
          <td align="center">Pamulapadu near Milk Dairy</td>
          <td align="center">Banumukula curve</td>
          <td align="center">Krishnapuram Anjineya Swamy curve</td>
          <td align="center">Krishnapuram curve</td>
          <td align="center">Venkatapuram curve</td>
        </tr>
        <tr>
          <td align="left">KM Post</td>
          <td align="right">61+780</td>
          <td align="right">66+180</td>
          <td align="right">78+780</td>
          <td align="right">79+580</td>
          <td align="right">80+980</td>
        </tr>
        <tr>
          <td align="left">Coordinates</td>
          <td align="right"><p>15°50'32.64&quot;N</p>
          <p>78°29'57.78&quot;E</p></td>
          <td align="right"><p>15°49'57.12&quot;N</p>
          <p>78°32'32.56&quot;E</p></td>
          <td align="right"><p>15°52'58.51&quot;N</p>
          <p>78°38'9.76&quot;E</p></td>
          <td align="right"><p>15°52'58.57&quot;N</p>
          <p>78°38'34.60&quot;E</p></td>
          <td align="right"><p>15°53'0.41&quot;N</p>
          <p>78°39'21.29&quot;E</p></td>
        </tr>
        <tr>
          <td align="left">Radius (m)</td>
          <td align="right">120</td>
          <td align="right">50</td>
          <td align="right">180</td>
          <td align="right">100</td>
          <td align="right">110</td>
        </tr>
        <tr>
          <td align="left">Superelevation (%)</td>
          <td align="right">5.0</td>
          <td align="right">7.0</td>
          <td align="right">4.0</td>
          <td align="right">5.0</td>
          <td align="right">5.0</td>
        </tr>
        <tr>
          <td align="left">Carriageway width (m)</td>
          <td align="right">7.20</td>
          <td align="right">7.10</td>
          <td align="right">7.00</td>
          <td align="right">7.10</td>
          <td align="right">7.20</td>
        </tr>
        <tr>
          <td align="left">Design speed (Kmph)</td>
          <td align="right">80</td>
          <td align="right">80</td>
          <td align="right">80</td>
          <td align="right">80</td>
          <td align="right">80</td>
        </tr>
        <tr>
          <td align="left">R1 -Inner radius(m)</td>
          <td align="right">118.20</td>
          <td align="right">48.25</td>
          <td align="right">178.25</td>
          <td align="right">98.25</td>
          <td align="right">108.20</td>
        </tr>
        <tr>
          <td align="left">R2 - outer radius (m)</td>
          <td align="right">121.80</td>
          <td align="right">51.75</td>
          <td align="right">181.75</td>
          <td align="right">101.75</td>
          <td align="right">111.80</td>
        </tr>
        <tr>
          <td align="left">Deflection Angle (Δ<sub>1</sub>)</td>
          <td align="right">12 30'</td>
          <td align="right">50° 58′ 27″</td>
          <td align="right">26° 31′ 00″</td>
          <td align="right">15° 01′ 00″</td>
          <td align="right">14° 47′ 00″</td>
        </tr>
        <tr>
          <td align="left">Deflection Angle (Δ<sub>2</sub>)</td>
          <td align="right">12 30'</td>
          <td align="right">50° 00′ 44″</td>
          <td align="right">26° 31′ 00″</td>
          <td align="right">16° 09′ 00″</td>
          <td align="right">14° 48′ 00″</td>
        </tr>
        <tr>
          <td align="left">MTD<sub>middle</sub> Average (mm)</td>
          <td align="right">0.391</td>
          <td align="right">0.993</td>
          <td align="right">0.570</td>
          <td align="right">1.123</td>
          <td align="right">0.380</td>
        </tr>
      </tbody>
    </table>
  </table-wrap>
  <sec sec-type="level-B" id="traffic-survey-procedure">
    <title>3.1 Traffic survey procedure</title>
    <p>In addition to geometry, traffic flow data were collected to assess exposure risk. Traffic flow and speed surveys were conducted during multiple field campaigns between January 2024 and May 2025 under representative traffic conditions. Classified traffic volume counts were conducted at representative sections along the study corridor using manual observations supported by video recording for verification. The observed traffic volumes were converted into Passenger Car Units (PCU/day) using <xref ref-type="bibr" rid="IRC2023">IRC:73-2023</xref> recommended equivalency factors for heterogeneous traffic streams. The annual average daily traffic was estimated at 3,975 passenger car units per day, with a peak hour factor of 0.74, indicating marked variations between peak and off-peak periods. The traffic stream was dominated by two-wheelers (50.7%), followed by three-wheelers (15.0%) and cars/jeeps (14.3%). Heavy vehicles contributed smaller proportions, including light commercial vehicles (9.4%), buses (5.1%), trucks (4.1%), and non-motorized vehicles (1.4%).</p>
    <p>To investigate vehicle operating behaviour under mixed traffic-flow conditions, spot speed surveys were conducted along the five selected horizontal curves on NH-340C during multiple field campaigns between January 2024 and May 2025. Surveys were performed during daylight hours under dry pavement conditions to minimize the influence of rainfall, poor visibility, and wet-surface friction variations. Observations were collected during representative morning, mid-day, and evening periods under free-flow traffic conditions while avoiding abnormal situations such as severe congestion, maintenance activity, or traffic diversions. High-precision radar speed guns were used to record spot speeds at four locations along each curve at, tangent (pre-start), entry, midpoint, and exit sections. The collected data covered major vehicle categories including two-wheelers, three-wheelers, cars, buses, and trucks. Each survey session continued until stable speed distributions were obtained for all vehicle classes.</p>
    <p>The speed profiles (<xref ref-type="fig" rid="fig2">Figure 2</xref>) indicate a consistent trend of deceleration from Pre-Start to the Middle of the curve, followed by partial recovery toward the Exit. Cars exhibited the highest operating speeds across all sections, with average speeds of 79.1 km/h at the tangent and 61.4 km/h at the curve midpoint, whereas three-wheelers recorded the lowest midpoint speeds (24.6 km/h). Standard deviation analyses revealed greater heterogeneity in car speeds (16.9–18.0 km/h) compared to the more uniform speeds of three-wheelers (7.6–8.4 km/h). Maximum speed observations reinforced these differences, with cars reaching up to 119 km/h at tangents while three-wheelers rarely exceeded 38 km/h at the curve midpoint. These results demonstrate that vehicle type and maneuverability strongly influence curve negotiation behavior in mixed traffic-flow conditions.</p>
<fig id="fig2">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 2</label>
<caption>
<title>The speed profiles of vehicles travelling on circular horizontal curves</title>
</caption>
<graphic xlink:href="e000139-Figure2.png" />
</fig>
    <p>Analysis of longitudinal dynamics highlighted distinct behavioral phases: sharp deceleration during curve entry (up to –6.3 m/s<sup>2</sup> for cars), sustained negative acceleration through the curve midpoint (–0.8 to –1.3 m/s<sup>2</sup> across vehicle types), and mild positive acceleration during exit (+0.5 to +1.2 m/s<sup>2</sup>). Although drivers attempted to regain speed at the exit, recovery seldom matched tangent speeds, emphasizing the restrictive influence of horizontal curvature. The consistent patterns across vehicle types, coupled with variation in acceleration magnitudes, reflect differences in stability, driver response, and traffic heterogeneity.</p>
    <p>The observed speed patterns, characterized by cautious entry speeds, controlled mid-curve travel, and partial exit recovery, reflect adaptive driver behavior under friction constraints. Repeated deceleration and limited speed recovery highlight potential instability when demand friction nears availability, reinforcing the necessity of safety indices to quantitatively relate driver response to curve geometry in mixed traffic environments.</p>
  </sec>
</sec>
<sec sec-type="level-A" id="development-of-fmax-for-indian-conditions">
  <title>4. Development of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> for Indian conditions</title>
  <p>A maximum side friction factor (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>) model for Indian conditions was developed using field-measured parameters, including Mean Texture Depth (MTD), midpoint operating speed, and superelevation (e). These inputs were analyzed using a nonlinear regression framework to ensure that the resulting model realistically captures pavement–tire interactions and mixed-traffic dynamics.</p>
  <p>The general model structure was specified as:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq2">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(2)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>a</mml:mi><mml:mo>*</mml:mo><mml:msup><mml:mrow><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow><mml:mi>b</mml:mi></mml:msup><mml:mo>*</mml:mo><mml:msup><mml:mrow><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>−</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msup><mml:mo>*</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:msup></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>where constants <italic>a, b, c,</italic> and <italic>d</italic> were estimated through optimization.</p>
  <p>For each observation <italic>i</italic>, the model prediction is:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq3">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(3)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>a</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>b</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>c</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>d</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mo>×</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mrow><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow><mml:mi>i</mml:mi><mml:mi>b</mml:mi></mml:msubsup><mml:mo>×</mml:mo><mml:msubsup><mml:mrow><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mi>−</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>×</mml:mo><mml:msubsup><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:msubsup></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Residual at row i (what the optimizer tries to drive toward zero):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq4">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(4)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>a</mml:mi><mml:mo>,</mml:mo><mml:mi>b</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo>,</mml:mo><mml:mi>d</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>A</mml:mi><mml:mi>A</mml:mi><mml:mi>S</mml:mi><mml:mi>H</mml:mi><mml:mi>T</mml:mi><mml:mi>O</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>−</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>a</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>b</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>c</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>d</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Stacking residuals into a vector <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:mi>r</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>&#x3F4;</mml:mi><mml:mo stretchy="false" form="prefix">)</mml:mo><mml:mo>&#x2208;</mml:mo><mml:msup><mml:mi>&#x211D;</mml:mi><mml:mi>n</mml:mi></mml:msup></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>, the objective minimized is the sum of squared errors (SSE):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq5">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(5)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:mi>E</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>&#x3F4;</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:msubsup><mml:msub><mml:mi>r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msup><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>&#x3F4;</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Conceptually, (Virtanen et al., 2020) start from an initial guess ϴ<sup>(0)</sup>. Later, Updates and to reduce SSE using a Gauss-Newton-type step with damping (Levenberg-Marquardt) or a trust-region method when bounds are used (<xref ref-type="bibr" rid="Levenberg1944">Levenberg, 1944</xref>; <xref ref-type="bibr" rid="Marquardt1963">Marquardt, 1963</xref>). Finally, stops when the gradient is small or the maximum number of iterations is reached. When no bounds are imposed, <italic>SciPy</italic> historically uses a Levenberg-Marquardt approach for performing nonlinear least-squares minimization (More et al., 1980). It uses the Jacobian to numerically approximate curvature (via finite differences) and take curvature-aware steps by solving (Moré, 1978; Virtanen et al., 2020).</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq6">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(6)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:msup><mml:mi>J</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mo>&#xB7;</mml:mo><mml:mi>J</mml:mi><mml:mo>+</mml:mo><mml:mi>&#x3BB;</mml:mi><mml:mi>I</mml:mi><mml:mo stretchy="false" form="prefix">)</mml:mo><mml:msub><mml:mi>&#x394;</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mi>J</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mo>&#xB7;</mml:mo><mml:mi>r</mml:mi></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>where <italic>J</italic> is the Jacobian of residuals and λ is the damping; J<sup>T</sup>: Transpose of the Jacobian; I: Identity matrix, Δp: Parameter update vector. Added to the current parameter vector in each iteration:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq7">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(7)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>w</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>u</mml:mi><mml:mi>r</mml:mi><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x394;</mml:mi><mml:mi>p</mml:mi></mml:msub></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Residual vector <italic>= r<sub>i</sub> = y<sub>i</sub> − f(x<sub>i</sub>;p)</italic></p>
  <p>For each residual r<sub>i</sub>, the solver uses sensitivities <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:mi>&#x2202;</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mi>&#x2202;</mml:mi><mml:mi>a</mml:mi><mml:mi>,</mml:mi><mml:mi>&#x2202;</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mi>&#x2202;</mml:mi><mml:mi>b</mml:mi><mml:mi>,</mml:mi><mml:mi>&#x2202;</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mi>&#x2202;</mml:mi><mml:mi>c</mml:mi><mml:mi>,</mml:mi><mml:mi>&#x2202;</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mi>&#x2202;</mml:mi><mml:mi>d</mml:mi><mml:mi>,</mml:mi></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> which for the model follow from and <italic>SciPy</italic> approximates these numerically.</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq8">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(8)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mi>.</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mi>X</mml:mi></mml:mrow><mml:mi>i</mml:mi></mml:msub><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mrow><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow><mml:mi>i</mml:mi><mml:mi>b</mml:mi></mml:msubsup><mml:mi>.</mml:mi><mml:msubsup><mml:mrow><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mi>−</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mi>.</mml:mi><mml:msubsup><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:msubsup></mml:mrow></mml:msub>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Hence,</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd><mml:mrow><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:mfrac><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:mfrac><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mrow><mml:mi mathvariant="normal">ln</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi><mml:msub><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:mfrac><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mrow><mml:mi mathvariant="normal">ln</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>s</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace></mml:mrow><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msub><mml:mover><mml:mi>F</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:mfrac><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mi>l</mml:mi><mml:mi>n</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>differences. Conceptually, it is expected that <italic>b &gt; 0</italic> (greater texture yields more friction), <italic>c &gt; 0</italic> (larger <italic>c</italic> indicates stronger reduction of friction with increasing speed), and <italic>d</italic> may be small but positive (superelevation marginally enhances available friction). Also, takes a &gt; 0.</p>
  <p>To improve numerical stability, a log-linear transformation was applied for initial parameter estimation:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq9">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(9)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>a</mml:mi><mml:mo>+</mml:mo><mml:mi>b</mml:mi><mml:mi>,</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi mathvariant="normal">ln</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi><mml:mo stretchy="false" form="prefix">)</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mi>c</mml:mi><mml:mo>&#xB7;</mml:mo><mml:mi mathvariant="normal">ln</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi><mml:mo stretchy="false" form="prefix">)</mml:mo><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:mo>&#xB7;</mml:mo><mml:mi mathvariant="normal">ln</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>e</mml:mi><mml:mo stretchy="false" form="prefix">)</mml:mo>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Let, <italic>y = ln F<sub>max</sub></italic>, <italic>x<sub>1</sub> = ln MTD</italic>, <italic>x<sub>2</sub> = ln Speed</italic>, <italic>x<sub>3</sub> = ln e</italic></p>
  <p>Runs an ordinary least squares on <italic>y = β0 + β.x<sub>1</sub> + β2.(−x<sub>2</sub>​) + β3.x<sub>3</sub></italic></p>
  <p>Ordinary least squares on this linearized form provided initial guesses and set: <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msup><mml:mi>a</mml:mi><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>O</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:msub><mml:mi>β</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:msup><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:msup><mml:mi>b</mml:mi><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0</mml:mn><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msub><mml:mi>β</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msup><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0</mml:mn><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mi>−</mml:mi><mml:msub><mml:mi>β</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:msup><mml:mi>d</mml:mi><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0</mml:mn><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msub><mml:mi>β</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></p>
  <p>These values were supplied to the nonlinear solver to accelerate convergence. The parameter vector was optimized as:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq10">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(10)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mover><mml:mi>&#x3F4;</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mover><mml:mo>∑</mml:mo><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>𝐜</mml:mi><mml:mi>u</mml:mi><mml:mi>r</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi><mml:mi>_</mml:mi><mml:mi>f</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>f</mml:mi><mml:mo>,</mml:mo><mml:mi>X</mml:mi><mml:mo>,</mml:mo><mml:mi>y</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mover><mml:mi>&#x3F4;</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mover><mml:mi>a</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo>,</mml:mo><mml:mover><mml:mi>b</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo>,</mml:mo><mml:mover><mml:mi>c</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo>,</mml:mo><mml:mover><mml:mi>d</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Where, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mover><mml:mi>&#x3F4;</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> are the best-fit parameters that minimize the sum of squared errors (SSE), <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mover><mml:mo>∑</mml:mo><mml:mo accent="true">̂</mml:mo></mml:mover></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> is the estimated covariance matrix derived from the Jacobian at the solution. From this, the standard error for parameter k was obtained as:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq11">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(11)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>S</mml:mi><mml:mi>E</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:msub><mml:mover><mml:mi>&#x3F4;</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mi>k</mml:mi></mml:msub><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msqrt><mml:msub><mml:mover><mml:mo>∑</mml:mo><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mrow><mml:mi>k</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:msqrt></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>And the 95% confidence intervals were computed as,</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq12">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(12)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mover><mml:mi>&#x3F4;</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mi>k</mml:mi></mml:msub><mml:mspace width="0.222em"></mml:mspace><mml:mo>±</mml:mo><mml:mn>1.96</mml:mn><mml:mspace width="0.222em"></mml:mspace><mml:mi>S</mml:mi><mml:mi>E</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:msub><mml:mover><mml:mi>&#x3F4;</mml:mi><mml:mo accent="true">̂</mml:mo></mml:mover><mml:mi>k</mml:mi></mml:msub><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>The above ensures the statistical reliability of the fitted parameters. The calibrated Indian-specific predictive model was obtained as:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq13">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(13)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mn>1.5721</mml:mn><mml:mo>*</mml:mo><mml:msup><mml:mrow><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mi>−</mml:mi><mml:mn>0.5508</mml:mn></mml:mrow></mml:msup><mml:mo>*</mml:mo><mml:msup><mml:mrow><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>−</mml:mi><mml:mn>0.583</mml:mn></mml:mrow></mml:msup><mml:mo>*</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mn>0.033</mml:mn></mml:msup></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>This functional form captures the nonlinear influence of pavement surface texture, speed, and geometry on the tire-road interaction. Predicted values ranged between 0.16 and 0.28, closely matching the AASHTO reference range of 0.15–0.30 while explicitly capturing Indian roadway conditions.</p>
</sec>
<sec sec-type="level-A" id="uncertainty-analysis">
  <title>5. Uncertainty analysis</title>
  <p>The uncertainty in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> is evaluated by the following methods.</p>
  <sec sec-type="level-B" id="delta-method-variance-estimation">
    <title>5.1 Delta-method variance estimation</title>
    <p>Uncertainty in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>​ due to parameter sensitivity is estimated using the first-order delta method and estimated through the propagation of parameter uncertainty (<xref ref-type="bibr" rid="Beutner2024">Beutner, 2024</xref>). Let <italic>&#x3F4;</italic> = [<italic>a, b, c, d</italic>]<italic><sup>T</sup></italic> and the covariance matrix Σ<italic><sub>θ</sub></italic>. Then:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq14">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(14)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow><mml:mo>≈</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>∇</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:msup><mml:mo>*</mml:mo><mml:mi>T</mml:mi></mml:msup></mml:msubsup><mml:msub><mml:mi mathvariant="normal">Σ</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msub><mml:mi>∇</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>Where <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msub><mml:mi>∇</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> is the gradient vector with respect to parameters:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq15">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(15)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>∇</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:mtable><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:mtable><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:mfrac><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mi>a</mml:mi></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mrow><mml:mi mathvariant="normal">ln</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mrow><mml:mi mathvariant="normal">ln</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>V</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd columnalign="center" style="text-align: center"><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mrow><mml:mi mathvariant="normal">ln</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>e</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mtd></mml:mtr></mml:mtable><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>Assuming independent parameter uncertainties,</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="normal">Σ</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>a</mml:mi><mml:mi>g</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0.05</mml:mn><mml:mi>a</mml:mi><mml:msup><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0.05</mml:mn><mml:mi>b</mml:mi><mml:msup><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mo>,</mml:mo><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0.05</mml:mn><mml:mi>c</mml:mi><mml:msup><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mo>,</mml:mo><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>0.05</mml:mn><mml:mi>d</mml:mi><mml:msup><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>A coefficient of variation of 5% was adopted for parameter uncertainty following common practice in transportation reliability analysis, where geometric and pavement measurement errors typically range between 3–7%. This assumption provides a conservative estimate of model variance while avoiding unrealistic uncertainty amplification (<xref ref-type="bibr" rid="Melchers2018">Melchers and Beck, 2018</xref>)</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq16">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(16)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mtext>The standard deviation of</mml:mtext><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mtext>is:</mml:mtext><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>σ</mml:mi><mml:mi>F</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow></mml:mrow></mml:msqrt></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>and the 95% confidence bounds are given by:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq17">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(17)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>−</mml:mo><mml:mn>1.96</mml:mn><mml:msub><mml:mi>σ</mml:mi><mml:mi>F</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>u</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo>+</mml:mo><mml:mn>1.96</mml:mn><mml:msub><mml:mi>σ</mml:mi><mml:mi>F</mml:mi></mml:msub></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p><xref ref-type="table" rid="tab2">Table 2</xref> presents the numerical outcomes of the delta-method uncertainty analysis. These reveal that <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mspace width="0.222em"></mml:mspace></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> sharply decreases with increasing operating speed, increases moderately with surface texture depth, and is only weakly enhanced by superelevation. The sensitivity magnitudes confirm that pavement texture and speed dominate frictional response, while the geometric term <italic>e</italic> exerts a positive influence.</p>
<table-wrap id="tab2">
<object-id pub-id-type="publisher-id"></object-id>
<label>Table 2</label>
<caption>
<title>Delta-Method Variance and Elasticities Summary</title>
</caption>
<table>
<thead>
          <tr>
            <th style="background-color:rgb(204,204,204)" align="center">Quantity</th>
            <th style="background-color:rgb(204,204,204)" align="center">Symbol</th>
            <th style="background-color:rgb(204,204,204)" align="center">Equation</th>
            <th style="background-color:rgb(204,204,204)" align="center">Typical Value (Mean ± SD)</th>
            <th style="background-color:rgb(204,204,204)" align="center">Interpretation</th>
          </tr>
        </thead>
        <tbody>
          <tr>
            <td align="left">Predicted friction</td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:mn>1.5721</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mrow><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mi>−</mml:mi><mml:mn>0.5508</mml:mn></mml:mrow></mml:msup><mml:mo>×</mml:mo><mml:msup><mml:mrow><mml:mi>S</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>−</mml:mi><mml:mn>0.583</mml:mn></mml:mrow></mml:msup><mml:mo>×</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mn>0.033</mml:mn></mml:msup></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="right">0.18 ± 0.02</td>
            <td align="left">Central model prediction</td>
          </tr>
          <tr>
            <td align="left">Variance (Delta)</td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msub><mml:mi>∇</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:msup><mml:mo>*</mml:mo><mml:mi>T</mml:mi></mml:msup></mml:msubsup><mml:msub><mml:mi mathvariant="normal">Σ</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msub><mml:mi>∇</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="right">1.1 × 10⁻⁴</td>
            <td align="left">Parameter-induced variability</td>
          </tr>
          <tr>
            <td align="left">Std. deviation</td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msub><mml:mi>σ</mml:mi><mml:mi>F</mml:mi></mml:msub></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msqrt><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow></mml:mrow></mml:msqrt></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="right">0.010–0.020</td>
            <td align="left">68 % confidence width</td>
          </tr>
          <tr>
            <td align="left">Elasticity wrt MTD</td>
            <td align="center"><italic>E<sub>MTD</sub> = b</italic></td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="right">−0.5508</td>
            <td align="left">Texture moderately increases friction</td>
          </tr>
          <tr>
            <td align="left">Elasticity wrt Speed</td>
            <td align="center"><italic>E<sub>V</sub> = c</italic></td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>V</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="right">−0.583</td>
            <td align="left">Friction decreases strongly with speed</td>
          </tr>
          <tr>
            <td align="left">Elasticity wrt e</td>
            <td align="center"><italic>E<sub>e</sub> = d</italic></td>
            <td align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>e</mml:mi></mml:mrow></mml:mfrac></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></td>
            <td align="right">+0.033</td>
            <td align="left">Minor geometric contribution</td>
          </tr>
        </tbody>
      </table>
    </table-wrap>
  </sec>
  <sec sec-type="level-B" id="local-elasticity-analysis">
    <title>5.2 Local elasticity analysis</title>
    <p>The elasticity of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>with respect to each input, its percentage sensitivity (De Kroon et al., 2000):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq18">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(18)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>∂</mml:mi><mml:msubsup><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow><mml:mrow><mml:mi>∂</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mi mathvariant="normal">ln</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>Given Eq. (13), the elasticities are considered as:</p>
    <p>E<sub>MTD</sub> = b = −0.5508, E<sub>V</sub> = c = −0.583, E<sub>e</sub> = d = 0.033</p>
    <p>A tornado-style plot (<xref ref-type="fig" rid="fig3">Figure 3</xref>) illustrates the relative elasticities of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>​ with respect to MTD, Speed, and <italic>e</italic>. Speed exhibits the highest negative sensitivity, followed by MTD, confirming its dominant influence on available side friction, whereas <italic>e</italic> contributes comparatively less positive influence.</p>
<fig id="fig3">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 3</label>
<caption>
<title>Elasticity tornado plot of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> model</title>
</caption>
<graphic xlink:href="e000139-Figure3.png" />
</fig>
    <p>Bootstrap confidence intervals: To evaluate the statistical robustness of parameters, nonparametric bootstrap resampling is applied to the log-linearized model in Eq. (9). The model was resampled, <italic>N<sub>boot</sub></italic> = 200 times with replacement, refitted in each iteration, and the empirical distribution of the resulting parameter vectors <italic>θ<sub>i</sub>​</italic> = (<italic>a, b, c, d</italic>) was used to derive mean estimates and 95 % confidence bounds as (Mokhtar et al., 2023):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq19">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(19)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msubsup><mml:mi>θ</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mi>θ</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>h</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mn>2.5</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mn>97.5</mml:mn><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="true" form="prefix">[</mml:mo><mml:msubsup><mml:mi>θ</mml:mi><mml:mi>i</mml:mi><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="true" form="postfix">]</mml:mo></mml:mrow></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>The results summarized in <xref ref-type="table" rid="tab3">Table 3</xref> demonstrate narrow confidence ranges for all four parameters, indicating that the nonlinear calibration is statistically stable and not overly sensitive to random sampling variation.</p>
<table-wrap id="tab3">
<object-id pub-id-type="publisher-id"></object-id>
<label>Table 3</label>
<caption>
<title>Confidence Intervals</title>
</caption>
<table>
<thead>
          <tr>
            <th style="background-color:rgb(204,204,204)" align="left">Parameter (log-linearized)</th>
            <th style="background-color:rgb(204,204,204)" align="center">Mean</th>
            <th style="background-color:rgb(204,204,204)" align="center">Confidence Interval low</th>
            <th style="background-color:rgb(204,204,204)" align="center">Confidence Interval high</th>
          </tr>
        </thead>
        <tbody>
          <tr>
            <td align="left">Scale coefficient</td>
            <td align="right">0.095358</td>
            <td align="right">0.02909</td>
            <td align="right">0.156274</td>
          </tr>
          <tr>
            <td align="left">Texture exponent</td>
            <td align="right">0.001336</td>
            <td align="right">&#x2212;0.00682</td>
            <td align="right">0.010113</td>
          </tr>
          <tr>
            <td align="left">Speed exponent</td>
            <td align="right">0.653567</td>
            <td align="right">0.638691</td>
            <td align="right">0.667525</td>
          </tr>
          <tr>
            <td align="left">Superelevation exponent</td>
            <td align="right">&#x2212;0.00083</td>
            <td align="right">&#x2212;0.02038</td>
            <td align="right">0.019295</td>
          </tr>
        </tbody>
      </table>
    </table-wrap>
    <p><xref ref-type="fig" rid="fig4">Figure 4</xref> interprets conceptually in the model calibration. It demonstrates the inverse power-law dependency between available friction and speed, confirming that as vehicle speed increases, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>​ declines exponentially.</p>
<fig id="fig4">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 4</label>
<caption>
<title><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> vs speed variation</title>
</caption>
<graphic xlink:href="e000139-Figure4.png" />
</fig>
  </sec>
  <sec sec-type="level-B" id="reliability-analysis">
    <title>5.3 Reliability analysis</title>
    <p>Under the First-Order Reliability Method (FORM) framework, the limit-state function separates safe and failure domains as (<xref ref-type="bibr" rid="Zhou2017">Zhou et al., 2017</xref>):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq20">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(20)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>g</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>V</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>R</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>e</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mo>=</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>V</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>R</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>M</mml:mi><mml:mi>T</mml:mi><mml:mi>D</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>e</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo><mml:mo>−</mml:mo><mml:msub><mml:mi>F</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>V</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>R</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>e</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>Failure occurs when g &lt; 0.</p>
    <p>Assuming normal uncertainty in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>, and the corresponding reliability index (β) and failure probability (<italic>P<sub>f</sub>​</italic> ​) are derived from the probabilistic relationship</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq21">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(21)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>β</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mi>g</mml:mi><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mrow><mml:mo stretchy="true" form="prefix">|</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:mo stretchy="true" form="postfix">|</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="normal">Φ</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>−</mml:mi><mml:mi>β</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>where <italic>Φ</italic>(⋅) denotes the cumulative normal distribution function.</p>
    <p>Advisory speeds <italic>V<sub>adv </sub></italic>(<italic>β</italic>) are computed for target reliability thresholds (<italic>β<sub>t</sub></italic> = 0.10,0.05,0.02,0.01) by iteratively solving:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq22">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(22)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>d</mml:mi><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>β</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>The reliability curves plot <italic>P</italic>(<italic>SI</italic>&lt;1) as a function of speed for each vehicle type, and Cars exhibited the highest probability of failure (<xref ref-type="fig" rid="fig4">Figure 4</xref>). The intersections of the probability curve with the dashed horizontal lines (<italic>β<sub>t</sub></italic>​) yield the advisory operating speeds (e.g., 22.56–23.16 m/s for the car case). The <italic>FORM-derived reliability surface</italic>, where smaller <italic>P<sub>f</sub>​</italic> values correspond to safer operating conditions and higher margins of frictional stability.</p>
<fig id="fig5">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 5</label>
<caption>
<title>Reliability curve for the car plot of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>– speed relationship</title>
</caption>
<graphic xlink:href="e000139-Figure5.png" />
</fig>
  </sec>
  <sec sec-type="level-B" id="global-sensitivity-analysis">
    <title>5.4 Global sensitivity analysis </title>
    <p>To assess overall influence ranges, global variance-based methods are applied to quantify the relative influence of texture (MTD), speed, and superelevation (e) on the predicted <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>:</p>
    <p><bold>(i) Sobol’s indices</bold> by Saltelli sampling (<xref ref-type="bibr" rid="Sobol2001">Sobol, 2001</xref>; <xref ref-type="bibr" rid="Saltelli2010">Saltelli et al., 2010</xref>):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq23">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(23)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi></mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:msub><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mi>Y</mml:mi><mml:mo stretchy="true" form="infix">|</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>Y</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:mfrac><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>S</mml:mi><mml:msub><mml:mi>T</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:msub><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi></mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo>∼</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:msub><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mi>Y</mml:mi><mml:mo stretchy="true" form="infix">|</mml:mo><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo>∼</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>V</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mi>Y</mml:mi><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>where <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mspace width="0.222em"></mml:mspace></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msub><mml:mi>S</mml:mi><mml:msub><mml:mi>T</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:msub></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> represent first-order and total effects, respectively.</p>
    <p>Sobol’s indices (<xref ref-type="table" rid="tab4">Table 4</xref>), the global sensitivity indices confirm that <italic>S<sub>T</sub>(MTD)</italic> = 0.84 &gt; <italic>S<sub>T</sub></italic>(<italic>e</italic>) = 0.67 &gt; <italic>S<sub>T</sub></italic>(<italic>Speed</italic>) = 0.51, indicating that pavement texture contributes most strongly to the total variance in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>​<bold>​</bold>, while speed and superelevation have secondary but consistent effects.</p>
<table-wrap id="tab4">
<object-id pub-id-type="publisher-id"></object-id>
<label>Table 4</label>
<caption>
<title>Sobol indices</title>
</caption>
<table>
<thead>
          <tr>
            <th style="background-color:rgb(204,204,204)" align="left">Parameter</th>
            <th style="background-color:rgb(204,204,204)" align="center"><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msub><mml:mi>S</mml:mi><mml:msub><mml:mi>T</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:msub></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></th>
          </tr>
        </thead>
        <tbody>
          <tr>
            <td align="left">MTD</td>
            <td align="center">0.839672</td>
          </tr>
          <tr>
            <td align="left">Speed</td>
            <td align="center">0.511073</td>
          </tr>
          <tr>
            <td align="left">e</td>
            <td align="center">0.672759</td>
          </tr>
        </tbody>
      </table>
    </table-wrap>
    <p>(ii) Morris screening (<xref ref-type="bibr" rid="Morris1991">Morris, 1991</xref>):</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq24">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(24)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msubsup><mml:mi>μ</mml:mi><mml:mi>i</mml:mi><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mrow><mml:mo stretchy="true" form="prefix">|</mml:mo><mml:msub><mml:mrow><mml:mi>E</mml:mi><mml:mi>E</mml:mi></mml:mrow><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="true" form="postfix">|</mml:mo></mml:mrow><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>σ</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>d</mml:mi><mml:mo stretchy="false" form="prefix">(</mml:mo><mml:mrow><mml:mo stretchy="true" form="prefix">|</mml:mo><mml:msub><mml:mrow><mml:mi>E</mml:mi><mml:mi>E</mml:mi></mml:mrow><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="true" form="postfix">|</mml:mo></mml:mrow><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>where <italic>EE<sub>i</sub>​</italic> are elementary effects from randomized trajectories.</p>
    <p>Morris screening (<xref ref-type="table" rid="tab5">Table 5</xref>) further supports this trend: the ranking <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>μ</mml:mi><mml:mi>i</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>(MTD) &gt; <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>μ</mml:mi><mml:mi>i</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>(Speed) ≫ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>μ</mml:mi><mml:mi>i</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>(e) highlights texture as the most influential factor, with superelevation showing minimal impact.</p>
<table-wrap id="tab5">
<object-id pub-id-type="publisher-id"></object-id>
<label>Table 5</label>
<caption>
<title>Morris Indices</title>
</caption>
<table>
<thead>
          <tr>
            <th style="background-color:rgb(204,204,204)" align="left">Parameter</th>
            <th style="background-color:rgb(204,204,204)" align="center">Mean (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>μ</mml:mi><mml:mi>i</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>)</th>
            <th style="background-color:rgb(204,204,204)" align="center">Standard deviation (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msub><mml:mi>σ</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>)</th>
          </tr>
        </thead>
        <tbody>
          <tr>
            <td align="left">MTD</td>
            <td align="right">0.332865</td>
            <td align="right">0.172659</td>
          </tr>
          <tr>
            <td align="left">Speed</td>
            <td align="right">0.294931</td>
            <td align="right">0.151806</td>
          </tr>
          <tr>
            <td align="left">e</td>
            <td align="right">0.020894</td>
            <td align="right">0.010989</td>
          </tr>
        </tbody>
      </table>
    </table-wrap>
    <p>Hence, the models are well-conditioned, with sensitivity concentrated in physically meaningful parameters, ensuring that the predicted friction behavior remains stable and interpretable under realistic operating conditions.</p>
  </sec>
</sec>
<sec sec-type="level-A" id="validation-and-comparative-evaluation-of-fmax-models">
  <title>6. Validation and comparative evaluation of <italic>F<sub>MAX</sub></italic> models </title>
  <sec sec-type="level-B" id="validation-against-aashto-and-international-models">
    <title>6.1 Validation against AASHTO and international models</title>
    <p>The model (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>) was compared with <xref ref-type="bibr" rid="AASHTO2018">AASHTO (2018)</xref> and other international standards, including South Africa (<xref ref-type="bibr" rid="CSIR2000">CSIR, 2000</xref>), <xref ref-type="bibr" rid="Austroads2021">Austroads (2021)</xref>, Sri Lanka (RDA, 1998), Malaysia (REAM, 2002), Abu Dhabi (<xref ref-type="bibr" rid="ADQCC2023">ADQCC, 2023</xref>), Korea (<xref ref-type="bibr" rid="KOTI2013">KOTI, 2013</xref>), <xref ref-type="bibr" rid="Malawi2020">Malawi (2020)</xref>, and <xref ref-type="bibr" rid="IRC2023">IRC:73 (2023)</xref> (<xref ref-type="fig" rid="fig6">Figure 6</xref>). However, AASHTO was used as a benchmark because it provides a well-established speed-dependent friction envelope. In contrast, the IRC standard adopts a constant friction value (0.15) independent of speed or surface texture. Therefore, AASHTO provides a suitable comparative framework for validating the speed-dependent model developed in this study.</p>
    <p>The Indian friction model <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> Eq. (13), aligns closely with the AASHTO friction envelope at moderate operating speeds (30–60 km/h). At lower speeds (&lt;20 km/h), the model predicts slightly higher friction availability due to the influence of pavement surface texture and superelevation. At higher speeds (&gt;80 km/h), the predicted values remain higher than those suggested by AASHTO, reflecting the influence of mixed-traffic behaviour and pavement heterogeneity typical of Indian highways. This difference arises because the proposed model incorporates field-measured pavement texture (MTD) and local operating speeds, which are not explicitly considered in conventional design guidelines.</p>
    <p>Across all models, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> declines consistently with speed (20–140 km/h). For example, AASHTO decreases from 0.36 at 20 km/h to 0.07 at 140 km/h, South Africa from 0.19 at 30 km/h to 0.11 at 120 km/h, and the present study from 0.346 to 0.111. At 60 km/h, the present study predicted 0.183 compared with 0.157 for Malaysia, 0.160 for Sri Lanka, and 0.18 for AASHTO. At 100 km/h, Malaysia and Abu Dhabi estimated 0.115 and 0.12, respectively, while the present study predicted 0.136. Austroads and Korea generally yielded lower values, and Malawi aligned (between 0.11 and 0.14) closely with IRC:73. In contrast, the present study derives friction availability directly from field measurements of pavement texture and operating speeds. Therefore, higher values at lower speeds reflect realistic tire-pavement interaction under Indian roadway conditions. Overall, the present study consistently produced higher friction availability at higher speeds, highlighting the need for region-specific calibration under Indian conditions. This trend reflects the influence of surface texture and road design on safety margins at higher speeds and underscores the need for cautious design at elevated velocities.</p>
<fig id="fig6">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 6</label>
<caption>
<title>Variation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:mspace width="0.222em"></mml:mspace><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> with speed for different models and the present study</title>
</caption>
<graphic xlink:href="e000139-Figure6.png" />
</fig>
    <p>In summary, the Present Study aligns with global models at low speeds but predicts higher side friction availability at high speeds compared to AASHTO, Korea, and Austroads. At 140 km/h, the Present Study (0.111) is well above AASHTO (0.07) and close to the conservative IRC:73 (0.15). While adopting conservative standards, such as AASHTO, ensures robust safety margins, the Present Study provides region-specific insights tailored to Indian pavement conditions and mixed-traffic behavior.</p>
  </sec>
  <sec sec-type="level-B" id="influence-of-curve-radius">
    <title>6.2 Influence of curve radius</title>
    <p>The effect of the curve radius on <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> and speed adaptation was examined using both the Indian model and AASHTO data. In the Indian study, distinct patterns emerged: sharp curves (60–80 m radius) were associated with reduced midpoint speeds (<italic>V<sub>middle</sub></italic>) but higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>; Medium-radius curves (100–140 m) showed wide variability in both speed and friction, and flatter curves (&gt;160 m) supported higher speeds with lower friction demand. The scatter within each group reflects heterogeneous driver adaptation strategies in mixed-traffic conditions.</p>
    <p>AASHTO results showed similar clustering but with more conservative and uniform friction estimates. At small radii (≈60–80 m, 20–50 km/h), the Indian model predicted <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> up to 0.6 compared with AASHTO’s 0.4. For medium radii (≈100–140 m, 40–80 km/h), peaks above 0.5 were observed, while AASHTO values remained ≤0.35. Even at large radii (≈160–180 m, 80–100 km/h), the Indian model predicted 0.2–0.3 compared with AASHTO’s 0.1–0.2. These results confirm that Indian drivers impose higher friction demands than those assumed in international guidelines.</p>
    <p>Overall, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> declines with increasing speed, with the curve radius significantly influencing the speed–friction relationship. Region-specific calibration is therefore vital, as global models tend to underestimate friction demand under Indian conditions, affecting reliable safety evaluation and design interventions.</p>
  </sec>
</sec>
<sec sec-type="level-A" id="the-safety-indices">
  <title>7. The safety indices</title>
  <p>After developing the <italic>F<sub>max</sub></italic> model specific to the measured field conditions and speed profiles of vehicles under mixed traffic-flow conditions, the safety conditions at three different segments of each curve at start, middle, and end locations were evaluated by using the following three safety indices:</p>
  <p>The safety index (SI) is defined as:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq25">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(25)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mfrac><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup><mml:msub><mml:mi>F</mml:mi><mml:mi>d</mml:mi></mml:msub></mml:mfrac></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>The above index represents the immediate safety margin between available and required friction. SI values &gt;1.0 indicate safe operating conditions with surplus friction, SI = 1.0 denotes a critical threshold with no margin of safety, and SI&lt;1.0 indicates unsafe operation with a risk of instability (<xref ref-type="bibr" rid="Donnell2016">Donnell et al., 2016</xref>).</p>
  <p>To capture how safety evolves along the curve, the change in safety index <bold>(</bold>ΔSI<bold>)</bold> was introduced. The variation from:</p>
  <p>Approach to midpoint:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq26">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(26)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>d</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="true" form="prefix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>Midpoint to exit:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq27">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(27)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>d</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>The above two indices distinguish deteriorating conditions (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>S</mml:mi><mml:mi>I</mml:mi><mml:mo>&lt;</mml:mo><mml:mn>0</mml:mn><mml:mo stretchy="false" form="postfix">)</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> from improving conditions (ΔSI&gt;0). Negative values indicate heightened risk zones, typically at entry or middle points, while positive values reflect recovery at exits.</p>
  <p>The dynamic curve safety index (DCSI) is defined as the ratio:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq28">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(28)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:mi>D</mml:mi><mml:mi>C</mml:mi><mml:mi>S</mml:mi><mml:mi>I</mml:mi><mml:mspace width="0.222em"></mml:mspace><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mtext mathvariant="normal">S</mml:mtext><mml:msub><mml:mtext mathvariant="normal">I</mml:mtext><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mtext mathvariant="normal">S</mml:mtext><mml:msub><mml:mtext mathvariant="normal">I</mml:mtext><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:mspace width="0.222em"></mml:mspace><mml:mfrac><mml:msub><mml:mtext mathvariant="normal">a</mml:mtext><mml:mn>1</mml:mn></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mfrac><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>d</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>−</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>−</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>d</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
  <p>The relationship quantifies curve risk balance values near 1.0 denote stability, &gt; 1.0 indicate entry instability, and &lt; 1.0 signal exit instability, while negative DCSI marks critical deceleration–acceleration patterns. Together, SI<bold>,</bold> ΔSI, and DCSI provide an integrated tool for assessing, tracking, and managing curve safety by combining analysis of speed, texture, and geometry.</p>
</sec>
<sec sec-type="level-A" id="curve-safety-evaluation-using-safety-indices">
  <title>8. Curve safety evaluation using safety indices</title>
  <p>Curve safety was assessed by linking operating speeds at approach, midpoint, and exit to side-friction demand and supply, and synthesizing these into three indicators: the Safety Index (SI), the Change in Safety Index (ΔSI), and the Dynamic Curve Safety Index (DCSI). Vehicle-class midpoint speeds formed the basis for estimating friction demand, which was compared with the calibrated <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> model (Eq. (13)) to compute the indices.</p>
  <sec sec-type="level-B" id="safety-index-si">
    <title>8.1 Safety index (SI)</title>
    <p>The SI expresses the balance between available friction (<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math>) and imposed demand (<italic>f<sub>d</sub></italic>) (<xref ref-type="fig" rid="fig7">Figure 7</xref>). Trucks and buses showed the greatest stability, with 76–80% of cases classified as safe (SI &gt; 1.0). Two-wheelers performed even better, with 94.3% of cases safe despite their lower speeds. By contrast, cars showed mixed performance, with 52.3% unsafe, while three-wheelers were the most vulnerable, with 56.9% unsafe at midpoints. These results emphasize vehicle-class disparities: heavy vehicles and two-wheelers were relatively stable due to their lower operating speeds, resulting in lower friction demand, whereas cars often operated near critical friction limits, underscoring the need for targeted speed management. Three-wheelers exhibit higher instability due to their asymmetric three-wheel configuration, narrow wheelbase, and elevated center of gravity, which reduces lateral stability during curve negotiation. In contrast, two-wheelers generally operate at lower speeds and riders tend to adopt conservative maneuvering strategies, resulting in lower friction demand.</p>
<fig id="fig7">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 7</label>
<caption>
<title>Risk condition vehicles in percent observed based on SI Values</title>
</caption>
<graphic xlink:href="e000139-Figure7.png" />
</fig>
  </sec>
  <sec sec-type="level-B" id="si-variation-with-radius-and-speed">
    <title>8.2 SI variation with radius and speed</title>
    <p>SI trends with curve radius reveal that a small radius (≈60–80 m) yielded highly scattered results in the Indian model, while AASHTO predicted predominantly unsafe or near-zero SI. At medium radii (≈100–140 m), Indian outcomes remained widely spread, reflecting heterogeneous driver behavior, whereas AASHTO values were more uniform. At large radii (≈160–180 m), the Indian model stabilized SI in the positive range (≈0.2–0.4) but with higher scatter, while AASHTO produced tightly clustered values (≈0.2–0.3). Thus, the Indian model better reflects mixed-traffic variability, whereas AASHTO offers more conservative, uniform estimates.</p>
  </sec>
  <sec sec-type="level-B" id="change-in-safety-index-δsi">
    <title>8.3 Change in safety index (ΔSI)</title>
    <p>ΔSI highlights safety progression along curves. Cars and two-wheelers largely remained in moderate risk (0 &lt; ΔSI &lt; 1), with 12–15% in high risk (ΔSI &lt; 0). Three-wheelers showed marked deterioration, with 17% unsafe at midpoints (ΔSI<sub>1</sub> &lt; 0) and nearly 48% unsafe at exits (ΔSI<sub>2</sub> &lt; 0). Trucks and buses performed well initially, with 59% safe from approach to midpoint (ΔSI<sub>1</sub> &gt; 1), but 46% deteriorated to unsafe conditions at exits (ΔSI<sub>2</sub> &lt; 0). Overall, safety tended to decline from midpoint to exit, with three-wheelers and heavy vehicles most affected. </p>
  </sec>
  <sec sec-type="level-B" id="dynamic-curve-safety-index-dcsi">
    <title>8.4 Dynamic curve safety index (DCSI)</title>
    <p>DCSI assesses captures entry–exit imbalances (<xref ref-type="fig" rid="fig8">Figure 8</xref>). Most vehicles fell into moderate risk (DCSI &lt; 1.0), but a large proportion faced high risk (DCSI &lt; 0), including all cars (n = 193) and all two-wheelers (n = 35). Severe instability (DCSI &lt;-1 to &lt;- 3) was particularly evident among three-wheelers and heavy vehicles. Compared with SI and ΔSI, DCSI more strongly revealed exit-segment deficits, where recovery from midpoint deceleration was insufficient.</p>
<fig id="fig8">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 8</label>
<caption>
<title>Risk condition vehicles in percent observed based on DCSI values</title>
</caption>
<graphic xlink:href="e000139-Figure8.png" />
</fig>
    <p>Together, SI, ΔSI, and DCSI form a comprehensive diagnostic framework: SI evaluates local stability, ΔSI captures its variation along the curve, and DCSI detects entry-exit imbalances. Integrating operating speed, friction, texture, and geometry, this framework effectively identifies hazardous curves and supports targeted safety interventions such as speed management, alignment correction, and pavement treatment.</p>
  </sec>
</sec>
<sec sec-type="level-A" id="results-and-discussion">
  <title>9. Results and discussion</title>
  <p><italic>Vehicle Speed Profiles:</italic> Spot speed surveys on five horizontal curves of NH-340C revealed marked differences across vehicle classes. Cars averaged 79.1 km/h at tangents and 61.4 km/h at midpoints, while three-wheelers were the slowest (≈24.6 km/h). Speed variability was highest for cars (SD: 16.9–18.0 km/h) and lowest for three-wheelers (7.6–8.4 km/h), reflecting heterogeneous driving behavior. Acceleration analysis showed sharp entry decelerations (up to −6.3 m/s<sup>2</sup> for cars), sustained negative accelerations at midpoints, and partial recovery at exits (+0.5 to +1.2 m/s<sup>2</sup>). These patterns indicate cautious entry, controlled mid-curve maneuvering, and incomplete exit recovery.</p>
  <p><italic>Validation of F<sub>max</sub> and International Comparisons:</italic> The Indian <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> model, derived from midpoint speeds, Mean Texture Depth (MTD), and superelevation, closely matched AASHTO trends at moderate speeds (30–60 km/h) but predicted higher values at low (&lt;20 km/h) and high (&gt;80 km/h) speeds, reflecting local pavement and mixed-traffic conditions (<xref ref-type="fig" rid="fig1">Figure 1</xref>). At low speeds (&lt;20 km/h), friction demand becomes negligible compared with available friction. Consequently, safety indices remain well above the critical threshold (SI&gt;1), resulting in limited differentiation between models. The proposed framework, therefore, provides greater analytical sensitivity at moderate and high operating speeds where friction demand approaches available limits. Comparisons with international models (South Africa, Austroads, Sri Lanka, Malaysia, Abu Dhabi, Korea, Malawi, and IRC) showed that the present study generally yields higher friction availability at higher speeds, highlighting the need for region-specific calibration (<xref ref-type="fig" rid="fig6">Figure 6</xref>).</p>
  <p><italic>Influence of Curve Radius on Speed and Friction:</italic> Curve radius strongly influenced both <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> and operating speeds. The Indian model showed scattered SI and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> values at small radii (≈60–80 m), wide variability at medium radii (≈100–140 m), and partial stabilization at large radii (&gt;160 m) (<xref ref-type="fig" rid="fig6">Figures 6</xref> and <xref ref-type="fig" rid="fig7">7</xref>). In contrast, AASHTO produced more conservative and tightly clustered predictions (<xref ref-type="fig" rid="fig6">Figure 6</xref>). Small-radius curves required <italic>F<sub>max</sub></italic> up to 0.6 under Indian conditions, exceeding the AASHTO limit of 0.4, while medium- and large-radius curves also showed higher and more variable friction demands. These results indicate that Indian drivers impose greater friction demands than international guidelines assume.</p>
  <p><italic>Curve Safety Assessment Using SI, ΔSI, and DCSI:</italic> Safety evaluation highlighted distinct class-specific patterns (<xref ref-type="fig" rid="fig7">Figures 7</xref> and <xref ref-type="fig" rid="fig8">8</xref>). SI indicated that trucks and buses were largely safe (76–80%), two-wheelers were highly stable (94.3%), while cars and three-wheelers often operated near critical limits (52.3% and 56.9% unsafe, respectively). ΔSI showed deteriorating safety toward exits, especially for three-wheelers (ΔSI<sub>1</sub> &lt; 0 for 23 cases, ΔSI<sub>2</sub> &lt; 0 for 113) and heavy vehicles (ΔSI<sub>2</sub> &lt; 0 for 229). DCSI revealed entry–exit imbalances, with all cars and two-wheelers in high-risk zones (DCSI &lt; 0) and severe instability among three-wheelers and heavy vehicles. Collectively, these indices form a comprehensive framework that integrates speed, friction, pavement texture, and geometry to identify hazardous curves and guide targeted interventions.</p>
  <p><italic>Implications for Highway Safety:</italic> The analysis demonstrates three critical points: (i) side friction demand decreases as speed increases, (ii) curve radius strongly moderates the speed–friction relationship, and (iii) international models often underestimate friction requirements under Indian conditions. Adoption of a region-specific model of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> improves the accuracy of safety assessments and supports interventions such as speed regulation, curve realignment, and pavement treatments, thereby enhancing safety in mixed-traffic environments.</p>
  <p><italic>Generalizability of the proposed model</italic>: Although the present study uses five curves from NH-340C, the selected curves cover a broad range of geometric and pavement conditions typical of Indian rural highways. The observed radii range between 50 m and 180 m, superelevation between 4 and 7%, and Mean Texture Depth between 0.38 and 1.12 mm. These ranges represent typical design parameters recommended in IRC standards for two-lane rural highways. Therefore, the calibrated friction model reflects realistic operating conditions under heterogeneous traffic environments. Future research will extend the dataset to multiple highways across different climatic regions of India.</p>
</sec>
<sec sec-type="level-A" id="speed-prediction-model">
  <title>10. Speed prediction model</title>
  <p>Field-observed speed profiles reveal a consistent deceleration trend from the curve approach to its midpoint, followed by partial acceleration toward the exit. The mid-curve represents the most critical control point, where three-wheelers record the lowest operating speeds. At the same time, cars exhibit greater variability and more abrupt acceleration patterns, leading to higher differential speeds among vehicle classes. Such fluctuations highlight the necessity of real-time driver advisory systems to maintain safe, uniform speed behavior on horizontal curves.</p>
  <p>To quantify these dynamics, a regression-based model was developed to predict mid-curve operating speeds (<italic>V<sub>middle</sub></italic>) under field conditions, using measured approach speeds (<italic>V<sub>start</sub></italic>) and geometric parameters. The calibrated model demonstrates a strong linear relationship, expressed as:</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq29">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(29)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>d</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.9821</mml:mn><mml:mspace width="0.222em"></mml:mspace><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mn>4.271</mml:mn></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
<fig id="fig9">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 9</label>
<caption>
<title>Predicted and measured mid-curve speeds (<italic>V<sub>middle</sub></italic>) using the regression model</title>
</caption>
<graphic xlink:href="e000139-Figure9.png" />
</fig>
  <p>Predicted and observed speeds showed close agreement (<xref ref-type="fig" rid="fig9">Figure 9</xref>). From training results, the statistical metrics found as Mean Absolute Percentage Error (MAPE) = 2.72%, Root Mean Square Error (RMSE) = 1.37 km/h, and Coefficient of Determination (R<sup>2</sup>) = 0.975, while testing results yielded MAPE = 2.96%, RMSE = 1.68 km/h, and R<sup>2</sup> = 0.962. These metrics confirm the model’s reliability and minimal error across datasets.</p>
  <p>These results confirm that the model accurately captures the link between approach and midpoint speeds, making it suitable for real-time curve-speed estimation. Trained with field data, it can generate safe speed alerts to improve vehicle stability and reduce excessive friction demand. Such models support safety measures such as dynamic speed display signs (Monsere et al., 2005) and audible warnings (Xu et al., 2023), helping drivers adjust their behavior and reducing crash risk</p>
</sec>
<sec sec-type="level-A" id="proposed-application">
  <title>11. Proposed application</title>
  <p>Advisory speed plaques are widely used for curve safety management (Traffic, 2009). Past studies linked advisory speeds to factors such as curve radius, geometry, side friction, superelevation, and warning signs (Bonneson et al., 2007; <xref ref-type="bibr" rid="Dixon2008">Dixon &amp; Rohani, 2008</xref>; <xref ref-type="bibr" rid="Montella2015">Montella &amp; Imbriani, 2015</xref>; Pratt et al., 2019). However, static advisory speeds often prove inadequate due to inconsistent design, varying driver behavior, and differences in decision-making (<xref ref-type="bibr" rid="Nicholson1998">Nicholson, 1998</xref>; Soole et al., 2013). Recent approaches emphasize adaptive systems that monitor real-time speeds and provide automated warnings through visual or audio feedback (<xref ref-type="bibr" rid="Davis2018">Davis et al., 2018</xref>; <xref ref-type="bibr" rid="Mahmud2023">Mahmud et al., 2023</xref>).</p>
  <p>The safety indices (SI, ΔSI, DCSI) proposed in this study offer a foundation for dynamic curve monitoring. Speed data, obtained via roadside sensors or GPS, can be used to calculate SI and identify unsafe conditions. Risk alerts may then be issued through vehicle-actuated flashing beacons, adaptive electronic signage, or integrated ADAS systems (<xref ref-type="bibr" rid="Antony2021">Antony &amp; Whenish, 2021</xref>). This approach supports immediate feedback, particularly under adverse conditions, and allows adjustments by vehicle type and site-specific characteristics.</p>
  <sec sec-type="level-B" id="practical-application-and-feasibility-of-the-proposed-warning-system">
    <title>11.1 Practical application and feasibility of the proposed warning system</title>
    <p>The conceptual framework is illustrated in <xref ref-type="fig" rid="fig10">Figures 10</xref> and <xref ref-type="fig" rid="fig11">11</xref>. Approach speed data are collected, side friction demand is estimated, and predicted midpoint speeds (Eq. 29) are applied to compute SI. If SI &lt; 1.0, the system triggers a speed reduction warning. Advisory speeds are dynamically calculated using Eq. (30), which incorporates friction demand, superelevation, curve radius, and gravitational effects to define safe operating limits.</p>
<p><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="block" id="eq30">
<mml:mtable displaystyle="true">
<mml:mlabeledtr>
<mml:mtd>
<mml:mtext>(30)</mml:mtext>
</mml:mtd>
<mml:mtd><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>d</mml:mi><mml:mi>v</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mfrac><mml:mrow><mml:mi>g</mml:mi><mml:mi>R</mml:mi><mml:mrow><mml:mo stretchy="true" form="prefix">(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mi>d</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>×</mml:mo><mml:mi>S</mml:mi><mml:mi>I</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">sin</mml:mi><mml:mi>θ</mml:mi><mml:mo stretchy="true" form="postfix">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">cos</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mi>θ</mml:mi></mml:mrow></mml:mfrac></mml:msqrt></mml:mrow>
</mml:mtd>
</mml:mlabeledtr>
</mml:mtable>
</mml:math></p>
    <p>Where the term (<italic>f</italic><sub>dmiddle</sub>×SI) represents <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi mathvariant="normal">max</mml:mi><mml:mo>&#8289;</mml:mo></mml:mrow><mml:mo>*</mml:mo></mml:msubsup></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> reflecting the actual friction demand at the predicted speed and the associated risky margin (SI); <italic>θ</italic> denotes the maximum superelevation, g is the acceleration due to gravity, and R is the curve radius.</p>
    <p>Drivers are alerted in real time through flashing beacons, adaptive signs, or in-vehicle notifications. By continuously updating safe speed thresholds along curve segments, the system ensures proactive interventions, reduces friction-related risks, and enhances safety in mixed-traffic environments.</p>
<fig id="fig10">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 10</label>
<caption>
<title>Control logic framework for applying safety index criteria on circular horizontal curves</title>
</caption>
<graphic xlink:href="e000139-Figure10.png" />
</fig>
<fig id="fig11">
<object-id pub-id-type="publisher-id"></object-id>
<label>Figure 11</label>
<caption>
<title>Conceptual framework for implementation <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mtable displaystyle="false"><mml:mlabeledtr><mml:mtd><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>d</mml:mi><mml:mi>v</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mspace width="0.222em"></mml:mspace></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math> on horizontal curves</title>
</caption>
<graphic xlink:href="e000139-Figure11.png" />
</fig>
  </sec>
</sec>
<sec sec-type="level-A" id="conclusions">
  <title>12. Conclusions</title>
  <p>This study developed a field-based framework for evaluating horizontal curve safety by integrating measured traffic speeds, pavement surface characteristics, and probabilistic reliability analysis. The framework was applied to five horizontal curves located along NH-340C to examine the relationship between operating speeds, pavement texture, and available side friction. The major findings are summarized as follows:</p>
  <list list-type="order">
    <list-item>
      <p>A locally calibrated model for estimating maximum available side friction (Fmax) was developed using field measurements of vehicle speed, Mean Texture Depth (MTD), and superelevation from the study curves. The predicted friction values ranged from 0.16 to 0.28 and showed reasonable agreement with established international guidelines while reflecting the observed pavement and mixed-traffic conditions at the selected sites.</p>
    </list-item>
    <list-item>
      <p>Uncertainty analysis using the Delta method, bootstrap resampling, and First-Order Reliability Method (FORM) indicated that the calibrated friction model remains statistically stable and provides consistent safety margins across the observed operating speed range.</p>
    </list-item>
    <list-item>
      <p>Global sensitivity analysis based on Sobol and Morris methods showed that pavement texture (MTD) has the strongest influence on friction variability, followed by operating speed, while superelevation contributes comparatively less to the variation in available friction.</p>
    </list-item>
    <list-item>
      <p>The safety indices provided a useful means of assessing safety conditions along the studied curves. The results suggest that cars and three-wheelers often operate closer to critical friction limits, whereas trucks, buses, and two-wheelers generally exhibit greater stability due to their lower operating speeds.</p>
    </list-item>
    <list-item>
      <p>A regression-based speed prediction model was developed to estimate mid-curve operating speeds with good accuracy. The model shows potential for supporting advisory speed estimation and safety monitoring on curves with similar geometric and traffic characteristics.</p>
    </list-item>
  </list>
  <p>Overall, the results demonstrate that integrating friction modelling, safety indices, and speed prediction provides a useful analytical approach for evaluating curve safety at the site level and may support future development of adaptive warning systems.</p>
</sec>
<sec sec-type="level-A" id="limitations-and-future-work">
  <title>13. Limitations and future work</title>
  <p>This study presents a field-based framework for assessing horizontal curve safety by relating vehicle operating speeds, pavement texture, and available side friction. However, several limitations should be acknowledged. The analysis was based on observations from five horizontal curves located along NH-340C. Although the selected curves represent a range of geometric and pavement conditions typical of rural highways, the limited number of study sites restricts the extent to which the results can be generalized beyond the investigated corridor. In addition, all field measurements were conducted under dry pavement conditions in order to isolate the influence of geometry and surface texture on friction availability.</p>
  <p>In real-world conditions, environmental factors such as rainfall, pavement wetness, and drainage characteristics can significantly affect tire-pavement interaction and reduce available friction levels. Future studies should therefore incorporate weather-related variables, including rainfall intensity and pavement surface moisture, possibly through continuous monitoring systems or sensor-based measurements. Expanding the dataset to include additional curves from other highway corridors and varying traffic compositions would also improve the robustness and broader applicability of the proposed framework. Further work may also investigate the integration of real-time traffic data and intelligent transportation technologies to support dynamic speed advisory systems for hazardous curves.</p>
 </sec>
<sec sec-type="level-A" id="credit">
  <title>CRediT contribution</title>
  <p><bold>R. Srinivasa Kumar:</bold> Conceptualization, Methodology, Software, Supervision, Writing – original draft, Writing – review &amp; editing. <bold>Vemula Rajesh Kumar:</bold> Data curation, Formal analysis, Investigation.</p>
 </sec>
<sec sec-type="level-A" id="ack">
  <title>Acknowledgements</title>
  <p>The authors are thankful to the Department of Civil Engineering for continuous support and the Traffic Police Department for permission and assistance during field data collection.</p>
 </sec>
<sec sec-type="level-A" id="coi">
  <title>Declaration of competing interests</title>
  <p>The authors have no conflicts of interest with other entities or researchers.</p>
 </sec>
<sec sec-type="level-A">
  <title>Declaration of generative AI use</title>
  <p>The authors declare that no generative AI or AI-assisted tools were used in the preparation of this manuscript.</p>
 </sec>
<sec sec-type="level-A">
  <title>Ethics statement</title>
  <p>The research reported in this paper involved field observations of vehicle speeds and roadway characteristics on public highways. No personal or identifiable data were collected from individuals. Therefore, the study did not involve human subjects and did not require formal ethical approval.</p>
 </sec>
<sec sec-type="level-A" id="fund">
  <title>Funding statement</title>
  <p>No agency was involved in funding this study.</p>
 </sec>
<sec sec-type="level-A">
  <title>Data availability statement</title>
  <p>Data is available upon request.</p>
 </sec>
<sec sec-type="level-A">
  <title>Editorial information</title>
  <p>Handling editor: <bold>Jiří Ambros</bold>, Transport Research Centre (CDV), Czechia.</p>
  <p>Reviewer: <bold>Alyssa Ryan</bold>, Michigan State University, the United States of America.</p>
  <p>Submitted: 13 March 2026; Accepted: 2 June 2026; Published: 13 June 2026.</p>
 </sec>
 </body>
 <back>
 <ref-list>
<ref id="AASHTO2018">
  <element-citation publication-type="book">
    <person-group person-group-type="author">
      <name><surname>AASHTO</surname></name>
    </person-group>
    <article-title>A policy on geometric design of highways and streets</article-title>
    <publisher-name>American Association of State Highway and Transportation Officials</publisher-name>
    <publisher-loc>Washington, DC</publisher-loc>
    <year>2018</year>
  </element-citation>
</ref>
<ref id="Abohassan2022">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Abohassan</surname><given-names>A.</given-names></name>
      <name><surname>El-Basyouny</surname><given-names>K.</given-names></name>
      <name><surname>Kwon</surname><given-names>T. J.</given-names></name>
    </person-group>
    <article-title>Effects of inclement weather events on road surface conditions and traffic safety: An event-based empirical analysis framework</article-title>
    <source>Transportation Research Record</source>
    <year>2022</year>
    <volume>2676</volume>
    <fpage>51</fpage>
    <lpage>62</lpage>
    <object-id pub-id-type="doi">10.1177/03611981221088588</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1177/03611981221088588">https://doi.org/10.1177/03611981221088588</ext-link>
  </element-citation>
</ref>
<ref id="ADQCC2023">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>ADQCC</surname></name>
    </person-group>
    <article-title>TR-514 road design standards</article-title>
    <publisher-name>Abu Dhabi Quality and Conformity Council</publisher-name>
    <year>2023</year>
  </element-citation>
</ref>
<ref id="Antony2021">
  <element-citation publication-type="book">
    <person-group person-group-type="author">
      <name><surname>Antony</surname><given-names>M. M.</given-names></name>
      <name><surname>Whenish</surname><given-names>R.</given-names></name>
    </person-group>
    <article-title>Advanced driver assistance systems (ADAS)</article-title>
    <source>Automotive embedded systems</source>
    <year>2021</year>
    <fpage>165</fpage>
    <lpage>181</lpage>
    <publisher-name>Springer</publisher-name>
    <person-group person-group-type="editor">
      <name><surname>Kathiresh</surname><given-names>M.</given-names></name>
      <name><surname>Neelaveni</surname><given-names>R.</given-names></name>
    </person-group>
  </element-citation>
</ref>
<ref id="Austroads2021">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Austroads</surname></name>
    </person-group>
    <article-title>Guide to road design part 3: Geometric design (Edition 3.4)</article-title>
    <year>2021</year>
  </element-citation>
</ref>
<ref id="Barnett1938">
  <element-citation publication-type="book">
    <person-group person-group-type="author">
      <name><surname>Barnett</surname><given-names>J.</given-names></name>
    </person-group>
    <article-title>Transition curves for highways</article-title>
    <publisher-name>U.S. Government Printing Office</publisher-name>
    <publisher-loc>Washington, DC</publisher-loc>
    <year>1938</year>
  </element-citation>
</ref>
<ref id="Bennett1994">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Bennett</surname><given-names>C. R.</given-names></name>
    </person-group>
    <article-title>A speed prediction model for rural two-lane highways</article-title>
    <year>1994</year>
  </element-citation>
</ref>
<ref id="Beutner2024">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Beutner</surname><given-names>E.</given-names></name>
    </person-group>
    <article-title>Delta method: Asymptotic distribution</article-title>
    <source>Wiley Interdisciplinary Reviews: Computational Statistics</source>
    <year>2024</year>
    <volume>16</volume>
    <issue>1</issue>
    <fpage>e1634</fpage>
    <object-id pub-id-type="doi">10.1002/wics.1634</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1002/wics.1634">https://doi.org/10.1002/wics.1634</ext-link>
  </element-citation>
</ref>
<ref id="Bonneson1999">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Bonneson</surname><given-names>J. A.</given-names></name>
    </person-group>
    <article-title>Side friction and speed as controls for horizontal curve design</article-title>
    <source>Journal of Transportation Engineering</source>
    <year>1999</year>
    <volume>125</volume>
    <issue>6</issue>
    <fpage>473</fpage>
    <lpage>480</lpage>
    <object-id pub-id-type="doi">10.1061/(ASCE)0733-947X(1999)125:6(473)</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1061/(ASCE)0733-947X(1999)125:6(473)">https://doi.org/10.1061/(ASCE)0733-947X(1999)125:6(473)</ext-link>
  </element-citation>
</ref>
<ref id="Bonneson2000">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Bonneson</surname><given-names>J. A.</given-names></name>
    </person-group>
    <article-title>Kinematic approach to horizontal curve transition design</article-title>
    <source>Transportation Research Record</source>
    <year>2000</year>
    <volume>1737</volume>
    <fpage>1</fpage>
    <lpage>8</lpage>
    <object-id pub-id-type="doi">10.3141/1737-01</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3141/1737-01">https://doi.org/10.3141/1737-01</ext-link>
  </element-citation>
</ref>
<ref id="Bonneson2009">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Bonneson</surname><given-names>J. A.</given-names></name>
      <name><surname>Pratt</surname><given-names>M. P.</given-names></name>
    </person-group>
    <article-title>Model for predicting speed along horizontal curves on two-lane highways</article-title>
    <source>Transportation Research Record</source>
    <year>2009</year>
    <volume>2092</volume>
    <fpage>19</fpage>
    <lpage>27</lpage>
    <object-id pub-id-type="doi">10.3141/2092-03</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3141/2092-03">https://doi.org/10.3141/2092-03</ext-link>
  </element-citation>
</ref>
<ref id="Chen2007">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Chen</surname><given-names>S. S.</given-names></name>
      <name><surname>Rakotonirainy</surname><given-names>A.</given-names></name>
      <name><surname>Loke</surname><given-names>S.</given-names></name>
      <name><surname>Krishnaswamy</surname><given-names>S.</given-names></name>
    </person-group>
    <article-title>A crash risk assessment model for road curves</article-title>
    <source>Proceedings of the 20th International Technical Conference on Enhanced Safety of Vehicles</source>
    <year>2007</year>
  </element-citation>
</ref>
<ref id="Chen2021">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Chen</surname><given-names>Z.</given-names></name>
      <name><surname>Fan</surname><given-names>W.</given-names></name>
    </person-group>
    <article-title>A freeway travel time prediction method based on an XGBoost model</article-title>
    <source>Sustainability</source>
    <year>2021</year>
    <volume>13</volume>
    <issue>15</issue>
    <fpage>8577</fpage>
  </element-citation>
</ref>
<ref id="CSIR2000">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>CSIR</surname></name>
    </person-group>
    <article-title>Guidelines for human settlement planning and design: Roads—Geometric design and layout planning</article-title>
    <year>2000</year>
  </element-citation>
</ref>
<ref id="Dash2020">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Dash</surname><given-names>D. P.</given-names></name>
      <name><surname>Sethi</surname><given-names>N.</given-names></name>
      <name><surname>Dash</surname><given-names>A. K.</given-names></name>
    </person-group>
    <article-title>Identifying the causes of road traffic accidents in India: An empirical investigation</article-title>
    <source>Journal of Public Affairs</source>
    <year>2020</year>
    <volume>20</volume>
    <fpage>e2038</fpage>
    <object-id pub-id-type="doi">10.1002/pa.2038</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1002/pa.2038">https://doi.org/10.1002/pa.2038</ext-link>
  </element-citation>
</ref>
<ref id="Davis2018">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Davis</surname><given-names>B.</given-names></name>
      <name><surname>Morris</surname><given-names>N. L.</given-names></name>
      <name><surname>Achtemeier</surname><given-names>J. D.</given-names></name>
      <name><surname>Patzer</surname><given-names>B.</given-names></name>
    </person-group>
    <article-title>In-vehicle dynamic curve-speed warnings at high-risk rural curves</article-title>
    <year>2018</year>
  </element-citation>
</ref>
<ref id="Dixon2008">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Dixon</surname><given-names>K.</given-names></name>
      <name><surname>Rohani</surname><given-names>J.</given-names></name>
    </person-group>
    <article-title>Methodologies for estimating advisory curve speeds on Oregon highways</article-title>
    <publisher-name>Oregon Department of Transportation</publisher-name>
    <year>2008</year>
  </element-citation>
</ref>
<ref id="Donnell2016">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Donnell</surname><given-names>E.</given-names></name>
      <name><surname>Wood</surname><given-names>J.</given-names></name>
      <name><surname>Himes</surname><given-names>S.</given-names></name>
      <name><surname>Torbic</surname><given-names>D.</given-names></name>
    </person-group>
    <article-title>Use of side friction in horizontal curve design: A margin of safety assessment</article-title>
    <source>Transportation Research Record</source>
    <year>2016</year>
    <volume>2588</volume>
    <fpage>61</fpage>
    <lpage>70</lpage>
    <object-id pub-id-type="doi">10.3141/2588-07</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3141/2588-07">https://doi.org/10.3141/2588-07</ext-link>
  </element-citation>
</ref>
<ref id="FigueroaMedina2007">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Figueroa Medina</surname><given-names>A. M.</given-names></name>
      <name><surname>Tarko</surname><given-names>A. P.</given-names></name>
    </person-group>
    <article-title>Speed changes in the vicinity of horizontal curves on two-lane rural roads</article-title>
    <source>Journal of Transportation Engineering</source>
    <year>2007</year>
    <volume>133</volume>
    <issue>4</issue>
    <fpage>215</fpage>
    <lpage>222</lpage>
    <object-id pub-id-type="doi">10.1061/(ASCE)0733-947X(2007)133:4(215)</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1061/(ASCE)0733-947X(2007)133:4(215)">https://doi.org/10.1061/(ASCE)0733-947X(2007)133:4(215)</ext-link>
  </element-citation>
</ref>
<ref id="Fitzpatrick2000">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Fitzpatrick</surname><given-names>K.</given-names></name>
      <etal/>
    </person-group>
    <article-title>Design factors that affect driver speed on suburban arterials</article-title>
    <publisher-name>Texas Transportation Institute</publisher-name>
    <year>2000</year>
  </element-citation>
</ref>
<ref id="Fleiter2005">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Fleiter</surname><given-names>J.</given-names></name>
      <name><surname>Watson</surname><given-names>B.</given-names></name>
    </person-group>
    <article-title>The speed paradox: The misalignment between driver attitudes and speeding behaviour</article-title>
    <source>Australasian Road Safety Research, Policing and Education Conference</source>
    <year>2005</year>
  </element-citation>
</ref>
<ref id="Gattis2005">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Gattis</surname><given-names>J. L.</given-names></name>
      <name><surname>Vinson</surname><given-names>B. F.</given-names></name>
      <name><surname>Duncan</surname><given-names>L. K.</given-names></name>
    </person-group>
    <article-title>Low-speed horizontal curve friction factors</article-title>
    <source>Journal of Transportation Engineering</source>
    <year>2005</year>
    <volume>131</volume>
    <issue>2</issue>
    <fpage>112</fpage>
    <lpage>119</lpage>
    <object-id pub-id-type="doi">10.1061/(ASCE)0733-947X(2005)131:2(112)</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1061/(ASCE)0733-947X(2005)131:2(112)">https://doi.org/10.1061/(ASCE)0733-947X(2005)131:2(112)</ext-link>
  </element-citation>
</ref>
<ref id="Gong2008">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Gong</surname><given-names>H.</given-names></name>
      <name><surname>Stamatiadis</surname><given-names>N.</given-names></name>
    </person-group>
    <article-title>Operating speed prediction models for horizontal curves on rural four-lane highways</article-title>
    <source>Transportation Research Record</source>
    <year>2008</year>
    <volume>2075</volume>
    <fpage>1</fpage>
    <lpage>7</lpage>
    <object-id pub-id-type="doi">10.3141/2075-01</object-id>
    <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3141/2075-01">https://doi.org/10.3141/2075-01</ext-link>
  </element-citation>
</ref>
<ref id="Han2023">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Han</surname><given-names>L.</given-names></name>
      <name><surname>Du</surname><given-names>Z.</given-names></name>
      <name><surname>Zheng</surname><given-names>H.</given-names></name>
      <name><surname>Xu</surname><given-names>F.</given-names></name>
      <name><surname>Mei</surname><given-names>J.</given-names></name>
    </person-group>
    <article-title>Reviews and prospects of human factors research on curve driving</article-title>
    <source>Journal of Traffic and Transportation Engineering</source>
    <year>2023</year>
  </element-citation>
</ref>
<ref id="Hassen2011">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Hassen</surname><given-names>A.</given-names></name>
      <name><surname>Godesso</surname><given-names>A.</given-names></name>
      <name><surname>Abebe</surname><given-names>L.</given-names></name>
      <name><surname>Girma</surname><given-names>E.</given-names></name>
    </person-group>
    <article-title>Risky driving behaviours for road traffic accidents among drivers in Mekele city, Northern Ethiopia</article-title>
    <source>BMC Research Notes</source>
    <year>2011</year>
    <volume>4</volume>
    <fpage>535</fpage>
  </element-citation>
</ref>
<ref id="Himes2019">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Himes</surname><given-names>S.</given-names></name>
      <name><surname>Porter</surname><given-names>R. J.</given-names></name>
      <name><surname>Hamilton</surname><given-names>I.</given-names></name>
      <name><surname>Donnell</surname><given-names>E.</given-names></name>
    </person-group>
    <article-title>Safety evaluation of geometric design criteria: Horizontal curve radius and side friction demand on rural two-lane highways</article-title>
    <source>Transportation Research Record</source>
    <year>2019</year>
    <volume>2673</volume>
    <fpage>516</fpage>
    <lpage>525</lpage>
  </element-citation>
</ref>
<ref id="Hummer2010">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Hummer</surname><given-names>J. E.</given-names></name>
      <name><surname>Rasdorf</surname><given-names>W. J.</given-names></name>
      <name><surname>Findley</surname><given-names>D.</given-names></name>
      <name><surname>Zegeer</surname><given-names>C. V.</given-names></name>
    </person-group>
    <article-title>Procedure for curve warning signing, delineation, and advisory speeds for horizontal curves</article-title>
    <publisher-name>North Carolina Department of Transportation</publisher-name>
    <year>2010</year>
  </element-citation>
</ref>
<ref id="IRC2023">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>IRC</surname></name>
    </person-group>
    <article-title>IRC:73-2023 geometric design standards for non-urban roads</article-title>
    <publisher-name>Indian Roads Congress</publisher-name>
    <year>2023</year>
  </element-citation>
</ref>
<ref id="KOTI2013">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>KOTI</surname></name>
    </person-group>
    <article-title>Ten lessons from road transport policy in Korea</article-title>
    <publisher-name>Korea Transport Institute</publisher-name>
    <year>2013</year>
  </element-citation>
</ref>
<ref id="Lee2023">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Lee</surname><given-names>E. H.</given-names></name>
    </person-group>
    <article-title>Traffic speed prediction of urban road network based on high-importance links using XGBoost and SHAP</article-title>
    <source>IEEE Access</source>
    <year>2023</year>
  </element-citation>
</ref>
<ref id="Levenberg1944">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Levenberg</surname><given-names>K.</given-names></name>
    </person-group>
    <article-title>A method for the solution of certain nonlinear problems in least squares</article-title>
    <source>Quarterly of Applied Mathematics</source>
    <year>1944</year>
    <volume>2</volume>
    <issue>2</issue>
    <fpage>164</fpage>
    <lpage>168</lpage>
  </element-citation>
</ref>
<ref id="Llopis-Castello2018">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Llopis-Castelló</surname><given-names>D.</given-names></name>
      <name><surname>González-Hernández</surname><given-names>B.</given-names></name>
      <name><surname>Pérez-Zuriaga</surname><given-names>A. M.</given-names></name>
      <name><surname>García</surname><given-names>A.</given-names></name>
    </person-group>
    <article-title>Speed prediction models for trucks on horizontal curves of two-lane rural roads</article-title>
    <source>Transportation Research Record</source>
    <year>2018</year>
    <volume>2672</volume>
    <fpage>72</fpage>
    <lpage>82</lpage>
  </element-citation>
</ref>
<ref id="Lu2023">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Lu</surname><given-names>X.</given-names></name>
      <name><surname>Chen</surname><given-names>C.</given-names></name>
      <name><surname>Gao</surname><given-names>R.</given-names></name>
      <name><surname>Xing</surname><given-names>Z.</given-names></name>
    </person-group>
    <article-title>Prediction of high-speed traffic flow around cities based on a BO-XGBoost model</article-title>
    <source>Symmetry</source>
    <year>2023</year>
    <volume>15</volume>
    <issue>7</issue>
    <fpage>1453</fpage>
  </element-citation>
</ref>
<ref id="Mahmud2023">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Mahmud</surname><given-names>M. S.</given-names></name>
      <name><surname>Bamney</surname><given-names>A.</given-names></name>
      <name><surname>Megat Johari</surname><given-names>M. U.</given-names></name>
      <name><surname>Jashami</surname><given-names>H.</given-names></name>
      <name><surname>Gates</surname><given-names>T. J.</given-names></name>
      <name><surname>Savolainen</surname><given-names>P. T.</given-names></name>
    </person-group>
    <article-title>Evaluating driver response to a dynamic speed feedback sign at rural highway curves</article-title>
    <source>Transportation Research Record</source>
    <year>2023</year>
    <volume>2677</volume>
    <fpage>1103</fpage>
    <lpage>1114</lpage>
  </element-citation>
</ref>
<ref id="Malawi2020">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Malawi Roads Authority</surname></name>
    </person-group>
    <article-title>Low volume roads manual: Geometric design and road safety</article-title>
    <year>2020</year>
  </element-citation>
</ref>
<ref id="Marquardt1963">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Marquardt</surname><given-names>D. W.</given-names></name>
    </person-group>
    <article-title>An algorithm for least-squares estimation of nonlinear parameters</article-title>
    <source>SIAM Journal on Applied Mathematics</source>
    <year>1963</year>
    <volume>11</volume>
    <issue>2</issue>
    <fpage>431</fpage>
    <lpage>441</lpage>
  </element-citation>
</ref>
<ref id="Melchers2018">
  <element-citation publication-type="book">
    <person-group person-group-type="author">
      <name><surname>Melchers</surname><given-names>R. E.</given-names></name>
      <name><surname>Beck</surname><given-names>A. T.</given-names></name>
    </person-group>
    <article-title>Structural reliability analysis and prediction</article-title>
    <publisher-name>Wiley</publisher-name>
    <year>2018</year>
  </element-citation>
</ref>
<ref id="McLean1981">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>McLean</surname><given-names>J.</given-names></name>
    </person-group>
    <article-title>Driver speed behaviour and rural road alignment design</article-title>
    <source>Traffic Engineering and Control</source>
    <year>1981</year>
    <volume>22</volume>
  </element-citation>
</ref>
<ref id="Milstead2011">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>Milstead</surname><given-names>R.</given-names></name>
      <name><surname>Qin</surname><given-names>X.</given-names></name>
      <name><surname>Katz</surname><given-names>B.</given-names></name>
      <name><surname>Bonneson</surname><given-names>J. A.</given-names></name>
      <name><surname>Pratt</surname><given-names>M.</given-names></name>
      <name><surname>Miles</surname><given-names>J.</given-names></name>
      <name><surname>Carlson</surname><given-names>P. J.</given-names></name>
    </person-group>
    <article-title>Procedures for setting advisory speeds on curves</article-title>
    <publisher-name>Federal Highway Administration</publisher-name>
    <year>2011</year>
  </element-citation>
</ref>
<ref id="Montella2015">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Montella</surname><given-names>A.</given-names></name>
      <name><surname>Imbriani</surname><given-names>L. L.</given-names></name>
    </person-group>
    <article-title>Safety performance functions incorporating design consistency variables</article-title>
    <source>Accident Analysis &amp; Prevention</source>
    <year>2015</year>
    <volume>74</volume>
    <fpage>133</fpage>
    <lpage>144</lpage>
  </element-citation>
</ref>
<ref id="MoRTH2022">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>MoRTH</surname></name>
    </person-group>
    <article-title>Road accidents in India 2022</article-title>
    <publisher-name>Ministry of Road Transport and Highways, Government of India</publisher-name>
    <year>2022</year>
  </element-citation>
</ref>
<ref id="MoRTH2023">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>MoRTH</surname></name>
    </person-group>
    <article-title>Annual report 2023–24</article-title>
    <publisher-name>Ministry of Road Transport and Highways, Government of India</publisher-name>
    <year>2023</year>
  </element-citation>
</ref>
<ref id="MoRTH2024">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>MoRTH</surname></name>
    </person-group>
    <article-title>Annual report 2024–25</article-title>
    <publisher-name>Ministry of Road Transport and Highways, Government of India</publisher-name>
    <year>2024</year>
  </element-citation>
</ref>
<ref id="Morris1991">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Morris</surname><given-names>M. D.</given-names></name>
    </person-group>
    <article-title>Factorial sampling plans for preliminary computational experiments</article-title>
    <source>Technometrics</source>
    <year>1991</year>
    <volume>33</volume>
    <issue>2</issue>
    <fpage>161</fpage>
    <lpage>174</lpage>
  </element-citation>
</ref>
<ref id="Nicholson1998">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Nicholson</surname><given-names>A.</given-names></name>
    </person-group>
    <article-title>Superelevation, side friction, and roadway consistency</article-title>
    <source>Journal of Transportation Engineering</source>
    <year>1998</year>
    <volume>124</volume>
    <issue>5</issue>
    <fpage>411</fpage>
    <lpage>418</lpage>
  </element-citation>
</ref>
<ref id="Parsa2020">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Parsa</surname><given-names>A. B.</given-names></name>
      <name><surname>Movahedi</surname><given-names>A.</given-names></name>
      <name><surname>Taghipour</surname><given-names>H.</given-names></name>
      <name><surname>Derrible</surname><given-names>S.</given-names></name>
      <name><surname>Mohammadian</surname><given-names>A. K.</given-names></name>
    </person-group>
    <article-title>Toward safer highways: Application of XGBoost and SHAP for real-time accident detection</article-title>
    <source>Accident Analysis &amp; Prevention</source>
    <year>2020</year>
    <volume>136</volume>
    <fpage>105405</fpage>
  </element-citation>
</ref>
<ref id="Pratt2015">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Pratt</surname><given-names>M. P.</given-names></name>
      <name><surname>Geedipally</surname><given-names>S. R.</given-names></name>
      <name><surname>Pike</surname><given-names>A. M.</given-names></name>
    </person-group>
    <article-title>Analysis of vehicle speeds and speed differentials in curves</article-title>
    <source>Transportation Research Record</source>
    <year>2015</year>
    <volume>2486</volume>
    <fpage>28</fpage>
    <lpage>36</lpage>
  </element-citation>
</ref>
<ref id="Saltelli2010">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Saltelli</surname><given-names>A.</given-names></name>
      <name><surname>Annoni</surname><given-names>P.</given-names></name>
      <name><surname>Azzini</surname><given-names>I.</given-names></name>
      <name><surname>Campolongo</surname><given-names>F.</given-names></name>
      <name><surname>Ratto</surname><given-names>M.</given-names></name>
      <name><surname>Tarantola</surname><given-names>S.</given-names></name>
    </person-group>
    <article-title>Variance-based sensitivity analysis of model output</article-title>
    <source>Computer Physics Communications</source>
    <year>2010</year>
    <volume>181</volume>
    <issue>2</issue>
    <fpage>259</fpage>
    <lpage>270</lpage>
  </element-citation>
</ref>
<ref id="Semeida2014">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Semeida</surname><given-names>A. M.</given-names></name>
    </person-group>
    <article-title>Application of artificial neural networks for operating speed prediction at horizontal curves</article-title>
    <source>Journal of Modern Transportation</source>
    <year>2014</year>
    <volume>22</volume>
    <fpage>20</fpage>
    <lpage>29</lpage>
  </element-citation>
</ref>
<ref id="Shalkamy2021">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Shalkamy</surname><given-names>A.</given-names></name>
      <name><surname>Gargoum</surname><given-names>S.</given-names></name>
      <name><surname>El-Basyouny</surname><given-names>K.</given-names></name>
    </person-group>
    <article-title>Towards a more inclusive and safe design of horizontal curves</article-title>
    <source>Accident Analysis &amp; Prevention</source>
    <year>2021</year>
    <volume>153</volume>
    <fpage>106009</fpage>
  </element-citation>
</ref>
<ref id="Singh2017">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Singh</surname><given-names>S. K.</given-names></name>
    </person-group>
    <article-title>Road traffic accidents in India: Issues and challenges</article-title>
    <source>Transportation Research Procedia</source>
    <year>2017</year>
    <volume>25</volume>
    <fpage>4708</fpage>
    <lpage>4719</lpage>
  </element-citation>
</ref>
<ref id="Sobol2001">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Sobol</surname><given-names>I. M.</given-names></name>
    </person-group>
    <article-title>Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates</article-title>
    <source>Mathematics and Computers in Simulation</source>
    <year>2001</year>
    <volume>55</volume>
    <fpage>271</fpage>
    <lpage>280</lpage>
  </element-citation>
</ref>
<ref id="WHO2023">
  <element-citation publication-type="report">
    <person-group person-group-type="author">
      <name><surname>WHO</surname></name>
    </person-group>
    <article-title>World health statistics 2023: Monitoring health for the SDGs</article-title>
    <publisher-name>World Health Organization</publisher-name>
    <year>2023</year>
  </element-citation>
</ref>
<ref id="Yang2022">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Yang</surname><given-names>Y.</given-names></name>
      <name><surname>Wang</surname><given-names>K.</given-names></name>
      <name><surname>Yuan</surname><given-names>Z.</given-names></name>
      <name><surname>Liu</surname><given-names>D.</given-names></name>
    </person-group>
    <article-title>Predicting freeway traffic crash severity using an XGBoost-Bayesian network model</article-title>
    <source>Journal of Advanced Transportation</source>
    <year>2022</year>
  </element-citation>
</ref>
<ref id="Zhao2022">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Zhao</surname><given-names>C.</given-names></name>
      <name><surname>Zhao</surname><given-names>X.</given-names></name>
      <name><surname>Li</surname><given-names>Z.</given-names></name>
      <name><surname>Zhang</surname><given-names>Q.</given-names></name>
    </person-group>
    <article-title>XGBoost-DNN mixed model for predicting driver behaviour during lane-changing decisions</article-title>
    <source>Sustainability</source>
    <year>2022</year>
    <volume>14</volume>
    <fpage>6829</fpage>
  </element-citation>
</ref>
<ref id="Zhou2017">
  <element-citation publication-type="journal">
    <person-group person-group-type="author">
      <name><surname>Zhou</surname><given-names>W.</given-names></name>
      <name><surname>Gong</surname><given-names>C.</given-names></name>
      <name><surname>Hong</surname><given-names>H. P.</given-names></name>
    </person-group>
    <article-title>New perspective on application of first-order reliability method for estimating system reliability</article-title>
    <source>Journal of Engineering Mechanics</source>
    <year>2017</year>
    <volume>143</volume>
    <issue>9</issue>
  </element-citation>
</ref>
</ref-list>
</back>
</article>