A Sustainable Land Utilization Pattern for Confirming Integrity of Economic and Ecological Objectives under Uncertainties

In this study, an integrated crop–forest system with market approach (ICFM) associated with recovering forest and withdrawing cultivation was developed for confirming regional integrity of economic and ecological objectives under uncertainties. A mixed quadratic stochastic-fuzzy programming method (QSF) was proposed for planning an ICFM issue under uncertainties. QSF can not only deal with spatial and temporal variations of meteorological condition, but also handle uncertainties expressed in terms of probability distributions and fuzzy sets. Meanwhile, it can also tackle nonlinear relationships between land resource plan and economic data. The developed QSF was applied to an ICFM issue in Xixian county, China. The results of adverse effects from irrigation, ecological effects from forest, land utilization with market approach and optimal system benefits were obtained. It can facilitate policymakers to adjust current land utilization with market approach to improve the productivities of land resources. The tradeoff between crop irrigation and forest protection can prompt generation of optimized plans with consideration of economic and ecological objectives, which can be availed to generate strategies for confirming integrity of socio-economic and eco-environmental sustainability.


Introduction
Over the past half century, the high-speed development of agriculture has greatly contributed to food safety, which could improve human lives.More land resources such as forest have been over-exploited into irrigated land to satisfy the accelerated demand for food production; meanwhile, woods have been felled for urban construction.However, over reclamation and cultivation could disturb the balance of groundcover; meanwhile, over-fertilization and irrational planting schedule have brought great stresses to water and environment, leading to land degradation [1,2].In recent years, land degradation has caused frequent detrimental phenomena such as sand storms, detritus streams, water deficits and water pollutions, which have challenged policymakers worldwide.To regain the frail ecological function of land resources, an integrated crop-forest system (ICF) project associated with recovering forest and withdrawing cultivation has been advocated in the context of rapid development of agriculture.However, the tradeoff between irrigative development and forest protection has aggravated the difficulty of ICF project promotion.On the one hand, ICF can be beneficial to regain ecological services from land (e.g., erosion control, recreation opportunity, soil and water conservation), leading to more intangible ecological returns to human beings [3,4].On the other hand, in the short run, an ICF project may reduce the direct benefit of economy due to withdrawing reclamation, which results in it being disputed.Therefore, an effective method such as market approach can be joined into a land resources reallocation issue to improve economic productivity of land usage [5].In an ICF with market approach (ICFM), land users may face the higher penalties or losses of land degradation; thus, they would consider withdrawing cultivation into forest land based on being given maximum environmental allowance, with aim to reduce these penalties or losses.Through a market approach, the land resources can be encouraged from lower value (or higher loss) to higher one (or lower loss), which can improve the productivity of land resources.Meanwhile, the land users may regulate their fertilizations and planting schedules to reduce the penalties for pollution discharges or losses of land degradation.It can also facilitate policymakers to adjust current irrigative structure and land utilization pattern in a more sustainable manner with consideration of economic and ecological objectives.
However, a practical ICFM issue is often complicated with numbers of uncertainties associated with economic development, environmental protection, marketable regulation, effective governance, and technical factors [6][7][8][9].For example, regional meteorological condition can impact ecological function of land resources, which can be deemed as spatial and temporal variations (e.g., water flows for plantation) existing in system components.Meanwhile, intricate factors such as market discipline and social-economic development would influence the processes of land planning, leading to vagueness [10].Moreover, the interactions between uncertain system components such as recourse actions between expected target and available land resource can be embodied into an implication of an ICFM issue, which would increase the complexity of ICFM system [11].The above uncertainties and their interactions require policymakers to obtain a comprehensive, complex and ambitious plan.
Previously, many researchers worked on planning land resources with consideration of uncertainties associated with precise environmental load, stochastic water flow, fuzzy ecological effect, complex economic development and uncertain market regulation [12][13][14][15][16][17].For instance, Perez-Garcia et al. [12] used Global Trade Model (GTM) and the Terrestrial Ecosystem Model (TEM) to estimate the potential effects of climate change on the scale of global forest land, which can provide ecological-economic aspects with uncertain economic, ecological and climate scenarios for policymakers.Hauk et al. [13] qualified the economic risk of short rotation woody crops (SRWC) to evaluate the economic diversification effects through Modern Portfolio theory, which can help to find an effective risk-management instrument and profitable crop in a farm scale.However, the above risk analysis methods have difficulty in reflecting the influence of random events and their corresponding effects.Therefore, a stochastic method is introduced to handle spatial and temporal uncertainties [14].For example, Whelan et al. [15] described a stochastic phosphorus (P) transferring model in an agricultural land and water resources issue, with aim to quantify the effects of spatial variability in an effective manner.Han et al. [16] planned land and water resources through an interval-parameter linear optimization model with stochastic vertices, which can arrange irrigative schedule against risk in the Hetao irrigation region.Djanibekov and Khamzina [17] concentrated on the economic viability of afforestation on marginal irrigated croplands with stochastic method in irrigated drylands.Djanibekov and Villamor [18] developed a farm-level dynamic mean-variance model combined with a real options approach to identify the relationship between agroforest conservation with market-based instruments (MBIs) and risk aversion of farmers.However, the above stochastic methods cannot reflect the recourse actions between expected target and random events (such as spatio-temporal water flow), which can be introduced to stochastic programming (SP) to handle this issue in a practical land planning issue.For instance, Guo et al. [19] developed an interval fuzzy-stochastic programming for irrigation planning within an agriculture system under multiple uncertainties.Zeng et al. [20] proposed a two-stage stochastic programming to reflect recourse actions between expected irrigative targets and random variables, which is effective to deal with stochastic uncertainties expressed as probability distributions.However, SP has difficulty in handling information with vagueness and fuzziness, which requires fuzzy programming method (FP) to regarded ambiguous parameters [21,22].
For instance, Maqsood et al. [23] mixed fuzzy programming, interval-parameter programming and two-stage stochastic programming into a framework to plan water resources/ecological environment management under uncertainty.Altunkaynak and Sen [24] introduced fuzzy membership functions to evaluate the dynamic valuation of ecosystem services in Lake Van, eastern Turkey.Deng et al. [25] improved entropy-based fuzzy matter-element model to reflect objective fact in river health assessment, which can improve feasibility of assessment.Nevertheless, traditional FP can present confidence degree of event occurrence regarded as possibility distribution, which cannot satisfy the higher requirement of fuzziness expression, particularly in the situation of vagueness superposition.Therefore, credibility constrained programming (CCP) can be jointed into FP to improve the abilities of fuzziness expression through mixed possibility and necessity measures [22].In addition, there are nonlinear relationships existing in land resource and economic parameters, which cannot be tackled by SP, FP, or CCP.Thus, quadratic programming (QP) can be introduced to reflect such nonlinearity to support for generation of global optimization under multiple uncertainties [6,11,20,[26][27][28][29]. Unfortunately, few studies focus on mixing SP, FP, CCP and QP methods for planning land resources.
Therefore, this study developed a sustainable land utilization pattern with market approach (ICFM) for confirming regional integrity of economic and ecological objectives under uncertainties.A mixed quadratic stochastic-fuzzy programming method (QSF) was proposed for an ICFM issue to handle uncertainties and systemic complexities, which can not only handle uncertainties expressed as probability and possibility distributions, but also quantify nonlinearity among economic parameters and land resources.The proposed model was applied to a real case of a crop-forest system in Xixian county, China.Sound land utilization plans associated with reclamation withdraw and forest recovery were designed to satisfy regional economic and ecological objectives.Meanwhile, a market approach was considered for a land utilization plan to increase the economic productivity of land resources.The obtained results can facilitate policymakers to identify optimized irrigated crop and forest protection policies with a sustainable manner, achieving integrity of socio-economy and eco-environmental sustainability

Study Area
Figure 1 presents Xixian county, part of Henan Province located in the upper Huaihe River Basin (113 • 15" E-114 • 46" E, 31 • 31" N-32 • 43" N), with area of 10,191 km 2 .The study region is composed of 65% mountain area and 35% flat depression.Since it is situated in the transition zone between the northern subtropical region and the warm temperate zone, the area is suitable for crop plantation with the average annual precipitation being 645 mm.In the study region, endemic plants such as Taxodium, Ginko and Metasequoia grow well to support regional wood-processing industry.Meanwhile, irrigation for crops such as grains, oil plants, vegetables and fruits has been developed rapidly in recent years, which can improve the quality of human living in the study region.
On the context of population growth and urbanization, forest lands have been exploited as irrigated lands, which could satisfy increasing food demand and accelerated urbanization, but bringing about great water and environment stresses.For example, from 2000 to 2015, irrigated land area has increased 3.56 × 10 3 km 2 , while water deficit of irrigating was 38.13 × 10 6 m 3 at its peak [30][31][32][33].Meanwhile, excessive fertilization can aggravate severe soil loss and pollution emission, leading to the destruction of land function.

Construction of an Integrated Crop-Forest System
To remedy the damage from crop irrigation, an integrated crop-forest system (ICF) can be built (Figure 2).

Construction of an Integrated Crop-Forest System
To remedy the damage from crop irrigation, an integrated crop-forest system (ICF) can be built (Figure 2).

Construction of an Integrated Crop-Forest System
To remedy the damage from crop irrigation, an integrated crop-forest system (ICF) can be built (Figure 2).Framework of land utilization pattern with consideration of ecological effect under uncertainties.In an ICF system, forest can play an important role in regaining ecological function, such as flood control, aquifer replenishment, sediment retention, and water filtration.A proper land utilization plan associated with recovering forest and withdrawing cultivation is required, where the adverse effects from crop irrigation can be transformed by forest system.For instance, excessive nitrogen/phosphorus discharges from excessive fertilization can be purified by the forest system through its purification capacity [17].Meanwhile, economic forest plantation can bring about direct and indirect economic returns.Based on eco-compensation mechanisms, an ICF issue can be considered an effective manner to regain ecological function and achieve regional sustainability in Xixian county.However, limited land resources can restrict the implementation of crop-forest strategy.Therefore, a market-based approach can be introduced to improve the economic productivity of land exchange from low to high values.

Modeling Formulation
A policymaker has the following challenges: (a) A land allocation to crop irrigation can generate economic benefits but can also bring about adverse effects on ecological environment.How to balance the relationship between economy and environment can be deemed as a key point in regional sustainability in study region.(b) Since a land plan in ICF is associated with withdrawing reclamation and recovering forest, the tradeoff between the losses of damages to ecosystem from irrigation and the cost of withdrawing reclamation (including decreased crop planting scale and increased forest recovery) can facilitate policymakers to adjust current land plans.How to identify an optimal strategy for coordination of economic and ecological objectives can be an important issue in ICF.(c) Several factors such as trading cost or reclamation cost can influence the exertion of land reallocation with a market approach in a ICFM issue.(d) Many uncertainties and their interactions existing in a ICFM issue, which could aggravate the complexity of ICFM.How to develop an effective QSF method to identify uncertainties and handle their interactions is a necessary step to generate an irrational result against systemic risk. max In this objective function, Model (1a) presents the maximum of system benefit in a ICFM issue.It equals the benefit from irrigated land and ecological effect minus the loss of water shortage, penalty of water pollution and water/soil loss (US$).BC t is the expected benefits from irrigative activity and forest protection (US$), which is shown in Model (1b).BD t is the losses of water deficit areas for crop irrigation and forest protection (US$), which is detailed in Model (1c).LE t is the losses of water/soil erosion due to irrigative activity (US$) [Model (1d)].PS t is the penalties for environmental pollution [Model (1e)].BW t is the benefits of ecological effect from land plan with market approach, which is presented in Model (1f).In this study, t is the subscript for planning periods, where three planning periods can be considered (t =1 Period 1, t = 2 Period 2, and t = 3 Period 3).Among them, one planning period is three years.Meanwhile, since crop irrigation and forest protection can be incorporated into a framework, three main types of crops [m = 1 grain plant (ha), m = 2 oil plants (ha), and m = 3 vegetables (ha)] and forests [n = 1 economic forest (ha), n = 2 shelter forest (ha), and n = 3 forest park (ha)] can be shown in an ICFM system.Therefore, m and n are the subscripts for main types of crops and forests.
Model (1b) presents expected benefit from a land plan incorporated with a forest system and an irrigation system.Since expected planning areas (or expected land demands) for crops and forests have been regulated in the beginning of planning period, expected benefit equals the net benefit parameter per area being allocated within the expected target area.W M tm and WEA tn are expected land demands for crop and forest area in period t (ha), respectively, which are the first-stage decision variables in an ICFM issue (ha).However, in a practical ICFM issue, there are various complex relationships between economic parameters and land plans, which can be expressed as nonlinearity.Therefore, a mixed quadratic stochastic-fuzzy programming method (QSF) can be introduced to handle nonlinear uncertain relationships between economic data and expected targets quadratically (e.g., (bm tm × W M tm + cm tm ) × W M tm ) (US$), where bm tm and cm tm are the quadratic coefficients for W M tm (US$/ha).The nonlinear relationship between economic data and expected target (W M tm ) can be calculated as a function (e.g., (bm tm × W M tm + cm tm ) × W M tm ).ba tn and ca tn .are quadratic coefficients for WEA tn .Therefore, bm tm , cm tm , ba tn and ca tn are quadratic coefficients for first-stage decision variables (i.e., W M tm and WEA tn ).Similar principles are expressed in Models (1c)-(1f).
Model (1c) shows the losses of water deficit areas (SM tm and SEA tn ) for crop irrigation (grain plant, oil plant and vegetables) and forest protection (economic forest, shelter forest and forest park) (US$).SM tm and SEA tn are second-stage decision variables to rectify the first-stage decision variables through QSF method (as shown in the Appendix A).In general, after the land resources are pre-regulated to each user, if the water has not been delivered to the regulated land, water deficit occurs, leading to economic losses.In the study region, the requirement of water per area equals the planting scale multiplied by the coefficient of water utilization for irrigation and forest (α tm and β tn ).α tm is the coefficient of water utilization for irrigative activity (ton/ha), and β tn is the coefficient of water utilization for forest protection (ton/ha).However, in a practical ICFM issue, water flows are often regarded as uneven levels, which can influence the exertions of economic and ecological functions.The seasonal water deficits can generate adverse effects on the regional development, reducing economic benefit.In Model (1c), p is the subscript for water flow levels, where five water levels (p = 1 low, p = 2 low-medium, p = 3 medium, p = 4 medium-high, and p = 5 high) are considered.dm tm , em tm , da tn and ea tn are the quadratic coefficients for SM tm and SEA tn , which indicate the nonlinear relationship between net loss of water deficit area and expected land demand (US$/ton).Therefore, dm tm , em tm , da tn and ea tn are the quadratic coefficients for second-stage decision variables.
With high runoff, water and soil loss lead to losses (or reduced benefits).Models (1d) and (1e) present loss for water/soil erosion [Model (1d)].In Model (5d), sl tm is the coefficient of soil loss from crop irrigation in period t (ton/ha) and wr tm is the coefficient of water loss from crop irrigation in period t (ton/ha).eLS tm and eLW tm are the quadratic coefficients of reduced benefit for soil and water erosion (US$/ton).
Moreover, with high runoff and soil loss, excess total nitrogen (TN) and total phosphorus (TP) from crop irrigation are discharged into water body directly, leading to pollution.Model (1e) presents the penalty for environmental pollution (US$).EN N tm is the TN content in lost soil due to irrigative activities in period t (%).EP P tm is the TP content in lost soil (%).eLS tm , gLS tm , eLW tm and gLW tm are quadratic coefficients for soil and water erosion, which can present nonlinear relationships between economic loss of water/soil erosion and land reallocation (US$/ton).Similarly, eCA N tm , gCA N tm , eCA P tm and gCA P tm are quadratic coefficients of penalties for TN and TP discharges (US$/ton), which can reflect nonlinearity issues in an ICFM issue.
Since forest systems play an important part of regaining ecological functions (such as flood control, aquifer replenishment, sediment retention, and water filtration), benefit from forest recovery can prompt achieving ecological objectives in the study region.Model (1f) presents the benefit of ecological effect from land plan with market approach.In an ICFM issue, the policymakers (or farmers) have to choose to pay the penalty for environmental pollution or reduce the crop planting scale to reduce pollutant discharge [25,26].In an ICFM issue, the market approach can support land from lower value to higher value, which is adapted to reclamation to reduce damage to environment.WRA tm and WRE tn are reallocated land resources for crop and forest area with a market approach, where the reallocated land resources can generate benefit due to withdrawing reclamation (ha).DA N tn , DA p tn , SI A tn and W I A tn are reduced amount of water/soil erosion and pollution from reclamation (ton/ha).Meanwhile, the reduced irrigative land can be transformed as forest land, which can generate indirect benefits due to ecological effect.DE N tn , DE p tn , SIE tn and W IE tn are increasing ecological effects due to forest recoveries (ton/ha).eBA tn , gBA tn , eBE tn and gBE tn are quadratic coefficients for reallocated land resources (WRA tm and WRE tn ), which can reflect the nonlinear relationships between benefit from ecological effects (US$/ton) and reallocated land resources.
In addition, f constraints associated with relocated land resource, water quantity, water supply capacity, pollutant discharge and technique constraint can be formulated as follows.

1.
Reallocated land resource with market approach: In a traditional land utilization plan, the actual land resources can be reallocated to each plant by the proportion based on expected target.The market approach can be introduced to prompt land resources from lower value to higher value by the law of value, which can support land reallocation optimally.Model (2a) shows land resources reallocation based on market approach, where the productivities of land resources can be improved by reallocated actions (WRA tm and WRE tn ) based on total land resources (TA t ) (ha).TA t is the total land resources in study region (ha).

2.
Water quantity and water supply capacity for irrigative activities: Model (2b) presents available water for irrigation and forest without market approach, where a land plan associated with water quantity based on regional water resource load can be expressed.If water cannot satisfy the expected land targets, water deficits occur, which are caused by uncertain water availabilities.Water availability equals available water from surface and underground (QF tp and QG tp ) minus evaporation (E tp ), watercourse loss (H tp ) and minimum ecological requirement ). QF tp is the available water from surface (m 3 ); QG tp is the available water from underground (m 3 ); E tp is the total evaporation in study region (m 3 ); H tp is the watercourse loss (m 3 ); REC t is water conservation ability of forest per ha (m 3 /ha); and RO t is rainfall runoff coefficient (%).In Model (6b), since available water can be deemed as stochastic and random variables impacted by spatio-temporal factors, fuzzy measure Cr can be advocated to express such fuzziness, where µ is the credibility level through the QSF method (as shown in the Appendix A).Model (2c) shows the water supply capacity for irrigative activities in period t under probability p tp (m 3 ).The model presents that maximum supply capacity (i.e., CSM tp ) can be restricted by water availability.A market approach can prompt the efficiency of land plan; limited water resources can restrict the development of crop planting and environmental protecting.Thus, Model (2d) presents available water for irrigation and forest protection through a market approach, where water deficits occur when water cannot be delivered to the reallocated land.SRA tm and SRE tn are the water shortage area (ha).

3.
Pollution purification capacity through market approach: Total nitrogen allowance: Total phosphorus allowance: Soil and water conservation capacity: Model (2e) presents that capacity of purification from forest system (through ecological effect) with market approach hinges on the coefficient of purification (i.e., NP t , NN t ) under probability p tp (m 3 ) in period t.DA N tm and DA P tm are the actual pollution purification capacities through ecological effect with a market approach (ton).NP t and NN t are the coefficient of purification with consideration of ecological effect, which can be obtained based on previous research works.Models (2f) and (2g) present that pollutant discharges from crop irrigation would impose restrictions on discharge allowance (TNP tp and TPP tp ).TNP tp and TPP tp are maximum allowable TN and TP discharges from irrigation in period t (ton), which have been expressed as credibility fuzzy manners.In fact, actual nitrogen and phosphorus discharges would be original discharge from irrigative activities minus the values that being purification through ecological effect [as shown in Models (2f) and (2g)].Since the capacity of water conservation and soil erosion of forest can relieve the adverse effect from irrigation, Models (2h) and (2i) present that the total capacities of water and soil conservation would be restricted by their maximum capacities (LWT tp and LST tp ) in study region.LMS tm and LMW tn are the coefficients of ecological effect for soil and water conservation from a forest system in period t (%).LWT tp and LST tp are maximum allowances for water and soil erosion in period t (ton).7. Non-negativity: Model (2j) is non-negativity restrictions.

Data Acquisition
Table 1 presents the economic data related to different expected land demand targets and water deficits, which can be expressed quadratically.The net benefit of land utilization and the loss of water deficit can be calculated according to regional statistical yearbook indirectly [30][31][32][33].

Adverse Effects from Crop Irrigation without Market Approach
Figure 3 presents the water deficits without market approach when η are 0.6 and 0.99, which displays that water deficits would vary under various expected land demand targets..In the study region, since expected land demand has been regulated by policymakers, when available water cannot satisfy the expected land demand, irrigative production and forest protection cannot be conducted in the expected area, which can result in the reduction of output and economic losses.Results demonstrate that, if the pre-regulated land demand target is greater than the actual land that can be satisfied by water availability, water deficit can occur.Since approximately 60% of the total annual precipitation concentrates in the wet season, when water inflows are middle or low, the deficit would be strengthened.For example, in Period 1, water deficits of vegetable would be 3.17, 2.71, 1.41 0.41 and 0 × 10 6 m 3 when water inflows are low, low-medium, medium, medium-high, and high levels, respectively, when η is 0.6.Meanwhile, the results display that deficits are influenced by various η-levels.The higher η-level (i.e., η = 0.99) would result in higher water deficits; by decreasing η-level, water deficits would drop (i.e., η = 0.6).For example, when water availability is low level in Period 3, water deficits of shelter forest would be 3.16 × 10 6 m 3 (η = 0.6) and 3.75 × 10 6 m 3 (η = 0.99).In comparison, the highest water deficit area would occur in grain plant due to its excessive planting scale.On the contrary, economic forest has a lowest water deficit, which indicates that expected land demand targets for economic forest would be rational in study region.Meanwhile, it implies that economic forest has a higher potential to be expanded due to its rational scale and higher ecological effect.In study region, the overextended farming and irrational fertilization would generate excessive TN and TP pollutant discharges, which would bring about great stresses on regional environment.Thus, policymakers have regulated the maximum pollutant allowance according to regional environmental carrying capacity. Figure 4 shows excessive pollutant discharges without market approach when µ is 0.6.Several sensitive parameters of pollution discharge without market approach can be obtained, where water flow level, crop planation scale and corresponding pollutant discharge rate would influence the excessive pollution discharges.(a) Various water levels can generate changed capacities of pollutant diffusion.In general, a higher water availability level can produce a higher diffusion capacity.However, after more available water is allocated to irrigative activities, more pollutant would be discharged due to over-fertilization.For instance, when water levels are low and high in Period 3, the excessive TN from grain plant would be 12.01 and 14.32 × 10 6 ton, respectively.(b) The results determine that the crop planation scale and corresponding pollutant discharge rate of varied crop planation can generate different amount of pollutant Meanwhile, various policy scenarios associated with changed allowable discharges can be considered.In the study region, based on calculation of environmental loads, policymakers have restricted maximum allowances (TNP tp and TPP tp ) for TN and TP discharges.The results indicate that looser allowance would bring about a lower excess discharge, and vice versa.Excessive pollutant discharges without market approach when µ is 0.6.(total nitrogen denoted as TN, total phosphorus denoted as TP).

Land Trade between Forest and Irrigation
Figure 5 displays the land transaction between forest and irrigation when µ is 0.6.The results present that, although irrigative land can produce economic benefit, it can also bring about environmental damage due to over-cultivation and over-fertilization, leading to environmental loss and penalty.Thus, a crop-forest system with market approach can be introduced to maximize the regional benefits with a sustainable and natural circulation mode, where the adverse impacts from irrigated production can be converted by forest reserve contribution.Based on the calculation from QSF, the following results are obtained: (a) In an ICFM issue, the crop irrigation deemed as a seller withdraws irrigative land due to its higher penalty of environmental damage.The results present that the greatest areas of withdrawing cultivation would occur in grain plant planting (achieving 1.14 × 10 6 ha at maximum) due to its lower net economic benefit and higher emission compared to the other irrigated productions.(b) Forest can be considered a buyer to recover forest land to improve the system benefits from ecological effects.Among the three types of forest recovery methods, the economic forest has a better economic benefit and ecological return than shelter forest and ecological park.Under this situation, the highest recovered forest would be economic forest, which would reach 0.89 × 10 6 ha at maximum.(c) The trend of land transition indicates that the current oversized planting scale in grain plant based on food security policy would generate a higher risk of environmental pollution, which easily leads to system failure.(d) The results present that the Excessive pollutant discharges without market approach when µ is 0.6.(total nitrogen denoted as TN, total phosphorus denoted as TP).

Land Trade between Forest and Irrigation
Figure 5 displays the land transaction between forest and irrigation when µ is 0.6.The results present that, although irrigative land can produce economic benefit, it can also bring about environmental damage due to over-cultivation and over-fertilization, leading to environmental loss and penalty.Thus, a crop-forest system with market approach can be introduced to maximize the regional benefits with a sustainable and natural circulation mode, where the adverse impacts from irrigated production can be converted by forest reserve contribution.Based on the calculation from QSF, the following results are obtained: (a) In an ICFM issue, the crop irrigation deemed as a seller withdraws irrigative land due to its higher penalty of environmental damage.The results present that the greatest areas of withdrawing cultivation would occur in grain plant planting (achieving 1.14 × 10 6 ha at maximum) due to its lower net economic benefit and higher emission compared to the other irrigated productions.(b) Forest can be considered a buyer to recover forest land to improve the system benefits from ecological effects.Among the three types of forest recovery methods, the economic forest has a better economic benefit and ecological return than shelter forest and ecological park.Under this situation, the highest recovered forest would be economic forest, which would reach 0.89 × 10 6 ha at maximum.(c) The trend of land transition indicates that the current oversized planting scale in grain plant based on food security policy would generate a higher risk of environmental pollution, which easily leads to system failure.(d) The results present that the water availability levels would influence the land transactions between irrigated production and forest recovery.For instance, in Case 1, the area of withdrawing cultivation in oil planting would be 0.46, 0.22 and 0 × 10 6 ha when water levels are low, medium and high in Period 2, respectively.It indicates that land resources and water resources, both deemed important impact factors, could be incorporated into an ICFM project to adjust current policy in study region.Land trades between forest and agriculture when µ is 0.6 (grain plant is denoted as CP, oil plants is denoted as OP, vegetable is denoted as VP, economic forest is denoted as EP, shelter forest is denoted as SP, and forest park is denoted as FP).

Ecological Effects and Corresponding Benefit from Market Approach
Figure 6 presents the ecological effects with market approach when µ levels are varied.The results indicate the forest system (including economic forest, shelter forest and forest park) has played important roles in soil intention, water conversation and pollution purification, leading to economic returns to human activities.Among them, shelter forest has a higher contribution to soil loss reduction, while the forest park has played an important role in water conversation.Moreover, economic plant has a higher contribution to pollution purification by self-purification capacity, which can generate the highest benefit.Figure 5. Land trades between forest and agriculture when µ is 0.6 (grain plant is denoted as CP, oil plants is denoted as OP, vegetable is denoted as VP, economic forest is denoted as EP, shelter forest is denoted as SP, and forest park is denoted as FP).

Ecological Effects and Corresponding Benefit from Market Approach
Figure 6 presents the ecological effects with market approach when µ levels are varied.The results indicate the forest system (including economic forest, shelter forest and forest park) has played important roles in soil intention, water conversation and pollution purification, leading to economic returns to human activities.Among them, shelter forest has a higher contribution to soil loss reduction, while the forest park has played an important role in water conversation.Moreover, economic plant has a higher contribution to pollution purification by self-purification capacity, which can generate the highest benefit.
Figure 7 presents the amounts of land trading and corresponding ecological benefits when µ is 0.6.In this study, the ecological benefits can be direct or indirect.Among them, withdrawing reclamation (i.e., seller) can bring about the direct benefit of ICFM, which can reduce penalty of water/soil loss and pollution discharge.For instance, withdrawing reclamation of grain plant (denoted as GP) would reduce the damage to environment, which would bring about 49.6% of the direct benefit at maximum.The indirect benefit would be from the improvement of ecological function (e.g., improvement of water quality and elevation of water/soil conservation) due to forest recovery.For example, in Period 3, forest park can bring about the highest indirect benefits from water conservation, achieving 46.7% of the total benefits.The synergy of direct and indirect environmental impacts can improve the regional ecological effects such as pollution purification, water and soil conservation.
is denoted as SP, and forest park is denoted as FP).

Ecological Effects and Corresponding Benefit from Market Approach
Figure 6 presents the ecological effects with market approach when µ levels are varied.The results indicate the forest system (including economic forest, shelter forest and forest park) has played important roles in soil intention, water conversation and pollution purification, leading to economic returns to human activities.Among them, shelter forest has a higher contribution to soil loss reduction, while the forest park has played an important role in water conversation.Moreover, economic plant has a higher contribution to pollution purification by self-purification capacity, which can generate the highest benefit.Figure 7 presents the amounts of land trading and corresponding ecological benefits when µ is 0.6.In this study, the ecological benefits can be direct or indirect.Among them, withdrawing reclamation (i.e., seller) can bring about the direct benefit of ICFM, which can reduce penalty of water/soil loss and pollution discharge.For instance, withdrawing reclamation of grain plant (denoted as GP) would reduce the damage to environment, which would bring about 49.6% of the direct benefit at maximum.The indirect benefit would be from the improvement of ecological function (e.g., improvement of water quality and elevation of water/soil conservation) due to forest recovery.For example, in Period 3, forest park can bring about the highest indirect benefits from water conservation, achieving 46.7% of the total benefits.The synergy of direct and indirect environmental impacts can improve the regional ecological effects such as pollution purification, water and soil conservation .

Figure 7.
The amounts of land trading and corresponding ecological benefits when µ is 0.6 (note: pollutant purification is denoted as PP, soil conservation is denoted as SI and water conservation is denoted as WC).

System Benefit with and without Market Approach
Figure 8 presents the system benefits between market and non-market approaches.The results show that the system benefits can vary with η levels.For example, the system benefits would decrease from $1.01 × 10 9 to $1.24 × 10 9 , with increasing η-necessity levels from 0.60 to 0.995.It means that a lower η level and increased uncertainty for the imprecise objective would correspond to an optimistic attitude on the expected system benefit, and vice versa.Meanwhile, the market Figure 7.The amounts of land trading and corresponding ecological benefits when µ is 0.6 (note: pollutant purification is denoted as PP, soil conservation is denoted as SI and water conservation is denoted as WC).

System Benefit with and without Market Approach
Figure 8 presents the system benefits between market and non-market approaches.The results show that the system benefits can vary with η levels.For example, the system benefits would decrease from $1.01 × 10 9 to $1.24 × 10 9 , with increasing η-necessity levels from 0.60 to 0.995.It means that a lower η level and increased uncertainty for the imprecise objective would correspond to an optimistic attitude on the expected system benefit, and vice versa.Meanwhile, the market approach is a more effective method for land planning, which can generate a higher benefit than that with non-market approach.For example, when η is 0.9, system benefits would be $1.08 × 10 9 with market approach, while system benefits are $0.94 × 10 9 without market approach.

Conclusions
This study developed a mixed quadratic stochastic-fuzzy programming method (QSF) for handling uncertainties in an integrated crop-forest system with market approach (ICFM) issue for coordinating irrigative activities and forest protection.QSF can not only deal with uncertainties expressed as probability and possibility distributions, but also handle nonlinearity with an effective manner.The proposed model was applied to a real case in Xixian county, China.Various results can be obtained as follows: (a) The results show that the current land policy has generated a higher water deficit due to increased irrigative exploitation.The highest water deficit for irrigation would achieve 23.43 × 10 6 m 3 (α = 0.6, ε = 0.99), leading to a loss.It means that the accelerated irrigative expansion has exceeded what water resources can afford.Thus, the adjusting the rate of agricultural development is an important issue for regional policymaker.(b) The results display that existing land utilization would produce excessive pollution to increase ecological risk in study region, where the highest excessive TP and TN discharges would reach 432.24 and 14.32  10 3 ton.(c) The results indicate that strategies for ecological expansions (convert cultivated land into forest) can decrease the excessive pollutant discharges and water deficits, which can generate a higher system benefit than that without withdrawing farmland and recovering forest.
Meanwhile, through the calculation from a practical ICFM issue in Xixian county (China), the obtained results determine that market approach can prompt land change between crop and forest systems, with aim to avoid higher retreatment costs from over-expanded cultivation and irrational fertilization, while achieving higher benefits from ecological returns.Through application of Xixian county, we found that the largest withdrawing reclamation is grain plantation, which would bring about 49.6% of the direct ecological benefit at maximum.Meanwhile, various forest types would generate varied ecological function, which can bring about indirect benefit.Among them, forest park can bring about the highest indirect benefits from water conservation, achieving 46.7% of the total benefits.Shelter forest plays an important role in soil conservation.In an integrated cropforest system with market approach (ICFM), market approach can grade land resources from lower to higher value, which can improve the productivities of land resources.Compared to a crop-forest system without market approach, it is found that the system benefit under market approach would be higher than that with non-market approach, which indicated that market would be a more effective method for land planning in study region.
Moreover, several discoveries have been made: (a) The tradeoff between the benefits of

Conclusions
This study developed a mixed quadratic stochastic-fuzzy programming method (QSF) for handling uncertainties in an integrated crop-forest system with market approach (ICFM) issue for coordinating irrigative activities and forest protection.QSF can not only deal with uncertainties expressed as probability and possibility distributions, but also handle nonlinearity with an effective manner.The proposed model was applied to a real case in Xixian county, China.Various results can be obtained as follows: (a) The results show that the current land policy has generated a higher water deficit due to increased irrigative exploitation.The highest water deficit for irrigation would achieve 23.43 × 10 6 m 3 (α = 0.6, ε = 0.99), leading to a loss.It means that the accelerated irrigative expansion has exceeded what water resources can afford.Thus, the adjusting the rate of agricultural development is an important issue for regional policymaker.(b) The results display that existing land utilization would produce excessive pollution to increase ecological risk in study region, where the highest excessive TP and TN discharges would reach 432.24 and 14.32 × 10 3 ton.(c) The results indicate that strategies for ecological expansions (convert cultivated land into forest) can decrease the excessive pollutant discharges and water deficits, which can generate a higher system benefit than that without withdrawing farmland and recovering forest.
Meanwhile, through the calculation from a practical ICFM issue in Xixian county (China), the obtained results determine that market approach can prompt land change between crop and forest systems, with aim to avoid higher retreatment costs from over-expanded cultivation and irrational fertilization, while achieving higher benefits from ecological returns.Through application of Xixian county, we found that the largest withdrawing reclamation is grain plantation, which would bring about 49.6% of the direct ecological benefit at maximum.Meanwhile, various forest types would generate varied ecological function, which can bring about indirect benefit.Among them, forest park can bring about the highest indirect benefits from water conservation, achieving 46.7% of the total benefits.Shelter forest plays an important role in soil conservation.In an integrated crop-forest system with market approach (ICFM), market approach can grade land resources from lower to higher value, which can improve the productivities of land resources.Compared to a crop-forest system without market approach, it is found that the system benefit under market approach would be higher than that with non-market approach, which indicated that market would be a more effective method for land planning in study region.
Moreover, several discoveries have been made: (a) The tradeoff between the benefits of ecological function and economic objectives can support producing several results, with aim to achieve harmony between economic development and ecological protection.(b) Although market approach is an effective method to increase the economic productivity of land resources, the occurrences of market failure would impede the process of withdrawing reclamation and recovering forest.(c) Risk preferences of policymakers can influence the outcome against risk-violation in an ICFM issue.Therefore, many corresponding suggestions can be summarized as follows: (a) a proper strategy associated with irrigative development and forest protection can not only protect regional land functions, but also facilitate minimizing the adverse impacts from human activities; (b) the improvement of marketable productivity and regulation of market order can promote the efficiency of land plan with market approach; and (c) the consideration of risk can prompt the reliability of land utilization plans, which should be accounted for in a policy-making process.

Figure 2 .Figure 1 .
Figure 2. Framework of land utilization pattern with consideration of ecological effect under uncertainties.

Figure 2 .
Figure 2. Framework of land utilization pattern with consideration of ecological effect under uncertainties.

Figure 2 .
Figure 2.Framework of land utilization pattern with consideration of ecological effect under uncertainties.
Sustainability 2018, 10, x FOR PEER REVIEW 11 of 19 policymakers have restricted maximum allowances ( tp TNP and tp TPP ) for TN and TP discharges.The results indicate that looser allowance would bring about a lower excess discharge, and vice versa.

Figure 3 .
Figure 3. Water deficits without market approach when µ are 0.6 and 0.99.(low water flow denoted as L, low-medium water flow denoted as LM, medium water flow denoted as M, medium water flow denoted as MH).

Figure 3 .
Figure 3. Water deficits without market approach when µ are 0.6 and 0.99.(low water flow denoted as L, low-medium water flow denoted as LM, medium water flow denoted as M, medium water flow denoted as MH).

Figure 3 .
Figure 3. Water deficits without market approach when µ are 0.6 and 0.99.(low water flow denoted as L, low-medium water flow denoted as LM, medium water flow denoted as M, medium water flow denoted as MH).

Figure 4 .
Figure 4. Excessive pollutant discharges without market approach when µ is 0.6.(total nitrogen denoted as TN, total phosphorus denoted as TP).

Figure 4 .
Figure 4. Excessive pollutant discharges without market approach when µ is 0.6.(total nitrogen denoted as TN, total phosphorus denoted as TP).

Sustainability 2018 , 19 Figure 5 .
Figure 5. Land trades between forest and agriculture when µ is 0.6 (grain plant is denoted as CP, oil plants is denoted as OP, vegetable is denoted as VP, economic forest is denoted as EP, shelter forest is denoted as SP, and forest park is denoted as FP).

Figure 6 .
Figure6.The ecological effects with market approach when µ is 0.6.

Figure 6 .
Figure 6.The ecological effects with market approach when µ is 0.6.

Sustainability 2018 , 19 Figure 8 .
Figure 8. System benefits between market and non-market approach.

Figure 8 .
Figure 8. System benefits between market and non-market approach.
Note: W M tm and WEA tn are expected land demands for crop and forest in period t (ha); SM tm and SEA tn are water deficit areas due to uncertain water availability in period t (ha).

Table 2
presents modeling input such as allowance TN/TP permit, water and soil erosion rate and TN/TP discharge rate in study region according to previous research works.
Note: TP is total phosphorus, and TN is total nitrogen.