Land Suitability and Insurance Premiums: A GIS-based Multicriteria Analysis Approach for Sustainable Rice Production

: The purpose of this research is to develop a land suitability model for rice production based on suitability levels and to propose insurance premiums to obtain maximum returns based on the harvest index and subsidy dependence factor for the marginal and moderately suitable lands in the northern part of Bangladesh. A multicriteria analysis was undertaken and a rice land suitability map was developed using geographical information system and analytical hierarchy process. The analysis identiﬁed that 22.74% of the area was highly suitable, while 14.86% was marginally suitable, and 28.54% was moderately suitable for rice production. However, 32.67% of the area, which was occupied by water bodies, rivers, forests, and settlements, is permanently not suitable; 1.19% is presently not suitable. To motivate low-quality land owners to produce rice, there is no alternative but to provide protection through crop insurance. We suggest producing rice up to marginally suitable lands to obtain support from insurance. The minimum coverage is marginal coverage (70%) to cover the production costs, while the maximum coverage is high coverage (90%) to enable a maximum return. This new crop insurance model, based on land suitability can be a rational support for owners of different quality land to increase production.


Introduction
Bangladesh is ranked fifth among the countries most exposed to natural disasters, including floods, cyclones, and droughts. The production losses of major cereal crops due to natural disasters over the past 29 years have been equivalent to an average of 6.4% of the national crop production every year [1]. The Bangladesh Climate Change Strategy and Action Plan 2009 has identified "insurance" as an effective disaster management tool and one of the 37 programs that various ministries/agencies plan to undertake to combat the adverse effects of climate change in Bangladesh [2]. A key challenge for the poor farmers of Bangladesh is the management of crop damage risk. One potential strategy to cope with the devastating effects of crop damage is to introduce a crop insurance scheme. The continuing growth of the population and consumption mean that the universal demand for food will increase for at least another 40 years [3]. Agriculture in Bangladesh is already under pressure, both from large and increasing demands for food as well as from obstacles related to the degradation of agricultural land [4]. Therefore, an increase in production with less land has become a global challenge [5]. In this regard, a land suitability analysis (LSA) could be the most appropriate approach to evaluate and fulfill the required criteria for these facilities.
owners suffer more when production is lower, prices fall, or crops are damaged by natural calamities in comparison to the other farmers. There are significant numbers of crop insurance products globally, but none of them coincide directly with the Bangladesh agriculture context, and no insurance products consider the quality or suitability of land for a specific crop. Almost all crop insurance considers the weather and, at times, natural calamities. Considering the limitations of data, knowledge, the awareness of the farmers, costs, and experience, this research attempts to propose a new insurance premium model using land suitability, harvest index, and the subsidy dependence factor through community involvement. To the best of our knowledge, no research has been undertaken on crop insurance premiums based on land suitability, harvest index, and subsidy dependence factors.
Therefore, the objective of this study was to first identify the suitable areas for rice production in poverty and risk prone areas, and second, to propose a crop insurance premium model based on land suitability, harvest index, and subsidy dependence factors.

Study Area and Data Collection
The study was conducted in the Rangpur, Lalmonirhat, and Kurigram districts of Rangpur Division, which are in the Teesta and Dharla river plains and are known as poverty-prone risk areas ( Figure 1). Approximately 47.23% of the population lives below the upper poverty line, while the national average is 24.3% [33]. This region has an area of 5879.60 km 2 and is located in the northern part of Bangladesh between 25°18′ and 26°33′ N latitudes and 88°56′ and 89°54′ E longitudes.
The area consists of 19 administrative units in 3 districts, with Rangpur, Lalmonirhat, and Kurigram having an overall population of 6.3 million [34]. The elevations range from 5 to 30 m above sea level. Agriculture is the main source of income for the population. Based on weather data, the minimum and maximum mean annual temperatures vary between 8.47 and 36.3 °C. The annual average rainfall recorded is 1932-2931 mm, with high humidity at 41%-77% [35]. This work attempted to develop an agricultural insurance premium model considering suitable lands for rice cultivation in the northern part of Bangladesh. Basically, the area is rice intensive and The area consists of 19 administrative units in 3 districts, with Rangpur, Lalmonirhat, and Kurigram having an overall population of 6.3 million [34]. The elevations range from 5 to 30 m above sea level. Agriculture is the main source of income for the population. Based on weather data, the minimum and maximum mean annual temperatures vary between 8.47 and 36.3 • C. The annual average rainfall recorded is 1932-2931 mm, with high humidity at 41%-77% [35].
This work attempted to develop an agricultural insurance premium model considering suitable lands for rice cultivation in the northern part of Bangladesh. Basically, the area is rice intensive and most of the farmers cultivate three types of rice (Aus, Aman, and Boro) all year round ( Figure 2). Major production influencing factors were used for rice land suitability analysis in the geospatial platform. Subsequently, considering the importance of the factors, the availability of data, the literature review, and expert opinions, the study identified eight criteria (slope, elevation, topsoil texture, land type, soil pH, temperature, rainfall, and flood prone) and one constraint (land use) for the analysis. Land use was identified with settlements, rivers, forests, and waterbodies as constraints. Different statistical data were obtained from the Bangladesh Bureau of Statistics (BBS). Rainfall and temperature data were collected from the regional meteorological stations, and data for agricultural production were collected from the Department of Agricultural Extension (DAE). The base maps of different factors were collected from the Local Government and Engineering Division (LGED) and the Bangladesh Country Almanac (BCA) to carry out this research work. most of the farmers cultivate three types of rice (Aus, Aman, and Boro) all year round ( Figure 2). Major production influencing factors were used for rice land suitability analysis in the geospatial platform. Subsequently, considering the importance of the factors, the availability of data, the literature review, and expert opinions, the study identified eight criteria (slope, elevation, topsoil texture, land type, soil pH, temperature, rainfall, and flood prone) and one constraint (land use) for the analysis. Land use was identified with settlements, rivers, forests, and waterbodies as constraints. Different statistical data were obtained from the Bangladesh Bureau of Statistics (BBS). Rainfall and temperature data were collected from the regional meteorological stations, and data for agricultural production were collected from the Department of Agricultural Extension (DAE). The base maps of different factors were collected from the Local Government and Engineering Division (LGED) and the Bangladesh Country Almanac (BCA) to carry out this research work. The research was conducted in two phases. In the first phase, a land suitability analysis was conducted to identify the areas suitable for rice production ( Figure 3). A land suitability map was prepared in ArcGIS 10.3 ® based on the rice production requirement criteria. An MCDA technique was adopted to carry out the suitability analysis. AHP and MCDA were used with eight factors and one constraint to incorporate the preferences of the experts and assign weights to the factors. On the other hand, GIS was applied to manage the thematic spatial dataset of the factors and constraints' layer maps using spatial analysis tools. A decision rule was applied to integrate GIS and AHP to develop a suitability map. In the second phase, the insurance premium model was developed based on the land suitability classifications. The following sections discuss land suitability using multicriteria analysis and insurance premiums. The research was conducted in two phases. In the first phase, a land suitability analysis was conducted to identify the areas suitable for rice production ( Figure 3). A land suitability map was prepared in ArcGIS 10.3 ® based on the rice production requirement criteria. An MCDA technique was adopted to carry out the suitability analysis. AHP and MCDA were used with eight factors and one constraint to incorporate the preferences of the experts and assign weights to the factors. On the other hand, GIS was applied to manage the thematic spatial dataset of the factors and constraints' layer maps using spatial analysis tools. A decision rule was applied to integrate GIS and AHP to develop a suitability map. In the second phase, the insurance premium model was developed based on the land Sustainability 2018, 10, 1759 5 of 28 suitability classifications. The following sections discuss land suitability using multicriteria analysis and insurance premiums.

Land Suitability with Multicriteria Analysis
This study applied MCDA in several steps in the spatial environment. The following nine factors were selected to conduct the MCDA: (1) slope; (2) elevation; (3) land type; (4) topsoil texture; (5) soil pH; (6) flood prone; (7) rainfall; (8) temperature; and (9) land use ( Figure 4). Only the land use factor was considered as a restriction on the study. A factor is a criterion that enhances or detracts from suitable alternatives for the activity under consideration, and a constraint serves to limit any alternative [36].

Land Suitability with Multicriteria Analysis
This study applied MCDA in several steps in the spatial environment. The following nine factors were selected to conduct the MCDA: (1) slope; (2) elevation; (3) land type; (4) topsoil texture; (5) soil pH; (6) flood prone; (7) rainfall; (8) temperature; and (9) land use ( Figure 4). Only the land use factor was considered as a restriction on the study. A factor is a criterion that enhances or detracts from suitable alternatives for the activity under consideration, and a constraint serves to limit any alternative [36].

Slope
The slope refers to the topography of the area or indicators of the geomorphological properties of the soil. The thickness of the soil layer decreases with the increase of the slope [37]. The original Shuttle Radar Topography Mission (SRTM) and digital elevation models (DEM) were used to produce the slope and elevation for the study area. ArcGIS was used to rectify the topographical maps using the Universal Transverse Mercator (UTM) projection and WGS 84 datum. The datasets prepared in this process were compared to the SRTM dataset to identify the accuracy of both datasets. The slope was calculated from the maximum rate of change between each cell and its neighbors. Every cell in the output raster had a slope value. A lower slope value indicates a flatter terrain and a higher slope value indicates a steeper terrain. Flat fields had a smooth surface, which was better for rice cultivation as it facilitates the even and equal distribution of water. Ninety-five percent of the study area has less than 5% slope (Table 1); almost all the study areas are suitable for rice production in the context of slope.

Slope
The slope refers to the topography of the area or indicators of the geomorphological properties of the soil. The thickness of the soil layer decreases with the increase of the slope [37]. The original Shuttle Radar Topography Mission (SRTM) and digital elevation models (DEM) were used to produce the slope and elevation for the study area. ArcGIS was used to rectify the topographical maps using the Universal Transverse Mercator (UTM) projection and WGS 84 datum. The datasets prepared in this process were compared to the SRTM dataset to identify the accuracy of both datasets. The slope was calculated from the maximum rate of change between each cell and its neighbors. Every cell in the output raster had a slope value. A lower slope value indicates a flatter terrain and a higher slope value indicates a steeper terrain. Flat fields had a smooth surface, which was better for rice cultivation as it facilitates the even and equal distribution of water. Ninety-five percent of the study area has less than 5% slope (Table 1); almost all the study areas are suitable for rice production in the context of slope.

Elevation
Elevation is an important factor that plays a role in the variation of plant cover, causing temperature changes, particularly in highlands. Rice can grow in extremely diverse elevations, including fields that are leveled, gently rolling, or steep, at altitudes of up to 2000 m [38]. However, when the elevation is less than 375 m, rice production is better [39]. Most of the study areas are plains land and the elevation is less than 70 m. Considering elevation, all the lands of the region are suitable for rice production ( Table 1).

Top Soil Texture
Rice grain yield is significantly affected by topsoil (bulk density, water holding capacity, permeability) [40]. The rice grain yield in clay soil was 46% higher than in sandy loam soil averaged across cultivars and water areas [41]. The effective soil depth is defined as the thickness of soil above a layer restricting root growth (e.g., consolidated rock or cemented materials, such as gravel). Most annual crops had a rooting depth of approximately 50 cm, while for tree crops, the rooting system can reach beyond 150 cm. However, most crops produce good yields in soils with an effective soil depth of approximately 100 cm, and this value has been used as an upper limit. There were seven categories of topsoil format available in the area: predominant clay, predominant silt clay, predominant silt clay loam, predominant silt loam, predominant clay loam, predominant silt loam, and predominant sandy loam (Table 1). These categories were converted into a land suitability class according to their characteristics for rice cultivation based on the United States Department of Agriculture (USDA) soil texture suitability rating for rice and are observed as sands, loamy sands, sandy loams, sandy clay loam, and silt (Table 2). Table 2. USDA topsoil texture suitability rating for rice.

Land Type
Good water retention capacity and a high amount of clay with organic matter are considered ideal for rice cultivation. Highlands are used to grow broadcast Aman and Aus rice, medium highlands are suitable for broadcast Aus and Aman, and medium lowlands are used to grow transplanted Aman and mixed broadcast Aus and transplanted Aman. Broadcast Aman is typically grown in lowland areas [42]. Based on seasonal flooding, the government of Bangladesh (GoB) has divided the land into five categories: highland, medium highland, medium lowland, lowland, and very lowland [43]. Most of the area is medium lowland (54%), followed by highland (28%), medium highland (12%), and lowland and very lowland (<6%). Medium highland and lowland are considered to be highly suitable, highland to be moderately suitable, lowland to be marginally suitable, and very lowland is considered to be presently not suitable for rice cultivation (Table 1).

Soil pH
The pH of the soil is defined as the negative logarithm of the hydrogen ion concentration of the soil solution. The pH is an important factor in the quality assessment for paddy soils [44]. Slightly acidic soils having a pH value of 6 to 7 are better for rice cultivation. However, rice can grow in a wide range of pH varying from values of 4 to 8 [45][46][47]. Considering the availability of the data format, the study ranked pH 5.5-7.3 as highly suitable (79%), pH 7.3-8.4 as moderately suitability (8%), and pH 4.5-5.5 as marginally suitable (<3%) ( Table 1).

Temperature
Temperature and rainfall are two climatic factors that have a favorable influence and, in some cases, unfavorable influence on the development, growth, and yield of rice. Being a tropical and subtropical crop, rice is normally grown at a fairly high temperature-high rainfall regime ranging from 20 to 40 • C. The capacity of soil is the amount of water held in the soil between tensions corresponding to field capacity and the permanent wilting point [45,47]. Although the research area had three administrative districts, the temperature was almost similar in all of the areas. The average temperature was 23.5 • C (Table 1).

Rainfall
Rainfall, in the critical stages of paddy growth, increases crop yield through the rapid dissolution of nutrients for uptake by plants [48]. Three districts in the study area receive 1250-2000 mm of annual rainfall. The annual average rainfall was recorded between 1932 and 2931 mm (Table 1).

Flood Prone
Bangladesh ranks as the sixth most flood-prone country in the world [49]. Flood-prone environments include deep water areas submerged under more than 100 cm of water from 10 days to a few months, and areas that are affected by flash floods of longer than 10 days [43,50]. Due to the availability of the data on flood prone areas, the area was categorized into five types: not flood prone (54%), low river flooding (16.6%), moderate river flooding (8.6%), moderate tidal surge (0.13%), and severe river flooding (20%). According to the depth and duration of water stagnation, not flood prone is considered as highly suitable, low river flooding is considered as moderately suitable, moderate river flooding is considered as marginally suitable, and moderate tidal surge and severe river flooding are considered as not suitable for rice farming [51] (Table 1).

Land Use
Land use is an important factor to identify the suitable sites for any crop production. In this study, we considered rivers, forests, waterbodies, and settlements as land use and they are considered as restrictions in the analysis. Approximately 33% of the total lands were occupied by settlements, rivers, and forests. As such, only 67% of the land was considered as suitable for rice production (Table 1).

Preferences of Criteria in Decision Analysis
The preferences of criteria can be referred by the weights. The purpose of weighting is to express the importance or preference of each criteria related to other criteria on crop yield and growth rate. A literature review of various references and interviews with local agronomists and researchers in this field were selected to identify the critical requirements for suitable rice growing areas. Factors were identified for selecting suitable sites for rice. The suitability levels for each of the factors were defined; these levels were used as a base to construct the maps (one for each factor). The suitability levels were: highly suitable-S1, moderately suitable-S2, marginally suitable-S3, and not suitable-N, based on the structure of FAO land suitability classification [52].
The procedure of the MCDA included several stages. First, the relevant factors and constraints were established [53]. Using the abovementioned factors, a pairwise comparison matrix was constructed. Although there are a variety of existing procedures for the development of weights, a pairwise comparison matrix in the context of a decision-making process known as the AHP was used [30]. The comparison concerns the relative importance of the two factors involved in determining the suitability of the stated objective [53]. In a land suitability analysis, a map represents each evaluation criteria with ordinal values (such as S1, S2, S3, and N), indicating the degree of suitability with respect to a sub-criterion and based on the crop requirements [54].

Analytical Hierarchy Process (AHP)
The AHP is one of the most widely accepted procedures and is considered as the most reliable multicriteria decision-making technique [55]. The process was applied to a set of factors to establish a hierarchical structure by providing a weight for each factor in the complete decision-making process [56]. Therefore, many decision-making methods attempt to determine the relative importance, or weight, of the alternatives in terms of each factor involved in each decision-making problem. The AHP provides a structural ground for quantifying the strong comparison of design factors in a pairwise technique, and thus decreases the complexity of the decision-making process [57]. Weights were used to determine the priorities of factors (elevation, slope, land type, soil pH, topsoil textures, flood prone, temperature, and rainfall) and to identify the suitability of different land uses for rice production. The resultant AHP weights were used to determine the priority of each factor for weighted overlay applications using GIS.
In the first stage of the analysis, we organized the factors of the decision model into a hierarchy for land suitability and an insurance premium ( Figure 3). The second stage involved scoring the factors via pairwise comparisons and scoring scales of relative importance (Table 3). Questionnaires were used to gather expert opinions on the relative importance of the factors considered (Appendix A). The comparative results (for each factor pair) were described as integer values of 1 (equal value) to 9 (extremely different), where a higher number denotes that the chosen factor was more important than the other factor to which it was compared. A rating of 9 indicates that in relation to the column factor, the row factor is more important. On the other hand, a rating of 1/9 indicates that relative to the column factor, the row factor is less important [54]. In cases where the column and row factors are equally important, they have a rating value of 1. For example, when comparing elevation and slope angle factors, a score of 1 indicates that both were equally relevant to evaluate suitability, and a score of 9 indicates that elevation is more important than the slope angle. All scores were assembled in a pairwise comparison matrix with diagonal and reciprocal scores located in the lower left-hand triangle. Reciprocal values (1/3, 1/5, 1/7, and 1/9) were used where the row factor was found to be less important than the column factor (Table 3). Table 3. AHP pairwise comparison scale (Saaty, 1990).

Equally important 3
Moderately important to one activity over another 5 Strongly important to one activity over another 7 Very strong important to one activity over another 9 Extremely important to one activity over another 2, 4, 6, 8 When compromise is needed Intermediate values Reciprocals Values for inverse comparison Third, we calculated the matrix and ensured the consistency of the pairwise comparison factors. The AHP also provided the measurements to calculate the normalized values of each factor and to determine the normalized principal eigenvalue and priority vectors. The pairwise matrix was calculated and is given by the following expression: The sum of each column of the pairwise matrix was denoted as follows: We then divided each element of the matrix by its column total to generate a normalized pairwise matrix as follows: Finally, we divided the sum of the normalized matrix column by the number of factors used (n) to generate the weighted matrix of priority factors as follows: The initial consistency vectors were derived by multiplying the pairwise matrix by the vector of weights as follows: The principal eigenvector (λ max ) was then calculated by averaging the values of the consistency vector as follows: Eigenvalues were calculated by averaging the rows of each matrix. Eigenvalues were also referred to as relative weights. The largest eigenvalue was equal to the number of factors, and when λ max = n, judgments were consistent. Normalized eigenvalues were generated as weights of priority factors. The principle value suggests that the eight criteria were consistent, as the calculation results reveal a maximum value of 8.34 (Table 4). The judgments were also checked to determine the consistency index (CI), which was calculated as follows:  Here, n is the total number of criteria. Saaty introduced the consistency ratio (CR) and compared it to the CI and random index (RI) values, which is the calculated value for matrices of different sizes [58] ( Table 5). The consistency ratio was calculated as follows: A lower CR ratio indicates more consistency. If the value of CR is 0.10 or more, then the weight values of the matrix indicate inconsistencies and the method (AHP) may not provide a meaningful result [30].

Weighted Overlay Analysis for Land Suitability
The rice crop suitability map was identified by the weighted sum overlay using the spatial analyst tools in ArcGIS 10.3 ® with both vector and raster layers ( Figure 5). Initially, nine thematic vector layers of slope, elevation, soil pH, topsoil texture, land type, temperature, rainfall, flood prone, and land use (rivers, water bodies, settlements, and forests) were taken in a base geographical coordinate system, WGS 1984 (Table 6). These vector layers were projected in WGS 1984 UTM Zone 45N to obtain the same geographic extent. The thematic layers were then converted to raster layers to conduct the spatial analysis. The AHP weights were obtained from the experts' opinions using a designed questionnaire for the priority weights of the factors. The constraints were stated as "restrictions". To aggregate the weights of the factors and constraints, a decision rule, WLC, was applied to determine suitable locations in the GIS environments.
The rice crop suitability map was identified by the weighted sum overlay using the spatial analyst tools in ArcGIS 10.3 ® with both vector and raster layers ( Figure 5). Initially, nine thematic vector layers of slope, elevation, soil pH, topsoil texture, land type, temperature, rainfall, flood prone, and land use (rivers, water bodies, settlements, and forests) were taken in a base geographical coordinate system, WGS 1984 (Table 6). These vector layers were projected in WGS 1984 UTM Zone 45N to obtain the same geographic extent. The thematic layers were then converted to raster layers to conduct the spatial analysis. The AHP weights were obtained from the experts' opinions using a designed questionnaire for the priority weights of the factors. The constraints were stated as "restrictions". To aggregate the weights of the factors and constraints, a decision rule, WLC, was applied to determine suitable locations in the GIS environments.

Crop Insurance
Farmers face a variety of substantial risks in the production and marketing phases, yet agricultural insurance has rarely been available outside of subsidized government offerings, which typically reach only a small fraction of small-scale farmers. A large number of newly developed products were piloted in developing countries including Bangladesh, however, none were perfect considering the socioeconomic characteristics of the farmers in Bangladesh [59]. In this study, the insurance premium mainly considers the quality of the land or suitability of the land. An AHP survey was first conducted to identify the criteria for the insurance premium. After literature reviews and a pilot survey, we identified some major factors that could be highly influential to take decisions about agricultural insurance. In different perspectives, we tried to identify the factors and causes of the agricultural insurance market. Subsequent to the careful analysis of all the factors, eight major factors were selected, which were highly influential in different stakeholder's decisions. The major factors identified were product price, government subsidy, awareness, willingness of the government, belief of the benefits, trust in the company, rules and regulations, and the collection and distribution of money. To determine the crop insurance premium, the factors were ranked by the experts and the influence in the crop insurance program was obtained; most of the experts mentioned product price as the main factor for participating in a crop insurance program (Table 7). Finally, the study proposed an insurance premium model based on land suitability with the consideration of the harvest index and subsidy dependence factor to increase rice production domestically. The product price depends on the production, which was influenced by the harvest index (HI). The HI can be calculated from the crop information of past years based on the biological yield, and yield of dry matter or yield of grain and straw [60,61]; the subsidy dependence factor was considered because it is difficult to find any crop insurance that was successful without a government subsidy. Table 7. The judgements of experts in prioritizing the weights of the different factors for introducing crop insurance premium model (A-J, experts from different fields). Most of the farmers in Bangladesh have limited incomes and consider insurance premiums as another cost of production. During the 1980s and 1990s, agricultural insurance was exclusively provided through the state-owned insurance company, Shadharan Bima Corporation (SBC). SBC offered an individual grower multiple peril crop insurance product as well as livestock mortality and aquaculture insurance. Insurance demand is moderately price elastic [62], while crop insurance demand varies across the household head's primary occupation, land ownership, and the size of the farm land [63].

Awareness
Farmers had limited awareness of agricultural insurance. In most cases, no field-level practice of crop insurance is the main cause for this. The lack of knowledge and awareness in the design, rating, and implementation of agricultural insurance need to be considered.

Belief of the Benefits
Almost all the farmers had some bad experiences, not directly owing to the insurance company, but some other microfinance organizations. Due to the bad experiences, farmers do not have sufficient confidence for a payment for any loss.

Collection and Distribution of Money
There are some fixed costs and high administrative costs for any operation, as rural infrastructure is not well developed and most of the farms are small in size. The sizes of many farms were less than one hectare, which causes the costs of insurance delivery, underwriting, and claims administration to be prohibitively high [64]. The main constraints and opportunities for crop and livestock insurance in developing countries are considered, such as insurance product types, hazards, vulnerability, and rural institutions that can support organization and distribution [65].

Trust in the Company
There is no clear policy framework for agricultural insurance in Bangladesh, including a lack of clarity on the role of government in supporting agricultural insurance through the private insurance sector. The commercial insurance companies are not willing to take a lead in developing crop insurance products.

Rules and Regulations
Establishing a business and continuing it in an appropriate manner needs some appropriate rules and regulations. How a company can work, what are the terms and conditions, and how a farmer participates in a crop insurance policy and gains benefit should be clear and understandable. While a premium subsidy is the most common intervention, other enabling measures are important, such as the legal and regulatory framework, reinsurance, technical and administrative assistance, and linkages to government extension services in agriculture, animal health, or meteorology [65].

Willingness of the Government
In developing countries, most of the decisions are made through government intervention. Agricultural insurance is normally undertaken as a market-based activity by private or state sector insurance companies, often with support measures from the government. In high and some middle-income countries, such as the United States, Canada, Spain, Italy, Japan, and South Korea, governments have significantly intervened, created sophisticated market structures, and accepted significant annual budgets for agricultural insurance [65].

Government Subsidy
Private commercial insurance companies in Bangladesh have limited financial capacity, and are, in general, reluctant to take a lead in investing in agricultural insurance staff and the design of products and policies, systems, and procedures, as this is a high-risk class of insurance. Thus, introducing crop insurance in a vast area that requires a good budget from the government for a subsidy. The estimated discrepancies between an expected indemnity and expected premium appear too large to be financed by government subsidy on a continuous basis [63].

Insurance Coverage and Premium Rate
Agricultural insurance is generally undertaken as a market-based activity by private or state sector insurance companies, often with support measures from the government [65]. The demand for crop insurance is said to be moderately price elastic [62,66]. Government-subsidized crop insurance programs have been established and developed for decades in developed countries, particularly in the United States and Canada [67]. In many developing countries including Bangladesh, crop insurance was introduced approximately three decades ago. Some of these countries are continuing to offer such insurance, while in some others, it has stopped functioning owing to heavy losses. Insurance in the United States is provided by 17 private companies. They work in agreement with the USDA Risk Management Agency and approximately 45% of the field crops' production value are insured (23% in the European Union (EU)). The average premium rate is close to 9%, which is much higher than in Europe (4%), mainly because they offer a wider coverage-revenue or yield insurances versus mainly single-peril or combined-risk insurances. The premium subsidies amount to $1900 million (58% of the total premiums). The United States government also provides funds for the administrative costs of insurance companies and provides reinsurance. Thus, the total support provided to insurance amounts to 72% of the total premiums (in the EU, approximately 500 M€ = 32% support) [68].
In this research, the insurance premium mainly considers the quality of land or suitability of the land to cover the crop revenue of farmers. An AHP survey was conducted to first identify the criteria involved for insurance premiums (Appendix B). To determine the premium rate, different low interest rates were considered that could have the ability or "willingness to pay." The range was 2%-6% of the insured revenue, which was considered as the willingness to pay the premium for their crops. The premium rate varies due to the land quality and risks [52]. The expected revenue was calculated with the help of the regional average production of time series data and the government floor price.
The insurance coverage was categorized as high coverage (90%), moderate coverage (80%), and marginal coverage (70%) considering the average production costs of the production of rice [69]. The minimum insurance coverage was considered marginal coverage (70%) to overcome the production costs, while the maximum coverage was estimated as high coverage (90%) to ensure the production efficiency of the factors of productions. Insurance coverage was the same for different land suitability, harvest index, and subsidy dependence factors, however, the premiums were different. Thus, it could benefit land owners with moderately suitable or marginally suitable lands. No insurance was covered for the presently not suitable areas due to the physical and economic viability. A farmer could buy a product of insurance from the three categories and could pay the premium in 12 installments (every month). The future value of annuity (V) was calculated using regional yield and the price of rice varieties as follows [70]: where P is the premium for different coverages and can be expressed by combining the harvest index and subsidy dependence index. The future expected revenue from rice production can be expressed as follows: where V j is the cumulative return from the yield for different coverages (j) of three seasonal of rice practices in one calendar year. The PPV is the average production capacity for seasonal rice (ton/ha; i is the seasonal variety: Aus, Aman, and Boro; and n is the number of years) and PRT is the average price of rice ($/ton). Therefore, the premium for different coverages (j) are based on the expected revenue, production function, and subsidy which can be expressed as follows: where P j is the premium for the specific coverage (j) of three seasonal rice practices in one calendar year and HI is the harvest index and ratio between the observed yield and biological yield. The SDF is the subsidy dependence factor and was assumed to be 0.50 as a welfare consideration. The parameter C f is the conversion factor for the pixel resolution for Landsat 8, OLI. The government of Bangladesh provides agricultural subsidies not only in the form of cash but also on price of agricultural inputs, such as seeds, fertilizer, irrigation, and fuel. Among production inputs, the cost of irrigation holds the highest share (nearly 20%) as a single cost item followed by fertilizer (approximately 18%) because of the price increase of these items at both the national and international markets [71]. Electricity in agriculture enjoyed a 20% subsidy for a long time, however, it was discontinued later to offset the impact of a rise in oil prices. Presently, the policy is proposed for at least a 25% subsidy on electricity and fuel needed to continue to meet future challenges, particularly, the Boro crop for which the irrigation is powered 40% by electricity and 60% by fossil fuel [72]. If the price of electricity and fuel is within the capacity of farmers, they will not lose interest in Boro cultivation. The subsidy expenditure as a total subsidy is presently offered at 17.30%, the proposed budget (2017-2018) by the government is 22.13%, and electricity or fuel for irrigation is proposed to increase from 20% to 25% [73].
The parameter C j is the coverage rate (highly, moderately, and marginally) and C f is the conversion factor to scale down the pixel level (30 m × 30 m, Landsat 8, OLI).

Land Suitability
We have found that highly suitable land (S1) covers an area of 1337.02 sq. km (22.74%) of the total land. Hence, a significant part of the study area was highly suitable for rice cultivation ( Table 8). The S1 (highly suitable) areas were characterized using a slope level of 0%-10%, elevation less than 50 m, soil pH level between 5.5 and 7.3, the topsoil texture class that included predominant clay, predominant silt clay, and predominant clay loam, the land type of medium high and medium lowland, a rainfall range between 1000 and 2000 mm, and temperatures between 20 and 30 • C. The result showed that 1678.04 sq. km (28.54%) of the area were moderately suitable (S2) for rice production, and 873.71 sq. km (14.86%) of the area were marginally suitable (S3) for the cultivation of rice with limitations of land type, soil texture, and the soil pH. This class was used less for rice production as its limitations and profitably corrected the suitability class to a moderate level. Currently, not suitable areas (N) had also covered a relatively insignificant portion of the study area of 69.97 sq. km. (1.19%). Generally, not suitable areas (N) were in very lowland areas with a sandy loam texture. In addition, the total lands that were permanently not suitable for rice production were 1920.87 sq. km (32.67%), occupied by settlements, rivers, forests, and water bodies. The results show that the majority (up to 67%) of the total land in the study area had varying degrees of suitability for rice production ( Figure 5). The highly and moderately suitable areas, taken altogether, which have, according to the analysis, covered 51% of the total land was commendable for rice production. Overall, these results indicate that the study area has enormous potential for rice production. Therefore, economic levels of agricultural production can be achieved by cultivating rice crops in highly (S1) and moderately (S2) suitable areas and practicing the diversification of marginally (S3) suitable areas to crops other than rice.

Insurance Premium and Coverage
The insurance coverages were categorized as high coverage (90%), moderate coverage (80%), and marginal coverage (70%) ( Table 9). Insurance coverage was the same for different land suitability, but the premiums were different. Thus, it could benefit land owners with moderately suitable or marginally suitable lands. Presently, not suitable areas are not recommended for insurance coverage because of cost effectiveness.
A farmer could buy a product of insurance from the three categories and could pay the premium in 12 installments of three rice crops (Aus, Aman, Boro) during crop plantation, in the middle of cultivation, and after harvest. The average weights for each factor were derived using individual judgments. Considering the fact that farmers want to produce rice with secured and minimum premiums, the experts' opinions were listed to recognize the criteria to introduce an insurance premium according to land suitability. Product price was the main factor (29%) to participate in crop insurance, which was highly involved with the production factor as mentioned in the harvest index. Farmers had limited awareness of agricultural insurance. There were 17% of the experts who mentioned that there was a low level of awareness and knowledge among Bangladeshi farmers on the role of agricultural insurance. Although the villagers had no understanding of crop insurance, when they understood the theme, they were rather worried (9%) about the benefits. In the case of collection and distribution of money, the collection of the premium and distribution of the insurer's money had a large involvement in the cost issue (3%). Trust in the company was also a challenge to prove (5%). In the rules and regulations criteria, the lack of a national framework (4.5%) for agricultural insurance and limited agricultural insurance provisions were mentioned by the experts. The willingness of the government was the second (18%) most important criteria for the establishment of the crop insurance market. In the case of the government subsidy, private commercial insurance companies in Bangladesh have limited financial capacity and are generally reluctant to take a lead in investing in agricultural insurance. Therefore, introducing crop insurance in the poverty and risk-prone areas requires the government to provide a subsidy. In the survey, we found that 16% of the experts' opinions consider a government subsidy as an important factor for farmer participation in crop insurance.

Discussion
In this research, the suitability mapping of rice production for the poverty prone and risk regions of Bangladesh was conducted to increase production and to ensure food security. From the expert analysis information, a key challenge for poor farmers of the study area was the management of the activities to reduce crop damage risk. One potential strategy for coping with the devastating effects of crop damage is to identify the suitability areas and introduce a crop insurance scheme on a cooperative basis. Agricultural insurance along with the cooperation of the community could help farmers to stabilize farm income and investment and guard against disastrous effects of losses due to natural hazards, low production, or low market prices. Crop insurance not only stabilizes farm incomes but also helps farmers to initiate production activity after a bad agricultural year, shock of crop losses, or low prices by providing farmers with a minimum amount of protection and helping farmers invest more in agriculture.
The socioeconomic identities of farmers, poor rural infrastructure, very low household incomes, large family sizes, low productivity, and insufficient training on agriculture made the regions more risk prone. To achieve sustainable livelihoods, a crop insurance policy needs to be introduced with a cooperative procedure, as most of the farmers do not deal with banks and nongovernmental organizations. In this research, a procedure to calculate premiums was introduced, based on land suitability areas, a harvest index, and a subsidy dependence factor. Considering the above factors, the government of Bangladesh needs to take initiatives by subsidizing this sector initially, which suggests a premium rate of 2%-6% of the expected outcome based on land suitability, harvest index, and subsidy dependence factor. Considering land suitability, there are three types of insurance coverage: high coverage (90%), moderate coverage (80%), and marginal coverage (70%). We suggest producing rice up to marginally suitable land; a low harvest index could gain support through insurance. The minimum coverage was marginal coverage (70%), which covered the production costs, while maximum coverage was designed using high coverage (90%) to utilize the maximum benefit against the risk-prone marginal areas. The expected outcome can be calculated from the time-series production data. The insurance premium introduced in India has offered a monumental and path-breaking bounty of crop insurance at a nominal premium-as low as 1.5% and 2% of the sum assured-which is unprecedented, to boost the farm sector and derisk the vast farming community. Under this new scheme, the government of India has substantially reduced the premium to be paid by the farmers-1.5% on Rabi crops (dry season), 2% on Kharif crops (wet season), and 5% on commercial/horticultural crops. To reduce delays and leakages, the payment of compensation was undertaken through direct transfers into the bank accounts of farmers [74]. Further research needs to be undertaken considering community acceptance to justify the premium as proposed in the research considering land suitability, harvest index, and subsidy in the risk-prone areas of Bangladesh.

Conclusions
To meet the demand of food for the increasing population with decreasing land resources, land suitability identification for crop production is very important. The study conducted an LSA analysis on the northern part of Bangladesh and found that approximately 33% of the total land, which was occupied by water bodies, rivers, forests, and settlements, was permanently not suitable, and 1.19% was presently not suitable. However, 22.74% of the study area was highly suitable, 28.54% was moderately suitable, and 14.86% was marginally suitable for rice production. Considering the demand for rice by the vast population, the regions need to secure rice production up to the marginally suitable lands. In the marginally suitable lands, the risks are higher than the moderately suitable areas, and the costs of production are also higher; the production could be lower than the moderately suitable lands. To motivate owners of low-quality land to produce rice, protection through crop insurance support must be provided. This new crop insurance model, "crop insurance based on land suitability", can be a very effective tool to manage production risks and provide rational support to owners of land of different quality, which will increase production to ensure food security. The socioeconomic characteristics of farmers in Bangladesh are far different from those in developed countries, and if we merely follow the latter's products and policies, rice production may not be sustainable. This study suggests a premium rate of 2%-6% of the expected outcome based on land suitability. The premium calculation that is based on costs and risks does not provide a suitable environment for the introduction of crop insurance. The strategy of the government of Bangladesh to provide subsidies to the farmers at times of shocks needs to be changed to subsidies to pay insurance premiums with appropriate regulations.
We suggest producing rice up to marginally suitable land, and this can be supported by insurance. The minimum coverage is marginal coverage (70%) to cover the production costs, while maximum coverage is ensured through high coverage (90%) to enable maximum efficiency. Based on these judgments, the highest priority was the product price (26.72%), followed by willingness of the government (20.27%), government subsidy (19.18%), awareness (11.92%), belief of the benefits (7.12%), rules and regulations (6.70), trust in the company (4.76%), and collection and distribution of money (3.34%). These results are helpful for policymakers and agricultural extension workers who assist the agriculture sectors to increase the production and ensure the welfare of the farmers. The information provided, at a local level, could be used by the growers to select suitable lands and can be relevant to other crops to increase production. The government could encourage private participation to provide varied products with premium rates for commercial farmers on an individual basis. It would be wise to integrate microinsurance with microfinance in order to manage the risk of default and enterprise losses in a more proactive manner.
Author Contributions: All the authors made substantial contributions to this research. M.M.I., the first author and PhD candidate, collected data and analyzed the data with satellite remote sensing datasets in GIS, AHP environment. In addition, he completed this manuscript as a part of his requirement for his doctoral dissertation. T.A. developed the original research concept and supervised the research to develop the land suitability and geospatial insurance premium. R.N. contributed to the AHP analysis and provided scholarly comments to develop this research. All the authors consulted one another while preparing the manuscript at all stages, interpreted the results, and agreed about the conclusions and further research directions.

Example: In taking the decision to develop the crop insurance how important this criterion:
Instruction: Please add a checklist (✓) on the column scale criteria (A) or at the column scale criteria (B) based on experts' opinion.

Definition of scale:
1: the second equally important criteria 3: criterion (A) is little more important than (B) 5: criteria (A) is more important than (B) 7: criterion (A) is much more important than (B) 9: criterion (A) is absolutely more important than (B) * and that goes for the reverse Example: In taking the decision to develop the crop insurance how important this criterion: If you give a sign checklist (✓) on a scale of 7 in column A, then it means that the criteria An in this example "Land Type" is more important than criteria B in this example "Topsoil Texture".
If you give a sign checklist () on a scale of 7 in column A, then it means that the criteria An in this example "Land Type" is more important than criteria B in this example "Topsoil Texture".