Process Optimisation of Steam Explosion Parameters on Multiple Lignocellulosic Biomass using Taguchi Method – A Critical Appraisal

Xylitol is a low calorie sweetener that can be produced through a bioconversion approach from lignocellulosic biomass that requires pre-treatment prior to the bioconversion of xylose to xylitol. Steam explosion (SE) is an industrially scalable pre-treatment (PT) process with the potential to liberate xylose monomers, however SE-PT has not been optimized for xylose release from multiple feedstock. The effect of pressure, substrate weight, phosphoric acid loading concentration and residence time on four feedstock [wheat straw (WS), corn stover (CS), Miscanthus (M), and willow (W)] for xylose release and minimal fermentation inhibitor production [furfural and 5-hydroxymethylfurfural (HMF)] was investigated using the Taguchi methodology for design of experiment (DoE) with variation at four levels (4 4 ). An L 16 orthogonal array design was utilized and all factors indicated inﬂuence on xylose release and inhibitor formation and the resulting xylose rich hydrolysate assessed for bioconversion to xylitol. The L 16 DoE gave hydrolysates containing 75–95% of xylose content in the original biomass, whilst retaining cellulose and lignin components in the ﬁber. The level of inhibitors were within boundary limits to enable microbial fermentation of the hydrolysates to xylitol. Fine tuning of the overall evaluation criteria (OEC) model imbibing 1.5 kg feedstock in 1.2% w/v orthophosphoric acid, 12 bar(g) and 6 min residence time resulted in 90% xylose recovery and production of > 1,000 L of wheat straw hydrolysate for bioconversion to xylitol. The advantages and limitations of the Taguchi OEC model and further improvements to this process are discussed in a bioreﬁning context.


INTRODUCTION
Use of fossil fuels has impacted national energy security and global warming, leading to an increased investigation into the use of biomass as a renewable substrate to generate multiple products (Maity, 2014). Key biomass resources include lignocellulosic agricultural residues and dedicated energy crops (2G resources), neither of which impact on food security, unlike 1G plant resources which enter the food v fuel debate. Major global biomass resources include; wheat straw, corn stover, short rotation coppice willow, and Miscanthus. O'Donohue (2014) reported that within Europe, 35 Mt DM wheat straw, 15 Mt DM corn stover, 2.5-5.5 Mt DM surplus hardwood are available for use as biorefining feedstock. In addition, Miscanthus a perennial C4 grass is currently receiving much research attention for use as a dedicated energy crop, bioremediation of soils, and enhancement of grade 3-5 (marginal) land meaning that in the future this could become a widely used crop across the UK and the world. (Heaton et al., 2008;Haughton et al., 2009;O'Donohue, 2014;Rettenmaier et al., 2014;Purdy et al., 2015;Hodgson et al., 2016;Clifton-Brown et al., 2017).
Biomass has predominantly been used for the production of energy to reduce dependence on fossil fuel imports and reduce carbon output, often utilizing biomass for a single product. A European report by O'Donohue (2014) highlighted the need to couple ethanol generation from cellulose with production of multiple products from the hemicellulose and lignin components. An analysis of multiple product portfolios indicated that generation of xylitol from the hemicellulosic fraction made significant improvement to the economic returns of a lignocellulosic based process (O'Donohue, 2014).
Xylitol is a sugar alcohol estimated to achieve a market value of $1 billion by 2020 and is of interest to pharmaceutical and nutraceutical industries due to its multiple health benefits; low glycaemic index, and reduction of dental carries (Venkateswar . The current production of this sweetener is performed via an energy intensive thermo-chemical synthesis route using corn cobs and birch wood hydrolysate. A sustainable production system through biomass pre-treatment to liberate xylan and subsequent bioconversion of xylose to xylitol has previously been considered (Mohamad et al., 2015). To the best of the author's knowledge, an optimized pre-treatment process to release xylose from multiple feedstock without undue release of fermentation inhibitors has not been carried out.
A myriad of pre-treatment techniques to facilitate deconstruction of lignocellulosic feedstock have been developed (e.g., biological, thermo-chemical, and mechanical) but few processes have been shown to be scalable (Alvira et al., 2010;Guerrero et al., 2017). Steam explosion is a simple hydrothermal, scalable, pre-treatment technology developed in 1926 that is industrially deployed to process tones of material per hour and has been used on a variety of feedstock (Mason, 1926;Hu and Ragauskas, 2012;Bryant et al., 2013;Martijn, 2015;Rabemanolontsoa and Saka, 2016). Much of the research utilizing steam explosion pre-treatment has focused on imbibing feedstock with sulfuric acid to serve as a catalyst (Eklund et al., 1995;Zimbardi et al., 2007;Sørensen et al., 2008;Martijn, 2015). Use of this harsh acid requires expensive capital costs for full-scale facilities due to the requirement of equipment to be constructed from exotic steel alloys, such as hastelloy, to withstand acid corrosion (Harmsen et al., 2010;Castro et al., 2014). Employing phosphoric acid as a catalyst enables conventional 316 stainless steel to be used in the reactor design and readily solubilizes the more labile hemicellulose fraction without releasing cellulose and further degrading the sugars to the fermentation inhibitors furfural and 5-Hydroxymethylfurfural (HMF) . For each of the feedstock wheat straw has had a substantial amount of SE-PT attention, as has corn stover, with few on Miscanthus or willow (Horn et al., 2011;Behera et al., 2014;Ding et al., 2018). As both the release of xylose and generation of inhibitors via SE-PT were unknown, multiple feedstock were evaluated in order to, (1) produce hydrolysates compatible for maximum xylitol yield and productivity and (2) address the issue of realistic feedstock supply for a commercial biorefinery.
To produce an economically sustainable SE-PT process, many individual process steps require optimization. Hundreds of experiments are necessary when following a full factorial design methodology leading to extensive labor input, costs, and time, particularly when performed at pilot scale. An alternative is to utilize design of experiment methodology such as the Taguchi method (Rao et al., 2008) that employs orthogonal array experimental designs, thereby reducing the number of runs required to obtain an optimum set of levels. The Taguchi method has been reviewed and critically appraised against other DoE methods such as Response Surface Methodology (RSM) and been shown to require less experiments to generate similar results to RSM (Aggarwal et al., 2008;Rao et al., 2008). Alongside RSM requiring significantly more experiments to achieve similar optimized results, the statistical design relates to number of variables not the statistical factorials that results in Taguchi DoE giving a more fully developed process (Box and Hunter, 1957;Rao et al., 2004). Selecting the key SE-PT factors and dividing these into several levels the Taguchi method provides the experimental conditions required before using statistical analysis to determine, from the results, the optimum level for each parameter to achieve the best xylose release.
The approach adopted in investigating and optimizing phosphoric acid SE-PT at pilot scale to release maximal levels of thermo-labile xylose for higher value xylitol production, differs from current research that has tended to focus on lab scale SE-PT for carbohydrate release for bioethanol production Negro et al., 2014;Sharma et al., 2015). To deliver a successful industrial biotechnological (IB) xylitol production process at commercial scale, multiple obstacles need to be overcome within the IB process that are then optimized and validated at both pilot and demonstration scale. The data presented here address the hurdles of feedstock selection, microbial biocompatibility of the crude hydrolysate, process parameter selection and design space stability of SE-PT. Recently, O'Donohue (2014) highlighted the main issues facing the development of a xylitol IB process were the availability of hydrolysate bioconversion, along with purification and crystallization.
The aim of this study was to identify the effect of four SE-PT parameters; acid loading, substrate weight, pressure, and time, on four feedstock; wheat straw, corn stover, willow, and Miscanthus to attain an optimal set of SE-PT conditions for xylose release and subsequent bioconversion to xylitol at pilot scale. With optimized conditions more than 1,000 l of hydrolysate was produced, concentrated and stored for bio conversions. A critique of the Taguchi method as applied to this pretreatment process is presented and future optimization scenarios are discussed.

Design of Experiment and Analysis
The Taguchi approach used Qualitek 4 DOE software with process factors and levels selected combined into an L 16 orthogonal array defining 16 unique experimental conditions and gave account of all process factor combinations ( Table 1).The 16 experiments were applied to each of the four feedstock with the four factors being selected based on preliminary experimental procedure (data not published) and literature review (Thomsen and Schmidt, 1999;Castro et al., 2014;Guerrero et al., 2017). Data handling and analysis were performed via "bigger is better" performance characteristic utilizing Qualitek 4 DOE software (USA), Microsoft Excel (Microsoft, USA), and SPSS (IBM, USA).

Steam Explosion Pre-treatment
Prior to SE-PT each feedstock was impregnated with 1/5 w/v water imbibed with phosphoric acid at a range of concentrations at 50 • C for 60 min, separated by draining the excess liquid through a crude sieve and fed into the Cambi pilot scale 30 L steam explosion rig (Cambi, Norway) (Figure 1). The reactor was charged to the requisite pressure with steam, maintained for the desired residence time followed by rapidly raising the pressure to 15 bar(g) prior to opening the release valve to maintain consistent explosive decompression. Following SE-PT, the material was cooled to below 50 • C, water-soluble carbohydrates (WSC) were recovered from the solids following addition of 1 L water at 50 • C. The resulting pulp was separation into solid and liquid fractions by filtration through a muslin cloth. The severity factor was calculated using the Overend et al. (1987) method Equation (1).
Where t is the residence time in minutes, T is the temperature of pre-treatment, 100 is the reference temperature and 14.75 is the arbitrary constant ω which is the activation energy from first order kinetics (Pedersen and Meyer, 2010).

Concentration of Steam Exploded Hydrolysate
The wet exploded material was then screw pressed (CP4 Screw press, Vincent Corporation, Tampa, Florida) at 2 bar(g) to separate solids and hydrolysate. Subsequently fiber was washed through a counter current system to extract majority of water soluble carbohydrates. The entire process was repeated until ∼1,000 L of hydrolysate was available, which was subsequently concentrated up. Concentration was adapted from Murthy et al. (2005) briefly; centrifugation through disk stack (Centrimax, Cologne, Germany), before being filtered through 100 kD membrane filter (Axium Process Ltd, Swansea, UK), before finally being concentrated through reverse osmosis (Reverse osmosis plant, Axium Process Ltd, Swansea, UK). Concentration required minimal pH adjustment (NH 4 OH) to ensure hydrolysate remained above pH 2 to allow RO membrane to operate within bounds.

Carbohydrate Analysis
SE-PT hydrolysate was cooled to room temperature and stored until required at −20 • C. Using NREL laboratory analytical procedures carbohydrate content was quantified in both the fiber and hydrolysates fractions (Sluiter et al., 2006(Sluiter et al., , 2008. Briefly, 300 mg (DW) freeze dried feedstock was added to 3 mL 72% H 2 SO 4 and mixed at 30 • C for 60 min, diluted to 4% v/v H 2 SO 4 by the addition of 84 mL deionized water, autoclaved at 121 • C for 60 min, cooled to room temperature and neutralized to pH 6 with CaCO 3 . To quantify the monomer to polymer ratio, liquid hydrolysate from SE-PT was acid hydrolyzed by adding 34.8 µL of 72% H 2 SO 4 to 1 mL of hydrolysate, followed by autoclaving and neutralization as described above.

High Performance Liquid Chromatography
Polymeric, oligomeric and monomeric carbohydrates, and inhibitory compound concentrations were quantified using three HPLC analytical techniques. Monomeric carbohydrates were quantified using calibrated standard curves of arabinose, galactose, glucose, xylose, mannose, fructose ranging from 100 to 0.6 µg mL −1 using a Dionex High Performance Anion-Exchange Chromatography with Pulsed Amperometer Detection (HPAE-PAD) HPLC system. Samples were centrifuged for 5 min at 14,000 rpm, the supernatant removed, diluted to within the range of the standard curve, filter sterilized through Sartorius Ministart 0.2 µm syringe filter (Fisher scientific) into 0.6 mL vials with pre-cut septum (Thermofisher) with 25 µL of each sample injected onto a CarboPac SA10 column, 250 mm long × 4 mm internal diameter and eluted with an isocratic mobile phase of 1 mM KOH at 1.5 mL min −1 (Basumallick and Rohrer, 2015). Quantification of organic acids was performed as in Bryant et al. (2011) using a JASCO HPLC system with a Aminex HPX-87C column, 300 mm long × 7.8 mm internal diameter, connected to a refractive index detector (RID) (Shodex RI-71, Showa Denko Europe GmbH, Munich, Germany) with analytes eluted at 0.6 mL min −1 using a 5 mM H 2 SO 4 isocratic mobile phase and 5 mM crotonic acid as an internal standard. Samples were prepared as above with 25 µL of sample quantified against calibrated standard curves formic acid and acetic acid ranging from 100 to 0.625 mg mL −1 .
Furfural and 5-hydroxymethylfurfural degradation products were quantified on the JASCO HPLC as above by injecting 20 µL of sample onto Aminex HPX-87H, 300 mm long with 7.8 mm internal diameter. Furfural and HMF were run as calibrated standard curves from 5 to 0.3 mg mL −1 (Sluiter et al., 2006).
Xylitol was analyzed on HPLC using JASCO HPLC system with a Aminex HPX-87H column, 300 mm long × 7.8 mm internal diameter, connected to a refractive index detector (RID) (Shodex RI-71, Showa Denko Europe GmbH, Munich, Germany) with analytes eluted at 0.6 mL min −1 using a 5 mM H 2 SO 4 isocratic mobile phase. Samples were quantified against calibrated standard curve for xylitol ranging from 10 to 0.625 mg mL −1 .

Compositional Analysis
The chemical composition of each feedstock prior to SE-PT are described in Table 2 and were found to be in broad agreement with previously published data (Thomsen and Schmidt, 1999;Sun et al., 2005;Maity, 2014;Swiatek et al., 2014). Lignin was highest in Miscanthus and willow similar to previously published reports (Fahmi et al., 2008;Hu and Ragauskas, 2012) while the xylose concentration was lowest in willow due to its higher lignin content (Fahmi et al., 2008), and highest in wheat straw.

Carbohydrate Release Following SE-PT
Pressure, time, acid loading and substrate weight were all shown to alter the concentration of xylose released from each feedstock ( Table 3). Due to the scale at which experiments were being investigated, duplicates were run to ascertain the variation in soluble xylose release from each feedstock following each treatment. A wide range of levels for each factor was investigated aiming to achieve minimal inhibitor production whilst still maximizing xylose release. The highest extraction of soluble xylose for W, CS, M, WS were; 133 mg g −1 DM, 265 mg g −1 DM, 197 mg g −1 DM, and 284 mg g −1 DM respectively (Figure 2). The xylose yield from corn stover and wheat straw was in agreement with the literature (Ballesteros et al., 2006;Fang et al., 2011;Pang et al., 2013), however in these reports, values were obtained using sulphuric acid to catalyze the release of carbohydrates or integrated microwave assistance, both of which increase the capital cost of processing. The variation seen across the experiments (Figures 2, 3)as well as that within duplicates indicate inherent variation of different process parameters experienced, with Auxenfans et al. (2017) finding a 3-26% variation in hemicellulose release within replicates across different combined severity factor experiments. This is exemplified in the optimal conditions of 12 bar(g), 3 min, 500 g, 1.2% w/v for experiment 13 where the error bars are smaller than many of the other high releasing experiments. For industrial scale production the process is required to be robust and consistent with (Figure 2) data highlighting the limitations for industrial use outside of the narrow operating window for SE-PT.
Frontiers in Energy Research | www.frontiersin.org Despite there being a wealth of research on steam explosion, the vast majority of research uses steam explosion as a pretreatment prior to enzymatic hydrolysis to produce both C5 and C6 sugars for bioethanol production. Reports solely focussing on hemicellulose extraction into hydrolysates is not highly reported possibly due to the question of addressing the resulting cellulosic pulp. However, an abundance of literature exists on uses for the remaining fibrous material, for example as substrate for biogas production or enzymatic hydrolysis for bioethanol generation (Horn et al., 2011;López-Linares et al., 2015;Arenas-Cárdenas et al., 2017;Lizasoain et al., 2017). Production of biogas with the remaining pulp could be integrated with a commercial xylitol production plant with the advantages of having a rapid integrated energy supply as SE-PT has been shown to increase the rate of biogas production (Horn et al., 2011).
The broad range of levels resulting in a wide array of severity factors (Table 1) were chosen to assess maximum xylose release whilst decreasing the fermentation inhibitors, several publications have discussed how high severity factors of SE-PT will lead to increased inhibitor generation (Iroba et al., 2014;Guerrero et al., 2017). The data in Figure 2 indicate that low severity factors result in low xylose release. Castro et al. (2014) achieved ∼66% xylose release from Eucalyptus chips utilizing 0.8% w/w phosphoric acid, 200 • C for 5 minu residence time whilst Negro et al. (2014) used 1% w/w phosphoric acid at 175 • C for 10 min to achieve 80% monomeric xylose extraction of olive tree pruning's. Castro et al. (2014) had high levels of inhibitors formed for high xylose release and Negro et al. (2014) generated double the furfural concentration by increasing the temperature from 175 to 195 • C with only moderate xylose increases. However, both were aimed at glucose extraction for fuel ethanol production not maximum xylose release. Zhang et al. (2014) has investigated steam explosion for hemicellulose release using a screw steam explosive extruder to extract almost 90% hemicellulose from corn cobs using sulphuric acid at 9 mg g −1 dried corn cobs, 1.55 MPa (201 • C) for 5.5 min. This work is informative however it used sulphuric acid which increases capital cost and requires the hydrolysate to be detoxified with activated carbon, cation and anion resins. These papers indicate a range of feedstock and the SE-PT parameters required to hydrolyze hemicellulose confirming that optimization of factor levels is required and that each feedstock is independent.

Inhibitor Generation
The production of the inhibitors HMF (from glucose) and furfural (from xylose) was investigated and the combined level of furfural and HMF released through steam explosion of biomass is presented in Figure 3. HMF was not detected within WS, low levels were detected in M and W (<6 mg g −1 DM), and appreciable levels within the CS hydrolysate (∼15 mg g −1 , ∼1 g L −1 ). The highest accumulation of furfural and HMF was found to be 110 mg g −1 DM in CS corresponding to ∼9.6 g L −1 which would be detrimental to bioconversion (van der Pol et al., 2014; Figure 3C). For WS, M and W the maximum individual inhibitor accumulation were below 4 g L −1 for HMF and furfural, which is important as up to 4 g L −1 of furfural or HMF has been shown to not inhibit microbial growth (Geddes et al., 2011;Mattam et al., 2016;Wang et al., 2016). Previous publications already mentioned have a range of inhibitors concentrations from 7.5 to 45.8 g kg −1  3.62 to 5.65 g L −1 (Negro et al., 2014)    used activated carbon with cation and anion resin before xylose crystallization. Optimization of both xylose release and reduced inhibitor production is required to ensure that scale-up with reduced process steps could be achieved.

Taguchi DoE L 16 Analysis
The influence of the selected factors on the pre-treatment process is given in Table 3 along with the average effect of a factor at each level.   the variation in the levels, the more the factor is influencing the process. Through Taguchi DoE analysis, each feedstock was shown to be influenced to a differing degree (percent influence) by factors affecting xylose release ( Table 4). The factor exerting the greatest influence on the process can be identified in the last column of the ANOVA table ( Table 4). The indication from ANOVA is that each feedstock was influenced to a different extent. For WS, acid concentration had a much higher influence (59.8%) as an individual factor on xylose release while pressure only had a 22% influence. This is in contrast to the other three feedstock where pressure had an effect of 66.2, 70.2, and 77.1% respectively as an individual factor. The individual factors of time, acid and substrate load had minor influence on the xylose release ( Table 4).
Interaction of the factors for each feedstock is shown in Table 5. The interaction between factors illustrates how in combination the factors involved can influence the pre-treatment process through either synergistic or antagonistic effects. Factor interactions can therefore be considered to have a similar effect to those of individual factors, thereby contributing positively or negatively on the net effect of the factors in the pre-treatment process. The magnitude of these interactions are calculated as a severity index (SI) percentage, which should not be confused with the severity factor used to describe pre-treatment processes. Time x substrate load and acid x substrate load were found to be most influential exerting the highest effect on the pre-treatment process for each feedstock. Pressure had high influence as an individual factor was not influential in combination with any other factor.
The quantity of lignin acting as a barrier to lignocellulose hydrolysis and the ability to disrupt the lignin structure are key factors recognized for reducing biomass recalcitrance (Behera et al., 2014;Rabemanolontsoa and Saka, 2016). Compositional analysis of each feedstock found that wheat straw had the highest levels of carbohydrates and the lowest level of lignin ( Table 2). The low lignin level was likely to be a key reason why acid exerted the highest influence on xylose release from WS (Davison et al., 2006), whereas pressure played an important role in the case for CS, M and W. The degradation of the lignin matrix is also dependent upon the lignin monomeric composition with previous studies indicating that hardwoods (e.g., willow) have a higher coniferyl and sinapyl residue ratio, whereas monocot's have a greater level of p-coumaryl residue (Davison et al., 2006;Cavka, 2013;Gu, 2013). The additional pressure being applied to CS, M and W may enhance the catalytic effect by ensuring deeper pore penetration of each feedstock, however the resulting increase in temperature can also lead to sugar degradation. There are few previous publications exclusively discussing xylose release from steam explosion for subsequent bioconversion to renewable commodities (Tomás-Pejó et al., 2008;Misra et al., 2013;Wang et al., 2013Wang et al., , 2015Fan et al., 2014;Zhang et al., 2014), with all but Tomás-Pejó et al. (2008) using sulphuric acid for acid impregnation. Tomás-Pejó et al. (2008) only used water as the catalyst for SE-PT for subsequent enzymatic hydrolysis. Studies by Wang et al. (2013) most closely matched the current work in scale using 1 kg dried material and a 10 L batch reactor. These authors found soaking time of 6 h, 0.3% sulphuric acid concentration at a liquid to corncob ratio of 3:1 with SE-PT conditions of 200 • C and 5 min residence time gave the optimal conditions for xylose release, however this achieved only 23.2% xylose mass recovery. Furthermore, this material required detoxification prior to bioconversion due to the level of inhibitors released through SE-PT.
Analysis of variance revealed that all four factors had a significant effect on the soluble xylose released from each feedstock (P < 0.01). ANOVA was used to determine which levels of each factor in combination will yield the highest xylose concentration and predict the xylose concentration achieved  from each of these levels. All four ANOVA analyses determined that the optimum conditions were 12 bar(g), 12 min, 1.2% acid and 0.4 kg substrate ( Table 6A).
The data presented in the L 16 data set show that 1.2% w/v phosphoric acid was effective in releasing maximum xylose from all feedstock and was specified in the model (Table 6). However, compared to the L 16 data the model predicted 12 min residence time with a severity factor of 3.8 as optimum. The L 16 used 6, 3, 3, and 3 min residence times for maximum xylose release with associated severity factors of 3.5, 3.2, 3.2, and 3.2 for W, CS, WS and MG respectively. Some of the observed variation in xylose release would appear to be due to the inherent variation of steam explosion at pilot scale, despite attempts to minimize this by warming up the reaction vessel prior to and cleaning in between each treatment. The error bars displayed in Figure 2 highlight the variation (standard deviation) between duplicate explosions, suggesting that day to day variations of steam explosion at this scale may be the result.

Overall Evaluation Criteria Modeling
Knowing that inhibitors are present in the hydrolysate an alternative model was assessed using a combination of maximum xylose release and minimum total HMF and furfural production, an overall evaluation criteria (OEC) analysis using Qualitek 4 was performed. This analysis with a 75% weighting to maximum xylose and 25% weighting to minimum total inhibitor provided a further set of optimum explosion parameters ( The OEC analysis (Table 8) indicated that the monocots; M, WS, and CS all required the same optimum conditions as for maximum xylose release (Table 6) whilst willow, an arboreal dicot, would require lower residence time to reduce the amount of furfural produced through the SE-PT. This was surprising considering the fact that W produced the lowest levels of inhibitor observed in the L16 data (Figure 3), but equally produced the lowest levels of xylose, Figure 2. The OEC model parameters were used in SE-PT of all feedstock to assess the models accuracy.

Taguchi Validation Experiments
The concentration of monomeric carbohydrates, released into the hydrolysate during steam explosion, revealed that the predicted model parameters gave lower values of all carbohydrates compared to the L 16 best performing SE-PT (Table 8 and Figure 2). A maximum xylose release (up to 94%) achieved within L 16 experimental conditions outperforms that observed in the OEC model conditions, with considerable losses apparent from each feedstock, up to 50% lower xylose release than L 16 . The indication here is that the increase in severity factor from the L 16 to the model has resulted in a decrease in xylose released. Table 6 illustrates the increase in severity factor from L 16 to model with the lower severity experiments yielding greater xylose levels than higher severities. Previous authors have noted that increased severity factors led to reduction in carbohydrate concentrations and greater inhibitor levels Iroba et al., 2014;Sharma et al., 2015). It may be unsurprising that the model prediction of a 4-fold increase in residence time for CS, M and WS resulted in lowering xylose recovery as the increased length of exposure of labile pentose monomers to acid at high temperature is known to result in degradation to furfural and subsequent decomposition to formic acid (Larsson et al., 1999). Indeed taking together the xylose released and inhibitor formed from WS, 65% of xylose was recovered and almost 300% more inhibitor was produced when compared to run 13 of the L 16 experiment (Figures 2-4). In isolation, time was shown to exert a minor influence on xylose recovery from all feedstock (Table 4). However, the severity index of interacting factors revealed that time x substrate weight applied an influence of between 80 and 96.4% on the generation of xylose for each feedstock (Table 5) demonstrating a strong interaction with biomass load. Zhao et al. (2013) also noticed an interactive influence in assessment of sewage sludge for solid fuel production utilizing an L 16 orthogonal array with steam explosion. The orthogonal array could not account for the interaction of factors; pressure and residence time, and the fluctuation in moisture content for a 15 s change in residence time was likely due to the interaction of these factors. In light of these findings the OEC model was modified by reducing the residence time to 6 min and increasing the biomass loading to 1.5 kg of WS, resulting in 90% recovery of the original xylose in the hydrolysate.
The lower xylose release from all feedstock is attributed to the increased residence time offered by the OEC model and the lower levels of inhibitors in CS and M suggest that the increase in residence time has led to further decomposition of furfural and HMF to levulinic acid and/or formic acid and that this decomposition occurred in a feedstock dependent manner (Mussatto and Roberto, 2004;van der Pol et al., 2014). The inhibitor data for wheat straw supported those of xylose release in that the SE-PT severity was too harsh resulting in degradation into inhibitory compounds. Specifically this appears to have resulted from the increased residence time and reduction in biomass loading as dictated by the Taguchi OEC model. The Taguchi methodology has enabled 256 experiments, if performed through full factorial assessment, to be reduced down to 16 experiments. Using this methodology each of the four feedstock gave high xylose yields (>75% total xylose), all of which successfully had preliminary single flask assessments for bioconversion to produce xylitol ( Figure 5). However, for process optimization the wide range of levels on each factor resulted in a broad range of xylose release ∼10 mg g −1 DM to ≥197 mg g −1 DM for all monocots, Figure 2. The xylose release from willow was substantial ∼75% (Figure 2) however, the amount in the initial material was much lower at 15.7% ( Table 2) than found in monocots, consequently further assessment of willow was discounted. For the three monocots the furfural and HMF produced also showed considerable variation (Figure 3) from ∼5 mg g −1 DM to ∼110 mg g −1 DM. The extensive range of both xylose and inhibitor levels suggest that the factor levels were too broad affecting the reliability of the Taguchi OEC model.
Taguchi DoE requires a systematic approach to understanding: measuring inputs and outputs, identifying and understanding the noise factors within the system, choosing appropriate control factors and levels, and selection of a suitable design (Rao et al., 2008). Options to improve the current approach would be to reduce the range of levels for each factor and to narrow the parameter range, increase the number of replicates to decrease the variation within the experiments, or to conduct a fine tuning Taguchi methodology model. Previous research has indicated that utilizing an L 16 to do an initial optimization of methane production can be further improved by fine tuning these results via an L 4 orthogonal array (Adu-Gyamfi et al., 2012). Teng and Xu (2008) utilized Taguchi for initial optimization with RSM used subsequently to further enhance lipase yield. Had the current investigation been carried out via RSM many more experiments would have been required, for example 4 factors at 3 levels requires 27 experiments, an increase of 50% on Taguchi DoE. Both of these examples indicate the importance that Taguchi methodology can have for enhancing production whilst reducing experiments required.  In future a L 4 Taguchi orthogonal array could be performed to finalize maximal xylose release from WS whilst minimizing formation of fermentation inhibitors. However, the increase in yield is anticipated to be modest compared to that presented here.
Based upon these data and subsequent modification of the OEC model an optimized scalable process was developed that enabled 90% xylose recovery from WS using the parameters of 1.5 kg DM, 1.2% w/v H 3 PO 4 , 12 bar(g) and a residence time of 6 min. All four feedstock generated hydrolysates indicated bioconversion of xylose to xylitol (data not shown) however to assess a pilot scale bioconversion system a single feedstock was chosen; CS was discounted due to the high inhibitor loadings seen in Figure 3. While Miscanthus performed well in terms of xylose release and inhibitor production, the low availability of this feedstock globally precluded use in subsequent studies. Table 6 illustrates that the monocotyledonous material all required the same SE-PT parameters to yield the highest concentration of xylose, indicating that a biorefinery processing plant could equally utilize any of these biomass for SE-PT toward Frontiers in Energy Research | www.frontiersin.org xylitol bioconversion. These results enable a broader spectrum of locals for a biorefinery to be situated for supply of enough biomass for production. SE-PT was performed on 10's of kg of wheat straw to produce a total yield of ∼1,000 L of hydrolysate for scale-up assessment of a bioconversion system for xylitol production up to 170 L fermentation.

CONCLUSION
This work helps fill in the knowledge void in the production of xylose following pre-treatment of lignocellulosic material at pilot scale. Agricultural residues, wheat straw, corn stover, and dedicated energy crops, Miscanthus, willow, could all be used as biomass resources for xylose release for xylitol production. Using a steam explosion pre-treatment of 12 bar(g), 3 min, 1.2% phosphoric acid, and 500 g substrate attained a yield of up to 94% release of initial xylose into hydrolysates. The use of Taguchi method DoE reduced the number of explosions from 256 per feedstock to 16 and gave a high yielding set of parameters. Steam explosion pre-treatment of wheat straw for xylose release gave high yields with low levels of fermentation inhibitors and subsequent fermentation assessment of the hydrolysates has proved these are non-toxic and allow for xylitol generation (Figure 5). However, the OEC Taguchi model predicted optimum conditions that were not confirmed by experimental data, yielding lower xylose levels than those seen via L 16 experiments. The utilization of Taguchi methodology DoE is based upon knowledge of the system being investigated. Methods to improve model outputs include narrow levels for factors via in-depth knowledge or preliminary investigation and use of multiple orthogonal arrays to initially limit the levels and factors before fine tuning the model with a tapered orthogonal array and reduced levels.
The results presented indicated that the Taguchi orthogonal experimental design did reduce experimental numbers, costs and time and that modification of the OEC model produced a scalable SE-PT for recovering 90% of the original xylose in a fermentable form.

AUTHOR CONTRIBUTIONS
DB, JG, SR, and AW conceived the project. SR and DB designed the experiments. DW performed the experiments. DW, SR, and DB performed analysis and data interpretation. DW and AS wrote the manuscript with all authors contributing to manuscript editing and revisions.