**Farming Influence on Physical-Mechanical Properties and Microstructural Characteristics of Backfilled Loess Farmland in Yan'an, China**

**Lina Ma 1,2,3, Shengwen Qi 1,2,3,\*, Bowen Zheng 1,2,3, Songfeng Guo 1,2,3,\*, Qiangbing Huang <sup>4</sup> and Xinbao Yu <sup>5</sup>**


Received: 27 May 2020; Accepted: 30 June 2020; Published: 8 July 2020

**Abstract:** A gigantic project named Gully Land Consolidation (GLC) was launched in the hill-gully region of the Chinese Loess Plateau in 2011 to cope with land degradation and create new farmlands for cultivation. However, as a particular kind of remolded loess, the newly created and backfilled farmland may bring new engineering and environmental problems because the soil structure was disturbed and destroyed. In this study, current situations and characteristics of GLC are introduced. Test results show that physical-mechanical properties and microstructural characteristics of backfilled loess of one-year and five-year farmland are significantly affected by the Gully Land Consolidation project. Compared to natural loess, the moisture content, density, and internal friction angle of backfilled loess increase. On the contrary, the porosity, plasticity index, particle size index, and cohesion index decrease. Through SEM tests, it is observed that the particles of backfilled loess are rounded, with large pores filled with crushed fine particles, which results in skeleton strength weakness among particles and pores. The pore size distribution (PSD) of the four types of loess (Q3 loess, Q2 loess, one-year farmland, and five-year farmland) was measured using mercury intrusion porosimetry (MIP) tests, showing that the pore size of Q3 loess is mainly mesopores 4000–20,000 nm in size, accounting for 67.5%. The Q2, five-year, and one-year farmland loess have mainly small pores 100–4000 nm in size, accounting for 52.5%, 51.7%, and 71.7%, respectively. The microscopic analysis shows that backfill action degrades the macropores and mesopores into small pores and micropores, leading to weak connection strength among soil particles, which further affects the physical-mechanical properties of loess. The disturbance of backfilled loess leads to an obvious decrease in cohesion and a slight increase in internal friction compared to natural loess. The farming effect becomes prominent with increased backfill time, while the loess soil moisture content increases gradually. Both the cohesion and internal friction of the backfilled loess soil decrease to different degrees. This study is helpful to investigate sustainable land use in the Chinese Loess Plateau and similar areas.

**Keywords:** Gully Land Consolidation; backfilled loess; physical-mechanical property; microstructural characteristic; pore size distribution

#### **1. Introduction**

The Chinese Loess Plateau (CLP) covers an area of 640,000 km2 in the upper and middle reaches of the Yellow River. It also has multiple crisscrossing hills and gullies with fractured geological structures [1]. Meanwhile, the semiarid climate, with only 464.1 mm of average annual precipitation, contributes to the fragility of the ecosystem in the CLP [2]. Owing to aridity, sparse vegetation, and concentrated rainstorms, the CLP is facing serious problems of soil erosion and water loss, becoming the most vulnerable ecological environment in China [3–5].

Since 1949, the ecological management of the CLP has gone through several important stages, such as "slope management," "integration management of gully and slope," "small watershed management," and the "Grain for Green (GFG) project," which increased vegetation coverage and improved ecological quality [6,7]. In particular, since the Chinese government implemented the GFG project in 1999 to convert farmland to forests, shrubland, and grassland, the vegetation coverage rate of the CLP increased significantly, from 31.6% in 1999 to about 65% in 2017 [8], thereby reducing soil erosion and water loss to a certain extent [9,10]. Most notably, the land consolidation measures dramatically transformed the ecology and landscape of the CLP, turning it from yellow to green [2,11]. However, in recent years, Chinese scientists and government officials have noticed problems: There have been significant reductions in arable land induced by the large-scale implementation of the GFP on the Loess Plateau. Especially in Yan'an, the total area of arable land has reduced by more than half, creating an urgent need for new farmland and a shortage in grain production, which are the main obstacles in sustainable agricultural development [12,13].

On the other hand, Yan'an city has 44,000 gullies that are 500 m or more in length and 20,900 gullies that are 1 km or more in length. These numerous gullies have a land creation capacity in Yan'an of about 1.5 million acres [11]. To offset the loss of cropland and cope with land degradation, Shaanxi Province, which has 13% of the area of the CLP, launched a megaproject called Gully Land Consolidation (GLC) in 2011 [2,14–16], aimed at creating 26.67 <sup>×</sup> 104 ha of farmland in creek valleys from 2011 to 2020, where the total investment has reached ¥30 billion [17]. The GLC project was officially listed as a major national land improvement project in 2013 [18]. Yan'an City created 3.33 × 104 ha of farmland from 2013 to 2017 [19]. The GLC project has become an important measure in increasing areas of cultivated land and to expand land resources. Meanwhile, it also effectively reduces the loss of soil and water in trenches, which can largely alleviate the problem of land shortage and support agricultural modernization [20].

As new comprehensive management, the main approaches of GLC include reshaping valleys by incising foot slopes, filling gullies and stream channels, constructing or rebuilding drainage canals, dams, and reservoirs, and creating flat farmlands [2,11,19]. The backfilled loess from farmland in the GLC is a special type of remolded soil. After field investigations, although layered excavation and mechanical compaction were used in the construction process, these operations and management still lack unified and standard guides. Some essential indicators that affect the arability of farmland soils, such as the layered thickness of backfilled soils, tamping methods, and soil compaction, were not adequately demonstrated and implemented. This results in uneven compaction of newly constructed farmland and vast differences of microstructure among backfilled loess. The microstructure is the fundamental structural unit of soil, including grain morphology, grain contact patterns, pore characteristics, and pore size distribution (PSD) [21–24]. Different types of loess have significantly different microstructural characteristics, even if the same type also shows distinctly different microstructure under different pressure, moisture content, dry density, and other conditions. Furthermore, the macroscopic and mechanical properties of soils, such as shear behavior, compressibility, collapsibility, and permeability, are controlled by its microstructural characteristics [25–27]. Under the influence of rainfall and irrigation, collapse subsidence is likely to occur, causing some new engineering and environmental problems.

However, there are few studies on the engineering characteristics of backfilled loess in the GLC in China and also around the world. Based on this, this paper selected two gullies of Chunshuyaozi (CSYZ) and Shijiagou (SJG) as study areas, which were typical GLC sites in the "National Soil and Water Conservation Demonstration Park of South Gully." Four types of loess samples were collected with minimal disturbance, including natural samples of Middle Pleistocene Lishi loess (Q2) and Upper Pleistocene Malan loess (Q3) from the excavated slope and backfilled samples from two farmlands located in the reshaped CSYZ and SJG gullies. Farmland in the CSYZ gully has been backfilled for five years, and that in the SJG gully for one year. Moreover, a series of routine soil tests were performed on the four types of loess samples, including basic physical property tests, direct shear tests, and compression tests. Meanwhile, this study also investigated the microstructural characteristics of the loess by conducting scanning electron microscope (SEM) and mercury intrusion porosimetry (MIP) tests. By comparing the physical-mechanical properties and microstructural characteristics of four types of loess samples, this paper shows the evolution of backfilled loess and provides a basis for sustainable land use in the CLP.

#### **2. Materials and Methods**

#### *2.1. Study Site*

Chunshuyaozi (CSYZ) gully is located in Chunshuyaozi village, Louping town, Ansai district, Yan'an city, Shaanxi Province, and is where the GLC project was implemented in 2013. After the farmland construction was completed in one year, the farmland was immediately planted. So far, it has been five years since CSYZ farmland was backfilled. Figure 1a,b shows satellite images of CSYZ before and after the GLC, taken on 21 July 2010 and 12 March 2019, respectively. Through comparison, it was found that the morphology of these slopes and trenches changed significantly after excavation and backfill.

**Figure 1.** Satellite images of Chunshuyaozi (CSYZ) gully (**a**) before and (**b**) after the Gully Land Consolidation (GLC) project. (**c**) Morphology of CSYZ gully.

The main direction of the CSYZ gully is north–south, with the head on the south side (Figure 1c). The farmlands are composed of three platforms and three sedimentary dams. The first-level platform is elongated and irregular, located at the south of the gully, and is 40 m wide in the south and 100 m wide in the north. A field survey showed that the major plant was alfalfa. Close to the north side of the first-level platform is the first-level dam, which is about 2.5 m high and nearly 50 m wide. The second-level platform is irregular and polygonal, with a length of 300 m in the east–west direction and about 130 m in the north–south direction. The main crop in this platform is black wolfberry. The secondary dam on the north side of the secondary platform has a height of about 15 m and a length of about 230 m in the east–west direction and consists of four steps. The third-level platform is an irregular square and resembles a dustpan. It is about 200 m long from north to south and about 170 m from east to west, and is mainly planted with canola. The gully mouth area is the third dam with a length of about 300 m from east to west and a height of about 5 m. Intersecting with the CSYZ gully is the Shendao gully in a nearly NW direction.

The Chang'an University team revealed the main stratigraphy of the CSYZ gully from new to old: Q3 loess, Q2 loess, Hipparion red soil of the Pliocene Series, and underlying sandy mudstone of Zaoyuan section, Yan'an Formation in the lower Jurassic (J1). The thickness of Q3 loess was mainly between 15 and 25 m, and some parts can reach 30 m. The thickness of Q2 loess ranged from 30 to 60 m. Loess in gullies mainly included backfilled loess of Q3 and Q2, between 10 and 25 m depth. The groundwater level of the original slope in this study area was 72 m deep. The thickness of backfilled farmland was between 12 and 22 m due to the undulations of the original channel [28].

Shijiagou (SJG) gully, located in the east of South Gully village, Louping township, Ansai district, Yan'an city, Shaanxi Province, began to implement the GLC project in 2019. Figure 2a,b shows satellite images of SJG before and after the GLC, taken on 13 February 2010 and 12 March 2019, respectively. The direction trend of SJG gully is NE 45◦, with a length of 600 m. The backfilled farmlands are composed of five platforms and five sedimentary dams. A field survey showed that the major plant in SJG was alfalfa. The SJG gully is just 1.5 km away from the CSYZ gully in the southeast direction, and thus the strata are similar.

**Figure 2.** Satellite images of Shijiagou (SJG) gully (**a**) before and (**b**) after GLC.

#### *2.2. Soil Sampling*

Loess soil samples were taken from 14 sampling sites of the excavated slopes and the backfilled farmlands in the study area, eight sites in the CSYZ gully, and six sites in the SJG gully. Two parallel samples were taken from each site. For the CSYZ gully, the sampling site labeled Q3 is Q3 loess located in the middle-upper part of the southern slope, which is about 15 m from the crest of the slope in the vertical direction and 1 m from the slope surface in the horizontal direction. The sampling site labeled Q2 is Q2 loess located in the lower part of the western slope, which is about 55 m from the crest of the slope in the vertical direction and 1 m from the slope surface in the horizontal direction. The sampling sites labeled F-C-1 to F-C-6 are backfilled loess, and samples were taken from about 0.5 m depth. Among them, F-C-1 is on the northern side of the second-level farmland, F-C-2 and F-C-3 are on the middle-upper and lower part of the secondary dam, respectively, and F-C-4 to F-C-6 are on the southern side of the three-level platform, from east to west (Figure 3).

**Figure 3.** Plan of CSYZ gully and sampling location.

For the SJG gully, sample sites labeled F-S-1, F-S-2, and F-S-3 are located on the west side of the fifth-level platform, respectively, from south to north, and sample sites labeled F-S-4 to F-S-6 are on the south side, from west to east. The sampling depth was about 0.5 m (Figure 4).

**Figure 4.** Plan of SJG gully and sampling location.

All soils were sampled by an undisturbed sampling method. Through experience, we found that the advantages of PVC sampling cylinders are that they are sealable, impervious, and lightweight compared to traditional iron cylinders. PVC pipes were cut into short tubes 16 cm in height and 11 cm in diameter in advance. Cardboard plates 11 cm in diameter were made from cardboard boxes to be used as caps for the PVC tubes. The sampling process and method are as follows (Figure 5): First, about the top 0.5 m surface layer was removed at the sampling site, and the soil was excavated to make an intact soil block about 25 × 20 × 20 cm (Figure 5a,b). A PVC tube was placed on top of the soil block and gently pushed into the block while the surrounding soil was trimmed off (Figure 5c,d). Then, multiple layers of preservative film were used to completely seal the PVC tube in order to preserve the soil moisture from evaporation loss. Cardboard plates were taped to both ends of the tube to secure the soil sample and prevent disturbance (Figure 5e). The tube was marked with the sample number and placed in a box to be transported to the laboratory.

**Figure 5.** Sampling process: (**a**) Excavation of soil block from sampling site, (**b**) cutting to suitable size, (**c**) pressing the soil sample, (**d**) sampling by PVC pipe, (**e**) sealing the soil sample.

#### *2.3. Experimental Method and Design*

Tests of basic physical properties on 28 soil samples from 14 sampling sites were conducted, including moisture content, density, specific gravity, liquid limit, plastic limit, and particle analysis, and the following properties were obtained through calculation: Dry density, porosity, saturation, plasticity index, and characteristic index of granularity. All test methods were based on the Geotechnical Test Method Standard GB/T 50123-1999 [29]. Specific gravity was tested with the pycnometer method, and the liquid–plastic limit was tested using the liquid–plastic limit combined determination method. The particle analysis tests were performed by a Malvern Mastersizer 3000 laser particle size analyzer. Samples of Q3, Q2, F-C-4, and F-S-2 were selected for consolidated and undrained direct shear tests. A total of 14 soil samples of natural and backfilled loess were tested for consolidation and compression. The porosity ratios were obtained, and coefficients of compressibility and modulus of compressibility were calculated. According to the regulations of the test standard, each site should have two specimens, and if the difference between the two values is within the error range, the average of the two is taken as the result. According to the tests, properties of natural Q2 and Q3 loess are relatively stable, and the values of the two specimens are very close. For SEM and MIP tests, samples of F-C-4 and F-S-2 were chosen to represent CSYZ and SJG loess.

A compression experiment is generally used to study the compressibility of soil. The compression tests of all samples were carried out with GZQ-1 pneumatic automatic consolidometer. Coefficient of compressibility (α) and modulus of compressibility (*E*s) were respectively calculated by Equations (1) and (2):

$$\alpha = \frac{\varepsilon\_1 - \varepsilon\_2}{p\_2 - p\_1},$$

where *p*<sup>1</sup> is the original pressure, defined as 100 kPa; *p*<sup>2</sup> is the summation pressure, defined as 200 kPa; and *e*<sup>1</sup> and *e*<sup>2</sup> are corresponding void ratios.

$$E\_s = \frac{1 + c\_1}{\alpha},\tag{2}$$

Scanning electron microscopy (SEM) and mercury intrusion porosimetry (MIP) tests are two well-accepted and effective methods of studying the microstructure characteristics of soils. SEM can obtain microscopic images of the soil surface, which are generally used for qualitative analysis of soil structures. MIP is generally used for quantitative analysis of pore characteristics by measuring the size and distribution of micropores in the soil. In this paper, SEM and MIP tests were performed on loess of Q3, Q2, CSYZ, and SJG. The micrographs of particle morphology were analyzed by means of SEM, and the pore size distribution characteristics were measured by MIP. The characteristics of micromorphology, soil grain morphology, grain contact patterns, types of porosity, and pore size distribution were comprehensively analyzed, and the evolution process of microstructure characteristics of backfilled farmland loess with increased backfilled time was discussed.

The apparatus used in SEM tests was an FEI Nova Nano SEM field emission scanning electron microscope. The test steps were as follows: (1) The sample was cut into a cuboid with dimensions of 2 cm × 1 cm × 1 cm, air-dried in a light-shielded place, then split into two parts by hand along the middle position, and a flat and fresh side was selected as the observation surface. (2) After the sample was evacuated under vacuum pumping, the surfaces were treated with gold spray. Then the sample was fixed on the target plate with conductive tape and transferred to the observation room. (3) The instrument software was opened to find the right position to observe the sample. The focus was adjusted to obtain optimal images with appropriate brightness at selected degrees of magnification. This experiment was observed at four magnifications: 200×, 400×, 800×, and 1600×. In this paper, representative images of 200× and 800× were selected for analysis.

The test principle of MIP is as follows. This method is based on the assumption that the measured pores are cylindrical. As a noninvasive liquid with liquid surface tension, mercury cannot enter pores without pressure. During the mercury intrusion process, the work of moving mercury into the pores is converted to overcome the work done by the surface tension. Hence the Washburn equation for the external pressure and pore radius was deduced [30] as Equation (3):

$$P = -\frac{2\gamma\cos\theta}{r},\tag{3}$$

where *P* is the applied pressure (kPa); *r* is the pore radius (nm); γ is the surface tension of mercury, 0.485 N/m; and θ is the contact angle of the mercury and the surface of the soil grain, 140◦. It can be known from the formula that pressure is inversely proportional to pore size. The greater the applied pressure, the smaller the pore radius that mercury can enter. Mercury intrusion starts from large pores, then mesopores, and finally small pores. By the volume of mercury intrusion at a certain pressure, the volume of the corresponding pore sizes can be determined [31]. The instrument used in this test was a Micromeritics Auto Pore IV 9500, which can provide a pressure range of 3–4 <sup>×</sup> 105 kPa, with a corresponding pore size range of 3 nm–450 μm, and the accuracy is 1 μL. The samples were cut into cuboids of 1 × 1 × 1 cm and air-dried in a light-shielded place.

#### **3. Results**

#### *3.1. Tests of Basic Physical Properties*

Basic physical test results can be summarized as follows (Table 1). The results (Figure 6a) show that the moisture content of F-C is between 16.5% and 18.4%, and F-S is between 19.1% and 24.4%. Compared with natural loess, there is a significant increase in moisture content of backfilled loess. At the same time, the moisture content of F-C is generally lower than that of F-S, and the scope of the former is less than that of the latter. The results of density tests show a similar tendency (Figure 6b). The density of F-C is about 1.77–1.95 g/cm3, and that of F-S is about 1.82–2.08 g/cm3. On average, the density of backfilled loess shows a significant increase compared with the natural loess, and the increase in SJG is greater than that in CSYZ.


**Table 1.** Basic physical indices of 14 loess groups.

**Figure 6.** Comparison of basic physical properties: (**a**) moisture content, (**b**) density, (**c**) porosity, and (**d**) plasticity index.

On the contrary, the porosity of backfilled loess is significantly lower than that of natural loess (Figure 6c). The porosity ratio of Q3 and Q2 is 0.97 and 0.74, respectively. The porosity ratio of F-C varies from 0.63–0.81, and that of F-S ranges from 0.58 to 0.79. The mean porosity of F-C increased compared with that of F-S. Similarly, the plasticity index of backfilled loess is also lower than that of natural loess from slopes (Figure 6d). The values of F-C samples are mainly distributed in the range of 8.4–9.8, while the values of F-S samples have an obvious decrease and are mainly distributed in the range of 6.3–8.5.

#### *3.2. Grain Size Distribution*

Figure 7 shows the results of grain analysis tests. It can be seen that the grain size accumulation curves of 14 samples is smooth and continuous, with a granularity distribution ranging from 0.523 μm to 111 μm. Among all curves, Q3 is located on the far left of the graph, Q2 and F-C are mostly concentrated in the middle part, and F-S is mostly distributed on the right. It can be seen from Table 2 that Q3 has four maximums of characteristic grain size indices, d60, d50, d30, and d10. Overall, although there are no significant differences, the indices of SJG are slightly higher than those of CSYZ.

**Figure 7.** Grain size accumulation curves of 14 samples.


**Table 2.** Characteristic indices of granularity of 14 loess groups.

**Note:** d10, d30, d50, d60 are values corresponding to 10%, 30%, 50%, and 60% finer by weight, respectively; Cu is uniformity coefficient, *<sup>C</sup>*<sup>u</sup> = *<sup>d</sup>*<sup>60</sup> *<sup>d</sup>*<sup>10</sup> ; and Cc is curvature coefficient, *<sup>C</sup>*<sup>c</sup> <sup>=</sup> *<sup>d</sup>*<sup>2</sup> 30 *d*10×*d*<sup>60</sup>

In addition, the uniformity coefficient (Cu) and curvature coefficient (Cc) are two important indices for evaluating grain grading. In the tests, all soil samples satisfied Cu ≥ 5 and Cc = 1–3, which indicates well-graded particle sizes of these soils. Comparing the indices of natural and backfilled loess, the latter decreased, and the decrease of CSYZ is greater (Table 2).

#### *3.3. Direct Shear Tests*

The shear strength indices of the four types of samples were obtained through direct shear tests (Table 3 and Figure 8). The results indicate that the cohesion of Q3, Q2, F-C, and F-S is 20.11 kPa, 34.51 kPa, 13.61 kPa, and 18.64 kPa and the internal friction angle is 23.7◦, 28.8◦, 31.0◦, and 33.8◦, respectively. Compared with natural loess, the cohesion of backfilled loess decreased, and the decrease of F-C is much larger than that of F-S, the internal friction increased, and the increase of F-C is less than that of F-S.


**Table 3.** Indices of shear strength for Q3, Q2, F-C, F-S samples.

**Figure 8.** Normal stress–shear strength curves of Q3, Q2, F-C, F-S samples.

#### *3.4. Compression Tests*

Using coefficients of compressibility and modulus of compressibility as criteria to evaluate the compressibility of soils, it can be found that loess of Q3, Q2, and F-S has an intermediate level of compressibility, and F-C is at a low level (Table 4). The porosity ratio of F-S is lower than that of Q3, but compressibility is greater (Figure 9), indicating that the compressibility of soil is related not only to the total pore volume, but also to the age and water content of the soil, the shape of soil particles, the strength of soil skeletons, and the pore distribution. Therefore, it is necessary to further study the microstructural characteristics of loess.


**Table 4.** Indices of compressibility for 14 loess groups.

**Note:** α is the coefficient of compressibility and ES is the modulus of compressibility.

**Figure 9.** (**a**) Coefficient of compressibility, (**b**) modulus of compressibility.

#### *3.5. Scanning Electron Microscope Tests*

The SEM tests indicate that the microstructure of Q3 loess is loose, with clear grains and outlines (Figure 10a,b). The shape of the grains is granular, and surfaces presented with friable minerals. The particles mainly had point-to-point contact with each other. The intergranular pores are developed, in which scaffold pores are predominant, and the connectivity among pores is pretty good (Figure 10a,b). Compared to the structure of Q3, the structure of Q2 loess is much denser, and the particle outlines are blurred. The particle forms mainly include aggregates and clots, and surfaces are mainly cemented face-to-face. The pores of Q2 are reduced, and most are mosaic and intra-particle pores (Figure 10c,d). The grain distribution of the F-C backfilled loess is more uniform, and its structure is denser than that of the Q3 loess. The aggregated clay particles are mainly attached on the surfaces, showing that soil particles tended to be rounded when softened by water. The particles have both point-to-point and face-to-face contact patterns. F-C is dominated by mosaic pores and a small number of scaffold pores. It can be seen that some small particles filled in the large gaps (Figure 10e,f). The structure of the F-S backfilled loess is also dense, while the grain distribution is not uniform. Similarly, the aggregated clay particles are mainly attached on the surfaces, and the particles have both point-to-point and face-to-face contact patterns. There are more mosaic and intra-particle pores. High magnification shows that there is structural damage among the soil particles, and some flocculent particles caused by water softening filled into the pores (Figure 10g,h).

**Figure 10.** SEM images of Q3, Q2, F-C, and F-S samples: (**a**) Q3 under 200×, (**b**) Q3 under 800×, (**c**) Q2 under 200×, (**d**) Q2 under 800×, (**e**) F-C under 200×, (**f**) F-C under 800×, (**g**) F-S under 200×, and (**h**) F-S under 800×.

#### *3.6. Mercury Intrusion Porosimetry Tests*

Figure 11 shows the curves of cumulative intrusion vs. pore size, and the volume of cumulative intrusion represents the cumulative volume of the pore. It can be seen that the largest volume of mercury intrusion is 0.3041 mL/g for the Q3 loess, followed by 0.2420 mL/g for the Q2 loess, and 0.1883 mL/g and 0.1785 mL/g for F-C and F-S. The corresponding porosity ratio of the four groups is calculated as 0.77, 0.61, 0.48, and 0.46. It can be found that the porosity of the backfilled loess is significantly reduced, which is consistent with the SEM test results.

**Figure 11.** Cumulative intrusion vs. pore size of Q3, Q2, F-C, and F-S samples.

Moreover, the slope of the cumulative intrusion vs. pore size curve of Q3 loess is abrupt in the range of pore diameters from 6000 to 15000 nm, indicating that as the pore diameter decreased, the amount of mercury intrusion changed significantly. Therefore, it can be deduced that the pore diameter of sample Q3 is mainly distributed in the range of 6000–15,000 nm, and this range is called the dominant pore diameter. In the same way, the curve of F-S loess also has a stage of steep slope, and its dominant pore diameter ranges from 1500 to 2500 nm. However, the curves of Q2 loess and F-C loess are relatively gentle, and the dominant pore diameter is distributed in the ranges of 500–10,000 and 300–10,000 nm, respectively. The characteristics of the dominant pore diameter are not significant, and the pore distribution is relatively uniform. At the same time, it can be found that compared with natural loess, the range of dominant pore diameters for backfilled loess moves in the direction of the smaller pore size.

Figure 12 shows curves of pore size distribution (PSD) of four types of loess. Since PSD curves for Q3 and F-S have a significant peak, they are named single-peak curves [32,33]. The peak diameter of Q3 is 10,062 nm, and F-S is 1613 nm. However, for loess of Q2 and F-C, there are two peaks on the PSD curves, which are named double-peak curves. The two peak values of Q2 are 1615 and 6582 nm, and those of F-C are 832 and 8053 nm.

**Figure 12.** Pore size distribution curves of Q3, Q2, F-C, and F-S samples.

#### **4. Discussion**

SJG gully is about 1.5 km away from CSYZ gully, and the main strata are basically similar. There is not much difference in the elevations at the sampling sites between CSYZ and SJG, which are 1170–1185 m and 1168–1172 m. respectively. Hence, the assumption of this study is that the contribution of Q2 and Q3 loess in backfilled sites for CSYZ and SJG would be similar. A field survey showed that the major plant in SJG was alfalfa, and found that some crops are not growing well in the backfilled loess. Alfalfa is a plant with strong vitality in northern China that can adapt to a newly backfilled farmland environment well. Therefore, farmers always plant alfalfa on newly backfilled farmland to improve the soil. Based on the above, the effects of farming and backfilled time on soil properties were comparatively studied.

#### *4.1. Analysis of Physical-Mechanical Properties*

From Table 2 and Figure 6, the test results show that the moisture content and density of backfilled loess increased compared with natural loess. On the contrary, the porosity and plasticity indices decreased. With increased backfill time, almost all of the above-mentioned physical properties converged toward the natural loess. This reflects that backfilled soils tend to recover to their natural state when subjected to environmental processes such as wet/dry and seasonal cycles.

The analysis of grain size distribution (Figure 7) shows that there was breakage of a few large soil particles in the GLC process, and the mean grain size of backfilled farmland loess decreased (Table 2). Since the sampling depth of 0.5 m is within the influence range of plant cultivation [34–36], farming will affect the properties of soils. Some studies have shown the importance of soil functioning and soil–plant interactions [37–39]. For this study, under the influence of farm operations such as cultivation, irrigation, and digging, loess particle distribution tended to be homogenized.

Compared with natural loess, the cohesion of backfilled loess decreases, and the decrease of F-C is much larger than that of F-S, while the internal friction increases, and the increase of F-C is less than that of F-S. The reason may be that backfilled loess was disturbed in the GLC process, and the cementation and electrostatic attraction among soil particles were destroyed while the occlusal effects between particles were strengthened, so the cohesion decreased a little, and the internal friction angles increased a little. It should be noted that compared with F-C (with five years of farming), F-S (with one year of farming) has larger cohesion and internal friction as well, which indicates that with rainfall and irrigation, the water content gradually increased (Table 1), and the cohesion and internal friction angles gradually decreased under farming.

#### *4.2. Analysis of Microstructural Characteristics*

The SEM tests indicate that the grains of backfilled loess (F-C and F-S) are tightly compacted, and the distribution of grains and pores is much denser than that of natural loess (Q2 and Q3), thus backfilled loess has a relatively higher internal friction angle than natural loess. However, the cement between the particles has been damaged by disturbance, and the backfilled loess has relatively lower cohesion than natural loess.

There are many studies on the classification of pores in loess at home and abroad [40–43], but there is still no unified conclusion. In this study, we divided pores into four groups according to the PSD characteristics (Figure 12): Macropores (pore diameter >20 μm), mesopores (4–20 μm), small pores (0.1–4 μm), and micropores (<0.1 μm). It can be seen from Figure 12 that the four types of loess have similar volume contents of macropores and micropores, but obvious differences in the distribution of mesopores and small pores. Calculating the percentages of the four groups of pores according to the ratio of certain pore volume to total pore volume (Table 5), it can be seen that the Q3 loess is mainly dominated by mesopores of 4000–20,000 nm, accounting for 67.5%, and the Q2, F-C, and F-S loess have mainly small pores of 100–4000 nm, accounting for 52.5%, 51.7%, and 71.7%, respectively. This indicates that some large and medium pores in the soil were destroyed and turned into small pores and micropores as a result of the engineering of Gully Land Consolidation.


**Table 5.** Volume percentages of four categories of pores for Q3, Q2, F-C, and F-S samples.

#### *4.3. Suggestion for Gully Land Consolidation*

GLC changes the macroscopic properties and microstructure of loess, especially the moisture content, density, and pore size distribution. Soil moisture is a key factor influencing soil nutrient movement and soil quality in the semiarid Loess Plateau [44–46]. Meanwhile, the surface loess of backfilled farmland undergoes mechanical action, resulting in increased density and pore ratio. However, dense soil is not conducive to the growth of plant roots. Field surveys showed that most crops do not grow well in backfilled loess in the first few years. Alfalfa, as a kind of plant with strong vitality in northern China, can adapt to this new land environment well. Therefore, farmers switched to planting alfalfa on backfilled farmland in subsequent years. Under the influence of farmers' operations such as cultivation, irrigation, and digging, characteristics of backfilled loess tend to resemble those of natural loess, which will become loose and fertile. Hence, this study suggests that newly constructed farmland may not be suitable for planting crops. Farmers can make the soil fertile by planting vigorous grasses or shrubs, like alfalfa and robinia pseudoacacia. At the same time, it is important to control the compactness of soil during construction.

#### **5. Conclusions**

Test results show that physical-mechanical properties and microstructural characteristics of backfilled loess of one-year and five-year farmland are significantly affected by Gully Land Consolidation.

(1) Compared with natural loess, the moisture content, and density of backfilled loess increase. On the contrary, the pore ratio, plasticity index, and particle size index decrease. Additionally, with the development of filling time, the physical indices of backfilled loess tend to resemble those of natural loess, and the fluctuation amplitudes gradually decrease.

(2) The microstructure of soil samples observed by SEM tests indicates that natural loess has a certain skeleton strength with a relatively stable structure between the grains and pores. However, the distribution of particles and pores in backfilled loess (F-C, F-S) becomes denser, but the skeleton strength between the particles is destroyed, and the structure is more unstable. It was observed that the particles of backfilled loess are rounded with large pores filled with crushed fine particles, which results in weakness of the skeleton and cement strength among particles and pores and strengthening of the internal friction angle.

(3) By MIP tests and microscopic analysis, it can be concluded that some larger and medium pores in backfilled soils were destroyed and turned into smaller pores and micropores, and the cement strength between particles was damaged, which essentially affects the physical-mechanical properties of loess.

(4) Disturbance of backfilled loess leads to an obvious decrease of cohesion, and a slight increase in internal friction compared with natural loess. With increased backfilled time, the farming effect becomes prominent over time, the loess soil moisture content increases gradually, and both the cohesion and internal friction of the backfilled loess soil decrease to different degrees.

(5) Newly constructed farmland may not be suitable for planting crops in the first few years. Farmers can make the soil fertile by planting vigorous grasses or shrubs, like alfalfa and robinia pseudoacacia. At the same time, it is important to control the compactness of soil during construction.

**Author Contributions:** Conceptualization, S.Q. and S.G.; Data curation, L.M., B.Z., and S.G.; Formal analysis, L.M, S.Q., B.Z., S.G., and X.Y.; Funding acquisition, S.Q., S.G. and Q.H.; Investigation, L.M., S.Q., S.G., and Q.H.; Methodology, S.Q., S.G. and Q.H.; Project administration, S.Q., B.Z. and S.G.; Resources, S.Q.; Supervision, S.Q. and S.G.; Validation S.Q.; Writing—original draft, L.M.; Writing—review and editing, S.Q., S.G., Q.H., and X.Y. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the State Key Research Development Program of China under grant no. 2017YFD0800501, the Natural Science Foundation of China under grants nos. 41790442 and 41825018, and the Q.H.Strategic Priority Research Program of the Chinese Academy of Sciences under grant no. XDA23090402.

**Acknowledgments:** Special thanks to Zhang Bin at the China University of Geosciences, Beijing, and Yang Jijin at the Institute of Geology and Geophysics, Chinese Academy of Sciences, for their help with the experiments.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
