Dessication cracking may be perceived as a two-fold process that entails the gradual loss of moisture from an initially saturated porous body, in this case a soil and after surpassing a threshold tensile strength, culminates into germination, propagation and widening of cracks with subsequent reduction in soil volume. [1] [2] describe this as a three-stage process: initial, primary and steady states. Several intrinsic factors like: particle size distribution, degree of saturation, presence of occluded void spaces, tensile and inter-aggregate cohesion forces are attributable to initiation of cracking in soils. Despite much research work on soil cracking phenomenon over the last decade, little has been reported on accurately predicting and localizing crack offset points as well as on geometry and direction of propagation. Other studies on crack growth and propagation have applied stochastic models [3] or physically based model [4] - [6] However, much of these investigations [7] - [9] have been qualitative and qualitative in nature and largely restricted to describing desiccation cracking as a phenomenon and not a dynamic process. Soil cracking therefore, has been described in terms of its functional properties: enhancing preferential flow for soil nutrients as well as pollutants in clay liners [10] or as in controlling both slope stability and those of built structures.

The complexity of dessication cracking in soils can better be understood through the Fractal Dimension. This approach seeks to characterize complex patterns at the Representative Elementary Volume (REV) and quantifying them as ratios of change from the micro-to mega-scale. It is based on the premise that dessication cracks do not exhibit typically Euclidean or topological dimensions in space and so fractal indices may only be expressed as non-integer values. Cracks are more or less curvi-linear and would necessarily have a dimension greater than 1 and less than a 2-dimensional square and this is the in-between Hausdorff-Besicovitch dimension.

The paper presents an approach based on the fractal dimension fd, that describes the spatio-temporal and dynamic change of crack development in a Black Cotton clay soil (Eutric Verticsol) extruded from along the banks of River Nile in Juba, South Sudan. The objective of the paper was to describe the cracking phenomenon during dessication in a Eutric Vertisol sample using: Crack Intensity Factor (CIF), Crack Area (CA), density functions f(l) and f(W) as indicators for both crack length and width respectively. The results of this study are expected to deepen our understanding on crack development as a dynamic process.

Visual observation (Figure 1) showed that crack patterns generally had a “T” shape at intersection nodes with secondary cracks perpendicular to primary propagating cracks, especially during the first crack segments and early stages of dessication. This is attributable to growth of cracks in the direction perpendicular to the local maximum tensile stress [19] - [21] . Although the crack pattern at the intersection nodes gradually assumed a more “Y” shape in the second crack segments, such crack patterns did not necessarily follow a strict “T” and “Y” sequence during crack propagation.

Thinned and skeletonized images of cracks during dessication are shown in Figure 2. In day 1 most thinned images of cracks are shown as simple discontinuous that gradually merged with each subsequent day of dessication. The contrast between the black binary image with its immediate background for example in days 4, 6, 8 and 10 were a result of poor quality of the original photo images captured that could not be adequately thinned and skeletonized.

From Table 1, it was found out that the fd prior and after thinng procedure positively correleted with the Cover or Bright/Difference and best expressed using an exponential function: y = 0.0841e^0.0029x at r² = 0.69 suggesting that increasing fd was directly associated with increment of black coverage and therefore of crack length in the digitized image during dessication. This was also corroborated with the increasing value of CIF especially between day 1 to day 8.

The log of crack area versus SMC (q) during drying of chosen soil specimen as a function time is shown in Figure 3. The results indicate that log CA increased quickly with decreasing SMC during the early stages of cracking. The CA of 967.5 cm² on day 1increased at an estimated rate of 2 cm²/day till day 5 and decreasing at an increasing rate of 1 cm²/day to CA 53.334 cm² on day 6. However, upon further drying the rate of CA increase started to decline and gradually reached stabilization on day 7. This observation is consistent with earlier

studies on desiccation cracking behavior of polypropylene fiber-reinforced clayey soil [19] [22] . This stabilization in crack area was approximated as the intersection between the log of crack area and soil moisture content lines. However, visual inspection of cracking during dessication showed that stabilization was attained on day 4 and no new cracks developed, rather the widening of pre-existing cracks as from day 8 till end of experiment on day 32. It can be mentioned that, crack initiation and propagation stopped at about SMC 37% enhancing only the widening of pre-existing cracks.

Dessication time versus CIF or Cover/Brightness for the chosen soil specimen is presented in Figure 4. It showed that, the CIF developed gradually between 0% to 1% during the first 4 days and experienced a rapid growth between day 4 to 8 and attained peak value of about 6% and about 11% (Cover/Brightness) during day 8. This accounted for about 98.6% of all cracks that developed between 45% SMC at the onset of drying to about 35% 8 days later. After day 8, there was a rapid CIF decline to about 2% and Cover/Brightness (about 5%) at SMC 28% and gradually decomposed over the entire experimental period to a value between 1% and 2% for the

CIF and between 3% and 4% for Cover/Brightness respectively. Our study also showed that, estimating the CIF from thinned binary images through the pixel count as well as the Cover/Brightness generated by the software were adequate and showed similar decomposition over time, although with an order of magnitude of 1 - 2 between them. The SMCs that enhanced initiation, rapid propagation and development of cracks during the first 1 -

4 days are herein referred to as critical soil moisture contents analogous to the critical suctions employed by [22] . As aforementioned, the Crack Intensity Factor (CIF) of a digitized binary image was determined and expressed as the ratio of the surface area of the cracks to the total surface area of the binary image of the soil specimen.

Variations of the SMC with the fd are shown in Figure 5. Low fd values between1.2 to 1.3 were observed at high critical SMC value of about 45% during day1 on the onset of cracking. At day 2, the SMC at about 42%, gradually decreased with concurrent increase of the fd to between1.4 to 1.5. With each proceeding day at day 4 to 8, the SMC decreased to about 35% with fd increase to about 1.7 as dessication cracking proceeded and the number of cracks increased [23] until a peak fd value between day 8 to 10 was attained. With further decrease in the SMC during day 10 to day 32 from 28% to16%, the fd correspondingly decreased from about 1.7 to 1.3. Fractal dimension can be treated as a measure of soil cracking during drying of the soil sample and is therefore not invariable. It increased to a peak value once maximum threshold value of cracking was attained and subsequently decreased. Our results showed that fd had both spatial and temporal variability implying that the fractal

behavior of cracks at time (t₁) was constant over a given crack length which corresponded to a fractal

dimension (fd₁) and changed at time (t₂) over yet another crack length that corresponded to a fractal

dimension (fd₂) and so on. Lower fd values closer to 1.0 during the first days of dessication cracking would in this case suggest crack pattern homogeneity with gradual transition to a more heterogenous crack pattern with fd values closer to 2.0. The fd downturn from 1.7 to 1.3 at day 10 to 32 suggested crack pattern homogenization as the soil dried up with widening of pre-existing cracks. Obviously, the low quality digitized crack images that were captured at day 10 to 32 showed more grey-black shades and without thresholding (thresholding is a procedure that transforms grey or color image objects into black and white), where further filter and skeletonization procedure tended to generate poor images. Continued filtering and skeletonization only led to more detached and isolated black spots along the crack length (Fournier dust) which reduced the fd value. Our results confirm that low fd value at the onset of dessication cracking characterize a more homogenous, simple curvi-linear crack patterns that with increasing fd value toward a more heterogenous and differentiated crack patterns.

According to the fractal theory, the relationship of log N(s) is linearly related to log (1/s) so that the correlation coefficient denoted as the slope is the fractal dimension fd. As illustrated in Figure 6, the fd before and after thinning of the binary image in the thirty days of dessication cracking was shown with fd values prior to thinning slightly higher than those after, although with no significant differences. Slight fd differences were evident at lower log (1/s) values which then merged at higher log (1/s) values. This “bifractal” behavior may be attributable to the presence of multiple black spots that necessitated a higher pixel count and so higher fd values prior to thinning. Both fd values prior and after thinning tended to merge at higher log (1/s). This could be as a result of widening of pre-existing cracks that yielded higher fd values and so were equal in magnitude to those fd values prior to thinning. All correlation coefficients estimated by the Fractal3e software and based on the fractal theory

are represented in Table 2. The results showed that fd values were between 1 and 2, i.e. (2 ³ fd ³ 1) and that the dessication cracks were of fractal nature.

Figure 7 generally showed a positive relationship between the fd and the CIF during dessication with the CIF varying between 0.0% to 6.0% whereas the fd between 1.2 to 1.7. Interestingly, the CIF increased during dessication in the days 1 to 8 attaining maximum value of 6% and fd 1.5 while it gradually reduced in the subsequent days 10 to 32 to as low as 1.2%.

Figure 8(a) showed that increase in crack length was time-dependent and attained a peak value at about f(l) = 0.015with a decrease of f(l) below 0.02 during the subsequent days of dessication. Such increase and decrease in crack length during dessication showed a log normal distribution pattern with parameters s = 91.9 and m = 80.59. Although both parameters had no physical significance in soils and could not empirically be corroborated in terms of their magnitudes, they seemed to significantly influence both waterflow and solute transport [24] . This unimodal distribution of f(l) suggested the concurrent increase in crack length in both primary and secondary cracks. It appears that any further dessication after day 8 had no significant effect on crack length increase rather than the widening of pre-exisiting cracks.

The crack number for the different crack widths (Figure 8(b)) also showed positive correlation with fd. The crack number was highest between 0.4 - 1.6 for crack width between 0 - 10 mm and fd between 0.2 - 1.6 as expressed by the polynomial function (0 - 10 mm = −0.0139x² + 0.6151x + 0.8921; r² = 0.89). However, low fd £ 1 values are unusual and would suggest the inadequacy of associating fd with crack number when characterizing crack geometry. Cracks with average width between 10 - 20 (10 - 20 mm = −0.0082x² + 0.4137x − 0.9383; r² = 0.88) and 20 - 30 mm (20 - 3 0 mm= −0.0001x² + 0.0135x − 0.0348; r² = 0.90) showed low crack numbers between 0.2 - 0.4 as well as fd values between 0.2 - 0.4. The cut-off values fd £ 1 also showed that the crack number as a parameter was inadequate in describing crack geometry since the relative numbers of widened cracks 10 - 20 and 20 - 30 mm was low.

The density functions of crack width f(W) during dessication cracking of a remolded vertisol are shown in Figure 9. The maximum value f(W) of a crack width may be perceived as most probable value (MVP) or the probability distribution of f(W) that is close to maximum [19] . The cracks with different widths showed typical bimodal distribution with maximum f(W) values at 6E-09 and 2E-09 for 0 - 10 mm; f(W) = 0.012 for 10 - 20 mm; and f(W) = 0.016 and 0.012 for 20 - 30 mm. This bimodality is a lump up of two sets of cracks; the main primary cracks and the secondary sub-cracks that developed from the former. The results indicated that the MVP and so the f(W) as well as the fd increased with dessication time to day 8 - 10 especially within the primary cracks, whereas this was between day 28 - 30 in the secondary cracks. The maximum fd = 1.3 for both crack widths 0 - 10 and 10 - 20 mm respectively while this was fd about 1.44 for 20 - 30 mm crack width suggesting a positive correlation between crack width and fd. This bimodality also suggested that maximum increase in crack widths for all tested crack sizes during dessication occurred in two phases: at the start (approximately at day 8 - 10) for predominantly primary cracks and at the end of dessication cracking (day 30) for the secondary cracks.

We analyzed and characterized a series of cracking processes in time as a function of water loss. The cracks that developed during dessication of a remolded Eutric Vertisol clearly showed fractal properties between 2 ³ fd ³ 1 During the initial and the subsequent stages, the crack length and number increased and showed that, the fd was not a constant but rather a spatio-temporal variable characterizing a type of cracking behavior that was conditioned by soil water loss during drying. From the conducted experiment, the following conclusions could be drawn:

1) The fd positively correlated with both CIF and time during dessication.

2) Crack length expressed as log of Crack Area CA, increased with decrease of SMC till a cracking stability phase was attained with no further crack length increase except the expansion or widening of pre-existing cracks.

3) Crack width increased in two phases whereof increase in the first 1 to 10 days during dessication was predominantly within the primary cracks, and between 20 - 30 days was within the secondary cracks.

The authors are very grateful to the Norwegian-funded project NUCOOP (Norwegian Universities Cooperation Project) under the title “Post-war Livelihood and Environmental Studies” Project No. 2000/10003 hosted at College of Natural Resources and Environmental Studies (CNRES), University of Juba for funding the purchase of both penetrologger and Theta moisture sensor. Similarly, the authors would like to thank the Japanese National Agriculture and Food Research Organization (NARO), for availing the Fractal3e software.

David Lomeling,Mandlena C. Kenyi,Modi A. Lodiong,Moti S. Kenyi,George M. Silvestro,Juma L. L. Yieb, (2016) Characterizing Dessication Cracking of a Remolded Clay (Eutric Vertisol) Using the Fractal Dimension Approach. Open Journal of Soil Science,06,68-80. doi: 10.4236/ojss.2016.64008