Development of a practical aggregate spatial road freight modal demand model system for truck and commodity movements with an application of a distance-based charging regime

We promote a view that more attention should be given to the freight sector in order to recognise that many initiatives designed to impact on passenger travel do also impact on the performance of the movement of freight vehicles and hence the ability to distribute commodities from the locations in which they are generated. This paper develops a practical freight demand model system and applies the models within an integrated passenger and freight model system for the Greater Sydney Metropolitan Area using a distance based charge for trucks and cars as a way of highlighting the importance of not ignoring truck traffic when assessing road pricing reform in the car passenger sector.


Introduction
The distribution of freight is a crucial feature of modern society. Between 2007 and 2016, the Australian domestic freight task grew by 50%, with 726 billion tonne kilometres of freight moved in 2015-2016 representing an annual average of 30,000 tonne kilometres for every person in Australia (Bureau of Infrastructure Transport and Regional Economics 2018). Close to 34% of all freight tonne kilometres, bulk plus non-bulk, in Australia are carried by trucks with trucks dominating in the non-bulk sector. The freight task is forecast to grow by an additional 26% between 2016 and 2026, 1 adding to congestion in the road network, especially approaching and within metropolitan areas.
The current interest in freight distribution by road stems from a plan to find a practical way of incorporating the freight task into an integrated transport and land use model system, be it with a metropolitan or regional/rural focus, which historically, with rare exception, has focussed on passenger travel activity. We promote a view that attention should be given to the freight sector within strategic transport model systems in order to recognise that many initiatives designed to impact on passenger travel do also impact on the performance of the movement of freight vehicles, and hence the ability to distribute commodities from the locations from which they are generated. De  has recently reviewed the broader literature on this theme, reinforcing the integration objective of this paper. The focus of this paper is in showing how an analyst can develop a practical policy sensitive freight model system that can be integrated into a strategic transport model system together with passenger models using available aggregate data sources.
In studying freight distribution activity, we need to identify the volume of freight being moved by commodity class between the spatially defined origins and destinations of a study area, and the amount and share of truck movements by type of road vehicle (namely rigid and articulated trucks) associated with each commodity class. Allowance for empty trucks is also of great relevance. While this may seem like a straightforward specification of the essential elements of a practical freight movement demand model system, the identification of appropriate data in many countries is quite a challenge. In the Australian context, there are a number of candidate data sources, but no one source provides all the data required. For example, while at least two data sources provide volumes of commodities transported between an agreed spatial representation of origins and destinations, the commodity classes are not strictly comparable; and one data set does not split kilo-tonnes by rigid and articulated trucks, and only one data set separates demand volume by port and domestic movement, the latter important to account for the import and export setting of freight movement.
In this paper, we share the journey of data sourcing, data incompatibility and ways in which data can be merged in order to take advantage of the strengths of specific data sets. All this data preparation is guided by the models that have to be estimated in a format that makes them spatially and output-wise compatible with two integrated transport and land use model systems into which they will be positioned; a metropolitan focussed system called Metroscan, developed at The University of Sydney (Ho et al. 2017;Hensher et al. 2020a, b, also Appendix A in the current paper), and a national and Statewide model system called RTresis 2 that incorporates metropolitan, rural and regional contexts. Importantly, the metropolitan model must recognise that some freight flows begin or end in origins or destinations that are outside of the metropolitan area and hence external travel zones must be included to capture these features of freight movements. MetroScan and RTresis feed demand forecasts into a benefit-cost and economic analysis framework as shown schematically in Appendix A and more detail in .
An important interaction between truck movements and passenger activity is through the volumes of vehicle movements on the road network. Within the traffic assignment procedure in MetroScan, we have to load all passenger and freight truck movements onto the road network, and together with the service levels on each link, obtain equilibrium levels of traffic movement and the associated performance metrics such as travel times and delays associated with congestion. Accounting for freight movement,, 3 is of great relevance and failure to do so results in an inability in benefit-cost and economic impact analyses to capture the fuller productivity benefits from improved transport infrastructure investment and other policy initiatives, a major source of wider economic impact.
The US Quick Response Freight Methods (QRFM) approach (Beagan et al. 2019) emphasises freight usage and performance, and the need to integrate freight demand models into a broader integrated transport and land use planning and policy system like MetroScan or the Oregon State Transportation and Land Use Model Integration Program (TLUMIP) (Donnelly et al. 2018), or the number of systems developed in Europe (summarised in , all of which recognise that other vehicles share networks with freight such as light commercial and passenger vehicles. Freight truck travel times reflect the congestion times in many localities, resulting from those other demands. As Beagan et al. (2019, page 10) state "Unless planners include the usage and performance of these other components of travel, it may not be possible to forecast freight demand and performance." This paper and its contribution is organised as follows. We begin by setting out the model framework within which to capture the key components of freight and truck movement. 4 This is followed by a detailed presentation of the data items required and where they were sourced from, and what assumptions were required to merge the data sources in order to have a data set compliant with the modelling requirements. A descriptive profile of the data is then provided, followed by the estimated models and the interpretation of the main findings. To ensure that the models are able to be built into MetroScan, they have to be calibrated to the base year vehicle type shares and overall frequency of movements by commodity class. We explain how this is done, and the final models are presented ready for implementation. A case study is for the Greater Sydney Metropolitan area is set out, with the freight model system integrated into MetroScan with an application of a distance-based road user charge with separate rates per kilometre applied to trucks and passenger cars.

The Road Freight Movement Model Framework
Four models are proposed as the representation of the main behavioural features of truck and commodity movements between spatially defined locations ( Fig. 1). We need to know the total numbers of vehicles by class moving between each origin and destination, and the amount of each commodity within a defined set of classes that is carried by each vehicle type. Truck movements often do not carry freight 5 and hence we have to allow for empty truck movements. This will enable us to predict the amount of freight vehicle movement activity as well as volumes of commodities being transported under a reference setting and a scenario setting, where the latter enables changes in freight distribution due in particular to exogenous socio-demographic influences such as population and employment forecasted changes, and changes in the transport network (e.g., improvements in the road network). A diagrammatic overview of the linkages between these models is given in Fig. 1 and also in the Appendix A schematic diagram for the interface with the overall Metroscan model system (Hensher et al. 2020a, b) where the freight component is embedded with passenger and land use models.
Formally, the quantity, X, of total truck movements for each of type k commodities transported, including empty truck movements (and hence zero kilo-tonnes transported), from origin zone i to destination j by truck class m, X k ijm , can be expressed as Eq. (1).
X k ij is the quantity of overall truck movements associated with commodities of type k produced in origin zone i and transported to destination zone j (by the distribution model), as well as zero commodity movement in empty trucks, and U k ijm is the observed utility of transporting commodity k from origin zone i to destination zone j by truck class m.
(1) Fig. 1 The structure of the model system and linkages Aggregate truck type share model The representative utility associated with the aggregate truck type share specification in (1) can be decomposed into the additive observed and unobserved influences of relative utility as in Eq. (2).
V k ijm represents the observed influences such as the volume of freight moved (in kilotonnes), the incidence of empty trucks, and level of service variables such as travel time between each ij pair, as well socio-economic influences such as population size that are specific to location i or j; and k ijm is an index of the aggregate unobserved influences, assumed to be represented in the sampled population as an extreme value type 1 distribution which, under the independent and identically distributed (IID) assumption, results in the popular multinomial logit (MNL) model form. Importantly, to ensure full clarity of the form in which the MNL model is implemented, given the aggregate nature of the frequency data, we do not refer to a 'truck type choice' model, but to an 'aggregate truck share' model. In Nlogit6, the software used (Greene 2016, Chapter N-321), frequency data consists of a set of frequency counts for the outcomes. Such frequencies are non-negative integers for the outcomes in the choice set and may be zero.

Aggregate truck movement model
(1) represents the annual frequency of aggregate truck movements for each commodity class, whose functional form is now discussed. Movement frequency is a positive number compliant with a count model such as zero inflation Poisson (ZIP) 6 and heterogeneity. Using standard regression methods, as is common in truck generation models, is not strictly correct. Truck movements are strictly count variables and should not be considered as continuous variables. We do, however, understand that some authors assimilate them as such in large scale applications with large counts. Estimation of a Poisson model instead of a linear regression is a very straightforward approach. Not only does it avoid an irrelevant estimation method based on a somewhat weak law of large numbers when it is not an empirical case, it is a correct scientific statistical approach.
The frequency model is connected to the aggregate truck type share model by the expected maximum utility (also referred to as logsum and defined as ln* ∑ Hensher et al.2015 ) associated with the aggregate truck type share model.
The total frequency of truck movements is a non-negative count value. For a discrete random variable, Y, observed over a period of length T n and observed frequencies, y n , (n observations), and explanatory variables x n (e.g., logsum), the Poisson regression model shown as Eq.
(3) as given in Greene (2000). (2) The selection of Poisson or other models is judged by the data. Poisson is typically used with discrete count data but can be used to model a continuous variable if the variable is not negative and satisfies the Poisson requirement. The assumption an equal mean and variance is not so strict; just like regression can be used if the distribution is not exactly normal. If data is highly dispersed, then a negative binomial (NB) may be considered (Poisson is a special case of NB), but our model results predict well and we support the ZIP specification in modelling the kilo-tonne model. Adding an extra dispersion term of (log of) for distributed heterogeneity in the NB did not change the results.
1 3 λ n is both the mean and variance of y n ; E[y n |x n ] = λ n . We allow for unobserved heterogeneity to recognise the possibility of partial observability if data on observed trip activity at the commodity and OD level exhibits zero trips. Specifically, the answer 'zero' could arise from two underlying responses. If we were unable to capture any truck movements between a specific origin-destination (OD) pair within a commodity class, we would only observe a zero; however, the zero may be due to the measurement period and the response might be some positive number in other periods. We define z = 0 if the response would always be 0, 1 if a Poisson model applies; y = the response from the Poisson model; then zy = the observed response. The probabilities of the various outcomes in the ZIP model are shown in (4).
The ZIP model is given as (Greene 2012) Y n = 0 with probability q n ; Y i ~ Poisson (λ n ) with probability 1q n so that Prob[Y n = 0] = q n + [1q n ]R n (0) and Prob[Y n = r > 0] = [1 q n ]R n (r). R n (y) = the Poisson probability = e −λn λ n yn /y n ! and λ n = e β′xn . We assume that the ancillary, state probability, q n , is distributed logistically; q n ~ Logistic [v n ]. Let F[v n ] denote the logistic CDF. Then, ʋ i can be defined by the form in Eq. (5) labelled the ZIP(τ) model (Greene 2012). τ may be positive or negative and if there is evidence of zero trips in any observations, the τ parameter is expected to be statistically significant; otherwise the Poisson form with normal latent heterogeneity is adopted.

The kilo-tonnes model
The annual kilo-tonnes per OD pair and commodity class 7 model is designed to recognise that exogenous changes, such as the growth in population and employment, will change the volume of commodities by class that need to be distributed. A zero inflated Poisson regression model is estimated of the same form as presented above for truck movement frequency. Any change in kilo-tonnes due to an exogenous shock will influence the vehicle type shares, and through the logsum obtained from the aggregate truck type share model, will impact on aggregate truck movements between each OD pair for a commodity class. There remains the possibility that kilo-tonnes is an endogenous influence on the aggregate truck share, as in the presence of empty trucks, and so we have to test for this possibility. There are a number of ways to set up a discrete choice share model that embeds the presence of endogeneity associated with a specific inclusion in the representative component of a utility expression (see Train and Wilson (2009), Guevara et al. (2019), and Wooldridge (2015)).
Control functions are statistical methods to correct for endogeneity by modelling the endogeneity in the relevant random components (Wooldridge 2015;Hensher et al. 2021). This involves two stages: (i) the endogenous variable is regressed on exogenous instruments; then, (ii) the residual (or a function of it) is incorporated into the utility function of the aggregate truck type share model as an explanatory variable denoted the control. An advantage of the control function approach is that the test that the parameter on the control function is zero is equivalent to a test of exogeneity. For the Poisson regression model, the generalized residual is r qi = Y(i) − Exp(beta'x(i)) = Observed(i) − Expected(i).

Aggregate truck movement frequency
To predict changes in aggregate truck movement frequency, the logsum (or expected maximum utility) associated with each OD pair and commodity class as an output from the aggregate truck type share model can be used as a predictor variable in a frequency (truck movement) model for each commodity class, as a way of predicting the overall movement frequency.
The logsum from the aggregate truck type share model provides the mechanism for linking two models in a nested structure that aligns with an underlying theory of random utility maximisation (see Hensher et al. 2015). Specifically, our model system enables us to link the truck movement shares with the absolute frequency of all truck movements between each OD pair for a given commodity class. The movement frequency is a positive number which is compliant with a count model. We estimated the frequency count model with the logsum as a predictor.

Model estimation order
The order of model estimation is as follows (see Fig. 1). Firstly, the kilo-tonne models are estimated for each commodity class and truck type at the OD pair level. Given that we want to establish whether kilo-tonnes in the aggregate truck type share model is endogenous, the residuals from the Poisson regressions, which are now control functions, are saved and are included together with the kilo-tonne variable in the truck type model.
Secondly, a separate model is estimated to distinguish empty and commodity-laden aggregate truck movements at the OD pair level. If the residuals are statistically significant, this suggests that there are correlated errors between the random errors of the aggregate truck type share models and the errors of the kilo-tonne model, which are now purged through the inclusion of these residuals. 8 If the residuals are not statistically significant, then the exogeneity assumption is assumed to hold.
Thirdly, a multinomial logit model using frequency counts for each truck type is estimated. The dependent variable is aggregated truck movements per annum by type between each OD pair for each commodity class (including zero commodity on board) for each of rigid and articulated trucks. The previous two models are linked using the predicted kilotonnes and the probability of empty truck movements. Residuals are tested for the presence or otherwise of endogeneity associated with the prior two models.
Finally, the logsum obtained from the aggregate truck type share model becomes an explanatory variable in the aggregate truck movement Poisson regression model. We now move to discuss the data sources used to construct the data set used in model estimation.

Sourcing and constructing a road freight database for Australia
Developing the aggregated spatial freight modal demand model system requires that we have several key elements in the dataset. First, the dataset should cover freight movement regions within Australia with details provided for both origins and destinations such as population and other socio-demographics for characterising the location that may have an impact on freight movement. Second, the dataset should also apply a widely adopted list of commodity types with which freights movements, kilo tonnes of goods, kilometres travelled and other freight information can be classified. Third, the dataset should allow filtering information by alternative modes, in this instance, rigid and articulated trucks. With the dataset containing the above three types of data, it can be used to model the relationship between frequency of movement and amount of freight by commodity type, truck type, and other information to inform how commodities are moved from one location to another.
No single available data source can provide all the required information. We were able to source and combine three data sets to construct the required freight dataset: (1) the road freight data set from the Australian Bureau of Statistics (ABS, 2014); (2) the "Strategic Freight Movement" (SFM) data set from Transport Performance and Analytics (TPA, 2016) of Transport for New South Wales; and (3) the number of movements by vehicle for both rigid and articulated trucks from TPA (TPA, 2016). The key characteristics of the three data sets are summarised in Table 1 The ABS data is based on a sample survey of 16,000 articulated and rigid vehicles that were registered with an Australian motor vehicle registry during the collection period. Respondents were asked to provide information about their vehicle configuration, load carried and origin-destination for all trips undertaken within a randomly allocated reference week. See https:// www. abs. gov. au/ ausst ats/ abs@. nsf/ mf/ 9223.0. SFM and TPA used ABS data and obtained unit record data to use in obtaining estimates of the variables not provided by ABS: https:// opend ata. trans port. nsw. gov. au/ datas et/ freig ht-data https:// data. datah ub. freig htaus tralia. gov. au/ it/ datas et/ strat egic-freig ht-movem ent-model.

3
The ABS data set contains road freight data obtained from a survey conducted in 2014. The data include total kilometres travelled, total tonnes carried, and total kilometre tonnes for both rigid and articulated trucks for 23 commodity items, across 351 SA3 10 level origins and destinations (see maps in Appendix B) for an entire year. It does not contain truck movement frequency between each OD pair by either commodity type or truck type.
The SFM data set contains forecasted freight demand, at the SA3 level, for the 40-year period between 2016 and 2056 for road and other modes. We extracted 2016 data from this data sets to have a close match with the other two data sets. The data includes annual kilo tonnes for each OD pair for a total of 49 commodities. 11 This commodity list is different from the one used by ABS but can be reclassified to the ABS types. This dataset does not provide separate freight kilo tonnes for each truck type, providing only the total kilo tonnes for each commodity class. The TPA data set provides the movement frequency by vehicle type which was missing in the ABS and SFM databases. The remaining fields in the data set correspond to the SFM database matching by origin, destination and commodity type.

Descriptive profile of the combined data
A descriptive profile of the data ready for model estimation, as a merging of the three data sets, is summarised in Tables 2 and 3. It is the most complete data set currently available in Australia for estimating the aggregate spatial models presented in a previous section.
With the exception of coal, all commodities are transported by a mix of rigid and articulated trucks, with food, followed closely by cork and wood, and crude materials, contributing the most to the truck movement task, with the exception of empty trucks. There is a good overall share of truck movements between the two classes of vehicles for all commodities, except coal, with rigid trucks dominating, in particular, for other manufactured articles, chemicals, machinery and transport equipment. Excluding coal, general freight and animal and vegetable oils, fats and waxes are mainly moved by articulated vehicles. 24.74% of all truck movements do not carry any commodities, 12 indicating vividly the opportunity to use this spare capacity through improved scheduling and greater cooperation amongst trucking companies. Two categories, "natural and manufactured gases" and "manufactured fertilisers" did not record movements by rigid and articulated trucks.
The volume of commodities carried per vehicle in kilo-tonnes per annum per vehicle per OD pair is shown in Table 3, noting that the product of the Table 3 OD volumes and  Table 2 truck movements can be used to obtain the annual total volume of commodities 11 The kilo-tonne estimates are derived from the SFM and as such are forecasted estimates. We went to great pains to check the estimates based on a few available primary data on kilo-tonnes given that the primary input data was too sparse for all OD pairs. 12 The percentage of empty trips is often much higher, approaching 50%, in some urban areas in the US. We acknowledge this advice from Kouros (Abolfazl Mohammadian). 10 Statistical area level 3 (SA3) is one of the spatial units defined under the Australian Statistical Geography Standard (ASGS). This level of spatial details provides a very good strategic level representation of areas within metropolitan areas; for example the Greater Sydney Metropolitan Area has 80 zones, from the north in Newcastle to the south in Wollongong of which 46 are in the Sydney Area. The ASGS provides a common framework of statistical geography and thereby enables the production of statistics which are comparable and can be spatially integrated. There are 358 spatial SA3 regions covering the whole of Australia without gaps or overlaps (Appendix B). For details see.
https:// www. abs. gov. au/ ausst ats/ abs@. nsf/ Lookup/ by% 20Sub ject/ 1270.0. 55. 001~July% 20201 6~Main% 20Fea tures ~Stati stical% 20Area% 20Lev el% 203% 20(SA3) ~10015. transported in each commodity class and vehicle type. The standard deviation numbers relative to the mean are of particular interest since they show the variation in volumes carried per truck type, with the greatest variation (defined by the coefficient of variation (CV), which is the Standard Deviation/Mean*100) being for cork and wood for rigid vehicles and cereal and grains for articulated trucks. There is a significant amount of variability in the volumes carried by each vehicle class within each commodity class. The profiles in Table 2 and 3 describe the dependent variables used in the models to be estimated that we now turn to.

The zero-inflated poisson model for annual kilo-tonnes per vehicle per OD pair
The first model in the system describes annual kilo-tonnes per vehicle between each OD pair for rigid and articulated trucks using a Zero-Inflated Poisson (ZIP) model, given the distribution of the dependent variable. There is a highly skewed distribution of kilotonnes between each OD pair for each commodity class; being towards zero (empty truck movements 13 ) and covering rigid and articulated trucks combined, with a mean of 8.05 per vehicle, a standard deviation of 51.32 and 1, 1, 5, and 17 respectively for 25, 50, 75 and 90 percentiles. As the dependent variable, kilo-tonnes can neither be modelled using ordinary least squares regression or an unaltered Poisson model, which would demand either a normal or a Poisson distribution. The resulting ZIP models for both rigid and articulated trucks are shown in Table 4. The Vuong statistics of 26.13 and 24.87 suggest that the ZIP model is strongly favoured over an unaltered Poisson model for both truck types. The volume of kilo-tonnes per vehicle per annum transported between each OD pair is statistically influenced by the population at the destination, which is a good indicator of the demand for commodities, with a larger population demanding more kilo-tonnes of commodities. In addition, we have a statistically significant effect that is positive for the kilo-tonnes distributed to locations in New South Wales compared to the other States. We do not have a measure of jobs at each origin, but as a proxy we use the number of people in the population who work, and although they may not work in the same zone in which they live, we assume that the statistical significance of the variable is a proxy for the catchment area in which individuals work. We know that there is a very high correlation between where workers live and where the jobs are even if it is not the same individuals in both definitions.

Accounting for empty truck movements
As shown in Table 2, there are many truck trips where no commodities are being moved. Although empty truck movements are treated as zero volumes of freight between each OD pair inTable 4, 14 we need to develop an additional model to identify the incidence of the number of empty trucks as influenced by the volumes of commodities, as well as the production and attraction factors such as population as a good proxy from economic activity and location specific effects. The probability of empty truck trips can be modelled as a function of commodity flows of various commodities, and other influences, following the suggestions in the literature (e.g., Holguín-Veras & Thorson 2003;Holguín-Veras & Patil 2008). Table 5 presents the two aggregate logit share models, given the aggregate frequency of truck movements for each truck type, in predicting the probability of a truck being empty versus carrying commodities, between each OD pair, with one model for the rigid truck and the other model for the articulated truck. The model findings show that the flows of some commodity classes are statistically significant in predicting the probability of trucks being empty versus non-empty between each OD pair. For example, the positive parameter estimate for live animal flows suggests that the more movement of live animals occurs, the higher is the probability that some trips between specific OD pairs for rigid and articulated trucks will move empty (often associated with the backhaul). Conversely, where petroleum and petroleum products are being transported, a higher volume of distribution tends to reduce the probability of empty rigid and articulated trucks. These findings highlight the importance of accounting for the commodity movement profile at the OD level in conditioning the likelihood of having vehicles on the road that are not transporting any product.
We also find that there are some State specific impacts associated with the incidence of empty trucks, noting that the positive parameter estimates for five States (New South Wales, Queensland, Western Australia, South Australia, Victoria), compared to Tasmania and the Northern Territory, have systematic differential impacts on the likelihood of empty truck volumes. We know that there are a lot of directional impacts associated with the lack of backhaul business, and this is greater in the larger populated States, with evidence that there is a higher tendency in the lesser populated States to ensure there is a backhaul load before agreeing to undertake a distribution activity which often is long distance, as in the Northern Territory. But also, there are far more trucks available in the network in the larger States which adds to the risk of empty truck movements. The results also indicate that OD pairs associated with greater populations are less likely to be associated with empty truck movements. This is plausible since the amount of demand ensures that vehicles are better utilised, even if there is a high level of competing business activity available to move commodities. The predicted probability of empty truck movements, rigid and articulated, at an OD level for each commodity, is introduced into the truck type choice model as discussed below. This should be thought of as similar to a kilo tonnes trip except that they are zero kilo tonnes, and in that sense it is like a 'commodity class', and introduced into the aggregate truck type share model to remove the bias attributed to otherwise failing to account for movements that do not carry any commodities. Another way of saying this is that it is a class of movement which, as we have shown in Table 5, is influenced by exogenous factors, and if not accounted for would result in all predicted truck movements being associated with non-zero commodity flows.

Aggregate truck type share model
The next model in the system is the multinomial logit model describing aggregate truck type shares, with results summarised in Table 6. The dependent variable is aggregated truck movements per annum by type between each OD pair for each commodity class. 15 Frequencies are a transformation of proportions or shares, and hence are equivalent. The log likelihood function associated with the model is given as Eq. (6).
W i is a weight, and F ij is a frequency. For frequency data W i = 1 and F ij = the frequencies (or W i = F i1 + F i2 ) and F ij = Freq(i,j)/W i equal to a proportion. The model estimation process converts these frequencies to predicted choice probabilities or shares, which are calibrated by additional commodity-specific constants (see below) to reproduce the observed vehicle type shares at a commodity level.
The overall goodness-of-fit is in line with what most non-linear discrete choice models obtained (Hensher et al. 2015), which in our case is 0.173. The signs of the estimated coefficients match a priori expectations. The statistically significant variables, in addition to commodity class-specific dummy variables, include total annual kilo-tonnes carried by each truck type between each OD pair, the average travel time between the centroids of each origin and destination as a generic parameter, 16 a Sydney-specific origin dummy variable in the utility expression for rigid trucks, which is positive; and two destination-specific dummy variables for Western Australia and South Australia included in the utility expression for articulated trucks, which are negative influences of the probability of truck movements to these two States compared to other States.
We initially included the trip operating cost, but it was found to be highly correlated with travel time (given little congestion in many jurisdictions) and hence time and cost are both correlated with distance. We then used a generalised cost (or generalised time) specification converting cost or time to common units using a value of travel time savings. This also did not provide a meaningful result, being statistically insignificant, and so we have selected to use travel time only, which is essential when the model system is built into Metroscan with travel time the equilibrium assignment criterion. For subsequent use of the model system in Metroscan, where a scenario is a cost change (e.g., a carbon tax, congestion pricing, fuel changes), we developed a conversion between operating cost and travel 15 We are using very straightforward multinomial logit models, but with aggregated modal shares (or frequencies). The popular aggregate 'choice' model developed by Berry, Levinsohn and Pakes (BLP 1995) is not required here since users cannot influence the levels of the statistically significant explanatory variables, in aggregate. Several other papers have developed road vehicle type choice (or aggregate truck type share) models, such as Bok et al. (2020aBok et al. ( , 2020b and Abate et al. (2014). 16 We tested alternative specific parameters and they were not statistically different. time for each class of vehicle using the standard relationship between operating cost (in cents/km) and speed (in kilometres per hour). 17 See Appendix C. The control function correction is used to test for and correct if appropriate, endogeneity associated with kilo-tonnes (Table 4) and empty truck movements (Table 5) in the aggregate truck type share model. It is simply an auxiliary variable to correct for the endogeneity (Guevara et. al., 2015), so it is important to also include the exogenous variables that will be considered in the aggregate truck type share model, otherwise the control function correction may confound the effect of the exogenous variable with that of the residual. For the Poisson regression the generalized residual for kilo tonnes is r qi = Y(i)-Exp(beta'x(i)) = Observed(i)-Expected(i), and for the aggregate logit share model for empty vs non-empty truck movements, the residuals are Y (1,0) minus the predicted choice probability.
To test for endogeneity associated with kilo-tonnes, the residuals of both ZIP models for vehicle kilo-tonnes were also used as an explanatory variable for the aggregate truck type share model as a way of purging the model of endogeneity bias associated with the specific variables. The squared residuals were included in the aggregate truck type share model. The residual variables are statistically significant but with relatively small effect sizes; so while it is appropriate to include them to correct for endogeneity, the expected impact is likely to be minimal, suggesting that the exogeneity assumption may not be problematic. 18 The variable representing the probability of a truck being empty, obtained from the model reported in Table 5, is an estimate. The sign is positive, suggesting that, all other influences held constant, when the probability that a truck movement is empty compared to not being so, the probability has a differential impact on the share of truck movements that are rigid and articulated. The residual associated with the empty truck movement model was not statistically significant. Together with the estimated kilo-tonnes moved, this is an estimate, and hence to correct the estimated asymptotic covariance matrix for the randomness of the estimators carried forward, the standard Murphy and Topel (1985) correction is implemented, so that the standard errors are asymptotically efficient.
The dummy variables for commodity classes represent all classes; however, some commodities from the classes in Table 3 had statistically insignificant parameter estimates, due mainly to the relatively small amount of commodity being transported. The commodity class dummy variables were included in the utility expressions for either the rigid or articulated truck alternatives as a way of allowing for the role, on average, that the underlying characteristics of commodities in a class have on the preference for rigid or articulated trucks. They are effectively a decomposition of the alternative-specific constant to account for unobserved heterogeneity differences amongst the commodity classes that influence the predicted share of truck movements by vehicle type.
Elasticity estimates are obtained for two of the key variables of particular interest, travel time and annual kilo-tonnes. The direct elasticity of the aggregate share of truck movement by rigid trucks with respect to travel time between each OD pair across all commodity classes is -1.30. The equivalent estimate for articulated trucks is -0.62. These estimates are quite plausible, and we would expect a greater sensitivity for rigid vehicles because of the typical value of commodities. For kilo-tonnes carried between each OD pair across all commodity classes, the direct elasticities are -2.44 for rigid trucks and -0.62 for articulated vehicles. We have relatively elastic responses for rigid trucks and relatively inelastic responses for articulated vehicles. There is a greater sensitivity for rigid vehicles to changes in kilo-tonnes carried than the actual travel times, which suggests that volume is especially important to the truck movement business.
The next task is to calibrate the aggregate truck type share model to reproduce, at an OD level, the share of truck movements associated with each commodity class. The calibration procedure is straightforward (see Train 2009, p. 37) and initially involves adjusting the overall constant for each truck type (initially set at zero for articulated trucks), after accounting for the influence of all the variables, by a new parameter which is the original estimated alternative-specific constant plus, for each commodity class, the natural logarithm of the ratio of the observed share of truck movements in the reference year to the predicted share. The predicted shares are compared to the actual shares, 19 and the constants are adjusted up or down until there is a match of the actual and predicted shares. Once the aggregate shares are reproduced through calibration, and checked through the assignment mechanism, we then looked at each and every OD pair and added an additional calibration constant where the observed and predicted shares deviate by more than 5%. There were very few OD pairs that required additional calibration. The final calibrated constants at the aggregated level (i.e., for all OD pairs) are summarised in Table 7 with the OD specific additional constants available on request.
We compared the model predicted shares over the actual shares of the two types of trucks before and after the calibration of the constants. The average OD based prediction accuracy has been greatly improved, with the average probability prediction error reduced from − 0.091 to − 0.038 for articulated trucks and the average probability prediction error reduced from 0.117 to 0.041 for rigid truck by each OD pair.
It should be noted that the accuracy of predicted aggregate shares of rigid and articulated trucks are related to the available sample. For those commodities with very large sample size (e.g., food with 4651 observations), the accuracy of the predicted aggregate share is very high (i.e., the difference is as small as 1% to the actual share in the data). Some of the commodity types with a small sample of around 1000 also performed well (i.e., the accuracy is as close as 0.4% in the difference for a sample of 1294). The less accurate commodity types are often related to the limited sample size especially those below 1000 records in the dataset, such as iron and steel. We combined some of these commodities with very small sample sizes for the truck frequency models.

Aggregate truck frequency model
A Poisson model for aggregate truck movements was estimated for each commodity type (19 models) with a constant term, a parameter estimate for logsum and the sigma (σ) term, where the latter parameter is the standard deviation of heterogeneity. The Poisson model without zero-inflation for each commodity class produced both good model fits and a high level of prediction accuracy, as shown in Table 8 for each model, since the ZIP parameter 19 The assignment module in MetroScan utilises the open-source traffic assignment platform PLANit (https:// github. sydney. edu. au/ PLANit), developed at ITLS (University of Sydney). The assignment configuration is hard coded and conducts a traditional static traffic assignment with the following properties: Route choice and network loading is Deterministic User Equilibrium (DUE) following Beckman et al. (1956); the shortest path algorithm is Dijkstra one-to-all; smoothing uses the Method of Successive Averages (MSA); the number of iterations is user configurable; when set to 1 (default), DUE collapses to an All-Or-Nothing (AON) assignment; the duality gap threshold is user configurable; and the approach overall is based on Bovy and Jansen (1983).

Table 7
Calibrated constants for commodity types for the aggregate truck type share model  was not statistically significant. All parameters are statistically significant at the 1 percent level. The predicted frequency and actual frequency for each commodity type shows a high level of prediction accuracy with a regression R 2 ranging from 0.81 to 0.94 for all commodity types, including commodity types with a relatively small sample size (e.g., general freight, see Table 2). The predicted truck movement frequencies for each row can almost perfectly predict the original truck movement frequency (R 2 = 0.93 for actual and predicted frequencies). Including other candidate variables such as population, origin or destination dummies and interaction terms of logsums by commodity types, did not perform well. The prediction of the frequencies became different from the actual counts or negative and accuracy significantly dropped. We suggest this is because the tested variables are already significant explanatory influences on aggregate truck type shares and kilo-tonnes and hence are accommodated elsewhere in the overall model system. We also are aware of the risk of over-fitting and also of ensuring identification.
The "other commodity" category combines the categories with very small sample sizes including coal, animal and vegetable oils, iron and steel and the original 'other' commodity categories. Figure 2 shows the distribution of predicted logsums across the commodity classes. The expected maximum utility (or logsum) varies a lot, as might be expected, between each commodity class, suggesting the influences on the choice between rigid and articulated trucks varies a great deal. 20 Fig. 2 Logsums obtained by the aggregate truck type share model for each commodity class 20 For the application in Metroscan, we want to ensure that, for each OD pair and commodity-specific calculation, that the sign of the product of Logsum and its parameter estimate in the frequency model provides a directionally valid response to any change in one of more of the calculated estimates for Beta*Logsum. We found that the great majority of observations (i.e., OD pair by commodity class) result in the correct sign, and a negligible few OD pairs and commodity-specific observations had a different sign to that

Application of the Freight models within MetroScan
In this section, we present and discuss simulated results and predictions made for three related transport scenarios, using the MetroScan system for the Greater Sydney Metropolitan Area (GSMA) in which the freight models are implemented 21 (See Appendix A for the overall structure of MetroScan and Appendix B for the location of the GSMA). We illustrate the use of Metroscan with a distance based charge on trucks and cars, with a particular interest in identifying the importance of including freight activity with passenger activity in an integrated transport and land use model system. Appendix D provides a summary of other studies that have investigated distance based charges; although a word of caution is required in attempting the compare the evidence since contexts vary a great deal.
Scenario 1 assesses the impact in the transport sector and the wider economy of a distance based charge (DBC) introduced to passenger cars and trucks from the year 2021 onwards in the GSMA, at 5 cents per kilometre for cars and 20 cents per kilometre for trucks. Scenario 2 is related to Scenario 1, but only imposes a DBC on passenger cars at 5 cents per kilometre without a congestion charge for trucks. Scenario 3 is the opposite of Scenario 2, with a DBC only for trucks. These three scenarios enable us to identify the impact of a DBC when only passenger cars or trucks are assessed (i.e., we set the DBC to zero for the other mode), and when the DBC is imposed on both cars and trucks. We are unaware of this being considered in previous research. In the application through Metroscan, we obtain many outputs associated with efficiency, equity and environmental impact, including changes to passenger mode shares, truck movements, total kilometres and other measures.

Background for the three scenarios
The idea of road pricing reform in transport is not new but has been a challenge politically for many years (Hensher and Puckett 2008;Hensher and Mulley 2014a, b;Bliemer 2014, Bok et al. 2020a). While the debate has not gone away, there remains support from industry and academics for a revision to the way that cars and trucks are charged for the use of the road network, to reflect the true cost of usage, including the broader set of externalities such as emissions, safety and congestion. We propose a distance based charge for cars and trucks on all roads and at all times of the day in the GSMA. While variations could also be considered such as a peak period DBC and even a cordon-based charge in some locations, our interest herein is identifying the implications of a DBC on both freight and passenger vehicles in order to recognise that the costs and benefits of a DBC in one Footnote 20 (continued) expected, largely associated with very low frequencies. Note also that the zero-inflated feature associated with Tau accommodates the probability of Poisson being non-zero versus the probability of zero occurring. 21 We have not discussed the Time of Day (ToD) component of the freight module in MetroScan. Currently, unlike the passenger modules, we do not have a ToD model and instead 'allocate' truck trips based on exogenously observed truck movements from traffic counts from an OD matrix in the GSMA. sector, passenger or freight, is impacted by policies implemented in the other sector; simply put because they share the same road network, and hence are affected by behavioural responses in both sectors (Moridpour et al. 2015).
Using Metroscan and introducing a DBC of 20 cents for trucks and 5 cents for passenger vehicles, we are able to predict the potential impact on traffic volumes through changes in freight and passenger trips, as well as changes in mode shares in both sectors. An important objective of road pricing reform is to improve the efficiency of the road network through reductions in traffic congestion as well as to achieve broader objectives linked to user benefits (firms and individuals) and emissions. An extensive number of outputs are obtained and reported as a way to identify the many ways in which a DBC impacts on users and the broader economy. Appendix C discusses how a DBC is introduced into the freight model system where there is no travel cost attribute.
The introduction of a distance based charge of road freight and passenger vehicles has varying impacts throughout the model system as is shown in the various output results in Table 9 (Scenarios 1-3). These changes occur throughout the interrelated model system in MetroScan; however, for example, when we introduce the DBC on trucks only, the passenger models greatest impacted are the vehicle kilometres model, linked to changes in travel times that are fed up from the mobility investment module (essentially mode and time of day for work and non-work trips). There are noticeable switches between car and public transport (see Scenario 3). The change in the freight model system is linked back to the change in travel time (and the conversion formula for generalised cost) which impacts on the aggregate truck type shares and the linked (via a logsum) aggregate truck movements. In addition, we have the traffic assignment that adjusts travel times etc., until we obtain an equilibrium outcome in respect of levels of freight movement for each commodity type and passenger vehicle movements for each OD pair. Table 9 summarises, for the three scenarios, a number of informative output results for 2030, noting that Metroscan generates forecasts of all such outputs for each year from 2021 onwards up to 2056. The DBC is set to start from 2021 for the three scenarios, with 2030 being an appealing year in which to show ongoing and consistent impact of the DBC policy.

Impact on the transport system
The total distance travelled of rigid trucks (TVKMR), and articulated trucks (TVKMA) under a DBC has the largest percentage reduction at -6.20% and -3.09% respectively for Scenario 3 when the DBC is applied for trucks only. This results in a 2.29% increase in car kilometres (TVKM), clearly due to improved travel conditions where there are less truck kilometres in the network. In contrast, for a DBC associated with cars only (Scenario 2), we see a 13.74% reduction in total annual car kilometres but a 1.59% increase in rigid truck kilometres and a 1.33% increase for articulated truck kilometres. When we assess the impact of a DBC applied to all road modes (Scenario 1), we see a 12.17% decrease in car kilometres, a 4.54% reduction in rigid truck kilometres, and a 2.8% decrease for articulated trucks. All of these findings reinforce the importance of not focussing on one modal segment since the flow-on effects across all road activity are significant and will have implications for the economic appraisal and impact of the road pricing reform agenda. Even if we only focussed on a DBC for cars, the impact on truck movements needs to be recognised (at 1.59% under Scenario 2). The generalised cost per person trip per for car (GCpersT) and generalised cost per person trip for public transport (GCPubT) are given as follows: GCpersT = VTTS*in-vehicle time + VoR*buffer time + operating cost ($/trip) + tollcost ($/trip) for all purpose of trips (peak/offpeak); and GCPubT = invt VTTS *invehicle travel time + out-of-vehicle VTTS *out of vehicle travel time + PT fare ($/trip) for all purpose of trips(peak/offpeak) All trip purposes = commuting, non-commuting and business Base year $2019 A related impact is on travel time, defined as the total end-use travel time in present value terms (TEUCPV_TTC). For a DBC on only cars (Scenario 2), there is a 4.38% improvement in travel time, which when combined with a DBC for trucks (Scenario 1), improves further to 5.23%, given a gain of 3.42% for trucks. This is an expected result given the dominance of cars in the road network. If only trucks were subject to a DBC (Scenario 3), the travel time benefit for cars improves by only 0.39%; however, trucks obtain a 4.19% increase in benefit.
While the impact of a DBC delivers improved travel times and contributes to reducing traffic congestion, it also has desirable environmental benefits. We see an 11.29% reduction in CO 2 emissions when the DBC is imposed on both cars and trucks, with 12.17% reduction associated with cars and 3.47% reduction for trucks. These are sizeable reductions in enhanced greenhouse gas emissions. With a single sector DBC, we would see CO 2 increasing in the sector not subject to a DBC, with car increasing by 2.29% when a DBC is introduced for trucks only and 1.43% for trucks when the DBC is applied to cars only. The overall impact on total energy consumed (in litres of fuel) is also significant and aligned with CO 2 changes given their common link to fuel efficiency, with a 12.99% reduction for cars under a DBC applied only to cars, which reduces to 11.64% when a DBC is also imposed on trucks.
The implications on truck movements are of special interest given the model framework presented above. When a DBC is applied to cars and trucks, we see a 4.54% and a 2.8% reduction respectively in rigid (movefreqR) and articulated (movefreqA) truck movements, which is adjusted up to 6.2% and 3.09% respectively when the DBC is only imposed on trucks. This again shows the interdependencies between cars and trucks regardless of whether one or both are subject to the DBC policy instrument.
The implications of this evidence flow through to a number of other output measures. For example, when the DBC is imposed on both sectors, we see a decrease in government fuel excise of 12.17% from cars and 3.47% from trucks; but this is compensated by a significant increase in revenue from the DBC, being $11.972 m per day (or approximately $3.424 billion per annum) for cars and 1.197 billion annual revenue from trucks. This a sizeable revenue stream, close to $4.6 billion per annum. Another informative output indicator is the consumer surplus change consequent on DBC. Holding residential and firm location fixed, and allowing for modal and time of day switching responses for all passengers, we see an average loss of consumer surplus benefit of 15.92% per person trip for passenger trips by all modes under a DBC in both sectors, but an increase on average of 21.25% when the DBC is applied only to trucks, and a 33.74% reduction when only applied to cars. It should be noted that these percentage changes are associated with numerically small absolute levels of consumer surplus before and after the introduction of a DBC.
In the passenger sector, when we have a DBC only on car kilometres, the mode shares for drive alone (TDA) and cars with passengers (TRS) reduce by 0.63% and 8.41% respectively, while the shares for train and bus increase by 28.21 and 21.33 respectively for train and bus, noting that public transport has a relatively small share (7.74%) of the overall travel movements. When we impose a DBC on cars and trucks, the reduction in percentage changes is 0.28, 7.87 respectively for car drive alone and car with passengers, which is lower than for a car only DBC, which is expected given the improvement associated with fewer truck kilometres. This also results in a drop in the percentage share for train and bus, now 23.18 and 17.02 percent respectively. This translates into a healthy increase in public transport revenue of 19.67% under a DBC for both sectors and 24.27% when the DBC is applied to cars only (where there is a greater switch into public transport). This is another important finding indicating that the switch into public transport is tempered when we account for pricing reforms in both the passenger and freight sectors.
There are significantly different contributions to local air pollution for passenger cars and trucks, especially heavy articulated trucks. While the suggested cost per kilometre recommended by TfNSW (2020) for passenger cars is only 3.37 cents/km, the suggested air pollution costs are as high as 16.5 cents/km for rigid trucks and 65.82 cents/km for articulated trucks. The impact of changes in total kilometres for cars and trucks can result in different patterns for CO 2 emissions and local air pollution. The combined change in CO 2 from cars and trucks under the various DBC scenarios is closer to the level of change in cars attributed to the smaller difference in CO 2 emissions from cars and trucks and much larger total kilometres by cars (i.e., 113 g/km for cars and 160 g/km for trucks on average for CO 2 emission). In contrast, for air pollution, the air pollution from trucks accounts for more than 50% of air pollution in the land transport sector, and hence the relatively small change in truck kilometres can dominate a much more sizeable change in car kilometres. For example, when DBC is imposed on trucks only (Scenario 3), the reduction in air pollution for trucks is disproportionately greater than the increase in air pollution for cars, which is not the case for CO 2 .

Conclusions
This paper has set out to develop a very practical and easy way to apply a model system for freight demand truck movements and to place the model system within a setting in which freight movement activity competes side by side with passenger transportation in the road network. The appeal is that it is simple to use, and captures the key linkages between the demand for commodities, the role of truck types in moving this volume, including zero commodity volume, and what this means for the overall amount of truck movement on the road network. Commodity demand and land use effects that change over time (given exogenous forecasts of influences such as population and employment) can be used to inform the freight movement task, and how it impacts on the overall performance of the transport system, beyond only freight activity.
The final set of calibrated freight demand models are inter-dependent, which is necessary for an integration into a unified, integrated transport and land use model system such as MetroScan that also recognises the interdependencies, hence endogeneity and feedback, between passenger and freight movement activity. Specifically, the kilo-tonnes model is a structural representation of the generation of volumes of commodities, including zero volume that have to be distributed and is influenced by population and employment in particular, which are proxies for economic activity. For a given quantum of kilo-tonnes which enters the aggregate truck type share model, we can predict the share of truck movements between each OD pair within a commodity class and for circumstances when truck movements do not carry any freight. The predicted aggregate truck type shares have an associated logsum variable that represents the expected maximum utility associated with the underlying preferences of the truck share model. The logsum is the key influence on the aggregate truck movements between each OD pair within each commodity class.
This aggregate truck movement in the base and application scenario feeds into the road traffic assignment module of MetroScan, and together with other road traffic, mainly passenger cars and light commercial vehicles, is used in establishing the equilibrium performance level in the road network.
Through iterations and feedback into the freight model system, where adjustments in travel times influence the changing truck shares and frequencies in the network, we arrive at a solution given the convergence criteria in the traffic assignment algorithm. The solution represents the travel times and the volume of truck movements between each OD pair within commodity class on the roads linking these geographical jurisdictions. We have demonstrated the rich policy capability of the Metroscan framework that embeds the freight model system with passenger and locational modules in order to study the implications of introducing a distance-based charging regime for trucks and cars. The key policy take away is that a pricing strategy be it on cars alone, trucks alone or both impacts of the performance and cost of each mode in the road network, and that both should be allowed for in an studying the impact of road priding reform.
In ongoing research we are enhancing the model system in two ways. Firstly, we will develop appropriate ways to recognise the changes in freight and passenger movements associated with the growth in e-(or online) shopping. At the time of model development and application, e-shopping in Australia was very small (less than 3%-see https:// www. finder. com. au/ online-shopp ing-stati stics). These purchases will still have to be transported to warehouses or integrated logistics hubs or other locations as part of the supply chain and this will retain a large amount of freight distribution by trucks; however some freight will then be distributed by light commercial vehicles (LCVs) as part of the last mile. We need also to recognise and account for the growing collection of online orders from shops as part of click and collect, and hence the change in passenger trips will not be fully associated with online shopping. Currently, the growth of online shopping can be accommodated with the current version of MetroScan through scenario analysis by adjusting the amount of freight being delivered by light commercial vehicles (LCVs) as well as reductions in passenger trips. Secondly, we will investigate whether an additional model that considers the role of rail over road distribution for non-bulk commodities is warranted, given the dominance of non-bulk movements being by truck; however the initial focus on road freight distribution is important not only as the main way in which freight movement interacts with passenger movements in the land transport context, but trucks provide an important role in distributing the goods moved by coastal shipping, rail and air to and from the initial origins and final destinations, commonly referred to as the first and last mile.

Appendix C: Travel time and costs for rigid and articulated trucks
In the aggregate truck type share model, we have included the trip time as one of the variables to predict the aggregate truck type shares ("BTIME" in Table 6). In practice, we often need to examine the impact of transport cost instead of using trip time on truck movement. These changes may be related to changes in fuel, toll, excise rate or any other types of operation and maintenance costs. The impact for potential new policies such as extra charges due to emission control or congestion control are also highly likely for the future. Travel time and travel cost are closely related in many measures such as the Value of Time (VOT). Using widely accepted values for VOT such as the ones used in Benefit Cost Analysis (BCA) recommended by government departments, we have derived the coefficient for cost per hour for both truck types, using Eq. 7.
We have used the commonly accepted VOT in Australia of $42.63/hour for rigid trucks and $64.72/hour for articulated trucks (TfNSW 2020) after also checking sources such as the Australian Transport Assessment and Planning (ATAP) and our past research. The parameters for the rigid and articulated trucks are -0.0296 and -0.0195 relatively, for transport costs in the dollar per hour unit.
To test a scenario with cost change, we can alter the part of a utility expression contributed by travel time for aggregate truck type share from β time *Time to β cost *VOT*(New/Current cost)*Time without changing the utility. To be specific, for rigid and articulated trucks, the related utility component related to time and cost can be given in two forms as below.
Using time only: Using cost and time: When we implement the freight model in the MetroScan system, the car and truck interactions, such as travel time, were also linked, based on past research. For example, the travel time associated with trucks was linked with the travel time for cars using the following equation in Kim and Mahmassani (1978, p 37): In their model, β 0 was -4.78 and β 1 was 1.075, respectively. With these links built in, freight models could work cohesively with other models such as a passenger mode choice model within an integrated system to reflect changes in both passenger car and freight vehicle movements.

Appendix D: Summary of a number of studies that have investigated distance-based charging
The research on either DBC or time and distance combined charges have been discussed by researchers from different countries such as Europe, US, Singapore and Australia. DBC policies have been implemented in some EU countries such as the Netherlands, with a main focus of research on optimising the DBC design or quantum (e.g., . In countries like Singapore where the transport authority has implemented cordon-based charges for the Central area, research has investigated combining distance-based and timebased charges into a simple cordon-based charge (e.g., Gu et. al., 2018;Liu et. al., 2014 andMeng et. al., 2012). In the USA and Australia where governments are facing inadequate funds for maintaining existing transport and road systems, studies have looked into ways of establishing an equitable and progressive DBC schemes to achieve multiple goals including raising adequate funds for road investment and maintenance while improving equity in these charges for both passengers and businesses (e.g., Yang et. al., 2014).
A summary of a number of studies that have investigated distance-based charging. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.