Stochastic analysis of transit route segments' passenger load variation for capacity and quality of service assessment

This paper investigates stochastic analysis of transit segment hourly passenger load factor variation for transit capacity and quality of service (QoS) analysis using Automatic Fare Collection data for a premium radial bus route in Brisbane, Australia. It compares stochastic analysis to traditional peak hour factor (PHF) analysis to gain further insight into variability of transit route segments’ passenger loading during a study hour. It demonstrates that hourly design load factor is a useful method of modeling a route segment’s capacity and QoS time history across the study weekday. This analysis method is readily adaptable to different passenger load standards by adjusting design percentile, reflecting either a more relaxed or more stringent condition. This paper also considers hourly coefficient of variation of load factor as a capacity and QoS assessment measure, in particular through its relationships with hourly average and design load factors. Smaller value reflects uniform passenger loading, which is generally indicative of well dispersed passenger boarding demands and good schedule maintenance. Conversely, higher value may be indicative of pulsed or uneven passenger boarding demands, poor schedule maintenance, and/or bus bunching. An assessment table based on hourly coefficient of variation of load factor is developed and applied to this case study. Inferences are drawn for a selection of study hours across the weekday studied.

Measuring route passenger load diversity for capacity and quality of service assessment

This poster introduces Passenger Relative Load Factor for a route or individual bus service as a capacity and quality of service measure, distinguishing it from Occupancy Load Factor. It introduces Load Diversity Coefficient as the ratio of Passenger Relative Load Factor to Occupancy Load Factor, and relates Load Diversity Coefficient to Coefficient of Variation in Occupancy Load Factor. It qualifies the operator’s and passengers’ perspectives of load factor based on Coefficient of Variation in Occupancy Load Factor along a route. A case study using weekday Automatic Fare Collection (AFC) data on a premium bus line in Brisbane, Australia illustrates the methodology. The compendium paper also qualifies the operator’s and passengers’ perspectives of these load factors along with Passengers’ Average Travel Time for capacity and quality of service assessment.

Ultimate flexural limit states analysis of prestressed concrete sleeper

Railway is one of the most important, reliable and widely used means of transportation, carrying freight, passengers, minerals, grains, etc. Thus, research on railway tracks is extremely important for the development of railway engineering and technologies. The safe operation of a railway track is based on the railway track structure that includes rails, fasteners, pads, sleepers, ballast, subballast and formation. Sleepers are very important components of the entire structure and may be made of timber, concrete, steel or synthetic materials. Concrete sleepers were first installed around the middle of last century and currently are installed in great numbers around the world. Consequently, the design of concrete sleepers has a direct impact on the safe operation of railways. The "permissible stress" method is currently most commonly used to design sleepers. However, the permissible stress principle does not consider the ultimate strength of materials, probabilities of actual loads, and the risks associated with failure, all of which could lead to the conclusion of cost-ineffectiveness and over design of current prestressed concrete sleepers. Recently the limit states design method, which appeared in the last century and has been already applied in the design of buildings, bridges, etc, is proposed as a better method for the design of prestressed concrete sleepers. The limit states design has significant advantages compared to the permissible stress design, such as the utilisation of the full strength of the member, and a rational analysis of the probabilities related to sleeper strength and applied loads. This research aims to apply the ultimate limit states design to the prestressed concrete sleeper, namely to obtain the load factors of both static and dynamic loads for the ultimate limit states design equations. However, the sleepers in rail tracks require different safety levels for different types of tracks, which mean the different types of tracks have different load factors of limit states design equations. Therefore, the core tasks of this research are to find the load factors of the static component and dynamic component of loads on track and the strength reduction factor of the sleeper bending strength for the ultimate limit states design equations for four main types of tracks, i.e., heavy haul, freight, medium speed passenger and high speed passenger tracks. To find those factors, the multiple samples of static loads, dynamic loads and their distributions are needed. In the four types of tracks, the heavy haul track has the measured data from Braeside Line (A heavy haul line in Central Queensland), and the distributions of both static and dynamic loads can be found from these data. The other three types of tracks have no measured data from sites and the experimental data are hardly available. In order to generate the data samples and obtain their distributions, the computer based simulations were employed and assumed the wheel-track impacts as induced by different sizes of wheel flats. A valid simulation package named DTrack was firstly employed to generate the dynamic loads for the freight and medium speed passenger tracks. However, DTrack is only valid for the tracks which carry low or medium speed vehicles. Therefore, a 3-D finite element (FE) model was then established for the wheel-track impact analysis of the high speed track. This FE model has been validated by comparing its simulation results with the DTrack simulation results, and with the results from traditional theoretical calculations based on the case of heavy haul track. Furthermore, the dynamic load data of the high speed track were obtained from the FE model and the distributions of both static and dynamic loads were extracted accordingly. All derived distributions of loads were fitted by appropriate functions. Through extrapolating those distributions, the important parameters of distributions for the static load induced sleeper bending moment and the extreme wheel-rail impact force induced sleeper dynamic bending moments and finally, the load factors, were obtained. Eventually, the load factors were obtained by the limit states design calibration based on reliability analyses with the derived distributions. After that, a sensitivity analysis was performed and the reliability of the achieved limit states design equations was confirmed. It has been found that the limit states design can be effectively applied to railway concrete sleepers. This research significantly contributes to railway engineering and the track safety area. It helps to decrease the failure and risks of track structure and accidents; better determines the load range for existing sleepers in track; better rates the strength of concrete sleepers to support bigger impact and loads on railway track; increases the reliability of the concrete sleepers and hugely saves investments on railway industries. Based on this research, many other bodies of research can be promoted in the future. Firstly, it has been found that the 3-D FE model is suitable for the study of track loadings and track structure vibrations. Secondly, the equations for serviceability and damageability limit states can be developed based on the concepts of limit states design equations of concrete sleepers obtained in this research, which are for the ultimate limit states.

Glycaemic load is associated with insulin resistance in older Australian women

Background: Diets with a high postprandial glycemic response may contribute to long-term development of insulin resistance and diabetes, however previous epidemiological studies are conflicting on whether glycemic index (GI) or glycemic load (GL) are dietary factors associated with the progression. Our objectives were to estimate GI and GL in a group of older women, and evaluate cross-sectional associations with insulin resistance. Subjects and Methods: Subjects were 329 Australian women aged 42-81 years participating in year three of the Longitudinal Assessment of Ageing in Women (LAW). Dietary intakes were assessed by diet history interviews and analysed using a customised GI database. Insulin resistance was defined as a homeostasis model assessment (HOMA) value of >3.99, based on fasting blood glucose and insulin concentrations. Results: GL was significantly higher in the 26 subjects who were classified as insulin resistant compared to subjects who were not (134±33 versus 114±24, P<0.001). In a logistic regression model, an increment of 15 GL units increased the odds of insulin resistance by 2.09 (95%CI 1.55, 2.80, P<0.001) independently of potential confounding variables. No significant associations were found when insulin resistance was assessed as a continuous variable. Conclusions: Results of this cross-sectional study support the concept that diets with a higher GL are associated with increased risk of insulin resistance. Further studies are required to investigate whether reducing glycemic intake, by either consuming lower GI foods and/or smaller serves of carbohydrate, can contribute to a reduction in development of insulin resistance and long-term risk of type 2 diabetes.

Post-tensioning and its effect on multi-level formwork load distribution

Multi-level concrete buildings requrre substantial temporary formwork structures to support the slabs during construction. The primary function of this formwork is to safely disperse the applied loads so that the slab being constructed, or the portion of the permanent structure already constructed, is not overloaded. Multi-level formwork is a procedure in which a limited number of formwork and shoring sets are cycled up the building as construction progresses. In this process, each new slab is supported by a number of lower level slabs. The new slab load is, essentially, distributed to these supporting slabs in direct proportion to their relative stiffness. When a slab is post-tensioned using draped tendons, slab lift occurs as a portion of the slab self-weight is balanced. The formwork and shores supporting that slab are unloaded by an amount equivalent to the load balanced by the post-tensioning. This produces a load distribution inherently different from that of a conventionally reinforced slab. Through , theoretical modelling and extensive on-site shore load measurement, this research examines the effects of post-tensioning on multilevel formwork load distribution. The research demonstrates that the load distribution process for post-tensioned slabs allows for improvements to current construction practice. These enhancements include a shortening of the construction period; an improvement in the safety of multi-level form work operations; and a reduction in the quantity of form work materials required for a project. These enhancements are achieved through the general improvement in safety offered by post-tensioning during the various formwork operations. The research demonstrates that there is generally a significant improvement in the factors of safety over those for conventionally reinforced slabs. This improvement in the factor of safety occurs at all stages of the multi-level formwork operation. The general improvement in the factors of safety with post-tensioned slabs allows for a shortening of the slab construction cycle time. Further, the low level of load redistribution that occurs during the stripping operations makes post-tensioned slabs ideally suited to reshoring procedures. Provided the overall number of interconnected levels remains unaltered, it is possible to increase the number of reshored levels while reducing the number of undisturbed shoring levels without altering the factors of safety, thereby, reducing the overall quantity of formwork and shoring materials.

Property damage crash equivalency factors to solve crash frequency-severity dilemma: Case study on South Korean rural roads

Safety interventions (e.g., median barriers, photo enforcement) and road features (e.g., median type and width) can influence crash severity, crash frequency, or both. Both dimensions—crash frequency and crash severity—are needed to obtain a full accounting of road safety. Extensive literature and common sense both dictate that crashes are not created equal, with fatalities costing society more than 1,000 times the cost of property damage crashes on average. Despite this glaring disparity, the profession has not unanimously embraced or successfully defended a nonarbitrary severity weighting approach for analyzing safety data and conducting safety analyses. It is argued here that the two dimensions (frequency and severity) are made available by intelligently and reliably weighting crash frequencies and converting all crashes to property-damage-only crash equivalents (PDOEs) by using comprehensive societal unit crash costs. This approach is analogous to calculating axle load equivalents in the prediction of pavement damage: for instance, a 40,000-lb truck causes 4,025 times more stress than does a 4,000-lb car and so simply counting axles is not sufficient. Calculating PDOEs using unit crash costs is the most defensible and nonarbitrary weighting scheme, allows for the simple incorporation of severity and frequency, and leads to crash models that are sensitive to factors that affect crash severity. Moreover, using PDOEs diminishes the errors introduced by underreporting of less severe crashes—an added benefit of the PDOE analysis approach. The method is illustrated with rural road segment data from South Korea (which in practice would develop PDOEs with Korean crash cost data).

Factors of process model comprehension : findings from a series of experiments

In order to make good decisions about the design of information systems, an essential skill is to understand process models of the business domain the system is intended to support. Yet, little knowledge to date has been established about the factors that affect how model users comprehend the content of process models. In this study, we use theories of semiotics and cognitive load to theorize how model and personal factors influence how model viewers comprehend the syntactical information of process models. We then report on a four-part series of experiments, in which we examined these factors. Our results show that additional semantical information impedes syntax comprehension, and that theoretical knowledge eases syntax comprehension. Modeling experience further contributes positively to comprehension efficiency, measured as the ratio of correct answers to the time taken to provide answers. We discuss implications for practice and research.

A tissue engineering solution for segmental defect regeneration in load-bearing long bones

The reconstruction of large defects (>10 mm) in humans usually relies on bone graft transplantation. Limiting factors include availability of graft material, comorbidity, and insufficient integration into the damaged bone. We compare the gold standard autograft with biodegradable composite scaffolds consisting of medical-grade polycaprolactone and tricalcium phosphate combined with autologous bone marrow-derived mesenchymal stem cells (MSCs) or recombinant human bone morphogenetic protein 7 (rhBMP-7). Critical-sized defects in sheep - a model closely resembling human bone formation and structure - were treated with autograft, rhBMP-7, or MSCs. Bridging was observed within 3 months for both the autograft and the rhBMP-7 treatment. After 12 months, biomechanical analysis and microcomputed tomography imaging showed significantly greater bone formation and superior strength for the biomaterial scaffolds loaded with rhBMP-7 compared to the autograft. Axial bone distribution was greater at the interfaces. With rhBMP-7, at 3 months, the radial bone distribution within the scaffolds was homogeneous. At 12 months, however, significantly more bone was found in the scaffold architecture, indicating bone remodeling. Scaffolds alone or with MSC inclusion did not induce levels of bone formation comparable to those of the autograft and rhBMP-7 groups. Applied clinically, this approach using rhBMP-7 could overcome autograft-associated limitations.

Human papillomavirus prevalence, viral load and pre-cancerous lesions of the cervix in women initiating highly active antiretroviral therapy in South Africa : a cross-sectional study

Background Cervical cancer and infection with human immunodeficiency virus (HIV) are both important public health problems in South Africa (SA). The aim of this study was to determine the prevalence of cervical squamous intraepithelial lesions (SILs), high-risk human papillomavirus (HR-HPV), HPV viral load and HPV genotypes in HIV positive women initiating anti-retroviral (ARV) therapy. Methods A cross-sectional survey was conducted at an anti-retroviral (ARV) treatment clinic in Cape Town, SA in 2007. Cervical specimens were taken for cytological analysis and HPV testing. The Digene Hybrid Capture 2 (HC2) test was used to detect HR-HPV. Relative light units (RLU) were used as a measure of HPV viral load. HPV types were determined using the Roche Linear Array HPV Genotyping test. Crude associations with abnormal cytology were tested and multiple logistic regression was used to determine independent risk factors for abnormal cytology. Results The median age of the 109 participants was 31 years, the median CD4 count was 125/mm3, 66.3% had an abnormal Pap smear, the HR-HPV prevalence was 78.9% (Digene), the median HPV viral load was 181.1 RLU (HC2 positive samples only) and 78.4% had multiple genotypes. Among women with abnormal smears the most prevalent HR-HPV types were HPV types 16, 58 and 51, all with a prevalence of 28.5%. On univariate analysis HR-HPV, multiple HPV types and HPV viral load were significantly associated with the presence of low and high-grade SILs (LSIL/HSIL). The multivariate logistic regression showed that HPV viral load was associated with an increased odds of LSIL/HSIL, odds ratio of 10.7 (95% CI 2.0 – 57.7) for those that were HC2 positive and had a viral load of ≤ 181.1 RLU (the median HPV viral load), and 33.8 (95% CI 6.4 – 178.9) for those that were HC2 positive with a HPV viral load > 181.1 RLU. Conclusion Women initiating ARVs have a high prevalence of abnormal Pap smears and HR-HPV. Our results underscore the need for locally relevant, rigorous screening protocols for the increasing numbers of women accessing ARV therapy so that the benefits of ARVs are not partially offset by an excess risk in cervical cancer.

Risk factors for moderate and severe microbial keratitis in daily wear contact lens users

Objective: To establish risk factors for moderate and severe microbial keratitis among daily contact lens (CL) wearers in Australia. Design: A prospective, 12-month, population-based, case-control study. Participants: New cases of moderate and severe microbial keratitis in daily wear CL users presenting in Australia over a 12-month period were identified through surveillance of all ophthalmic practitioners. Case detection was augmented by record audits at major ophthalmic centers. Controls were users of daily wear CLs in the community identified using a national telephone survey. Testing: Cases and controls were interviewed by telephone to determine subject demographics and CL wear history. Multiple binary logistic regression was used to determine independent risk factors and univariate population attributable risk percentage (PAR%) was estimated for each risk factor.; Main Outcome Measures: Independent risk factors, relative risk (with 95% confidence intervals [CIs]), and PAR%. Results: There were 90 eligible moderate and severe cases related to daily wear of CLs reported during the study period. We identified 1090 community controls using daily wear CLs. Independent risk factors for moderate and severe keratitis while adjusting for age, gender, and lens material type included poor storage case hygiene 6.4× (95% CI, 1.9-21.8; PAR, 49%), infrequent storage case replacement 5.4× (95% CI, 1.5-18.9; PAR, 27%), solution type 7.2× (95% CI, 2.3-22.5; PAR, 35%), occasional overnight lens use (<1 night per week) 6.5× (95% CI, 1.3-31.7; PAR, 23%), high socioeconomic status 4.1× (95% CI, 1.2-14.4; PAR, 31%), and smoking 3.7× (95% CI, 1.1-12.8; PAR, 31%). Conclusions: Moderate and severe microbial keratitis associated with daily use of CLs was independently associated with factors likely to cause contamination of CL storage cases (frequency of storage case replacement, hygiene, and solution type). Other factors included occasional overnight use of CLs, smoking, and socioeconomic class. Disease load may be considerably reduced by attention to modifiable risk factors related to CL storage case practice.

Dynamic load-sharing of longitudinal-connected air suspensions of a tri-axle semi-trailer

The effects of suspension parameters and driving conditions on dynamic load-sharing of longitudinal-connected air suspensions of a tri-axle semi-trailer are investigated in this study. A novel nonlinear model of a multi-axle semi-trailer with longitudinal-connected air suspensions is formulated based on fluid mechanics and thermodynamics and validated through test results. The effects of road surface conditions, driving speeds, air line inside diameter and connector inside diameter on dynamic load-sharing capability of the semi-trailer were analyzed in terms of load-sharing criteria. Simulation results indicate that, when larger air lines and connectors are employed, the DLSC (Dynamic Load-Sharing Coefficient) optimization ratio reaches its peak value when the road roughness is medium. The optimization ratio fluctuates in a complex manner as driving speed increases. The results also indicate that if the air line inside diameter is always assumed to be larger than the connector inside diameter, the influence of air line inside diameter on load-sharing is more significant than that of the connector inside diameter. The proposed approach can be used for further study of the influence of additional factors (such as vehicle load, static absolute air pressure and static height of air spring) on load-sharing and the control methods for multi-axle air suspensions with longitudinal air line.

Transport and health in China : recent collaborative research relating to vehicle factors and pollution

The nature of the transport system contributes to public health outcomes in a range of ways. The clearest contribution to public health is in the area of traffic crashes, because of their direct impact on individual death and disability and their direct costs to the health system. Other papers in this conference address these issues. This paper outlines some collaborative research between the Centre for Accident Research and Road Safety - Queensland (CARRS-Q) at QUT and Chinese researchers in areas that have indirect health impacts. Heavy vehicle dynamics: The integrity of the road surface influences crash risk, with ruts, pot-holes and other forms of road damage contributing to increased crash risks. The great majority of damage to the road surface from vehicles is caused by heavy trucks and buses, rather than cars or smaller vehicles. In some cases this damage is due to deliberate overloading, but in other cases it is due to vehicle suspension characteristics that lead to occasional high loads on particular wheels. Together with a visiting researcher and his colleagues, we have used both Queensland and Chinese data to model vehicle suspension systems that reduce the level of load, and hence the level of road damage and resulting crash risk(1-5). Toll worker exposure to vehicle emissions: The increasing construction of highways in China has also involved construction of a large number of toll roads. Tollbooth workers are potentially exposed to high levels of pollutants from vehicles, however the extent of this exposure and how it relates to standards for exposure are not well known. In a study led by a visiting researcher, we conducted a study to model these levels of exposure for a tollbooth in China(6). Noise pollution: The increasing presence of high speed roads in China has contributed to an increase in noise levels. In this collaborative study we modelled noise levels associated with a freeway widening near a university campus, and measures to reduce the noise(7). Along with these areas of research, there are many other areas of transport with health implications that are worthy of exploration. Traffic, noise and pollution contribute to a difficult environment for pedestrians, especially in an ageing society where there are health benefits to increasing physical activity. By building on collaborations such as those outlined, there is potential for a contribution to improved public health by addressing transport issues such as vehicle factors and pollution, and extending the research to other areas of travel activity. 1. Chen, Y., He, J., King, M., Chen, W. and Zhang, W. (2014). Stiffness-damping matching method of an ECAS system based on LQG control. Journal of Central South University, 21:439-446. DOI: 10.1007/s1177101419579 2. Chen, Y., He, J., King, M., Feng, Z. and Chang, W. (2013). Comparison of two suspension control strategies for multi-axle heavy truck. Journal of Central South University, 20(2): 550-562. 3. Chen, Y., He, J., King, M., Chen, W. and Zhang, W. (2013). Effect of driving conditions and suspension parameters on dynamic load-sharing of longitudinal-connected air suspensions. Science China Technological Sciences, 56(3): 666-676. DOI: 10.1007/s11431-012-5091-3 4. Chen, Y., He., J., King, M., Chen, W. and Zhang, W. (2013). Model development and dynamic load-sharing analysis of longitudinal-connected air suspensions. Strojniški Vestnik - Journal of Mechanical Engineering, 59(1):14-24. 5. Chen, Y., He, J., King, M., Liu, H. and Zhang, W. (2013). Dynamic load-sharing of longitudinal-connected air suspensions of a tri-axle semi-trailer. Proceedings of Transportation Research Board Annual Conference, Washington DC, 13-17 January 2013, paper no. 13-1117. 6. He, J., Qi, Z., Hang, W., King, M., and Zhao, C. (2011). Numerical evaluation of pollutant dispersion at a toll plaza based on system dynamics and Computational Fluid Dynamics models. Transportation Research Part C, 19(2011):510-520. 7. Zhang, C., He, J., Wang, Z., Yin, R. and King, M. (2013). Assessment of traffic noise level before and after freeway widening using traffic microsimulation and a refined classic noise prediction method. Proceedings of Transportation Research Board Annual Conference, Washington DC, 13-17 January 2013, paper no. 13-2016.

Static load bearing exercises during rehabilitation of individuals with transfemoral amputation fitted with osseointegrated implant: Kinetic analysis

The desire to solve problems caused by socket prostheses in transfemoral amputees and the acquired success of osseointegration in the dental application has led to the introduction of osseointegration in the orthopedic surgery. Since its first introduction in 1990 in Gothenburg Sweden the osseointegrated (OI) orthopedic fixation has proven several benefits[1]. The surgery consists of two surgical procedures followed by a lengthy rehabilitation program. The rehabilitation program after an OI implant includes a specific training period with a short training prosthesis. Since mechanical loading is considered to be one of the key factors that influence bone mass and the osseointegration of bone-anchored implants, the rehabilitation program will also need to include some form of load bearing exercises (LBE). To date there are two frequently used commercially available human implants. We can find proof in the literature that load bearing exercises are performed by patients with both types of OI implants. We refer to two articles, a first one written by Dr. Aschoff and all and published in 2010 in the Journal of Bone and Joint Surgery.[2] The second one presented by Hagberg et al in 2009 gives a very thorough description of the rehabilitation program of TFA fitted with an OPRA implant. The progression of the load however is determined individually according to the residual skeleton’s quality, pain level and body weight of the participant.[1] Patients are using a classical bathroom weighing scale to control the load on the implant during the course of their rehabilitation. The bathroom scale is an affordable and easy-to-use device but it has some important shortcomings. The scale provides instantaneous feedback to the patient only on the magnitude of the vertical component of the applied force. The forces and moments applied along and around the three axes of the implant are unknown. Although there are different ways to assess the load on the implant for instance through inverse dynamics in a motion analysis laboratory [3-6] this assessment is challenging. A recent proof- of-concept study by Frossard et al (2009) showed that the shortcomings of the weighing scale can be overcome by a portable kinetic system based on a commercial transducer[7].

Numerical study on the surface depression of the concrete runway pavement under impact load

Airport runway pavement always subjected to huge impact loading due to the hard landing of aircraft on the pavement surface. Therefore runway pavements should have sufficient impact resistance capability to avoid damage causing by hard impact like surface deflection in downward or penetration since the repair works is cumbersome within the operating condition of airport and also increases the service life cost of the pavement structure. Several research works have been carried out on airport runway pavement to measure the present condition of pavement and also to predict future performance of it. However, most of the works are confined by pavement response under moving aircraft loading. Nevertheless, no comprehensive research work is yet conducted to identify the controlling factors which might have significant effect in changing the common pavements damage like surface penetration depth under impact of aircraft. Therefore, a 3D FE study is conducted to determine some effective factors in controlling the top surface penetration depth of runway pavement. Among the exterior factors, mass of the impactor, velocity of the impactor, impact angle and boundary conditions are selected and as interior factors, thickness of the runway pavement, compressive strength and density of materials used in the runway pavement are selected.

Taxonomy of factors which influence heavy metal build-up on urban road surfaces

Heavy metals build-up on urban road surfaces is a complex process and influenced by a diverse range of factors. Although numerous research studies have been conducted in the area of heavy metals build-up, limited research has been undertaken to rank these factors in terms of their influence on the build-up process. This results in limitations in the identification of the most critical factor/s for accurately estimating heavy metal loads and for designing effective stormwater treatment measures. The research study undertook an in-depth analysis of the factors which influence heavy metals build-up based on data generated from a number of different geographical locations around the world. Traffic volume was found to be the highest ranked factor in terms of influencing heavy metals build-up while land use was ranked the second. Proximity to arterial roads, antecedent dry days and road surface roughness has a relatively lower ranking. Furthermore, the study outcomes advances the conceptual understanding of heavy metals build-up based on the finding that with increasing traffic volume, total heavy metal build-up load increases while the variability decreases. The outcomes from this research study are expected to contribute to more accurate estimation of heavy metals build-up loads leading to more effective stormwater treatment design.