Understanding the Mechanisms of Earth Fissuring for Hazard Mitigation in Najran, Saudi Arabia

: Being a fast-growing city with a high rate of urbanization and agricultural development, the city of Najran, situated in the southwest of the Kingdom of Saudi Arabia, has witnessed a series of earth ﬁssuring events and some other geo-environmental hazards in recent times. These ﬁssures have posed a signiﬁcant threat to inhabitants and infrastructure in the area. A few studies suggest that excessive groundwater withdrawal is responsible for ﬁssuring activities. Because of the intensity of this geo-hazard, this article presupposes that groundwater extraction alone cannot be responsible for the magnitude of ﬁssuring activity in the area and discusses other severe factors that could be responsible for the earth ﬁssures. The study proposes that the cause of the problem is multifaceted and synergistic, and outlines threatening factors that can inherently trigger more ﬁssures in the region, based on the geologic history of the area and a critical review of investigative studies conducted in the area and beyond. Predicated on the region’s structural history, some undiscovered elements that can potentially cause ﬁssuring in the region were identiﬁed and discussed. Some of these include the pre-existence of a fault system, a crack from the bedrock ridge, the existence of paleochannels, the collapsibility of loess, the tectonic (earthquake) history of the area, and differential compaction due to heterogeneity. The use of a metaheuristic and a combined application integrating other optimization algorithms can be utilized to determine optimum hyperparameters and present their statistical importance, thereby improving accuracy and dependability in ﬁssure prediction in Najran. Reliable models would primarily be used to monitor active ﬁssures and identify key factors utilizing spatial information, subsidence, groundwater-related data sets, etc.


Introduction
Recently, earth fissures, which are deformations or cracks at the Earth's surface that occur in regions with little or no rain, mainly caused by the immoderate withdrawal or pumping of groundwater [1,2], have become one of the disturbing environmental problems across different cities in the world. It causes damage to homes and businesses, highways, bridges, waterways, drains, powerlines, farmland, pipeline infrastructure, animals, ecosystems, biodiversity, and individual lives, as well as providing a pathway for toxins to enter groundwater. According to [2][3][4][5][6][7][8], several researchers have confirmed that groundwater is vulnerable to impurity through earth fissures, and cleaning these contaminants can be an arduous and time-consuming task. As a result, this geological risk has had a huge impact on business activities, social security, and environmental conservation, garnering more attention in recent decades [9]. The occurrence of and damage caused by earth fissures have been recorded in different places around the world, including the southeast corner of the North China Plain, the southwest of the American Basin, the North Indian River Plain, The formation of earth fissures in Qinglong Graben is divided into three steps, namely (1) regional extension creating normal faults under the surface, (2) inadequate groundwater extraction causing normal faults to rupture the surface, and (3) erosion widening the cracks [23] 3 Haram-to-Haram Highway, Iran Extensive field mapping, on-site experiments, sampling, and evaluation of soil in the laboratory.
Fissure development and expansion in the studied area were aided by a combination of shallow and deep mechanisms as well as other environmental agents, with soil features playing a key role. [24] 4 Central and northern Arabia Integration of field, interferometric synthetic-aperture radar (InSAR), groundwater geology, geodetic studies, and geospatial analysis were conducted to decipher the cause of fissuring.
The fissures were mostly constrained by a graben fault system which played a role in the development of the fissures. Fissuring activity was triggered by aquifer depletion. [25] 5 Weihe Basin, China Field mapping, trench excavation, geophysics, and geotechnical testing methods were applied The ground fissures mirrored the deeply buried faults on the surface. Loess collapsibility and seepage were involved in the mechanism of the fissures. Differential settlement occurred in the collapsible and non-collapsible portions of the loess, with ground fissures serving as the borderline. [26] 6 Su-Xi-Chang, southeast corner of northern China The geologic history and piezometric features of the area were studied. Based on this, a physics model was used to evaluate the mechanism of ground fissure creation.
The creation of fissures in the location is caused by repeated groundwater impoundment and pumping, as well as peculiar geological conditions. [18] 7 Yingxian area, Datong Basin, China Trench excavation analysis, spudding, seismic exploration techniques, and the geologic history of the area were employed.
The findings and interpretations show that (1) the fissures are connected to two large active faults that provide the background for the fissure formation; (2) intensified groundwater harvesting has revived the underlying faults and exacerbated fissure growth; and (3) assuagement is best achieved by ceasing excess groundwater exploitation, suitable engineering design, and avoiding major structure location.
[27] To determine the origins of fissure formation, an integrated ground survey with trenches and geophysical prospecting profiles was conducted.
The pre-existence of buried channels and faults offered a favorable condition for the origination of fissures, and groundwater pumping further accelerated the expansion and the fissures. [28] 9 Shuanghuaishu, Shaanxi Province, China An on-site water immersion test was performed to study the cause of fissures The creation of collapsible earth fissures resulted from tensile tension and strain caused by the unequal collapse of loess. The loess encountered tensile failure and developed an earth fissure when the force was greater than the eventual deformation [29] 10 Wuxi, China Stress displacement assessment using an interface infinite model that was developed from an impermeable and non-compressible rock cliff. Aquifer depletion and the existence of pre-existing sedimentary deposits such as lake deposits promote friable soil compaction. [32] From the reviews and ongoing discussion among scientists about fissure formation processes, it can be summarized that (a) deep dynamic tectonism controls ground fissure location, (b) the internal tectonic force of the earth's crust leads to a collection of ground fissures, (c) stresses at the fault planes induces the formation of ground fissures, (d) groundwater overmining results in the reactivation and expansion of subsurface faults, leading to fissure formation, and (e) buried paleo geomorphology aids the collapsibility of soils and formation of fissures. The assemblage and appearance of ground fissures are thus majorly influenced by the geological environment. More crucially, the combined impacts of interior geological forces (such as crustal instability, tectonic forces within a craton, and motions at fault zones) and human-instigated stresses, such as excessive groundwater extraction, heavy engineering structures, etc., cause ground fissure propagation.
A more detailed understanding of the mechanism and causative factors of fissures would assist in designing suitable ways to mitigate the fissure hazards. This paper (1) identifies and discusses five major factors that cause ground fissures by reviewing scientific studies conducted by different researchers, (2) relates these factors to suggested causes of fissures in the city of Najran for a better understanding of its development mechanism, (3) shares insightful conclusions about other potential factors that could contribute to the development of fissures in in the city of Najran, (4) presents honest opinions about fissure occurrences in Najran, and (5) makes recommendations on potential approaches that could be taken to study, monitor, and mitigate fissuring in the city.

Geological Background
Najran is the administrative center of the Najran region in the southwest of Saudi Arabia, close to the Yemeni border (Figure 1a,b). Najran is one of the kingdom's extremely fast-developing municipalities, with a populace that has increased nearly tenfold across an area of around 360,000 km 2 . In the last four decades, both the government and the private sector have invested significantly in Najran. The area is primarily an agricultural zone in Wadi Najran's floodplain. The elevation rises to 1447 m, with elevations ranging from 883 m east to 2330 m west [2]. Three distinct geomorphic groups exist, namely (a) elevated zones in the western part, (b) lowland regions with sediments deposited across the Wadi, and (c) up the east dominated by sandy dunes. Najran is composed of volcanic Precambrian rocks, as well as Wajeed sandstone deposited during Cambrian-Ordovician age, and the later deposited tertiary bedrock [33,34]. Wadi Najran is covered in alluvial quaternary sediments, especially dunes prevailing eastward near the Empty Quarter fields atop alluvial soils.
In the study area, Wajid sandstone is found in a sequence of hill and rock formations at scattered outcrops. It lays unconformably atop the roughly flat basement rocks of the Shield of Arabia. The topmost exposure of the Wajid sandstone formations does not quite extend into the research region, but it does unconformably underlie the Permian Khuff Formation to the northeast. The Wadi Najran region is made up of Arabian-Nubian basement massif rocks such as granites, metamorphosed andesite, basalt, gneisses, granitoid formed during orogeny, and gabbro and diorite volcanic rocks. The Arabian plate's N-S trending structures were created by extensional tectonic movements of the Arabian-Nubian Shield [35,36]. The NW-SE trending dextral fault system of Najid was created in the very last phase of its rifting tectonic system, after which the Arabian Shield was made almost plain, and the Wajid Group was accumulated during the Paleozoic in the southern region of Saudi Arabia [37]. Numerous tectonic events, mostly in the Phanerozoic, reactivated the Precambrian underlying structures [38][39][40][41]. The Wajid Group's regional-scale fissures are most likely the result of the reactivation of those structures [42]. Sedimentary rocks cover the eastern section of the Najran region as well (Figure 1b).
The Wajid Group is a siliciclastic sequence that was deposited during the Cambrian and Permian ages and underlies the southwestern zone of Saudi Arabia [43,44]. Its exposures extend from Wadi Al-Dawasir in the south of Najran, with just a small exceptional case in some areas [44] (Figure 1b). The Wajid Group in Najran is divided into five stratigraphic units, from the earliest to the most recent formation, and they include Dibsiyah, Sanamah, Qalibah, Khusayyayn, and Juwayl [43,45], as shown in Figure 2. Based on the previous lithologic description, the lower Dibsiyah Formation is a reddish cross-stratified sandstone with some conglomerate layers deposited in some kind of fluvial channels. The topmost of Dibsiyah Formation's intercalated massive and bioturbated sandstones were accumulated in transitional shallow marine settings. The Sanamah Formation is composed of coarsegrained conglomeratic sandstones, considered a mix of glacial and fluvial sediments. The Khusayyayn Formation is composed of river-and wind-deposited sandstone with grain sizes ranging from medium to coarse. The Juwayl Formation is massive and/or layered, consisting of grains that range from fine to coarse, deposited in a glaciofluvial setting [42,[46][47][48]. The Group's sandstones range in texture from fine to coarse, are slightly sorted, and are partly rounded to sub-angular in form. Quartz arenites that are composed of limited percentages of feldspars form the basal units; however, Khusayyayn Formation sandstone is a subarkose sandstone having a highly disproportionate amount of feldspar grains [42], as illustrated in Figure 2. In the study area, Wajid sandstone is found in a sequence of hill and rock formations at scattered outcrops. It lays unconformably atop the roughly flat basement rocks of the Shield of Arabia. The topmost exposure of the Wajid sandstone formations does not quite extend into the research region, but it does unconformably underlie the Permian Khuff Formation to the northeast. The Wadi Najran region is made up of Arabian-Nubian basement massif rocks such as granites, metamorphosed andesite, basalt, gneisses, granitoid formed during orogeny, and gabbro and diorite volcanic rocks. The Arabian

The Causes and Mechanisms of Fissures Formation
To better comprehend the fissuring process and establish a plan to mitigate earth fissuring and, consequently, reduce hazards, it is essential to identify the elements that influence earth fissuring. However, it is unclear how these factors affect the development and evolution of earth cracks. In the intervening decades, several theories have been put forth to elucidate how groundwater loss causes surface faulting as well as earth fissuring. Some include localized vertical differential compaction, regional variable compaction, and capillary stress related to a lowering water table [51]. It has also been proposed by [52] and [53] that aquifer displacement, particularly lateral movement, is the cause of earth fissures associated with groundwater depletion. The complicated processes of earth fissuring resulting from inordinate groundwater drawdown, as perhaps the major influential anthropogenic factor, are quite well discussed in the literature [54][55][56][57][58][59][60][61][62]. By pressurizing susceptible aquifer systems, groundwater exploitation significantly contributes to soil slope failure and fissure emergence. The subsidence associated with the extraction of liquids such as water from subterranean strata is likely the most well-studied of all anthropogenic and natural sources of subsidence [63]. In aquifer systems, lateral deformation is caused by removal and occasional recharge/discharge pressures [17,[64][65][66][67]. Earth fissure is most often provincial and widespread in aquifer systems undergoing compaction, with provincial lateral stresses rarely exceeding 2 ppm and regional misalignments rarely surpassing 240 arcseconds. Linear strains can be substantial locally, especially near pumping boreholes with higher pressure gradients, along with inflections in the bedrock landforms of underground water systems, and near the borders of hydrogeologic units with differing hydraulic and (or) mechanical characteristics [63]. Enough scientific work is necessary to clarify the role of lateral strains in the formation of geological fissures in basins predisposed to aquifer system compression. Research has shown that not only horizontal deformation but shear on vertical planes and rotating stresses also play a part in the development of certain earth cracks [20]. However, scientists have argued that lateral displacement of the earth's surface was only detected in a few spots in areas of subsidence related to groundwater depletion. They maintain that lithologic differences, such as the amount of consolidation or cementation in alluvial deposits [68], and subsurface structures such as faults, both influence the formation and placement of fissures [58], whereas the inordinate exploitation of groundwater only activates them. Based on previous studies, this paper categorizes the mechanisms of fissure formation into three, namely tectonic, non-tectonic, and anthropogenic (Table 2). Figure 3 illustrates how these mechanisms are interrelated and dependent on each other, with groundwater pumping as the central factor that triggers other processes.  The following subsections discuss these processes and how they are triggered by disequilibrium caused by groundwater over-pumping and other natural hydrological conditions.

Crack from the Bedrock Ridge
The strata structure of the aquifer system is influenced by the bedrock ridge. The large differential settlement arises because of sediment consolidation under higher effective pressures from groundwater drawdown when groundwater is removed disproportionately for a long time. Bending can cause tensile deformation in sediments above the bedrock ridge, as shown in Figure 4. A fissure occurs when tension surpasses the tensile strain limit of soils. The elastic deformation of sediments under plane strain circumstances can be calculated using Equations (1)-(3), as the fissure develops longitudinally along the bedrock ridge.
where u = lateral displacement (m), v = vertical displacement (m), Ex = perpendicular strain in the horizontal direction, Ey = perpendicular strain in the vertical direction, and Ƿxy = shear deformation (strain) The following subsections discuss these processes and how they are triggered by disequilibrium caused by groundwater over-pumping and other natural hydrological conditions.

Crack from the Bedrock Ridge
The strata structure of the aquifer system is influenced by the bedrock ridge. The large differential settlement arises because of sediment consolidation under higher effective pressures from groundwater drawdown when groundwater is removed disproportionately for a long time. Bending can cause tensile deformation in sediments above the bedrock ridge, as shown in Figure 4. A fissure occurs when tension surpasses the tensile strain limit of soils. The elastic deformation of sediments under plane strain circumstances can be calculated using Equations (1)-(3), as the fissure develops longitudinally along the bedrock ridge.  The following subsections discuss these processes and how they are triggered by disequilibrium caused by groundwater over-pumping and other natural hydrological conditions.

Crack from the Bedrock Ridge
The strata structure of the aquifer system is influenced by the bedrock ridge. The large differential settlement arises because of sediment consolidation under higher effective pressures from groundwater drawdown when groundwater is removed disproportionately for a long time. Bending can cause tensile deformation in sediments above the bedrock ridge, as shown in Figure 4. A fissure occurs when tension surpasses the tensile strain limit of soils. The elastic deformation of sediments under plane strain circumstances can be calculated using Equations (1)-(3), as the fissure develops longitudinally along the bedrock ridge.
where u = lateral displacement (m), v = vertical displacement (m), Ex = perpendicular strain in the horizontal direction, Ey = perpendicular strain in the vertical direction, and Ƿxy = shear deformation (strain) where u = lateral displacement (m), v = vertical displacement (m), Ex = perpendicular strain in the horizontal direction, Ey = perpendicular strain in the vertical direction, and PEER REVIEW 10 of 23 The following subsections discuss these processes and how they are triggered by disequilibrium caused by groundwater over-pumping and other natural hydrological conditions.

Crack from the Bedrock Ridge
The strata structure of the aquifer system is influenced by the bedrock ridge. The large differential settlement arises because of sediment consolidation under higher effective pressures from groundwater drawdown when groundwater is removed disproportionately for a long time. Bending can cause tensile deformation in sediments above the bedrock ridge, as shown in Figure 4. A fissure occurs when tension surpasses the tensile strain limit of soils. The elastic deformation of sediments under plane strain circumstances can be calculated using Equations (1)-(3), as the fissure develops longitudinally along the bedrock ridge.
Ex = dv/dy; Ey = du/dx; (1) where u = lateral displacement (m), v = vertical displacement (m), Ex = perpendicular strain in the horizontal direction, Ey = perpendicular strain in the vertical direction, and Ƿxy = shear deformation (strain) = shear deformation (strain)  For illustration, [30] studied land subsidence and the accompanying ground fissures in the Suzhou, Wuxi, and Changzhou (Su-Xi-Chang) areas of southern Jiangsu Province in China in 2010. Excessive groundwater drawdown, in combination with the underlying bedrock ridge, was identified as the causal element in previous site investigations and geophysical surveys. As per new findings, earth fissures may have ruptured from the underlying bedrock cliffs or ridges to the surface of the ground. Ref. [68] asserted that previously postulated mechanisms to produce fissures from groundwater were predicated on assumptions and abstract justifications which might or might not be coherent with the physics of fissure formation from water table decrease. In the research, the authors aimed to better understand the mechanical mechanism which culminates in the creation of earth fissures in an unconfined aquifer due to fluctuating groundwater level fall, as well as the important elements that regulate the commencement of earth fissures. Ref. [68] found that earth cracks in aquifers emerging from groundwater overpumping are created by the effect of shearing on vertical planes and spinning, rather than tension, as is usually assumed. It demonstrates how the gradient of the subduction bowl, which is composed of simple shear deformation and spin can be used to anticipate the onset of earth fissures. This finding supported the idea that applied stress could cause earth fissures to fracture from the bedrock cliff face to the outer surface of the earth. Ref. [56] shed more light by outlining two distinct mechanisms that are involved in the formation of fissures influenced by bedrock ridge: first, bending or draping of horizontal bedding planes over the ridge by differential vertical compaction and, second, tensile failure at depth along the top of the ridge caused by horizontal movement of the aquifer in at least one direction away from the ridge. According to the study by [69], the upper layers on either one or both sides of the ridge drape or spin in a reverse orientation as a result of the vertical compression at a range from the ridge. Aquifer fabric may be stretched horizontally over a rigid bedrock ridge as a consequence of uneven vertical compression close to the ridge, creating a tension region in which the aquifer material is thinnest. The latter process is predicated on a groundwater aquifer that is moving horizontally and is partially or completely terminated by the bedrock. This denotes a fresh impetus for the beginning and propagation of an earth fissure. Additionally, this horizontal displacement may result in a deep tensile failure. If there are two pumping centers on each side, this failure might happen in the material next to the ridge in the axis For illustration, [30] studied land subsidence and the accompanying ground fissures in the Suzhou, Wuxi, and Changzhou (Su-Xi-Chang) areas of southern Jiangsu Province in China in 2010. Excessive groundwater drawdown, in combination with the underlying bedrock ridge, was identified as the causal element in previous site investigations and geophysical surveys. As per new findings, earth fissures may have ruptured from the underlying bedrock cliffs or ridges to the surface of the ground. Ref. [68] asserted that previously postulated mechanisms to produce fissures from groundwater were predicated on assumptions and abstract justifications which might or might not be coherent with the physics of fissure formation from water table decrease. In the research, the authors aimed to better understand the mechanical mechanism which culminates in the creation of earth fissures in an unconfined aquifer due to fluctuating groundwater level fall, as well as the important elements that regulate the commencement of earth fissures. Ref. [68] found that earth cracks in aquifers emerging from groundwater over-pumping are created by the effect of shearing on vertical planes and spinning, rather than tension, as is usually assumed. It demonstrates how the gradient of the subduction bowl, which is composed of simple shear deformation and spin can be used to anticipate the onset of earth fissures. This finding supported the idea that applied stress could cause earth fissures to fracture from the bedrock cliff face to the outer surface of the earth. Ref. [56] shed more light by outlining two distinct mechanisms that are involved in the formation of fissures influenced by bedrock ridge: first, bending or draping of horizontal bedding planes over the ridge by differential vertical compaction and, second, tensile failure at depth along the top of the ridge caused by horizontal movement of the aquifer in at least one direction away from the ridge. According to the study by [69], the upper layers on either one or both sides of the ridge drape or spin in a reverse orientation as a result of the vertical compression at a range from the ridge. Aquifer fabric may be stretched horizontally over a rigid bedrock ridge as a consequence of uneven vertical compression close to the ridge, creating a tension region in which the aquifer material is thinnest. The latter process is predicated on a groundwater aquifer that is moving horizontally and is partially or completely terminated by the bedrock. This denotes a fresh impetus for the beginning and propagation of an earth fissure. Additionally, this horizontal displacement may result in a deep tensile failure. If there are two pumping centers on each side, this failure might happen in the material next to the ridge in the axis of pumpage, but it may also happen right over the is horizontal displacement which is a function of S, and H and L are vertical and horizontal thickness, respectively.
For illustration, [30] studied land subsidence and the accompanying ground fissures in the Suzhou, Wuxi, and Changzhou (Su-Xi-Chang) areas of southern Jiangsu Province in China in 2010. Excessive groundwater drawdown, in combination with the underlying bedrock ridge, was identified as the causal element in previous site investigations and geophysical surveys. As per new findings, earth fissures may have ruptured from the underlying bedrock cliffs or ridges to the surface of the ground. Ref. [68] asserted that previously postulated mechanisms to produce fissures from groundwater were predicated on assumptions and abstract justifications which might or might not be coherent with the physics of fissure formation from water table decrease. In the research, the authors aimed to better understand the mechanical mechanism which culminates in the creation of earth fissures in an unconfined aquifer due to fluctuating groundwater level fall, as well as the important elements that regulate the commencement of earth fissures. Ref. [68] found that earth cracks in aquifers emerging from groundwater over-pumping are created by the effect of shearing on vertical planes and spinning, rather than tension, as is usually assumed. It demonstrates how the gradient of the subduction bowl, which is composed of simple shear deformation and spin can be used to anticipate the onset of earth fissures. This finding supported the idea that applied stress could cause earth fissures to fracture from the bedrock cliff face to the outer surface of the earth. Ref. [56] shed more light by outlining two distinct mechanisms that are involved in the formation of fissures influenced by bedrock ridge: first, bending or draping of horizontal bedding planes over the ridge by differential vertical compaction and, second, tensile failure at depth along the top of the ridge caused by horizontal movement of the aquifer in at least one direction away from the ridge. According to the study by [69], the upper layers on either one or both sides of the ridge drape or spin in a reverse orientation as a result of the vertical compression at a range from the ridge. Aquifer fabric may be stretched horizontally over a rigid bedrock ridge as a consequence of uneven vertical compression close to the ridge, creating a tension region in which the aquifer material is thinnest. The latter process is predicated on a groundwater aquifer that is moving horizontally and is partially or completely terminated by the bedrock. This denotes a fresh impetus for the beginning and propagation of an earth fissure. Additionally, this horizontal displacement may result in a deep tensile failure. If there are two pumping centers on each side, this failure might happen in the material next to the ridge in the axis of pumpage, but it may also happen right over the ridge's highest point where the aquifer material is moving away from the ridge on both flanks. As the horizontal aquifer migration intensifies, so any collapse might ultimately migrate upward and start a fracture in the layers above.

Pre-Existence of Faults
Existing faults in constant shear states may become more active and stretch to the ground surface, causing cracks in the hanging fault block. Besides that, sustained aquifer depletion in such an area can cause the faults' hanging wall to consolidate and settle, increasing the unequal settlement on both blocks of the faults and hastening the mobility and formation of the ground fissure [70]. Fault zones and buried fault scarps can have an impact on the accumulation of stresses over time. The faulted zones have a considerable influence on extensional and compressional patterns, whereas the buried fault scarp has a greater impact on shear stress patterns. Ref. [58] discovered that a fissure caused by tension is near the faulted block at the sliding wall, affirming that fissures begin near the surface and make their way downward, terminating near the boundary between the saturated and unsaturated column. The presence of a fault zone has been established as a determinant in a fissure's propagation course. This was further confirmed by the authors of [28], who used site evaluation, drilling, and trench unearthing to define and assess the basis for fissure formation in Qinglong Gaben, China. Syn-sedimentary faults, as well as an excessive drawdown of groundwater and erosion processes, were shown to be major drivers in the creation and multiplication of fissures in the region, according to the research. Recent research has revealed that earth fissure cluster formation is polygenetic and the result of the earth's internal forces and anthropogenic actions. The origination of fault-induced earth cracks can be subdivided into three developmental phases: first, initial extension resulting from fault planes underground, followed by the excessive pumping of groundwater, and then, washouts aided by the development of the fissure. Deep crustal movements underneath the basin influence earth fissuring, which is driven by the local geodynamic stress of the nearby fault zones [71]. The fault displacement increases as the depth of the fault grows, which is typical of synthetic and other syn-depositional faults. Vertical displacement characterizes these fissures, which can have a yearly growth rate that ranges from 1 to 3 cm. The tectonic structure of a basin's hidden faults is thought to have an influence on the formation of large-scale fissures, which may also be influenced by the basin's regional extensional stress. Over-exploitation of groundwater further increases the level of activities leading to fissure formation [72]. As a result, groundwater exploitation is an important contributor to fault-induced earth fissure development. Intensive groundwater extraction results in the establishment of an inner multi-directional (radial) compression zone near pumping wells and an external radial tensional zone outside them. Lateral tensile stress (σ) minimizes the resisting force on the fissure surface and increases its activation in fissures in the external radial tensional zone. If a well is bored in one of the fault compartments, for example, the rate of the water table decrease in the sliding (or hanging) wall will be faster than it is in the footwall, resulting in more shear stress (τ) in the vertical plane of the fissure. Unequal surface settling (∆H > ∆F) causes tensional stress around the ground fissure, which encourages fissure migration and expansion as shown in Figure 5. The sediment thickness in the hanging wall is generally greater than thickness of the accumulated sediment in the footwall due to the build-up of pre-existing faults. Fissuring can form aberrant layer placements and subsequently stretch to pre-existing faults as groundwater pumping intensifies.

Existence of Paleochannels
Paleochannels affect the creation of rock fissures, which are commonly caused by intense rains, squeezing, or earthquakes, and inadequate groundwater pumping [73,74]. It creates an ideal environment for fissure formation [71]. Because of the inhomogeneity, anisotropy, and friability of shallowly buried ancient channel deposits, instabilities and uneven settling are common, resulting in fissure development, piping, and surface erosion, among other things, and human activities hasten this process. Earth fissures linked with paleochannels have short lengths and superficial depths overall, as well as the forming trend of earth fissures on the edges of paleochannels, which reflects the channel's meandering model [73,75]. Due to gravitational force and extra pressures, the top has strata curves and breaks above the empty spaces when unequal vertical strain occurs in the compressible layer. According to a theoretical examination of their characteristics, channel controlled fissures are generated by gravity, initial joint stresses, and groundwater withdrawal operations, and are driven by paleo-fluvial factors [74]. The base of paleochannel-induced fissures usually coincides with curving or meandering surfaces or gradational zones between valley slopes and terraces, implying that they are generated by the self-weight of overlying strata [75]

Collapsibility of Loess
In recent years, experts have become increasingly interested in nontectonic ground cracks. Nontectonic ground fissures were explored, and their causes were investigated in [6,29]. They used field study, trenching, and geotechnical testing to conclude that nontectonic ground cracks are mostly caused by loess collapsibility [6], determined that earth cracks are produced by loess collapsibility and erosion is caused by infiltration, and that loess cracks create a favorable environment for the creation of earth fissures. The studies reached these conclusions by studying groundwater variability, topology, soil physico-mechanical indexes, as well as the lithological backdrop. Ref. [26] looked at the relationship connecting four soft soil properties and earth fissuring: density of dry loess, its porosity, water saturation, and self-gravity collapsibility coefficient. They described how earth fissures are generated by employing the spinning mechanism of a loess cantilever beam, as well as the computation equation of tensile force. Although the following research on fissures caused by factors other than tectonic forces has distinct foci

Existence of Paleochannels
Paleochannels affect the creation of rock fissures, which are commonly caused by intense rains, squeezing, or earthquakes, and inadequate groundwater pumping [73,74]. It creates an ideal environment for fissure formation [71]. Because of the inhomogeneity, anisotropy, and friability of shallowly buried ancient channel deposits, instabilities and uneven settling are common, resulting in fissure development, piping, and surface erosion, among other things, and human activities hasten this process. Earth fissures linked with paleochannels have short lengths and superficial depths overall, as well as the forming trend of earth fissures on the edges of paleochannels, which reflects the channel's meandering model [73,75]. Due to gravitational force and extra pressures, the top has strata curves and breaks above the empty spaces when unequal vertical strain occurs in the compressible layer. According to a theoretical examination of their characteristics, channel controlled fissures are generated by gravity, initial joint stresses, and groundwater withdrawal operations, and are driven by paleo-fluvial factors [74]. The base of paleochannel-induced fissures usually coincides with curving or meandering surfaces or gradational zones between valley slopes and terraces, implying that they are generated by the self-weight of overlying strata [75].

Collapsibility of Loess
In recent years, experts have become increasingly interested in nontectonic ground cracks. Nontectonic ground fissures were explored, and their causes were investigated in [6,29]. They used field study, trenching, and geotechnical testing to conclude that nontectonic ground cracks are mostly caused by loess collapsibility [6], determined that earth cracks are produced by loess collapsibility and erosion is caused by infiltration, and that loess cracks create a favorable environment for the creation of earth fissures. The studies reached these conclusions by studying groundwater variability, topology, soil physico-mechanical indexes, as well as the lithological backdrop. Ref. [26] looked at the relationship connecting four soft soil properties and earth fissuring: density of dry loess, its porosity, water saturation, and self-gravity collapsibility coefficient. They described how earth fissures are generated by employing the spinning mechanism of a loess cantilever beam, as well as the computation equation of tensile force. Although the following research on fissures caused by factors other than tectonic forces has distinct foci and creation mechanisms, they all agree that nontectonic fissure formation is linked to the ability of loess to collapse and allow fluids to infiltrate. According to [24], all fissures formed by soil collapsibility go through a phase from emergence to maturity, starting with microcracks and ending with eroded fissures as broad as gullies. It is said that as surface water accumulates and infiltrates into the ground, the wetness of the soil increases while its tensile strength falls, causing fissures to form. The saturation of collapsible soil occurs because of water infiltration, which is followed by hydro-compaction and the creation of new nascent hairline fissures. Additional water infiltrates inside the subsurface through tiny cracks, eroding telescopic soils and forming underground passageways that are not obvious on the earth's surface in the initial state of growth and can only be observed through some cave-ins along the fissures. More water erosion broadens and increases the subway, causing the roofs to collapse and forming pits along the fissures. Finally, owing to piping erosion, these potholes connect, producing more collapse. Continuous erosion enlarges the collapsed structures and transforms them into eroded gullies, which are the mature stage of the fissures.

Differential Compaction Due to Heterogeneity
Aquifer heterogeneity results in a sudden shift in an aquifer's thickness or in its bed composition. According to [53], such thickness fluctuations could lead to rotation, vertical shear, or even lateral expansion at depth. For instance, the author in [56] illustrated that the vertical compaction strength varies on either side of a geometric heterogeneity. He explained that localized differential vertical movement in the underlying, non-compacting interval is caused by a geometric anomaly that migrates upward, causing tilting and shearing at the surface of the terrain. On the other hand, localized differential horizontal distortions may cause fissures at the ground surface as well as a deep extensional zone and induce an opposite direction of rotation (see Figure 6). and creation mechanisms, they all agree that nontectonic fissure formation is linked to the ability of loess to collapse and allow fluids to infiltrate. According to [24], all fissures formed by soil collapsibility go through a phase from emergence to maturity, starting with microcracks and ending with eroded fissures as broad as gullies. It is said that as surface water accumulates and infiltrates into the ground, the wetness of the soil increases while its tensile strength falls, causing fissures to form. The saturation of collapsible soil occurs because of water infiltration, which is followed by hydro-compaction and the creation of new nascent hairline fissures. Additional water infiltrates inside the subsurface through tiny cracks, eroding telescopic soils and forming underground passageways that are not obvious on the earth's surface in the initial state of growth and can only be observed through some cave-ins along the fissures. More water erosion broadens and increases the subway, causing the roofs to collapse and forming pits along the fissures. Finally, owing to piping erosion, these potholes connect, producing more collapse. Continuous erosion enlarges the collapsed structures and transforms them into eroded gullies, which are the mature stage of the fissures.

Differential Compaction Due to Heterogeneity
Aquifer heterogeneity results in a sudden shift in an aquifer's thickness or in its bed composition. According to [53], such thickness fluctuations could lead to rotation, vertical shear, or even lateral expansion at depth. For instance, the author in [56] illustrated that the vertical compaction strength varies on either side of a geometric heterogeneity. He explained that localized differential vertical movement in the underlying, noncompacting interval is caused by a geometric anomaly that migrates upward, causing tilting and shearing at the surface of the terrain. On the other hand, localized differential horizontal distortions may cause fissures at the ground surface as well as a deep extensional zone and induce an opposite direction of rotation (see Figure 6). Figure 6. Heterogeneous layer of a highly compressible clay causes differential compaction, leading to vertical shear, rotation, and horizontal extension [56]. Figure 6. Heterogeneous layer of a highly compressible clay causes differential compaction, leading to vertical shear, rotation, and horizontal extension [56].
The ongoing discussion is summarized in Table 3, which highlights the major factors and processes of fissure formation as discussed above. Other anthropogenic factors and processes such as groundwater over-pumping is an integral part of the processes discussed above since they facilitate them.

Earth Fissures in Najran
From the previous discussion, we can see that it is a popular opinion among researchers that groundwater depletion makes a significant contribution to the formation of fissures. Researchers have attributed inadequate groundwater removal as the primary cause of land subsidence and earth fissures in several locations throughout the Kingdom of Saudi Arabia [2,32,34,76]. Ref. [34] highlighted the diverse kinds of ground fissures (Figure 7 A study conducted in [2] estimated that the city of Najran witnesses only about 31.3 × 106 m 3 rainfall per year. This is the reason for the constant reduction in water table levels, which caused the drying of more than a thousand wells in the city. Rapid groundwater depletion remains one of the most profound sustainability challenges facing many cities. In Najran, the authors in [2] indicated that the Najran fissures have varying lengths and widths; one of the fissures elongated to about 600 m in a northeast direction and has a varying breadth of thirty or fifty centimeters (see Figure 8). Although these reported initial widths of the fissures might seem small, the authors in [77] asserted from their investigation of fissures in south-central Arizona that similar fissures being initially narrow and having a width of about 2 cm wide and 1-2 km long, widen and elongate as time passes. As discussed above, in Section 2, the length of the fissures or the extent of deformations resulting from fissuring activity typically depends on the underlying geology of the area, with causes such as pre-existing faults, the degree of groundwater depletion, and the A study conducted in [2] estimated that the city of Najran witnesses only about 31.3 × 106 m 3 rainfall per year. This is the reason for the constant reduction in water table levels, which caused the drying of more than a thousand wells in the city. Rapid groundwater depletion remains one of the most profound sustainability challenges facing many cities. In Najran, the authors in [2] indicated that the Najran fissures have varying lengths and widths; one of the fissures elongated to about 600 m in a northeast direction and has a varying breadth of thirty or fifty centimeters (see Figure 8). Although these reported initial widths of the fissures might seem small, the authors in [77] asserted from their investigation of fissures in south-central Arizona that similar fissures being initially narrow and having a width of about 2 cm wide and 1-2 km long, widen and elongate as time passes. As discussed above, in Section 2, the length of the fissures or the extent of deformations resulting from fissuring activity typically depends on the underlying geology of the area, with causes such as pre-existing faults, the degree of groundwater depletion, and the resultant pressure imbalance. As reported from a survey carried out in Su-Xi-Chang Plain, in the southeastern part of the North China Plain, earth fissures can be up to two thousand meters long [20]. It is therefore acceptable to presume that earth fissures in Najran might become worse than they are presently.
Su-Xi-Chang Plain, in the southeastern part of the North China Plain, earth fissures ca be up to two thousand meters long [20]. It is therefore acceptable to presume that ear fissures in Najran might become worse than they are presently. In Najran, most exposed ground fissures caused by groundwater depletion appe on the surface as linear ruptures with some vertical displacements, and the resultin surface scarps result in myriad destruction on roadways, farms, and residences with range of influence of 30 to 400 m. Due to tensile forces caused by excessive groundwat extraction, earth cracks in the city of Najran can stretch for tens of meters to kilomete [51]. Studies have shown that Najran's earth fissures are linked to groundwater depletio ((Interdisciplinary Earthquake Hazard Research in Gulf of Aqaba and Strait of Tira (GAST)-NASA/ADS, n.d.) [34,78]. The rapid drop in the water table level cause subsidence as well as ground failure, which was evidenced by earth rifts, owing to th aquifer's composition of poorly consolidated sands with high porosity and aquitards wi low pore connectivity and high compressibility [2,34]. If increased water exploitation the city of Najran persists, fissuring and some other environmental threats includin sinkholes, seismic events, and structural failure to engineering constructions will al persevere, causing the size of existing subsidence and fissures throughout the location grow even larger [15].
Research has proven that factors causing the formation of fissures are synergist implying that, although groundwater depletion triggers or enhances the process, oth pre-existing geologic tectonic and nontectonic factors play major roles in enabling the formation [23]. Apart from the excessive drawdown of the groundwater table in Najra the subsurface topography of the study area has been also suggested as a major fact contributing to the formation of the fissures in that area. As discussed earlier, subsurfa slopes such as bedrock ridges can influence the formation of fissures [21,58]. Th secondary cause of the Najran fissure was proposed by the authors in [2], who use advanced electrical resistivity tomography (ERT) to evaluate the subsurface exten topological information, and the spread of fissures in the Najran metropolis.
Several regional and local active seismic zones have been reported in the weste region of Saudi Arabia, including the Red Sea, the Gulf of Aqaba, and the Najid Fau In Najran, most exposed ground fissures caused by groundwater depletion appear on the surface as linear ruptures with some vertical displacements, and the resulting surface scarps result in myriad destruction on roadways, farms, and residences with a range of influence of 30 to 400 m. Due to tensile forces caused by excessive groundwater extraction, earth cracks in the city of Najran can stretch for tens of meters to kilometers [51]. Studies have shown that Najran's earth fissures are linked to groundwater depletion ((Interdisciplinary Earthquake Hazard Research in Gulf of Aqaba and Strait of Tiran (GAST)-NASA/ADS, n.d.) [34,78]. The rapid drop in the water table level caused subsidence as well as ground failure, which was evidenced by earth rifts, owing to the aquifer's composition of poorly consolidated sands with high porosity and aquitards with low pore connectivity and high compressibility [2,34]. If increased water exploitation in the city of Najran persists, fissuring and some other environmental threats including sinkholes, seismic events, and structural failure to engineering constructions will also persevere, causing the size of existing subsidence and fissures throughout the location to grow even larger [15].
Research has proven that factors causing the formation of fissures are synergistic, implying that, although groundwater depletion triggers or enhances the process, other pre-existing geologic tectonic and nontectonic factors play major roles in enabling their formation [23]. Apart from the excessive drawdown of the groundwater table in Najran, the subsurface topography of the study area has been also suggested as a major factor contributing to the formation of the fissures in that area. As discussed earlier, subsurface slopes such as bedrock ridges can influence the formation of fissures [21,58]. This secondary cause of the Najran fissure was proposed by the authors in [2], who used advanced electrical resistivity tomography (ERT) to evaluate the subsurface extent, topological information, and the spread of fissures in the Najran metropolis.
Several regional and local active seismic zones have been reported in the western region of Saudi Arabia, including the Red Sea, the Gulf of Aqaba, and the Najid Fault system in addition to many others that can potentially accommodate different magnitudes of earthquakes such as the destructive 1993 (Mw 6.1) and 1995 (Mw 7.3) earthquakes that occurred in the Gulf of Aqaba [79,80]. These seismic zones activate tectonic deformation of the earth and consequent faulting across the region, and the city of Najran, being in the western region, is not entirely free from their impacts. The Wajid Group which underlies the city of Najran is reported to be fractured on a regional scale as a result of the paleo reactivation of the Arabian Shield structures [42]. These structures can be easily reactivated due to the stress imbalance accruing from excessive groundwater pumping. Although it has been suggested that subsurface topography is a key factor in the fissure's formation, it is not exhaustive and can be an insignificant contributor to the fissuring process. It is, therefore, required to conduct fault or fracture mapping studies across the fissuring zones. An in-depth understanding of the structural architecture and fault history of the Najran Basin would be an immense guide toward understanding other unknown factors responsible for fissuring in the area. There is a need to conduct further research into the other potential causes of this hazard. A robust understanding of this would enable the government to develop holistic mitigation and control measures to protect life and properties in the area.

Najran Fissure Monitoring Recommendation
So many arithmetical, theoretic, geomechanical, and physics models have earlier been established to explore the configurations, tensile forces, and resumption of fissuring predicated on underground water lateral motion, drainage, agrarian operations, subsidence, continuous soil cracking, bedrock characteristics, and aquifer/aquitard layer thickness, as evidenced by the ongoing discussions in the literature [29,[81][82][83][84][85]. We have shown that the factors responsible for fissuring in Najran are multifaceted and may be more serious than previously thought. If all the probable sources of hazard are identified, they can be combined to create more robust conceptual and/or predictive models of the fissure development mechanism in the area. It is especially important to have a holistic model of the earth fissure hazards in the area to recognize susceptible groundwater areas for the proper enactment of water management policies and efficacious implementation of groundwater recharge strategies to achieve environmentally friendly preservation plans and protect existing groundwater resources [86]. Predicting and comprehending hazardous locations can assist policymakers in managing susceptible areas and reducing the likelihood of damage. The findings would also be useful for water resource managers in formulating judgments on how to properly regulate groundwater extraction. To investigate the creation mechanism of fissuring in Najran, statistically based models such as time series models or artificial intelligence (AI) can be applied. These models would primarily be used to monitor active fissures and identify key factors, utilizing spatial information, subsidence, groundwater-related data sets, etc. [85][86][87][88][89]. As discussed by the authors of [90] for landslide prediction, metaheuristic algorithms such as artificial bee colony (ABC), genetic algorithms (GAs), [90][91][92][93], gray wolf optimization (GWO) algorithms [94,95], particle swarm optimization (PSO) algorithms [96,97], and water cycle algorithms (WCAs) [97] can also be implemented for fissure predictions. Ref. [90] further demonstrates that a combined application integrating k-fold cross-validation, metaheuristic support vector regression, and nonparametric Friedman tests can be utilized to determine optimum hyperparameters and present their statistical significance, thereby improving the accuracy and dependability of AI-based forecasts, and with the help of these methods, different machine-learningbased geo-hazard models might be analyzed and contrasted to ensure high prediction precision [98].

Future Fissure Studies in Najran
It is suggested that fissure-related studies in Najran and the surrounding cities focus on identifying all the major potential contributors to their formation processes, as well as assessing the contribution of an ephemeral variability of predictors such as rainfall, groundwater drawdown, groundwater variability, and so on, in order to develop a reliable model that can be utilized to forecast and monitor fissure formation in the area. This would be a great opportunity to learn more about the earth fissures. The primary determinants in fissure formation can be revealed through a sensitivity analysis of the predictors. Future findings could serve as a foundation for fissure research in other sections of Saudi Arabia and beyond. The increased demand for the prediction of possible hazards and susceptibility mapping can be met by constructing a credible fissuring model [7,[99][100][101][102]. Irregular bedrock topography, low groundwater recharge, excessive aquifer depletion, water table decline, a high density of groundwater wells, a high density of constructed roads, collapsible sediment distribution, fault zones, and other factors could make these areas more vulnerable to fissure hazards, and an accurate predictive model could help the government monitor these areas.

Conclusions
The factors that cause fissures are complex and synergistic. The investigation of this danger in many areas has revealed that multiple elements contribute to its creation. Because of the stress imbalance caused by groundwater extraction, pre-existing geologic features such as faults, fractures, joints, or sloppy terrain are mostly activated. The likelihood of pre-existing underground faults being activated by excessive groundwater extraction is quite likely, given the fracture history of the underlying geology of the city of Najran. Because prior studies of the fissures in this area failed to account for pre-existing faults or channels, future research should concentrate on this and other potential variables. Moreover, additional research might focus on finding all probable fissure predictors in this area and constructing holistic, realistic, and efficient predictive models that could be used to monitor areas with significant fissuring potentials. Such models would aid the government in making critical choices on facility development, city planning, agricultural expansion, and groundwater extraction regulation for both irrigation and domestic use.