Assessment of Climatic and Anthropogenic Controls on Bridge Deck Drainage and Sediment Removal

Bridge deck drainage is essential to prevent hydroplaning and maintain safety along major roadways. With projected changes in climate, current designs may not be sufficient and a better understanding of the primary controls (climate, bridge deck, and inlet design) on the hydraulic efficiency and sediment removal of drainage systems is needed to maintain public safety. To evaluate the controls on hydraulic drainage efficiency, 576 controlled laboratory experiments were conducted testing grate type (rectangular bar vs. curved vane) and downspout configuration (square vs. circular and 20 cm vs. 25 cm) across a range of flow rates, cross slopes, and longitudinal slopes. An additional 144 sediment erosion experiments were performed to identify controls on the removal of sediment. Hydraulic testing indicated that inflow driven by climate is a primary control on drainage efficiency and spread of water on a roadway. For anthropogenic controls, downspout opening size was found to be the primary control followed by longitudinal slope. Sediment removal results indicated that inflow regime and grate type were the primary controls on the sediment removal rate. Given that inflow, driven by climate, is a control on both hydraulic and sediment removal performance, hydraulic engineers should consider forecasted changes in rainfall intensity in their present-day drainage designs. We provide design guidance and discussion for developing a proactive approach to hydraulic infrastructure in the face of future climate uncertainty.


Introduction
Roadway systems are important aspects of the modern world and allow for the interaction between people and communities that prompts economic growth and stability [1]. To protect this function and ensure trust within the public, traffic safety is a top priority for engineers. Traffic safety can be compromised by large precipitation events as the accumulation of water within a roadway increases the risk of hydroplaning, potentially causing the loss of vehicle control and human life [2]. To prevent hydroplaning, drainage systems are designed to convey stormwater runoff away from the roadway surface and limit the spread of accumulated water into traffic lanes, highlighting the importance of drainage system design when considering roadway infrastructure.
Roadways and bridge systems are in crucial need for repair with the condition of over 43% of public roadways in America rated as poor or mediocre [3]. The adequacy of current infrastructure systems to safely convey traffic during rainfall events will only worsen as changes to climate, such as extreme precipitation, temperature, and windstorm events, increase infrastructure stress [4][5][6][7][8]. In particular, climate change will stress roadway drainage by increasing rainfall intensities in much of the world, resulting in more runoff and debris flow [8] and increasing the susceptibility for clogging. There is potential that current drainage system designs lack adequate capacity to handle these forecasted stressors, lowering the safety for commuters on roadways [4,5]. Preparing for these issues now is stressors, lowering the safety for commuters on roadways [4,5]. Preparing for these issues now is crucial to reducing future risks due to climate change [9]. With the current need for updated infrastructure, the opportunity is present to review and update current roadway drainage systems and designs to better prepare for the future.
The appropriate choices regarding how to update roadway or bridge infrastructure depend upon the variables that control ponding of water. For example, roadway design, drainage inlet design, and inlet clogging status all affect the potential for water buildup and hydroplaning. Roadway drainage first conveys sheet flow laterally to the curb where it channelizes as gutter flow and is conveyed longitudinally toward drainage inlets (Figure 1). By capturing and removing gutter flow, drainage inlets limit the spread of water (T) across the deck. Drainage system design is governed by the goal to contain the spread of water within the shoulder area and out of driving lanes [10]. To determine the number and spacing of inlets required for necessary drainage, the projected magnitude of gutter flow and capacity of the inlets to intercept flow are considered [11,12]. Additionally, bridge drainage inlets are often plagued by debris accumulation (Supplementary Materials Table S1), which reduce their ability to convey runoff from the roadway making it also important to understand the impact of design parameters on debris movement. Recent literature [2,13] on bridge deck drainage systems has focused on the implementation and assessment of experimentally derived design equations. We aim to understand the controls on drainage design from a broader perspective by comparing different combination of anthropogenic controls (bridge deck and inlet design parameters) and climate controls to provide holistic design guidance. Additionally, we analyze the impact of anthropogenic and climatic controls on sediment removal rate of drainage inlets which, to our knowledge, has not been examined before. To achieve our aim, we use a scaled physical model to determine the primary controls on bridge deck drainage efficiency and drainage inlet sediment removal by varying climate drivers, bridge deck design parameters, and drainage inlet design parameters. We develop two primary experimental setups. The first experimental setup examines hydraulic efficiency and spread for multiple configurations of gutter flow, roadway slope, and inlet designs to compare standard design performance metrics for design guidance. The second experimental setup tested sediment removal rates with multiple configurations of gutter flows, bridge slopes, and grate designs in order to develop guidelines for preventing the commonly observed field issue of inlet clogging (Supplementary Materials Table S1). In each set of experiments, design factors were examined through related experimental variables. Climate is examined through inflow settings, bridge deck design by model slope, and inlet design by downspout opening size and grate type. The results of this study provide guidance toward effective and Recent literature [2,13] on bridge deck drainage systems has focused on the implementation and assessment of experimentally derived design equations. We aim to understand the controls on drainage design from a broader perspective by comparing different combination of anthropogenic controls (bridge deck and inlet design parameters) and climate controls to provide holistic design guidance. Additionally, we analyze the impact of anthropogenic and climatic controls on sediment removal rate of drainage inlets which, to our knowledge, has not been examined before. To achieve our aim, we use a scaled physical model to determine the primary controls on bridge deck drainage efficiency and drainage inlet sediment removal by varying climate drivers, bridge deck design parameters, and drainage inlet design parameters. We develop two primary experimental setups. The first experimental setup examines hydraulic efficiency and spread for multiple configurations of gutter flow, roadway slope, and inlet designs to compare standard design performance metrics for design guidance. The second experimental setup tested sediment removal rates with multiple configurations of gutter flows, bridge slopes, and grate designs in order to develop guidelines for preventing the commonly observed field issue of inlet clogging (Supplementary Materials Table S1). In each set of experiments, design factors were examined through related experimental variables. Climate is examined through inflow settings, bridge deck design by model slope, and inlet design by downspout opening size and grate type. The results of this study provide guidance toward effective and proactive approaches to ensuring new and existing bridge infrastructure and drainage systems continue to meet traffic safety mandates in the face of a changing climate.

Laboratory Experimental Setup
The primary method to test the goals of the study was through the use of a physical scale model as is commonly used in studies [2,[10][11][12][13] related to roadway drainage. The use of experimental model allowed for directly testing the impact of longitudinal slope, cross slope, inflow, and inlet type on the removal of gutter flow from a roadway surface and width of water on the roadway. Each of these parameters are important for drainage design. An experimental model provides flexibility as multiple combinations of variables can be tested and it allows us to control the quality and quantity of data needed for our analysis. This was especially crucial for analysis of sediment removal rates related to grate design where limited literature exists.

Hydraulic Efficiency Experiment
The first set of experiments was conducted to quantify hydraulic efficiency, E, or the fraction of total gutter flow captured by an inlet as defined by the following equation [10]: where Q c is the flow captured by the inlet (m 3 /s) and Q is the total gutter flow in the channel (m 3 /s). Gutter flow, as described in Figure 1, is defined by the following relationship in design [11,12]: where Q = total gutter flow rate, k g = gutter flow unit conversion factor (0.375 for SI units), n = Manning's roughness coefficient, T = width of flow, S x = cross slope of the roadway, and S = longitudinal slope of the roadway. This parameter is measured in the experiment but shown here for background on the relationship between experimental variables. A single opening in the deck gutter was created to test grate inlet designs (Figure 2b). The inlets and grates were scaled using a modeling software and 3-D printed for precision and control of design (Figure 2c). The four inlet types printed for this study simulated 0.61  A single opening in the deck gutter was created to test grate inlet designs (Figure 2b). The inlets and grates were scaled using a modeling software and 3-D printed for precision and control of design ( Figure 2c). The four inlet types printed for this study simulated 0.61 m square openings. Bar and vane grates were scaled and printed following design drawings provided by Kansas (KDOT) [14] and Illinois (ILDOT) [15] Departments of Transportation. Square and circular downspout shapes were tested for two widths/diameters (20 cm and 25 cm). Upstream inflow to the model was provided by a hose with a calibrated flow meter (Tuthill TT10PN) (Figure 2d) with an accuracy of +/-1.00%. The hose was placed against the curb of the deck facing downslope to immediately simulate channelized gutter flow. Measurement tapes were attached 0.76 m upstream of the inlet to measure the spread of water. To represent the roughness of bridge deck surfaces, the standard Manning's n (0.016) used by KDOT [14] was scaled to 0.012 by application of a length ratio relationship. Sand particles (0.5 mm silica) with an equivalent grain roughness were adhered to the deck. Flow captured by the inlet (Figure 2e) was measured using a scale (Ohaus SD75) and collection bin placed beneath the deck. The scales used in the experiment are professionally calibrated every year and a check was performed before experimentation. The tolerance of each scale was 50 g. The system was automated with solenoid valves to drain the collection bins between experiments. A full breakdown of experimental model and materials can be found in Supplementary Materials Table S2.
The procedure for measuring E dependence on cross slope, longitudinal slope, grate type, and downspout configuration was to run the experimental series listed in Figure 3 at low (0.7 m 3 /h), medium (1.0 m 3 /h), and high (1.6 m 3 /h) inflow regimes with three replicates per configuration. For each trial, the deck surface was wetted by inflow for 30 s prior to data collection to remove the influence of surface material on water spread. Each experimental run lasted three minutes. Starting and ending captured weights and spreads were recorded for each run. Analysis of the results from the hydraulic efficiency experiment compared parameters of efficiency and spread.

Sediment Removal Experiment
The second set of experiments was conducted to measure inlet sediment removal rate based on grate type (vane or bar) and flow intensity. The experimental setup for hydraulic efficiency was modified to extend experiments from 3 to 10 min, allowing for measurable erosion to occur. At the inlet, a cohesive sediment mixture was packed underneath the grate to simulate field conditions ( Figure 4).

Sediment Removal Experiment
The second set of experiments was conducted to measure inlet sediment removal rate based on grate type (vane or bar) and flow intensity. The experimental setup for hydraulic efficiency was modified to extend experiments from 3 to 10 min, allowing for measurable erosion to occur. At the inlet, a cohesive sediment mixture was packed underneath the grate to simulate field conditions ( Figure 4).

Sediment Removal Experiment
The second set of experiments was conducted to measure inlet sediment removal rate based on grate type (vane or bar) and flow intensity. The experimental setup for hydraulic efficiency was modified to extend experiments from 3 to 10 min, allowing for measurable erosion to occur. At the inlet, a cohesive sediment mixture was packed underneath the grate to simulate field conditions ( Figure 4). The sediment mixture consisted of modeling clay, 0.5 mm silica, and water with composition percentages of 30%, 46%, and 24%, respectively. The wet density of the mixture was 2.04 g/cm 3 and the dry density was 1.99 g/cm 3 . The experiments were conducted as described in Figure 3 with three trials per combination. The sediment removal rate (g/s) was measured by recording the loss of sediment mass over the 10 min trials. Spread measurements were recorded at a location upstream and downstream of the inlet. The sediment mixture consisted of modeling clay, 0.5 mm silica, and water with composition percentages of 30%, 46%, and 24%, respectively. The wet density of the mixture was 2.04 g/cm 3 and the dry density was 1.99 g/cm 3 . The experiments were conducted as described in Figure 3 with three trials per combination. The sediment removal rate (g/s) was measured by recording the loss of sediment mass over the 10 min trials. Spread measurements were recorded at a location upstream and downstream of the inlet.

Sediment Mixture Comparison
Measured experimental erosion was contextualized within the literature by analyzing the relationship between dimensionless shear stress and transport. We present our results alongside a collection of natural mix types found in Walder [16]. For cohesive sediment, a dimensionless parameter for shear stress, φ, and a dimensionless transport parameter, R, are defined as in Walder [16]: where ε = the sediment entrainment rate, ρ s = the density of the sediment, ρ = the density of water, τ c = the critical shear stress, u * = the shear velocity, and u cr = the critical shear velocity. Sediment entrainment was determined by van Rijn [17] as ε = M/At where M = the total mass of sediment lost, A = surface area, and ∆t = the measurement period.
The shear velocity was determined as u * = ghS, where g = gravity, h = the hydraulic radius, and S = the longitudinal slope. A critical shear stress of 0.014 Pa was used as determined from experimental testing of the upstream bed shear with a specific gravity of 2.04 and a dimensionless critical shear stress of 0.14. Experimental erosion data by grate type were compared to meta-analysis by Walder [16] for a range of sediment types (Supplementary Materials Figure S1). The results of our created mixture were similar to other natural cohesive sediments as indicated in Walder [16].

Precipitation Analysis
To inform real-world decision making regarding optimal grate type as a function of flow regime, experimental results were applied to environmental conditions from eight sites across Illinois (IL) and Kansas (KS) by analyzing 15 min precipitation data over a 10-year period (2000-2010) [18]. The Kansas gauges were in the towns of Iola, Lawrence, Smolan, and Wilson and the Illinois gauges were in Chicago, Danville, Illinois City, and Quincy. The rainfall data were used to categorize actual events into three experimental inflow regimes (low, medium, and high). Cutoffs were determined through transforming the experimental inflow regimes of 0.7, 1.0, and 1.6 m 3 /h to rainfall intensities (mm/h) by dividing by the deck area resulting in equivalent intensities of 94, 141, and 220 mm/h representing low, medium, and high intensities, respectively. Low-to-medium flow occurrences were defined as the number of events in the range of 51 to 220 mm/h and high flow occurrences were events greater than 220 mm/h.

Statistical Analysis
Statistical analysis of data from both experiments were analyzed using a combination of IBM SPSS Statistics 27 software [19] and SciPy statistical package in Python [20]. All statistical tests used in this paper were performed with a significance level of α = 0.05. Statistical testing for this study incorporated all experimental trials. Results reflect all measurement error and variability; represented as the standard deviation within the results.
For the hydraulic efficiency experiment, our response variables were efficiency and spread. These results were non-normally distributed, as indicated by the Shapiro-Wilks's test; thus, hypothesis testing was performed using the non-parametric Mann-Whitney U-test [21]. One-Way Analysis of Variance (ANOVA) with the Tukey Post hoc test was used to examine the influence of inflow, longitudinal slope, grate type, downspout shape, and downspout size on the performance metrics of efficiency and spread. ANOVA is typically used to difference in means when more than two groups are present as in this study [22,23]. For the sediment removal experiment, the response variable was sediment removal rate. Distribution of erosion results was not normally distributed and statistical significance was evaluated with the Mann-Whitney U-test. One-Way ANOVA with the Tukey post hoc test was applied to analyze the influence of inflow, bridge deck slope, and grate type on sediment removal rate. Data were split by inflow regime to be consistent compared to the analysis of the hydraulic efficiency experiment.

Results
This section examines the results of the two experimental setups: hydraulic efficiency and sediment removal. These two results are explored within the context of three drivers: flow regime (a proxy for climate), bridge deck design, and inlet design. Significant experimental results are highlighted and discussed in each subsection.

Hydraulic Efficiency
Efficiency and spread were compared for 6% cross slope only as measurements at 2% cross slope showed no response to discharge ( Figure 5). The application of the one-way ANOVA indicated inflow regime significantly influenced hydraulic efficiency (F 2,285 = 150.02, p < 0.01). Analysis with the Tukey post hoc test indicated efficiency was significantly higher for the low flow condition (µ = 0.93, σ = 0.06) compared to the medium (µ = 0.85, σ = 0.09) and high (µ = 0.69, σ = 0.13) flow conditions. Inflow regime was also a significant influence (F 2,285 = 154.92, p < 0.01) on the spread of water into the roadway. Significant difference in average spread (cm) was observed for all comparisons of low (µ = 15.33, σ = 1.82), medium (µ = 17.57, σ = 1.82), and high (µ = 19.99, σ = 1.87) inflow groupings. Our results indicate that for bridge decks with large cross slopes, the intensity of gutter flow, which we infer to be a proxy for precipitation intensity and climate, is a primary control on the performance of water capture efficiency and spread. 0.85, σ = 0.09) and high (µ = 0.69, σ = 0.13) flow conditions. Inflow regime was also a significant influence (F2,285 = 154.92, p < 0.01) on the spread of water into the roadway. Significant difference in average spread (cm) was observed for all comparisons of low (µ = 15.33, σ = 1.82), medium (µ = 17.57, σ = 1.82), and high (µ = 19.99, σ = 1.87) inflow groupings. Our results indicate that for bridge decks with large cross slopes, the intensity of gutter flow, which we infer to be a proxy for precipitation intensity and climate, is a primary control on the performance of water capture efficiency and spread.

Sediment Removal
The magnitude of inflow significantly influenced the rate of sediment removal within drainage inlets (F2,141 = 3.93, p = 0.02). Mean sediment removal rate (g/s) at the low inflow regime (µ = 0.047, σ = 0.019) were significantly different from erosion at the medium (µ = 0.065, σ = 0.045) and high (µ = 0.063, σ = 0.035) inflows. Surprisingly, no significant differences were observed between sediment removal rates at the medium and high inflow regimes, which we expected to occur given that high flows produce larger, more erosive

Sediment Removal
The magnitude of inflow significantly influenced the rate of sediment removal within drainage inlets (F 2,141 = 3.93, p = 0.02). Mean sediment removal rate (g/s) at the low inflow regime (µ = 0.047, σ = 0.019) were significantly different from erosion at the medium (µ = 0.065, σ = 0.045) and high (µ = 0.063, σ = 0.035) inflows. Surprisingly, no significant differences were observed between sediment removal rates at the medium and high inflow regimes, which we expected to occur given that high flows produce larger, more erosive shear stresses [16,17]. Potential explanations include non-linear erosion dynamics between flow rate and entrainment, uncertainty with mixture uniformity, measurement error, or the influence of other experimental variables as described below. Our results indicate that greater flow intensities, such as what may be experienced with climate change, cause reduced hydraulic efficiency of water capture; however, flow impact on sediment removal is less certain.

Hydraulic Efficiency
The relationship between efficiency and spread was markedly different for the two cross slopes (Figure 5). At a 6% cross slope, efficiency co-varied with spread, whereas at the 2% cross slope, there was no consistent relation of efficiency in response to spread. The lack of response at 2% cross slope was attributed to a low energy condition, where friction or surface roughness may have controlled the transport of water downslope due to small flow depths [24]. With less gutter flow directly over the inlet openings, the amount of water removed from the opening is less-controlled by the inlet design [25,26]. The design implication of shallow cross slopes is that they perform poorly when it comes to channelizing roadway runoff into gutter flow, limiting the ability of drainage inlets to remove water, thus causing greater spread into the traffic lanes [27]. For shallow cross slopes, such as the 2% tested in this study, the measured spread is large, and the design of shoulder widths must consider whether they will be large enough to contain these high-water depths.
Across most experiments, hydraulic efficiency generally decreased as longitudinal slope increased ( Figure 6). The exception to this observation was at the largest longitudinal slope where we observed an uptick in hydraulic efficiency. With respect to flow, hydraulic efficiency decreased as flow increased irrespective of longitudinal slope. This relationship was significant for all inflow rates: low (F 3,92 = 102.64, p < 0.01), medium (F 3,92 = 21.17, p < 0.01), and high (F 3,92 = 5.57, p < 0.01). Regarding spread of water onto the roadway, we observed significantly lower spread as longitudinal slope increased (Figure 7). Irrespective of longitudinal slope, we found that spread significantly increased with increasing inflow across all flow regimes: low (F 3,92 = 25.06, p < 0.01), medium (F 3,92 = 89.96, p < 0.01), and high (F 3,92 = 65.61, p < 0.01).

Sediment Removal
No statistical difference in mean sediment removal rates (g/s) between the cross slopes of 2% (µ = 0.054, σ = 0.027) and 6% (µ = 0.063, σ = 0.042) was observed. This was Overall, efficiency and spread co-varied for longitudinal slope across all inflow regimes, as shown in previous studies [2,[25][26][27]. However, within our experiment the largest longitudinal slope produced an increase on average in efficiency for the low and medium inflow regimes ( Figure 6) compared to the other studies in which the efficiency decreases as longitudinal slope increased. It is hypothesized that this is attributed to the placement of the inflow opening at the low end of the cross slope where it channelized water more than would be represented in the real world causing less flow bypass around the side of the inlet. These results confirm that the longitudinal slope of the bridge deck is a significant control on the overall performance of bridge deck drainage systems.

Sediment Removal
No statistical difference in mean sediment removal rates (g/s) between the cross slopes of 2% (µ = 0.054, σ = 0.027) and 6% (µ = 0.063, σ = 0.042) was observed. This was also true when statistical testing was applied within inflow regimes. The mean sediment removal rates at the 2% cross slope for low, medium, and high inflow were 0.045 (σ = 0.014), 0.058 (σ = 0.038), and 0.060 (σ = 0.024) compared to rates at 6% cross slope of 0.050 (σ = 0.023), 0.073 (σ = 0.050), and 0.066 (σ = 0.046). Uncertainties in our estimates range from 30 to 70%, and thus no statistically significant differences could be determined. Increasing the number of experiments could help alleviate this issue. Despite this, cross slope should still be considered to have an impact on sediment scour as larger cross-slopes assist in greater channelization of gutter flow, larger longitudinal velocities, and higher shear stresses in the drainage inlet [28].
While longitudinal slope was not a significant influence on sediment removal rates for the low inflow regime (F 3,44 = 5.57, p = 0.78), it was significant at the medium (F 3,44 = 6.04, p < 0.01) and high (F 3,44 = 3.54, p = 0.02) inflow regimes. For the medium inflow regime, mean sediment removal rates (g/s) at the steepest longitudinal slope of 4% (µ = 0.106, σ = 0.070) were significantly different compared to the slopes of 0.5% (µ = 0.060, σ = 0.027), 1% (µ = 0.046, σ = 0.016), and 2% (µ = 0.050, σ = 0.014). However, no significant differences in mean erosion were found comparing within the smaller slopes. The results for the longitudinal slope indicated that, on average, the sediment removal rates at a slope of 4% were 50% greater compared to the other three tested slopes. Since the flow would mainly be in the streamwise direction due to the steep slope gradient, a greater increase in bed shear over the sediment element occurs to remove particles from the inlet at a greater rate [29]. The use of steeper longitudinal slopes, when possible, could alleviate issues related to debris in bridge systems.

Hydraulic Efficiency
Hydraulic efficiency was similar for both grate types (bar and vane) at low and medium flow regimes, but performance differed at the highest inflow rate. Average efficiency for the bar grate across low, medium, and high inflows was 0.93 (σ =0.07), 0.86 (σ = 0.07), and 0.74 (σ = 0.10), respectively, whereas vane grate values were 0.93 (σ = 0.06), 0.84 (σ = 0.09), and 0.65 (σ = 0.14), respectively. The bar grate design was significantly more efficient at the high inflow rate compared to the vane design. Statistical comparison of spread between grate types were similar with average values at low, medium, and high inflows of 15.28 (σ = 1.67), 17.54 (σ = 1.69), and 20.11 (σ =1.86) for the bar grate and 15.38 (σ =1.97), 17.59 (σ = 1.97), and 19.89 (σ = 1.90) for the vane grate, respectively. The larger vertical opening area of the bar grate compared to the vane rate allowed for greater flow interception at the high inflow regime explaining the difference in efficiency. This indicates that designers must consider grate type as a possible control on water drainage, particularly at high flows.
The downspout shape and width were also identified as controls on efficiency for medium and high flow, but not for low flow (Figure 8). For square downspouts, average efficiencies for low, medium, and high inflow were 0.94 (σ = 0.05), 0.87 (σ = 0.07), and 0.73 (σ = 0.11), respectively, whereas for circular downspouts they were 0.92 (σ = 0.07), 0.82 (σ = 0.09), and 0.66 (σ = 0.13), respectively. Similar relationships were observed for downspout widths (20 cm vs. 25 cm). Average efficiencies for 20 cm downspouts were 0.92 (σ = 0.07), 0.82 (σ = 0.09), and 0.61 (σ = 0.11) for low, medium, and high inflow, respectively, compared to 0.94 (σ = 0.06), 0.87 (σ = 0.07), and 0.73 (σ = 0.11) for 25 cm downspouts. In each case the design with the larger average cross-sectional area (square shape and/or 25 cm), performed better indicating downspout area was an important control. At the high inflow regime, mean efficiency was significantly different between all opening areas ( Figure 8). Application of larger downspout areas lessens the gap in efficiency between inflow regimes indicating that inclusions of larger opening areas in design would be beneficial for maintaining efficiency across all precipitation intensities. Downspout size and shape were a significant influence for the low (F3,92 = 10.67, p < 0.01), medium (F3,92 = 3.66, p = 0.02), and high (F3,92 = 4.67, p < 0.01) inflow regimes ( Figure  9). For low inflow regimes, the smallest downspout area had significantly larger spreads compared to the three larger areas. For the medium inflow, spread gradually decreases as the opening area increases with significant difference between areas of 3.88 and 7.72 cm 2 . A transition region is likely occurring where the spread is partially controlled by downspout design, but the largest size alleviates the control condition. Similarly, for the high inflow regime in Figure 9, significant difference is only observed between the largest opening area and the two smallest sizes. Overall, the larger downspout design is beneficial to controlling spread across all conditions, and thus reducing the risk of hydroplaning. Downspout size and shape were a significant influence for the low (F 3,92 = 10.67, p < 0.01), medium (F 3,92 = 3.66, p = 0.02), and high (F 3,92 = 4.67, p < 0.01) inflow regimes (Figure 9). For low inflow regimes, the smallest downspout area had significantly larger spreads compared to the three larger areas. For the medium inflow, spread gradually decreases as the opening area increases with significant difference between areas of 3.88 and 7.72 cm 2 . A transition region is likely occurring where the spread is partially controlled by downspout design, but the largest size alleviates the control condition. Similarly, for the high inflow regime in Figure 9, significant difference is only observed between the largest opening area and the two smallest sizes. Overall, the larger downspout design is beneficial to controlling spread across all conditions, and thus reducing the risk of hydroplaning.

Sediment Removal
Sediment removal rates between bar and vane grates differed based on inflow regime (Table 1). In general, sediment removal rate tends to increase with flow rate. This was true for all bar grate flow scenarios; however, this was not the case for vane grates as erosion increased initially from low to medium inflows, but then decreased at high flows (Table 1). No significant differences between grate types at each inflow regime were found, but we mainly attributed this to the small experimental sample size per inflow (n = 24).

FOR PEER REVIEW
12 of 17 Figure 9. Spread measurement boxplot based on inflow category comparing across downspout opening area. (*) Indicates statistical similarity (p > 0.05) between means with numbers corresponding to opening size and shape given on the x-axis.

Sediment Removal
Sediment removal rates between bar and vane grates differed based on inflow regime (Table 1). In general, sediment removal rate tends to increase with flow rate. This was true for all bar grate flow scenarios; however, this was not the case for vane grates as erosion increased initially from low to medium inflows, but then decreased at high flows (Table  1). No significant differences between grate types at each inflow regime were found, but we mainly attributed this to the small experimental sample size per inflow (n = 24). Table 1. Sediment removal rates (± 1 σ) by grate design and average inflow (n = 24 for each inflow).

Sediment Removal Rate (g/s) Bar Grate
Vane Grate 0.7 0.045 ± 0.015 0.049 ± 0.022 1.0 0.056 ± 0.023 0.075 ± 0.059 1.6 0.067 ± 0.045 0.059 ± 0.023 At low-to-medium flows, it is hypothesized that the vane grate is able to guide the water into the inlet as its design has curved bars, which transition from horizontal to vertical, aiding gravity in changing the direction of gutter flow to improve sediment removal rates. Observations from video recordings of the experiments found that, at low-to-medium flows, water could be seen traveling along the vane and diverted into the inlet along the vertical side of the cohesive sediment more efficiently compared to the bar grate that  At low-to-medium flows, it is hypothesized that the vane grate is able to guide the water into the inlet as its design has curved bars, which transition from horizontal to vertical, aiding gravity in changing the direction of gutter flow to improve sediment removal rates. Observations from video recordings of the experiments found that, at low-to-medium flows, water could be seen traveling along the vane and diverted into the inlet along the vertical side of the cohesive sediment more efficiently compared to the bar grate that allowed more flow to bypass horizontally over the sediment material, supporting the results in Table 1. However, at higher flows, when grates are submerged, it is hypothesized that the curvature of the vanes may cause an obstruction in the vertical flow path of water toward the downspout pipe, limiting erosion potential. Grate-splashover velocity relationships found in Johnson and Chang [30] and Brown et al. [31] support this idea as a vane grate was shown to have more splashover compared to a bar grate for large gutter flows. This explanation is plausible as the rectangular bar grate, with its bars parallel to potential vertical flow paths, had improved cleanout over the curved vane grate at high flows. The experiment shows how bridge configuration and inflow can affect the movement of debris based on grate type.

Controls and Implications
Previous works on roadway drainage [2,[11][12][13]25,26] focus on the empirical relationship between efficiency and spread based on the controls tested in this study. We considered three types of controls for hydraulic performance of water and sediment removal: climate, bridge design, and drainage design to define the importance of each control for a more robust discussion. Inflow intensity, a proxy for climate, had significant impacts to hydraulic efficiency and spread for nearly all experiments, regardless of the evaluated parameter, e.g., longitudinal slope (Figures 6 and 7), and downspout area (Figures 8 and 9). However, inflow intensity is not a variable that can be controlled, but rather must be designed for while considering climate uncertainties. The two variables hydraulic engineers have some control over include bridge (cross and longitudinal deck slope) and drainage (grate and downspout configuration) design. Bridge design is typically not selected by hydraulic engineers, but rather by structural engineers. However, this variable was the primary control of measured spreads and efficiencies ( Figure 5). In situations where structural engineers have flexibility with deck design, increased cross slope could help with improving hydraulic efficiency and reducing roadway spread. Our study indicates that low cross slopes drive the design based on surface roughness control of the water. Current literature is lacking in this regard as no known studies investigate this scenario indicating that intermediate cross slopes should also be investigated to understand the threshold in which this condition subsides. Additionally, longitudinal slope influenced design across all inflow regimes, which is consistent with previous studies [2,13] and should be examined alongside cross slope when feasible. On the other hand, the shape and size of the downspout (drainage design) has significant performance implications and is specified by hydraulic engineers. We found that hydraulic engineers should choose square downspout rather than circular downspout designs. A square downspout of a given width will perform better, particularly during medium and high flows (Figures 8 and 9), as a circular design of an equivalent diameter, predominantly due to larger opening area despite having a similar installation footprint. Typically, related studies [2,11,13,25,26] solely focus on a single inlet design, but we suggest further investigation should be coupled with our experimental results across more designs as various types are used in the practice.
Understanding the interactions between bridge deck design controls is important for designing safe and resilient drainage systems as there is the potential for increased rainfall intensities in the future due to climate change [4]. To illustrate, a study by Moustakis et al. [32] examined precipitation extremes and indicated that in the Midwest and Eastern United States high intensity storm events are expected to occur more frequently, i.e., what is now a 20-year event will become a 7-year event for 75% of the US. Designers must understand how their area will respond to future climate scenarios to ensure standard safety factors are met. Implementation of proactive design such as the examination of bridge slopes in the context of drainage and increasing the typical downspout size will help to account for these changes in rainfall intensity. Alongside implementing the design relationships found in this study, designers should also consider different measures for determination of design storm events. Currently, designs within the context of water engineering applies intensity-duration-frequency curves derived from series of annual maximum and partial duration series [6]. However, this approach relies on historical data to assess a return period based on stationarity [6]. A study by Mailhot and Duchesne [9] determined that due to climate change, stationarity is not satisfied. Instead, for the context of water resource design, they [6] proposed a simplified approach where a critical return period is based on design parameters to account for future climate extremes acting as a safety factor to increase the standard return period. Applying regional specific information with newer guidelines will allow for designers to account for future climate uncertainty.
The results of our study not only provide strong justification for changes in drainage design but highlight important metrics to consider when analyzing current infrastructure resiliency and sustainability. Current literature [33] finds that poor road drainage systems have devastating impacts the economy, as structure failures due to poor drainage increasingly impact society through increased travel time and decrease roadway safety. Barreiro et al. (2021) [34] analyze the resilience to flooding scenarios in cities through 1-D and 2-D drainage models with respect to volume, flooding duration, flooded area, and affected infrastructure. Incorporation of our experimental results for cross slope, longitudinal slope, inflow, and inlet type can help better model these systems through portraying the important parameters for runoff removal. For instance, since inflow is an important parameter for drainage efficiency, our results would highlight that multiple simulation types for each inflow regime should be conducted to understand the resilience of systems. Additionally, our findings indicate how downspout configurations that impact drainage efficiency, highlight how models must incorporate the designs found in roadway infrastructure when examining flooding resilience for cities. Understanding the construction and design aspects of when a system was constructed is also crucial to understanding the impact of each system component on the sustainability of a network. Based on our findings, it is important to understand how inlet design and bridge deck design have changed overtime to be able to better understand where changes in roadway drainage need to be made to keep the drainage systems operating at necessary levels.

Sediment-Precipitation Analysis
One of the main applications of this work is that the choice of drainage design (i.e., grate and downspout specifications) impacts the flow and sediment dynamics, which are influenced by the geographic location of the bridge. Our results highlight this idea as sediment removal was shown to change based on grate type and inflow regime (Table 1). To illustrate this point, precipitation analysis across eight sites in Kansas and Illinois were used to select a recommended grate type based on expected storm intensities. Each site's unique storm intensity characteristics resulted in different recommendations for the type of grate to use in design. For the four sites in Kansas, events falling within the high-flow regime are the predominant storm types (Table 2), thus suggesting that bar grates are more suitable for effective design given that our experiments showed they performed better under these conditions (Table 1). For the Illinois sites of Chicago, Danville, and Illinois City, low-to-medium intensity events are temporally more common; thus, designers should utilize curved vane grates as they perform best under low and medium inflow rates. For Quincy, IL, the breakdown of events indicate that it would be advantageous to use the bar grate as it could more efficiently remove sediment as high intensity events are common. However, designers should note that surface runoff generation is directly associated with temporal patterns in rainfall intensity, which drives sediment yield through rainfall impact and scour [35]. It should also be noted that precipitation might not be the sole driver of sediment removal as sediment mixture is location dependent. Further investigation into debris-grate relationships should examine the impact of location specific soil type, land use, and traffic patterns to address other factors that may control erosion potential. Overall, with the creation of debris removal data, designers could easily integrate precipitation analysis coupled with climate change adjustment [6] to address issue of clogging in the field.
Further studies are needed to understand relationships between grate types and sediment transport to establish specific design guidance as well as examine future trends related to climate change. Sediment buildup and transport processes are highly variable with spatial heterogeneities across drainage areas [36], which need to be explored further for developed structures related to precipitation extremes. Along with determining sediment removal relationships, similar experiments can also be used to quantify a clogging coefficient to provide optimized safety factors for design based on debris patterns as performed by Guo and Mackenzie [37] and Gomez et al. (2019) [38]. This would allow designers to better predict field values of efficiency especially when coupled with specific design efficiency curves. Overall, accounting for debris removal would alleviate costs due to maintenance and bridge deck degradation, which can help to prolong bridge structures for the future and allow for better bridge performance [39,40].

Conclusions
In bridge deck drainage design, the primary goal is to reduce hydroplaning by removing runoff generated during storm events. Our results showed that downspout size and shape are the most practical design elements to improve flow removal, particularly during high flow rates. Incorporating larger-than-presently-necessary downspout designs is a simple solution that can combat expected future increases in precipitation intensity and frequency. With design of new infrastructure, the inclusion of hydraulic design within the structural and transportation design process could be important for future design policy as experiments results find that under low-medium inflow conditions design is controlled by bridge deck slope. Additionally, sediment removal experiments indicate that larger bridge deck slopes have more potential to erode debris build-up. Furthermore, curved vane grates had greater sediment cleanout rates at lower flow rates, whereas rectangular bar grates performed better at high inflows. Comparing the experimental inflows with rainfall over a ten-year period at eight sites showed how the distribution of precipitation events can be compared to help decide between grate types if debris removal performance is known. Finally, our data and analysis on controls of drainage efficiency and sediment removal rates provide researchers and practitioners guidance on metrics important to modeling roadway drainage, planning for future infrastructure needs, and assessing current resilience to flood of roadway structures.
The primary limitations identified in our study stem from our experimental model. The first limitation is we only examined two cross-slopes and four longitudinal slopes as model scale, materials, and available time did not allow to test more combinations of slopes. We recommend the testing of a wider range of cross slopes to determine when it becomes a primary control on intercepted flow compared to inflow. A greater density of longitudinal slopes should be tested to cover all available design conditions to better quantify the relationship between intercepted flow and spread width. Further analysis of sediment removal needs to be conducted to develop sediment removal design tables as well as examine future precipitation projections. We only analyzed a single cohesive mixture, but mixture types based on regional materials should be analyzed to understand how sediment removal rates relate to regional climates. Further research on debris in drainage inlets should couple sediment removal with sediment buildup to understand the process holistically. Overall, this study identifies the primary controls on the bridge deck drainage design through experimental results and analysis to facilitate proactive infrastructure upgrades in face of changing precipitation regimes.