Effects of Grass and Forests and the Infiltration Amount on Preferential Flow in Karst Regions of China

Preferential flow is an important water infiltration phenomenon in karst regions. The response of preferential flow to vegetation restoration requires urgent investigation due to the special soil structure of karst regions. In order to study the effect of vegetation restoration on water movement in karst regions, four kinds of ponded water infiltration experiments were carried out in Pinus Yunnanensis plantation forestland, secondary forestland, and natural grassland. A brilliant blue dyeing experiment was conducted to visualize the distribution of water infiltration in soil (a total of 150 stained images from vertical soil slices). Results showed that the average depth of matrix flow in natural grassland was approximately six times those in plantation and secondary forestlands. An increase in matrix flow will have a negative effect on the development of preferential flow. Water transported in preferential flow paths affects the distribution of nutrients and organic matter in the soil. However, preferential flow in grassland can promote the accumulation of available nutrients, and preferential flow in plantations can inhibit the loss of organic matter. Preferential flow in grasslands and forest plantations is less than that in native forests soils. Preferential flow increases the percolation of water in soils. The effect is that preferential flow can obstructs water uptake by the roots under low rainfall conditions, and decreases surface runoff before soil saturation under high rainfall conditions. In the process of nutrient element migration, preferential flow has a good contribution, which is conducive to the migration and accumulation of elements required for surface vegetation growth. The contribution of preferential flow needs to be considered in studies on vegetation restoration planning and land degradation.


Introduction
The concept of preferential flow is different from that of matrix flow [1].Water infiltration initially shows a wetting front parallel to the soil surface.This part is defined as matrix flow, which is classified as uniform flow [2].Then, the wetting front becomes irregular with different water inflow velocities.This part is called preferential flow, which is classified as non-uniform flow [3,4].Preferential flow is a ubiquitous soil water movement phenomenon [5].Preferential flow refers to the uneven and rapid movement of water and solutes through porous media (typically soil).It is an important factor that increases soil water conservation [6] and groundwater recharge [7] because of spatial connectivity of soil [8], which has been studied by certain mathematical models [9].Moreover, preferential flow, as an important part of soil water infiltration, can affect the transport of pollutants [10], thus influencing the capacity of the soil to buffer and filter [11].Some studies have found that the concentration of water-borne pollutants is not uniform when they diffuse [12]; however, few studies have focused on whether the distribution of soil nutrients will change with water infiltration.The special soil structure in karst area will affect the soil macropores, thus affecting preferential flow and soil permeability [13].Therefore, preferential flow can affect the transport of elements and conveyance system for groundwater [14], which are the key factors in vegetation restoration and land degradation [15].In this work, we used Brilliant Blue dyed tracer [16] and digital image processing techniques to determine the marks and analyses the characteristics of preferential flow, respectively [17].
The rapid process of preferential flow reduces the water retention time [18]; in order to make a profound study about the law of water cycle and the practical significance of natural forest and artificial forest on bottom of slope, the geochemical processes and the mode of vertical leakage need to be investigated [19].At present, most studies that focus on different water inputs have lots of effects on preferential flow under different soil architectures [1], land cover [20], and land use [21].Several studies show that a large rainfall of short duration contributes to the occurrence of preferential flow [22].Few studies on the transport of metal ion pollutants by preferential flow after rainfall [23] and snowpack [24] have been reported.Some mathematical models simulate the vertical movement process of solutes carried by water [25].However, insufficient attention has been given to preferential flow and its influence on water and nutrient transport under different precipitation levels in situ simulation experiments.The reason for the importance of studying on preferential flow is that it can easily form karst rock desertification due to the fragile ecological environment in karst regions; moreover, its environmental capacity is small, its anti-interference ability is low, and its resistance to pollutant diffusion is poor [26], thereby seriously restricting the development of agroforestry.Therefore, the effect of preferential flow on soil nutrient elements for vegetation restoration requires urgent study.Preferential flow has attracted wide attention because of its widespread existence, and its significant impact on the environment.Preferential flow has been studied based on dye tracer [27], CT scan [28,29], penetration curve [30], permeameter [31], and mathematical model [32].Among them, the dyeing tracer method has become the most widely used preferential flow research method because it can visually and clearly reflect the region through the dyeing water flows.As a dye tracer, Brilliant Blue can visually show the distribution of non-uniform flow (preferential flow) in soil [33].
The ecological environment in karst regions needs to be governed and restored immediately [34].However, the high variability in the soil spatial structures increases the difficulty of vegetation restoration.The pore structure and water infiltration characteristics of the soil can be obtained indirectly by studying preferential flow dyeing areas.In order to study the response of different forest sources (plantation restoration forestland and secondary forest or grassland) to the soil environment in this area, and discuss the influence of plantation on soil, experiments were carried out on Pinus Yunnanensis plantation forestland, secondary forestland and grassland in Jianshui, Yunnan Province, and the secondary forestland and grassland were selected for the study of primitive vegetation communities.A man-made forest (e.g., Pinus Yunnanensis plantation forestland) comprises excellent native tree species because of its extremely drought-resistant, barrenness resistance, and fast-growing under any conditions.
On this basis, the characteristics of soil erosion in karst areas are further investigated, and the advantages and disadvantages of vegetation restoration in karst areas are evaluated.These practices can provide theoretical support for the prevention and control of rock desertification.Therefore, this study aims (1) to analyze the distribution of preferential flow in a typical stand origin and the characteristics of water preferential infiltration movement, (2) to study the response of preferential flow to precipitation amounts and (3) to determine the relationship between nutrient accumulation and preferential flow in different typical stand origins.

Site Description
Field experiments were conducted in forest soil from July 2018 to August of 2018 in Yangjie Catchment (102 • 55 E, 23 • 37 -23 • 44 N) in the southwestern karst area of Jianshui Town in Yunnan Province, China.The karst regions in Southwest China are the largest continuous karst regions in the world [35].We selected three typical forest lands (e.g., Pinus Yunnanensis plantation, secondary forestland, and natural grassland) at the foot of the slope.The landform is that of a typical karst basin.The altitude is 1350-1700 m.It is affected by the monsoons of the Southwest Indian Ocean, which can form distinct dry and wet seasons.The seasonal drought of this area is severe, indicating that the rainy season (May to October) is humid, and rainfall during the dry season (November to April of the following year) is scarce.The annual average temperature is 19.8 • C, and the annual average temperature of the land surface is 20.8 • C. The average annual sunshine hours are 2322 h, and the average annual precipitation is 805 mm.The plantation forestry was planted in 1985, and its recovery time has been more than 30 years.Thus, the forest soil structure is stable.Secondary forestland has been restored since 1970, and no human disturbance was observed in the secondary forestland and grassland.Basic descriptions of the experimental sites are provided in Tables 1 and 2.

Site Description
Field experiments were conducted in forest soil from July 2018 to August of 2018 in Yangjie Catchment (102°55′ E, 23°37′-23°44′ N) in the southwestern karst area of Jianshui Town in Yunnan Province, China.The karst regions in Southwest China are the largest continuous karst regions in the world [35].We selected three typical forest lands (e.g., Pinus Yunnanensis plantation, secondary forestland, and natural grassland) at the foot of the slope.The landform is that of a typical karst basin.The altitude is 1350-1700 m.It is affected by the monsoons of the Southwest Indian Ocean, which can form distinct dry and wet seasons.The seasonal drought of this area is severe, indicating that the rainy season (May to October) is humid, and rainfall during the dry season (November to April of the following year) is scarce.The annual average temperature is 19.8 °C, and the annual average temperature of the land surface is 20.8 °C.The average annual sunshine hours are 2322 h, and the average annual precipitation is 805 mm.The plantation forestry was planted in 1985, and its recovery time has been more than 30 years.Thus, the forest soil structure is stable.Secondary forestland has been restored since 1970, and no human disturbance was observed in the secondary forestland and grassland.Basic descriptions of the experimental sites are provided in Tables 1 and 2.

Site Description
Field experiments were conducted in forest soil from July 2018 to August of 2018 in Yangjie Catchment (102°55′ E, 23°37′-23°44′ N) in the southwestern karst area of Jianshui Town in Yunnan Province, China.The karst regions in Southwest China are the largest continuous karst regions in the world [35].We selected three typical forest lands (e.g., Pinus Yunnanensis plantation, secondary forestland, and natural grassland) at the foot of the slope.The landform is that of a typical karst basin.The altitude is 1350-1700 m.It is affected by the monsoons of the Southwest Indian Ocean, which can form distinct dry and wet seasons.The seasonal drought of this area is severe, indicating that the rainy season (May to October) is humid, and rainfall during the dry season (November to April of the following year) is scarce.The annual average temperature is 19.8 °C, and the annual average temperature of the land surface is 20.8 °C.The average annual sunshine hours are 2322 h, and the average annual precipitation is 805 mm.The plantation forestry was planted in 1985, and its recovery time has been more than 30 years.Thus, the forest soil structure is stable.Secondary forestland has been restored since 1970, and no human disturbance was observed in the secondary forestland and grassland.Basic descriptions of the experimental sites are provided in Tables 1 and 2.

Site Description
Field experiments were conducted in forest soil from July 2018 to August of 2018 in Yangjie Catchment (102°55′ E, 23°37′-23°44′ N) in the southwestern karst area of Jianshui Town in Yunnan Province, China.The karst regions in Southwest China are the largest continuous karst regions in the world [35].We selected three typical forest lands (e.g., Pinus Yunnanensis plantation, secondary forestland, and natural grassland) at the foot of the slope.The landform is that of a typical karst basin.The altitude is 1350-1700 m.It is affected by the monsoons of the Southwest Indian Ocean, which can form distinct dry and wet seasons.The seasonal drought of this area is severe, indicating that the rainy season (May to October) is humid, and rainfall during the dry season (November to April of the following year) is scarce.The annual average temperature is 19.8 °C, and the annual average temperature of the land surface is 20.8 °C.The average annual sunshine hours are 2322 h, and the average annual precipitation is 805 mm.The plantation forestry was planted in 1985, and its recovery time has been more than 30 years.Thus, the forest soil structure is stable.Secondary forestland has been restored since 1970, and no human disturbance was observed in the secondary forestland and grassland.Basic descriptions of the experimental sites are provided in Tables 1 and 2.

Dyeing Experiment
The sample plots were selected in the karst fault basin in Jianshui County, Yunnan Province.Sample plots were selected to represent local vegetation types (Pinus Yunnanensis plantation forestland, secondary forestland, and natural grassland) and growth conditions and have relatively consistent site conditions (e.g., elevation, slope position, and slope).Secondary and plantation forests with a stable soil structure and reduced disturbance were selected [33].The recovery time of the three kinds of woodlands was more than 20 years.A natural secondary forest occurs from secondary bare land through a series of successive stages of the local plant community.And local tree species Pinus Yunnanensis in plantation accounted for a large proportion.
To study the preferential flow characteristics for different stand origins under varied precipitation conditions, we selected flat fields in plantation forestland, secondary forestland, and natural grassland with two replicated transects, under similar environmental characteristics at least three weeks after rainfall to remain dry prior to the dye experiment.Each transect included five vertical sections, and we requested the sample sites to be located at least 10 m from each other to avoid interferences in the experiments [36].Before the test, we removed fallen leaves and other coverings from the surface to prevent water from being absorbed by elements other than the soil.At each site, plastic sheets with a length × width × height of 0.6 m × 0.6 m × 0.3 m were used to surround the field.Subsequently, the plastic sheets were embedded 10 cm into the soil, and the soil around the plastic sheets was uniformly compacted at 5 cm.
As a dye tracer, Brilliant Blue can visually show the distribution of non-uniform flow in soil, that is, preferential flow [16].Therefore, a 4 g/L dye solution was used in 36 tracing experiments [37].According to rainfall data for the years 2008-2017 in our study area, the distribution of rainfall amounts was: 69% of the rainfall events for <10 mm rain (light rain); 23% of events for 10-25 mm (moderate rain); 7% of events for 25.1-50 mm rain (heavy rain), and 2% of events for >50.1 mm rain (a storm).In order to simulate and describe water infiltration phenomena under different conditions of surface water accumulation using an in-situ test, we used the one-dimensional ponded water infiltration method to simulate the change of ponded water in four kinds of rainfall conditions, which can reduce water loss and minimize the wastage.Our simulation experiments included rainfall amounts of 5 mm (a small rainfall event), 15 mm (median), 35 mm (large), and 55 mm (extreme), because the probability of these precipitation events occurring in four kinds of rainfall amounts is higher.The maximum rainfall depth recorded at the study site was 67 mm during 24 h.Therefore, the applied amounts, especially 35 mm and 55 mm, should be reasonable.The fields, images and soil samples were numbered G5, G15, G35, and G55, according to the amount of Brilliant Blue solution used.Then, each field was covered with a waterproof cloth that was removed 24 h after the cessation of infiltration to avoid evaporation of the solution [36].We removed the waterproof cloth and left a 10 cm buffer zone at the inner edge of each sheet to prevent the cutting process from affecting the pore structure.The final data were collected in the 50 cm × 50 cm staining area in the middle of the quadrant, and each quadrant was cut into five vertical sections with 10 cm spaces between each other.We photographed each vertical section using a digital camera (500D, Canon, Japan) with a 35 mm focal length.In cases of bright sunlight, we used a parasol to reduce shadow interference [38].
Soil samples were taken at depths of 0-10 cm, 10-20 cm, and 20-30 cm.The samples were taken back to the laboratory of graphite digestion instrument to test the chemical properties of soil after digestion with concentrated sulfuric acid and perchloric acid.The samples were assayed by automatic chemical analyzers (the model is Smartchem).Available potassium (AK) used flame photometry method, nitrate nitrogen (NN) used hydrazine sulfate method, total nitrogen (TN) used salicylic acid-hypochlorite spectrophotometry method, available phosphorus (AP) used NaHCO 3 method, total phosphorus (TP) used the ammonium molybdate method, organic matter (Org) used potassium dichromate oxidation method.

Image Analysis
The photos were processed in Photoshop CS3 (Adobe Systems Inc., San Jose, CA, USA) and Image Pro Plus 6.0 (Media Cybernetics Inc., Rockville, MD, USA) to quantify the stained areas [39].Image analysis geometrically corrected the captured images by using Photoshop CS3.The sizes of the vertical sections of the pictures were cropped to 50 cm × 50 cm by referring to the image ruler.Then, the saturation, brightness, greyscale, and threshold values were adjusted.The processed images were digitized using Image Pro Plus 6.0, and the pixel of each image was set to 500 × 500 to convert the image into binary information, classified into black (pixel value 0, represents the dye position) and white (pixel value 255, represents the undyed position).We calculated preferential flow indices by using Sigmaplot to calculate these binary matrices (see Figure 1 for visual examples).
Water 2019, 11, X FOR PEER REVIEW 5 of 20 the vertical sections of the pictures were cropped to 50 cm × 50 cm by referring to the image ruler.Then, the saturation, brightness, greyscale, and threshold values were adjusted.The processed images were digitized using Image Pro Plus 6.0, and the pixel of each image was set to 500 × 500 to convert the image into binary information, classified into black (pixel value 0, represents the dye position) and white (pixel value 255, represents the undyed position).We calculated preferential flow indices by using Sigmaplot to calculate these binary matrices (see Figure 1 for visual examples).

Data Analysis
Dye coverage (DC/%) refers to the percentage ratio of the dyed area [40].The dyeing area ratio is mainly obtained by the ratio of the number of dyed pixels at different depths in each soil layer to the total number of pixels in the image to reflect the degree of preferential flow, which is calculated as: where DC corresponds to DC (%), D is the stained area of the soil profile (cm 2 ) and ND is the unstained area of the soil slices (cm 2 ).Preferential flow fraction (PF) is used to reflect the distribution and degree of preferential flow [41].A high PF value indicates a high degree of preferential flow [42], which is calculated as: where PF is the soil preferential flow fraction (%); UF is the depth of uniform flow (cm), that is, matrix flow; W is soil profile horizontal width (cm); DT is the total dyed area of the soil sections (cm 2 ).The length index (Li) of preferential flow has a positive correlation with the degree of heterogeneity of preferential flow [43] and is calculated as: where Li is length index of preferential flow (%), DCi is DC area ratio to the i-th layer (%), and n is the number of pixels in the vertical sections (500 in this study).
The variation coefficient of the DC in the preferential flow area (CV) is used to quantitatively describe the change in dyeing morphology after the removal of the matrix flow areas [44] and is calculated as:

Data Analysis
Dye coverage (DC/%) refers to the percentage ratio of the dyed area [40].The dyeing area ratio is mainly obtained by the ratio of the number of dyed pixels at different depths in each soil layer to the total number of pixels in the image to reflect the degree of preferential flow, which is calculated as: where DC corresponds to DC (%), D is the stained area of the soil profile (cm 2 ) and ND is the unstained area of the soil slices (cm 2 ).Preferential flow fraction (PF) is used to reflect the distribution and degree of preferential flow [41].A high PF value indicates a high degree of preferential flow [42], which is calculated as: where PF is the soil preferential flow fraction (%); UF is the depth of uniform flow (cm), that is, matrix flow; W is soil profile horizontal width (cm); DT is the total dyed area of the soil sections (cm 2 ).The length index (L i ) of preferential flow has a positive correlation with the degree of heterogeneity of preferential flow [43] and is calculated as: where L i is length index of preferential flow (%), DC i is DC area ratio to the i-th layer (%), and n is the number of pixels in the vertical sections (500 in this study).
The variation coefficient of the DC in the preferential flow area (C V ) is used to quantitatively describe the change in dyeing morphology after the removal of the matrix flow areas [44] and is calculated as: where C V is the variation coefficient of DC in the preferential flow area (-), N is the number of soil sections of preferential flow in vertical sections (-), x is the average value of the staining area ratio (-), and x i is the dye area ratio of the soil profile in layer i-th (-).
The variation coefficients of the maximum infiltration depths (C µ ) are used to describe the non-uniform characteristics of the soil water flow in vertical slices [45] and are calculated as: where C µ is the variation coefficient of maximum infiltration depths (-), h max (x i ,y i ) is the maximum depth of dyeing located at (x i , y i ) of the vertical section (cm), h is the average of all maximum depths (cm), and m is the number of pixels in the horizontal direction of the vertical section (500 in this study).

Structural Equation Model (SEM)
SEM is widely used due to its empirical analysis capabilities and suitability for latent variable analysis [46].Compared with traditional statistical methods, SEM introduces latent variables that can simultaneously consider endogenous variables, which can allow errors of measurement in parameter estimation whilst analyzing the direct and indirect effects between variables [47].The SEM method consists of structural equations and measurement models that are generally represented by three matrix equations, as follows [48].
Equation ( 6) represents SEM, which reflects the structural relationships between latent variables, which cannot be measured in practice.Where η is the endogenous latent variable, B is the correlation coefficient matrix between endogenous latent variables, Γ is the structural coefficient matrix of exogenous latent variables on endogenous latent variables, ξ is an exogenous latent variable and ζ is the measurement error.
Equations ( 7) and ( 8) represent the measurement models, which can reflect the relationship between observed and latent variables, where X is the exogenous observation variable vector; Y is the endogenous observation variable vector; Λ x is the factor loadings of the exogenous observation variable vector; Λ y is the factor loadings of endogenous observation variable vector; ε and σ are the measurements of the exogenous and the endogenous observation variables, respectively.Latent variables can be reflected by measurable variables by the use of a measurement model.We can obtain the relationship between each latent variable and measurable variable by solving Equations ( 6)- (8).After building the initial model, the path coefficients represent the extent of the relationship between the variables.The path coefficients are usually calculated by using the maximum likelihood method.

Spatial Variation of Dyeing Distribution
Blue Brilliant FCF dye tracer was used to visualize the trace of water at the initial stage of infiltration (24 h), which is a common and effective dyeing reagent [40,49].We compared and calculated the data extracted from 150 stained vertical sections in three vegetational types, and nine examples are shown in Figure 2. The distribution of dye area was not uniform or similar in the three experimental fields (Figure 2).The uniform infiltration depth (cm), which was calculated as the depth in which DC was 80% and belongs to the matrix flow, can represent the depth of the uniform infiltration front [43].In Figure 2, the images in first line are taken from Pinus Yunnanensis forestland fields, the images in second line are taken from secondary forestland fields, and the images in third line are taken from grassland fields; a, e, i correspond to the images under G5 condition; b, f, j correspond to the images under G15 condition; c, g, k correspond to the images under G35 condition; d, h, m correspond to the images under G55 condition.The dyeing area ratio decreased as the soil depth increased [17,39].The degree of image change indicated the presence of preferential flow in this area and the dominance of matrix flow in the soil slices of grassland.Preferential flow was apparent at Pinus Yunnanensis plantation forestland and secondary forestland, but the spatial distribution of the secondary forestland was noticeably more non-uniform.At the same time, Figure 2, from a to d, e to h, i to m, showed that the infiltration depth increases with the increase of infiltration water volume.We can discuss the characteristic of spatial distribution of preferential flow by drawing the standard deviation images of DC (Figure 3).In general, the standard deviation of DC in Pinus Yunnanensis plantation forestland was more normal and concentrated; the standard deviation of DC in secondary forestland fluctuated significantly; a predominant downward trend in DC in grassland was observed mostly, and where the standard deviation was also the lowest.
Water 2019, 11, X FOR PEER REVIEW 7 of 20 line are taken from grassland fields; a, e, i correspond to the images under G5 condition; b, f, j correspond to the images under G15 condition; c, g, k correspond to the images under G35 condition; d, h, m correspond to the images under G55 condition.The dyeing area ratio decreased as the soil depth increased [17,39].The degree of image change indicated the presence of preferential flow in this area and the dominance of matrix flow in the soil slices of grassland.Preferential flow was apparent at Pinus Yunnanensis plantation forestland and secondary forestland, but the spatial distribution of the secondary forestland was noticeably more non-uniform.At the same time, Figure 2, from a to d, e to h, i to m, showed that the infiltration depth increases with the increase of infiltration water volume.We can discuss the characteristic of spatial distribution of preferential flow by drawing the standard deviation images of DC (Figure 3).In general, the standard deviation of DC in Pinus Yunnanensis plantation forestland was more normal and concentrated; the standard deviation of DC in secondary forestland fluctuated significantly; a predominant downward trend in DC in grassland was observed mostly, and where the standard deviation was also the lowest.We calculated the experimental data using the aforementioned formulas (Table 3).The dyeing area distinctly declined under different conditions.DC rapidly declined in the soil layer from 0 to 5 cm.The dyeing morphological parameters of preferential flow under different infiltration amounts in Table 3 showed that the depth of matrix flow was deepest in grassland, followed by those in Pinus Yunnanensis plantation forestland, and secondary forestland.This phenomenon indicated that the depth of matrix flow may be related to the near-surface vegetational types, water infiltration into the soil through fissures, macropores and pore networks formed by roots, and the higher abundance and uniformity of grass roots in surface soil than in forest land.
However, in the deep soil layer (i.e., 5-15 cm), a significant difference was observed in the DC of vertical sections under different field types; these large variations may be related to the pore structure of the soil, that is, the effect of integrated environmental factors, such as soil properties [50].Table 3 shows that the characteristic parameters of the preferential flow of secondary forestland were bigger than those of Pinus Yunnanensis plantation forestland and grassland.The reason for this was that the soil layer in Pinus Yunnanensis plantation forestland is thicker than in secondary forestland, but the main preferential flow path is the root system of trees, which influences the connectivity of macropores; thus, the standard deviation of Pinus Yunnanensis plantation forestland is higher than that of grassland, and this finding is consistent with that of Alaoui et al. [37].For the same reason, the forestland macropores can efficiently transport vertically downwards compared to those of the grassland soil.This phenomenon is important in soil hydrological processes (e.g., the soil water cycle and water infiltration rate).The preferential flow pathway of secondary forestland is extremely developed because of the following: (1) the special morphology and high heterogeneity of the soil structure in karst regions [35], (2) the fact that the fracture between rock and soil has become a preferential pathway for transporting surface water to deep soil layers [51], and (3) the non-uniform distribution of fissures in surface soil, which is caused by the lack of uniform vegetation cover.We calculated the experimental data using the aforementioned formulas (Table 3).The dyeing area distinctly declined under different conditions.DC rapidly declined in the soil layer from 0 to 5 cm.The dyeing morphological parameters of preferential flow under different infiltration amounts in Table 3 showed that the depth of matrix flow was deepest in grassland, followed by those in Pinus Yunnanensis plantation forestland, and secondary forestland.This phenomenon indicated that the depth of matrix flow may be related to the near-surface vegetational types, water infiltration into the soil through fissures, macropores and pore networks formed by roots, and the higher abundance and uniformity of grass roots in surface soil than in forest land.
However, in the deep soil layer (i.e., 5-15 cm), a significant difference was observed in the DC of vertical sections under different field types; these large variations may be related to the pore structure of the soil, that is, the effect of integrated environmental factors, such as soil properties [50].Table 3 shows that the characteristic parameters of the preferential flow of secondary forestland were bigger than those of Pinus Yunnanensis plantation forestland and grassland.The reason for this was that the soil layer in Pinus Yunnanensis plantation forestland is thicker than in secondary forestland, but the main preferential flow path is the root system of trees, which influences the connectivity of macropores; thus, the standard deviation of Pinus Yunnanensis plantation forestland is higher than that of grassland, and this finding is consistent with that of Alaoui et al. [37].For the same reason, the forestland macropores can efficiently transport vertically downwards compared to those of the grassland soil.This phenomenon is important in soil hydrological processes (e.g., the soil water cycle and water infiltration rate).The preferential flow pathway of secondary forestland is extremely developed because of the following: (1) the special morphology and high heterogeneity of the soil structure in karst regions [35], (2) the fact that the fracture between rock and soil has become a preferential pathway for transporting surface water to deep soil layers [51], and (3) the non-uniform distribution of fissures in surface soil, which is caused by the lack of uniform vegetation cover.The characteristic parameters of the preferential flow (e.g., ID (infiltration depth), UF, DC, PF, and L i ) can be calculated using the image data of stained vertical profiles, which can quantitatively evaluate the characteristics of dyeing morphology and the differences under various infiltration conditions and reveal the non-uniform characteristics of preferential flow under different infiltration conditions in karst regions (Table 3).The characteristic parameters under different vegetational types vary.Most indicators showed that G > PY > SF regardless of infiltration conditions, such as UF, DC and L i .However, exceptions were observed, such as ID and PF.The UF values of Pinus Yunnanensis plantation forestland, secondary forestland, and grassland were 17%, 3%, and 14% of ID, respectively, indicating that the infiltration depth can be affected by the types of vegetation on the surface when using the same infiltration amount.The C V and C µ of Pinus Yunnanensis plantation forestland were larger than those of secondary forestland and grassland.The surface coverage of grassland was nearly 100%; grassland can preserve moisture in an environment with a very large differences in the humidity values, such as in Yunnan, due to the shallow root systems.However, water infiltration in forestland, such as Pinus Yunnanensis plantation forestland and secondary forestland, is largely dependent on the root system and fractures, such as the crevices between rocks and soil, which can preferentially transport water to deep soil layers.Thus, preferential flow in forestlands is highly developed.The standard deviations of UF under G55 and ID under G15 and G55 are large, indicating that the process of water flow, including non-uniform flow, can confirm the presence of preferential flow [39].Relative to non-karst regions, ID is extremely deep.The reason may be due to the unique dual structure of karst landform (soil-rock) [51].Preferential flow becomes the major form of soil moisture transport in this area due to high heterogeneity and fissure structure.L i can also reflect the difference between preferential and matrix flows [52].Although PF and L i are parameters used for measuring preferential flow, the former emphasizes the dyeing region in the whole stained vertical section and the other focuses on different soil layers.
Figure 4 is based on the average DC of Pinus Yunnanensis plantation forestland, secondary forestland, and grassland after treatment with different infiltration amounts.It is evident that there is a downward trend in the average DC.Regardless of the infiltration changes amongst G5, G15, and G35, the regular infiltration of DC was similar and uniform, thereby indicating that the soil has a good water storage capacity.However, the fluctuation of forestland in storm conditions (G55) was conspicuous.The changing trend of grassland was similar to those in other infiltration treatments, but the fluctuation of Pinus Yunnanensis plantation forestland was less than that of secondary forestland, which meant that Pinus Yunnanensis plantation forestland can buffer vertical water infiltration better than the secondary forestland.Table 3 shows that the characteristic parameters of stained vertical sections increased with the amount of infiltration.The UF of G35 was similar to that of G55, which was almost twice those of both G5 and G15.The ID of G35 was similar to that of G15, but the ID of G55 was 1.60 times higher than that of G35, and those of G15 and G35 were 2.65 times higher than that of G5.The DC of G55 was significantly higher than those of the other infiltration conditions, which were 3.19, 1.86, and 1.54 times those of G5, G15, and G35, respectively.The PF under the four infiltration conditions was similar but still followed the rule of G5 (97%) < G15 (97%) < G35 (97%) < G55 (98%), and PY < SF < G, which meant that preferential flow can occupy a large proportion during water infiltration and occupy a small proportion in Pinus Yunnanensis plantation forestland incidents with less ponded water (G5).The L i of G55 was 2.40 times that of G5 and 1.30 times those of G15 and G35.The C V was similar under the four infiltration conditions, but that of G15 was significantly larger than the others.The C µ of G55 was 21 times those in the other three conditions, indicating that, among the extreme events, large water accumulation will affect the variability of the maximum water infiltration depth.The increase in infiltration amounts thus had a certain promoting effect on UF and ID.Moreover, the UF of G5, G15, G35 and G55 were 34%, 33%, 39% and 28% of ID, respectively, indicating that, in typical ponded water events, an increase in infiltration water has a certain promoting effect on the depth of matrix flow and its proportion.When the extreme events of large water accumulation occur, preferential flow dominates the process of water infiltration due to the special soil structure of karst landforms.This group of data demonstrated that the matrix flow under the different infiltration conditions showed no remarkable variation.L i increased with precipitation, and G55 was significantly larger than G35 and G15.G5 was significantly smaller than the three other infiltration conditions.Most parameters under G55 and its standard deviations were far greater than those in the three other infiltration conditions, indicating that water infiltration (e.g., preferential flow) will remarkably fluctuate during large ponded water events.Thus, we need to investigate the process of preferential flow under different conditions, its influencing factors and its impact on the environment.precipitation, and G55 was significantly larger than G35 and G15.G5 was significantly smaller than the three other infiltration conditions.Most parameters under G55 and its standard deviations were far greater than those in the three other infiltration conditions, indicating that water infiltration (e.g., preferential flow) will remarkably fluctuate during large ponded water events.Thus, we need to investigate the process of preferential flow under different conditions, its influencing factors and its impact on the environment.The changes in CV and Cμ can reveal the permeation mode of spatial variability in karst regions [53].CV does not decrease with the increase in the infiltration amount (Table 3).This finding conflicts with the results of Yao et al. ( 2018) [39].The reason may be due to the different spatial heterogeneities of soil in karst regions, which can strongly affect the depth of soil moisture and the development of preferential flow.As shown in Table 3, Cμ was negatively correlated with the amount of infiltration, and the water volume had a remarkable influence on the non-uniformity of the maximum infiltration depth.The Cμ of Pinus Yunnanensis plantation forestland was significantly greater than those in secondary forestland and natural grassland, indicating that the depth of infiltration had a strong spatial dependence in each forestland.This phenomenon illustrated substantial water translocation between the preferential flow path and the surrounding soil matrix.Transferring soil moisture to compacted soil or rock is more difficult than from the soil matrix to the surrounding environment [54].However, cracks caused by rocks or roots can enhance the development of preferential flow and increase the spatial variation in the morphological characteristics of the infiltration.Therefore, artificial land preparation and the rational allocation of surface vegetation can change the spatial heterogeneity of vertical preferential infiltration, and these results are consistent with those of Shinohara and Otsuki [55].The spatial variation (CV and Cμ) of Pinus Yunnanensis plantation forestland was greater than that of secondary forestland, and grassland was the smallest, indicating that artificial land preparation and rational allocation of surface vegetation can increase the variability of the water infiltration process.

Relationship between Preferential Flow and Infiltration
Different from loess plateaus [36] and wetlands [56], for which pollutant diffusion models have been established, karst areas require further research.Thus, we need to study the relationship between surface water infiltration caused by ponded water and preferential flow.Many latent variables cannot be directly observed during the infiltration process because of the imperceptibility of soil, and these latent variables cause errors.SEM has been shown to be verifiable [48].To accurately study preferential flow, we need to consider not only the correlation between variables but also the errors.Hence, we designed a relationship model between latent variables and estimated the fit of the The changes in C V and C µ can reveal the permeation mode of spatial variability in karst regions [53].C V does not decrease with the increase in the infiltration amount (Table 3).This finding conflicts with the results of Yao et al. ( 2018) [39].The reason may be due to the different spatial heterogeneities of soil in karst regions, which can strongly affect the depth of soil moisture and the development of preferential flow.As shown in Table 3, C µ was negatively correlated with the amount of infiltration, and the water volume had a remarkable influence on the non-uniformity of the maximum infiltration depth.The C µ of Pinus Yunnanensis plantation forestland was significantly greater than those in secondary forestland and natural grassland, indicating that the depth of infiltration had a strong spatial dependence in each forestland.This phenomenon illustrated substantial water translocation between the preferential flow path and the surrounding soil matrix.Transferring soil moisture to compacted soil or rock is more difficult than from the soil matrix to the surrounding environment [54].However, cracks caused by rocks or roots can enhance the development of preferential flow and increase the spatial variation in the morphological characteristics of the infiltration.Therefore, artificial land preparation and the rational allocation of surface vegetation can change the spatial heterogeneity of vertical preferential infiltration, and these results are consistent with those of Shinohara and Otsuki [55].The spatial variation (C V and C µ ) of Pinus Yunnanensis plantation forestland was greater than that of secondary forestland, and grassland was the smallest, indicating that artificial land preparation and rational allocation of surface vegetation can increase the variability of the water infiltration process.

Relationship between Preferential Flow and Infiltration
Different from loess plateaus [36] and wetlands [56], for which pollutant diffusion models have been established, karst areas require further research.Thus, we need to study the relationship between surface water infiltration caused by ponded water and preferential flow.Many latent variables cannot be directly observed during the infiltration process because of the imperceptibility of soil, and these latent variables cause errors.SEM has been shown to be verifiable [48].To accurately study preferential flow, we need to consider not only the correlation between variables but also the errors.Hence, we designed a relationship model between latent variables and estimated the fit of the initial model with the measured data.This study used SEM to construct and evaluate the index system, which can reflect the interrelationship between latent variables (exogenous and endogenous latent variables) [48].In this study, we used the software Amos 7.0 (Analysis of Moment Structures, developed by James L. Arbuckle.Amos 7.0 is a plugin for SPSS.) to construct the SEM equation.
The difference between the correlation coefficient calculated by SEM and the Pearson correlation coefficient is that Pearson will assume no limit due to measurement error, but SEM can provide not only an overall model test but also an independent parameter estimation test, thereby reducing the errors.Thus, the interference in correlation estimates between latent variables caused by measurement errors can be reduced.This method can explain the measurement errors.According to the multivariate statistics of James Stevens' social science application, at least 15 samples per factor are observed in the standard multiple regression analysis.We had 150 samples of vertical sections; therefore, this method was appropriate.
Initially, we measured the reliability of all the data with split-half reliability to display the Cronbach's alpha coefficient of each characteristic parameter.The Cronbach's alpha was calculated by SPSS (Statistical Product and Service Solutions, developed by IBM) as 0.72, which is greater than 0.70, indicating good data reliability.The path diagram of SEM was established after several revisions (Figure 5).Tables 4 and 5 show the path coefficients, errors, critical ratios, and significance of the influencing factor of the model between each characteristic parameter of preferential flow [46], Table 6 show the evaluation and fitting results of overall SEM fitness.We used maximum likelihood to perform simulation operations and corrected the model according to critical ration (CR); the path coefficient is significantly different from 0 at 95% confidence when CR > 2. The CR of each observed variable was greater than 2, indicating that we can use this model to study the correlation of preferential flow characteristic parameters.
In this study, the absolute fit index and the relative fit index were used to test the model (Table 5).Each of them conformed to the standard, indicating that the available of the model.Figure 5 and Tables 4 and 5 show that UF was significantly negatively correlated with PF (standardized estimate is −0.95), and DC was significantly positively correlated with L i (standardized estimate is 0.73).PF was calculated by UF (Equation ( 2)).Thus, the correlation coefficient between these parameters was the largest.The increment in UF also increased DC, and DC would promote L i , which indicated that the increase in matrix flow would not only inhibit the development of preferential flow but also increase the difference between the matrix and preferential flows.Ponded water was significantly positively correlated with ID (standardized estimate is 0.68) and DC (standardized estimate is 0.62), indicating that water infiltration caused by ponded water promoted increases in DC and ID, that is, the development degree of preferential flow.Our results are consistent with those of Yao et al. (2017) [39].When the pore structure of soil remains unchanged, DC will increase with the infiltration amount, and the increasing DC will increase L i .C V will decrease as L i increases because of the significant negative correlation between C V and L i .Thus, the increase in the amount of ponded water may decrease C V .Dyeing morphology is more standardized if L i is large, and the vertical migration of water is more significant than other forms of migration.Table 6 shows that the results in Tables 4 and 5 are credible, the results of SEM is feasible, because the statistical test quantity of the model conform the standards.Therefore, a large amount of infiltration will reach deep soil rather than remain in shallow soil.x 2 /df corresponds to the ratio of chi-square and freedom. 2 GFI corresponds to the goodness of fit index. 3RMSEA corresponds to the root meant square error of approximation. 4 NFI corresponds to the normative fit index. 5TLI corresponds to Tucker-Lewis index. 6CFI is the comparative fit index.

Relationship between Preferential Flow and Soil Nutrients
The occurrence and development of preferential flow in the soil are closely related to water, transport, and the accumulation of nutrients in the soil [57,58].In this study, we used SPSS to calculate the soil nutrients and characteristic parameters of preferential flow to study the correlation between preferential flow and the accumulation of soil nutrients (Table 7).Some of the measured nutrients showed large heterogeneity between fields [23], such as AK, AP and Org.However, the correlation between TP, TN, NN, and preferential flow is low, which may be due to the relative stability of TN, TP, and NN; the accumulation of these nutrients may be related to the biomass of surface vegetation [1].However, available nutrients are more easily absorbed by plants, and because of strong solubility of available nutrients [59], we should pay attention to the correlation between preferential flow and available nutrients.The correlations amongst AK and ID, DC and L i were positive, and the correlations between AK and C V and C µ were negative.The correlations between AP and UF and DC were positive, and the correlations between AP and C V and C µ were negative.The correlations between Org and PF and C V were positive, and the correlations between Org and UF was negative.Therefore, the preferential flow is closely related to the transport and distribution of important nutrients.The relationship between preferential flow and nutrient accumulation and distribution requires further study.This experiment cannot be processed by SEM due to the limitations in the sample size.Thus, the correlation coefficient calculated by SPSS is also reliable.
Table 7 indicates that the accumulation of available nutrients increases with DC and decreases with the increase of C V and C µ .The correlation between the distribution of AP and DC is evidently strong.This indicates that preferential flow can transport AP to deeper soil and is consistent with Julich et al. [60].By combining the results detailed in Table 7 and the analysis in Section 3.1, DC of grassland is higher than the Pinus Yunnanensis plantation forestland and secondary forestland, and C V and C µ are lower than Pinus Yunnanensis plantation forestland and secondary forestland, which indicates that the increase of preferential flow variability in grasslands may delay the accumulation of available nutrients.The accumulation of organic matters increases with PF and C V and decreases with the increase of UF.By combining the results detailed in Table 7 and the analysis in Section 3.1, UF of secondary forestland is lower than Pinus Yunnanensis plantation forestland and grassland, PF of secondary forestland is the highest, and C V of Pinus Yunnanensis plantation forestland is also the highest, which indicates that natural forestland contained more Org.However, C V is strongly correlated with available nutrients (negative correlation) and Org (positive correlation), and C V of Pinus Yunnanensis plantation forestland is higher than secondary forestland and grassland.Thus, the preferential flow of plantation forestland may promote the accumulation of organic matters.The accumulation of organic matter in karst should be studied because organic matters can improve the physical properties and softness of soil.Table 7 shows that the variability of preferential flow in the plantation can promote the accumulation of organic matters.
C V decreases with the increase of infiltration, according to the conclusions of Section 3.2.Thus, the increase in the infiltration leads to the loss of organic matter.Owing to the special geological conditions, soil formation speed is slow and the underground connectivity is high, thereby developing preferential flow.However, the risk that nutrients in the soil enter underground with preferential flow before being absorbed and utilized in the growth of surface plants is high due to the frequent rainfall in karst areas.Therefore, studies on soil erosion under karst landforms should be carried out, not only to include hydraulic erosion but also soil nutrient loss.
Plantation forestland is actively being planned to control rocky desertification in karst areas.The rational allocation of vegetation restoration should be considered for the sustainable development of planted forests.The soil can provide enough nutrients and water resources for afforestation by improving the accumulation of available nutrients and water conservation.The results showed that the accumulation of available nutrients in shallow soil could be promoted by increasing the planting of near-surface vegetation (herbaceous plants).It can also promote the accumulation of organic matter by reducing UF or increasing PF and C V .The contribution of preferential flow to groundwater and solute transport needs to be studied under more kinds of vegetation types because of the limitations in the limited sample size at the present.In future studies on karst land degradation, preferential flow needs to be considered and researched thoroughly. 1 AK corresponds to available potassium/ug•mL −1 . 2 NN corresponds to nitrate nitrogen/mg•kg −1 . 3TN corresponds to total nitrogen/g•kg −1 . 4AP corresponds to available phosphorus/mg•kg −1 . 5TP corresponds to total phosphorus/g•kg −1 . 6Org corresponds to organic matter/%. 7* correlation is significant at a 0.05 level (two-tailed). 8** correlation is significant at a 0.01 level (two-tailed).

Conclusions
Preferential flow can quickly transport surface water to deep soil.Because of the high surface temperature and evaporation in this area, preferential flow reduces surface water retention and increases available water for underground root growth.The preferential flow characteristics and accumulative distribution of nutrients in Pinus Yunnanensis plantation forestland, secondary forestland, and grassland fields were investigated using dye tracer experiments, which can reveal water infiltration patterns.This study reflects the response of preferential flow to the amount of infiltration (and some nutrients) under different water infiltration conditions.The precipitation amounts of nutrients can be significantly affected by the preferential flow in karst areas.Moreover, preferential flow can be effectively studied by SEM.The variation in preferential flow with depth results in spatial variability, which is influenced by the characteristic parameters of preferential flow.The matrix flow of grassland was superior to that of secondary forestland, whereas the preferential flow of the secondary forestland was better than other stand origins, and plantation forestland can effectively reduce the degree of preferential flow.Preferential and matrix flows have a large influence on the accumulation of available nutrients and organic matter under plantation forestland, secondary forestland, and natural grassland, or infiltration water volumes, indicating that preferential flow may affect the distribution of soil nutrients and the growth of surface vegetation by influencing the distribution of soil nutrients.Reasonable allocation of plantation forests has a certain mitigation effect on soil erosion in Karst areas, and preferential flow under this special geomorphological type is worth studying.Preferential flow can transport nutrients to deeper soil for roots according to the data of this study.Therefore, plantation is feasible under karst landform conditions, but it is better to combine herbaceous plants in a plantation.The results could provide suggestions for the restoration of rocky desertification and the advantages or disadvantages of vegetation restoration engineering in karst areas.

Figure 1 .
Figure 1.Example of stained profiles (cm) and image processing.

Figure 1 .
Figure 1.Example of stained profiles (cm) and image processing.

Figure 2 .
Figure 2. Example dye images for vertical soil sections in the three types of sites under different precipitation amounts.(a-d) correspond to the vertical sections of PY at G5, G15, G35, and G55, respectively; (e-h) refer to the vertical sections of SF at G5, G15, G35, and G55, respectively;(i-m) correspond to the vertical sections of G at G5, G15, G35, and G55, respectively.Among them, G5, G15, G35 and G55 correspond to 5 mm, 15 mm, 35 mm, and 55 mm ponded water infiltration simulations, which simulate the change in ponded water in light rain, medium rain, heavy rain, and storm, respectively.

Figure 2 .
Figure 2. Example dye images for vertical soil sections in the three types of sites under different precipitation amounts.(a-d) correspond to the vertical sections of PY at G5, G15, G35, and G55, respectively; (e-h) refer to the vertical sections of SF at G5, G15, G35, and G55, respectively;(i-m) correspond to the vertical sections of G at G5, G15, G35, and G55, respectively.Among them, G5, G15, G35 and G55 correspond to 5 mm, 15 mm, 35 mm, and 55 mm ponded water infiltration simulations, which simulate the change in ponded water in light rain, medium rain, heavy rain, and storm, respectively.

Figure 3 .
Figure 3. Average dye coverage (DC) (yellow curve) and standard deviation (black part) of three site types.PY corresponds to Pinus Yunnanensis, SF corresponds to secondary forestland, and G corresponds to natural grassland.

Figure 3 .
Figure 3. Average dye coverage (DC) (yellow curve) and standard deviation (black part) of three site types.PY corresponds to Pinus Yunnanensis, SF corresponds to secondary forestland, and G corresponds to natural grassland.

Figure 4 .
Figure 4. Average DC curve under different precipitation amounts ((a) is G5, (b) is G15, (c) is G35, (d) is G55).PY corresponds to Pinus Yunnanensis forestland, SF corresponds to secondary forestland, and G corresponds to natural grassland.Among them, G5, G15, G35, and G55 correspond to 5 mm, 15 mm, 35 mm, and 55 mm ponded water infiltration simulation, which simulate the change of ponded water in light rain, medium rain, heavy rain, and storm, respectively.

Figure 4 .
Figure 4. Average DC curve under different precipitation amounts ((a) is G5, (b) is G15, (c) is G35, (d) is G55).PY corresponds to Pinus Yunnanensis forestland, SF corresponds to secondary forestland, and G corresponds to natural grassland.Among them, G5, G15, G35, and G55 correspond to 5 mm, 15 mm, 35 mm, and 55 mm ponded water infiltration simulation, which simulate the change of ponded water in light rain, medium rain, heavy rain, and storm, respectively.

Figure 5 .
Figure 5. SEM of preferential flow characteristic parameters (e represents the errors).

Figure 5 .
Figure 5. SEM of preferential flow characteristic parameters (e represents the errors).

Table 1 .
Basic situation of the experimental sites.

Table 1 .
Basic situation of the experimental sites.

Table 2 .
Basic soil characteristics of the experimental sites.

Table 1 .
Basic situation of the experimental sites.

Table 2 .
Basic soil characteristics of the experimental sites.

Table 1 .
Basic situation of the experimental sites.

Table 2 .
Basic soil characteristics of the experimental sites.

Table 2 .
Basic soil characteristics of the experimental sites.

Table 3 .
Characteristic parameters of preferential flow variables.

Table 4 .
Regression results of SEM variables.
1S.E.corresponds to the smallest eigenvalue.2CRcorresponds to critical ration.3Pcorresponds to the probability of obtaining the critical ratio in absolute value.4***corresponds to the probability of obtaining the critical ratio in absolute value is less than 0.001.In other words, the regression weight is significantly different from zero at the 0.001 level (twotailed).

Table 5 .
Covariances results of SEM variables.
1Rf corresponds to the simulated amount of ponded water infiltration.

Table 4 .
Regression results of SEM variables.
1S.E.corresponds to the smallest eigenvalue.2CRcorresponds to critical ration.3Pcorresponds to the probability of obtaining the critical ratio in absolute value.4***corresponds to the probability of obtaining the critical ratio in absolute value is less than 0.001.In other words, the regression weight is significantly different from zero at the 0.001 level (two-tailed).

Table 5 .
Covariances results of SEM variables.

Table 6 .
Evaluation and fitting results of SEM overall fitness.

Table 7 .
Correlation analysis between cumulative content of nutrients in soil and the characteristic parameters of preferential flow.