An Integrated Assessment of GIS-MCA with Logistics Analysis for an Assessment of a Potential Decentralized Bioethanol Production System Using Distributed Agricultural Residues in Thailand

: Crop residues derived from post-harvesting process have been problematic due to an on-ﬁeld incineration, which caused air pollutants and greenhouse gas (GHG) emissions. An appropriate utilization of those biomasses can improve the environmental situation and provide a substitute for fossil fuels. Therefore, this study intends to analyze how left-over agricultural residues should be valorized in the decentralized bioethanol production configuration. With integrated techniques of geographical information system and multi-criteria analysis (GIS-MCA), we identify suitable locations for exhibiting decentralized sites matching the geographical backgrounds in each region. Under the precondition of a complete utilization of the agricultural residues, we found optimal installation numbers 71 units of decentralized production in total through suitability analysis. Conducting the location– allocation model, it is possible to determine production scales from the collectable spatially distributed biomass and transportation distances. Under the presumed conditions of installing 1 to 25 units, the logistics cost and total capital investment can reach USD 1.17–2.46 L − 1 and USD 1.17–6.93 L − 1 , respectively. The results from examining the technical potential and economic feasibility aspects are key to designing decentralized bioethanol production facilities and maximizing the utilization of agricultural residues in Thailand.


Introduction
Biofuel from energy crops has been adopted as a competent alternative for replacing fossil fuels such as gasoline and diesel. Bioethanol from sugar and starch crops plays a significant role with a global consumption of 106 billion L per year in 2019 [1]. On the scale of Thailand, commercialized bioethanol production has been promoted from sugarcanes and cassavas, aiming to achieve a 30% renewable energy share of the total final energy consumption [2] and to increase the production of bioethanol of 240% between 2010 and 2020 [3]. Nonetheless, conflicts around using bioethanol as sugar-and starch-based crops are still happening and need to be dealt with, such as the requirement of land use expansion and the depletion of soil nutrients from crop plantation [4,5].
Concerning environmental impacts from the conventional supply of energy crops, lignocellulosic biomass, also known as second-generation bioethanol, has received more attention as an alternative to food crop-based feedstocks, also called first-generation bioethanol. The advantages of this type of biomass include the abundant quantity and the reduced conflicts associated with using land use for food productions [6]. Utilizing cellulosic biomass in bioenergy production potentially contributes to reducing waste from the agricultural sector as well as providing additional income to farmers [7].
Bringing cellulosic biomass derived from agricultures into the bioenergy industry is not a novel technology since bioethanol from lignocellulosic feedstocks has reached a commercial level [8][9][10]. Nevertheless, as the production capacity increased, the higher exploitation of the disperse availability of residues can lead to a higher financial burden from the process of feedstock transportation, which influences the ethanol factory-gate price [11,12]. Due to the low energy density of agricultural residues, biofuel production requires the transportation of a high volume of input materials. Therefore, it is necessary to assess the feasibility of biomass logistics along with proposing optimal locations for biorefinery installation.
Previously, Thailand has been verified for the potential of regional lignocellulosic biomass derived from agricultural cultivation process in terms of quantity and bioethanol convertibility [13,14]. Moreover, advocating the utilization of this potential source of biomass may help to minimize the open-field burning of post-harvest residues in cultivation fields as a common treatment method in Thailand [15,16], which causes a serious issue of air pollution nowadays [17]. Although the post-harvest biomass in Thailand has played roles in private sectors and local communities, the utilization of residues has not been promoted in terms of liquid biofuel policy that encourages the commercialization of residual biomass on a broader scale [14].
The current bioethanol production chain in Thailand is considered following a centralized distribution. They have been implemented in facilities that are mainly located in the central and northeastern regions based on the determination of the proximity from conventional bioenergy crops, including cassava and sugarcane plantation [18]. Thus, in order to encourage the valorization of distributed agricultural residues, decentralized configuration is a proposed concept in this study.

Decentralization Designs for Biorefineries
Traditionally, biofuel production facilities or biorefineries rely on the centralized supply chain with large-scale production. However, it is becoming more challenging when the technology needs to be compatible to the volatility of biomass quantity [19]. On the other hand, decentralized facilities initiate the idea to exploiting localized resources to be utilized for on-site production [20]. According to the description by Fleischmann et al., centralization is defined as a scheme consisting of a number of refineries operating in the same concentrated areas [21]. On the contrary, decentralization refers to the structure of the same manufacturing configurations that are sited and operated simultaneously in several distributed areas [21].
However, traditionally, there has been a reluctance of decentralized production systems, mainly due to the economy of scale which is potentially achievable by larger production facilities. To pursue economic feasibility, large-scale production can be more favorable for the lower operational costs and capital investment [22,23]. For the conventional sugar-based ethanol production, small-scale bioethanol in Brazil can possibly be feasible at 1000 L day −1 as a micro-distillery [24]. However, regarding the current global status of lignocellulosic bioethanol production in 2015-2020, the commercial production capacity ranges from 40 to 200 mL year −1 [8,25,26]. In terms of decentralized refineries, there have not been specific thresholds to define the productivity of the decentralized cellulosic bioethanol facility.
Nevertheless, the decentralized production scheme could improve sustainable development. For instance, it is possible to encourage the collection of waste into different forms of energy production, and higher valued products [26,27]. An outstanding advantage of practicing localized biorefineries concerns the accessibility of dispersed biomass, resulting in a potential reduction in the logistics costs [11,28,29], reduced material and energy loss from transportation [29], as well as greenhouse gas reduction from transporting shorter distances [30].
In the context of Thailand, the development of crop residue management has been mainly motivated by the aim to restrict field burning incidents [31,32]. The utilization Sustainability 2022, 14,9885 3 of 24 of crop residues is commonly targeted on the post-harvest of rice and sugarcane, which has thus far been proposed in light of the use of electricity generation through biomass combustion or gasification technology [33,34]. Although residues derived from palms have been paid attention to, deploying palm residues for bioethanol production has not been approached in Thailand [35]. The concept of employing crop residues in lignocellulosic bioethanol has been theoretically contemplated to assess its potential [13,14]. Furthermore, in the previous study by Delivand et al., the logistics costs of rice residues have been analyzed via a collecting and transporting technical approach, resulting in a logistics cost of USD 19-20 for each ton of biomass −1 [33]. However, integrating geographical information system (GIS)-based analysis for decentralization configuration still represents a gap in the research which is believed to possibly enhance biomass assemble efficiency.

Bioenergy Supply Chain and Logistics Analysis
In the generic bioenergy supply chain from agricultural residues, the process comprises biomass collection after harvesting, preprocessing, and transporting from sources to the production facility [33,36]. Therefore, optimizing the biomass delivery as a part of the biorefinery operating costs is considered imperative in enhancing the techno-economic feasibility [36], as the logistics cost dominates around 25% of total production costs [37]. In some cases, the logistics cost was found to make up a major proportion of the total feedstock price, i.e., almost 20-30% [37,38].
With an attempt to reduce operating costs, there are several potential solutions to minimize logistics costs. One of the important solutions is to lessen the transportation costs associated with delivery distances. Although several distances are recommended under different conditions, road transportation for lignocellulosic biomass is recommended for short distances under 100 km [39]. In the cases without depot hubs, the approximate optimal distance for a biomass collection radius was reported to be within 80 km of the direct collection to the refinery [40].
Optimizing the logistics cost can also be approached from biorefinery site selection. The GIS application provides useful instruments for dealing with the scattered spatial allocation of biomass through network analysis to connect biomass sources with the refinery sites [41][42][43]. Several studies identified the optimal plant locations and production scale based on high density of potential resources through the GIS-based location-allocation model [20,44,45]. The site selection can be determined for the heterogeneous resources, as well as recommending optimal unit numbers of biorefinery facilities for minimizing delivery distances [46].

Goal and Scope of the Study
The main goal of this study is to propose the conceptual decentralization based on the notion to maximize the collection and utilization of spatially distributed agricultural residues in cellulosic bioethanol refineries. The present study aims to determine optimal locations for plant installation with units and production sizes concerning suitable geographical conditions in different areas of Thailand. The goal is also settled to investigate the effects of logistics costs and capital investment costs based on the proposed decentralized designs in each region as a contribution to techno-economic assessment. Lastly, decentralized models, as proposed in this study, are also contextualized by demonstrating an environmental approach. The amount of collectable crop residues is used to investigate how much atmospheric emissions can be abated from preventing open-field burning.
Thus far, the identification of biomass for biorefineries has been mainly focused on a single kind of biomass associated with only some specific locations, or with the purpose of combusting biomass for electric power generation [47,48]. Therefore, comparing to relocating the heterogeneous agricultural post-harvest from the national distribution would essentially become a novelty of the study.

Materials and Methods
The methodologies in the current study are constructed from GIS-based modeling and mathematic computation for logistic analysis. Figure 1 illustrates a flowchart indicating the input dataset, modelling procedures used in this study, and the expectable outputs. Spatial analysis and network analysis are conducted through the ArcGIS 10.8.1 ® software by ESRI (Environmental Systems Research Institute, Redlands, CA, USA). Input spatial information used in the modeling analysis was collected both in vector and raster format. For nonspatial inputs, data were collected from statistical reports and related publications, which are elaborated upon in the following sections. However, with respect to the calculation of logistics and investment costs based on the same baseline, the current work did not consider the modification process. of combusting biomass for electric power generation [47,48]. Therefore, comparing to relocating the heterogeneous agricultural post-harvest from the national distribution would essentially become a novelty of the study.

Materials and Methods
The methodologies in the current study are constructed from GIS-based modeling and mathematic computation for logistic analysis. Figure 1 illustrates a flowchart indicating the input dataset, modelling procedures used in this study, and the expectable outputs. Spatial analysis and network analysis are conducted through the ArcGIS 10.8.1 ® software by ESRI (Environmental Systems Research Institute, Redlands, CA, USA). Input spatial information used in the modeling analysis was collected both in vector and raster format. For non-spatial inputs, data were collected from statistical reports and related publications, which are elaborated upon in the following sections. However, with respect to the calculation of logistics and investment costs based on the same baseline, the current work did not consider the modification process.

Regional Case Studiues and Gricultural Residues Allocation
As previously clarified by Jusakulvijit et al., the most generated residues in Thailand could be obtained from crop plantation comprising sugarcane residues (SRs), cassava residues (CRs), rice residues (RRs), and palm residues (PRs), providing availability and den-

Regional Case Studiues and Gricultural Residues Allocation
As previously clarified by Jusakulvijit et al., the most generated residues in Thailand could be obtained from crop plantation comprising sugarcane residues (SRs), cassava residues (CRs), rice residues (RRs), and palm residues (PRs), providing availability and density distributed differently in each region of Thailand [13]. These types of biomass are selected as the targeted feedstocks in this study with a summary of the total production yields in Table 1. The residue productivity, as obtained from the previous study, refers to the conversion from crop productivities based on the national statistics database according to the Office of Agricultural Economics and residue-to-crop ratio derived from the literature [13,49]. Table 1. Characteristics of residual biomass from agricultural sector in Thailand and the potential bioethanol conversion (adapted with permission from Jusakulvijit et al. [13]). In this study, we divided the country level into 4 sub-boundaries on the basis of administrative boundaries and the regional characteristics of biomass availability according to previous schematic findings of crop residues [13]. The context of existing bioethanol facilities also differentiated the backgrounds of each targeted region. Case studies R1_South and R2_North aim to depict the construction of decentralization due to non-located facilities in this area. Meanwhile, the case of R3_Northeast accommodates the current infrastructures for 8 production units combined with new candidate plants. For the case R4_Central, biomass allocation and logistics analysis were performed by considering only currently operational facilities, with 17 production facilities in this region.

Regions
Based on the provincial scale of crop residue availability, residue distribution was created into 1 × 1 km 2 pixels of biomass distribution by incorporating a layer of economic crop zoning, a district scale of residue production, and fishnet layers from the fishnet tool of ESRI ArcGIS 10.8.1, following the procedure by Wochele et al. and Perpiñá et al. [50,51]. The economic crop zoning in the shapefile format was acquired from the Land Development Department [52]. After that, only one type of residue that contributed the maximum production was selected in order to represent the production per one grid cell. The precise scales of resource allocation could be obtained to be further used in the suitability analysis and the location-allocation model.

Geophical Indicators as Input Multi-Criteria
In the initial process of identifying suitable locations, geographical criteria comprising 'residue density', 'type of residue', 'land use type', 'proximity to road network', and 'slope' were accommodated in this study. The weights of important input attributes used in the suitability analysis were normalized from the Delphi-AHP process and applied in the input layers [53]. The scores of criteria prioritizations after the normalization from the publication are summarized in Table 2. Residue density [49] Preferable residue density depends on regional distribution. Due to uneven distribution, reclassification is created based on Jenks' natural breaks.
[54] 23 Types of residue [13] Preferential levels are determined from the previous study, classified from the highest values of residues-to-ethanol convertibility to the lowest. Cassava residues > palm residues > sugarcane residues > rice residues [13] 19 Land use type [55] (paddy fields, field crops, perennial, orchard, horticulture, pasture and farm house, rangeland, industrial land, other built-up land, integrated farm/diversified farms).

Suitability Analysis
In the current work, the analysis for determining suitable locations was conducted by using integrated tools combining GIS applications and multi-criteria decision analysis (MCA) [58]. The integrated methods of GIS-MCA, as a tool for suitability analysis, are broadly used in many studies related to location identification. To incorporate the opinions in decision-making methods, multi-criteria evaluation can account for the opinions from actual players [28,62]. Dealing with MCA in this study, the analytic hierarchy process (AHP) was adopted in the combined tools to consider the important levels of regarded parameters [44,45].
In order to propose specific locations for decentralized biorefinery development, suitability analysis was applied in this study. Firstly, the values of input spatial layers were reclassified into 5 levels since each parameter has a different range of values. Then, the value of ith pixel in the kth level of preference of criteria layers (C ik ) could be derived. Following that, the normalized weights of each parameter (w k ), as derived from Section 2.2, were applied to respective meshed polygons to determine the suitability index, as described in Equation (1) through the weighted overlay function in ArcGIS. The output map illustrates the suitability levels which indicate indexed scores representing preferential levels 1 to 5 (S i ).
Since suitable locations are suggested under the 1 × 1 km 2 resolutions, there are several suitable points that share the same district (Amphoe level: the second sub-divided administrative boundary of Thailand). In order to prevent the supply competition between adjacent bioethanol producers, we summarized multiple points to a single point to be used as a candidate production plant per district by dissolving the centroids among suitability level 5, which share the same district boundary through the functionality on ArcGIS.

Location-Allocation Model
The location-allocation model was applied to allocate input materials to production facilities within the catchment areas. From preselected suitable candidate locations, the location-allocation model enabled the selection of optimal plants (n) that could connect the shortest distances from the source point i to the candidate facility j (d ij ) via p-median problem solving [54]. Equation (2) describes the target of the computation based on the input vector layer of the road network, locations of the candidate facilities, and the residue weights. The impedance cut-off at 80 km was established as a threshold distance restricting the maximum radius from the facility to the resource supply. Equations (3)-(5) refer to the condition that the biomass from each source is limited to assign to only one facility in order to avoid overlapping between demand facilities. Then, biomass resource weights at point i (w i ) could be supplied to facility j only when the facility is still open in the respective area. Finally, when assigning the numbers of facilities to p, the total numbers of selected candidate facilities are denoted to be equal to p.
In this study, we run the location-allocation model for V = {1, 2, . . . , n}. The setup numbers (n) for R1_North, R2_South, R3_Northeast, and R4_Northeast were 25, 25, 45, and 17, respectively. By increasing the set-up numbers of installing units, it was possible to illustrate the coverage rates of residual biomasses. The allocated residues to the proposed facilities could then be used to determine the potential capacity from the supplied feedstocks.
where x ij is an indicator of the assigned residue. If it is equal to 1 when there is agricultural residue assigned to facility j and 0 if not, y ij indicates the candidate points when it is selected as facility point y ij = 1 or y ij = 0.

Collectable Residue Quantity and Bioethanol Production Capacity
From the location-allocation model, the quantity of available wet biomass in the catchment area could be obtained (AWB). However, to pursue the biomass quantity that could reflect realistic values, the usable biomass per biorefinery facility (UB) was calculated through Equation (8) by dividing it by n numbers of optimal plants. The AWB was applied with a usable biomass rate (UBR) and a collectable residue ratio (CRR), which were taken as 0.85 and 0.50, respectively [36,63]. This collectable ratio was also considered based on the functionality of soil erosion protection from left-over crop residues. Meanwhile, considering the quantities of residues with a bioenergy purpose, farm accessibility rates indicating the farmland factor (FF) and the transportation and storage loss (TSL) were also applied in the calculation [36]. For the moisture content (M), the dataset for each type of biomass was acquired from the previous reports in the literature [13]. Then, through Equations (8) and (9), the annual dry biomass production based on the average per facility (DUB) could be formulated.
Then, the theoretical bioethanol production capacity (Ecap) was estimated based on the usable dry biomass quantity (DUB) fed per biorefinery and the biomass to bioethanol convertibility (BEC ), as presented in Equation (10). The residue-to-ethanol conversion factors (BEC) have been retrieved in the previous findings by Jusakulvijit et al. [13]. The value of parameters used in the calculation can be referred in Supplementary Table S1 with input data sorted from the literature [36,63,64].

Logistics Cost Analysis
Logistics cost estimation conducted in this study comprised the costs from biomass handling, collecting, loading, and delivering from resources to biorefineries. On the other hand, we did not include the biomass pre-processing procedures and the in-plant processing cost at the biomass storage process. The costs associated with nutrient replacement from residue utilization in this study are also not included since the calculation already concerns the collection fraction for bioenergy purposes at 50% of the available biomass. The logistics costs were calculated separately through spreadsheet model development under each regional area using spatial analysis, referring to the analysis of Sultana et al. [65]. Input monetary data for cost calculation in this study were collected from the literature and adjusted with an inflation index to the year 2020. For the price regarding fuel consumption in biomass transportation, values were collected from the diesel selling price in 2020 and converted from THB to USD [66]. Regarding the parameters for the transportation cost calculation, variable values were collected from the literature, as summarized in Supplementary Table S1. (detailed input parameters see Refs. [67][68][69][70][71][72][73]).
In this study, the collecting methods were generalized under the same baling method for every type of biomass since the bale processing costs are lower than the pelletized biomass [66]. Equation (11) elaborates the formula of the biomass collecting and handling cost (HC), constructed from the processing costs of biomass shredding (HC shredding ), biomass raking (HC raking ), biomass baling (HC bale_collection ), biomass twining (HC twining ), biomass storage (HC storage ), farmer premium (HC f armer_prem ), biomass grinding to size 10 mm (HC grinder ), and tub-grinding operation (HC tub_grinding ). HC = (HC shredding + HC raking + HC bale collection + HC twining + HC storage + HC f armer_prem + HC grinder Next, the transportation cost regarding biomass delivery (DC) was carried out following Equations (12) and (13). The cost was estimated from the portion of truck delivery and the biomass loading-unloading cost from trucks. The delivery cost was formulated from multiplying d with the tortuosity factor (Tr) [39], the variable of the truck transport per weight per distance in a trip (Tv), and the delivery rounds calculated from usable biomass weight divided by the truck capacity (Truck cap ). For the biomass loading-unloading cost, the loading-unloading factor (BM load ) was applied in the equation.
Finally, the logistics cost per production facility (LC) and the unitary logistics cost (L uc ) per production volume were computed in the study, following Equations (14) and (15). Then, the logistics cost results per delivered feedstock were used for comparing the transportation costs at different scales of bioethanol production capacity in the later section.
where d is an average round trip distance from residue collection delivered to bioethanol refineries and N Line represents the number of connected routes derived from the locationallocation model.

Total Capital Investment (TCI) Estimation
TCI is one of the most important factors for determining the optimal scale when designing a decentralized scheme for comparing the production size with an estimated logistics cost [74]. Technically, to address the precise investment cost, it is necessary to break down the cost of associated factory elements. In this study, however, the cost estimation process was considered from the baseline non-linear scale factor (n), as mentioned in the technical report by Humbird et al. [75]. The exponent factor n = 0.6 was applied in none-linear regression, as described in Equation (16). The new TCI could be computed with an input base plant capacity (P 0 ) and base cost (C 0 ) of 230.58 mL year −1 and USD 594.28 million, respectively, which are values reported by Humbird et al. [75]. Thus, the expectable new cost (C) could be obtained depending on the new input plant capacity (P), as estimated in this study.

Emission Abatement from Open Field Burning
GHG emissions and atmospheric emissions in Thailand mainly occurred from sugarcane and rice cultivation due to the process of post-harvest biomass [15,16], whereas in general the post-harvest residues from palm plantation are being left on plantation areas without the incinerating process [76]. The open-field burning process emits GHGs, as we considered CO 2 , CH 4 , and N 2 O in this study, and health-threatening particulate matters (PMs) such as PM 10 and PM 2.5 [15,77].
Firstly, the total emission factor (EF total ) for each type of residue was formulated from the summation of emission factors of CO 2 , CH 4 , and N 2 O, as written in EF CO 2 , EF CH 4 , and EF N 2 O , respectively (in Equation (17)). For EF CH 4 and EF N 2 O , global warming potential (GWP) values of 25 and 298 were multiplied in order to convert the equivalent of kg CO 2 [78].
Then, in order to calculate the GHG abatement potential and air pollutants, as presented in Equation (18), the EF total was multiplied with a fraction of biomass subjected to open burn (FR), and the combustion factors (CR) which correspond to each kind of agricultural residues and the collectable dry residue quantity (W c ). The methods performed by Sornpoon et al. and Gadde et al. are referenced in the current study [77,79]. Meanwhile, the amount of PM reduction potential was estimated from multiplication with emission factors of PM 2.5 (EF PM 2.5 ) and PM 10 [80][81][82][83][84][85][86][87][88][89][90].

Agricultural Residues Allocation
High resolutions of biomass allocation were successfully illustrated, indicating available yields per pixel (kg km −2 ) of SRs, CRs, RRs, and PRs. The distribution of four types of agricultural residue was also classified in a gridded map on ArcGIS into different categorized production densities, as shown in Supplementary Figure S1. Different colors indicate a variety of biomass types and the maximum density of residue generation per square kilometer.
A summary of the distributed biomass quantifications is categorized in each region in Table 3. It can be observed that R3_Northeast produced the highest biomass availability rates, whereas the lowest available biomass was found in R1_North. In terms of biomass type categorization, the most generated residues were SRs, followed by RRs. The least generated biomass was PRs, which was mostly collectable in R2_South. In terms of cultivation areas, the largest cultivation farms were found in the R3_Northeast, occupying an area of 139,377 km 2 . The dominant residues were from sugarcane's post-harvest SRs. Following that, the second largest farm areas were found in the R1_North, with SRs representing a majority biomass allocation. A similar size was allocated to the total cultivation areas in R4_Central, i.e., the most available biomass, dominated by RRs. Finally, the R2_South constituted the smallest cultivation area with palm plantation and paddy fields. In total, for Thailand's potential farms as a whole nation, we identified 296,858 km 2 with collectable biomass in this study.

Site Selection from Suitability Analysis
An exemplary output map (see Figure 2) of suitable areas was classified into five levels, as illustrated in Figure 2b. The most suitable areas at level 5 were illustrated in red raster layers. We could identify the suitable level 5 into one specific candidate location per district boundary by dissolving the points of extensive layers, as presented in Figure 2c. As a result, R3_Northeast found the greatest numbers of suitable locations, resulting in 221 areas. R3_Northeastern was distributed as the most widespread area because paddy fields were the largest arable land in this region, and had the highest RR production. Following that, R1_North and R2_South had 152 and 121 areas. In sum, 494 candidate areas were proposed in total for simulating the decentralized installation. Suitable areas in R2_South were favorably located in the same spatial patterns as PR availability. For the suitability analysis, results of every region can be inferred in Supplementary Figures S2 and S3.

Site Selection from Suitability Analysis
An exemplary output map (see Figure 2) of suitable areas was classified into five levels, as illustrated in Figure 2b. The most suitable areas at level 5 were illustrated in red raster layers. We could identify the suitable level 5 into one specific candidate location per district boundary by dissolving the points of extensive layers, as presented in Figure 2c. As a result, R3_Northeast found the greatest numbers of suitable locations, resulting in 221 areas. R3_Northeastern was distributed as the most widespread area because paddy fields were the largest arable land in this region, and had the highest RR production. Following that, R1_North and R2_South had 152 and 121 areas. In sum, 494 candidate areas were proposed in total for simulating the decentralized installation. Suitable areas in R2_South were favorably located in the same spatial patterns as PR availability. For the suitability analysis, results of every region can be inferred in Supplementary Figure S2,3.

Decentralized Bioethanol Production Design from the Location-Allocation Model
With the location-allocation model, optimal biorefineries in four regions were selected from candidate suitable plants in addition to the existing operational plants. Visualizations of the connected routes from defined biorefineries to the captured residue resources are presented in Figure 2d, whereas changes in residue coverages from increasing biorefinery numbers in the regional boundaries are shown in Supplementary Figure S4. It can be noted that bioethanol facilities are likely to be chosen in high-density biomass areas as the first priority as a significant indication can be observed in R1_North. The model succeeded in allocating the feedstock biomass supply without conflicts between refineries.

Decentralized Bioethanol Production Design from the Location-Allocation Model
With the location-allocation model, optimal biorefineries in four regions were selected from candidate suitable plants in addition to the existing operational plants. Visualizations of the connected routes from defined biorefineries to the captured residue resources are presented in Figure 2d, whereas changes in residue coverages from increasing biorefinery numbers in the regional boundaries are shown in Supplementary Figure S4. It can be noted that bioethanol facilities are likely to be chosen in high-density biomass areas as the first priority as a significant indication can be observed in R1_North. The model succeeded in allocating the feedstock biomass supply without conflicts between refineries. Figure 3 demonstrated relations of biomass utilization rates and the average potential bioethanol capacities per plant corresponding to the increasing numbers of set-up installation facilities. Although increasing production units led to incremental residue utilization rates, there were peak points which represented the maximum biomass coverage (utilization rates). Thus, raising installation numbers after the maximum coverages also caused declines in biomass collection rates because a certain amount of the overall supply per region has to be shared to more producers. The case of R2_South illustrates the collectability which declined after increasing the number of facilities over the maximum due to condensed installation limited by accessible biomass supply.
As a result, we could identify the greatest collectable amounts in R1_North, R2_South, R3_Northeast, and R4_Central as 13,568, 37,304, 321,085, and 11,010 kTons year −1 , respectively. The maximum biomass coverage rates in respective regions were 91.8, 90.8, 90.5, and 79.7%. According to the outcomes of maximum collection, the optimal installation numbers were suggested to be 18, 7, 30, and 16 units in R1_North, R2_South, R3_Northeast, and R4_Central, respectively. Biomass utilization rates also created different results in different cases, both with and without considering current plants. In the case of considering only existing locations, R4_Central has showed that the biomass collectable curve almost reached the plateau after installing to approximately five units. It is possibly due to the condition of fixed locations which controlled the modeling opportunities to expand to broader areas.
R2_South, R3_Northeast, and R4_Central as 13,568, 37,304, 321,085, and 11,010 kTons year −1 , respectively. The maximum biomass coverage rates in respective regions were 91.8, 90.8, 90.5, and 79.7%. According to the outcomes of maximum collection, the optimal installation numbers were suggested to be 18, 7, 30, and 16 units in R1_North, R2_South, R3_Northeast, and R4_Central, respectively. Biomass utilization rates also created different results in different cases, both with and without considering current plants. In the case of considering only existing locations, R4_Central has showed that the biomass collectable curve almost reached the plateau after installing to approximately five units. It is possibly due to the condition of fixed locations which controlled the modeling opportunities to expand to broader areas.

Logistics Analysis and Estimated Capital Investment Cost
Delivery distances were derived from the location-allocation model. Figure 4 displays the analysis outcomes of the unitary logistics cost per bioethanol productivity associated with installation units. Every regional case study showed a lower logistics cost in line with more numbers of installation. Among four cases, the unitary logistics cost in R2_South was found to be the highest in the range of USD 1.73-2.46 L −1 . Nonetheless, R2_South displayed the largest drop in the logistics cost (almost 30%) by increasing

Logistics Analysis and Estimated Capital Investment Cost
Delivery distances were derived from the location-allocation model. Figure 4 displays the analysis outcomes of the unitary logistics cost per bioethanol productivity associated with installation units. Every regional case study showed a lower logistics cost in line with more numbers of installation. Among four cases, the unitary logistics cost in R2_South was found to be the highest in the range of USD 1.73-2.46 L −1 . Nonetheless, R2_South displayed the largest drop in the logistics cost (almost 30%) by increasing facility numbers. The model succeeded in saving logistics costs due to an improvement in searching resources with higher conversion factors, specifically referring to RRs, in addition to PRs in this case. On the other hand, R1_North, R3_Northeast, and R4_Central could benefit from the CR availability, resulting in lower logistics costs. The R4_Central showed the lowest rate in the range of USD 1.18-1.32 L −1 , which did not change significantly when increasing the number of facilities. Figure 5 illustrates correlations between the unitary logistics cost per delivered feedstocks and bioethanol production scales in different regional case studies. In every case, it is noticeable that when the production capacity per unit is smaller due to more installed facilities, the logistics cost per unitary feedstock delivery became cheaper. The results demonstrated the benefits of locating more facility units for shortening distances from feedstock resources. The cost curves vary on each region, and are structured from different geographical features and residue availabilities. Rapid increases in unitary logistics costs, from the smallest to the flat curve at different points of the production capacity, were evident. For example, R1_North reached a stable transportation cost per feedstocks after a production capacity of around 50 mL year −1 . The curve changes from logistics analysis can be used as a development framework for optimizing the productivity if the logistics cost is reduced.
facility numbers. The model succeeded in saving logistics costs due to an improvement in searching resources with higher conversion factors, specifically referring to RRs, in addition to PRs in this case. On the other hand, R1_North, R3_Northeast, and R4_Central could benefit from the CR availability, resulting in lower logistics costs. The R4_Central showed the lowest rate in the range of USD 1.18-1.32 L −1 , which did not change significantly when increasing the number of facilities.  Figure 5 illustrates correlations between the unitary logistics cost per delivered feedstocks and bioethanol production scales in different regional case studies. In every case, it is noticeable that when the production capacity per unit is smaller due to more installed facilities, the logistics cost per unitary feedstock delivery became cheaper. The results demonstrated the benefits of locating more facility units for shortening distances from feedstock resources. The cost curves vary on each region, and are structured from different geographical features and residue availabilities. Rapid increases in unitary logistics costs, from the smallest to the flat curve at different points of the production capacity, were evident. For example, R1_North reached a stable transportation cost per feedstocks after a production capacity of around 50 mL year −1 . The curve changes from logistics analysis can be used as a development framework for optimizing the productivity if the logistics cost is reduced.

Estimated Capital Investment Costs
In light of the TCI computation based on production capacity, Figure 6 shows th changing curves of the TCI according to more plant installations. As consistent with ever scenario, the TCI was cheaper when there were fewer exhibited production units in th region, which refers to a larger capacity per plant. Therefore, with more numbers of plan installation, the production capacity per unit become smaller, leading to more expensiv TCIs. The expectable ranges of the TCI in R1_North were observed to range between USD 3.33 and 6.93 L −1 as the installation plants increased from 1 to 25 plants. For R2_South, th estimated TCI ranged between USD 1.67 and 4.98 L −1 . Following that, the price range o the TCI in R3_Northeast within 1-25 plants varied from USD 1.17 to 2.48 L −1 , making the most cost-effective solution due to the large volumes of productivity in this region Among four regions, R4_Central became the highest TCI curve since it has the lowest b oethanol production potential, ranging from USD 3.51 to 6.50 L −1 .

Estimated Capital Investment Costs
In light of the TCI computation based on production capacity, Figure 6 shows the changing curves of the TCI according to more plant installations. As consistent with every scenario, the TCI was cheaper when there were fewer exhibited production units in the region, which refers to a larger capacity per plant. Therefore, with more numbers of plant installation, the production capacity per unit become smaller, leading to more expensive TCIs. The expectable ranges of the TCI in R1_North were observed to range between USD 3.33 and 6.93 L −1 as the installation plants increased from 1 to 25 plants. For R2_South, the estimated TCI ranged between USD 1.67 and 4.98 L −1 . Following that, the price range of the TCI in R3_Northeast within 1-25 plants varied from USD 1.17 to 2.48 L −1 , making it the most cost-effective solution due to the large volumes of productivity in this region. Among four regions, R4_Central became the highest TCI curve since it has the lowest bioethanol production potential, ranging from USD 3.51 to 6.50 L −1 .
3.33 and 6.93 L −1 as the installation plants increased from 1 to 25 plants. For R2_South, the estimated TCI ranged between USD 1.67 and 4.98 L −1 . Following that, the price range of the TCI in R3_Northeast within 1-25 plants varied from USD 1.17 to 2.48 L −1 , making it the most cost-effective solution due to the large volumes of productivity in this region. Among four regions, R4_Central became the highest TCI curve since it has the lowest bioethanol production potential, ranging from USD 3.51 to 6.50 L −1 .  In accordance with the optimal installation numbers based on the maximum biomass utilization rates, the results presented the placement of diverse production scales in one region. The illustration in Figure 7a-d presents geographical mapping of bioethanol scales with the estimation of a capital cost per unit in the respective region. Corresponding to the classified scales, the bubble sizes and categorized colors on the right-hand graphs portrayed the potential scales of the bioethanol refinery. In this study, the proposal construction of bioethanol scales was grouped into four different scales, i.e., 'small', 'medium', 'large', and 'very large'.
Firstly, R1_North with 18 optimal units resulted in TCIs per bioethanol production unit ranging between USD 3.11 and 20.80 L −1 . Plant installations were mostly dominated by a small-size capacity. As the northern area of R1_North is constituted with mountain ranges and high elevated lands, the outcomes of decentralized small-scales seemed to be specifically appropriated to enhance resource accessibility. Next, the TCI results in R2_South clarified the costs ranging from USD 1.91 to 7.13 L −1 . In this region, the capacity was found under the category of large-scale production. For R3_Northeast, it was possible to identify all four ranges of production scales, showing the TCI ranges at USD 1.19-5.90 L −1 . However, most of production plants fell into large and very large size productions. Large-scale bioethanol productions are generally produced from SRs and CRs where the TCIs showed relatively low rates. Lastly, the outcomes of R4_Central indicate a small-scale production, which could effectively represent decentralized facilities. The TCIs in this region altered between USD 4.37 and 18.73 L −1 . In accordance with the optimal installation numbers based on the maximum biomass utilization rates, the results presented the placement of diverse production scales in one region. The illustration in Figure 7a-d presents geographical mapping of bioethanol scales with the estimation of a capital cost per unit in the respective region. Corresponding to the classified scales, the bubble sizes and categorized colors on the right-hand graphs portrayed the potential scales of the bioethanol refinery. In this study, the proposal construction of bioethanol scales was grouped into four different scales, i.e., 'small', 'medium', 'large', and 'very large'.
Firstly, R1_North with 18 optimal units resulted in TCIs per bioethanol production unit ranging between USD 3.11 and 20.80 L −1 . Plant installations were mostly dominated by a small-size capacity. As the northern area of R1_North is constituted with mountain ranges and high elevated lands, the outcomes of decentralized small-scales seemed to be specifically appropriated to enhance resource accessibility. Next, the TCI results in R2_South clarified the costs ranging from USD 1.91 to 7.13 L −1 . In this region, the capacity was found under the category of large-scale production. For R3_Northeast, it was possible to identify all four ranges of production scales, showing the TCI ranges at USD 1.19-5.90 L −1 . However, most of production plants fell into large and very large size productions. Large-scale bioethanol productions are generally produced from SRs and CRs where the TCIs showed relatively low rates. Lastly, the outcomes of R4_Central indicate a smallscale production, which could effectively represent decentralized facilities. The TCIs in this region altered between USD 4.37 and 18.73 L −1 .
In addition to formulating various scales of biorefineries, we projected the potential scenarios under different assumptions of biomass utilization rates (BURs). The reason for this scenario projection is based on assuming that the participation of biorefineries may not reach 100% of the proposal in reality, due to the fact that smaller scales of producers need to pay unitary TCIs at higher rates. If the biomass utilization may be prioritized, the shadowed areas of BUR at around 80% and 50% indicated the production facilities that are encouraged to be established. For instance, the coverage of BUR 85% in R1_North alluded to biorefinery No. 10, 11, 12, 14, 15 and 17, representing a set of small and medium production units, while BUR 50% covered only two units including No. 10 and 17.

Emissions and Pollutants Abatement from Biomass Collection
By increasing the numbers of bioethanol production facilities, the amount of residue collection is expected to contribute to expectable GHGs and pollutant eliminations. The result of GHG emission abatement curves, as compared to the emission rates without biomass utilization under different numbers of refineries, is presented in Supplementary  Figure S5. The maximum generated values of complete combustion from left-over biomass was found to emit 9.1, 0.5, 167.1, and 9.4 M ton CO2 eq of GHG emissions in R1_North, R2_South, R3_Northeast, and R4_Central, respectively. Nonetheless, when the biomass residues were collected, significant abatement could be observed. Table 4 summarized GHG and particulate matter (i.e., PM10 and PM2.5) reduction potential from the maximum biomass collection in each location. A comparison among four regions unveiled that R3_Northeast resulted in the highest GHG, with a saving potential of PM10 and PM2.5 at 153.7 M tons of CO2-eq year −1 reaching 605.3 and 552.8 k Tons year −1 , respectively. An estimation of this region is around 90% of a whole estimated national emission. In addition to formulating various scales of biorefineries, we projected the potential scenarios under different assumptions of biomass utilization rates (BURs). The reason for this scenario projection is based on assuming that the participation of biorefineries may not reach 100% of the proposal in reality, due to the fact that smaller scales of producers need to pay unitary TCIs at higher rates. If the biomass utilization may be prioritized, the shadowed areas of BUR at around 80% and 50% indicated the production facilities that are encouraged to be established. For instance, the coverage of BUR 85% in R1_North alluded to biorefinery No. 10, 11, 12, 14, 15 and 17, representing a set of small and medium production units, while BUR 50% covered only two units including No. 10 and 17.

Emissions and Pollutants Abatement from Biomass Collection
By increasing the numbers of bioethanol production facilities, the amount of residue collection is expected to contribute to expectable GHGs and pollutant eliminations. The result of GHG emission abatement curves, as compared to the emission rates without biomass utilization under different numbers of refineries, is presented in Supplementary  Figure S5. The maximum generated values of complete combustion from left-over biomass was found to emit 9.1, 0.5, 167.1, and 9.4 M ton CO 2 eq of GHG emissions in R1_North, R2_South, R3_Northeast, and R4_Central, respectively. Nonetheless, when the biomass residues were collected, significant abatement could be observed. Table 4 summarized GHG and particulate matter (i.e., PM 10 and PM 2.5 ) reduction potential from the maximum biomass collection in each location. A comparison among four regions unveiled that R3_Northeast resulted in the highest GHG, with a saving potential of PM 10 and PM 2.5 at 153.7 M tons of CO 2 -eq year −1 reaching 605.3 and 552.8 k Tons year −1 , respectively. An estimation of this region is around 90% of a whole estimated national emission.

Spatial Distribution and Biomass Allocation to Decentralized Facilities
The bioethanol potential convertibility from SRs, CRs, RRs, and PRs in the decentralized units was previously evaluated between 0.38 and 1.39 L bioethanol per kg of dry matter, based on the biomass production in 2018 [13]. As such, this study proved that the conversion of those figures can utilize scattered crop residues to potentially reach 13,444 million L of bioethanol from the maximum collectable residues.
The estimated amount of bioethanol can not only substitute fossil fuels but also reduces additional GHG emissions, when the agricultural residues are no longer burned on the field. We could prove these benefits in terms of the GHG emission reduction by transitioning from centralization to decentralization. Taken as the preliminary estimation, if the generated agricultural residues are completely utilized in bioethanol production, the CO 2 abatement from substituting the fossil-based gasoline can be expected to reach around 21.99 M tons of CO 2 (calculated based on emissions from direct combustion) on top of the abatement from avoiding the open-burning action [91]. However, in order to scrutinize the precise reduction effect, an overall supply chain of bioethanol production needs to be accounted for the life-cycle GHG emissions [92,93].
In this paper, the spatial features of input variables were created based on 1 × 1 km 2 in a high resolution. The constructed grid cells from geographic layers successfully enabled the function of suitability analysis and location-allocation modeling. By accommodating GIS-MCA, the characteristics of suitable spots in respective regions can effectively reflect the suitable areas with the high density of agro-residues. For example, R1_North highly recommended areas that were mostly in the lower part of the north, i.e., non-mountainous areas with a high density of biomass availability and cultivation, as well as industryoriented land use types. Similarly, suitable areas in R2_South are restricted by preserved forests and extended mountain ranges in spite of high densities of palm plantation.
Although accuracy suitability analysis can be improved from applying additional spatial information (e.g., water sources or power utility systems) [94], the geographic results from suitability analysis in this study are in line with the suitable locations assessed by Cheewaphongphan et al., who addressed the locations of biomass power plants from a rice straw for different utilization technologies [15]. Based on candidate points, when existing production has not been originally situated, there is more freedom for the model to seek out the most optimal sites and expand catchment areas to collect residues.

Logistics for Biomass Supply Chain
In this study, the connection between biorefineries and supply sources was assigned based on the actual spatial information of road network in Thailand. This approach effectively and accurately determined the adjacent resources allocated to producers. The approach in the study is distinctive from some of the literature [44,95] based on the assumption of limiting radial distances to avoid exceeding 80 km. However, according to the modeling execution, acquired output data ranged between 37 and 55 km. These figures demonstrate that the highest economic viability in biomass-to-biofuel industries is between 45 and 80 km [40,51].
With the restriction of catchment distances, it was possible to limit the rising production scales from expanding a longer collecting radius, and also allowed the study to expand biomass spatial coverage from multiplying production units instead of expanding catchment areas. With our simplified logistic analysis approach, we proved that increasing installation units can improve average logistics costs. This significant result could be observed in the case of R2_South, as presented in Figure 7. Although the case of R2_South resulted in the most expensive logistics cost, the result revealed that R2_South contributed the most significant change following expanding numbers of production refinery per region.
Although trend lines of unitary biomass logistics per delivered feedstock are relatively similar across regions, as presented in Figure 5, R3_Northeast seemed to gain the most advantageous outcomes owing to the high density of biomass availability and the high potential of the bioethanol conversion ratio. The results also imply that bioethanol convertibility and density per area are influential factors for estimating unitary logistics costs. Thus, the possibility of reducing logistics costs per capacity could be achieved by blending input materials with a high convertible biomass, leading to a higher bioethanol yield.

Opportunities of Decentralization and Prospect for the Future Development
Decentralized schemes, as proposed in this work, proposed 71 units of decentralized biomass refineries in total. The conceptual designs could be proposed by depicting production units which consolidate different scales (small, medium, and large), as shown in Figure 7. Since the study did not prioritize the crop type selection, the concept of mixed feedstock could be expected to prevail against the limitations from single-type feedstock, especially regarding the harvesting seasonality. In real practice, SRs are collectable around November to March, whereas CRs can be collected between March and May [96]. For RRs, the harvesting period is between October and December [49]. Therefore, the heterogeneous feedstock utilization can potentially compensate production yields between cultivation cycles.
Considering the regions without current invested facilities, R1_North could be used for initiating stand-alone and small-scale production facilities, as R1_North consists of many remote areas with obstacles including high elevation and land slope. Meanwhile, R2_South has more advantages in terms of PRs that can be used with the highest biomass-to-crop ratio. In some areas of R2_South, medium-to-very large scales can be used due to the high generated density of PRs, which require mass production facilities. For R3_Northeast and R4_Central which have operating plants, the existing locations could possibly be used for supplying alternative feedstocks.
In decentralized production scales, the TCI is also a vital factor for determining the project investment scales. Taking R4_Central as an example, TCIs from newly developing plants may cost USD 4.37-18.73 L −1 , which is considered to be around 2-4 times more expensive than some pioneer second-generation plants [97,98]. To deal with the economic limitation, integrating the first-and the second-generation technologies could encourage a reduction in capital investment. Another possible solution concerns the adjustment of radii for the catchment area to be larger than 80 km in order to increase the quantity of feedstocks and production scales.
Despite economic challenges, the study demonstrated opportunities of encouraging local producers and farmers to potentially participate in future project developments. To encourage the proposed scheme for policy implications, incentivized mechanisms may become necessary in order to persuade farmers to deliver post-harvested biomass into the bioethanol production chain. One of the important driving factors regards the feedstock selling price that can attract farmers, so they can benefit as value-added incomes along with the main crop production.

Sensitivity Analysis
The sensitivity analysis was conducted to observe the effects of variables involving biomass availability, delivery distances, delivery costs, and the fluctuating supplying feedstocks. Firstly, from the study of logistic cost sensitivities, as shown in Figure 8, the most sensitive parameter in every case study was indicated to the delivery distance. The similar trends across regions demonstrate that the costs increase as the delivery distances become longer, and vice versa. The next influential factor was the delivery cost per kilogram per kilometer. On the other hand, it is noticeable that the influence of varying quantities of biomass supply were hardly affected by the cost changes within ±30% from the base case. The TCI directly relies on the production scales. Hence, the densities of available biomass and bioethanol conversion factors from different types of biomass are important indicators for estimating the production capacity. We demonstrated the absence of each kind of biomass to detect the effects from different types of input feedstocks. Without each type of biomass, the TCI soared, particularly from the scenario without SRs in R1_North, R3_Northeast, and R4_Central. For R2_South, unavailable PRs affected the TCI the most, as it costs more than the double price of the TCI from the base case scenario to over USD 6.00 L −1 .  Another clarification is associated with the missing specific type of input feedstocks due to the seasonality or supply allocation. As a result, except for the R2_South, logistics costs in every region are sensitive to unavailable SRs, since SRs offer benefits including the high density of availability. On the contrary, the exclusion of RRs resulted in lower logistics costs, as displayed in Figure 8a,c. The exclusion of RRs positively affected the logistics costs due to the fact that it offered the lowest bioethanol convertibility. Interestingly, the most significant logistics cost reduction in R2_South was offered from factoring out PRs as input feedstocks owing to the massive quantity of collectable residues which increases the transportation cost. Hence, reducing the input PRs can also be beneficial in the logistics cost reduction.
In addition to the sensitivity of logistics costs, Figure 9 showed the sensitivity analysis results of the TCI according to changing variables. As a result of changing biomass availability, among four regions, R3_Northeast demonstrated the most sensitive TCI changes, as the price ranges showed the largest span of the lowest and the highest price, i.e., USD 1.89 L −1 to USD 2.42 L −1 . On the other hand, the least sensitive scenario depending on biomass availability was discovered in R2_South.  The TCI directly relies on the production scales. Hence, the densities of available biomass and bioethanol conversion factors from different types of biomass are important indicators for estimating the production capacity. We demonstrated the absence of each kind of biomass to detect the effects from different types of input feedstocks. Without each type of biomass, the TCI soared, particularly from the scenario without SRs in R1_North, R3_Northeast, and R4_Central. For R2_South, unavailable PRs affected the TCI the most, as it costs more than the double price of the TCI from the base case scenario to over USD 6.00 L −1 .

Conclusions
The agricultural residues from post-harvesting sugarcane, cassava, rice, and palm were proposed as feedstocks for bioethanol production. In this study, the geographical variables were created at 1 × 1 km 2 grid-scale resolutions, providing mapping information of spatial biomass availability. By implementing the integrated GIS-based suitability analysis and prioritizing the techno-economic viability, 494 candidate locations were suggested in the model. The goal to maximize biomass collection rates was achieved through the locationallocation model to determine the feasible numbers of installation units, resulting in 18, 7, 30, and 16 optimal facilities in the north, south, northeast, and center of Thailand, respectively. In total, 71 locations were chosen in this work, revealing the optimal numbers of installation facilities that cover 80-92% of spatially available residual biomass rates. The amount of collected residues could contribute to reductions of 170 M tons of CO 2 -eq GHGs and 1271 k tons of particulate matters, suggesting potential solutions to abolish open-field burning of crop residues.
By installing more decentralization units, cutting logistic costs represents a decline of 10-30%. Nonetheless, exhibiting more biorefineries causes trade-offs from investment costs along with benefits. The range of possible logistics costs in the scope of the nation within 1-25 plant installations was indicated at USD 1.17-2.46 L −1 , whereas the total capital investment was in the range of USD 1.17-6.93 L −1 . The results can be used as an economic framework for further analysis to verify the minimum bioethanol selling price in the future work.
The proposal of combining different biorefinery scales in respective regions in this study could be taken as a recommendation for policy planning or commercialization in terms of location selection and for determining favorable scales based on optimal installation units. In order to drive the scheme into national policy, plants with a high production capacity may be taken as an installation priority in order to cover 50-85% of available biomass per region, which could be economically attractive for encouraging bioethanol producers in relation to investment costs.
Even though comprehensive operational cost analysis from production process still needs to be assessed, the outcomes of logistics costs and capital costs per ethanol production units can be used as baselines to capture the minimum bioethanol selling price. Future research requires further scenario-based analysis by optimizing decentralized schemes based on crossing dimensions of techno-economic and environmental aspects, which is expected to sustainably enhance lignocellulosic bioethanol development.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/su14169885/s1. Figure S1: Spatial allocation of agricultural residues indicating per 1 km 2 grid cell in (a) R1_North, (b) R2_South, (c) R3_Northeast, and (d) R4_Central case study. Figure S2: Suitability level indicating in (a) R1_North, (b) R2_South, and (c) R3_Northeast. Figure S3: Identification of candidate plants to be used as newly proposed decentralized bioethanol production identifying in 3 regions in (a) R1_North, (b) R2_South, and (c) R3_Northeast. Figure S4: Chosen biorefineries according to location-allocation models under the different assumption of installation numbers resulting in connections to biomass allocation in the catchment areas. Figure S5: GHG and atmospheric pollutants emission abatement of (a) R1_North, (b) R2_South, (c) R3_Northeast, and (d) R4_Central case study. Table S1: Adjusted cost input for using the logistics analysis model. Table S2

Data Availability Statement:
The data that support the findings of this study are available upon reasonable request from the corresponding author, P.J.