Multi-Scale Analysis of the Evolution of Jiangsu’s Ecological Footprint Depth and Its Factor Decomposition

: The ecological footprint (EF), as a set of land-based ecological indicators, plays an important role in land ecology and evaluations of ecological pressure. Multi-scale levels of Jiangsu’s three-dimensional EF were analyzed, and 3D maps were presented to demonstrate the geographical distribution of the ecological footprint depth (EFD) of Jiangsu’s counties in 1995–2015 at the geographic scales of prefecture-level cities and counties. The results show that the overall EFD of Jiangsu gradually increased during the study period. The county-scale results show that the distribution of EFDs was high in the south and low in the north, and EFDs were mainly concentrated in urban areas of prefecture-level cities. The logarithmic mean Divisia index (LMDI) was used to decompose the factors in explaining the change in EFD. The LMDI analysis shows that the changes in factors every year differ among geographical units on different scales. Afﬂuence is the main factor that promotes EFD, and the change in the ratio between EFD and scientiﬁc and technological level is the main factor that suppresses EFD. Countermeasures and suggestions for balancing ecological pressure in speciﬁc regions and reducing the depth of the EF from various factors with multi-scale heterogeneity are suggested.


Introduction
The issue of sustainable ecological development is a major concern for the ecological security of human beings [1], and the relationships between ecology and society are challenging [2]. Humans have exerted more pressure on the biosphere and climate change than its capacity [3]. In particular, the problems of land ecological pressure [4] and carbon emissions [5] in relation to sustainability have received worldwide attention and research. How to accurately and objectively measure human pressure on the ecological environment is an important part of sustainable development research. The ecological footprint (EF) model has played an increasingly important role in accounting for global ecological boundaries and evaluating ecological pressure [6], land use [7], and carbon footprint [8,9]. The EF developed by [10,11], as a concept and model, has several main advantages: the ecological footprint model takes the global standard unit of land as the evaluation index, which can uniformly link the land use with the carbon footprint; it provides a good perspective for the research of other types of pollution; and it is extended to ecology-related research contents such as the water footprint [12], land footprint [7], and material footprint [13].
The EF concept is being improved to be applicable in different scenarios, e.g., the input-output method (IO) is used to study the flow of embodied footprints [14], LCA is used to study more specific individual ecological footprints, and the emergy-based EF is applied in situations where it is difficult to calculate the original ecological footprint. The three-dimensional ecological footprint (3D EF) was introduced [15,16] to explore the relationships between ecological flows and natural capital. The 3D EF is a measure of how much of the current biocapacity supplies the ecological consumption demand [17]. In comparison with other EF accounting improvements, 3D EF provides a new perspective Land 2022, 11, 1997 2 of 20 for EF accounting, especially in the relationship between ecological flows and natural capital [15,18]. The 3D EF model was measured by EF size (EFS) and EF depth (EFD). EFS is a concept related to biocapacity (BC), and EFD is a concept related to ecological pressure. Compared to ecological pressure with traditional EF, the measurement of EFD is a more dynamic process because the measurement of the flow and stock of natural capital in 3D EF provides a dynamic description of the relationship between ecological occupation and ecological carrying capacity [17]. Therefore, the 3D EF model is more efficient for measuring ecological security and ecological pressure [19], and many studies of principles, accounting, and explanations of the model have been conducted, e.g., the types of interactions between flows and stocks of energy and materials from the 3D EF perspective were discussed [17]. Table 1 summarizes the study areas, focused indicators, methods and key findings of recent literature on 3D EF. Most of them raise the advantage of 3D EF in multi angles, e.g., Tang et al. used 3D EF to explore the ecological carrying capacity of the Greater Bay Area of China and found the modified 3D EF is more significant than the traditional EF model in driving factor explanations [20]. The regional scales of most of the studies were national scales [18], provincial scales [21], and prefecture-level scales [22,23], while there are fewer studies on the regional 3D EF at smaller scales or at multiple scales. With multiple scales of 3D EF, research can better explore ecological spatial changes and spatial relationships. Due to the close relationship between the economic society and EF, exploring the relationship between socioeconomic factors and changes in EF is an important part of proposing countermeasures and suggestions for optimizing EF. The influencing factors of EF have been studied from multiple angles [1,[26][27][28][29][30], and the influencing factors of 3D EF have also been studied, especially regarding the footprint depth, due to the rarity of changes in 3D EF. Most of the studies implemented or extended the IPAT or STIRPAT model [31], and the explanatory model forms included the exploration of linear relationships, combined with the EKC model to explore the quadratic relationship of affluence [32] and with a spatial econometric model [33]. Some studies focus on research in a certain industry, e.g., Lv et al. [34] used STIRPAT and a spatial econometric model to study the spatial relationship between transportation networks' 3D EF and socioeconomic factors. P. Li, Zhang, and Xu [24] used partial regression analysis (PLS) to explore the influencing factors of Urumqi's 3D EF in China, and their study focused on the ecosystem service values in 3D EF. They concluded the main driving factors for EFD were built-up area, population, and per capita GDP. Compared with those regression models, factor decomposition can analyze the influence of factors in each time period in more detail. For factor decomposition models, the logarithmic mean Divisia index (LMDI) [35,36] provides full decomposition with no residuals. This index is widely used in the decomposition of carbon emissions [37], energy consumption [38], and PM2.5 [39]. However, LMDI it is rarely used in the decomposition EF or EF depth; nevertheless, it is applicable to EF, as the accounting mechanism and expression of ecological occupancy content are similar to those of carbon emissions.
For the driving factors (variables) of EF depth or EF from Table 1, most of the studies on EFD focused on the driving factors of GDP (per capita GDP) [20,21,23,24], population [24,25], income [18], industrial structure [22,25], etc. What is more, the process of urbanization is generally considered; population agglomeration, technological progress, and other aspects of human activities in the process of urbanization have both a demand and an impact on EF. Population and economic growth are generally the main positive influencing factors of EF and EF depth in many studies [21]. There are different views on the impact of science and technology improvement on ecology. Ecologist Commoner B. believes that the most important factor leading to the deterioration of the ecological environment is the use of industrial technology. In contrast, modern ecological theory contends that economic and social development and ecological environmental protection can be simultaneously achieved through scientific and technological progress and development to promote economic and environmentally friendly growth and upgrade the industrial structure. Industrial structure has a significant effect on EF depth [22].
From the above literature review, the following points could be concluded. (1) Most studies were carried out on a province or prefecture-level city scale; the EFD at a small scale or comparison with multi-scales are less common. As the ecological pressure has a great relationship with the urbanization population concentration and accounting of ecological pressure varies significantly at different scales, such as at the country scale and the city scale, it is of more ecological significance for EFD studies at multi-scale. (2) For the EF factor analysis, more studies paid more attention to EF values even though they employed 3D EF, which was close to the analysis framework of the traditional EF, and fewer studies were mainly concentrated on EFD, which was essential for 3D EF. In particular, there were few studies on ecological sustainability analysis from the concentrated perspective of EFD and ecological pressure. (3) Most of them used OLS or PLS regression or correlation analysis for EF driving factor explorations, while there was little research on the staged comparison of the factor decomposition. This kind of research can better describe the difference in stage change, while the former ignores comparing stage change differences.
Jiangsu Province of China was chosen as the study area, as it is representative of China's high ecological pressure and economic development; as measured, Jiangsu's EF exceeded its BC by 6.28 times in 2015.
Based on the above analysis for the EFD analysis improvement, this study tries to fill the gap and take Jiangsu Province as a demonstration for analyzing the EFD and decomposing its social and economic factors in multi-scales. Specifically, this paper analyzes the distribution, features, and evolution of the EFD of Jiangsu Province from provincial, prefecture-level city, and county-level scales, taking 5 years as a time interval in the period 1995-2015, corresponding to China's economic planning from the 9th through the 12th five-year plan. The main content is as follows: (1) Accounting and mapping of the EFDs of Jiangsu are performed. As mapping is important in ecosystem studies [40], this study used ggplot2 graphic tools to vividly illustrate the distribution and geographic variation in EFDs among Jiangsu's counties and prefecture-level cities. The changes in each factor and the changes in EFD are analyzed and the mechanistic relationships between them are explored. Suggestions for ecological balance and sustainable development from the 3D EF perspective are provided.

Study Area
Jiangsu Province is located in East China. As an economically developed city in China, Jiangsu has a large population density and low per capita BC. As a province in the forefront of China in terms of urbanization level, per capita GDP, and technological level, Jiangsu is considered to have a relatively high ecological footprint pressure. This study has a certain reference for research on other provinces in China. As shown in Figure 1, there are 13 prefecture-level cities in Jiangsu. Each prefecture-level city is composed of county-level cities, counties, new urban districts, and traditional urban areas. This paper divides Jiangsu Province into 73 units of districts and counties; for the convenience of description, these units are collectively referred to as counties.
(2) The LMDI is used to decompose the changes in EFDs with economic and so tors. The influencing factors of EFD changes are compared among different and in different time stages. This paper is based on the model of population fluence (A), technology level (T), and industrial structure (S) to construct th decomposition list. (3) The changes in each factor and the changes in EFD are analyzed and the mec relationships between them are explored. Suggestions for ecological balance a tainable development from the 3D EF perspective are provided.

Study Area
Jiangsu Province is located in East China. As an economically developed China, Jiangsu has a large population density and low per capita BC. As a provinc forefront of China in terms of urbanization level, per capita GDP, and technologic Jiangsu is considered to have a relatively high ecological footprint pressure. Thi has a certain reference for research on other provinces in China. As shown in F there are 13 prefecture-level cities in Jiangsu. Each prefecture-level city is comp county-level cities, counties, new urban districts, and traditional urban areas. Thi divides Jiangsu Province into 73 units of districts and counties; for the conveni description, these units are collectively referred to as counties.  The data are from the China County Statistical Yearbook, China Urban Statistical Yearbook, Jiangsu Province Statistical Yearbook, statistical yearbooks of prefecture-level cities and some existing districts and counties, as well as the regional national economic and Social Development Statistical Bulletin. The global average production data of each agricultural product come from FAO and Geo remote sensing data.

Three-Dimensional Ecological Footprint Model
EF is an indicator based on the consumption of EF-related products. There are many EF accounting methods, mainly from the perspective of production (EF = EF of local production + EF of import product conversion-EF of export product conversion) and from the perspective of consumption. The former is more suitable for ecological footprint accounting at a national scale and with detailed import and export data, while the latter is more suitable for EF accounting at small scales such as cities and counties. Based on this, the Land 2022, 11,1997 5 of 20 EF accounting in this paper is from the consumption of EF accounting; the EF of arable land, forestland, grassland, and fishing land can be accounted for by the following equation: where YF i is the yield factor for product i, used to convert the local area to the global average production area, and EQF i is the equivalence factor for product i, used to convert this land type to the global average land type area (global hectare). The built-up footprint is accounted for by transferring its direct land occupation (arable land) into EF. The carbon footprint is accounted for by transferring the necessary carbon uptake land (forestland) into EF. Different from the calculation of EF, the calculation of BC is based on how much productive land there is in the area, converted into the size of a unified unit. The calculation method is as follows [41]: The EF accounting step: firstly, A i is the bioproductive area for the production of product i. Figure 2 demonstrates the conversion from 2D EF to 3D EF. As designed by Niccolucci et al. [16], the 3D EF model is defined as follows: duction + EF of import product conversion-EF of export product conversion) and from the perspective of consumption. The former is more suitable for ecological footprint accounting at a national scale and with detailed import and export data, while the latter is more suitable for EF accounting at small scales such as cities and counties. Based on this, the EF accounting in this paper is from the consumption of EF accounting; the EF of arable land, forestland, grassland, and fishing land can be accounted for by the following equation: where is the yield factor for product , used to convert the local area to the global average production area, and is the equivalence factor for product , used to convert this land type to the global average land type area (global hectare). The built-up footprint is accounted for by transferring its direct land occupation (arable land) into EF. The carbon footprint is accounted for by transferring the necessary carbon uptake land (forestland) into EF. Different from the calculation of EF, the calculation of BC is based on how much productive land there is in the area, converted into the size of a unified unit. The calculation method is as follows [41]: The EF accounting step: firstly, is the bioproductive area for the production of product . Figure 2 demonstrates the conversion from 2D EF to 3D EF. As designed by Niccolucci et al. [16], the 3D EF model is defined as follows: Figure 2. From 2D EF to 3D EF. Figure 2. From 2D EF to 3D EF.

LMDI
The logarithmic mean Divisia index (LMDI) has two forms, i.e., the additive form and the multiple form. The additive form and the multiple form are equivalent, and the former is usually used in physical units while the latter is used in indexes for convenient expression [36]. Since this paper intends to compare the difference in the amount of EFD data decomposition between different years, the decomposition in the additive form is more suitable for the research. According to IPAT and its improvements, this paper selected P to represent the population size, G to represent the total GDP, G s to represent the total GDP of the tertiary industry, and T A to represent the total number of patents to use in the LMDI decomposition.A = G/P obtains the per capita GDP, indicating the affluence; S = Gs/G obtains the proportion of GDP of the tertiary industry, indicating the industrial structure; T s = T A /G s obtains the regional patent authorization per unit of GDP of the tertiary industry as a whole in that year, representing the technological intensity of the tertiary industry; and E t = EF depth /T A represents the size of the EF depth per unit of patent authorization, representing the intensity of scientific and technological demand for Land 2022, 11,1997 6 of 20 pollution control. The PyLMDI [42] was used to calculate the results in this study. The equation of the LMDI is as follows: The additive form decomposition for this article is as follows: For a general explanation of those sectors, ∆EFD X is expanded below: X could be P, A, S, T s , or E t in this article, and i is the economic department.

Result of Jiangsu's Ecological Footprint Depth
The results of the figures below were drawn by ggplot2 [43]. The calculation results showed that the EFD increased from 2.10 in 1995 to 6.28  To better reflect the evolution of the EFD in prefecture-level cities and districts and counties under their jurisdiction, the following analysis was carried out at two regional scales.
From the perspective of the prefecture-level city scale, the evolution of the EFD of 13 prefecture-level cities in Jiangsu Province from 1995 to 2015 is shown in Figure 3 and Table A1. The location heights in the figure represent the EFD values of the corresponding city. It can be seen that in the study years, the EFD presents a spatial distribution pattern of high values in the south and low values in the north. In 1995, the regional-level cities showed relatively low differential EFD distributions, among which the EFDs of Nanjing and Wuxi were significantly higher, reaching 4.22 and 4.49, respectively. From 1995 to 2005, the change was small, and the distribution did not change significantly. From 2000-2005, the EFD of the prefecture-level cities underwent relatively significant changes, and Suzhou especially experienced significant growth, increasing 5.71. The main feature of change in 2005-2010 was that the EFD in all cities increased significantly. In 2010-2015, the most obvious change in the EFD was in Lianyungang.
At the county scale, the three-dimensional map of the ecological footprint depth of each county in Jiangsu Province from 1995 to 2015 is shown in Figure 4 and Table A2. The overall spatial correlation also gradually increased in the study years. From the perspective of spatial agglomeration, most of the counties are low-EFD agglomeration areas, which are distributed in the northern and central parts of Jiangsu, and high-EFD agglomeration areas are mainly distributed in southern Jiangsu. Due to the process of urbanization, the population is gradually agglomerating in urban areas, and high values are mainly distributed in urban areas of prefecture-level cities, showing a very high EFD that was significantly higher than that of non-urban territorial units. At the county scale, the three-dimensional map of the ecological footprint depth of each county in Jiangsu Province from 1995 to 2015 is shown in Figure 4 and Table A2. The overall spatial correlation also gradually increased in the study years. From the perspective of spatial agglomeration, most of the counties are low-EFD agglomeration areas, which are distributed in the northern and central parts of Jiangsu, and high-EFD agglomeration areas are mainly distributed in southern Jiangsu. Due to the process of urbanization, the population is gradually agglomerating in urban areas, and high values are mainly distributed in urban areas of prefecture-level cities, showing a very high EFD that was significantly higher than that of non-urban territorial units.  From 1995 to 2000, the EFDs of a small number of counties decreased, namely, the EFDs of Jiangyan City, Taixing City, Changshu City, Jingjiang City, Peixian County, Ganyu County, and Liyang City. The EFDs of Donghai County, Xuyi County, and Sihong County were unchanged at 1 because their EFs were less than the corresponding BC for each county. During this period, 12 counties increased their EFD by 1, and 7 were prefecture-level urban areas, namely, Nanjing Urban, Nantong Urban, Taizhou Urban, Xuzhou Urban, Suzhou Urban, Wuxi Urban, and Huai'an Urban Areas. Others were economically developed districts and counties in southern Jiangsu: Zhangjiagang City, Jiangyin City, Taicang City, Wuzhong District, and Wujiang District. In addition, counties such as Jianhu County, Kunshan City, and Tongshan District had a faster EFD growth rate.
From 2000 to 2005, there were 34 counties whose EFD growth changed by more than From 1995 to 2000, the EFDs of a small number of counties decreased, namely, the EFDs of Jiangyan City, Taixing City, Changshu City, Jingjiang City, Peixian County, Ganyu County, and Liyang City. The EFDs of Donghai County, Xuyi County, and Sihong County were unchanged at 1 because their EFs were less than the corresponding BC for each county. During this period, 12 counties increased their EFD by 1, and 7 were prefecture-level urban areas, namely, Nanjing Urban, Nantong Urban, Taizhou Urban, Xuzhou Urban, Suzhou Urban, Wuxi Urban, and Huai'an Urban Areas. Others were economically developed districts and counties in southern Jiangsu: Zhangjiagang City, Jiangyin City, Taicang City, Wuzhong District, and Wujiang District. In addition, counties such as Jianhu County, Kunshan City, and Tongshan District had a faster EFD growth rate.
From 2000 to 2005, there were 34 counties whose EFD growth changed by more than 1, most of which were distributed in urban areas and counties in southern Jiangsu. The higher changes were mainly Zhenjiang urban area (+8.74), Jiangyin City (+10.53), Nanjing Urban Area (+13.97), Xuzhou Urban Area (+15.69), and Nantong Urban Area (+16.39). Some districts and counties experienced little change, such as Danyang City, Jurong City, Jinhu County, Xuyi County, and Donghai County.
From 2005 to 2010, the growth slowed significantly. Some counties that grew rapidly in the first half of the study period hardly changed or even decreased in this range. These include Xuzhou Urban Area and Nantong Urban Area. Most of the districts and counties with large EFD growth are mainly concentrated in the central and southern regions, including Liuhe District, Suzhou Urban Area, Ganyu County, Zhangjiagang City, Taicang City, and Lishui District.
From 2010 to 2015, the overall growth of the EFD in each county further slowed, and 12 counties experienced negative growth, e.g., Jiangdu District, Nanjing Urban Area, Jianhu County, Taixing City, and Liuhe District. Figure 5 shows the distribution of changes in the overall EFDs, and most of them were between 0 and 0.5 for this period. The counties with high values higher than 1 include Lianyungang Urban Area and Tongzhou District. Other high values are obviously smaller than the changes between 2005 and 2010.

LDMI Results of EF Depth
The following study analyses the influence of the five factors-, , , , andon the EFD changes in Jiangsu Province in each time period from the provincial, prefecture-level city, and county levels. We decompose the EFD changes under each influencing factor. Figure 6 shows the decomposition of the overall EFD change in Jiangsu Province with LMDI. It can be seen from the figure that in each study period had the main negative effect on the change in EFD, and the and factors both showed a positive effect on

LDMI Results of EF Depth
The following study analyses the influence of the five factors¯A, E t , P, S, and T s -on the EFD changes in Jiangsu Province in each time period from the provincial, prefecturelevel city, and county levels. We decompose the EFD changes under each influencing factor. The decompositions of EFD changes in 13 prefecture-level cities are shown in Figure  7. And the summary of factors' decomposition for 13 prefecture-level cities of Jiangsu is shown in Table A3. The decompositions of each factor from 1995 to 2000 were relatively balanced among the prefecture-level cities. It is obvious that the main negative influence of EFD at this stage was (mean: −3.186); although had a negative effect in Nanjing (−1.717), its overall positive effect was the largest among these factors, and a positive effect was also found for , , , which had a positive impact on most of these cities. From 2000 to 2005, the situation changed variously compared with that in the previous five years: Factor became the main positive influencing factor, and , , and showed a relatively uniform positive and negative influence on the number of cities affected; that is, the negative effect became stronger. The factor was still the most negative influencing factor, showing a positive impact in three cities, namely, Changzhou (0.337), Suzhou (1.638), and Nantong (0.03). It can be seen that on the whole, these 13 cities in the time period also had the largest change in EFD compared with other time periods, but the overall variation was small, mainly due to the overall decrease in the absolute value of negative factors. From 2005 to 2010, the absolute value of the influence of various factors increased, especially in Changzhou, Nanjing, Nantong, Suzhou, and Wuxi, cities with relatively developed economies in Jiangsu. The main positive influencing factors were , , and , and the influence of EFD was positive in all cities.
was the most important positive influencing factor, and P had a positive impact on Suzhou (4.082), Nanjing (1.371), and Wuxi (1.468), which was due to the positive population growth (immigration) of these The decompositions of EFD changes in 13 prefecture-level cities are shown in Figure 7. And the summary of factors' decomposition for 13 prefecture-level cities of Jiangsu is shown in Table A3. The decompositions of each factor from 1995 to 2000 were relatively balanced among the prefecture-level cities. It is obvious that the main negative influence of EFD at this stage was E t (mean: −3.186); although T s had a negative effect in Nanjing (−1.717), its overall positive effect was the largest among these factors, and a positive effect was also found for P, S, A, which had a positive impact on most of these cities. From 2000 to 2005, the situation changed variously compared with that in the previous five years: Factor A became the main positive influencing factor, and P, S, and T s showed a relatively uniform positive and negative influence on the number of cities affected; that is, the negative effect became stronger. The E t factor was still the most negative influencing factor, showing a positive impact in three cities, namely, Changzhou (0.337), Suzhou (1.638), and Nantong (0.03). It can be seen that on the whole, these 13 cities in the time period also had the largest change in EFD compared with other time periods, but the overall variation was small, mainly due to the overall decrease in the absolute value of negative factors. From 2005 to 2010, the absolute value of the influence of various factors increased, especially in Changzhou, Nanjing, Nantong, Suzhou, and Wuxi, cities with relatively developed economies in Jiangsu. The main positive influencing factors were Ts, A, and S, and the influence of EFD was positive in all cities. T s was the most important positive influencing factor, and P had a positive impact on Suzhou (4.082), Nanjing (1.371), and Wuxi (1.468), which was due to the positive population growth (immigration) of these cities during this period. The largest negative factor in each city was E t . From 2010 to 2015, the absolute value of the overall influencing factors decreased, and the most important positive influencing factor became A, followed by S. P was positive except in Yancheng (−0.015), but the influence became very small (mean: 0.108); T s and E t both showed positive and negative influences in 13 cities; the former mainly had a positive influence, and the largest negative influence factor was still E t .  From 1995-2000, there were three counties with EFDs equal to 1 because each EF was smaller than the corresponding BC, i.e., Donghai County, Xuyi County, and Sihong County. There was no change during this period for the three counties, and the figure shows that the factor change was 0. There were three change factors with a more obvious decomposition effect, i.e., A, S, and P. A showed a stronger positive effect, and the decomposition values were all positive. A had the largest average value of the decomposition amount of the EFD change at this stage, which is 1.422. S mainly showed a positive decomposition effect, with an average decomposition of EFD changes of 0.439. P generally had positive effects, with an average decomposition value of 0.225. At this stage, there were only 16 counties for P with a negative effect, and the absolute value of the number was small. In addition, T s showed more negative effects, and E t showed both positive and negative effects.

Discussion
A discussion of EFD differences and the factor decomposition effects at different scales is necessary. From the research results at the provincial, prefecture-city level, and district and county levels, it can be seen that when different scales are unfolded layer by layer, the finer the geographical units are, the greater the difference in EFD changes and factor decomposition. The finer the regional scale is, the greater the heterogeneity. The data at the provincial level and the prefecture-level city scale show little change, which is mainly reflected in the different distributions between the north and south of the prefecture-level cities. The county-level data not only showed the difference between the north and the south but also highlighted the significant differences between the urban areas and nonurban areas. In terms of the absolute value of the factor decomposition of each county, the urban areas are also significantly higher than the value of the non-urban areas. Due to the agglomeration effect of urbanization population, economy, and other factors in large cities, the population density of prefecture-level urban areas was significantly higher than that of surrounding districts and counties, and the population size gap between regions was large. There were some differences in the impact of population factors at different scales: from the perspective of 20-year temporal and spatial changes, the overall impact of population factors on EFD in Jiangsu Province decreased, while from the perspective of prefecture-level cities and counties, especially from 2005 to 2010, the factors affecting the population increased in some regions. This was due to the migration of the population between regions brought about by urbanization, especially the migration of the population from districts and counties to prefecture-level cities at the same level. This finding shows that urbanization will affect the ecological pressure of developed regions to a certain extent.
From the perspective of counties, the results of the regression analysis to study the relationship between EFD and explanatory variables can be explained from the regression coefficient, and the decomposition method mainly gives the decomposition value of the change in EFD caused by the change in decomposing factors. The correlation analysis of each decomposition factor and EFD is shown in Figure 9. It can be seen that in the research years, population size, affluence, industrial structure, and depth of ecological Land 2022, 11,1997 13 of 20 footprint mainly showed a significant positive correlation, while the relationships between the level of science and technology in the tertiary industry, the intensity of scientific and technological demand for pollution control, and EFD were diverse. As defined in Equation (7), for EF as a whole to input into the equation, the direction of change of positive and negative change in the same direction is shown in Figure 10, and the points are all located in the first and third quadrants. The relationship of the change in each factor and the change in EFD and the relationship between each factor and EFD show similar characteristics: the respective differences in population size, affluence, and industrial structure have a strong correlation with the corresponding decomposition of EFD in each time stage, while ∆T s and ∆E t have a weak relationship with the corresponding decomposition amount of EFD. This weak relationship presents a variety of relationships at the county scale. It plays a relatively important decomposition role in some counties, and the influence of the EFD was less stable. In some cases, EFD has a positive impact, and in some cases a negative impact.  From the comparative discussion of the literature findings, there is no EF literature from the perspective of EFD in Jiangsu Province, and there are findings in Jiangsu Province mainly from the provincial scale EF [44,45]. The results of the previous studies indicate the overall trend of increasing ecological pressure in Jiangsu province, and the results of this study are clearer about the ecological pressure distribution and driving factors. In analysis of the results of other related 3D EF studies [22,25], most of the EFD results have a larger scale of geographic units. Comparative analysis shows that this study is better able to show the 3D EF of each regional unit and to respond to spatial information, and it is easy to carry out the evaluation of ecological pressure and ecological balance within the city. The results of this study differ significantly from other studies in terms of the effects of demographic factors on EFD [25], partly due to model variability and partly due to variability between regions. The analysis of GDP (GDP per capita) results [20,21] has similar conclusions to most studies and serves as a major positive driver, showing more clearly the difference in the impact of GDP on EFD over the study year range than most studies. In conclusion, there are differences in other factors, and the overall difference in the study in terms of influence factors and interactions is obvious in time and space between regions, highlighting the degree of subdivision of the influencing factors in different regions and different periods.

Conclusions
In this article, the focus points are the distribution and changes in Jiangsu's EFDs and the decomposition of the EFD changes with influencing factors. The following conclusions are drawn: (1) Multi-scale EFD research in Jiangsu Province can mine ecological information at different geographic scales and compare differences and scales, which can provide support for the study of ecological balance between and within regions. At the county scale, the unbalanced north-south difference in the EFD distribution increases year by year, and the difference between the urban areas and the counties is obvious. (2) The factor decomposition method divides the factors of each time interval of EFD in Jiangsu Province. Affluence was the most important decomposing factor, which has always been shown to affect the change in EFDs at different scales, different regions, and different time periods. Other factors all showed the effect of diversity and quality in the above three conditions, among which the factor of the ratio of EFD Land 2022, 11, 1997 15 of 20 to technological level was the main inhibitory factor. The annual factor changes have certain differences in geographical units at different scales. Industrial structure and population factors mainly showed the promoting effect of EFD, while the technological intensity of the tertiary industry has a relatively large heterogeneous effect on EFD. (3) Finally, this paper discusses the relationship between factors in a targeted manner and proposes countermeasures and suggestions to reduce EFD and balance the ecological pressure among regions from multiple perspectives.

Policy Implications
In view of the results of EFD and factor zoning in Jiangsu Province, to reduce EFD and balance the ecological pressure among regions, the following countermeasures are suggested: (1) The ecological pressure analysis in Jiangsu Province should be conducted at multiple scales, focusing on the balanced distribution of the population and the balance of ecological pressure. To this end, it is necessary to pay attention to the balanced development among regions. The transfer and layout of industries from developed counties with high ecological pressure to economically underdeveloped districts and counties with low ecological pressure would help transfer the impact of population factors on EFD and reduce the original high EFD value. The specific implementation can be combined with the research results (Tables A1 and A2) in this paper for the reference of policymakers, e.g., in 2015, the ecological pressures of Suzhou, Nanjing, Wuxi, Changzhou, Nantong, and Zhenjiang, most of which were in the south of Jiangsu, were significantly higher than those of other prefecture-level cities. Therefore, the industries and population of those cities could be transferred to other cities with small EFDs. At the same time, attention should be given to the integrated urbanization development of prefecture-level cities and districts and counties to increase the city's overall ecological pressure resistance capability and avoid extreme increases in local ecological pressure caused by agglomeration. Therefore, it is necessary to reasonably optimize their urbanization rates, and to strengthen the construction of rural basic livelihoods, so as to alleviate urban ecological pressures. Within a prefecture-level city, taking Huai'an City for example in 2015, the EFD of the Huai'an Urban reached 24.44, while in the same period, the EFD of non-urban areas under the prefecturelevel city was only in the range of 1.68-3.71. Therefore, the high EFD industries in Huai'an Urban Area need to be balanced with other regions, while expanding the attractiveness of non-urban areas within the prefecture-level city and reducing the pressure on the urban area. (2) It is recommended to focus on different growth-influencing factors at different stages, because the influencing factors have a certain change pattern at different stages. For example, the population factor has a greater impact in the early stage but has less influence in the later stage. In addition, the heterogeneity of critical influencing factors of EFD among districts and counties should be highlighted to develop differentiated countermeasures to reduce EFD. From the results of the analysis, for example, in 2015, to effectively refrain the growth of EFD, Huai'an Urban should focus on the increase in T s , while Nanjing Urban should focus on the increase in E t . (3) The results show that, for most of the regions in Jiangsu, attention should be given to the development mode of the green economy and to optimizing the industrial structure to reduce the impact of economic growth and the economic structure on EFD. In addition, a focus should be on improving the GDP output of the tertiary industry per unit of science and technology and promoting the speed of science and technology development to exceed the growth rate of EFD, which will help reduce EFD at the scales of provinces, cities, districts, and counties.

Limitations of this Study and Future Study Recommendations
Limitations of this study are the following: (1) Since there is no district/county level input-output table, this paper does not use LMDI to decompose EFD influencing factors based on industry sectors. (2) Due to the fact that this study focuses on the decomposition of EFDs with LMDI, this paper did not conduct a comparative study among multiple decompositions or regression models. (3) There was no analysis on EFDs by land type.
Future study recommendations are the following: (1) It can be suggested that the factor research of EFD can be compared based on LMDI, regression, and other models.
(2) The spatial agglomeration effect and heterogeneity of these influencing factors of LMDI can be further studied. (3) The EFDs of different land type could be analyzed further.