2.2.1. Land Suitability Model
In the land suitability model, lands dedicated to the energy crop and ethanol processing were identified by considering land topography, current land use, soil temperature and moisture regimes, rainfall characteristics and other climate factors, etc. The model was built on several important factors and criteria identified by the Clean Fuels Development Coalition [
4,
13,
14,
16,
28,
29,
30]. The model addresses two objectives: (i) identify possible land locations for banagrass production and estimating corresponding banagrass yield at each location; and (ii) determine possible locations for ethanol processing plants.
Data preparation is critical in order to exclude areas that are not feasible for banagrass production by employing a series of overlay operations. Banagrass is known as a tropical grass found in most tropical or subtropical countries [
31]. High-altitude areas with elevation above 2000 m were found to be unsuitable for banagrass production and were not considered as land candidates. Areas affected by natural hazards such as volcanoes or reserved for different purposes such as forest reserves, sanctuaries, parks, historic sites, military installations, and others were also excluded. To avoid competing with lands currently growing agricultural crops and engaged in related activities, these were excluded. Land areas classified as pasture land and non-prime, non-unique agricultural lands, as categorized in ALUM and ALISH map layers, were selected as potential land candidate for banagrass production.
This study assumed mechanized banagrass production management and harvesting. Therefore, land types, slope gradient, and road conditions play important roles and were taken into consideration. All land areas categorized as sandy, very rocky or stony, extremely stony, and rough broken land in the soil survey [
20] were excluded from possible production areas. Slope gradient measured as percentage rise was calculated from the Shuttle Radar Topography Mission (SRTM) 90-m resolution Digital Elevation Model (DEM). Different harvesting techniques require different slope limits in order to operate machines [
32]. A slope gradient less than 20% was found to be suitable for grass crops [
13]. As a result, lands located in areas with slope gradient exceeding 20% were not considered. Biomass feedstock transportation requires frequent and heavy truck hauling operations. Therefore, road conditions play a critical role. Land areas adjacent to existing major highways classified as multiplane divided expressways or principal highways on Hawaii and Maui Islands were selected [
19].
The next important step in the land suitability model is the ground-truthing process to validate the potential production land areas resulting from the operations. Since the map layer data of ALUM and ALISH were produced in the 1970s, it is possible there could be more recent residential development areas mixed with other land projects. The resulting land areas were then cross-checked with TMK parcel maps (most recently updated in 2013) by the State Office of Planning to define specific land owners and land use purposes [
19]. The areas characterized as agricultural and rural lands by attribute were selected for further consideration. The final land selection was based on the process of ground-truthing to help improve the quality of the analysis by using high spatial resolution satellite imagery from Google Earth and Microsoft Aerial Bing Maps. All the shapefile layers were projected to the coordinate system NAD 1983 UTM Zone 4N and processed for Hawaii and Maui Islands. The result was the shapefile layer presenting possible land areas for growing banagrass feedstock. Data for these areas were retrieved for land suitability analysis and LP model.
The integration of AHP and GIS for the land suitability model has been widely used in various studies, such as landfill selection in Konya’s Lake Beyşehir catchment area, Turkey [
33], urban development in a secondary city of Bangladesh [
34], and site evaluation for ecotourism [
35]. AHP is a multi-criteria decision-making technique based on decision theory to determine which criteria are the most important. This paper uses a technique to evaluate land suitability levels for banagrass production. The AHP model reports the weight of each factor and the consistency of the results is validated by calculating a consistency ratio (CR). The threshold CR value is 0.1. A CR value less than 0.1 indicates an acceptable level of consistency in pairwise comparisons. A CR value greater than or equal to 0.1 indicates that the criteria should be reexamined until CR reaches the consistency level less than 0.1. The use of other studies, expert judgments, and experimental data are important to determine the ranking system for the pairwise comparison matrix in AHP.
The first step is to identify the critical factors affecting banagrass growth. Four critical factors were considered in this study, including soil temperature and moisture regimes, precipitation, minimum temperature, and solar radiation. These critical factors were selected based on a series of studies by Black and Veatch (2010), Ferraris (1978); Ferraris and Sinclair (1980); and Ferraris, Mahony and Wood (1986) [
1,
36,
37,
38]. Each critical factor has its own set of attributes that also has effects on banagrass growth. These attributes were retrieved according to the conditions on Hawaii and Maui Islands and categorized into different groups (
Table 1). Soil temperature and moisture regimes were divided into different groups by combining temperature conditions and moisture regimes such as hot and wet; hot and moist; hot and saturated; etc. Other factors—precipitation, minimum temperature, and solar radiation—were categorized into five groups by applying the commonly used Jenks natural breaks classification method [
39] in ArcGIS, which is a data clustering method designed to determine the best arrangement of values for different classes. This was done by minimizing each class’s average deviation from the class mean, while maximizing each class’s deviation from the means of the other groups. In other words, this method reduces the variance within classes and maximizes the variance between classes.
The next step was to determine the weights of the four critical factors and the weights of attributes within each critical factor. The primary objective for this step in AHP analysis is to determine the ranking of criteria and to compare their rankings by constructing a pairwise comparison matrix. The numerical values 1, 3, 5, 7, and 9 were assigned to indicate if one factor is thought to be equally important, moderately more important, strongly more important, very strongly more important, and strictly superior in comparison with the other criteria in a pairwise comparison, respectively. Given the four critical factors and attributes, there were five sets of AHP results to estimate the weight of each factor and the weight of each attribute within each factor. This study used both (i) experimental data and (ii) expert’s judgments for this step.
- (i)
The experimental data were collected from several banagrass field trials located on the Hawaiian Commercial & Sugar Company (HC&S) land on Maui Island. Monitoring of banagrass yields, soil characteristics, and daily weather conditions on the field trials was done from 2012 to 2015 as a part of the USDA/NIFA Biomass Research and Development Initiative (BRDI) project at the University of Hawaii at Manoa. These study results are also a part of the BRDI project.
- (ii)
Dr. Richard Ogoshi, Department of Tropical Plant & Soil Sciences, University of Hawaii at Manoa, is the task leader in charge of working on banagrass field trials. He is an agronomist/soil scientist who works with researchers from HC&S to monitor, collect, and analyze the effects of weather conditions and soil characteristics, etc. on banagrass yields over the years. For Hawaii, he is the lead expert in this field. We consulted with him regarding determination of the ranking system for the pairwise comparison matrix that is used in conducting AHP analysis for the four critical factors and their attributes.
Dr. Ogoshi is the leading expert who works with researchers (other experts) from a different organization (HC&S) on the same research project (BRDI). A recent AHP study by de FSM Russo and Camanho (2015) reviewed the literature (33 other AHP studies) thoroughly but used only one expert in their AHP analysis [
40]. Another study, Tsyganok et al. (2012), suggests that “…if the number of experts in a group is relatively small, expert competence should always be taken into consideration” [
41].
Of the four critical factors, soil temperature and moisture regimes appear to be the most important factors, followed by precipitation. Minimum temperature and solar radiation are less important as compared to other factors. Therefore, a value of 5 is assigned to indicate that soil temperature and moisture regimes are thought to be strongly more important than precipitation. A value of 3 meant that precipitation is moderately more important than minimum temperature or solar radiation. The same procedure was done to create the full pairwise comparison matrix for AHP analysis for the critical factors.
This procedure was also applied for AHP analyses to determine the weights of attributes for each critical factor. For example, higher banagrass yield likely occurred in the locations with higher soil temperature and moisture conditions. Hot and wet soil is thought to have a higher yield as compared to other attributes of soil temperature and moisture regimes. A value of 3 indicated that hot and wet soil is thought to be moderately more important than hot and moist soil; a value of 5 meant hot and wet soil is thought to be strongly more important hot and saturated soil; or a value of 9 showed that hot and wet soil is strictly superior to cold and moist soil, etc. For other critical factors, higher banagrass yield is thought to be associated with higher solar radiation, higher minimum temperature and precipitation ranging from 1000 mm to 5000 mm as compared to other attributes within the solar radiation, minimum temperature, and precipitation factors. All these conditions were used in building pairwise comparison matrices to estimate the weights of factors and their attributes that determine suitable areas for banagrass production.
The AHP results were then inputted to the related layers for weighted overlay operations processed in ArcGIS version 10.1. These overlay operations were used to locate suitable banagrass production areas for Hawaii and Maui Islands.
The estimation of banagrass yield in natural conditions was based on the study by [
1] and validated by data collected from several field trials on Maui during the period from 2012 to 2015 [
42]. The estimated yield data were then compiled to estimate potential biomass production on Hawaii and Maui Islands and the potential ethanol production for the state of Hawaii based on the banagrass biomass feedstock.
Proper siting of an ethanol plant is critical to economically develop renewable energy alternatives. To account for logistical, social, and environmental concerns, a guide was developed by the Clean Fuels Development Coalition (CFDC) and the Nebraska Ethanol Board in 2006. This guide [
28] suggests that the land area for an ethanol processing plant should be at least 40 acres. This land area should provide enough space for on-site transportation and future development of the ethanol plant. Properly located, this plant should not impose potential problems with neighboring community and residential areas. Following this requirement, the ethanol plant location should also be within the biomass feedstock production areas prescribed by the land suitability model. Land areas less than 40 acres and/or near to community and residential locations were excluded from further plant location consideration. Other criteria used in the plant location process include transportation and water requirements [
28]. Possible processing sites should include a 3200-m buffer zone around major highways [
16,
30] and availability of aquifer systems to meet sustainable yields for water demands in 2025 [
16]. The presence of sustainable water aquifers was stated in the Hawaii and Maui counties’ water use and development plans [
43,
44]. These data were processed and resulted in land locations that meet all requirements for the ethanol processing plant site selection. These processing ethanol plants were located at the center of each island’s banagrass production sites.
2.2.2. LP Model
The overall objective of the LP model is to minimize the cost of ethanol production in order to meet 20% of the state’s highway fuel demand by 2020. Results from the LP model include the lowest estimated unit cost for ethanol production, the land areas and their locations needed to produce the target amount of banagrass feedstock such that this feedstock can be processed into ethanol to meet the ethanol demand targets. The framework for the LP model is presented in
Figure 3. The program involves several cost components, including the costs of banagrass feedstock production, feedstock transportation, ethanol conversion, and ethanol transportation.
To obtain data for the LP model, 1-km
2 grid cells were created and used to delineate possible locations for banagrass production. The land areas and corresponding estimated banagrass yield that falls within the grid cells were retrieved as inputs for the LP model. The center point of each grid cell was generated in order to estimate distances from production areas to the ethanol processing plants and then to harbor ports on each island using network analysis. This analysis aims to minimize transportation costs by using the calculated distances from the center of each banagrass production grid cell to destination locations through the origin-destination (OD) cost matrix analysis integrated in ArcGIS version 10.1. The detailed street map for Hawaii and Maui Islands [
45] were used to create this network dataset. The OD cost matrix calculates least-cost paths from production areas to the ethanol processing plants, from production areas to harbor ports and from ethanol processing plants to harbor ports for each island. These data are the inputs to the cost analyses for feedstock production and transportation, and ethanol processing and transportation.
The first component of the cost analysis is the cost of banagrass feedstock production, assuming that feedstock production occurs on both Hawaii and Maui Islands. GIS data on possible locations of land available for banagrass production and the corresponding yields associated with each location are used as inputs in the program. Let x
i be the land area variable at each location i on the island of Hawaii, i ∈ {1, 2, …, n} and xj be the land area variable at each location j on the island of Maui, j ∈ {1, 2, …, m}. These variables represent the land areas selected for feedstock production, measured in km
2. The feedstock production equation is shown as:
where a
i and y
i denote the maximum land area availability and the estimated banagrass yield at production location i on Hawaii Island, respectively; a
j and y
j represent the maximum land area availability and estimated yield at production location
j on Maui Island, respectively; and Y
1 and Y
2 are the total amount of feedstock in dry tons produced on Hawaii and Maui Islands, respectively. The costs of feedstock production on Hawaii and Maui Islands are expressed as C
F1 and C
F2, respectively, and calculated via Equation (2):
where c is the cost of feedstock production per hectare (ha) per year. This includes the costs of land preparation, planting, field operations such as fertilizers, herbicides for weed control and other field operations, harvesting, other operations, and operating overhead [
1,
14].
The second component of the cost analysis involves banagrass feedstock transportation from production areas to the ethanol processing plant. Assuming that ethanol will be processed at one processing plant on either Hawaii or Maui Islands, the transportation logistics and costs are comprised of two parts—ground transportation and interisland shipping. Distances for feedstock ground transport are measured in kilometers (km) and are calculated based on a round-trip basis. The transportation cost calculation depends on where ethanol processing occurs. Equation (3) shows the case where ethanol is processed on Hawaii Island.
Feedstock produced on Hawaii Island is transported directly to the ethanol processing plant. Feedstock produced on Maui Island will be hauled by heavy truck-trailers from production areas to Kahului port, then shipped to Hilo or Kawaihae or both ports on Hawaii Island before being transported to the ethanol processing plant. Variables C
FT11, C
FT12, and C
FT13 are the costs of feedstock transportation from production areas on Hawaii Island to the processing plant; production areas on Maui Island to Kahului port; and inter-island shipping and ground transportation from Kahului port on Maui Island to the ethanol processing plant on Hawaii Island through Hilo or Kawaihae ports. Variables D
i and D
jk represent distances (measured in km) from feedstock production locations i on Hawaii Island to the ethanol processing plant and feedstock location j on Maui Island to Kahului port. The terms α and β are the fixed and variable ground transportation costs (measured in $/dry ton and $/dry ton/km); W
11 and W
12 denote the amount of feedstock shipped to Hilo and Kawaihae harbors on the island of Hawaii; D
h and D
k denote distances from Hilo and Kawaihae ports to the ethanol plant; and γ is the interisland shipping cost measured in $/dry ton.
Equation (4) presents the case for ethanol processing on Maui Island. For this case, inter-island shipping occurs in order to transport feedstock produced on Hawaii Island to Maui Island for ethanol processing. The transportation logistics include feedstock ground transportation from production areas to either Hilo or Kawaihae ports on Hawaii Island and inter-island shipping to the Kahului port on Maui Island and ground transportation to the ethanol processing plant. The transportation costs are featured as CFT22, CFT23, and CFT24. Feedstock produced on Maui Island is transported directly to the ethanol processing plant and the transportation cost is CFT21. Dj denote distances from feedstock production sites on Maui to the ethanol processing plant; Dih and Dik represent distances from feedstock production sites on Hawaii Island to Hilo and Kawaihae ports; and Dkp indicates distance from Kahului port to the ethanol processing plant on Maui Island.
The total feedstock transportation cost is given in Equation (5):
The third component is ethanol conversion costs. Since feedstock is transported to the ethanol plant immediately after harvesting, it is assumed that the moisture content remains the same as it was at the time of harvesting [
14,
46]. Cleaning and size reduction processes are accomplished at the ethanol plant. Ethanol production cost is calculated on a dry basis and given as:
where C
E is the ethanol processing cost (in $); δ is the rate of banagrass size reduction at the plant (in %); and c
ef and c
vf represent the fixed cost and variable cost of production per dry ton of feedstock ($/dry ton), Respectively.
The last component is the cost of ethanol transportation to Oahu for blending. Let ψ denote the ethanol conversion rate. Ethanol will be transported to either Hilo or Kawaihae ports if ethanol processing occurs on Hawaii Island or to Kahului if ethanol processing is on Maui Island. The produced ethanol will then be shipped to Oahu for blending with gasoline to produce E10 automobile fuel. Total ethanol produced and the transportation costs to the blending sites are calculated as:
where C
ET1 and C
ET2 denote the ethanol transportation cost on Hawaii and Maui Islands, respectively; c
et is the ethanol transportation cost per km; and c
e1 and c
e2 are the inter-island shipping cost per gallon ($/gallon/km) from Hawaii Island to Oahu and from Maui Island to Oahu.
From Equations (1)–(7), the cost optimization objective function is:
All variables are positive. Parameters for the analysis were proxied by using estimates from related and relevant studies. All monetary values are inflated to their 2014 value by using CPI indexes. These monetary variables include the estimated banagrass production cost of $1386.27/acre [
14], ethanol processing variable and fixed annual costs for a 1MGY facility with corresponding costs of $1.19/gallon and $600,700, respectively [
46], feedstock transportation variable and fixed costs of $0.18/dry ton/km and $6.79/dry ton, respectively [
47], ethanol distribution cost of $0.18/ton/km [
48], feedstock inter-island shipping cost of $136.67/wet ton [
49] and ethanol inter-island shipping costs of $0.15/gallon and $0.11/gallon from Hawaii and Maui Islands to Oahu Island, respectively [
50].