Identifying the Relationship between Soil Properties and Rice Growth for Improving Consolidated Land in the Yangtze River Delta, China

: China has widely implemented land consolidation, which was expected to increase the amount of cultivated land and enhance grain yields. Key components of land consolidation include ﬁlling mall waterbodies and leveling land, both of which have strong impacts on the environment in the Yangtze River Delta. The impacts of land consolidation on soil ecology and agricultural production are not yet clear. Here, we conducted a ﬁeld survey of soil properties and rice growth to detect the effects of land consolidation in the ﬁrst growing season. The normalized difference vegetation index (NDVI) was used to analyze the remote sensing data. We found signiﬁcant differences in the soil properties under different types of land leveling, with a general NDVI pattern of: control > borrowed topsoil area > ﬁlled waterbodies area > topsoil cutting area. We found signiﬁcant heterogeneity in rice NDVI after land consolidation. The NDVI of rice had extremely signiﬁcant positive correlations with soil organic matter and available zinc. The spatial variation in soil properties caused by land consolidation was a dominant factor leading to the heterogeneity of rice NDVI. Fertilizing soil and strengthening ﬁeld management should be adopted to provide more ecological services while increasing quantity.


Introduction
Land consolidation is a valuable tool for land-use management in many countries, such as for example China and Iran in Asia [1,2], Poland and Spain in Europe [3,4], and others [5,6]. China's land consolidation programs started in 2000, and were initially intended to increase cultivated land area and grain-yield capacity [7,8]. Chinese has invested about 300,000 million CNY in land consolidation in the past 10 years [1], accumulating 30 million hectares of consolidated farmland, and making substantial contributions to farmland protection and grain security [1,9,10]. As a result of a growing population and limited farmland in China, land consolidation has overemphasized the increase of farmland quantity and grain yield capacity for a long time [8,11], resulting in managers paying little attention to the negative effects of land consolidation. As the economy improves in China, more attention has turned to ecological security and environmental protection [1,12,13]. Land consolidation may now emphasize not only farmland quantity and grain yield, but also ecosystem services [12,13]. Farmland not only produces grain, it also serves as a multifunctional provider of ecological services, densely covered, which is typical of Yangtze River Delta. This study area is raining season in summer, with two crops each year, which is mainly a rice-wheat rotation. Between January-May 2017, land leveling and a water conservancy project were completed, and then farmed by big plantation owners, with unified management.
Sustainability 2018, 10, x FOR PEER REVIEW 3 of 14 densely covered, which is typical of Yangtze River Delta. This study area is raining season in summer, with two crops each year, which is mainly a rice-wheat rotation. Between January-May 2017, land leveling and a water conservancy project were completed, and then farmed by big plantation owners, with unified management.

Soil Sampling and Analysis
Soil sampling was collected on 1 September 2017, about four months after the land was leveled. The boundaries of each treatment area (filled waterbodies (FW), borrowed topsoil (BT), topsoilcutting (TC), and a control area (CA)) were determined using engineering maps, remote sensing images, and field investigations. We collected six samples of FW, 15 samples of TC, 16 samples of BT, and six samples of CA ( Figure 1). Soil samples-five pooled samples selected randomly from each plot [42] of (0-20 cm depth)-were collected for each of the four soil treatment categories. Coordinates were taken with a handheld Global Positioning System (GPS) (Qmini A5/A7, Hi-target, Shenzhen, China).
Soil samples were dried at room temperature for about one week, and then sieved through a 0.25-mm sieve for physicochemical analysis [42]. We chose to test soil properties that are known to affect rice growth, such as soil texture, pH, soil organic matter (SOM) content, total N (TN) [43], available P (AP), available K (AK) [44], available Zn (AZn), and available Si (ASi). Soil texture was measured using a particle size analyzer (BT-9300Z, Baxter, Dandong, China). Soil pH was measured by the potentiometric method (water: soil ratio of 1:2.5) [42]. SOM was determined using the potassium dichromate oxidation outer heating method [45]. TN was determined using the semimicro Kjeldahl method [46]. Soil AP was determined using the hydrochloric acid ammonium chloride method [47]. Soil AK was determined by caries by ammonium acetate extraction-flame photometry [46] (FP640, Jingke, Shanghai, China). Soil AZn was determined using diethylene triamine pentacetate acid (DTPA) extraction atomic absorption spectrophotometry [48] (TAS-990, PGENERAL, Beijing, China). Soil ASi was determined using the citric acid buffered leaching silicon molybdenum blue colorimetric method [46].

Interpretation and Processing of Remote Sensing Data
Three GeoEye images (Earth Eye Satellite, Inc., Dulles, VA, USA) were purchased for 31 December 2016, 6 April 2017, and 23 August 2017, respectively representing before, during, and after land consolidation ( Figure 2). The spatial resolution is 1.65 m [49]. We trimmed the three-phase images, and processed them with the conventional methods (geometric correction, atmospheric

Soil Sampling and Analysis
Soil sampling was collected on 1 September 2017, about four months after the land was leveled. The boundaries of each treatment area (filled waterbodies (FW), borrowed topsoil (BT), topsoil-cutting (TC), and a control area (CA)) were determined using engineering maps, remote sensing images, and field investigations. We collected six samples of FW, 15 samples of TC, 16 samples of BT, and six samples of CA ( Figure 1). Soil samples-five pooled samples selected randomly from each plot [42] of (0-20 cm depth)-were collected for each of the four soil treatment categories. Coordinates were taken with a handheld Global Positioning System (GPS) (Qmini A5/A7, Hi-target, Shenzhen, China).
Soil samples were dried at room temperature for about one week, and then sieved through a 0.25-mm sieve for physicochemical analysis [42]. We chose to test soil properties that are known to affect rice growth, such as soil texture, pH, soil organic matter (SOM) content, total N (TN) [43], available P (AP), available K (AK) [44], available Zn (AZn), and available Si (ASi). Soil texture was measured using a particle size analyzer (BT-9300Z, Baxter, Dandong, China). Soil pH was measured by the potentiometric method (water:soil ratio of 1:2.5) [42]. SOM was determined using the potassium dichromate oxidation outer heating method [45]. TN was determined using the semi-micro Kjeldahl method [46]. Soil AP was determined using the hydrochloric acid ammonium chloride method [47]. Soil AK was determined by caries by ammonium acetate extraction-flame photometry [46] (FP640, Jingke, Shanghai, China). Soil AZn was determined using diethylene triamine pentacetate acid (DTPA) extraction atomic absorption spectrophotometry [48] (TAS-990, PGENERAL, Beijing, China). Soil ASi was determined using the citric acid buffered leaching silicon molybdenum blue colorimetric method [46].

Interpretation and Processing of Remote Sensing Data
Three GeoEye images (Earth Eye Satellite, Inc., Dulles, VA, USA) were purchased for 31 December 2016, 6 April 2017, and 23 August 2017, respectively representing before, during, and after land consolidation ( Figure 2). The spatial resolution is 1.65 m [49]. We trimmed the three-phase images, and processed them with the conventional methods (geometric correction, atmospheric correction, radiation correction, and other pretreatment) [49,50]. Remote sensing is widely used in agricultural monitoring [43,[51][52][53][54][55][56][57][58][59]. The selected measure, NDVI, is sensitive to surface cover, and thus reflects the growth status and coverage of rice [52][53][54][59][60][61]. We used ENVI 5.0 software to process the GeoEye images, and combined engineering drawings of the field areas to extract the NDVI value of each field. NDVI is the difference between the values of the near-infrared and visible-infrared bands divided by the sum of the two bands and its ratio range is [−1, 1] [61]. The formula is as follows: Spatial clustering analysis can identify the characteristics of rice growth. Spatial clustering analysis was based on NDVI values [60]. We calculated the Local Moran's I index associated with the NDVI values of each pixel [62][63][64][65][66] using the following formula: In Equation (2), X i and X j are the NDVI values of the i and j pixels, respectively. X is the average of all of the fields' NDVI values, n is the total number of pixels, and W ij is the spatial weight between the i and j of the pixel. An I i > 0 indicates that the picture element has neighboring pixels containing NDVI values that are equally high or equally low, and thus, the picture element belongs to a high-value or low-value cluster. If I i is <0, it indicates that the picture element has neighboring pixels containing different NDVI values, and thus the picture element is an abnormal value. By setting the confidence interval to 95%, all of the pixels can be distinguished from statistically significant high-high (HH) clusters, low-low (LL) clusters, which are outliers whose high values are mainly surrounded by low-high (HL) outliers, and high-low (LH) outliers whose low values are mainly surrounded by high values. As one example, if a treatment area had a low value in the LL cluster of a high-low outlier (LH) region that was surrounded by high values that were significantly lower than the values of the other treatment types, it would indicate that this treatment area was less conducive to crop growth. correction, radiation correction, and other pretreatment) [49,50]. Remote sensing is widely used in agricultural monitoring [43,[51][52][53][54][55][56][57][58][59]. The selected measure, NDVI, is sensitive to surface cover, and thus reflects the growth status and coverage of rice [52][53][54][59][60][61]. We used ENVI 5.0 software to process the GeoEye images, and combined engineering drawings of the field areas to extract the NDVI value of each field. NDVI is the difference between the values of the near-infrared and visibleinfrared bands divided by the sum of the two bands and its ratio range is [−1, 1] [61]. The formula is as follows: Spatial clustering analysis can identify the characteristics of rice growth. Spatial clustering analysis was based on NDVI values [60]. We calculated the Local Moran's I index associated with the NDVI values of each pixel [62][63][64][65][66] using the following formula: In Equation (2), and are the NDVI values of the and pixels, respectively. ̅ is the average of all of the fields' NDVI values, is the total number of pixels, and is the spatial weight between the and of the pixel. An > 0 indicates that the picture element has neighboring pixels containing NDVI values that are equally high or equally low, and thus, the picture element belongs to a high-value or low-value cluster. If is <0, it indicates that the picture element has neighboring pixels containing different NDVI values, and thus the picture element is an abnormal value. By setting the confidence interval to 95%, all of the pixels can be distinguished from statistically significant high-high (HH) clusters, low-low (LL) clusters, which are outliers whose high values are mainly surrounded by low-high (HL) outliers, and high-low (LH) outliers whose low values are mainly surrounded by high values. As one example, if a treatment area had a low value in the LL cluster of a high-low outlier (LH) region that was surrounded by high values that were significantly lower than the values of the other treatment types, it would indicate that this treatment area was less conducive to crop growth. The NDVI value of rice was affected by many factors. Stepwise regression analysis can continuously enter each factor into the regression equation and screen out the key factors. The mathematical model of stepwise regression is: In Equation (3), " " is the unknown parameter in the model, and " " is the random error. The NDVI value of rice was affected by many factors. Stepwise regression analysis can continuously enter each factor into the regression equation and screen out the key factors. The mathematical model of stepwise regression is: y = a 0 + a 1 x 1 + a 2 x 2 + · · · + a n x n + B In Equation (3), "a" is the unknown parameter in the model, and "B" is the random error. Statistical analyses of experimental data were performed using SPSS 20.0 (IBM, New York, NY, USA) for the Tukey-Kramer method on the NDVI values and soil properties to determine whether different land-consolidation treatments differed significantly from the control. Correlation analysis was used to determine the relationship between NDVI values and each property in the soil, and line regression analysis was used to determine the relationship between NDVI values and SOM, AZn. Mapping was done using ArcGIS 10.2 (ESRI, San Diego, CA, USA) and Origin 9.0 (Origin Lab, Northampton, MA, USA) software.

Short-Term Interference of Land Consolidation on Soil Properties
In addition to AP, the soil texture, pH, SOM, TN, AK, AZn, and ASi in land-leveling treatments differed significantly from CA (Table 1). Soil clay and silt content in land-consolidation treatments was higher than in CA, and soil viscosity was also higher in the land-consolidation treatments. Land leveling had little effect on the soil surface pH in TC and BT, but in FW, the pH was significantly lower (p < 0.05). In the land-leveling treatments, the content of SOM, TN, AK, and AZn were generally lower, with the following trend by treatment type (highest to lowest): CA > BT > FW > TC. The difference in ASi content was opposite; the ASi content in TC was 327.35 mg kg −1 , which was significantly higher than all of the other treatments. Land consolidation had little effect on AP, with no significant change. In summary, land leveling had a significant influence on soil physicochemical properties, and the impact of BT was smaller than that of FW or TC.

Spatial Difference of Rice Growth in Different Land Leveling Areas
There was a significant difference in the growth potential of rice between CA and the two types of FW and TC (p < 0.05) ( Table 2). Among the treatments, the average value of NDVI in CA was the highest, and the average NDVI value in TC was 10.67% smaller. The growth of rice in CA was the highest and uniform, and it was consistent with the results of the field survey (Figure 3), suggesting that all of the land-consolidation treatments caused changes in soil properties that negatively impacted rice growth.
The spatial clustering of NDVI values in the study area is shown in Figure 4. In general, there were five types of clusters: not significant, LL cluster, HH cluster, LH outlier, and HL outlier. Most of the NDVI values fell in the not-significant area. There were also large areas of LL clusters in the study area, mainly in TC and FW. The HH cluster areas were mostly in BT with better rice growth, which was consistent with the results of the field observations. study area, mainly in TC and FW. The HH cluster areas were mostly in BT with better rice growth, which was consistent with the results of the field observations. Within each row, the difference is not significant when there is a letter with the same mark, and the difference is significant when there are different mark letters.

Relationship between Soil Properties and Rice Growth
After land consolidation, farmland was planted by big plantation owners with unified management. The difference in the rice growth was mainly attributed to the change of soil nutrients caused by land leveling. The NDVI values of all of the plots had a significant positive correlation with SOM, TN, AK, and AZn (Table 3). There was a significant negative correlation between ASi content and NDVI values. No significant correlation was found between NDVI values and soil clay, pH, or AP (p > 0.05). In this study, SOM, AZn, TN, and AK, which were significantly correlated with NDVI, were selected. Stepwise regression analysis was used to establish a stepwise regression equation. At p < 0.05, only SOM and AZn entered the equation, and the organic matter entered first. This indicates that SOM and AZn were the keys factors affecting rice growth.

Relationship between Soil Properties and Rice Growth
After land consolidation, farmland was planted by big plantation owners with unified management. The difference in the rice growth was mainly attributed to the change of soil nutrients caused by land leveling. The NDVI values of all of the plots had a significant positive correlation with SOM, TN, AK, and AZn (Table 3). There was a significant negative correlation between ASi content and NDVI values. No significant correlation was found between NDVI values and soil clay, pH, or AP (p > 0.05). In this study, SOM, AZn, TN, and AK, which were significantly correlated with NDVI, were selected. Stepwise regression analysis was used to establish a stepwise regression equation. At p < 0.05, only SOM and AZn entered the equation, and the organic matter entered first. This indicates that SOM and AZn were the keys factors affecting rice growth. There was a significant positive correlation between soil properties and the NDVI values of rice growth in different leveling types, with R 2 being the largest in CA ( Figure 5). In terms of individual indicators, the correlations among SOM, AZn, and NDVI were significant ( Figure 5), but the correlations among SOM, AZn, and NDVI were significantly different in different leveling types. For example, the correlation coefficients between SOM and NDVI were 0.465 and 0.796 in FW and CA. Although this aspect showed that the disturbances had caused the heterogeneous growth of rice after land consolidation, SOM positively affected rice growth. Our results also showed that the soil physicochemical properties in FW showed an overall decrease trend; soil properties and soil stability had been strongly disturbed. As a result, the correlation between soil properties and NDVI was reduced. AZn was also an important factor positively affecting rice growth, whereas there was a significant negative correlation between ASi and NDVI (p < 0.05). Rice depends on silicon to grow, and silicon can promote photosynthesis in rice growth. Due to the high content of ASi in the study area, which exceeds the suitable concentration range for rice growth, the correlation coefficient was only 0.238. From the field investigations, the NDVI value not only reflected rice growth, it also reflected the distribution of soil nutrients after land leveling. In summary, NDVI may be a convenient replacement for field investigations to monitor soil quality after large-scale land consolidation. physicochemical properties in FW showed an overall decrease trend; soil properties and soil stability had been strongly disturbed. As a result, the correlation between soil properties and NDVI was reduced. AZn was also an important factor positively affecting rice growth, whereas there was a significant negative correlation between ASi and NDVI (p < 0.05). Rice depends on silicon to grow, and silicon can promote photosynthesis in rice growth. Due to the high content of ASi in the study area, which exceeds the suitable concentration range for rice growth, the correlation coefficient was only 0.238. From the field investigations, the NDVI value not only reflected rice growth, it also reflected the distribution of soil nutrients after land leveling. In summary, NDVI may be a convenient replacement for field investigations to monitor soil quality after large-scale land consolidation.

Variability of Soil Properties under Different Land-Leveling Patterns
Land consolidation is an important tool for promoting sustainable use [10,24], but deep turning, excavation, and landfilling during the land consolidation process inevitably break the original soil configuration and change the vertical or horizontal distribution of soil nutrients [20,33,67]. In the next

Variability of Soil Properties under Different Land-Leveling Patterns
Land consolidation is an important tool for promoting sustainable use [10,24], but deep turning, excavation, and landfilling during the land consolidation process inevitably break the original soil configuration and change the vertical or horizontal distribution of soil nutrients [20,33,67]. In the next 10 years, China will carry out large-scale land consolidation. The state requires that 53.33 million hectares of high-standard permanent basic farmland be completed by 2020, and 80.0 million hectares of high-standard permanent basic farmland be completed by 2030. However, a slight change in the soil may cause major changes in the global environment [29,68,69]. In particular, large-scale land consolidation in China is bound to have an impact on the soil carbon and nitrogen cycle. When farmland is leveled, construction units generally adopt a program of soil balance within a small area to reduce the transportation costs. This program may increase the probability of borrowed topsoil and topsoil cutting.
As a complex project, land consolidation has caused large disturbances to soil properties. The variation properties of soil traits were different due to the different construction methods and different types land leveling. During land consolidation, a large number of deep soils are turned over to the surface, and water bodies that are filled are often filled with mud first, making it difficult for borrowed soil to be compacted. It is common for soils to experience subsidence after land consolidation. Land consolidation has caused the disturbances to soil physicochemical properties in FW, TC, and BT from Table 1. The pH content in FW was 5.75 to 5.98, which was significantly lower than the other three types (p < 0.05). We speculate that the main reason for the serious soil acidification was that the slurry sludge was derived from the silt of the river. The acidification of river bottom mud is a serious issue in the Yangtze River Delta, and high groundwater levels can easily cause acid regurgitation. This feature was very similar to the results of previous studies [70,71]. The thickness of the surface soil in FW and TC was generally thin, and the texture of the soil was heavy. This can easily form a barrier layer that negatively affects crop growth. The ASi mainly exists in the silicic acid solution of soil. After land consolidation, the topsoil is fully mixed with the underlying soil, and the viscosity of the new cultivation layer increased, so it can absorb more silicic acid. We suspect this absorption was the reason that FW, TC, and BT were all higher in ASi than in the CA.
Soil nutrients are influenced by both human and natural factors. It is generally considered that with the increase in land-leveling practices, SOM, N, P, and K will decrease at first, and then increase, so fertility will decline, but only in the short-term [28,31,72]. The field investigations in this study also found that land leveling had significant negative effects on soil properties and crop growth in the short term, which is consistent with previous research results [20,30,67]. This study found that the changes of AP were not significant, which differs from the results of Brey et al. [20]. Under conditions of consistent field management, AP are mainly controlled by the soil parent material [68]. Land leveling mainly disturbs the surface of the soil, and thus rarely affects the conditions of the parent materials or resulting AP. We found no significant changes in AP among our treatments. The significant decreases observed in previous studies were attributed the loss of AP to inconsistent fertilization and field management [2,67].

Factors Influencing of the Ndvi Value of Rice Growth
The surface layer of the soil contains abundant organic matter and nutrients. From the stepwise regression results in Table 3, SOM and AZn were the controlling factors for the NDVI value of rice growth. In this study, the NDVI value of CA was the largest, and the topsoil was considered a highly developed mature soil layer. In addition, with the same time of rain and heat in Yangtze River Delta, it is more suitable for rice planting and rice growth than most parts of China [73]. The NDVI values of rice in TC and BT decreased, and the variation among values became larger. There was clear spatial variation in the NDVI metrological characteristics of rice growth within the land consolidation project, and the NDVI value BT was significantly higher than that in TC (p < 0.05). The spatial cluster analysis showed that the LL cluster was mainly TC and FW, which was probably a result of different land-leveling types. When farmland is leveled, the method of cutting high and filling low is generally used to flatten the surface of the farmland. The surface matured soil of TC was borrowed to BT, and the surface soil in BT is fertile. However, the NDVI of the FWs was lower, and there was an LL cluster; this was significantly different from CA, because the FWs were filled with mud and then borrowed topsoil of TC, which cannot be compacted. After the summer rainfall, uneven settlement occurred in the loose soil layer, which disturbed the borrowed soil, and resulted in the variation of soil nutrients in the cultivated horizon. As a result, rice growth in the FWs was not only poor, it was also uneven. In the short term, land consolidation disturbed the nutrient status of the surface soil and thus negatively affected rice growth [2,44]. Previous studies have shown that soil physicochemical properties can have extremely significant effects in the short term after land consolidation [21,74], which include the severe compaction of soil, increased viscosity, decreased fertility, and decreased plant growth. Moreover, these adverse effects mainly occurred in the first growing season after land consolidation [30,31], when the soil layer was disturbed, especially the surface layer, and the topsoil physicochemical properties changed [20,30]. It is generally accepted that soils are affected by mechanical mixing in the short term, because the quality decreases [6,28]. This is similar to the results of TC and BT in this study.

Mitigating the Disturbance of Land Leveling to Soil
Ecological civilization has become the mainstream of development in China [22,75]. Land consolidation should change from adding the quantity of cultivated land and gain yield capacity to ecological rehabilitation, which focuses on the multifunctional ecological services of farmland [1,13]. This study revealed that land consolidation seriously disturbs the soil environment in the short term. This is consistent with the conclusions of previous studies [2,28]. The original soil configuration can be destroyed by land leveling, and the nutrient losses that occur when the topsoil is impacted are serious [21]. The current leveling method is unscientific; the regulations for land consolidation in China stipulate that first, the topsoil layer should be peeled off and set aside before leveling so that after the area is leveled, it can be backfilled with the same topsoil. However, rarely is the regulation strictly implemented. To minimize construction time and transportation costs, it is more common to cut directly high and fill low, resulting in a thinner topsoil and the loss of soil nutrients. Soil protection is neglected during leveling, and a large amount of topsoil is lost, causing soil quality declines in the short term. Our study suggests that it will be particularly important for FW and TC to implement land-leveling practices that strictly comply with the regulations in the future. For FW, it should be considered that there will be subsidence after the first irrigation, and that the thickness of borrowed soil should be thickened as needed according to the situation. For TC, the topsoil should first be removed (0-25 cm), the site shall be uniformly leveled by taking the high fill and low, and then, the original mature surface soil shall be backfilled. If necessary, the raw humus shall be used to cover the soil once more [2]. These measures are conducive to the internal balance of the topsoil and the equalization of yield capability.
National regulations in land consolidation are not comprehensive, and only specify requirements for soil thickness after land consolidation without requiring the maintenance of soil properties or quality. This work shows that there is value in regulations taking into account the value of multifunctional ecological services [13]. Land consolidation generally promotes the change of regional farmland ecological environment, but there is still a certain degree of quality improvement after land consolidation compared with high-quality farmland [30]. Follow-up management should include ecological control, scientific fertilization, and the application of green manure and humus fertilizer. Further, ecological management should be combined with the local soil, climate characteristics [14], and crop growth needs to formulate the best fertilization plan and usage. According to the actual conditions of different leveling areas, a long-term mechanism would be established to achieve a win-win situation for agricultural sustainable development and ecological protection [2,9].

Conclusions
As a powerful tool for land-use management, land consolidation plays a positive role in increasing the quantity of farmland, enhancing farmland infrastructure, and improving grain yield capacity. However, land consolidation has also strongly disturbed soil environment in the short term. It has also not yet received sufficient attention in China. This study was conducted in the land consolidation project of Gengzhuang Village in Yangtze River Delta. Remote sensing and field investigations were used to monitor the effects of land consolidation on soil properties and rice growth. The work yielded the following conclusions: (1) Land consolidation can seriously disturb soil nutrients in the short term. Land leveling significantly decreased SOM, TN, AK, and AZn in the short term, whereas ASi increased significantly. Soil nutrients differed among the different land-leveling types, which are generally represented as CA > BT > FW > TC.
(2) There was a significant difference in the growth status of rice in the first growing season under different land-leveling types, which are shown as CA > BT > FW > TC. The mean value of the NDVI of FW was 0.763 ± 0.03, TC was 0.747 ± 0.06, BT was 0.808 ± 0.05, and CA was 0.834 ± 0.02. NDVI low-low clusters areas were concentrated in TC and FW.
(3) The NDVI value of rice growth was significantly positively correlated with SOM and AZn, and the correlation coefficients were 0.831 and 0.810, respectively. The spatial variation of SOM and AZn caused by land leveling were the main factors leading to rice NDVI heterogeneity.
(4) NDVI values could replace field soil investigations and facilitate the large-scale monitoring of soil quality after land consolidation. Using remote sensing found that there was a strong disturbance after land consolidation in TC and FW. Follow-up management should give priority attention, such as through fertilizing soil and strengthening field management, which could help the soil retrieve its original state after a period of time. Remote sensing has the appropriate tools to identify the defective area of crop growth on a large scale, and provide a quick and convenient mean for the planning and construction of land consolidation, which can optimize land management.