Rupture and Migration Law of Disturbed Overburden during Slicing Mining of Steeply Dipping Thick Coal Seam

School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China Institute of Energy, Hefei Comprehensive National Science Center, Hefei 230031, China State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines, Anhui University of Science and Technology, Huainan 232001, China School of Civil and Hydraulic Engineering, Dalian University of Technology, Dalian 116024, China


Introduction
Coal seam mining causes overburden movements, which can jeopardize coal production safety. erefore, mining-caused overburden movement patterns are vital for mining technology [1]. However, due to the high cost of field research, limitations of monitoring equipment, and other factors, the overburden deformation and rupture laws are mainly determined by the physical analog modeling (PAM). is is a robust laboratory technique for reproducing, at small scale, geometric and kinematic features of natural geologic structures by a proper selection of deformation styles and analog materials. e latter include sand, claycake, silicon putty, and so forth and simulate rock strata table-top scale analog layers. Under the condition of satisfying the basic similarity principle, PAM can accurately reflect the engineering structure relationship. Because of its unique superiority, it is widely used in geotechnical engineering [2,3]. Many scholars have used PAM to analyze the gangue slip filling characteristics after rupture of steeply inclined coal seams overburden. ey studied the dip and antidip directions' pile structure in the process of steeply inclined coal seams workface gangue slip, obtained asymmetric "ruptured arch" of steeply inclined coal seams overburden, and revealed the asymmetric "arch shell" form of overburden space mining by steeply inclined coal seams [4,5,6,7,8,9]. However, due to the coal seam occurrence conditions and mining technology impacts, in steeply inclined and thick coal seams sliced mining, the overburden structure changes under the influence of multiple mining. Given this, the laws of disturbed overburden rupture and migration after overburden structural change in steeply inclined coal seams mining are very topical. e PAM approach envisages deformation monitoring via total stations, dial gauges, strain gauges, and so forth, whose testing accuracy and sensitivity need further improvement: most of them use point monitoring, while the layout of points and lines is too complicated. In general, deformation monitoring can only monitor the model's surface but fail to implement internal monitoring. With the rapid development of acoustic emission (AE) nondestructive detection technology, numerous researchers had successfully used it to detect the number of events and energy rates of the PAM when the chasm slips were unstable [10]. ey obtained the precursor information of chasm instability and revealed the characteristics and dynamic response laws of stope pressure before and after the activation of chasms [11]. Besides, the AE technology was also utilized in the assessment of characteristics of the spatial-temporal evolution of dynamic instability of large-inclination coal and rock masses [12,13,14,15,16]. Another hot field of research is fiber sensing technology, which can achieve long-distance distributed monitoring. Using the advantages of anticorrosion, antielectromagnetic interference, small size, lightweight, high sensitivity, high precision, and so forth, it ensures the stability monitoring for geotechnical [17,18], geology [19,20], water conservation [21], aerospace [22], defense [23], and other applications. References [24][25][26] incorporated the optical fiber sensing technology into the logical/ physical model to monitor the deformation and weighting of workface overburden, which yielded interesting findings.
is technology has been successfully applied to the monitoring of overburden migration in coal mine goaf [27,28].
In this study, the PAM, AE, and distributed fiber sensing (DBS) are jointly applied to steeply inclined coal seams for assessing the disturbed overburden secondary deformation and rupture patterns, as well as the characteristics of sound and light response in the process of rock deformation and rupture. is research can provide a theoretical reference for the prevention of hazards such as rib spalling, roof falling, and support instability. e research results are considered instrumental in the AE and DBS detection of rock and soil structure deformation and damage.

AE-Based Positioning by Time Differences.
Positioning by time differences is aimed to determine the sound source's coordinates or position by measuring the parameters of each AE channel signal TDOA, wave speed, and sensor spacing and their processing via a certain algorithm.
is accurate yet complicated positioning method has been widely applied in component detections. When performing planar positioning on a diamond array composed of four AE sensors, AE can obtain a real positioning sound source, as shown in Figure 1. In actual component AE monitoring, the number of AE sensors often exceeds four. If an individual AE sensor fails during monitoring, the remaining four AE sensors can still monitor a real positioning sound source.
Assume that the hyperbolic curve 1 is obtained from the time difference Δt X between S 1 and S 3 AE sensors, the hyperbolic curve 2 is obtained from the time difference Δt Y between S 2 and S 4 AE sensors, the positioning sound source is Q, the distance between S 1 and S 3 sensors is a, the distance between S 2 and S 4 sensors is b, and the wave velocity is V.
en, the positioning sound source is located at the intersection point Q (X, Y) of the two hyperbolae, and its coordinates can be expressed as follows: (1)

Optical Fiber Sensing
Measurements. DBS uses optical fibers as sensing elements. When the measured fiber is subjected to an external force or temperature variations, the frequency of the stimulated Brillouin scattered light in the fiber will change. e strain data transferred along the fiber can be obtained by conversion. e relationship between the Brillouin frequency shift, temperature, and strain is [29] Δε � where ΔV B is Brillouin frequency shift, MHz; C 1 is temperature sensitivity coefficient, MHz/°C; C 2 is strain sensitivity coefficient, MHz/με; ΔT is temperature variation,°C; and Δε is strain variation. is method is based on the principle of light scattering and measures the frequency change of Brillouin scattered light.
It is generally believed that when the temperature variation is within 5°C, the Brillouin frequency shift caused by the temperature can be ignored. In this study, this test is carried out indoors, and the measured temperature variation is less than 3°C. erefore, the temperature effect on the measurement results can be ignored.

Case Study.
e experiment uses 12123 working face of the Huainan Panbei Coal Mine in Anhui Province, China, as the engineering geological case study. e main coal seam is #13-1 coal seam, the average seam inclination angle is 40°, the average thickness is 4.4 m, and inclined mining in slices is adopted. e upper slice has a mining height of 2.0 m and a surface length of 120 m. Recovery was completed in 2015.
e lower slice has a mining height of 2.4 m and a surface length of 90 m. e recovery interval between the upper and lower slices is two years. e lower slice is arranged in the middle and lower parts of the goaf in the upper slice. Specifically, the lower slice's gate road is 5.0 m away from the gate road in the upper slice. Figure 2 is a schematic diagram of the location and production conditions of 12123 working face. e coal seam hardness coefficient (f ) is 0.3∼0.5, the immediate roof is the interslice of mudstone and sandy mudstone, and the average thickness is 9.3 m. e immediate floor is mudstone, and the average thickness is 3.0 m.

Physical Analog Model (PAM).
e test uses a rotating plane PAM frame; the device dimensions are length × width × height � 2.0 m × 0.2 m × 2.0 m. PAM requires that each part's size is enlarged or reduced according to a uniform ratio between the physical model and the prototype, which is l m /l p � C l . According to the 12123 working face engineering geological data and test device dimensions, we determined the geometric similarity constant l m /l p � C l � 100, bulk density similarity constant C c � c p /c m � 1.6, stress similarity constant C c � C c. C l � 160, and time similarity constant C τ � C l 1/2 � 10. In the above similar constant formulas, the subscript parameters p and m represent the prototype and model, respectively. e test uses river sand as aggregate, gypsum, lime as cementing material, and water as the bonding material. Due to the article's limited space, the detailed similar material ratio and modeling process are omitted and can be found elsewhere [30,31]. Based on the similarity constant and the physical and mechanical parameters of coal, the weighted mixing test between the above analog materials was used to adjust and determine the similarity ratio and then calculate the proportions of analog materials in each coal slice. Analog materials are also stirred, laid, and compacted slice by slice. Different colors are sprayed on different rock slices of the high-level floor of the workface to distinguish them. e physical and mechanical parameters of coal and similar materials are listed in Table 1.
e test laying height is 1.5 m, and the workface burial depth ranges from −375 to −495 m. After the overlying nonsimulated rock slice is converted by a load similar constant, the top of the physical model is loaded by a multifunctional adjustable PAM loading device. e complete physical model is depicted in Figure 3.

Monitoring System.
Coal pillars of 50 cm are left on each side of the upper slice. e upper and lower simulation slices' excavation distance is 10 cm, the upper slice is excavated 12 times, and the lower one is excavated nine times; each excavation duration is about 1 min. e interval of each excavation is 20 min.
Steeply inclined coal seams are generally cut from top to bottom, and the frame is moved from bottom to top. When the rock slice is broken, the slippage mainly occurs during the frameshift phase. erefore, the bottom-up simulation of coal cutting is adopted in the two-dimensional physical model. Six AE sensors, numbered P 1 -P 6 , are embedded beforehand in the physical model, to facilitate the picking of 3D positioning points. Among them, P 1 , P 3 , and P 5 are embedded in the front, while P 2 , P 4 , and P 6 are implanted in the back. e specific locations are shown in Figure 4. Before the workface simulation recovery, the AE parameters need to be calibrated, and the upper slice recovery process is used to check the accuracy of the AE parameters. e AE sensor's resonance frequency is finally determined to be 100-600 kHz, the sampling rate is 3 MHz, the preamplifier gain is 40 dB, and the monitoring threshold is 90 dB. With such AE parameters, the AE signal is not affected by the excavation action, which is convenient for later data screening. When the lower slice is excavated, the AE is turned on synchronously until the overburden structure is stable.
e distributed fiber uses a high-transmission low-rigidity responsive sensing optical cable with a 2 mm diameter, and the fiber placement is shown in Figure 4. e optical fiber is embedded in analog materials. To make the optical fiber tight, a certain pretightening force is applied to both optical fiber ends. A small amount of cement slurry is put to make it tightly coupled with analog materials. Like AE parameter calibration, optical fiber also needs to be calibrated for strain coefficient accuracy. For brevity sake, a detailed calibration process can be found elsewhere [32]. e spatial resolution of the optical fiber strain measurement is set to 50 mm. e optical fiber can monitor the rock slice's deformation and rupture before and after the overburden structure disturbance in the upper and lower slices on a real-time scale. e optical fiber strain's positive and negative values indicate that the analog material is in a stretched and compressed states, respectively. e physical model is installed in front of the camera to monitor the rock deformation and rupture.

Rock Slice Rupture and Migration Features.
After the upper slice excavation, the immediate roof breaks at the height of 9.3 m, and the gangue moves down along the floor of the workface. e gangue's nonuniform filling characteristics are enhanced, showing the goaf filling with an "upper empty, middle full, and down solid" feature. e immediate high-level roof has a macrofracture (crack line in A-A red frame), and the rock slice below the fracture surface loses its bearing capacity. In contrast, a stable cantilever beam structure (A-A red frame) is formed above the fracture surface, which has a certain bearing capacity. In the workface Advances in Civil Engineering lower area, the gangue piles up, the overburden's movement space in the mined-out area decreased, and overburden shows transverse and longitudinal cracks and delamination [33,34]. It is still not completely separated from the stable rock slice and has the form of a hinged rock beam structure.
Its load directly acts on gangue (or bracket). According to the time similarity constant, after the upper slice overburden is stabilized, continue to excavate the lower slice, as shown in Figure 5. After the lower slice's excavation, there is a free space height below the disturbed overburden, and the  gangue's ability to restrain the disturbed overburden is weakened. Under the superposition of overburden gravity and mining stress, the hinged rock beam formed by disturbed overburden breaks and slips after "activation." e disturbed overburden structure is changed.
Specifically, the overlying low-level disturbed overburden (B-B middle and lower parts sprayed with red and dark green) emerges during mining, which exposes the high-level disturbed overburden above it. Gangue slides along the workface floor to fill the lower part of the goaf, while the free space below the disturbed overburden in the middle and upper parts of the goaf (available for the rock slice to break and rotate) increases. e disturbed overburden breaks up to the main roof (B′-B′ part sprayed with dark blue paint). e breaking height is 16.5 m. After the simulated excavation reaches the upper end of the lower slice, the structure of the cantilever beam formed after the upper slice excavation is unstable and together with the disturbed overburden above the coal pillar has slip damage without any support (A′-A′ red frame), with the development of coal pillar cracks and severe gangling (A′-A′ yellow frame).
Comprehensive analysis shows that, during the lower slice recovery process, the disturbed overburden is severely broken under the secondary mining action, and it is very easy to induce the roof collapse, which causes the bracket to be improperly connected, and the bracket falls, slips, reverses, and becomes unstable. For the lower slice return airway, the first required measure is to prevent the roof from breaking, and the second one is to avoid the breaking impact of the cantilever beam formed by the high immediate roof above the coal pillar on the support after the upper slice recovery.

Rock Slice Rupture AE Response.
e physical analog model has a large scale, and the distance between the AE sensors is relatively large. According to a single AE sensor, it is difficult to obtain the AE response law of distorted overburden deformation and rupture. Hence, based on the time node after disturbed overburden structure mutation of lower slice mining, the data of six AE sensors are synthesized to extract the parameters, such as energy and ring count, and are recorded as the number of one event. e 3D positioning (AE events) in the gangue sliding process is relatively random. Only the 3D positioning AE data of the disturbance caused by high roof breaking, rupture of the disturbed main roof, and sliding of disturbed overburden are extracted. Energy is a relative value monitored by the AE analyzer per unit time and is an important indicator of the fracture rate and size of the rock mass; the ring count is the cumulative value monitored by the AE analyzer per unit time and is an important indicator of the rock mass destruction. Figure 6 depicts the relationship between the break and slip of the lower slice disturbed overburden and time. During the lower slice recovery process, primary cracks in the disturbed overburden continue to develop, a large amount of elastic energy accumulating at the crack tip and its surroundings. e AE energy gives rise to six higher-energy bundles. e first energy peak of the AE multimodal energy bundle is the largest. en, the peak value gradually decreases within the duration of each disturbed overburden break, or it appears as a single peak. It shows that the instability of the lower slice disturbed overburden has gone through six stages, and when the disturbed overburden breaks in each stage, the first break strength is the largest. e first energy bundle appears when the disturbed low-level immediate roof (B-B part painted in red in Figure 5) is broken. e elastic wave released by the broken overburden at this stage is low, and its duration is shorter. e peak value of AE energy is 778.3 mV·mS, and the cumulative ring count increases slowly. e second energy bundle appears when the disturbed high-level immediate roof (B-B part painted in dark green part Figure 5) is broken. At this stage, the peak value of AE Advances in Civil Engineering energy shows a "group shock" phenomenon. e highest peak value reaches 4973.1 mV·mS, the duration is 38 s, the cumulative ring count jumps from 2867 to 3631, and the total number of events in this stage is 2462 n. e third energy bundle appears when the gangue slides along the workface floor after the disturbed high-level immediate roof is broken. e lower slice's gangue is stacked, the free space height in the middle and upper regions increases, and the exposed length of the disturbed main roof increases. e fourth energy bundle appears when the exposed disturbed main roof (BB sprayed blue part in Figure 5) breaks. At this stage, the peak value of AE energy is "isolated." e isolated peak value of AE is 3914.3 mV·mS. e duration is 3s, the cumulative ring count jumps rapidly from 4541 to 5574, and the total number of events at this stage is 726 n. is shows that the rock breakage at this stage is "abrupt." After the rock slices at the first and second stages are broken, gangue has filled the lower and middle parts of the lower slice. e space available for gangue slippage after the fourth stage rupture of the disturbed main roof is reduced; that is, the fifth energy bundle caused by gangue slippage is weak and exhibits a "solitary shock" form. When the upward simulation excavation reaches the upper part of the lower slice, the disturbed overburden above the coal pillar breaks and slips under the joint action of overburden load and mining stress. e sixth energy bundle emerges. e highest energy peak is 3801.5 mV·mS, the duration is 11 s, the   Advances in Civil Engineering cumulative ring count jumps from 5854 to 7270, and the total number of events is 2084 n. e procedure requires breaking the disturbed high-level immediate roof, importing the 3D positioning AE data of the rupture of the disturbed main roof, and sliding the disturbed overburden into AutoCAD software to generate 3D positioning line drawings and perform postprocessing such as coloring. Figure 7 shows the spatial positioning pattern of the generated disturbed overburden breaking AE signals. As observed, the distribution of AE spatial positioning points is roughly corresponding to the broken position of the disturbed overburden (A′-A′, B′-B′ red box in Figure 5) in the lower slice. Disturbed high-level immediate roof breaking off locations of AE in the space are more scattered, mainly extending from the location of the disturbed high-level immediate roof to the surrounding area, showing the relative ductility destruction of the disturbed high-level immediate roof rock slice. Rupture of disturbed main roof AE positioning points is concentrated in space, showing the relative brittleness damage of disturbed main roof. After simulating the excavation of the lower slice's upper area, disturbed overburden damage, and slips, the AE space positioning points are distributed in a "flow-like" way. e comprehensive analysis revealed that macro-and microcracks in the disturbed high-level roof were more developed. Under the superimposed action of overburden load and mining stress, macro-and microcracks continued to store elastic energy. When the energy storage limit was reached, the disturbed high-level immediate roof breaking lasted for a long time, and the energy released was quite high (the second energy bundle). However, the amount of energy released per unit time was relatively small. e internal cracks in the rock slice had sufficient time to evolve and develop. e distribution of 3D positioning points was divergent, and the damage degree of the rock body was increased. e rock body was broken, and the broken gangue moved violently in the lower slice goaf (or support top beam). Disturbed main roof macro-and microcracks were not developed. Under the superimposed action of overburden load and mining stress, the macro-and microcracks released higher energy after reaching the energy storage limit in a short time (the fourth energy bundle). e cumulative ringing count per unit time exhibited a sharp jump, and the rock mass was quickly fractured. e internal structure of the rock slice had no time to evolve and was damaged. e distribution of 3D positioning points was aggregated, and the overall damage degree of the rock body was small. e disturbed main roof possessed the characteristics of sudden damage and could produce the impact-caused instability to the lower sliced support. When the simulated excavation reached the lower slice's upper part, the disturbed overburden slips, but its integrity after fracture was high. e distribution of 3D positioning points was streamlined. e cumulative value of the ringing count in the disturbed overburden slipping process per unit time was between the value of disturbed high-level immediate roof process and the value of rupture of disturbed main roof process.

Rock Rupture Fiber Response.
e first mining of the upper slice (advance of 10 m) is detected by optical fiber sensing, indicating that the fiber senses the weak deformation inside the rock slice. As shown in Figure 8(a), when the upper slice excavation is simulated, the fiber response shows a "double peak" shape, and peak sizes are 436 and 534 με, respectively. e corresponding physical model heights are 680 and 980 mm, respectively. e fiber response double peak position is projected onto the upper slice workface, which is mainly located in the middle and upper regions of the upper slice. After the steeply inclined coal seams mining, the overburden in the upper part of the workface breaks and slips. e gangue fills the lower part of the workface, hindering the continued breakage of the overburden in the lower part of the workface. ere is a high Advances in Civil Engineering peak fiber response in the upper part of the workface and a small peak response in the lower part, further explaining the asymmetric feature of steeply inclined coal seams workface overburden breakage.
Taking the fiber response before and after the disturbed high-level immediate roof breakage as an example, the optical fiber response change characteristics before and after the lower slice disturbed overburden breakage are analyzed, and the results are shown in Figure 8(b). As observed, before the lower slice disturbed high-level immediate roof breaks, the tensile stress concentration occurs in the disturbed overburden (the fiber response value is positive), and the elastic deformation in the rock gradually accumulates. e peak value of the optical fiber response is 2095 με, and the peak position is located above the breaking line of the disturbed high-level immediate roof, while the height of the corresponding physical model is 673 mm. After an immediate high-level roof of the disturbed lower slice is broken, the rock mass in the fiber-implanted slice shows a certain degree of stress release, and the accumulated elastic strain energy is also released. e optical fiber response's peak value drops from 2095 to 1717 με, and the peak position is the same as before the disturbed high-level immediate roof breakage. is is consistent with AE's response law before and after the breakage of the disturbed high-level immediate roof (the AE energy accumulation in the rock before the breakage, and the release of AE energy after the breakage). For the rupture of the disturbed main roof and disturbed overburden, the fiber response before and after the slip also follows the same pattern (more detailed analysis is omitted for brevity). e critical point on the optical fiber response curve in the whole process of the lower slice simulation excavation from monotonous increase to monotonous decrease is roof breaking. e instantaneous response before breaking reaches the maximum value. After breaking, the roof rock slice quickly releases pressure, and stress is transferred to the deep stable rock mass around the goaf. At this time, the roof rock slice's stress at the breaking position reaches the minimum value, and the overburden movement gradually stabilizes. After that, as the work surface advances, it monotonically increases again, and the process forms the work surface with a periodic mineral pressure trend. e cycling process forms the appearance of the workface's periodic mineral pressure along the inclined workface.
During the lower slice simulation excavation, the optical fiber response data after the disturbed high-level immediate roof rupture, rupture of the disturbed main roof, and disturbed overburden slipping are extracted, as shown in Figure 8(c). As observed, after the rupture of the disturbed high-level immediate roof, the rock slice is   Advances in Civil Engineering opened to form a separation crack, which causes the tensile stress concentration in the optical fiber. With the upward excavation of the PAM and the rupture of the disturbed main roof, the optical fiber response's peak value jumps from 1717 to 2547 με in steps. e peak position is located above the line of the rupture of the disturbed main roof. e specific height of the corresponding PAM is 796 mm. With the simulated excavation of the lower slice's upper area, the optical fiber senses the disturbed overburden slip in the upper area of the lower slice, and the peak value of the optical fiber strain also shows a stepped jump of 1575 με. e peak position corresponds to the PAM height of 1078 mm. e distributed fibers monitor the process of abscissa slice development, expansion, closure, rebirth, and so forth. e abscissa slice development length gradually expands, as the lower slice progresses upward along the inclination. e sandstone in the upper part of the lower slice (A′-A′ and B′-B′ sprayed green part in Figure 5) does not change significantly in the subsequent excavation process, indicating that the upper part sandstone constitutes the main key slice of the physical model.

Substantiation of ree Key Protection Areas in the Lower Slice by Fiber Sensing and AE Methods.
e AE and distributed fiber monitoring results inevitably have errors. PAM has a plane stress state from an experimental perspective, and constraints do not effectively control its surface. e brittle sand material affects the contact area of the sensor and the effective transmission of strain data. From the perspective of AE and fiber sensing technology, the broken rock slice leads to incomplete contact between the sensor and rock slice, causing slippage and shedding. All these factors may bring uncertainty to the measurement. Although different sound and light measurement methods have certain differences in specific values, the overall regularity is not changed. It is necessary to analyze the influencing factors in further research and eliminate them gradually.
Under the influence of secondary mining, the lower slice disturbed overburden slips, and the lower slice's return airway is seriously damaged. e disturbed main roof in the central area of the workface breaks "suddenly," which has a positive impact on the bracket and may cause the bracket to fall, slip, and become unstable. e disturbed high-level immediate roof in the lower area of the workface breaks, the gangue slides down along the floor of the workface, and the workface undergoes such disasters such as roof collapse. To control the mining system's stability in different regions, the so-called "triunity collaborative support system" is proposed for the protection of roadway sides, workface roof, and support brackets. e above support system implementation in the three key protection areas of the lower slice will make it possible to avoid the workface's large-scale roof collapse, side falling, bracket instability, or other violation of the workface mining safety. Engineering practice has verified the feasibility of the proposed subdivision of the three key protection areas in the lower slice.

Conclusions
e physical analog modeling of the particular case study (the Huainan Panbei Coal Mine in Anhui Province, China) was used to derive the rupture characteristics and migration of the disturbed overburden in steeply inclined coal seams. e joint application of AE and fiber response techniques made it possible to draw the following conclusions: (1) e secondary mining results in breakage and slip of the disturbed overburden hinged rock beam in the lower area and the cantilever beam in the upper area. e overburden structure is mutated, which causes the lower slice roof to break, and it is easy to induce the workface roof collapse, side falling, instability, and so forth.
(2) e lower slice's disturbed overburden experiences six stages of damage and instability, and at each stage of the broken overburden, the AE energy bundle is multimodal. e first energy peak in the AE multimodal energy bundle is the largest. (3) e optical fiber can sense a small-scale deformation of the rock mass. e critical point where the optical fiber response develops from monotonous increasing to decreasing is roof breaking. e instantaneous strain before breaking attains its maximum value, and the overburden reaches stability when the stress is released to a minimum value. After that, as the workface advances, it is again monotonically increasing. e cycling process creates periodic pressure fluctuations along the inclined workface. (4) During the lower slice's mining process, three key prevention and control areas, namely, the disturbed high-level immediate roof rupture, rupture of the disturbed main roof, and sliding of disturbed overburden, are observed. e AE energy per unit time in the three phases is higher than that in other breaking phases, and the rock mass fracture rate is higher, while the fiber sensing exhibits a sudden change in response value and a stepped jump increase. e research results are considered instrumental in the application of acoustic emission and fiber sensing to the deformation and damage monitoring of rock and soil structures.
Data Availability e data used for conducting classifications are available from the corresponding author upon request.

Conflicts of Interest
e authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.