Abstract
The Northwestern University Freight Rail Infrastructure & Energy Network Decarbonization (NUFRIEND) framework is a comprehensive industry-oriented tool for simulating the deployment of new energy technologies including biofuels, e-fuels, battery-electric, and hydrogen locomotives. By classifying fuel types into two categories based on deployment requirements, the associated optimal charging/fueling facility location and sizing problem are solved with a five-step framework. Life-cycle analysis (LCA) and techno-economic analysis (TEA) are used to estimate carbon reduction, capital investments, cost of carbon reduction, and operational impacts, enabling sensitivity analysis with operational and technological parameters. The framework is illustrated on lower-carbon drop-in fuels as well as battery-electric technology deployments for the US Eastern and Western Class I railroad networks. Drop-in fuel deployments are modeled as admixtures with diesel in existing locomotives, while battery-electric deployments are shown for varying technology penetration levels and locomotive ranges. When mixed in a 50% ratio with diesel, results show biodiesel’s capacity to reduce emissions at 36% with a cost of $ 0.13 per kilogram of CO2 reduced, while e-fuels offer a potential reduction of 50% of emissions at a cost of $ 0.22 per kilogram of CO2 reduced. Battery-electric results for 50% deployment over all ton-miles highlight the value of future innovations in battery energy densities as scenarios assuming 800-mi range locomotives show an estimated emissions reduction of 46% with a cost of $ 0.06 per kilogram of CO2 reduced, compared with 16% emissions reduction at a cost of $ 0.11 per kilogram of CO2 reduced for 400-mi range locomotives. The NUFRIEND framework provides a systematic method for comparing different alternative energy technologies and identifying potential challenges and benefits in their future deployments.
Keywords
The transportation sector is the largest contributor to greenhouse gas (GHG) emissions in the US, contributing 27% of the emissions in 2020 (
In 2019, the US freight rail sector accounted for approximately 40% of the national freight ton-miles and emitted nearly 40 megatons of CO2 into the atmosphere in the process, an amount equivalent to the emissions of all the passenger vehicles in Texas alone (
Diesel-electric locomotives have dominated US freight rail operations since the 1960s (
Several high-profile pilot studies have been conducted in partnership between multiple railroads, locomotive manufacturers, and local and state governments to test the viability of alternative technologies on revenue service (
Picking the right mix and schedules to invest and deploy the next generation of energy technologies is a challenging process. Technological uncertainties, network effects, regional economics, and economies of scale all render mathematical optimization formulations of the problem essentially intractable. Decarbonization decisions will no doubt have far-reaching environmental, operational, and financial impacts on railroads, shippers, regulators, and other stakeholders in the greater supply chain. While previous research focused on conventional fuel types and highly simplified railroad networks, there is a significant research gap in developing optimization models to support the deployment of infrastructure to support rail decarbonization.
Research Aims and Contributions
The Northwestern University Freight Rail Infrastructure & Energy Network Decarbonization (NUFRIEND) Framework was developed to assist the rail industry in planning and evaluating the adoption of alternative fuels for decarbonization efforts. Scenario-specific simulation and optimization modules provide estimates for emissions reduction, capital investments, cost of carbon reduction, and operational impacts for any deployment profile.
The framework relies on two different approaches to capture the characteristics and requirements of the two main groups of energy technologies:
Drop-in fuels (Figure 1): Lower-carbon drop-in fuels can directly replace diesel fuel in locomotives and refueling stations. We assume no significant changes to existing assets or infrastructure are required for their deployment. Figure 1 shows the flowchart of the framework developed to analyze drop-in fuel deployment scenarios.
Energy storage fuel technology (Figure 2): Hydrogen and battery-electric technology deployment poses a more complex problem as it requires significant investments to be made in the siting of refueling/charging facilities and the replacement of locomotive fleets. Thus, strategies for locating and sizing refueling/charging facilities on a railroad’s network to meet their energy demands must be developed to aid deployment decisions. However, jointly locating and sizing facilities quickly becomes a combinatorial problem because of the interconnectivity of the many potential facility locations seen on networks, and both regional economics and economies of scale affect the cost of a facility deployment strategy. To reduce the problem complexity, we decompose the facility location and facility sizing problems from each other and insert a flow routing module in between to assign freight flow that must be served by the selected facilities. Figure 2 depicts a flowchart highlighting the five-step framework developed to address the deployment of the refueling/charging infrastructure to support hydrogen or battery-electric locomotives.

Flowchart of framework to support the deployment of lower-carbon drop-in fuels on the rail network.

Flowchart of five-step framework to support the deployment of refueling/charging infrastructure for hydrogen and battery-electric technologies.
The NUFRIEND framework aids stakeholders in analyzing alternative fuel technology deployments for freight rail operations. We assess and classify potential energy technologies based on deployment requirements and provide practical alternatives to diesel locomotives. The otherwise intractable facility location and sizing problem is solved with a five-step framework consisting of nominal problems from graph theory. A key advantage is the flexibility to apply the framework for any railroad considering the specific network structure and freight demand, outputting evaluation metrics for the associated emissions and costs relative to diesel operations. Equipped with the capability to efficiently simulate technology adoption scenarios with life-cycle analysis (LCA) and techno-economic analysis (TEA), this framework supports sensitivity analyses for different operational and technological parameters through a transparent and flexible input module, thereby addressing the uncertainties surrounding technological developments.
This paper reviews the relevant literature in freight rail decarbonization and new energy technologies in the next section. The methodology of the five-step NUFRIEND framework, in particular LCA and TEA, is discussed in the third section, before introducing the data sources and parameters adopted in the fourth section. The results are then illustrated with examples and scenario comparisons in the US freight rail network in the fifth section, concluding with a discussion of limitations and future research directions in the final section.
Literature Review
Network design and network flow optimization problems have seen rich applications in rail, with a review of recent work on their applications in multimodal freight in Archetti et al. (
Early approaches to conducting freight rail decarbonization analyses date back to the 1970s with Schwarm (
The analysis of decarbonization pathways such as hydrogen or battery-electric technologies requires the deployment of refueling/charging facilities over the rail network. Taxonomies of facility location problems appear in Mirchandani and Francis (
Previous research efforts evaluate fuel efficiency and resulting emissions mainly through system-average efficiency factors (
As diesel-electric locomotives combust diesel to generate electricity onboard to power an electric motor, battery tender cars can be attached to power the electric motor directly with limited retrofitting (
Methodology
This framework simulates and evaluates the deployment of alternative fuel technologies on the rail network, and as such is dependent on the technological requirements of each energy technology. The energy technologies considered in this paper can be divided into two categories: (i) lower-carbon drop-in fuels such as biodiesels and e-fuels and (ii) new energy storage technologies such as battery-electric or hydrogen, as discussed in the introductory section.
Lower-carbon drop-in fuels are generally deployed as admixtures (e.g., 20% biodiesel and 80% diesel). Their deployment is considered uniform across the network at a desired penetration rate taken as the rate of the admixture. Thus, origin–destination (O-D) flows are simply routed by commodity group on the existing (baseline) network’s links. These flows are used to calculate the costs and emissions associated with their deployment by weighing the relevant fuel parameters with their corresponding admixture rates. Figure 1 demonstrates the essential steps and flow of information.
However, for alternative energy sources that require locomotive powertrains to be converted and new refueling infrastructure to be deployed, we present a five-step sequential framework to (a) select O-D paths that leverage economies of density, (b) locate facilities along these paths, (c) reroute flows on the existing and alternative technology rail networks, (d) size the facilities to serve the rerouted flows, and (e) evaluate the deployment in relation to their emissions, costs, and operational impacts. The flow of information between each of the steps can be seen in Figure 2.
Sequential Framework for Facility Location and Sizing
The joint facility location and sizing problem is a combinatorial problem to solve (
We formulate the facility location integer program as an extension of the set covering problem (
As flows are assumed to originate and terminate at any candidate facility location on the rail network, they may be assigned to the alternative technology enabled subnetwork. The actual assignment of flows depends on the specified routing policy, which may or may not allow for flows to be rerouted from their original routing on the baseline network, may have a maximum distance increase if flows are rerouted, or may only allow flows originating
As the energy intensities of moving different commodities are known to vary considerably (
LCA of Energy Technologies
In this study, we examine the GHG emissions of different energy technologies with a system boundary covering both the well-to-pump (WTP) and pump-to-wheel (PTW) stages. The functional unit for the emissions was set as gCO2/ton-mile. Together, WTP and PTW stages comprise well-to-wheel (WTW) analysis. While the WTP examines the environmental impact of production, transportation, and distribution of feedstock and fuels, PTW focuses on the vehicle operation. Note that the impacts from the vehicle manufacturing cycle, including stages such as material extraction, component manufacturing, assembly, and recycling of vehicle components, are beyond the scope of this analysis.
We use ANL’s GREET model (
For battery-electric locomotives, we consider the GHG emissions associated with the upstream processes of electricity generation and the losses in the transmission and distribution system. For current and projected (Mid-case Standard Scenario) generation mixes, we capture the state-wise variation of electricity generation mixes in gCO2/kWh supplied to the charger station based on the results from Burnham et al. (
TEA of Energy Technologies
For conventional diesel, biodiesel, and e-fuels, this study focuses on the levelized cost of refueling only, as these energy technologies do not require additional infrastructure investments. We estimate the levelized cost of refueling by multiplying the commodity-specific link flow, energy intensity, and fuel cost. For battery-electric and hydrogen technology, we consider the charging/refueling infrastructure cost in addition to the battery/hydrogen tender car capital investment and refueling costs. The sequential framework described above estimates the facility size based on the peak demand (in kWh/day for battery and kgH2/day for hydrogen) at each location. Based on a given locomotive battery/hydrogen tender car capacity, the peak and average demand is used to estimate the peak and average locomotive throughput at each facility. From the peak locomotive throughput, we estimate the number of chargers/pumps needed at each location to support the peak demand, such that a provided maximum station utilization is not exceeded. ANL’s bottom-up TEA tools, Heavy-duty Electric Vehicle Infrastructure Scenario Analysis Model (HEVISAM) for battery-electric, and Hydrogen Delivery Scenario Analysis Model (HDSAM) for hydrogen (
The levelized costs of operation are estimated in cost per quantity of energy (e.g., $/kWh) or fuel (e.g., $/kgH2, $/gallon) and are converted to cost per revenue ton-mile using energy intensity parameters from Table 1 and the results from the flow routing and facility sizing steps to determine the ton-miles served by a particular energy technology. The levelized cost per ton-mile provides a standardized way to represent alternative technology costs as an operational metric.
Parameters
The WTW GHG emissions and the levelized cost of operation could be synthesized into one metric to compare across energy technologies. The cost of avoided emissions of a particular technology is the ratio of the levelized cost of operations (in $/ton-mile) and the WTW GHG emissions intensity (in kgCO2/ton-mile), relative to the baseline diesel operations, as represented in Equation 4. The cost of avoided emissions serves as a key policy metric, as it provides the cost per unit of carbon reduced from emissions for a specific technology, a metric that can be compared with the social cost of carbon or carbon credit/tax schemes.
Application Data Sources and Parameters
The framework presented above is next illustrated through application to evaluate alternative decarbonization scenarios for the US Class I railroads network. In this application, the parameters involved in modeling rail operations and energy sources collected from multiple sources are summarized in Table 1.
Rail Network
The existing rail network was extracted from the North American Rail Network data set compiled by the Federal Railroad Administration and the Bureau of Transportation Statistics using work on rail facility classification done by Oak Ridge National Laboratory in WebTRAGIS (
Freight rail demand for 2019 was estimated from the STB’s annually compiled Carload Waybill Sample (CWS) (
Energy Requirement
As the flow routing from the CWS provides estimates of the ton-miles of goods by commodity on the network, these must be converted into terms of energy from which fuel or electricity consumption can be calculated. The energy required to move a ton of goods 1 mi varies by commodity because of differences in physical characteristics (e.g., density) and operational practices (e.g., shipment speed), by railroad because of topographical and fleet variations, and by energy technology because of differences in powertrain designs (
New Energy Sources
New energy sources are evaluated relative to diesel operations based primarily on differences in their cost and emissions. Baseline diesel cost and emissions data are extracted from the R-1 report (
Levelized Cost of Battery
For the case of battery-electric locomotive deployment, the sizeable capital cost of the batteries must be captured. Building on assumptions and data from Popovich et al. (
Levelized Cost of Hydrogen Tender Car
Having considerably lower energy density relative to diesel, hydrogen fuel locomotives will require the use of a tender car for fuel storage. The storage of hydrogen fuel requires advanced temperature and pressure control systems, making hydrogen tender cars capital-intensive investments (
Analysis of Results
To illustrate the functionality of the model, scenario simulation and evaluation results are shown for Western (BNSF, Canadian National, and Canadian Pacific Railways) and Eastern (CSX Transportation, Kansas City Southern, and Norfolk Southern Railways) railroad networks.
Figure 3 shows an example of the NUFRIEND dashboard. It allows users to model different scenarios based on inputs including railroads, energy sources, commodity groups, battery ranges, and target deployment percentages on the left, as well as other specific parameters on a separate pane. In the context of battery-electric and hydrogen deployment, the five-step sequential framework is applied to consider the railroad network and freight traffic, locate and size the charging facilities, route the rail traffic, and estimate the emissions and costs based on LCA and TEA, as outputted in the right. Metrics including WTW emissions, levelized costs of operation (LCO), the proportion of ton-miles served by each energy technology, the cost of avoided emissions, and other detailed operational results are shown. The scenario WTW emissions are the sum of emissions of diesel (blue) and battery/hydrogen (green) routes. The LCO of diesel is the current fuel cost, while the battery-electric LCO is composed of charging facility capital costs, battery capital, and operations and maintenance (O&M) cost. Electricity cost and the hydrogen LCO is composed of refueling facility capital costs, energy tender car capital cost, and fuel cost. Above the battery-electric/hydrogen LCO, the scenario average LCO is shown, which includes the cost of diesel refueling needed to serve the segments of the network not covered by battery-electric or hydrogen. Results are also shown down to the facility and track levels with green color denoting coverage of battery-electric or hydrogen technology. Users can hover over them to examine more granular information such as traffic volume, number of chargers or pumps, and station utilization.

Example of Northwestern University Freight Rail Infrastructure & Energy Network Decarbonization (NUFRIEND) dashboard.
Example Illustration
This section showcases the functionalities of the dashboard with scenario results for each of the two energy technology categories. Each of the figures shows aggregate emissions, cost, and deployment plots, as well as a network map that shows the specific alternative technology deployment. The WTW emissions bar plots show the emissions (in kton CO2) attributable to diesel (blue) and the alternative technology (green) operations, with the emissions intensity (in g CO2/ton-mile) of diesel and the deployment scenario in question overlaid (yellow diamond). The levelized cost of operation bar plots show the LCO (in cent/ton-mile) for diesel and the alternative technology in question. The key cost components are displayed, such as fuel (blue) for diesel, biofuel, e-fuel, and hydrogen; electricity (light blue) for battery; battery and hydrogen fuel tender car costs (orange); and charging/refueling station capital costs (red). The pie chart shows the deployment of the alternative technology as the share of ton-miles captured by diesel (blue) and by the alternative technology (green).
Battery-Electric and Hydrogen
Figures 4 and 5 show an example (hypothetical) deployment of a 400-mi range battery-electric locomotive technology for Western and Eastern railroads, respectively, needed to serve approximately 50% of their ton-mileage. Western rail networks are in general more expansive and require more charging facilities compared with Eastern railroads (57 versus 21). This also leads to difficulties in connecting the whole network and affects the overall cost of avoided emissions. The emissions and costs associated with battery-electric technologies are highly dependent on those of the electric grid and therefore sensitive toward the locations of the charging facilities. To reduce emissions by 1 kg of CO2, it costs Western railroads $ 0.11 and Eastern railroads $ 0.09 under the examined deployment scenario.

Example of battery-electric deployment for Western railroads.

Example of battery-electric deployment for Eastern railroads.
Figures 6 and 7 show a corresponding example deployment of a hydrogen locomotive technology (with an approximate 1,000-mi range) for Western and Eastern railroads, respectively, needed to serve approximately 55% of their ton-mileage. Though Western rail networks are in general more expansive than Eastern rail networks, the long range of this hydrogen locomotive diminishes the number of required refueling facilities compared with the battery case (nine for Western versus three for Eastern). Most notably, hydrogen locomotive operations are not considerably cleaner than diesel operations, as hydrogen fuel is primarily produced through natural gas reforming. Furthermore, liquid hydrogen fuel exhibits high costs because of the energy-intensive process of onsite liquefaction. The relatively low emissions reductions and high incremental costs lead to high costs of avoided emissions. To reduce emissions by 1 kg of CO2 through hydrogen operations, it costs Western railroads $ 0.9 and Eastern railroads $1.05 under the examined deployment scenario. More environmentally friendly (e.g., solar, renewable, nuclear powered) and economical processes for hydrogen production would be required to help hydrogen decarbonize the rail freight sector.

Example of hydrogen deployment for Western railroads.

Example of hydrogen deployment for Eastern railroads.
Biofuel and E-Fuel
For drop-in fuels such as biofuels and e-fuels, fuel blends are assumed to be applied uniformly across all locomotives on the network. Figures 8 and 9 show the results for 50% deployment of biofuels in Western and Eastern railroads, respectively. This deployment of biodiesels would contribute to a 36% reduction in emissions (for both railroad groups) relative to diesel, with a cost of $ 0.13 per kilogram of CO2 reduced. Figures 10 and 11 show the results for 50% deployment of e-fuels in Western and Eastern railroads, respectively. As e-fuels are nearly carbon neutral, they would provide a more promising environmental solution than biofuels, albeit at significantly greater cost (nearly double that of conventional diesel). Their deployment in this scenario would contribute to a 50% decrease in carbon emissions at a cost of $ 0.22 per kilogram of CO2 eliminated.

Example of 50% biodiesel deployment for Western railroads.

Example of 50% biodiesel deployment for Eastern railroads.

Example of 50% e-fuel deployment for Western railroads.

Example of 50% e-fuel deployment for Eastern railroads.
Scenario Comparison
This section evaluates the optimization and simulation results of the battery-electric deployments, demonstrating the framework’s potential for analyzing and comparing across scenarios.
Deployment Percentage
A key functionality of our framework implemented in the dashboard allows users to input different deployment percentages to simulate, optimize, and evaluate different intermediate stages of a technology’s roll-out.
Figures 12 to 14 show the results for Eastern railroads under (approximately) 30%, 50%, and 100% target deployments. As the deployment percentages increase, the network coverage increases drastically, with the number of facilities growing from 12, through 21 to 167 at 100%. This highlights the high number of facilities required to serve the “last miles” of the rail network. While emissions decrease rather proportionally with the extent of the roll-out, the initial LCO at 30% deployment is considerably high, because of starting costs on capital infrastructure projects. Higher deployment percentages enable economies of scale to reduce the costs of avoided emissions from $ 0.17 per kilogram of CO2 at 30% to $ 0.09 at 50%.

Results for Eastern railroads with 30% target deployment percentage.

Results for Eastern railroads with 50% target deployment percentage.

Results for Eastern railroads with 100% target deployment percentage.
Range
In the context of battery-electric deployment, locomotive range affects both the economics and environmental performance of a particular scenario. Figures 15 to 17 illustrate the key trade-off associated with locomotive range. Increasing the locomotive range increases the energy storage capacity per locomotive and thus increases the total battery purchase and operating costs. However, an increase in locomotive range allows for greater reach and fewer charging stations to be deployed, reducing the total facility capital costs, while increasing the network penetration of service that can be provided. This trade-off is illustrated by the costs of avoided emissions which are $ 0.10 per kilogram of CO2 at 200-mi range, $ 0.11 at 400-mi range, and $ 0.06 at 800-mi range. The stark decrease in the cost of avoided emissions between the 400-mi and 800-mi range cases comes from the consolidation of freight along key corridors that allows emissions reductions to go from 16% with 400-mi range locomotives to 46% with 800-mi range locomotives. Furthermore, locomotives with longer ranges (e.g., 800 mi) can significantly reduce emissions, as they can be used to decarbonize more energy-intensive commodities (i.e., intermodal) that are typically shipped over long distances. Note that with locomotives with a 200-mi range, though a 50% target deployment level was set, only 31% of ton-miles could be served by battery-electric locomotives, because of the insufficient range on the expansive Western railroad network. The flexibility of the range parameter supports the sensitivity analysis of optimal technology range values.

Results for Western railroads with 200-mi range.

Results for Western railroads with 400-mi range.

Results for Western railroads with 800-mi range.
Conclusion
This paper presents a simulation framework to aid stakeholders to optimize, plan, and analyze deployment strategies for alternative fuel propulsion technologies in freight rail. It is applied to evaluate different classes of energy technologies that can be used to decarbonize freight rail such as biofuels, e-fuels, hydrogen, and battery-electric pathways. The five-step approach developed in this work builds on nominal graph theory problems to simplify and solve otherwise intractable problems in the field of facility location and sizing. This framework can efficiently simulate a variety of technology adoption scenarios with built-in LCA and TEA tools and features input flexibility for evaluating a variety of deployment scenarios and performing sensitivity analysis of key technological and operational parameters. Key metrics covering emissions, costs, facility locations and sizes, and traffic volumes are outputted in aggregate and granular levels for detailed analysis.
The framework’s functionalities are demonstrated with US Class I railroad data for Eastern and Western railroad networks and technological parameters for battery-electric, biofuel, and e-fuel technologies, over scenarios with various deployment percentages and ranges. Drop-in fuel deployments are modeled as admixtures with diesel in existing locomotives, while battery-electric deployments are shown for varying technology penetration levels and locomotive ranges. A 50% admixture ratio of diesel with biodiesel is estimated to provide a 36% emissions reduction with a cost of $ 0.13 per kilogram of CO2 reduced, while a similar mixture with e-fuels would cut emissions by 50% at a cost of $ 0.22 per kilogram of CO2 reduced. Battery-electric technology deployments at 50% of all ton-miles highlight the need for technological developments in battery energy densities as scenarios for 800-mi range locomotives are estimated to provide a 46% emissions reduction at a cost of $ 0.06 per kilogram of CO2 reduced, compared with a 16% emissions reduction at a cost of $ 0.11 per kilogram of CO2 reduced for 400-mi range locomotives.
As the deployment of alternative fuel technologies would pose considerable operational impacts to the railroads, future work should incorporate the potential delay, congestion, additional fleet (locomotives and rail cars), and track infrastructure maintenance costs arising from possible requirements and re-routing options associated with a particular deployment strategy. Furthermore, this paper is limited in its discussion of future technological parameter uncertainties and the effects on technology roll-outs. Future work will conduct sensitivity analysis on key values for capital, operations and maintenance, locomotive and energy storage, and charging/refueling costs, technological parameters for energy storage tender car configurations, and emissions intensity values.
The presented framework provides a powerful tool for simulating, evaluating, and analyzing future alternative energy technology deployments for freight rail decarbonization. The current version of this tool is a first step in bridging the gap between scientific research and implementation of efforts to reduce the environmental impact of hard-to-decarbonize transportation modes.
Footnotes
Acknowledgements
We thank all ARPA-E personnel who have provided insight and expertise that greatly assisted the research, especially Robert Ledoux and Mirjana Marden. We also thank Industry Advisory Board members and Breton Johnson who were part of the overall project and provided comments on the work.
Author Contributions
The authors confirm contribution to the paper as follows: study conception and design: A. Hernandez, M. T. M. Ng, N. Siddique, P. L. Durango-Cohen, A. Elgowainy, H. S. Mahmassani, M. Wang, Y. Zhou; data collection: A. Hernandez, M. T. M. Ng, N. Siddique; analysis and interpretation of results: A. Hernandez, M. T. M. Ng, N. Siddique, P. L. Durango-Cohen, A. Elgowainy, H. S. Mahmassani, M. Wang, Y. Zhou; draft manuscript preparation: A. Hernandez, M. T. M. Ng, N. Siddique, P. L. Durango-Cohen, A. Elgowainy, H. S. Mahmassani, M. Wang, Y. Zhou. All authors reviewed the results and approved the submission of the manuscript.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This paper is based on research funded by Advanced Research Projects Agency-Energy (ARPA-E) as part of a project titled “LOwering CO2: Models to Optimize Train Infrastructure, Vehicles, and Energy Storage” (LOCOMOTIVES) conducted by the Northwestern University Transportation Center (NUTC) in collaboration with Argonne National Laboratory (ANL).
