Research on Expansion Characteristics of Aquaculture Ponds and Variations in Ecosystem Service Value from the Perspective of Protecting Cultivated Lands: A Case Study of Liyang City, China

Abstract

In the context of global food insecurity, a large amount of cultivated land in China has been occupied by aquaculture ponds, leading to a series of variations in the ecological environment. The Chinese government pays close attention to the problem. In order to achieve sustainable development and ensure the safety of China’s cultivated land, the paper uses Liyang City as an example to discuss the spatial characteristics of the expansion of aquaculture ponds through occupying cultivated lands and analyzes the variations in ecosystem service value and cultivated land function. The conclusions are as follows: (1) 2073.24 hectares of cultivated lands were occupied for expanding aquaculture ponds in Liyang from 2009 to 2019, and there was a small number of new aquaculture ponds in the ecological protection area, which shows that the aquaculture ponds in Liyang City are at the stage of disorderly expansion; (2) the total value of ecosystem services increased by 1.43%; supply and support services values decreased, but the increase in regulation and cultural services values was sufficient to more than compensate for the mentioned losses; and (3) the expansion of aquaculture ponds leads to a decrease in the carbon storage of cultivated land, which in turn has negative impacts such as an increase in atmospheric carbon concentration.

1. Introduction

In May 2021, the Food and Agriculture Organization (FAO) of the United Nations and the World Food Programme (WFP) issued the 2021 Global Report on Food Crisis. According to the report, the problem of food insecurity reached its maximum level in the previous five years and the increasing seriousness of the food crisis was influenced by COVID-19 and the international situation. Cultivated land is not only the main source of food production [1], but also an important guarantee for the steady development of the social economy [2,3]. Although cultivated land is a renewable resource [4], its recovery cycle is very long. Therefore, when a large amount of cultivated land is occupied by other land types, it will undoubtedly directly threaten national food security. With the development of urbanization and industrialization in China, cultivated lands experienced a rapid decline. Later, China successively issued the ecological civilization construction strategy and rural revitalization strategy. Under dual restraints of the ecological red line and the red line of cultivated land protection, the decline in the area of cultivated lands was reduced to a certain extent, but the amount of cultivated lands used for grain continued to decrease. The fact was that the phenomenon of non-grain production of cultivated land occurred with a structural adjustment in the agricultural industry [5,6]. The “non-grain production of cultivated land” refers to the behavior of planting non-food crops on cultivated lands, for example, planting economic crops, excavating aquaculture ponds and afforestation. All these behaviors will affect the quality of cultivated land and may threaten the security of cultivated land and food in cases of disordered development. The behavior of excavating aquaculture ponds has been widely considered as an important cause of the decrease in cultivated lands [7]. Increasingly more cultivated lands are occupied by aquaculture ponds [8]. It is estimated that globally, about 1.5 million hectares of land is occupied by aquaculture every year [9]. The rapid development of aquaculture ponds depends on occupying large areas of high-quality cultivated land. Although it has created huge economic benefits, it also has some negative influences on the ecological environment [10,11,12,13,14]. In this context, the paper quantitatively evaluates the spatial distribution and variations in the characteristics of aquaculture ponds, which is important and significant for the protection of regional cultivated lands and the sustainable development of agriculture.

As one of the main utilization types of terrestrial ecosystems, cultivated lands play an important role in ensuring human and food security, and variations in cultivated land have an important influence on the regional ecological environment [15]. After issuing the reform scheme of rural land circulation in 2008, China saw rapid development in rural land circulation, the utilization type of cultivated lands also emerged [16], and transformation from cultivated lands into aquaculture ponds gradually became a prominent utilization mode in the transformation of cultivated lands. Some studies indicate that the expansion of aquaculture ponds in the Yellow River Delta was achieved through occupying cultivated lands [17]. Some other scholars found that a notably large amount of cultivated land was transformed into aquaculture ponds through their analysis of the spatial variation in aquaculture ponds in the coastal area of China [18]. Due to its speed and low cost [19], remote sensing data is widely used in studying aquaculture ponds. Some scholars supervised and classified the land cover of aquaculture ponds through remote sensing data [20] and some other scholars plotted the distribution of aquaculture ponds in the coastal area of the East China Sea through TM imagery [21]. In addition, different scholars studied aquaculture ponds from different perspectives, for example, hazards of aquaculture ponds in Southeast Asia in closed water areas [22]; expansion of aquaculture ponds along the coastline of China [23]; spatial distribution of aquaculture ponds along the coastline in India [24]; the influence of aquaculture on the ecological environment along the coastline in Bangladesh [25]; etc. In terms of research methods, many scholars tend to recognize the spatial variation in aquaculture ponds based on the remote sensing satellite data [26] or automatically extract positions of aquaculture ponds from Google image data [27]. Some scholars studied the driving factors of occupying cultivated lands for the expansion of aquaculture ponds in Guangdong Province through questionnaire data [28]. Some scholars analyzed land use variations in the coastal area due to aquaculture ponds through remote sensing (RS) and geographic information system (GIS) [29]. Some scholars analyzed the degree of influence of the expansion of aquaculture ponds on changes in local land cover through the landscape pattern index [9].

Another emphasis in aquaculture pond research is problems in the ecological environment due to the expansion of aquaculture ponds. Many studies indicate that the expansion of aquaculture ponds can promote local economic development but, to a certain extent, will also have a negative influence on the local ecological environment [30,31]. The long-term negligence of ecological value and social value of cultivated land resources is the largest shortcoming and flaw in protecting cultivated lands in China [32]. The rapid development of aquaculture ponds through occupying cultivated lands leads to further degradation of the ecological environment. Some scholars studied the problem from the view of ecosystem service. The ecosystem service concept first came into view in the 1980s [33]. The ecosystem service refers to life support products and services that are obtained directly or indirectly through the structure, process and function of the ecosystem [34]. The ecosystem service is the material basis and basic condition for the survival and development of humans and is the key natural capital of humans [35]. The ecosystem can provide necessary food and raw materials for humans and can support non-material services such as biological and chemical circulation and carbon sequestration [36]. The decisions of policy-makers are often influenced by economic data. Since ecosystem services cannot be directly measured by economic data, they must be quantified in economic terms based on their price in the market or the price of alternative goods and services [37]. Ecosystem service value (ESV) is an important quantitative and evaluation method for the function of ecosystem services [38,39]. The change in ESV is related to factors such as land use change, climate change and biodiversity [40], among which land use change is the most direct reflection of human activities. Therefore, ESV measurement methods are widely used in related research on land use change, such as the impact of land use change on the ecological environment in the process of urbanization [41,42]. In addition, ESV has good practical application value in the selection of regional sustainable development strategies [43], urban planning and management [44] and the formulation of ecological compensation policies [45].

In the context of paying close attention to the ecological environment and protection of cultivated lands in China, evaluating the influence of transformation between aquaculture ponds and cultivated lands on ESV is important and significant for the resource protection and sustainable development of cultivated lands. The paper aims to understand the spatial characteristics of the expansion of aquaculture ponds through occupying cultivated lands in the countryside of China and to analyze the influence of the process on local ESV and the function of cultivated lands. The paper takes Liyang City (Jiangsu Province) as an example to investigate the objectives of: (1) analyze the spatial characteristics and quantitative relationship of the expansion of local aquaculture ponds through occupying cultivated lands from 2009 to 2019; (2) analyze local ESV change from the perspective of service function and land use type, and analyze the value flow of the ecosystem service due to the expansion of aquaculture ponds through occupying cultivated lands; and (3) use the carbon capacity to evaluate the influence of the expansion of aquaculture ponds through occupying cultivated lands on the function of cultivated lands to provide a reference for promoting the sustainable utilization of regional resources of cultivated lands.

2. Materials and Methods

2.1. Overview of Research Area

Located in eastern China, Liyang City is a county-level city in the southwest of Jiangsu Province, as shown in Figure 1. With the northern latitude 31°09′~31°41′ and eastern longitude 119°08′~119°36′, Liyang has a length of 59.06 km from north to south, width of 45.14 km from east to west and a total land area of 1535 km². Located at the edge of the northern subtropical region, Liyang has a maritime, moist and monsoon climate and significant monsoon characteristics. The prevailing wind direction throughout the year is the easterly wind, annual average temperature is 17.5 °C and annual average precipitation is 1149.7 mm. Liyang City is located in complex topography of mountainous land, hills, plains and polder areas. It has plains from west to east with flat terrain, and mountainous lands and hills in the west, north and south.

Liyang City is a nationally advanced county in food production. However, many cultivated lands are uncultivated or occupied by aquaculture ponds, and large areas of non-food crops have been cultivated in recent years (Figure 2a). The per capita area of cultivated land decreased from 1.05 mu to 0.96 mu, which is far lower than the national per capita level of 1.51 mu. Although Liyang has gained the honorable title of “National Ecological City”, it still faces many challenges in the ecological environment; the problem of organic matter and soil pollution is gradually increasing; the situation of environmental water protection is serious (Figure 2b) and direct drainage of agricultural wastewater is still prominent; and the ecological environment indicates a significant trend of environmental degradation, and fragmentation of ecological space tends to be serious.

2.2. Data Source and Processing

The 2009 land use data for Liyang City were taken from The Second China’s Territory Survey Data, while the 2019 land use data and boundary data of administrative districts were taken from The Third China’s Territory Survey Data. The vector data of the ecological protection area were provided by the Liyang Ecological Environment Bureau. The cultivated area of food crops, yield and dosage of chemical fertilizer and pesticide were taken from the Liyang Statistical Yearbook. The national average price of food crops was taken from the Collection of National Cost Benefits of Agricultural Products in 2019.

2.3. Research Methods

2.3.1. ESV Evaluation

The proposed value coefficient by Costanza et al. [38] has been widely used in research on ESV, but the value coefficient is not universally agreed. To better analyze the value variation in ecosystem services in the research area, the study assigns value equivalence for ecosystem services on cultivated lands and aquaculture ponds in the research area based on existing research findings [46] and the amended “Table of value equivalence for ecosystem service per unit area for China’s land ecosystem” by Xie Gaodi et al. [47] based on actual conditions in China (Table 1); and calculates ESV in the research area through Formulas (1) and (2):

(1)$\mathrm{ESV}={\displaystyle \sum }_{i=1}^n{A}_i\times V{C}_i$

(2)$V{C}_i={\displaystyle \sum }_{j=1}^{k}E{C}_j\times E_a$

where, ESV is the ecosystem service value, Yuan/a; i is the land use type; j is the service type of the ecosystem; $A_i$ is the area of type i land use, ha; VC_i is the ESV per unit area of type i land use, Yuan/hm⁻² a⁻¹; EC_j is the value equivalence of ecosystem service for type j land use; Ea is the economic value of ecosystem service per unit, Yuan/hm⁻² a⁻¹.

To accurately determine the economic value of ecosystem services per unit in the research area, the study corrects related parameters through Formula (3) [48]. The obtained the economic value of natural food production on farmland in the research area in 2019, which is the economic value of ecosystem serviced per unit, is RMB 4502.99 Yuan/hm⁻² a⁻¹.

(3)$E_a=1/7{\displaystyle \sum }_{i=1}^n\frac{m_i{q}_i{p}_i}{M}\times MCI$

where, Ea is the economic value of ecosystem services per unit in the research area, Yuan/hm⁻²a⁻¹; i is the type of food crop; $m_i$ is the national average price of type i food crop, Yuan/kg; $q_i$ is the cultivated area of type i food crop, hm²; $p_i$ is the per unit yield of type i food crop, kg/hm²; M is the total cultivated area of food crop, hm²; MCI is the multi-cropping index of farmland in the research area, set at 2.5.

2.3.2. Value Flow Analysis for Ecosystem Services

The value flow analysis of ecosystem services calculates profit and loss in ESV due to mutual transformation between different land use types through transformation data of land use types. The purpose is to analyze the influence of land use transformation on ESV [49]. The calculation is shown in Formula (4):

(4)$P{L}_{ij}=\left(V{C}_j-V{C}_i\right)\times A_{ij}$

where, $P{L}_{ij}$ is the profit and loss in ESV after transforming type i land use into type j land use; $V{C}_j$ and $V{C}_i$ are value coefficients of ecosystem services for type i land use and type j land use; $A_{ij}$ is the area of transforming type i land use into type j land use.

2.3.3. Evaluation Method for Functions of Cultivated Lands

The expansion of aquaculture ponds through occupying cultivated lands and drainage of aquaculture wastewater will, to a certain extent, have a negative influence on the quality of cultivated lands. Due to a lack of soil data, quantitative evaluation cannot be directly made for soil texture. Hence, the study uses the carbon capacity of cultivated lands to evaluate the functional variation in cultivated lands.

The soil organic carbon is an important index for representing the variation in soil fertility, which influences the physical, chemical and biological characteristics of soil and lays an important foundation for the stable production of food crops and sustainable development of agriculture. The carbon capacity of different types of land use is represented by the product between carbon density of soil and the area of related land use type, as shown in Formula (5):

(5)$C{S}_i=C{D}_i\times A_i$

where, $C{S}_i$ is the carbon capacity of type $i$ land use, kg; $C{D}_i$ is the carbon density of type $i$ land use, kg/m²; $A_i$ is the area of type $i$ land use, m². In Jiangsu Province, the carbon density (0~100 cm) is 9.29 kg/m² for cultivated lands and 8.11 kg/m² for aquaculture ponds [50,51].

3. Results

3.1. Characteristic Analysis for Cultivated Lands and Aquaculture Ponds

3.1.1. Quantitative Variation in Aquaculture Ponds

The agricultural area in Liyang City was 59,797.42 hectares in 2009 but decreased by 14,887.47 hectares to 44,909.95 hectares in 2019, where 6279.21 hectares of cultivated lands was occupied by aquaculture lands (new) and 4205.97 hectares of aquaculture ponds was transformed into cultivated lands (decrease in aquaculture ponds). In fact, 2073.24 hectares of cultivated lands was occupied for expanding aquaculture ponds from 2009 to 2019.

The increased and decreased areas of aquaculture ponds in various towns are different. See Table 2 for detailed changes. The area of aquaculture ponds increased by 3675.99 hectares in Licheng Town, Shanghuang Town, Tianmu Lake Town, Shangxing Town and Shezhu Town; but decreased by 1602.75 hectares in Daitou Town, Bieqiao Town, Zhuze Town, Nandu Town and Kunlun Street. The greatest increase in aquaculture pond area occurred in Shezhu Town with about 2478.99 hectares, which was 39.48% of the total increased areas of aquaculture ponds in the city. The smallest decreased in aquaculture pond area was in Bieqiao Town with about 890.00 hectares, being 21.16% of the total decreased areas of aquaculture ponds in the city.

3.1.2. Spatial Expansion of Aquaculture Ponds

As shown in Figure 3, new aquaculture ponds in Liyang City are mainly distributed in Shezhu Town and Shangxing Town in the southwest, as well as Shanghuang Town in the northeast corner of Liyang City. Although significant and concentrated distribution is not found in the other towns, some small and fragmented spots can be seen in each town. The decreased aquaculture ponds are mainly distributed in the north of Liyang City, including Daitou Town, Bieqiao Town, Zhuze Town and Nandu Town, and some aquaculture ponds can also be found in the southern Daibu Town and Shezhu Town.

Generally, the ecological protection area is the area that requires mandatory protection due to ecological sensitivity [52]. According to laws and regulations in China, development and productive construction activities are not allowed in ecological protection areas. As shown in Figure 3, however, some new aquaculture ponds are located in ecological protection areas. After calculation, there are 540.28 hectares of new aquaculture ponds in ecological protection areas, indicating a disordered expansion of aquaculture ponds in Liyang City.

3.2. ESV Change Analysis

3.2.1. ESV Change

The expansion of aquaculture ponds through occupying cultivated lands leads to ESV change in Liyang City. From 2009 to 2019, ESV increased from RMB 13.821 billion Yuan to RMB 14.019 billion Yuan and total EVS increased by RMB 198 million Yuan.

The dominant type of ecosystem service function is regulation services. From 2009 to 2019, the regulation service value was over 90% of the total ecosystem service value and increased by RMB 293 million Yuan; the cultural service value increased by RMB 4.19 million Yuan; the supply service value and support service value decreased by RMB 92 million Yuan and RMB 8 million Yuan, respectively. See Table 3 for detailed changes.

From the view of the ecological value of land use type (Table 4), the ESV of cultivated lands is higher than aquaculture ponds. In the research period, the ecological service value decreased by RMB 437 million Yuan for cultivated lands and increased by RMB 691 million Yuan for aquaculture ponds; but the regulation service value decreased the most for cultivated lands, up to RMB 307 million Yuan, followed by the supply service value, about RMB 87 million Yuan; the decreased cultural service value was lowest, about RMB 6 million Yuan. For aquaculture ponds, the regulation service value increased the most, up to RMB 653 million Yuan, but the supply service value increased the least, only RMB 8 million Yuan. Hence, the expansion of aquaculture ponds through occupying cultivated lands can lead to a decrease in supply service value and support service value but an increase in regulation service value and cultural service value in the ecosystem in Liyang City.

3.2.2. Value Flow of Ecosystem Service

From 2009 to 2019, due to the occupation of cultivated lands by aquaculture ponds, the ESV increased by RMB 1.03 billion Yuan in Liyang City, where the ESV increased by RMB 3.13 billion Yuan for the transformation of cultivated lands into aquaculture ponds and decreased by RMB 2.1 billion Yuan for the transformation of aquaculture ponds into cultivated lands. The main function leading to ESV change is the regulation function, of which there is a 95.92% change in total ESV in Liyang City.

From the view of ESV, expansion of aquaculture lands through occupying cultivated lands can improve the total ESV of the region. From the view of protection of cultivated lands and food security, the expansion of aquaculture lands through occupying cultivated lands will firstly lead to a quantitative decrease in cultivated lands and then lead to a quality decrease in cultivated lands. In the context of global food insecurity, the problem of protecting cultivated lands cannot be neglected.

3.3. Analysis of Functional Change in Cultivated Lands

Carbon emissions come mainly from human activity, including combustion of fossil fuels and land use change [53]. Different types of land use have different carbon storage capacities. The transformation of land use from a type with high carbon storage capacity into a land use type with low carbon storage capacity generally leads to carbon emissions into the atmosphere [54], which will threaten the regional ecosystem.

In Liyang City, the expansion of aquaculture ponds leads to the rapid discharge of organic carbon. The transformation from cultivated lands into aquaculture ponds leads to a decrease in soil carbon capacity of 8.74 × 10⁷ kg. The transformation from aquaculture ponds into cultivated lands leads to an increase in soil carbon capacity of 5.86 × 10⁷ kg. The expansion of aquaculture ponds leads to a decrease in soil carbon capacity of 2.89 × 10⁷ kg. The reduced emission volume of stored carbon storage into the atmosphere leads to further improvement in carbon concentration in the atmosphere.

4. Discussion

4.1. Causes of Expansion for Aquaculture Ponds

In recent years, the Chinese government has gradually paid close attention to the problem of protecting cultivated lands. In the early stage, many cultivated lands were occupied for construction due to the rapid development of industrialization and urbanization in China. With social and economic development, the disordered expansion of construction lands is effectively controlled and replaced by the expansion of aquaculture ponds. From 2009 to 2019, aquaculture lands rapidly expanded by occupying cultivated lands, and in Jiangsu Province there was a significant increase in the area of aquaculture ponds in Liyang City [18]. Aquaculture ponds have been found to be one of the land use types with rapid growth. Similarly, some other studies also indicate that the main land use type for promoting an increase in aquaculture ponds is cultivated land [17]. Although the expansion of aquaculture ponds can create enormous economic benefits, it also leads to a number of environmental problems [55].

In fact, the earnings of transforming cultivated lands into aquaculture ponds are far higher than grain earnings due to the low earnings of food crops on cultivated lands. According to an interview with local aquaculture households in Liyang City, the earnings per hectare of cultivated land is less than RMB 2895.97~3787.03 Yuan, while the earnings per hectare of aquaculture can reach RMB 119,402.99~149,253.73 Yuan. For this reason, the peasants are not willing to plant food crops and choose to transform cultivated lands into aquaculture ponds. Due to the influence of surrounding peasants [56], increasingly more peasants prefer to transform cultivated lands into aquaculture lands, which is the main reason for the rapid expansion of aquaculture ponds in Liyang City. Because there is unlimited control over cultivated lands in the Land Management Law of the People’s Republic of China, the first-level managers cannot effectively control the transformation from cultivated lands to aquaculture ponds by way of supervision and law enforcement, which is another important reason for the expansion of aquaculture ponds.

4.2. Influence of Expansion of Aquaculture Ponds

At present, two methods can be used to study ESV, including the method based on the equivalence factor of unit area value and that based on the price of unit service functions [48]. The method based on the price of unit service functions requires more accurate ecological data, while the method based on the equivalence factor of unit area value requires only area data of different land use types. In view of data availability, the equivalence factor method is widely used in China. The ESV difference has a strong correlation with the list of used equivalence factors. To obtain more accurate results, the paper corrects the ESV equivalence factor based on the actual conditions of Liyang City. LEI Jun-zheng et al. used the list of equivalence factors of XIE Gaodi and obtained consistent results with this study [49].

In 2009, the aquaculture of freshwater shrimps greatly promoted local economic development and solved employment, to a certain extent, in Shezhu Town. Later, Liyang City considered aquaculture as the dominant type of agricultural industry for eliminating poverty. From 2009 to 2019, the dominant land use transformation in Liyang City was from cultivated lands into aquaculture ponds. Related studies indicate that aquaculture ponds include large areas and rapid growth [57]. The expansion of aquaculture ponds leads to an increase in total ESVs in Liyang City due to a lower ESV per unit area of the cultivated area than aquaculture ponds. Some other studies indicate that an expansion in water domain land will lead to an increase in regional ESV [58]. Although the expansion of aquaculture ponds through occupying cultivated lands can improve regional ESV, it will also have some negative influences, which is neglected in previous studies. Most of the previous studies consider the ecosystem of a region as a whole and study the ecosystem of different land use types as an integral part but neglect differences among different land use types in the ecosystem. Some studies found that cultivated lands would be occupied for the expansion of aquaculture ponds [59], but few studies estimated the ESV change. In consideration of the urgency of the phenomenon, the study selects cultivated lands and aquaculture ponds to independently analyze their ESV changes, which will improve research in the field.

The example of expanding aquaculture ponds through occupying cultivated lands is not specific to Liyang City. Some studies indicate that aquaculture ponds are also expanded through occupying cultivated lands in the Yangtze River Delta of China [17], Sinaloa, Mexico [9] and Thailand [60]. Farmers need to regularly change the water in ponds in order to improve the breeding output. However, the breeding wastewater will often be discharged onto cultivated land near the pond, and the sediment in the wastewater will lead to the degradation of cultivated land, which is ultimately not suitable for grain planting [61]. The transformation from cultivated lands into aquaculture ponds will promote regional carbon emission, which conflicts with existing low-carbon policy and will not promote the long-term sustainable development of the local ecosystem. During the research period, although some aquaculture ponds were transformed into cultivated lands in Liyang City, newly cultivated lands had low grain production capacity [62] and could not maintain their original grain production capacity.

When studying the problem of expanding aquaculture ponds through occupying cultivated lands, we shall consider that it can improve regional ESV and threaten the production capacity of cultivated lands. On the practical level, regional managers shall comprehensively consider the ESV and the protection of cultivated lands based on regional development objectives.

4.3. Limitations and Innovations

The different classification of land use types may lead to different ESV results. Although Costanza [38] and XIE Gao-di [63] provided value equivalence of different land use types, their classifications did not include aquaculture ponds. Hence, the paper uses ESV equivalence of the water area only to substitute aquaculture ponds, but the water domain land includes multiple Grade 2 classifications such as lakes, aquaculture ponds, reservoirs and rivers, which may result in certain deviations from the actual conditions in Liyang City. Due to the lack of soil texture data, the change in cultivated lands can be analyzed only through carbon capacity. Despite the above problem, it is possible to determine the influence of the expansion of aquaculture ponds on ESV at the macroscopic level. Therefore, the paper aims to attract the attention of scholars to a series of problems due to the expansion of aquaculture ponds through occupying cultivated lands.

Most previous studies analyze the spatial distribution of aquaculture ponds along the coastline of a region based on macro-scale or medium-scale remote sensing image data, instead of studying aquaculture ponds and cultivated lands as a whole. In addition, ESV-related studies usually analyze all land use types as a whole, which may overlook regional characteristics. For the above problem, the paper analyzes the expansion of aquaculture ponds through occupying cultivated lands and its ESV change in Liyang City through vector data based on the perspective of a microcosmic city. The paper firstly refined previous studies in the research scale, then used the vector data with higher accuracy than raster data, and finally considered aquaculture ponds and cultivated lands as a whole to offset previous deficiencies.

5. Conclusions

Based on the land use data of Liyang City in 2009 and 2019, the paper analyzes variation in the ESV and function of cultivated lands due to occupation of cultivated lands by aquaculture ponds for expansion and draws conclusions as below:

Large areas of cultivated land is occupied for expanding aquaculture ponds in Liyang City. The increased aquaculture ponds are mainly distributed in the southwest and northeast corners, while decreased aquaculture ponds are mainly distributed in the north of Liyang City. It is noted that new aquaculture ponds in the ecological protection area indicate disordered expansion of aquaculture ponds.

The expansion of aquaculture ponds through occupying cultivated lands will lead to an increase in total ESV in Liyang City, where values of supply services and support services will decrease but values of regulation services and cultural services will increase. The regulatory function is the main factor that results in the value variation in ecosystem services.

Although the expansion of aquaculture ponds through occupying cultivated lands can improve Liyang City’s ESV, it will have a negative influence on cultivated lands. During the research period, the carbon storage of cultivated land in Liyang City significantly decreased, resulting in a large volume of carbon emissions into the atmosphere, which will have an adverse influence on the atmospheric environment and sustainable development of agriculture.

Author Contributions

Conceptualization, B.Z. and Y.L.; methodology, B.Z. and Y.L.; software, B.Z.; validation, B.Z., Y.W. and Y.L.; formal analysis, B.Z.; investigation, B.Z., Y.W. and Y.L.; resources, Y.W.; data curation, B.Z.; writing—original draft preparation, B.Z.; writing—review and editing, B.Z. and G.Z.; visualization, B.Z.; supervision, Y.L.; project administration, Y.L.; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Social and economic data are from the official website of Liyang municipal government.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Geographical position of Liyang City.

Figure 2. Status quo of aquaculture ponds and ecological environment. (a) Aquaculture ponds and (b) contaminated water.

Figure 3. Spatial expansion and distribution of aquaculture ponds in Liyang City.

Table 1

Value equivalence of ecosystem services.

Grade 1 Ecosystem Service	Grade 2 Ecosystem Service	Cultivated Lands	Aquaculture Ponds
Supply service	Food production	1.36	0.8
	Raw materials production	0.09	0.23
Regulation service	Gas regulation	1.11	0.77
	Climate regulation	0.57	2.29
	Environmental purification	0.17	5.55
	Hydrological regulation	2.72	102.24
Support service	Soil conservation	0.01	0.93
	Maintain nutrient cycling	0.19	0.07
	Biological diversity	0.21	2.55
Cultural service	Aesthetic landscape	0.09	1.89

Table 2

Change in aquaculture ponds in each town of Liyang City from 2009 to 2019.

Towns	Transformation fromAquaculture Ponds into Cultivated Lands/ha(Area Decrease of Ponds)	Transformation fromCultivated Lands intoAquaculture Lands/ha(Area Increase of Ponds)	Actual Change in Aquaculture Ponds/ha
Licheng Town	71.75	124.32	52.56
Shanghuang Town	114.52	617.89	503.37
Tianmu Lake Town	61.01	136.30	75.29
Shangxing Town	363.60	929.37	565.78
Shezhu Town	361.84	2840.83	2478.99
Daitou Town	378.53	186.40	−192.13
Daibu Town	173.36	104.37	−68.99
Bieqiao Town	1235.68	345.68	−890.00
Zhuze Town	325.55	267.95	−57.60
Nandu Town	995.29	621.41	−373.88
Kunlun Street	124.83	104.68	−20.15
Total	4205.97	6279.21	2073.24

Note: In China’s administrative divisions, the level of streets and towns is the same.

Table 3

ESV of different service functions in Liyang City.

Classification ofEcosystem Service	2009		2019		ESV Change/RMB 100 Million Yuan
	ESV/RMB 100Million Yuan	Proportion	ESV/RMB 100Million Yuan	Proportion
Supply service	4.96	3.59%	4.05	2.89%	−0.92
Regulation service	126	91.39%	129	92.19%	2.93
Support service	4.75	3.44%	4.67	3.33%	−0.083
Cultural service	2.19	1.58%	2.23	1.59%	0.042
Total	138.21	100.00%	140.19	100.00%	1.98

Table 4

ESV of cultivated lands and aquaculture ponds in Liyang City.

Grade 1	Grade 2	2009	2019	2009	2019
		ESV of Cultivated Lands/RMB100 Million	ESV of Cultivated Lands/RMB100 Million	ESV of Aquaculture Ponds/RMB100 Million	ESV of Aquaculture Ponds/RMB100 Million
Supply service	Food production	3.66	2.75	0.82	0.89
	Raw materials production	0.24	0.18	0.24	0.25
Regulation service	Gas regulation	2.99	2.24	0.79	0.83
	Climate regulation	1.53	1.15	2.36	2.48
	Environmental cleaning	0.46	0.34	5.71	6.01
	Hydrological regulation	7.32	5.5	105	111.07
Support service	Soil conservation	0.03	0.02	0.96	1.01
	Maintenance of nutrient circulation	0.51	0.38	0.07	0.08
	Biological diversity	0.56	0.43	2.62	2.76
Cultural service	Aesthetic landscape	0.24	0.18	1.94	2.05
Total		17.54	13.17	120.51	127.42

References

1. Song, W.; Pijanowski, B.C. The effects of China’s cultivated land balance program on potential land productivity at a national scale. Appl. Geogr.; 2014; 46, pp. 158-170. [DOI: https://dx.doi.org/10.1016/j.apgeog.2013.11.009]

2. Song, X.; Zhu, O. Connotation of Multifunctional Cultivated Land and Its Implications for Cultivated Land Protection. Prog. Geogr.; 2012; 31, pp. 859-868.

3. Sen, W.; Daochi, Q.; Qingyuan, Y.; Chunxia, W.; Chunyang, L. Study on the Index System of Infield Resources Safety Estimation. Chin. Agric. Sci. Bull.; 2007; 08, pp. 466-470.

4. Owen, O.S.; Chiras, D.D. Natural Resource Conservation: An Ecological Approach; Macmillan: New York, NY, USA, 1980; 593.

5. Li, Y.; Zhao, B.; Huang, A.; Xiong, B.; Song, C. Characteristics and driving forces of non-grain production of cultivated land from the perspective of food security. Sustainability; 2021; 13, 14047. [DOI: https://dx.doi.org/10.3390/su132414047]

6. Su, Y.; Qian, K.; Lin, L.; Wang, K.; Guan, T.; Gan, M. Identifying the driving forces of non-grain production expansion in rural China and its implications for policies on cultivated land protection. J. Land Use Policy; 2020; 92, 104435. [DOI: https://dx.doi.org/10.1016/j.landusepol.2019.104435]

7. Cheng, L.; Brown, G.; Liu, Y.; Searle, G. An evaluation of contemporary China’s land use policy–The Link Policy: A case study from Ezhou, Hubei Province. Land Use Policy; 2020; 91, 104423. [DOI: https://dx.doi.org/10.1016/j.landusepol.2019.104423]

8. Pan, P.P.; Yang, G.S.; Wang, X.M.; Wang, L.Y.; Wang, X.R.; Jing-Feng, G.E. Analysis and Forecast on Dynamic Response of Grain Production to Cropland Use Changes in Taihu Lake Basin. Resour. Environ. Yangtze Basin; 2019; 28, pp. 2364-2375.

9. Berlanga-Robles, C.A.; Ruiz-Luna, A.; Bocco, G.; Vekerdy, Z. Spatial analysis of the impact of shrimp culture on the coastal wetlands on the Northern coast of Sinaloa, Mexico. Ocean Coast. Manag.; 2011; 54, pp. 535-543. [DOI: https://dx.doi.org/10.1016/j.ocecoaman.2011.04.004]

10. Cao, L.; Wang, W.; Yang, Y.; Yang, C.; Yuan, Z.; Xiong, S.; Diana, J. Environmental impact of aquaculture and countermeasures to aquaculture pollution in China. Environ. Sci. Pollut. Res.-Int.; 2007; 14, pp. 452-462.

11. Costanzo, S.D.; O’Donohue, M.J.; Dennison, W.C. Assessing the influence and distribution of shrimp pond effluent in a tidal mangrove creek in north-east Australia. Mar. Pollut. Bull.; 2004; 48, pp. 514-525. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2003.09.006]

12. Yang, P.; Lai, D.Y.; Jin, B.; Bastviken, D.; Tan, L.; Tong, C. Dynamics of dissolved nutrients in the aquaculture shrimp ponds of the Min River estuary, China: Concentrations, fluxes and environmental loads. Sci. Total Environ.; 2017; 603, pp. 256-267. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2017.06.074] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28628817]

13. Morshed, M.M.; Islam, M.S.; Lohano, H.D.; Shyamsundar, P. Production externalities of shrimp aquaculture on paddy farming in coastal Bangladesh. Agric. Water Manag.; 2020; 238, 106213. [DOI: https://dx.doi.org/10.1016/j.agwat.2020.106213]

14. Páez-Osuna, F.; Álvarez-Borrego, S.; Ruiz-Fernández, A.C.; García-Hernández, J.; Jara-Marini, M.E.; Bergés-Tiznado, M.E.; Piñón-Gimate, A.; Alonso-Rodríguez, R.; Soto-Jiménez, M.F.; Frías-Espericueta, M.G. Environmental status of the Gulf of California: A pollution review. Earth-Sci. Rev.; 2017; 166, pp. 181-205. [DOI: https://dx.doi.org/10.1016/j.earscirev.2017.01.014]

15. Pan, G.X.; Zhao, Q.G. Research on evolution of soil carbon pool for farmland in China: Global change and national food security. J. Adv. Earth Sci.; 2005; 20, pp. 384-393.

16. Wang, Y.; Chen, Y.J.; Yi, X.Y.; Xiao, B.L. Problem of “non-food crop” in circulation of cultivated lands and countermeasures. J. China Agric. Resour. Reg. Plan.; 2011; 32, pp. 13-16.

17. Ren, C.; Wang, Z.; Zhang, B.; Li, L.; Chen, L.; Song, K.; Jia, M. Remote monitoring of expansion of aquaculture ponds along coastal region of the Yellow River Delta from 1983 to 2015. J. Chin. Geogr. Sci.; 2018; 28, pp. 430-442. [DOI: https://dx.doi.org/10.1007/s11769-017-0926-2]

18. Ren, C.; Wang, Z.; Zhang, Y.; Zhang, B.; Chen, L.; Xi, Y.; Xiao, X.; Doughty, R.B.; Liu, M.; Jia, M. et al. Rapid expansion of coastal aquaculture ponds in China from Landsat observations during 1984–2016. Int. J. Appl. Earth Obs. Geoinf.; 2019; 82, 101902. [DOI: https://dx.doi.org/10.1016/j.jag.2019.101902]

19. Jia, M.; Liu, M.; Wang, Z.; Mao, D.; Ren, C.; Cui, H. Evaluating the effectiveness of conservation on mangroves: A remote sensing-based comparison for two adjacent protected areas in Shenzhen and Hong Kong, China. Remote Sens.; 2016; 8, 627. [DOI: https://dx.doi.org/10.3390/rs8080627]

20. Virdis, S.G.P. An object-based image analysis approach for aquaculture ponds precise mapping and monitoring: A case study of Tam Giang-Cau Hai Lagoon, Vietnam. Environ. Monit. Assess.; 2014; 186, pp. 117-133. [DOI: https://dx.doi.org/10.1007/s10661-013-3360-7]

21. Zhang, T.; Li, Q.; Yang, X.; Zhou, C.; Su, F. Automatic mapping aquaculture in coastal zone from TM imagery with OBIA approach. Proceedings of the 2010 18th International Conference on Geoinformatics; Beijing, China, 18–20 June 2010; pp. 1-4.

22. Eng, C.T.; Paw, J.N.; Guarin, F.Y. The environmental impact of aquaculture and the effects of pollution on coastal aquaculture development in Southeast Asia. Mar. Pollut. Bull.; 1989; 20, pp. 335-343.

23. Duan, Y.; Li, X.; Zhang, L.; Chen, D.; Ji, H. Mapping national-scale aquaculture ponds based on the Google Earth Engine in the Chinese coastal zone. Aquaculture; 2020; 520, 734666. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2019.734666]

24. Prasad, K.A.; Ottinger, M.; Wei, C.; Leinenkugel, P. Assessment of coastal aquaculture for India from Sentinel-1 SAR time series. Remote Sens.; 2019; 11, 357. [DOI: https://dx.doi.org/10.3390/rs11030357]

25. Azad, A.; Jensen, K.R.; Lin, C.K. Coastal aquaculture development in Bangladesh: Unsustainable and sustainable experiences. Environ. Manag.; 2009; 44, pp. 800-809. [DOI: https://dx.doi.org/10.1007/s00267-009-9356-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19701597]

26. Ottinger, M.; Clauss, K.; Kuenzer, C. Large-scale assessment of coastal aquaculture ponds with Sentinel-1 time series data. Remote Sens.; 2017; 9, 440. [DOI: https://dx.doi.org/10.3390/rs9050440]

27. Xia, Z.; Guo, X.; Chen, R. Automatic extraction of aquaculture ponds based on Google Earth Engine. Ocean. Coast. Manag.; 2020; 198, 105348. [DOI: https://dx.doi.org/10.1016/j.ocecoaman.2020.105348]

28. Zou, B.L. The Phenomenon of Dig Ponds to Breed Fish and Cultivated Land Protection—Based on the Investigation in Yangchun City. South China Rural Area; 2013; 29, pp. 24-27.

29. Rajitha, K.; Mukherjee, C.; Chandran, R.V. Applications of remote sensing and GIS for sustainable management of shrimp culture in India. Aquacult. Eng.; 2007; 36, pp. 1-17. [DOI: https://dx.doi.org/10.1016/j.aquaeng.2006.05.003]

30. Wu, H.; Peng, R.; Yang, Y.; He, L.; Wang, W.; Zheng, T.; Lin, G. Mariculture pond influence on mangrove areas in south China: Significantly larger nitrogen and phosphorus loadings from sediment wash-out than from tidal water exchange. Aquaculture; 2014; 426, pp. 204-212. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2014.02.009]

31. Ottinger, M.; Clauss, K.; Kuenzer, C. Aquaculture: Relevance, distribution, impacts and spatial assessments—A review. Ocean Coast. Manag.; 2016; 119, pp. 244-266. [DOI: https://dx.doi.org/10.1016/j.ocecoaman.2015.10.015]

32. Ren, P.; Tao, W.U.; Zhou, J.M. Study on Value Loss Evaluation Model and Compensation Mechanism of Cultivated Land Conversion. Sci. Agric. Sin.; 2014; 47, pp. 786-795.

33. Gómez-Baggethun, E.; De Groot, R.; Lomas, P.L.; Montes, C. The history of ecosystem services in economic theory and practice: From early notions to markets and payment schemes. Ecol. Econ.; 2010; 69, pp. 1209-1218. [DOI: https://dx.doi.org/10.1016/j.ecolecon.2009.11.007]

34. De Groot, R.S.; Wilson, M.A.; Boumans, R.M. A typology for the classification, description and valuation of ecosystem functions, goods and services. Ecol. Econ.; 2002; 41, pp. 393-408. [DOI: https://dx.doi.org/10.1016/S0921-8009(02)00089-7]

35. Xie, G.D.; Lu, C.X.; Cheng, S.K. The progress of research on value of global ecosystem service. Resour. Sci.; 2001; 23, pp. 5-9.

36. Program, M.E.A. Ecosystems and Human Well-Being: Wetlands and Water Synthesis; World Resources Institute: Washington, DC, USA, 2005.

37. Hasan, S.S.; Zhen, L.; Miah, M.G.; Ahamed, T.; Samie, A. Impact of land use change on ecosystem services: A review. Environ. Dev.; 2020; 34, 100527. [DOI: https://dx.doi.org/10.1016/j.envdev.2020.100527]

38. Costanza, R.; d’Arge, R.; De Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’neill, R.V.; Paruelo, J. The value of the world’s ecosystem services and natural capital. Nature; 1997; 387, pp. 253-260. [DOI: https://dx.doi.org/10.1038/387253a0]

39. Reid, W.V.; Mooney, H.A.; Cropper, A.; Capistrano, D.; Carpenter, S.R.; Chopra, K.; Dasgupta, P.; Dietz, T.; Duraiappah, A.K.; Hassan, R. Ecosystems and Human Well-Being-Synthesis: A Report of the Millennium Ecosystem Assessment; Island Press: Washington, DC, USA, 2005.

40. Wang, Z.; Cao, J.; Zhu, C.; Yang, H. The impact of land use change on ecosystem service value in the upstream of Xiong’an new area. Sustainability; 2020; 12, 5707. [DOI: https://dx.doi.org/10.3390/su12145707]

41. Peng, J.; Tian, L.; Liu, Y.; Zhao, M.; Wu, J. Ecosystem services response to urbanization in metropolitan areas: Thresholds identification. Sci. Total Environ.; 2017; 607, pp. 706-714. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2017.06.218]

42. Su, S.; Xiao, R.; Jiang, Z.; Zhang, Y. Characterizing landscape pattern and ecosystem service value changes for urbanization impacts at an eco-regional scale. Appl. Geogr.; 2012; 34, pp. 295-305. [DOI: https://dx.doi.org/10.1016/j.apgeog.2011.12.001]

43. Matthias, S.; Koellner, T.; Alkemade, R.; Arnhold, S.; Bagstad, K.J.; Erb, K.H.; Frank, K.; Kastner, T.; Kissinger, M.; Liu, J. Interregional flows of ecosystem services: Concepts, typology and four cases. Ecosyst. Serv.; 2018; 31, pp. 231-241.

44. Burkhard, B.; Kroll, F.; Nedkov, S.; Müller, F. Mapping ecosystem service supply, demand and budgets. Ecol. Indic.; 2012; 21, pp. 17-29. [DOI: https://dx.doi.org/10.1016/j.ecolind.2011.06.019]

45. Wen, Y.; Ma, L.; Xie, J.; Ma, Y.; Zhu, Y.; Liu, G. Quantitative Research of Ecosystem Service Function Space Transfer-A Case of Guanting Reservoir Watershed Region. Environ. Prot. Sci; 2018; 44, pp. 95-102.

46. Shi, Y.Y.; Xiao, L.; Huang, X.J.; Miao, Y.U. Arable Land Use Transitions and Its Response of Ecosystem Services Value Change in Jiangsu Coastal Areas. J. Nat. Resour.; 2017; 32, pp. 961-976.

47. Xie, G.D.; Zang, C.X.; Zang, L.M.; Chen, W.H.; Li, S.M. Improvement of value method for ecosystem service based on unit area value equivalent factor. J. Nat. Resour.; 2015; 30, pp. 1243-1254.

48. Lei, J.C.; Liu, J.X.; Yong, F.; Liu, H.M.; Jun, W.U.; Ding, H.; Wang, J.M.; Song-Qin, W.U.; Cheng, S.; Cui, P. The influence of land use variation on ecosystem service value in Xunwu County. ACTA Ecol. Sin.; 2019; 39, pp. 3089-3099.

49. Ya-Ru, L.; Cong, W.; Li-Jiao, Y. Impacts of land use change on ecosystem services in the agricultural area of North China Plain: A case study of Shangqiu City, Henan Province, China. Chin. J. Appl. Ecol.; 2018; 29, pp. 1597-1606.

50. Xiaowei, C.; Xianjin, H.; Zeqing, Z.; Mei, Z.; Qilin, L.; Li, L.; Junyu, L.U. The influence of land use variation in Jiangsu Province on carbon capacity of terrestrial ecosystem. Resour. Sci.; 2011; 33, pp. 1932-1939.

51. Chuai, X.; Huang, X.; Wang, W.; Zhao, R.; Zhang, M.; Wu, C. Land use, total carbon emissions change and low carbon land management in Coastal Jiangsu, China. J. Clean. Prod.; 2015; 103, pp. 77-86. [DOI: https://dx.doi.org/10.1016/j.jclepro.2014.03.046]

52. Yang, S.S.; Zou, C.X.; Shen, W.S.; Shen, R.P.; Xu, D.L. Construction of ecological security patterns based on ecological red line: A case study of Jiangxi Province. Chin. J. Ecol.; 2016; 35, pp. 250-258.

53. Quéré, C.; Raupach, M.R.; Canadell, J.G.; Marland, G.; Bopp, L.; Ciais, P.; Conway, T.J.; Doney, S.C.; Feely, R.A.; Foster, P. Trends in the sources and sinks of carbon dioxide. Nat. Geosci.; 2009; 2, pp. 831-836. [DOI: https://dx.doi.org/10.1038/ngeo689]

54. Post, W.M.; Kwon, K.C. Soil carbon sequestration and land-use change: Processes and potential. Glob. Chang. Biol.; 2010; 6, pp. 317-327. [DOI: https://dx.doi.org/10.1046/j.1365-2486.2000.00308.x]

55. Duan, H.; Zhang, H.; Huang, Q.; Zhang, Y.; Hu, M.; Niu, Y.; Zhu, J. Characterization and environmental impact analysis of sea land reclamation activities in China. Ocean Coast. Manag.; 2016; 130, pp. 128-137. [DOI: https://dx.doi.org/10.1016/j.ocecoaman.2016.06.006]

56. Gellrich, M.; Baur, P.; Koch, B.; Zimmermann, N.E. Agricultural land abandonment and natural forest re-growth in the Swiss mountains: A spatially explicit economic analysis. Agric. Ecosyst. Environ.; 2007; 118, pp. 93-108. [DOI: https://dx.doi.org/10.1016/j.agee.2006.05.001]

57. Wu, Z.; Lin, C.; Shao, H.; Feng, X.; Wang, S. Ecological risk assessment and difference analysis of pit ponds under different ecological service functions-A case study of Jianghuai ecological Economic Zone. Ecol. Indic.; 2021; 129, 107860. [DOI: https://dx.doi.org/10.1016/j.ecolind.2021.107860]

58. Wei, S. Land-use/land-cover change and ecosystem service provision in China. Sci. Total Environ.; 2017; 576, pp. 705-719.

59. Pueppke, S.G.; Nurtazin, S.; Ou, W. Water and land as shared resources for agriculture and aquaculture: Insights from Asia. Water; 2020; 12, 2787. [DOI: https://dx.doi.org/10.3390/w12102787]

60. Flaherty, M.; Vandergeest, P.; Miller, P. Rice paddy or shrimp pond: Tough decisions in rural Thailand. World Devel.; 1999; 27, pp. 2045-2060. [DOI: https://dx.doi.org/10.1016/S0305-750X(99)00100-X]

61. Chowdhury, M.A.; Shivakoti, G.P.; Salaquzzaman, M.A. conceptual framework for the sustainability assessment procedures of the shrimp aquaculture industry in coastal Bangladesh. Int. J. Agric. Resour. Gov. Ecol.; 2006; 5, pp. 162-184. [DOI: https://dx.doi.org/10.1504/IJARGE.2006.009162]

62. Li, Y.Y.; Tan, M.H.; Zang, H.Q. The influence of variation in spatial pattern of China’s cultivated lands on food production. Strateg. Study CAE; 2018; 20, pp. 90-95.

63. Di, X.G.; Xia, L.C.; Kui, C.S.; Du, Z. Value evaluation for ecological assets of Qinghai-Tibet Plateau. J. Nat. Resour.; 2003; 18, pp. 189-196.