Montclair State University Digital Montclair State University Digital Commons Commons Assessment of Switchgrass-Based Bioenergy Supply Using GIS- Assessment of Switchgrass-Based Bioenergy Supply Using GIS-Based Fuzzy Logic and Network Optimization in Missouri (U.S.A.) Based Fuzzy Logic and Network Optimization in Missouri (U.S.A.)

Part of the Biogeochemistry Commons, Environmental Education Commons, Environmental Health and Protection Commons, Environmental Indicators and Impact Assessment Commons, Environmental Monitoring Commons, Geophysics and Seismology Commons, Natural Resource Economics Commons, Natural Resources and Conservation Commons, Natural Resources Management and Policy Commons, Oil, Gas, and Energy Commons, Other Environmental Sciences Commons, Paleontology Commons, Sustainability Commons, and the Water Resource Management Commons Abstract: Bioenergy has been globally recognized as one of the sustainable alternatives to fossil fuels. An assured supply of biomass feedstocks is a crucial bottleneck for the bioenergy industry emanating from uncertainties in land-use changes and future prices. Analytical approaches deriving from geographical information systems (GIS)-based analysis, mathematical modeling, optimization analyses, and empirical techniques have been widely used to evaluate the potential for bioenergy feedstock. In this study, we propose a three-phase methodology integrating fuzzy logic, network optimization, and ecosystem services assessment to estimate potential bioenergy supply. The fuzzy logic analysis uses multiple spatial criteria to identify suitable biomass cultivating regions. We extract spatial information based on favorable conditions and potential constraints, such as developed urban areas and croplands. Further, the network analysis uses the road network and existing bioreﬁneries to evaluate feedstock production locations. Our analysis extends previous studies by incorporating biodiversity and ecologically sensitive areas into the analysis, as well as incorporating ecosystem service beneﬁts as an additional driver for adoption, ensuring that biomass cultivation will minimize the negative consequences of large-scale land-use change. We apply the concept of assessing the potential for switchgrass-based bioenergy in Missouri to the proposed methodology.


Introduction
Bioenergy has been recognized as one of the sustainable alternatives to fossil fuels, and the quantification of bioenergy potential is an essential field of research. The design of the biomass supply chain is a major stumbling block for bioenergy and includes three spatially interlinked elements, including (i) resource assessment, (ii) logistic planning, and (iii) biorefinery design [1]. Perennial biomass crops, such as switchgrass (Panicum virgatum L.), face adoption challenges owing to the considerable time required to establish the plant and uncertainties associated with land-use change and price fluctuations [2][3][4]. However, an inadequate volume of potential biomass feedstocks also presents a barrier; Richard et al. suggested that ethanol production at the commercial level is challenged by the low density of biomass feedstock across the landscape [5,6]. This situation has led investments in the bioenergy industry to be described as a chicken-and-egg problem, wherein investment in feedstock processing and fuel production is not forthcoming due to uncertainty in feedstock supply, while farmers values for switchgrass land cover. Our analysis extends prior studies by incorporating biodiversity and ecologically sensitive areas into the fuzzy logic analysis, ensuring that biomass cultivation will minimize negative consequences of large-scale land-use change and how these changes impact ecosystem services relating to carbon sequestration, annual water yield, and soil erosion. We apply the following conceptual framework (Figure 1) to assess the potential for switchgrass-based bioenergy in Missouri.
Energies 2020, 13, x FOR PEER REVIEW 3 of 18 from bioenergy crop transitions, we integrated an ecosystem service assessment based on the known values for switchgrass land cover. Our analysis extends prior studies by incorporating biodiversity and ecologically sensitive areas into the fuzzy logic analysis, ensuring that biomass cultivation will minimize negative consequences of large-scale land-use change and how these changes impact ecosystem services relating to carbon sequestration, annual water yield, and soil erosion. We apply the following conceptual framework (Figure 1) to assess the potential for switchgrass-based bioenergy in Missouri.

Methods
The current study assesses bioenergy potential for biomass cultivation and spatial optimization of existing bioenergy facilities. Figure 1 illustrates the conceptual model structures used in this research. In the first step of this method, we applied multiple-criteria spatial factors (including environmental and socioeconomic) to locate suitable regions for growing energy crops. These spatial factors were described as geospatial rasters (30 m resolution) and transformed into binary images by reclassifying cells within the favored areas by "1" and cells outside the area by "0". These reclassified layers were then used as input data to a fuzzy overlay model which is expressed by the following equation [14]: where , was the value of the i-th cell in the final overlay maps; , was the Boolean value (0,1) assigned for the i-th cell in the j-th input layers; n is the number of factors considered in the study.
We reclassified the collected data of both preferred and constraint factors to create discrete variables based on the selected growing conditions. The reclassified layer was used as an input for the overlay analysis by combining favorable conditions for all variables using the "AND" combination method [23]. As such, if a location was identified as suitable under all data layers, then it would be classified as suitable. Consequently, if any aspects/variables were unsuitable, the entire cell was determined to be unsuitable. Therefore, our assessment represents a conservative estimate of suitable growing areas on which a bioenergy crop could be produced. We used the land use/cover

Methods
The current study assesses bioenergy potential for biomass cultivation and spatial optimization of existing bioenergy facilities. Figure 1 illustrates the conceptual model structures used in this research. In the first step of this method, we applied multiple-criteria spatial factors (including environmental and socioeconomic) to locate suitable regions for growing energy crops. These spatial factors were described as geospatial rasters (30 m resolution) and transformed into binary images by reclassifying cells within the favored areas by "1" and cells outside the area by "0". These reclassified layers were then used as input data to a fuzzy overlay model which is expressed by the following equation [14]: where C E,i was the value of the i-th cell in the final overlay maps; C i,j was the Boolean value (0,1) assigned for the i-th cell in the j-th input layers; n is the number of factors considered in the study. We reclassified the collected data of both preferred and constraint factors to create discrete variables based on the selected growing conditions. The reclassified layer was used as an input for the overlay analysis by combining favorable conditions for all variables using the "AND" combination method [23]. As such, if a location was identified as suitable under all data layers, then it would be classified as suitable. Consequently, if any aspects/variables were unsuitable, the entire cell was determined to be unsuitable. Therefore, our assessment represents a conservative estimate of suitable growing areas on which a bioenergy crop could be produced. We used the land use/cover raster to Energies 2020, 13, 4516 4 of 18 build scenarios by changing land use/cover attributes. For example, scenarios could be based on using marginal land for growing the crop in the research region, or land use/cover scenarios can be modified based on available conditions.
In the second phase, the study utilizes network location and sensitivity analyses to identify transportation frameworks to optimize logistics and supply chain aspects of biomass delivery. The network analysis is based on Dijkstra's algorithm to find the single-source, shortest-path connecting two nodes (a source and a destination) in a weighted graph G = (V, E), where V is a set of vertices or elements, including nodes, junctions, or intersections in a road network, and E is a set of directed edges or road segments that connect ordered pairs of vertices in V. The weight or cost of transport can be travel time and/or travel distance. The algorithm builds a set S of vertices, S V, which includes the final smallest weight w from the determined source s. Noting that w is the weight of an ordered pair (u, v) E, the algorithm repeatedly finds the vertex u V-S, adds u to S, and relaxes all edges different than u [24]. In this methodology, by placing the existing biorefinery at the center of assigned service areas as determined sources, s, we identify the potential biomass supply region for each refinery based on minimum distance computations. The outcome of the network location analyses is a set S of potential cultivation areas that connect to the biorefineries by the shortest paths. The ArcGIS Network Analyst extension is executed to achieve this step.
To estimate the feedstock area of existing bioethanol refineries, we used the Service-Area solver to create concentric serviced regions that encompass all open streets within a specified impedance. To estimate optimized transportation paths used to calculate the costs of transportation of biomass delivery, we applied Closet-Facility solver in the ArcGIS platform, which minimizes the distance from individual incidents and existing local biorefineries to their potential biomass sources. By doing this, we determined viability at each of these plants, though resources may not be allocated evenly, as many of these plants may not have sufficient supply depending on the scenario. We assumed that trucks depart from biorefineries to collect feedstock from potential areas identified through fuzzy logic analysis. We used two break conditions: break time and break distance for a professional truck driver [15]. This option allowed areas that do not have direct road access to be considered in the analysis. In both Service-Area and Closest-Facilities solvers, the impedances were travel time and travel distance. The Facilities-To-Find was set with the existing local biorefineries, while the Location-Type for both potential sites and facilities was established with an either-side-of-vehicle approaching type, to prevent the vehicles from having to approach from only one side of the road. The Search-Tolerance for street map features was 1.6 km (1-mile) to make sure that all the potential elements were captured. We assigned the Global-Turns to generate turn restrictions; the study assumed that all vehicles are allowed to perform U-turns only at intersections and dead-ends and that all potential sites and facilities could be approached from either side of the vehicle [15,22] Biomass transportation cost ($. Mg −1 ) is often calculated as a sum of total time cost and total distance transportation cost. While time cost accounts for the capital cost of the truck and labor cost, the distance transportation cost accounts for a fully loaded truck traveling to the facility and for empty trucks on the backhaul [9]. In this study, we calculated the transportation cost using the equation below (Equation (2) [9,25]).
where L, logistic cost, is the cost for loading and unloading biomass from a purpose-built bale truck, T, travel time (minutes) is time for a truck collecting biomass to travel both to and from the facility. This value is an output of Closest-Facility solver, while wage cost, W, is the average regional wage of a truck driver ($ per minute), fuel cost, F, is the average regional cost for fuel consumed during biomass transportation. Fixed cost, FC, is an average cost that represents vehicle maintenance, replacement cost, and other upkeep costs. Finally, we performed a sensitivity analysis by applying Monte Carlo simulation to identify possible outcomes and probabilities the transport cost will occur based on uncertainties of cost parameters.

Case Study
For the case study of growing switchgrass (Panicum virgatum L.) in Missouri, we located areas that are most suitable for switchgrass based on ideal growing conditions including soil type, annual rainfall, seasonal temperature, and the slope of the terrain. Many studies agree that switchgrass is considered as a high potential bioenergy feedstock as it can produce ample biomass, has sizable ethanol conversion potential, is native to the U.S., and can be grown on low-quality lands with limited use of fertilizers [26,27]. Switchgrass provides a range of ecosystem services and has the potential to be an excellent biomass feedstock. Mitchell et al., confirmed that growing switchgrass on marginal soil returned an ethanol yield equal to or higher than that of no-till corn.
The state of Missouri, located in the Midwestern United States, represents a varied landscape that consists of marginal and degraded lands, as well as lands that lie in floodplains, which researchers and policymakers have emphasized as necessary for bioenergy feedstock cultivation. Furthermore, switchgrass has a relatively high yield potential throughout the state of Missouri [28]. We assumed that new cellulosic-based biorefineries would be co-located with the existing ones in order to benefit from existing infrastructure. Six biorefineries are operating in the study region, whose details on location, feedstock, biofuel platform, and biofuel generation capacity can be found in Table 1 [29]. Using the biorefineries' locations, we created a 128-km (80-mile) buffer, which represents the longest one-way distance that a truck with a long-log trailer can haul biomass feedstocks effectively [25,30,31]. The buffered polygons were merged to generate one unified preliminary area. The final area of interest includes Missouri and neighboring states, including Kansas, Nebraska, Iowa, and Illinois, and covers a total area of 130,612 km 2 . We excluded wetlands, conservation regions, the US reservation areas, and developed areas (including street network, urban sections, and pipeline network) from the analysis, as they would be unsuitable for switchgrass cultivation. Our assumptions regarding favorable growing conditions were informed based on published literature. Table 2 lists the details, including values/ranges for specific variables, data types, and references for our assumptions.
The data used in this study were selected to represent the environmental conditions that determine feedstock suitability based on availability and economic feasibility of transportation. Spatial soil properties, land use types, road networks, temperature, rainfall, and other relevant inputs were collected in vector and raster format from various sources including geospatial data from the US Department of Agriculture (USDA), US Geological Survey (USGS), and ESRI ® Data & Maps. The maximum temperature data for 30 years were extracted from the National Oceanic and Atmospheric Administrator (NOAA) 1981-2010 Climate Normals database. Mean precipitation values were collected from the USDA Natural Resources Conservation Service (NRCS) from 1981-2010 [32]. The land use/cover data were extracted from the National Land Cover Database 2011 [33]. The spatial resolution was maintained at a 30 m resolution for all raster layers during suitable region assessment analyses.
To analyze the potential transportation of switchgrass to biorefineries, we used the 2005 detailed U.S. street map layer to plan routes from refineries to potential fields [34]. While it would be ideal to have a more recent road map, the required richness of data for our method necessitated using an older, but a complete road map.
Terrain slope has essential considerations when evaluating a site for potential agriculture due to difficulties associated with operating farming/harvesting equipment. The slope, in percentage, was represented in a raster using the slope spatial analysis tool in ArcGIS version 10.6. The input data for this tool used a digital elevation model (DEM) raster extracted from National Elevation Data 30m (NED_30m) compiled from USDA/NRCS-National Geospatial Center of Excellence [35].
We collected detailed soil texture data at the county level for the study region. The soil attribute database was downloaded from the Geospatial Data Gateway, USDA [36]. We used Soil Data Viewer 6.2 (SDV 6.2), an ArcGIS compatible toolkit supplied by the Natural Resources Conservation Service, to extract the necessary soil texture shapefiles to categorize data into broad soil-texture types and make our analysis tractable [37]. SDV 6.2 allows one to access soil interpretations and soil properties that are contained within the complex soil database. The SDV 6.2 generated spatial attributes of soil types and displays results as spatial maps. Since all spatial soil texture data at the county level were collected, we merged them and dissolved to create one soil layer. Soils were classified into 11 categories, including loamy sand, silt loam, silty clay loam, clay loam, fine sandy loam, loam, very gravelly silt loam, very gravelly silty clay loam, gravelly silt loam, moderately decomposed plant material, and silty clay [27,[37][38][39][40][41].
Falling temperatures can cause lower growth rates and ultimately lower yields. Therefore, we constructed a temperature raster that represented temperatures only during the growing season (May-September). Triangulated irregular network (TIN) interpolation is one of the most popular methods for representing digital models, which is used to estimate missing values of any location-dependent attributes such as precipitation, temperature, or elevation. We used the TIN interpolation tool in QGIS to generate the spatial temperature using monthly data from 150 adjacent national meteorological stations, resulting in a spatial temperature raster [42]. This raster was projected to WGS 1984/UTM Zone 15S and used as one of the inputs for the fuzzy overlay analysis. Critical species habitats, extracted from U.S. Fish & Wildlife Service, ECOS [43], were excluded from the potential area to mitigate unintended consequences. In the case study, based on the land use/cover classification, we investigated scenario (scenario 1) using marginal lands, which included areas of fallow/idle cropland, barren land, and switchgrass, as classified by Homer et al. [33]. Scenario 2 used marginal lands and non-food cropland for growing switchgrass. The non-food crop was defined as industrial crops, including grass, pasture, alfalfa, hay, and sod/grass seed. Table 3 describes land-use-land-cover (LULC) changes in different scenarios where fallow/idle cropland, alfalfa, grass/pasture, and sod/grass seed are theoretically converted to switchgrass. For applying network analyses in the case study, we used a detailed street map, which was extracted from ESRI_DM (2005), clipped to the preliminary potential area, and projected to the WGS 1984/UTM Zone. We used an eleven-hour break time and a break distance of 240 km (150 miles) as requested by the Federal Motor Carrier Safety Administration (FMCSA) for a truck carrying agricultural products [31,53]. We also set a 1.6-km (1-mile) search tolerance, which snapped the facilities to the closest street map features, allowing areas that do not have direct road access to be still considered in the analysis. While the biorefinery location is used as each facility point, a centroid pixel located at the center of the polygon, representing a suitable area, serves as a supply point (an incident) in the network [22]. Figure 2 describes the location and selected factors of the case study in Missouri.
Although we assumed that switchgrass could replace existing land cover where suitable areas are defined, these changes can cause negative impacts to their ecosystem services. To evaluate these potential impacts, we quantified carbon storage, annual water yield, and overland sedimentation to estimate potential feedbacks of the LULC changes in different scenarios using the Integrated Valuation of Ecosystem Services Tradeoffs (InVEST) 3.8.0 modeling software [54]. The annual water yield, a difference between the annual precipitation and real evapotranspiration, is estimated using the Budkyo curve method. The carbon storage model matches land cover to estimate carbon pools in soil, above-ground and below-ground biomass, and woody debris using a lookup table. The sediment delivery ratio model computes total sediment generated, sediment export, sediment retention using the universal soil loss equation (USLE) combined with a connectivity index [55]. We used LULC maps in three scenarios, including base, scenario 1, and scenario 2, as LULC inputs of models to compare between scenarios. Supplementary Materials describes the InVEST model's settings in detail.
To calculate transportation costs, we used the values presented in Table 4, and all values are adjusted to 2019 dollars. agricultural products [31,53]. We also set a 1.6-km (1-mile) search tolerance, which snapped the facilities to the closest street map features, allowing areas that do not have direct road access to be still considered in the analysis. While the biorefinery location is used as each facility point, a centroid pixel located at the center of the polygon, representing a suitable area, serves as a supply point (an incident) in the network [22]. Figure 2 describes the location and selected factors of the case study in Missouri.

Results
Switchgrass cultivation shows significant promise in the northern region of Missouri. If the entire proposed study areas in scenario 2 were cultivated, it would provide up to 8923 × 10 6 L of bioethanol at an average transport cost of $8.60 per Mg of switchgrass, or just over $0.16 per liter of ethanol. Also, switchgrass, due to its physiology, could reduce sedimentation by up to 2% and increase the carbon storage in the landscape by 3% by changing only a small portion of the land area. However, scenario 1 shows that at this spatial resolution, marginal land does not provide substantial feedstock for a switchgrass-based industry with current cropping patterns.

Potential Cultivating Regions of Switchgrass
In scenario 1, where only marginal land, fallow, and idle crops were selected, the suitable area accounted for approximately 4755.6 hectares, equivalent to a 74% increase in current switchgrass plantings. In scenario 2, wherein marginal land and non-food cropland were selected, the suitable area increased significantly to more than 1.7 million hectares, which accounts for upwards of 13.5% of the total study area (Figure 3). There were 116 and 3343 suitable locations found by the fuzzy overlay analysis in scenarios 1 and 2, respectively, that were suitable for switchgrass cultivation. The sizes of parcels were varied widely, ranging from 0.1 ha to 212.6 ha with a median of 9 ha in scenario 1, and 1 ha to 136,992 ha with a median of 150.8 ha in scenario 2. Figure 3 demonstrates the broad range of region sizes between scenarios 1 and 2. The dark blue pixels in Figure 3 each represent an aggregation of suitable 900 square meter parcels of land. Each operating biorefinery was shown with its production capacity, in millions of gallons per year (MMgy) and its feedstock type.

Switchgrass Feedstock Supply Areas
We used the Service-Area analyst to quantify accessibility from facility locations to the surrounding area. In our study, the service area estimated the ability of biorefineries to access their potential switchgrass suppliers by the local road network. Figure 4 presents the Service-Area analyst results. Although the network analyst was set with the maximum 11-h drive time, the results show that a haul truck departing from an existing biorefinery could reach any potential site in the study area within only 180 min, driving one-way. In one hour, the service areas covered about 24.3% and 25.9% of switchgrass suitable regions in scenario 1 and scenario 2, respectively. However, the coverage increased to 96.8% and 90.3% for an additional hour of driving time. Table 5 illustrates the size of service areas within 60 min, 120 min, and 180 min. Since we did not incorporate the cutoff costs at facilities such as loading/unloading time, the maximum travel time to collect biomass from a potential facility was approximately six hours. The dark blue pixels in Figure 3 each represent an aggregation of suitable 900 square meter parcels of land. Each operating biorefinery was shown with its production capacity, in millions of gallons per year (MMgy) and its feedstock type.

Switchgrass Feedstock Supply Areas
We used the Service-Area analyst to quantify accessibility from facility locations to the surrounding area. In our study, the service area estimated the ability of biorefineries to access their potential switchgrass suppliers by the local road network. Figure 4 presents the Service-Area analyst results. Although the network analyst was set with the maximum 11-h drive time, the results show that a haul truck departing from an existing biorefinery could reach any potential site in the study area within only 180 min, driving one-way. In one hour, the service areas covered about 24.3% and 25.9% of switchgrass suitable regions in scenario 1 and scenario 2, respectively.  The dark blue pixels in Figure 3 each represent an aggregation of suitable 900 square meter parcels of land. Each operating biorefinery was shown with its production capacity, in millions of gallons per year (MMgy) and its feedstock type.

Switchgrass Feedstock Supply Areas
We used the Service-Area analyst to quantify accessibility from facility locations to the surrounding area. In our study, the service area estimated the ability of biorefineries to access their potential switchgrass suppliers by the local road network. Figure 4 presents the Service-Area analyst results. Although the network analyst was set with the maximum 11-h drive time, the results show that a haul truck departing from an existing biorefinery could reach any potential site in the study area within only 180 min, driving one-way. In one hour, the service areas covered about 24.3% and 25.9% of switchgrass suitable regions in scenario 1 and scenario 2, respectively. However, the coverage increased to 96.8% and 90.3% for an additional hour of driving time. Table 5 illustrates the size of service areas within 60 min, 120 min, and 180 min. Since we did not incorporate the cutoff costs at facilities such as loading/unloading time, the maximum travel time to collect biomass from a potential facility was approximately six hours. However, the coverage increased to 96.8% and 90.3% for an additional hour of driving time. Table 5 illustrates the size of service areas within 60 min, 120 min, and 180 min. Since we did not incorporate the cutoff costs at facilities such as loading/unloading time, the maximum travel time to collect biomass from a potential facility was approximately six hours.  Figure 5 shows the shortest paths from the six biorefineries in Missouri to 116 potential sites in scenario 1 and 3343 potential sites in scenario 2, respectively. Based on the road system, the Closest-Facility apparatus determined the shortest path from each plant to each of the potential sites. Tables 6 and 7 present the results of the suitable selection for scenarios 1 and 2, respectively, where each biorefinery is presented with its potential growing site (facility), area (in hectares), and potential yield (in Mg).

Shortest Transportation Route
Energies 2020, 13, x FOR PEER REVIEW 10 of 18  Figure 5 shows the shortest paths from the six biorefineries in Missouri to 116 potential sites in scenario 1 and 3343 potential sites in scenario 2, respectively. Based on the road system, the Closest-Facility apparatus determined the shortest path from each plant to each of the potential sites. Tables 6 and 7 present the results of the suitable selection for scenarios 1 and 2, respectively, where each biorefinery is presented with its potential growing site (facility), area (in hectares), and potential yield (in Mg).     In scenario 2, the Lifeline Foods LLC and ICM Pilot Int Cell plant are estimated to gather the largest amount of biomass among biorefineries since there are 1387 potential sites potentially accessible from them. The accumulated area of these sites accounts for 51% of the total potential area. Following this plan are Golden Triangle Energy LLC and Poet Biorefining 2, which can access 755 and 605 sites, respectively.

Least Transportation Costs
Using the shortest path result ( Figure 5) and estimated regional transport cost per minute of travel, we estimated the transportation cost for collecting switchgrass as their feedstock. The cost in scenario 1 ranged from $3.62-$23.12 (Mg −1 ), and $4.0-$22.52 (Mg −1 ) in scenario 2. The wide range of transportation costs was caused by a broad range of distances carried by potential haul trucks. These distances were very dependent on the spatial distributions of potential sites. However, using a Monte Carlo simulation for 10,000 repetitions showed that 90% of transportation costs fell in the range of $3.70-$14.75 per Mg in scenario 1 and $2.75-$14.66 per Mg in scenario 2.
Further analysis of scenario 2 results showed that the least cost biorefinery point was Show Me Ethanol (3); this biorefinery also covered 6% of the potential switchgrass (see Figure 6). While it had the lowest transport cost, it covers very little of the study area effectively. Conversely, LF, LLC & ICM Pilot Int. Cell (2) covered 51% of the area while having the second-lowest transport cost. This makes it a prime candidate as a potential pilot plant for the region, with both low transportation costs and ideal land quality. In scenario 2, the Lifeline Foods LLC and ICM Pilot Int Cell plant are estimated to gather the largest amount of biomass among biorefineries since there are 1387 potential sites potentially accessible from them. The accumulated area of these sites accounts for 51% of the total potential area. Following this plan are Golden Triangle Energy LLC and Poet Biorefining 2, which can access 755 and 605 sites, respectively.

Least Transportation Costs
Using the shortest path result ( Figure 5) and estimated regional transport cost per minute of travel, we estimated the transportation cost for collecting switchgrass as their feedstock. The cost in scenario 1 ranged from $3.62-$23.12 (Mg −1 ), and $4.0-$22.52 (Mg −1 ) in scenario 2. The wide range of transportation costs was caused by a broad range of distances carried by potential haul trucks. These distances were very dependent on the spatial distributions of potential sites. However, using a Monte Carlo simulation for 10,000 repetitions showed that 90% of transportation costs fell in the range of $3.70-$14.75 per Mg in scenario 1 and $2.75-$14.66 per Mg in scenario 2.
Further analysis of scenario 2 results showed that the least cost biorefinery point was Show Me Ethanol (3); this biorefinery also covered 6% of the potential switchgrass (see Figure 6). While it had the lowest transport cost, it covers very little of the study area effectively. Conversely, LF, LLC & ICM Pilot Int. Cell (2) covered 51% of the area while having the second-lowest transport cost. This makes it a prime candidate as a potential pilot plant for the region, with both low transportation costs and ideal land quality.

Ecosystem Services
Broader impacts of these shifts are captured by ecosystem benefits given by land-use change to switchgrass. Namely, the InVEST models for Carbon Storage, Annual Water Yield, and Sediment Delivery Ratio were employed to quantify the feedback on carbon stocks, total water yield, surface erosion, and local streamflow caused by land-use conversion to switchgrass. Table 8 demonstrates changes in ecosystem services in each scenario in terms of carbon storage, water yield, sediment total, and sediment export. Table 8 also shows the results of scenarios 1 and 2 compared to the base case represented by the rate of change. Shifting to switchgrass showed a generally positive impact on both scenarios for the ecosystem. Although there was a small 0.004%, reduction in annual water yield in scenario 1, scenario 2 showed a strong positive effect on water yield. The effect is similar across the board, as scenario 1 simply identified too little available land for switchgrass growth for a substantial regional impact. However, scenario 2 identified far more available land, which, when cultivated, had a much more significant impact on ecosystem function. Carbon stocks could be increased by up to 3.325%, equivalent to 59,281,2967 Mg, by switching marginal land and non-food cropland to switchgrass due to the deep root system of switchgrass [51]. Soil erosion would decline substantially, with a reduction of about 2.14% in total sediment and a 5.06% decline in sediment export. Additionally, around 34.4 million m 3 of water would be allowed to infiltrate into the soil layer or drain to the stream network and account for a 1.36% increase in available water yield compared to the base case.
Subsidies have often been proposed for switchgrass cultivation [4,9], and ecosystem service benefits can be used to help justify these subsidies beyond end-use carbon offsets. Carbon stocks alone increase between 33.14-72.54 Mg per hectare for the land cover types studied. This increases the value of growing switchgrass as a crop over its utility as a bioenergy feedstock, which can be used

Ecosystem Services
Broader impacts of these shifts are captured by ecosystem benefits given by land-use change to switchgrass. Namely, the InVEST models for Carbon Storage, Annual Water Yield, and Sediment Delivery Ratio were employed to quantify the feedback on carbon stocks, total water yield, surface erosion, and local streamflow caused by land-use conversion to switchgrass. Table 8 demonstrates changes in ecosystem services in each scenario in terms of carbon storage, water yield, sediment total, and sediment export. Table 8 also shows the results of scenarios 1 and 2 compared to the base case represented by the rate of change. Shifting to switchgrass showed a generally positive impact on both scenarios for the ecosystem. Although there was a small 0.004%, reduction in annual water yield in scenario 1, scenario 2 showed a strong positive effect on water yield. The effect is similar across the board, as scenario 1 simply identified too little available land for switchgrass growth for a substantial regional impact. However, scenario 2 identified far more available land, which, when cultivated, had a much more significant impact on ecosystem function. Carbon stocks could be increased by up to 3.325%, equivalent to 59,281,2967 Mg, by switching marginal land and non-food cropland to switchgrass due to the deep root system of switchgrass [51]. Soil erosion would decline substantially, with a reduction of about 2.14% in total sediment and a 5.06% decline in sediment export. Additionally, around 34.4 million m 3 of water would be allowed to infiltrate into the soil layer or drain to the stream network and account for a 1.36% increase in available water yield compared to the base case.
Subsidies have often been proposed for switchgrass cultivation [4,9], and ecosystem service benefits can be used to help justify these subsidies beyond end-use carbon offsets. Carbon stocks alone increase between 33.14-72.54 Mg per hectare for the land cover types studied. This increases the value of growing switchgrass as a crop over its utility as a bioenergy feedstock, which can be used as a reason for additional subsidies based on an agreed carbon price divided out by years of cultivation.

Discussion
Expanding bioenergy feedstock supply while avoiding competition with food crops is a priority. This research estimates the potential for energy crop cultivation by considering consequences of land-use change and accounting for a range of ecosystem benefits. We used the combination of fuzzy logic and network optimization to evaluate potential switchgrass-based biomass supply and applied the InVEST ecosystem service approach to validate our selection for land cover transitions in Missouri. Furthermore, we count for the variations in yield, operation costs during growing and transporting biomass from potential sites to existing biorefineries within the study area.
Our results show that this methodology can be used to find the biorefineries with the highest potential for new feedstocks. In the case study, Show Me Ethanol LLC was the most efficient refinery in terms of biomass transport (see Tables 6 and 7). However, Lifeline Foods, LLC & ICM Pilot Int. Cell can potentially benefit from the highest land availability. Given that this plant is already collecting cellulosic materials to produce biofuel, it may be suitable to consider switchgrass-based biofuel production at this facility.
Switchgrass cultivation in Missouri shows some promise given the requisite infrastructure and demand. However, in scenario 1, where switchgrass cultivation was considered only on marginal land, we found that the potential regions could only cover about 4755 ha. The area of potential regions could increase to 1,766,454 ha in scenario 2, where non-food cropland such as grasslands, alfalfa, and pasture/hay were considered. The potential regions were spread mostly through the northwestern and the western regions of the research area; the unequal distribution can be explained by several factors, including the steep slope of the terrain and critical species' habitats in the southern and the southeastern sections of the study area preventing development. In scenario 2, using suitable regions to grow switchgrass could yield around 18 × 10 6 ÷ 39.6 × 10 6 Mg of biomass. As the 2016 Billion-Ton Report, 56 gal ton −1 (approximately 225.325 L Mg −1 ) can be used as the conversion rate for the transformation of biomass to biofuel [9]. Using this assumption, the biorefineries could produce about 4055 × 10 6 ÷ 8923 × 10 6 L of bioethanol depending on the level of fertilization and cultivation of the switchgrass (see Tables 2 and 3). Given the relatively small study area and the clustered nature of suitable land, the conversion plan in scenario 2 could be a potential approach to achieve large-scale production of bioenergy.
The ecosystem services provided by such an endeavor provide substantial benefit, including directly providing 59 million tons of carbon capture in addition to the carbon offset from ethanol use. The decrease in sedimentation and increased available water in a heavily agricultural area is significant as well.
The InVEST models applied in this research revealed the significant benefit switchgrass could have on ecosystem services, including annual water yield, carbon storage, and soil erosion. The outputs confirmed that converting suitable areas from original crops to switchgrass could yield positive feedback to the environment. In scenario 1, the impacts were indistinct because the land cover change is too small, less than 1% of the total area. In scenario 2, in contrast, the impacts on all studied services were sustainably remarkable. Besides yielding remarkable biomass, switchgrass is well-known for its ability to reduce soil erosion, increase soil organic carbon, and increase carbon sequestration [15]. In scenario 2, where almost 14% of the total area is employed growing switchgrass, we see a 3.33% increase in carbon storage, which accounts for around 59.3 × 10 6 Mg of carbon.
Similarly, alternative land-use can generate 0.28 Mg ha −1 less sedimentation than idle cropland and fallow land. Also, while no habitat risk assessment was conducted, there are studies on the benefit of growing switchgrass to species across the US [26]. For example, Werling et al. confirmed that low-input cultivation switchgrass enhances biodiversity in a study conducted in Michigan and Wisconsin [23]. Other studies, e.g., Bies and Hartman and Werling suggested that switchgrass provides multiple ecological niches for different species, including birds, insects, reptiles, and mammals by maintaining multi-dimension habitat structure [23,61,62]. The lack of adequate field tests for our region prevented analyzing the positive and negative externalities of biodiversity and increased fertilization on switchgrass cultivation.
While the extent of this research goes beyond what many have previously accomplished [12], that the study does not use yield as a normalizing factor for transportation cost, as costs are assumed to be incurred by optimized, fully loaded trucks. Additionally, transportation in our case only includes trucks as a potential vehicle of transport, as there is no existing infrastructure for other transportation vectors for most existing biorefineries in the study area. We also feel that analysis can be improved by using partial values in fuzzy logic rather than using a discrete, strict yes/no terminology. Furthermore, a full implementation might require additional considerations for the interactions between the various conditions that affect switchgrass growth, particularly the interplay between rainfall, soil type, and yield, where higher precipitation in conjunction with well-drained soils might produce favorable conditions [46]. Given the seasonal nature of the temperature data required for this study, we used recorded station data from the study region. Finally, as we only considered potentially expanding existing biorefineries for switchgrass-based biofuels, as opposed to building a new plant, our analysis can be extended by identifying the most suitable locations for cellulosic biorefineries.

Conclusions
Our study shows that existing biorefineries can relatively quickly increase their access to cellulosic feedstocks, focusing on using marginal land and non-food cropland to build towards greater bioenergy use ultimately. The analysis evaluated the accessibility of potential bioenergy based on locating suitable areas in which biomass can be produced and minimizing the transport cost from existing biorefineries to their potential suppliers. It identified potential biomass suitability zones based on favorable conditions as well as possible constraints such as developed urban areas, biodiversity, ecosystem services, and traditional food crop zones. Then, based on these cultivation zones, the study evaluated accessibility based on potential bioenergy crop yield and the least cost of transportation of locally existing biorefineries. This study can be used to provide data for developing a biomass supply system capable of operating cost-effectively at a commercial scale [9]. Based on the location of the biomass suitability zones, biorefineries can also consider partnering with local farmer associations or entering into contracts with farmers to ensure year-round biomass supply.
We used the methodology to assess the suitability of switchgrass-based bioenergy in Missouri and neighboring states as a case study. The case study results suggest that while that degraded or marginal land has been considered as a resource for growing switchgrass as an energy crop, the potential yield harvested from marginal land alone is not enough for proper development of the bioethanol industry in the state of Missouri. While switchgrass might not be the ideal crop in all non-food croplands, utilizing grasslands and pasture/hay lands for switchgrass cultivation might contribute significantly towards cellulosic-based biofuel expansion. Our results show new information for determining where this kind of analysis can be deployed and what biorefineries may be optimal by considering feedstock, transportation, and associated cost constraints. By following this method, one may find potential biorefinery locations in their research, which will lead to increased efficiency in feedstock transportation.
While our model addresses many factors that play a major role in the potential development of cellulosic-based biofuel production, future research should continue to expand on external factors to provide more details to inform effective decision-making. Our model focuses primarily on positive ecosystem services associated with bioenergy development, which could be expanded on to address the complex extensive feedback systems that industry development has on the environment. For example, the effects of monoculture expansion have been associated with adverse environmental consequences which alter the natural ecosystem even in systems utilizing native species [63]. Challenges associated with monoculture, such as reduced pollination and soil microbe biodiversity, should be explored for switchgrass development. Future work can incorporate additional transportation networks, ecosystem service externalities, and partial values in fuzzy analysis and variable interactions to add to the body of knowledge in this area and contribute towards policy planning, economic feasibility modeling, and best practices development and deployment.