Abstract
To promote urbanization in the next stage, it is of great significance to explore the ecological efficiency of green ecological regions and advance the sustainable development of a social economic system. However, spatial heterogeneity has not been fully considered in the existing evaluation models or methods for regional ecological efficiency (REE), and the corresponding decisions on sustainable development are not the optimal solutions. To solve the problems, this paper explores the evaluation of REE and intelligent decision support for sustainable development by analyzing environmental big data. Firstly, the spatiotemporal evolution of REE was examined based on environmental big data to clarify the spatial layout of REE and the sources of the spatial differences. Next, a multiobjective optimal decision-making model was established for the sustainable development of a regional ecosystem, and the solving method was presented for the model. The proposed model was proved valid through experiments.
1. Introduction
Urban development faces several problems, namely, the development model is extensive, the newly developed areas are underpopulated, and the resources are in short supply. With the development of the industry, these problems have intensified the contradiction between population, environment, and resources [1–5]. Rather than pursuing the single goal of economic interest, the sustainable development of the urbanization system aims to organize production for the composite goal of resources/energy, environment, economy, and society.
China is vigorously implementing a strategy of regional ecosystem development and environmental protection, which will continuously enhance the coordination between the regional environment, regional economy, and regional society [6–12]. Against this backdrop, in order to promote urbanization in the next stage, it is of great significance to explore the ecological efficiency of green ecological regions and advance the sustainable development of the social economic system [13–21].
Environmental pollution is a potential bottleneck of stable economic growth. As an effective measure of sustainable development, ecological efficiency can fully reflect the actual level of coordinated development between economy and environment. Liu and Sun [22] constructed a data envelopment analysis (DEA) model with environmental pollution and resource consumption as inputs and total economy as the output. Ratner [23] described the monitoring of ecological and economic efficiencies of the activities in the regional economic system, compared the applicability of several methods to the dynamic DEA model, and attempted to generate time series based on environmental and economic efficiency points. To evaluate the performance of the urban system, Giordano et al. [24] divided the urban system into traffic subsystem, built environment subsystem, and social economy subsystem, evaluated the efficiency of each subsystem by Takagi–Sugeno model, and recursively derived the overall efficiency of the entire system. Hoang and Alauddin [25] designed an input-oriented DEA framework, which allows the measurement and decomposition of economic, environmental, and ecological efficiencies in agricultural production of different countries.
On the sustainable development of the urban ecological system, Wang and Li [26] set up an overall framework for the sustainable urban spatial development model based on smart cities. The framework, involving such four dimensions as the monitoring, collection, interconnection, and sharing of urban data, provides an important guide for the spatial planning and construction of smart cities. Ma et al. [27] analyzed the evolution of an urban ecosystem and its dissipation structure and constructed an evaluation index system (EIS) for the sustainable development capability of the urban ecosystem.
The existing studies provide some useful references. However, the evaluation models and methods for regional ecological efficiency fail to fully consider spatial heterogeneity, take account of the dynamic evolution of ecological efficiency distribution, and thoroughly discuss the key factors affecting the ecological efficiency identification and animal habitat changes in different regions and at different levels. Concerning sustainable development, the available decision-making models overlook the balance of relative development for the regional ecosystem in spatiotemporal evolution and cannot converge to the optimal decisions. To solve the problems, this paper explores the evaluation of REE and intelligent decision support for sustainable development based on environmental big data. Section 2 analyzes the spatiotemporal evolution of REE based on environmental big data, identifies the spatial layout and spatial difference of REEs, and discusses the key environmental factors. Section 3 establishes a multiobjective optimal decision-making model for the sustainable development of the regional ecosystem and presents the solving method for the model. Experimental results demonstrate the effectiveness of our model.
This research discusses the inherent linking mechanism of system sustainable development with environmental, economic, and social impacts, and optimizes the absolute performance of sustainability and relative development simultaneously.
2. Spatiotemporal Evolution Analysis
This paper identifies the sources of environmental data for the spatiotemporal evolution of REE. Table 1 lists the REE evaluation items based on big data. Several environmental types are enumerated, namely, building area, traffic area, water area, and landscape area, and the corresponding factors that affect ecological efficiency are presented in detail.
As mentioned before, this paper aims to clarify the spatial layout of REE and the sources of the spatial differences and realize the REE evaluation and make intelligent decisions for sustainable development. For this purpose, the spatial distribution of ecological efficiency was described intuitively based on environmental big data. Then, the degree of spatial differences of REE was investigated, and the sources of these differences were revealed. Finally, the dynamic evolution law of REE was explored through kernel density estimation (KDE) and Markov chain.
2.1. Spatial Difference Analysis
Through Dagum’s decomposition of the Gini coefficient, this paper analyzes the spatial differences of global and local ecological efficiencies of the ecosystem in the study area. Dagum’s decomposition method was selected for its capability of considering the overlap between subsamples, solving heterogeneous sources of REE spatial differences, and measuring how much overall regional differences are affected by the intraregional difference, inter-regional difference, and transvariation intensity.
Let buv and bgs be the urban ecological efficiencies of the u-th and -th regions, respectively; let λ be the mean ecological efficiency of all cities in a region, m be the total number of samples (cities), l be the number of regions, and mu and be the number of cities in the u-th and -th regions, respectively. Then, the relative difference between all cities in a region in ecological efficiency can be measured by the regional overall Gini coefficient:
To decompose the Gini coefficient by a regional subsystem, the first step is to sort all cities in a region by mean ecological efficiency is as follows:
Ecological efficiency is the ratio of ecological outputs to ecological inputs. The outputs refer to the value of products and services provided by enterprises or economies, while the inputs refer to the resources and energies consumed by enterprises or economies, as well as the environmental load caused by enterprises or economies. The Gini coefficient of the u-th region can be calculated by
The intraregional contribution difference Hq can be calculated by
The inter-regional Gini coefficient GNug between the u-th and -th regions can be calculated by
The inter-regional contribution difference Hmr can be calculated by
The transvariation intensity contribution Hd can be calculated bywhere ei = miλi/(mλ) and tj = mj/m. Let cug be the ecological efficiency difference between regions; tug be the first moment of transvariation. Then, the relative influence Cug of the unit ecological efficiency of the u-th and -th regions can be calculated by
Let Ou and be the cumulative density distribution functions of the u-th and -th regions, respectively. Then, cug can be characterized by the mathematical expectation of the sum of the positive sample values for the difference between the ecological efficiencies, buv and bgs, of the cities in the u-th and -th regions:
In addition, tug can be calculated by
2.2. Dynamic Evolution of Distribution
The traditional KDE, a popular tool in spatial nonequilibrium analysis, is a nonparametric estimation approach. It can describe the form of each random variable with a continuous curve. In this way, virtually no statistical error will be incurred by the improper setting of ecological efficiencies in the region. Let M be the total number of cities in the region, Ai be the independent identically distributed observations, be the mean value, τ be the bandwidth, and Γ(·) be the kernel function. Then, the density function μ(a) of the random variable A can be calculated by
The greater the τ value, the larger the neighbourhood of a. If τ is too large, the kernel density function Γ(·) will be too smooth. Then, some important features of the function will be smoothened out, causing a large deviation. Hence, the τ value should be minimized if possible. Thus, Γ(·) needs to meet the following condition:
The dynamic evolution features of REE distribution can be estimated by the Gaussian kernel function:
2.3. Long-Term Transition Trend
This paper adopts the Markov chain to describe the ecological efficiency of each city in a region, that is, sets up the transition probability matrix of the Markov chain. Let {A(p), p∈ψ} be the stochastic process of the discrete events corresponding to the Markov chain, N be the finite set of the solutions, and ψ be the set of indices for different stages of the stochastic process. The Markov chain highlights that, historical dynamic actions share the same features as future dynamic actions. In our problem, the identicalness of the features manifests as follows: the state of urban ecological efficiency A in year p + 1 directly bears on the probability for A to belong to type j in year p. Let GRij be the probability for urban ecological efficiency belonging to type i in year p to transfer to type j in year p + 1; that is, the maximum likelihood estimation of the transition is GRij = mij/mi; let mij be the number of cities shifting from type k in year p to type j in year p + 1 and mi be the total number of cities belonging to type i in the study period. Then, we have
By dividing urban ecological efficiencies in the region into x states, it is possible to obtain an x × x state transition probability matrix of the Markov chain.
3. Multiobjective Optimal Decision-Making Model
3.1. Preliminaries
Figure 1 summarizes the evaluation attributes for the sustainable development of the regional ecosystem. The evaluation system for the sustainable development of the regional ecosystem involves multiple dimensions, such as environmental sustainability, economic sustainability, social policy sustainability, and technical sustainability. The evaluation dimensions can all be quantified. For the sustainable development of the regional ecosystem, decision makers can also set up individualized evaluation systems based on objective needs and subjective judgement.

To clarify the incremental costs and benefits of regional ecosystem development, this paper focuses on the incremental costs incurred during the optimization of plan design and measures implementation for regional ecosystem development. The costs of these two stages are mainly invested continuously by green project investors, developers, and implementers. The incremental benefits will continue to appear after the optimization measures are implemented in the regional ecosystem. Figure 2 details the costs and benefits of the regional ecosystem in a development cycle.

Our problem is to sustainably improve the regional ecosystem, that is, combining and optimizing multiple optional measures for the ecosystem optimization to realize multiple goals under specific constraints, in reference to multiple sets of environmental big data in different periods, such that the optimized regional ecosystem is much more sustainable than the original system. This paper builds a weighted multiobjective optimization model to improve the sustainability of the regional ecosystem in multiple dimensions, relying on optimization measures. On this basis, the authors strived to make sustainable intelligent decisions.
3.2. Model Construction and Solving
Targeting the regional ecosystem, our model tries to improve the ecological efficiency of the ecosystem by adopting the improved combination of optimization measures for sustainable development. Figure 3 provides the specific steps to improve sustainable development efficiency: characterization of the regional ecosystem, sustainable development evaluation of REE, sustainable development optimization of REE, and sustainable development improvement of REE.

The optimal decision making for the sustainable development of REE involves multiple steps and goals in four dimensions, including environmental sustainability, economic sustainability, social policy sustainability, and technical sustainability. The decision-making process is complicated by the diverse constraints on the different goals and the interference of uncertain decision-making factors. Hence, this paper presents a multiobjective optimal decision-making model (Figure 4) for the sustainable development of the regional ecosystem. The proposed model consists of two parts, namely, the characterization of the ecosystem, and construction and solving of sustainable development decision-making model and opens the path to improve the sustainable development of the regional ecosystem through spatiotemporal evolution.

The proposed multiobjective optimal decision-making model is targeted at regional ecosystems, whose ecological efficiencies obey different spatial distributions. The model was constructed through composite weighting of each subobjective, setting weighted objectives and sustainable objectives, and solving the model. Firstly, the interval analytic hierarchy process (IAHP) was adopted to assign composite weight to each objective, eliminating the need to consider the interaction between subobjectives. The composite weighting begins with the construction of an interval comparison matrix.
The subjective weight of each subobjective was quantified through expert judgement, which is uncertain to a certain extent. Therefore, the comparison matrix is composed of the following interval numbers:where m is the number of optimization subobjectives; is the relative importance of subobjective i to subobjective j, both of which are expressed as interval numbers. Next, the interval comparison matrix F was converted into the interval priority matrix S:where
After that, the subjective weights were calculated. Let , , , and be the positive and negative biases in objective programming, respectively; be the interval value range of the relative importance for each subobjective corresponding to the optimal solution. Then, the optimal solution of matrix S can be obtained by
The subjective weight of each objective can be obtained by
The overall objective for the sustainable development of the regional ecosystem combines multiple subobjectives, which are impossible to achieve all at once. To realize the overall objective, the traditional subobjective setting approach needs to adjust the subobjectives repeatedly. Based on multiple sets of environmental big data on the optimization system in different periods, this paper puts forward a novel weighted subobjective setting strategy.
According to the weights of subobjectives for the sustainable development of the regional ecosystem as well as the sustainable performance of each subobjective reflected by the multiple sets of environmental big data, the weighted subobjectives can be configured under a constraint as follows:where be the performance value of subobjective i reflected by the normalized environmental sample data j (If ψj = 1 or 0 in ψj(1 − ψj) = 0, then the j-th set of environmental sample data can or cannot serve as a reference); D is the constraint; T ∈ D is the requirement that any combination of optimization measures must satisfy D.
Through the above steps, all weighted subobjectives can be obtained as . Each weighted subobjective corresponds to a specific model. Each combination of optimization measures leads to the optimal solution to each specific model.
This paper likens the set of optimal solutions on all weighted subobjectives as the ideal solution to the objective of sustainable development, which can be vectorized as
The ideal solution above can be transformed to an actual solution under the condition . For any feasible solution, the corresponding combination of optimization measures Tb can be vectorized as
Considering the length and direction of vector functions, the degree of absolute improvement of sustainable development in the ideal solution IS can be given by
The degree of absolute improvement of sustainable development in the feasible solution b can be given by
The equilibrium of relative development of multiple weighted subobjectives in the ideal solution IS can be given by
The equilibrium of relative development of multiple weighted subobjectives in the feasible solution b can be given by
By the principle of vector projection, the equilibrium of relative development of weighted subobjectives can be unified with the absolute improvement of sustainable development. In ideal conditions, the optimal degree of improvement of regional sustainable development can be given by
The improvement effect of the feasible solution Tb on sustainable development can be quantified by
The intelligent decision making of the objective of sustainable development can be derived from formulas (27) and (28). To display the decision-making process more intuitively, an objective achievement rate was introduced to transform the two formulas:
The above formula shows that the objective achievement rate of IS is 1 and that of Tb falls in between (0, 1) (the closer the value is to 1, the better the optimization effect). Through the above analysis, the overall optimization objective of sustainable development for the regional ecosystem can be expressed as
The above formula also outputs an interval number. Under uncertainty conditions, the model can be solved by probability comparison:
Any interval numbers and can be compared in terms of probability. If Wa > Wt, i.e., S(a>t), S(a>t) > 0.5 means Wa is probably better than Wt; the inverse is also true. Similarly, if and , Wa must be better than Wt.
4. Experiments and Result Analysis
Figure 5 presents the evolution trends of the overall ecological efficiency of the study area in 2005–2020. The REEs exhibited obvious features of phased development, dropping from 0.6024 in 2006 to the lowest point of 0.5785 in 2013. From 2014 to 2020, the REEs rose significantly and peaked at 0.6549. The reasons for the phased development are as follows:

In the early phase of the sample period, the regional economy developed extensively in the traditional model, and the society had poor awareness of ecoenvironmental protection. That is why the REEs in the study area continued to decline. In 2013–2014, the cities in the study area quickly adjusted their development pattern and publicized ecological civilization, which contribute to the continuous growth of REEs.
There were some differences between regions in ecological efficiency. Region A, a transportation hub in a strategic location, remained the leader of ecological efficiency throughout the 15-year-long sample period. The ecological efficiencies of Regions A and B oscillated similarly as the overall ecological trend of the study area. After 2015, the ecological efficiency of Region B increased markedly each year, with increasingly strong momentum. Region C had the lowest ecological efficiency, which grew stably over the 15 years.
Figure 6 compares the REE evolution trends of each pair of the three regions. Through the sample period, the largest regional difference in ecological efficiency existed between Regions A and B, whose mean Gini coefficient was 0.274. The second-largest regional difference was observed between Regions A and C, whose mean Gini coefficient was 0.261. The smallest regional difference was observed between Regions B and C, whose mean Gini coefficient was 0.253.

In general, the ecological efficiency gap between Regions A and C continued to narrow at an annual change rate of −2.14%. The gap between Regions B and C widened after 2010, despite some fluctuations, with an annual change rate of 0.87%. Specifically, the gap between Regions A and C declined with fluctuations from 0.347 in 2006 to 0.206 in 2020, reaching the valley of the sample period. The gaps between Regions A and B, and between Regions B and C both decreased first and then increased. The gap between Regions A and B minimized at 0.223 in 2016, while that between Regions B and C minimized at 0.204 in 2010.
To sum up, the relevant government departments should strengthen the coordination between Regions B and C in green development and reduce the spatial difference between regions in ecological efficiency without ignoring the regional difference between A and B.
Through Gini coefficient decomposition, this paper looks for the sources of REE spatial differences, that is, analyzes the contributions of intraregion difference, inter-regional difference, and transvariation intensity to overall regional spatial differences. Table 2 presents the calculated results on intra and inter-regional contributions.
Through the sample period, the inter-regional difference of A and B contributed the greatest to the spatial gap of REE, followed in turn by that of B and C, and that of A and C; the intraregional difference of C contributed the greatest to the spatial gap of REE, followed in turn by that of A and that of B. It can be observed that the spatial imbalance of REEs can be effectively solved by reducing inter-regional differences. Regions A and B should be the focus of regulation because their inter-regional difference is the major contributor to the spatial gap of REE.
Tables 3–5 report the transition probability matrices of the Markov chain for Regions A–C, respectively. Note that period p1 has a one-year lag, and period p2 has a two-year lag. The long-term transition trend of REE can be derived from these matrices. In each matrix, the elements to the left of the diagonal were smaller than those to the right, indicating that the overall ecological efficiency of each region moves clearly from the low level to the high level. Meanwhile, the REE had a greater stable probability than transition probability. This means that the state transition of ecological efficiency is accompanied by apparent polarization. As period p1 changed to period p2, i.e., the lag increased from 1 year to 2 years, the overall ecological efficiency of the study area basically remained the same, with a very slight decline. In addition, there were many nonzero elements on both sides of the diagonal in each matrix, suggesting the difficulty for the study area to realize the leapfrog transition of ecological efficiency.
This paper employs the multistage method to measure the contributions of intraregion difference, inter-regional difference, and transvariation intensity in 2005–2020. The results are presented in Table 6 and Figure 7. The measured results of our method have a smaller variation than those of the traditional DEA. The following can be inferred from the data in Table 6. In 2005–2020, the mean contribution was 0.813, 0.990, and 0.821 for intraregional difference, inter-regional difference, and transvariation intensity, respectively. Overall, the study area had relatively high ecological efficiency since 2005. The regional economic development achieves a high cost effectiveness in terms of energy-saving technology and green environmental investment. In 2010 and 2013, the contribution of inter-regional difference was merely 0.758 and 0.761, respectively. In 2009–2011 and 2013–2016, the inter-region contributions were below the mean of the 15 years. Through the 15-year-long sample period, the transvariation intensity made relatively high contributions, reflecting the good effect of technical and management measures adopted by the regional ecosystem.

5. Conclusions
Based on the analysis of environmental big data, this paper investigates the evaluation of REE and intelligent decision support for sustainable development. Specifically, the spatiotemporal evolution of REE was analyzed based on environmental big data, and the spatial layout of REE and the sources of the spatial differences were both clarified. In addition, the authors constructed a multiobjective optimal decision-making model for the sustainable development of a regional ecosystem and presented the solving method for the model. After that, experiments were carried out to analyze the REE changes and inter-regional REE differences in 2005–2020. The experimental results demonstrated the feasibility of our analytical method. Furthermore, the authors discussed the contributions of intraregional difference, inter-regional difference, and transvariation intensity to the overall spatial difference of regional ecological efficiency, constructed the transition probability matrix of the Markov chain for the REE in each region, and provided the measured results on regional REEs.
The future research will try to develop a decision model capable of considering and solving different interest preferences and demands of different decision makers and realize multirole, multiattribute decision making as well as multirole, multiobjective decision making.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Humanities and Social Science Foundation of Education Department of Hubei Province (Grant no. 194101), the Young Scientists Fund of Humanities and Social Sciences of the Ministry of Education of China (Grant no. 19YJC630236), and the Philosophy and Social Science Research Project of the Hubei Provincial Department of Education (Grant no. 19Q093).