 Research
 Open Access
 Published:
Stochastic technoeconomic analysis of alcoholtojet fuel production
Biotechnology for Biofuelsvolume 10, Article number: 18 (2017)
Abstract
Background
Alcoholtojet (ATJ) is one of the technical feasible biofuel technologies. It produces jet fuel from sugary, starchy, and lignocellulosic biomass, such as sugarcane, corn grain, and switchgrass, via fermentation of sugars to ethanol or other alcohols. This study assesses the ATJ biofuel production pathway for these three biomass feedstocks, and advances existing technoeconomic analyses of biofuels in three ways. First, we incorporate technical uncertainty for all byproducts and coproducts though statistical linkages between conversion efficiencies and input and output levels. Second, future price uncertainty is based on casebycase timeseries estimation, and a local sensitivity analysis is conducted with respect to each uncertain variable. Third, breakeven price distributions are developed to communicate the inherent uncertainty in breakeven price. This research also considers uncertainties in utility input requirements, fuel and byproduct outputs, as well as price uncertainties for all major inputs, products, and coproducts. All analyses are done from the perspective of a private firm.
Results
The stochastic dominance results of net present values (NPV) and breakeven price distributions show that sugarcane is the lowest cost feedstock over the entire range of uncertainty with the least risks, followed by corn grain and switchgrass, with the mean breakeven jet fuel prices being $0.96/L ($3.65/gal), $1.01/L ($3.84/gal), and $1.38/L ($5.21/gal), respectively. The variation of revenues from byproducts in corn grain pathway can significantly impact its profitability. Sensitivity analyses show that technical uncertainty significantly impacts breakeven price and NPV distributions.
Conclusions
Technical uncertainty is critical in determining the economic performance of the ATJ fuel pathway. Technical uncertainty needs to be considered in future economic analyses. The variation of revenues from byproducts plays a significant role in profitability. With the distribution of breakeven prices, potential investors can apply whatever risk preferences they like to determine an appropriate bid or breakeven price that matches their risk profile.
Background
Aviation currently accounts for approximately 5% of total anthropogenic radiative forcing [1, 2]. In the absence of mitigation measures, total greenhouse gas (GHG) emissions associated with aviation are expected to be 400–600% higher in 2050 than in 2010, driven by an increase in global air traffic of up to seven times [3]. Against this backdrop, the International Air Transport Association [4] (IATA) has set a goal of carbonneutral growth of aviation by 2020, and a 50% reduction of CO_{2} emissions by 2050 compared to 2005 levels. Similarly, the United States (US) Federal Aviation Administration [5] (FAA) aims for carbonneutral growth of aviation by 2020. These goals are to be achieved by improvements in aircraft operations, airport and air traffic management, airframe and engine technologies, as well as through the largescale introduction of biofuels with significantly lower GHG emissions than petroleumderived jet fuel, on a life cycle basis [6]. To date, no mandate exists specifically for aviation biofuel usage; however, these fuels can qualify under the Renewable Fuel Standard (RFS). Moreover, the US FAA has set a shortterm goal of 1 billion gallons of alternative fuel consumption by 2018 for military and commercial applications [5].
Reduction in the climate impact of aviation may be achieved via the use of biofuels. However, unlike ground transportation which can transition to ethanol or electricity, aviation requires the use of energy dense, nonoxygenate, hydrocarbon, liquid fuels [7]. There are four major aviation biofuel technologies that are currently technically feasible: Fischer–Tropsch (F–T), hydroprocessed renewable esters and fatty acids (HEFA), sugar conversion (fermentation, thermochemical), and direct liquefaction (pyrolysis) [8]. In addition to the potential climate benefits, aviation biofuel production could help to meet the 36 million RFS targets by 2022, and could help reduce US dependence on energy imports and increase energy security [9]. More than twenty airlines have already used aviation biofuels blended with petroleumderived jet fuel on thousands of passenger flights [10].
The existing biofuels TEA literature focuses mainly on bioethanol and biodiesel production. Recent biodiesel TEA literature focuses on vegetable oils for carbon chain attributes similar to petroleum diesel [11–13]. Other existing biofuel TEA literature emphasizes bioethanol production from lignocellulosic biomass, because lignocellulosic feedstocks have lower expected feedstock costs and avoid direct competition with food [14–19]. Generally, the TEA literature calculates breakeven prices, internal rates of return (IRR), and net present values (NPV), and uncertainty has been incorporated in a number of studies in order to estimate distributions of these values. Bauer and Hulteberg [20] developed a probability distribution for production cost using Monte Carlo simulation when evaluating a new thermochemical production process for isobutanol. Abubakar et al. [21] graphed the variations of mean NPV with the increase of the sample size. Sensitivity analyses conducted by Reyes Valle et al. [22] estimated how breakeven prices respond to ±30% uncertainty in fixed capital costs. Zhu et al. [23] used a sample size of 100 experimental cases to derive a breakeven price distribution when evaluating a Benchscale woody biomass hydrothermal liquefaction (HTL) upgrading plant; however, their sample size is insufficient to estimate a breakeven price distribution without randomization and the authors did not consider how price projections and price uncertainties would influence the distribution results.
Very little existing TEA literature focuses specifically on aviation biofuel production, and most studies in the literature are deterministic. Atsonios et al. [24] modeled the ATJ process and evaluated five pathways of converting corn stover and wheat straw to aviation fuels deterministically. They obtained a $1.39/L breakeven price for an F–T plant, which is lower than for a mixed alcohols synthesis (MAS) plant. They concluded that the expected breakeven price of ATJ is higher, despite better performance in terms of carbon utilization and thermal efficiency, than the F–T Synthesis (FTS) route. Staples et al. [25] calculated breakeven prices of renewable middle distillate (diesel and jet) fuels from fermentation and advanced fermentation technologies, using sugarcane, corn grain, and switchgrass as feedstocks. The authors employed three scenarios and found that breakeven prices for sugarcane, corn grain, and switchgrass range from $0.61 to 2.63, $0.84 to 3.65, and $1.09 to 6.30/L of middle distillate fuel, respectively. Their analysis showed that breakeven prices are the most sensitive to feedstock type, fuel conversion efficiency, and feedstock costs. Pearlson et al. [26] estimated baseline breakeven prices for HEFA production ranging from $1.01 to 1.16/L. Maximizing jet fuel yield rather than total fuel yield in the HEFA process adds $0.07–0.08/L to the breakeven prices due to the increased hydrogen requirements and reduced middle distillate fuel yield. Similarly, Seber et al. [27] assessed the breakeven price of HEFA middle distillate fuel production from waste oils and tallow. The estimated breakeven prices were $0.88–$1.06/L for yellow grease (YG)derived HEFA and $1.05–1.25/L for tallowderived HEFA. The authors found that feedstock cost contributes the most to breakeven price, and that the breakeven price of middle distillate HEFA from YG and tallow was higher than petroleumderived diesel fuel prices, but lower than the breakeven price of soybean oil HEFA. de Jong et al. [28] compared six shortterm renewable jet fuel pathways by combining possible feedstocks and technologies, as well as ten greenfield, three retrofitting, and nine colocating strategies. Their results showed that HEFA is the most competitive pathway in the short term. However, none of the pathways can compete with petroleumderived jet fuels on a price basis. Their analyses pointed out that conversion efficiency in fermentation is critical in determining breakeven prices. The authors examined the breakeven price and NPV variation ranges in different scenarios of investments, yields, feedstock prices, and hydrogen consumption. However, they did not estimate the distribution patterns of breakeven prices and NPV.
To the best of our knowledge, only one other TEA study for aviation biofuels incorporates stochasticity into key input and output variables: Bittner et al. [29] carried out a stochastic TEA of aviation biofuel from corn stover using a fast pyrolysis process. They investigated policies of reverse auction and capital subsidies, and found that reverse auction is more risk reducing.
This study makes three contributions to the existing biofuel TEA literature. First, most existing stochastic TEA analyses do not integrate the individual uncertain variables with other related variables in the process. We evaluate uncertainty in the conversion efficiency of two steps of the ATJ process, and then link related model variables by statistical estimation to the random draws from distributions of the conversion efficiency factors. The linked variables include capital costs, utility requirements, feedstock quantity, fuel and byproducts output quantity, and natural gas costs, among others.
Second, we employ timeseries price projection based on historical casebycase patterns instead of conventional Brownian motion or mean reversion price assumption. Timeseries estimation captures the uniqueness of the motion processes of each product market, based on historical prices [30].
Third, TEA studies generally translate all the uncertainties into NPV distributions and only calculate the breakeven prices for most likely cases. In this study, we derive breakeven price distributions by considering all combinations of uncertainties. This approach also permits stochastic dominance comparison and gives a guidance of benchmark investing price at each uncertainty level for private investors.
The point of departure for this research is the previous analysis by Staples et al. [25] on renewable middle distillate production via fermentation and advanced fermentation technologies. We extend this work by considering future price projections and introducing technical uncertainties in ATJ production, thereby developing a deeper and more comprehensive understanding of the ATJ pathway.
Methods
Pathway and feedstock descriptions
ATJ involves upgrading of biomassderived alcohols to a dropin jet fuel or blendstock specification. Typically, ATJ technologies extract polymer sugars from a biomass feedstock via mechanical, chemical or biological means. The polymer sugars are then decomposed to monomer sugars, and metabolized (or fermented) by an engineered microorganism to an alcohol platform molecule (ethanol or isobutanol). Finally, the alcohol is dehydrated, oligomerized, and hydrogenated to a final fuel product slate which includes some proportion of dropin jet fuel or blendstock. A number of private corporations, such as Byogy Renewables, Inc. and Gevo, Inc., has been pursuing ASTM certification and commercialization of ATJ technologies. Gevo’s ATJ production has been approved by ASTM standard in March, 2016 and up to a 30% blend in conventional jet fuel is anticipated to be used for commercial flights [31]. The subject of this analysis is a subset of ATJ technologies, which includes sugars derived from sugarcane, corn grain or switchgrass, followed by fermentation to an ethanol platform molecule. These feedstocks are selected to represent the present and future of renewable fuel production: corn grain and sugarcane are commonly used for the production of ethanol in the US and Brazil, respectively, and herbaceous lignocellulosic crops, such as switchgrass, can be used for the production of secondgeneration renewable fuels such as cellulosic ethanol. The final fuel product slate includes diesel, jet, heavy fuel oil, and naphtha, and we also consider nonfuel coproducts from the ATJ process. ATJ derived from corn grain results in the coproduction of distiller dry grains and solubles (DDGS). Bagasse produced after juice extraction from sugarcane, and biomass residues generated after sugar extraction and fermentation from switchgrass can be cofired to meet the utility requirements of the biorefinery, and excess electricity can be exported to the grid [25]. A simplified schematic of the ATJ process is shown in Fig. 1.
Model framework and basic assumptions
Our model is designed to capture and quantify variation in ATJ process inputs, fuel outputs, and coproducts. Mass and energy balance relationships, the base case scenario and the range of feasible parameter values are derived from Staples et al. [25], where the base case is defined as the most likely or mode value. We present ATJ by two main process steps: feedstocktoethanol conversion and ethanoltofuel conversion. Both steps require water, electricity, and heat (generated from natural gas) inputs. Two conversion efficiency factors are developed corresponding to the two steps, denoted as C _{fset} and C _{etfl} in Fig. 2, and the product of the two conversion efficiency factors is overall conversion efficiency (see next section for more details). The two conversion efficiency factors link feedstock inputs in with fuel outputs and drive variation in the utility requirements, quantities of coproduct generated, and capital costs associated with the ATJ process. C _{fset}, C _{etfl} and other price variables have independent stochastic distributions, represented by ovals in Fig. 2. Each iteration of the Monte Carlo simulation yields a random value from each independent stochastic distribution, and drives the changes of variables shown as rectangles in Fig. 2. Variables shown as parallelograms, such as water, power, and other inputs (enzymes, yeast, and chemicals), are less than 0.01, 0.1, and 1% of the total costs for each feedstock, respectively. Their variations do not significantly impact calculated NPV and breakeven price distributions. We treat them as exogenous and deterministic. We use @Risk, an excel addin software, to perform Monte Carlo simulations [32].
All the price projections and breakeven price distributions are presented in real dollars. Financial analysis in this study first conducted in nominal terms and then converted to real. The deterministic assumptions in this analysis are taken from previous research by Staples et al. [25] and Seber et al. [27] assuming a facility size of 4000 bpd with 8400 operating hours per year. We assume a construction period of 3 years, followed by 20 years of production, and 8, 60, and 32% of the initial fixed capital investments are expended during the construction phase. We assume 20% equity and 80% of capital investment, financed through loans at a 5.5% interest rate for the first 10 years, and working capital is calculated as 20% of first production year (4th project year) operating costs. Since working capital is added back in the last production year, the only financial cost is the implicit interest cost of the working capital advance. We adopt the variable declining balance (VDB) depreciation method for the first 10 production years. The nominal discount rate is 15%; the income tax rate is 16.9%; and all values are presented in 2012 US dollars.
Technical uncertainty
Conversion efficiency
C _{fset} and C _{etfl} reflect the conversion efficiency of the feedstocktoethanol and ethanoltofuel processes, respectively. The three feedstocks considered in this analysis have different feedstocktoethanol conversion factors, but share the same ethanoltofuel conversion factor. The two conversion efficiency factors are expressed in units of kg feedstock per kg of ethanol and kg ethanol per MJ fuel, respectively, and the product of the two is the overall conversion efficiency factor in units of kg feedstock per MJ fuel.
Both the feedstocktoethanol and ethanoltofuel conversion factors are bounded and assumed to follow a PERT distribution. The PERT distribution shares the same parameters as a triangular distribution (defined by min, mode, and max values), but more of the probability density is located around the mode than a triangular distribution. The min, mode, and max values are obtained from Staples et al. [25] original technical estimation work, and the mean value of the PERT distribution is calculated as (min + 4*mode + max)/6. The min, mode, max and mean values of the low, base, and high cases are shown in Table 1.
We assume that the total final fuel output quantities are the same for all three feedstocks, and we use statistical regressions to link the two conversion efficiency factors with feedstock inputs, utility requirements, and the share of each fuel for total fuel output. Therefore, both inputs and outputs are varied based on random draws of the two conversion efficiency factors generated in the Monte Carlo simulation.
The feedstocktoethanol process includes preprocessing, saccharification, and fermentation process steps. In each of these three subprocesses, the electricity, water, and heat utility requirements and output fuel shares, are correlated to the two conversion factors, C _{fset} and C _{etfl}, as well as the interaction between the two conversion efficiency factors. In the interaction terms, C _{fset} takes either quadratic or linear form and all of the resulting regression equations are significant with R ^{2} values over 0.98. The resulting equations are:
The ethanoltofuel process consists of separation and postprocessing. In each of these subprocesses, utility inputs of electricity, water, and heat, and the output fuel product shares, are determined by a quadratic function of C _{etfl}:
A detailed list of regressions for each utility input in each subprocess is presented in the Additional file 1: Table A1. Feedstock inputs are calculated from the input–output mass balances, and determined by the overall conversion efficiency factor. Through the three equations shown above, all input and output quantities are subject to variations in the two conversion efficiency levels.
Capital cost
Uncertainty in capital investment presents another aspect of technical uncertainty. Capital cost consists of two components: preprocessing and fermenter costs, and dehydration, oligomerization, and hydrotreating costs. Dehydration, oligomerization, and hydrotreating costs are treated as a linear function of facility size. Feedstock preprocessing and fermenter costs are a function of feedstock input quantity and dollarsperunitmass of feedstock processing capacities estimated from Staples et al. [25]. For sugarcane, the range is from $20 to 30/kg of capacity [33, 34], for corn grain the range is $55–95/kg of capacity [35, 36], and for switchgrass, the range is $115–215/kg of capacity [36, 37]. Since these capital costs are also bounded, we again choose a PERT distribution for the stochastic analysis. The modes of the preprocessing and fermenter capital cost distributions for corn grain, sugarcane, and switchgrass are $300, $347, and $697 million, respectively. The total capital cost distribution for corn grain and sugarcane follow a Beta General distribution with 90% of the values falling into the range $261–341 and $305–390 million, respectively. The total capital costs for switchgrass follows a gamma distribution with 90% of the values falling into the range from $537 to 899 million. The capital costs of preprocessing and fermenter capacity are lowest for corn grain, followed by sugarcane and switchgrass. Corn grain preprocessing is well established and is feedstock intensive; sugarcane milling involves handling the bagasse coproduct; and switchgrass is a lignocellulosic process involving handling large volumes of feedstock material, as well as costly feedstock preprocessing steps.
Price uncertainty
The future prices of the three biomass feedstocks, natural gas inputs, and diesel are projected with uncertainty. We employ two major price estimation methods: casebycase timeseries estimation is used for corn, sugarcane feedstock prices, natural gas prices, and diesel prices and contractbased price estimation, indexed by yield, is used for switchgrass prices.
Timeseries price estimation
Future price projection is a central challenge for stochastic TEA, and in much of the literature either Brownian motion or mean reversion techniques are employed. However, neither approach is completely satisfactory: Meade [38] compared Brownian motion and mean reversion by examining daily Brent and West Texas Intermediate (WTI) crude oil prices via density forecasts. He found that Brownian motion is only accurate for one or two years, and that the addition of mean reversion does not improve the performance of the model. Postali and Picchetti [39] found that mean reversion is more accurate in representing the evolution of oil prices over time without considering structural breaks, and that geometric Brownian motion (GBM) had fewer evaluation errors with low mean reversion rate. GBM may be a better choice when no reverting trend is apparent, otherwise mean reversion is a superior choice [39]. Lucia and Schwartz [40] proposed three mean reversion models with jumps and spikes when studying energy commodity prices. He found that a price derived from a proper jumpdiffusion model is closer to market price data than the GBM model in the short term. Petter and Tyner [41] found that mean reversion is a more appropriate method for diesel and gasoline price projections. From this review of the literature, there is no consistent conclusion about which method is preferred for estimating future fossil fuel prices. In addition, the motion processes underlying price movements may be different for unique commodity markets.
Given the existence of mature markets for all of the nonswitchgrass inputs and outputs for the ATJ process, future prices can be projected using historical price data. Historical prices can also be used for the fuel products of the pathway, because the renewable fuels produced via the ATJ pathway have very similar performance characteristics to their petroleumderived analogues. Therefore, we assume ATJderived and petroleumderived fuels to be fungible products, up to a blend of 50% ATJ, with identical market prices. We go beyond previous analyses by employing timeseries estimation using historical price data for each commodity price series, in order to forecast future feedstock, natural gas, and fuel product prices. Historical data for each commodity price are tested in order to determine the timeseries process that best fits each commodity.
Corn grain and sugarcane are commodities with mature markets, and annual historical prices from 1980 to 2014 are available from the US Department of Agriculture [42, 43]. Based on Akaike information criterion (AIC), the secondorder moving average (MA2) turns out to be the best price projection for corn grain and sugarcane by following the form [44]:
where (1) P _{ t } is the corn grain or sugar prices in time t; (2) μ = E(P _{ t }); (3) ɛ _{ t } = σN _{ t }, σ is the volatility parameter, and \(N_{t} \,\sim\,{\text{Normal}} \left( {0,1} \right)\); (4) Var(P _{ t }) = σ ^{2}(1 + b _{1} ^{2} + b _{2} ^{2} ), and b _{1} and b _{2} are the moving average coefficients.
The upper bounds for corn grain and sugar price timeseries simulation are approximately identical to their maximum historical prices, while the lower bounds sometimes generate negative values. Since negative commodity prices are unrealistic, we truncate each year’s price distribution at 0.75 times their minimum historical prices, and the fraction of the lower bound tails generated by truncation is negligible. Sugar prices are converted to sugarcane prices assuming a yield of 1 kg raw sugar from 10 kg of sugarcane [25]. All parameter estimates are presented in Table 2.
DDGS is a byproduct of ATJ pathway from corn grain, and its prices are positively correlated with corn grain prices. It is an important revenue source in the corn grain ATJ case. We use a simple ordinary least square (OLS) regression to represent the relationship between prices of DDGS and corn grain prices with a R ^{2} of 0.87:
Natural gas accounts for over 90% of utility input costs in the base case ATJ for all three feedstocks. Natural gas is used for both heat and hydrogen production. Therefore, the variability in natural gas prices make the profitability of ATJ production more uncertain. Timeseries estimation based on historical prices since 1997 is used to project future natural gas prices. AIC criterion indicates that the firstorder moving average process (MA1) is the timeseries stochastic projection method with the best fit, following Eq. (6), defined by the parameters shown below:
where (1) P _{ t } is the natural gas prices in time t; (2) μ = E(P _{ t }); (3) ɛ _{ t } = σN _{ t }, σ is the volatility parameter, and \(N_{t} \,\sim\,{\text{Normal}} \left( {0,1} \right)\); (4) Var(P _{ t }) = σ ^{2}(1 + b _{1} ^{2} ), and b _{1} is the moving average coefficient.
Similar to corn grain and sugar prices, the natural gas price distributions are truncated on the low end at 0.75 times of the minimum historical prices in order to avoid negative prices, and are converted to units of 2012 US dollars per MJ.
Sale of fuel products is the major revenue stream for the ATJ pathway. In addition to variation in the quantity of fuel produced, driven by the two conversion efficiency factors, future prices for jet, diesel, naphtha, and heavy fuel oil are also uncertain. Diesel prices are forecasted using timeseries estimation, and jet, naphtha, and heavy fuel oil prices are calculated on the basis of their historical correlation with diesel prices.
Future diesel prices follow a firstorder autoregressive moving average (ARMA11) process shown in Eq. (7), following the parameter estimates shown in Table 2 [45].
where (1) P _{ t } is the diesel prices in time t; (2) μ = E(P _{ t }); (3) ɛ _{ t } = σN _{ t } and σ is the volatility parameter, and \(N_{t} \,\sim\,{\text{Normal}} \left( {0,1} \right)\); (4) Var(P _{ t }) = σ ^{2}(1 + b _{1} ^{2} + 2a _{1} b _{1})/(1 − a _{1} ^{2} ), a _{1} is the autoregressive coefficient, b _{1} is the moving average coefficient.
Historical data demonstrate that jet and diesel prices are almost identical, with correlations up to 0.996 in some periods. Ordinary least squares regression is used to regress diesel price on jet fuel, and the final regression relationship is
Our analysis also demonstrates that heavy fuel oil and naphtha prices are highly correlated with diesel prices. We link the prices of these products to diesel prices using their historical price ratios.
Contractbased price estimation indexed by yield for switchgrass
In contrast to corn grain and sugarcane feedstocks, switchgrass is not currently a traded commodity, and there are no historical price data to draw upon for price forecasting. Therefore, a different approach is required for this feedstock. The cultivation of switchgrass would require farmers to make a change in their land use for a period of at least 10 years. In order to mitigate risk associated with future revenues, switchgrass producers may choose to operate under longterm price contracts [46]. Significant research exists on contract design to effectively share risk between farmers and biofuel plants. For example, Yoder et al. [46] found that contracts based on dollarsperhectare prices, regardless of yield, were the best option for riskaverse farmers growing miscanthus, a herbaceous cellulosic crop similar to switchgrass. Therefore, this analysis assumes that switchgrass is planted and contracted using fixed dollarsperhectare contracts.
In addition, a number of studies have estimated switchgrass yields under different production conditions. To derive our switchgrass price uncertainty ($/kg), we combine the fixed annual payment ($/ha) with varying annual yield (kg/ha) to estimate the uncertainty in unit switchgrass cost ($/kg).
The yield of switchgrass varies according to the weather conditions each year, and the ecosystem in which the crop is cultivated: switchgrass yields in upland and lowland ecosystems are reported to be distributed with mean (±standard deviation) 8.7 ± 4.2 and 12.9 ± 5.9 1000 kg/ha, respectively [47]. The coefficient of variation (CV) for upland and lowland conditions are 0.483 and 0.457, respectively. The mean of the two CVs is 0.47, the average yield for upland and lowland switchgrass is 10.8 1000 kg/ha, and we use these values to calculate the standard deviation for the average yield, which is 5.08 1000 kg/ha.
We assume the abovederived values for mean and standard deviation of switchgrass yield in order to gauge switchgrass yield uncertainty. To capture a realistic range of real world yields, we assume a bounded PERT distribution that approximates a normal distribution with the above mean and standard deviation. We set the mode to the estimated mean (10.8 1000 kg/ha) and the minimum and maximum values to ±2 standard deviations, leading to a minimum value of 0.6 1000 kg/ha and a maximum of 21.0 1000 kg/ha, respectively. The resulting mean of the PERT distribution is exactly 10.8 1000 kg/ha, with a standard deviation of 3.8 1000 kg/ha.
To derive uncertainty in unit switchgrass feedstock prices ($/kg), we combine the payment from the fixed annual farmer contract ($/ha) with varying annual yield (kg/ha). The average cost of switchgrass is estimated as $116.5/1000 kg according to a report published by the National Academy of Sciences (NAS) [48]. We use this cost together with the yield to calculate the farmer payment ($1258.2/ha):
Using this procedure, we derive the stochastic feedstock price ($/kg) each year, which is the fixed farmer payment ($/ha) from Eq. (9) divided each year by a random draw from the switchgrass yield distribution.
Quantities of the base case for all inputs and outputs and associated prices are presented in Table 3 for an annual production of 232 million liters (61 million gallons), or approximately 4000 bpd (barrels per day), of total fuel production [25].
Breakeven jet price distributions
In addition to NPV distributions, we also develop a way to calculate and present breakeven jet price distributions. Breakeven jet price is the constant real jet price through the entire production period that makes NPV equal to zero. With the variation of the stochastic variables described previously, the diesel and jet prices that make the present value of accumulated revenues equal to the costs also changes. Breakeven price distributions permit potential investors to select any desired risk level, and then to determine the corresponding breakeven price. It also permits comparison among feedstocks.
The basic procedure is to run the standard Monte Carlo simulation and to save all the simulated values. Then the simulated values are used to calculate the breakeven price for each iteration using the Excel Goal Seek function. The breakeven prices are then fit to an appropriate standard distribution. This distribution then can be used to determine the probability for any breakeven price.
Results and discussion
NPV distributions
A summary of NPV distribution results is presented in Table 4. The mean NPV for corn grain, sugarcane, and switchgrassderived ATJ are all negative. Sugarcane has the highest NPV and smallest standard deviation, and switchgrass has the lowest NPV and largest standard deviation (Fig. 3). All three feedstocks’ probability of loss is higher than 85%, and there is more uncertainty in switchgrass ATJ fuel production. We apply stochastic dominance tests to the three distributions and found that sugarcane firstorder stochastic dominates (FSD) corn and corn FSD switchgrass. The definitions of firstorder and secondorder stochastic dominance relationship are introduced in Additional file 1: Section A2. These results imply that under current diesel, jet and feedstock prices, technology levels, and projected future product prices, incentives would be needed to stimulate investment in aviation biofuel production via ATJ (Additional file 1: A2, Figure A1).
The NPV results show that sugarcane is the leastcost option for the ATJ pathway among the three feedstocks considered, under all circumstances. Sugarcane ATJ production does not require heat and electricity utility inputs because cofiring of the coproduced sugarcane bagasse provides more than sufficient heat and power for fuel production, and permits 168 GWh of electricity to be exported to the grid annually, in the base case. Although the combustion of biomass residues generated during switchgrass ATJ production can also offset heat and electricity requirements, it still requires more natural gas and generates less power for export than sugarcane ATJ. In addition, the conversion efficiency of switchgrass ATJ is also lower than corn grain ATJ. In general, we find that the mean NPV of the different renewable jet fuel pathways are inversely proportional to the recalcitrance of simple sugars in the raw feedstock to be converted to ethanol; switchgrass is the most recalcitrant feedstock examined (requiring greater utility, energetic and feedstock inputs per unit of monomer sugar extracted) and has the lowest NPV and, in contrast, sugarcane is the least recalcitrant feedstock (requiring fewer utility, energetic and feedstock inputs per unit of monomer sugar extracted) and has the highest mean NPV.
Breakeven price distributions and policy implications
Fitted breakeven price distributions for corn grain, sugarcane, and switchgrassderived ATJ follow normal, Beta General, and PERT distributions, respectively. The statistics and quintiles of these distributions are presented in Table 5. We find that the breakeven price distribution for switchgrass ATJ has the largest standard deviation, which is because it is represented with higher technical uncertainty than the other two processes.
The stochastic dominance relationship is presented in Fig. 4. The distribution with higher probability to have lower breakeven ATJ fuel prices is more cost efficient. By definition, switchgrass ATJ FSD corn grain and sugarcane ATJ. While we find that switchgrassderived jet fuel firstorder stochastically dominates corn and sugarcanederived fuels, corn grain does not with regard to sugarcane, as the cumulative density functions intersect at the 90% probability level (sugarcane only secondorder stochastically dominates corn). This is because DDGS prices increase with corn grain prices, which generates additional revenue when corn grain prices are high. Therefore, at higher feedstock prices, corn grain ATJ is less costly than sugarcane ATJ. However, sugarcane is the best feedstock option in ATJ fuel production under 90% of circumstances analyzed.
The mean [5 percentile; 95 percentile] breakeven jet prices per liter of ATJ from corn grain, sugarcane, and switchgrass are $1.01 [$0.88; $1.15], $0.96 [$0.79; $1.17], and $1.38 [$1.10; $1.81], respectively. The mean values are the price for jet fuel at which investors have a 50% probability of earning more than their threshold discount rate. The breakeven price distributions are within the deterministic range calculated by Staples et al. [25] confirming that our results are consistent with this analysis. Our breakeven prices for corn and sugarcane ATJ are at the lower range of breakeven prices calculated by Pearlson et al. [26] and Seber et al. [27] for HEFA pathways.
We also conduct Welch’s t test to determine whether the three breakeven price distributions are statistically different from each other. The twosample test assuming unequal variances, conducted for three pairwise breakeven price samples, confirms that the mean values of the three breakeven price distributions are significantly different from each other (Additional file 1: A3, Table A2).
From a policyperspective, risk profiles as those developed in this paper can also be used to assess the impact of alternative policies such as loan guarantees, tax credits, crop insurance, end user offtake agreements, reverse auction based on offtake contract, and capital subsidy on reducing project risk [49]. This is especially important given that derisking investment has been shown to be one of the core levers for incentivizing a more rapid scaleup of the aviation biofuel industry [50].
Sensitivity analysis
Figure 5 presents the sensitivity summary for corn grain, sugarcane, and switchgrass ATJ. The results indicate the minimum and maximum values that the NPV can achieve with variation of each individual parameter with the uncertainty ranges assumed in this analysis [51, 52]. The base case NPV is the mean value of NPV distributions with all mode input values. We only report the sensitivity results for the feedstocktoethanol and ethanoltofuel conversion factors, and the feedstock preprocessing and fermentation capital costs. Price uncertainty is not included here because there is a stochastic price variable each year for each price, which cannot be simply aggregated to a single range. The results show that the two conversion factors cause the largest impacts on NPV variation. Corn grain and sugarcane ATJ are most sensitive to ethanoltofuel conversion factors, followed by feedstocktoethanol conversion factors. In contrast, switchgrass ATJ is more sensitive to feedstocktoethanol conversion factors, followed by ethanoltofuel conversion factors. The feedstocktoethanol conversion factor imposes greater uncertainty for switchgrass ATJ, as compared to corn grain and sugarcane ATJ. The feedstocktoethanol conversion factors’ effects on corn grain and sugarcane ATJ are very similar, while its impact on switchgrass ATJ is four times larger than the impacts on corn grain and sugarcane ATJ.
Conclusions
This study makes three contributions to current stochastic TEA: (1) we take technical uncertainty into account by linking conversion efficiency with input and output quantities through statistical methods; (2) in addition to NPV, we develop breakeven price distributions to provide potential investors the price level required to achieve their stipulated rate of return at each probability level; (3) price forecasts are based on casebycase historical timeseries analyses. Sugarcane is the lowest cost feedstock over the entire range of uncertainty with the least risks, followed by corn grain and switchgrass, with the mean breakeven jet fuel prices being $0.96/L ($3.65/gal), $1.01/L ($3.84/gal), and $1.38/L ($5.21/gal), respectively. The probability of loss given the future fuel market price projections for sugarcane, corn grain, and switchgrass ATJ are 88, 95, and 100%, respectively. Price support policies based on breakeven price distributions should be implemented to avoid potential losses and achieve targeted profitability.
Incorporating both technical and economic uncertainty is critical in characterizing the economic performance of any new technology and needs to be considered in future economic analyses. We find that the variation of revenues from byproducts can impact profitability differently at different probability levels.
Abbreviations
 ATJ:

alcoholtojet
 ARMA11:

firstorder autoregressive moving average
 DDGS:

distiller dry grains and solubles
 FAA:

Federal Aviation Administration
 F–T:

Fischer–Tropsch
 FTS:

F–T synthesis
 GBM:

geometric Brownian motion
 GHG:

greenhouse gas
 HEFA:

hydroprocessed renewable esters and fatty acids
 HTL:

hydrothermal liquefaction
 IATA:

International Air Transport Association
 IRR:

internal rate of return
 MA1:

firstorder moving average
 MA2:

secondorder moving average
 MAS:

mixed alcohols synthesis
 NPV:

net present value
 OLS:

ordinary least squares
 RFS:

renewable fuel standard
 TEA:

technoeconomic analysis
 YG:

yellow grease
References
 1.
Lee DS, Fahey DW, Forster PM, Newton PJ, Wit RC, Lim LL, Owen B, Sausen R. Aviation and global climate change in the 21st century. Atmos Environ. 2009;43(22):3520–37.
 2.
Dessens O, Köhler MO, Rogers HL, Jones RL, Pyle JA. Aviation and climate change. Transp Policy. 2014;34:14–20.
 3.
ICAO. Environmental report. Montreal: Destination Green; 2013.
 4.
IATA. A global approach to reducing aviation emissions. First stop: carbonneutral growth from 2009. http://www.iata.org/whatwedo/environment/Documents/globalapproachreducingemissions.pdf. Accessed 9 Sept 2015.
 5.
Federal Aviation Administration. Destination 2025. Federal Aviation Administration, U.S. Department of Transportation, Washington, DC, 2011. https://www.faa.gov/about/plans_reports/media/Destination2025.pdf. Accessed 9 Sept 2015.
 6.
Rahmes TF, Kinder JD, Henry TM, Crenfeldt G, LeDuc GF, Zombanakis GP, Abe Y, Lambert DM, Lewis C, Juenger JA. Sustainable bioderived synthetic paraffinic kerosene (BioSPK) jet fuel flights and engine tests program results. Report No AIAA 2009, 7002.
 7.
Air Transport Action Group. Beginner’s guide to aviation biofuels. Cary: Aviation’s Envrionmental Resource Center; 2009.
 8.
Faaij A, van Dijk M. White paper on sustainable jet fuel. Sky energy white paper 2013, 12.
 9.
Renewable Fuel Standard (RFS). http://www.epa.gov/oms/fuels/renewablefuels/.
 10.
Passenger Biofuel Flights. http://aviationbenefits.org/environmentalefficiency/sustainablefuels/passengerbiofuelflights/.
 11.
Apostolakou A, Kookos I, Marazioti C, Angelopoulos K. Technoeconomic analysis of a biodiesel production process from vegetable oils. Fuel Process Technol. 2009;90(7):1023–31.
 12.
Tang ZC, Zhenzhou L, Zhiwen L, Ningcong X. Uncertainty analysis and global sensitivity analysis of technoeconomic assessments for biodiesel production. Bioresour Technol. 2015;175:502–8.
 13.
Koutinas AA, Chatzifragkou A, Kopsahelis N, Papanikolaou S, Kookos IK. Design and technoeconomic evaluation of microbial oil production as a renewable resource for biodiesel and oleochemical production. Fuel. 2014;116:566–77.
 14.
Gubicza K, Nieves IU, Sagues WJ, Barta Z, Shanmugam KT, Ingram LO. Technoeconomic analysis of ethanol production from sugarcane bagasse using a Liquefaction plus Simultaneous Saccharification and coFermentation process. Bioresour Technol. 2016;208:428
 15.
Hamelinck CN, Van Hooijdonk G, Faaij AP. Ethanol from lignocellulosic biomass: technoeconomic performance in short, middleand longterm. Biomass Bioenergy. 2005;28(4):384–410.
 16.
Ljunggren M, Wallberg O, Zacchi G. Technoeconomic comparison of a biological hydrogen process and a 2nd generation ethanol process using barley straw as feedstock. Bioresour Technol. 2011;102(20):9524–31.
 17.
KleinMarcuschamer D, Simmons BA, Blanch HW. Technoeconomic analysis of a lignocellulosic ethanol biorefinery with ionic liquid pretreatment. Biofuels Bioprod Biorefin. 2011;5(5):562–9.
 18.
Gnansounou E, Dauriat A. Technoeconomic analysis of lignocellulosic ethanol: a review. Bioresour Technol. 2010;101(13):4980–91.
 19.
Sassner P, Galbe M, Zacchi G. Technoeconomic evaluation of bioethanol production from three different lignocellulosic materials. Biomass Bioenergy. 2008;32(5):422–30.
 20.
Bauer F, Hulteberg C. Isobutanol from glycerine—a technoeconomic evaluation of a new biofuel production process. Appl Energy. 2014;122:261–8.
 21.
Abubakar U, Sriramula S, Renton NC. Stochastic technoeconomic considerations in biodiesel production. Sustain Energy Technol Assess. 2015;9:1–11.
 22.
Reyes Valle C, Villanueva Perales AL, VidalBarrero F, GómezBarea A. Technoeconomic assessment of biomasstoethanol by indirect fluidized bed gasification: impact of reforming technologies and comparison with entrained flow gasification. Appl Energy. 2013;109:254–66.
 23.
Zhu Y, Biddy MJ, Jones SB, Elliott DC, Schmidt AJ. Technoeconomic analysis of liquid fuel production from woody biomass via hydrothermal liquefaction (HTL) and upgrading. Appl Energy. 2014;129:384–94.
 24.
Atsonios K, Kougioumtzis MA, Panopoulos KD, Kakaras E. Alternative thermochemical routes for aviation biofuels via alcohols synthesis: process modeling, technoeconomic assessment and comparison. Appl Energy. 2015;138:346–66.
 25.
Staples MD, Malina R, Olcay H, Pearlson MN, Hileman JI, Boies A, Barrett SRH. Lifecycle greenhouse gas footprint and minimum selling price of renewable diesel and jet fuel from fermentation and advanced fermentation production technologies. Energy Environ Sci. 2014;7(5):1545–54.
 26.
Pearlson M, Wollersheim C, Hileman J. A technoeconomic review of hydroprocessed renewable esters and fatty acids for jet fuel production. Biofuels Bioprod Biorefin. 2013;7(1):89–96.
 27.
Seber G, Malina R, Pearlson MN, Olcay H, Hileman JI, Barrett SRH. Environmental and economic assessment of producing hydroprocessed jet and diesel fuel from waste oils and tallow. Biomass Bioenergy. 2014;67:108–18.
 28.
de Jong S, Hoefnagels R, Faaij A, Slade R, Mawhood R, Junginger M. The feasibility of shortterm production strategies for renewable jet fuels—a comprehensive technoeconomic comparison. Biofuels Bioprod Biorefin. 2015;9(6):778–800.
 29.
Bittner A, Tyner WE, Zhao X. Field to flight: a technoeconomic analysis of the corn stover to aviation biofuels supply chain. Biofuels Bioprod Biorefin. 2015;9(2):201–10.
 30.
Box GE, Jenkins GM, Reinsel GC, Ljung GM. Time series analysis: forecasting and control. New York: Wiley; 2015.
 31.
Gevo Inc. Gevo’s alcoholtojet fuel meets approved ASTM standard. Englewood: Biomass Magazine; 2016.
 32.
Palisade Corporation. @Risk 6. 2015. http://www.palisade.com/. Accessed 25 July 2015.
 33.
AsiaPacific Economic Cooperation. Biofuel costs, technologies and economics in APEC economies. Singapore: APEC secretariat; 2010.
 34.
Goldemberg J. The Brazilian biofuels industry. Biotechnol Biofuels. 2008;1(6):1–7.
 35.
Entrix C: Current state of the US ethanol industry. Report for US Department of Energy Office of Business Programs Washington, DC Retrieved October 2010, 23:2012.
 36.
Wallace R, Ibsen K, McAloon A, Yee W. Feasibility study for colocating and integrating ethanol production plants from corn starch and lignocellulosic feedstocks (Revised). In: National Renewable Energy Lab., Golden, CO (US); 2005. http://www.nrel.gov/docs/fy05osti/37092.pdf. Accessed 9 Sept 2015.
 37.
Bain RL. World biofuels assessment; worldwide biomass potential: technology characterizations (milestone report). National Renewable Energy Laboratory (NREL), Golden, CO(USA); 2007. http://www.nrel.gov/docs/fy08osti/42467.pdf. Accessed 9 Sept 2015
 38.
Meade N. Oil prices—Brownian motion or mean reversion? A study using a one year ahead density forecast criterion. Energy Econ. 2010;32(6):1485–98.
 39.
Postali FAS, Picchetti P. Geometric Brownian motion and structural breaks in oil prices: a quantitative analysis. Energy Econ. 2006;28(4):506–22.
 40.
Lucia JJ, Schwartz ES. Electricity prices and power derivatives: evidence from the nordic power exchange. Rev Deriv Res. 2002;5(1):5–50.
 41.
Petter R, Tyner WE. Technoeconomic and policy analysis for corn stover biofuels. ISRN economics. 2014; 2014:13. doi:10.1155/2014/515898.
 42.
Sugar and Sweeteners Yearbook Tables. United States Department of Agriculture Economic Research Service. https://www.ers.usda.gov/dataproducts/sugarandsweetenersyearbooktables.aspx. Accessed 24 Oct 2014.
 43.
Statistics by Crops and Plants. United States Department of Agriculture National Agricultural Statistics Service. https://www.nass.usda.gov/Statistics_by_Subject/index.php?sector=CROPS. Accessed 24 Oct 2014
 44.
Palisade Corporation: RISK MA2. @Risk 6 help manual. 2015. http://www.palisade.com/. Accessed 25 July 2015.
 45.
Palisade Corporation: RISK ARMA11. @Risk 6 help manual. 2015. http://www.palisade.com/. Accessed 25 July 2015.
 46.
Yoder JR, Alexander C, Ivanic R, Rosch S, Tyner W, Wu SY. Risk Versus Reward, a Financial Analysis of Alternative Contract Specifications for the Miscanthus Lignocellulosic Supply Chain. Bio Energy Res. 2015;8(2):644–56.
 47.
Wullschleger S, Davis E, Borsuk M, Gunderson C, Lynd L. Biomass production in switchgrass across the United States: database description and determinants of yield. Agron J. 2010;102(4):1158–68.
 48.
Committee on Economic and Environmental Impacts of Increasing Biofuels Production. Renewable fuel standard: potential economic and environmental effects. Washington: National Research Council; 2011.
 49.
Tyner WE, Van Fossen K. Policy research needs relevant to alternative jet fuels. In: Research and development team white paper series: policy impact CAAFI: CAAFI; 2014.
 50.
Gegg P, Budd L, Ison S. The market development of aviation biofuel: drivers and constraints. J Air Transp Manag. 2014;39:34–40.
 51.
Palisade Corporation: Advanced AnalysisTornado. @Risk 6 help manual. 2015. http://www.palisade.com/. Accessed 25 July 2015.
 52.
Van Dael M, Van Passel S, Pelkmans L, Guisson R, Reumermann P, Luzardo NM, Witters N, Broeze J. A technoeconomic evaluation of a biomass energy conversion park. Appl Energy. 2013;104:611–22.
Authors’ contributions
GY designed the technoeconomic model, conducted all the economic analyses in this study, and drafted the paper. MDS developed the mass and energy balance relationships, provided the base case scenario data and the range of feasible parameter values. He also wrote the section of pathway and feedstock descriptions and prepared Fig. 1. RM offered the chemical and economic perspectives. WET supervised this study and provided economic and political insights. All authors edited and reviewed several versions of this paper. All authors read and approved the final manuscript.
Authors’ information
GY is a Ph.D. student in agricultural economics at Purdue University. She obtained her Master’s degree at Purdue as well with research focus on energy economics. MDS is a Ph.D. student in engineering systems at Massachusetts Institute of Technology (MIT). His research focuses on economic and environmental performance of alternative fuels for aviation. RM is a professor of Environmental Economics at Hasselt University in Belgium and a research affiliate at Laboratory for Aviation and the Environment at MIT. WET is James & Lois Ackerman Professor of Agricultural Economics at Purdue University. He is an expert in technoeconomic and energy policy analysis.
Acknowledgements
The authors are pleased to acknowledge funding support from the Federal Aviation Administration for this research. See the funding section.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analyzed during the current study were available from the corresponding author on reasonable request.
Funding
This work was partially funded by the US Federal Aviation Administration (FAA) Office of Environment and Energy as a part of ASCENT Project [107208] under FAA Award Number [13CAJFEPU]. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the FAA or other ASCENT Sponsors.
Author information
Additional file
13068_2017_702_MOESM1_ESM.docx
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Stochastic technoeconomic analysis
 Breakeven price distributions
 Stochastic dominance
 Aviation biofuel
 Alcoholtojet
 Monte Carlo simulation