Liquefaction in Palu: the cause of massive mudflows

The 7.5 Mw tectonic earthquake that hit Palu City on 28 September 2018 was followed by tsunami and liquefaction, triggered massive mudflows in Balaroa, Petobo, and Jono Oge areas. This study focuses on the generating factors of liquefaction such as the condition of soil lithology, depth of water table, the distance to the focal mechanism, and the thickness of soft sediment. Microtremor data, including the Horizontal Vertical Spectral Ratio (HVSR), geological condition, and borehole data, were examined to conduct the liquefaction analysis. The analysis results based on the microtremor data showed that the distribution of ground shear strain values in Palu City ranged from 0.75 × 10–4 to 2.56 × 10–4. The distribution of the locations of the liquefaction was correlated to the distribution of ground shear strain values. High ground shear strain values and a shallow groundwater level were discovered in Palu City valley, which indicates that liquefaction in Palu City will undoubtedly occur. The semi-empirical method confirmed that Balaroa, Petobo, and Jono Oge had undergone large-scale liquefaction at a maximum depth of 16 m below the ground level. The average peak of water runoff that generated the mudflow was estimated to be at 11.31 cm3/s. Since the soil has loose soil grain with high water content, the soil will turn into a massive amount of mud during the liquefaction.


Introduction
On September 28, 2018, an earthquake of M w 7.5 struck Central Sulawesi Province, Indonesia with the earthquake hypocenter located at a depth of 11 km and the geographic coordinate of 0.18°S and 119.85°E. the earthquake had prompted extensive damages ( Fig. 1) and mudflows in Palu City Valkaniotis et al. (2018), Mason et al. (2019). The mudflow covered the areas of Balaroa, Petobo, and Jono Oge, as shown in Fig. 1 b-d, respectively. The enormous life loss due to the liquefaction in Palu could probably be attributed to the three mudflows, which moved material within hundreds of meters downslope parallelly. Two of these mudflows originating as lateral spreads directly beneath the Gumbasa irrigation channel occurred at Petobo and Jono Oge (Fig. 1c,  d). Meanwhile, the mudflow in Balaroa was located close to the Palu Koro fault line, and it was interpreted that it fell along the fault as it was close to the fault line. The mudflow left severe damages on a large part of the soil on the gentle slopes of Palu City Socquet et al. (2019) within a distance of 1 km approximately. The damage of the Gumbasa irrigation channel triggered flash floods which turned into wave-like mudflows (Mason et al. 2019). The Gumbasa irrigation channel was an anthropogenic hazard of unexpected landslides in Palu that clearly demonstrated the need for more proactive assessment (Bradley et al. 2019).
The epicentral distance to the liquefaction site in Palu City was plotted in an image developed (Ambraseys 1988) based on the database regarding shallow earthquakes around the worldwide (Fig. 2). As shown in Fig. 2, the distance between Palu City and the epicenter of the 7.5 M w earthquake was 80 km, which indicated consistency with the findings.
The Meteorological, Climatological, and Geophysical Agency (BMKG) reported that the earthquake that hit Jalil et al. Geoenviron Disasters (2021) 8:21 the city and coastal area caused enormous damages to the settlements and buildings (Sadly, 2018). This circumstance also developed liquefaction, landslide, and mudflows Mason et al. (2019). Thus, it is required to carry out liquefaction and mudflows research to identify and reduce the disaster risk.
The Palu City soil investigation report mentioned that the area was composed of alluvium deposits with a sand  (Widyaningrum 2012). Widyaningrum (2012) also conducted the semi-empirical and liquefaction potential index (LPI) methods with the results reported a shallow groundwater level in Palu City. Although this study provides valuable information regarding the liquefaction susceptibility in the Palu area, the researchers only used semi-empirical and geological procedures. Hence, the effect of liquefaction that causes mudflow was not explained in this study. In general, previous studies analyzed the liquefaction potential based on the semi-empirical approach proposed by Idriss and Boulanger (2008). The liquefaction potential index (LPI) was determined by referring to the (Iwasaki et al. 1981) methods. The National Center for Earthquake Study released a collateral hazard map due to the 2018 earthquake in Palu, which included flow liquefaction and mudflow Irsyam et al. (2019).
Some researches on liquefaction studies apply microtremor wave propagation. Futhermore, the researchers also made numerical models of geotechnical effects of differences in seismic parameters (Araujo and Ledezma 2020). This study was carried out by correlating the wave characteristics, soil site classification, and liquefaction evidence during the 2018 earthquake. The objective of this study is to evaluate the soil susceptibility of Palu City to the liquefaction phenomenon, which might cause flash floods that resulted in the massive mudflows and soil dynamic effect. In addition to this, to support the regional infrastructure planning and development. The soil susceptibility resulting from the dynamic effects was analyzed by combining the velocity amplification values obtained from microtremor measurements with the SPT derived from the site classification by using an ArcGIS framework.

Geological condition of Palu
The geology of the Palu area is dominated by sedimentary and volcanic deposits from the Cretaceous and Paleogene (Fig. 3). Palu contains five formations: Formation of the Alluvium and Coastal Deposits, Formation of the Celebes Molasse, Formation of Granite, Metamorphic Complex, and Formation of Tinombo. The Geological Survey of Indonesia investigated the formations in the 1970s by Sukamto (1973).
The Palu City area is a quarter sedimentary expanse. Silty clay, silt alluvial sand and old river channel deposits were found along the Palu river at the lower elevations near the valley's center. Young and old alluvium fan deposits are surficial sediments in the low relief hills, which extend on the west and east sides of the valley. Deposits of colluvium debris from gravel sand often occur at higher altitudes Thein et al. (2014).
Balaroa morphology is approximately the similar to Petobo, although Balaroa is located near the western hills where the slope is more considerate than Petobo. There is no irrigation channel in Balaroa, as opposed to another area affected by soil flow. The most likely factor that triggered significant liquefaction in the region is the influence of the colluvial fan topographic, which is controlled by natural shallow groundwater Faris et al. (2019). Petobo and Jono Oge are located on a gentle slope in the eastern  (Sukamto 1973) hills. The morphology of the area is a colluvial deposition area originating from the hill.

Horizontal to vertical spectrum ratio (HVSR) method
Data resulted from the microtremor measurement was analyzed using the Horizontal to Vertical Spectral Ratio (HVSR) method. Subsequently, a comparison was carried out between H (horizontal) and V (vertical) spectrum data. The HVSR method was introduced by Nakamura (1989), to estimate the resonance frequency and local sediment amplification factor of microtremor data. Moreover, it is adopted to estimate the soil susceptibility index (Nakamura 1997) and the building susceptibility index (Sato et al. 2008). The critical parameters generated from the HVSR method are the resonance frequency and amplification factor. The HVSR method carried out for sediment data measurements should meet the criteria recommended by SESAME (2004).
The dominant frequency and amplification factor of the site can be determined by using the H/V ratio. Nakamura (1997) also proposed the ground shear strain (GSS) to be used as a mechanism to determine the level of liquefaction by comparing the value of ground shear strain (GSS) and seismic susceptibility index (K g ) to determine the location of the soil movement in the 2018 liquefaction. The value of GSS and K g can be used to determine the location of vulnerable liquefaction points. A point with a higher K g value indicates higher liquefaction potential.
The GSS values are determined from the soil structure strains, which can be defined as.
where γ = ground shear strain, α = peak ground acceleration at bedrock. The value of K g can be used to verify the predicted liquefaction points by comparing them to the damage resulted from the 2018 earthquake. The K g value is obtained by squaring the HVSR spectrum with the peak value of the resonant frequency, which is determined by Nakamura (2000) as: where K g = the seismic susceptibility index, A g = amplification factor of the HVSR spectral peaks, and f g = natural frequency of the ground surface (Hz).
In this study, the evidence of Balaroa liquefaction events was compared by applying the (Idriss and Boulanger 2008) method, which is more suitable for comparing the theory and field results of the liquefaction events. Many researchers have widely used this method, such as Mase et al. (2020), Jalil et al. (2020). The safety factor against liquefaction is determined by comparing the cyclic resistance ratio (CRR ) and cyclic stress ratio (CSR), which is expressed as: where FS liq = the safety factor against liquefaction. The values were predicted based on the N-SPT penetration resistance measured in various depths.

Field investigation
This paper uses the results of a microtremor survey in 150 measurement stations carried out by Thien (2015). The measuring instruments include the PIC, Portable Intelligent Collector, microtremor measurement device with three sensor components, and a data recording unit.
Microtremor was repeatedly recorded in 40.96 s (4,096 data in 100 Hz sampling) at each measured location, and 10.24 s of data were selected from the viewpoint of less artificial noise (Sato et al. 2008). Then the selected data was transformed into Fourier before the horizontal to vertical spectrum ratio was calculated for each component of each measurement. Finally, the H/V spectrum ratio was derived as the average spectrum ratio (Nakamura, 1989). The dominant frequency and amplification factor were obtained from the H/V spectral ratio. This procedure was performed in 1990. The data microtremor survey is presented in Fig. 3, and the seven boreholes (BH) used for the analysis are shown in Fig. 4 which are located in the Balaroa village and Palu coastal area. The boreholes BH-1 and BH-2 were at a depth of 20 m, and BH-3 to BH-7 were at a depth of 30 m below the ground surface. Figure 5 shows the ground slope plan-view from Balaroa village to the coastal area. The loose sand layer was found on the surface between 1 and 15 m depth below the ground surface, whereas a relatively rigid layer was found in the borehole BH-7. At a depth of 15 and 25 m, the silty gravel layer was found with an N-SPT value of 50. The average N-SPT value was higher than 50 at a depth of 25 and 30 m, indicating a dense layer. The dominant sedimentary soil that stretches in the Palu basin consists of gravely sand, sandy silt, sand, silt, sand and gravel, and silty gravel.
The groundwater level (GWL) in the command area was between 0.3 and 10 m below the ground surface. Near the coast of Palu City, the groundwater level was found at a depth of 0.3 m. While in the hills on the west of Balaroa, the groundwater was found at a depth of 10 m (Fig. 5). The boreholes BH-1 and BH-2 were located near the coastal area of Palu City, and BH-3 to BH-7 were on Balaroa gentle slope. A higher SPT value was found in BH-7 and BH-2, which were located on Balaroa steep slope and near the coastal area, respectively. The soil resistance of BH-3, BH-4, BH-5, and BH-6 was relatively weak in comparison to BH-7 and BH-2. The variation of (N 1 ) 60 values is, as shown in Fig. 5. The groundwater level in the Palu coast and Balaroa gentle slope was 0.3 m and at a depth of 10 m, respectively. The groundwater level in the Palu basin was shallow, with a deposit soil layer containing loose sandy soil, a confined aquifer, and low-level permeability (JICA 2019).

Flow liquefaction analysis
The flow liquefaction phenomenon may take place in either saturated loose sands, silts, or very sensitive clays as it is associated with the contractive behavior of loose sand materials to become liquefaction which possibly occurs during a strong earthquake (JICA, 2019). Li and Ming (2000) carried out a performance model to explain the liquefaction loading conditions. The simulations for undrained triaxial compression were established from anisotropic stress conditions. Laboratory test data and models which calculate the flow liquefaction and cyclic mobility illustrate the behavior of saturated loose sand, dense sand, and undrained soils during the shear loading (monotonic and repetitive), as have been published by several researchers (Robertson, 2010). Furthermore, by using SPT penetration resistance data, flow liquefaction can be predicted.
The corrected SPT blow counts (N 1 ) 60 shown in Fig. 5 is employed to predict flow liquefaction. The procedure to evaluate the contractive state susceptibility of the soil to flow liquefaction based on the relationship between yield strength ratio and penetration resistance was proposed by Olson and Stark (2003). The yield strength ratio in Fig. 9 presents the penetration resistance (N 1 ) 60 value that is limited to be less than 12 (Eq. 4). The boundary of the relationship between yield strength ratio and penetration resistance are described as: is yield strength ratio, (N 1 ) 60 is penetration resistance.
The liquefaction susceptibility to flow is highly dependent on the initial state of the soil. Fear and Robertson (1995) proposed a susceptibility boundary relationship employing the concept of a critical state soil mechanics and laboratory test results. Fear and Robertson (1995) recommended the boundary for practices that describe whether a field is susceptible to a flow failure based on the compatibility to the theory, laboratory test results, and field case history. This relationship is shown in Fig. 10. The recommended boundary relations can be expressed below:

Flash flood caused mudflows
Rapid landmass movement on gentle slopes is influenced by excess pore water pressure when the earthquake shocks trigger the water pressure in the aquifer to emerge to the ground surface. The generation of pore water pressures on loss soil shear strength has changed into the debris flow wave on the ground surface. Because there are no data available to reveal the peak runoff during the earthquake, this phenomenon is assumed to be a hypothesis similar to flash floods due to high rainfall events in the river basins. Rainfall-induced landslides are one of the disasters in complex areas of terrain, especially regions with heavy rainfall regularly. Landslides can occur in a catchment area of unconsolidated soil prone to landslides, as well as flash floods that occur caused by the catchment area broken watersheds. The landslides that occurred in Balaroa, Petobo, and Jono

Fig. 4 Locations of the boreholes in Balaroa area
Oge were caused by a broken catchment area after the Palu earthquake. The gentle slope had fissures/cracks in the catchment area of Balaroa, Petobo, and Jono Oge by 0.4 km 2 , 1.43 km 2 , and 1.35 km 2 , respectively (Mason et al. 2019). It is clear that flash flood is the primary cause of infrastructure and casualties during earthquakes in gentle slopes. Rainfalls flow from the catchment area into low-level slope areas. The surface runoff is referred to hydrograph unit and used to design the drainage system. The empirical methods were performed to represent the complex relationship between rainfall and peak surface runoff. The rational methods were conducted to analyze the runoff in a relatively small catchment area. The runoff is expressed as: where Q = the maximum runoff in (cm 3 /s), i = the designed mean intensity of rainfall (m/s), A = area of (6) Q = KiA catchment (m 2 ), K = surface runoff coefficient, for slope K is 1 (Geotechnical Control Office 1984).
Equation 6 is generally used to calculate the runoff, thus the water spray from the aquifer layer can be assume as a runoff. Furthermore, i is the mean intensity of rainfall, which is assume to the velocity of water flowing on the sand as proposed by Darcy's law as follows: By applying the Darcy equation, Eq. (6) becomes where K = coefficient of surface runoff, v = velocity (m/ sec), Δh = different in height (cm), k = coefficient of permeability (cm/s), l = distance (cm).

Ground shear strain (GSS) and soil deformation
The value of the ground shear strain is the parameter used to represent the soil strain in the event of earthquake. The higher the soil shear strain value, the more likely the surface sediment layer can stretch. These strains and shifts will cause deformations, such as soil fractures, landslides, and liquefaction (Nakamura, 1989). Conversely, the smaller the GSS value, the lesser soil layer deformation is likely to occur. Ground shearstrain at the surface can be obtained by using Eq. 1. The PGA values at the bedrock are within the range 100-550 gals (Thein et al. 2014). Figure 6 shows the distribution of ground shearstrain values depicted on the command area's Digital Elevation Model (DEM) map. The sufficient value of the shear strain was almost the same as the distribution of the Kg value, where the moderate to high value was found in the Palu valley area and was directed by the Palu Koro fault line. Soil layers will experience plastic conditions if the ground shear strain (GSS) values range from 1000 × 10 -6 , and the ground shear strain (GSS) values range higher than 10,000 × 10 -6 (Ishihara, 1996).
Such soil condition is vulnerable to landslide and liquefaction whenever an earthquake occurs.
The values of GSS show shallow sedimentary layers in Palu City, as presented in Fig. 6. The magnitude of GSS values ranged from 0.75 × 10 -4 to 2.56 × 10 -4 . Thus by having such a range, the soil began to be elastic-plastic that the soil was easy to experience fractures and surface soil settlement which is shown in Table 1 (Ishihara, 1996). In addition, high ground shear strain values were found in the southern, eastern, and western parts of Palu City surrounding Balaroa and following the trail of the Palu Koro fault. These results tend to be similar to GSS (BMKG, 2018). These findings in Balaroa are in line with the 2018 Palu earthquake geotechnical damage facts, especially in soil displacement caused by fault movements resulting in liquefaction and mudflows. Miyajima et al. (2019) also discovered that those fault movements induced large soil displacement in Palu City.
Equation 2 is used to determine the value of the seismic susceptibility index (K g ). The Seismic susceptibility index is also related to geomorphological condition resonant frequency (see Fig. 7). A high value of seismic susceptibility indicates that the area consisted of alluvial material similar to the coastal area. Areas with a high seismic susceptibility index (K g ) were found in the valley and coastal sediment areas. In hilly areas, it indicates a shallow seismic susceptibility index.

The liquefaction analysis under the empirical methods
The evaluation method of liquefaction using the N-SPT data is called the simplified procedure. It is a deterministic method using a seismic load by considering the cyclic stress ratio (CSR). The ability of the soil layer to resist liquefaction is expressed in the cyclic resistance ratio (CRR ). Liquefaction resistance of granular soils is also determined by using the CRR of the liquefaction potential assessment under the simplified shear stress technique (Idriss and Boulanger 2008).
Based on the geological characteristics, the investigation of the subsurface soil layer in the Balaroa areas shows that it was dominated by granular material. Sandy silt, silt and gravel, and dominant silt sand found in the  first 20 m depth were with FC (fines content) in the range of 15-40%. This soil layer was followed by silt and sand in the 30 m depth, with FC ranging from 5 to 45%. The soil classification was determined based on the average N-SPT of the top 30 m of the soil column. The procedure for calculating the average N-SPT follows the Indonesian Standards (SNI 03-1726-2019 2019). It indicates that soil in Balaroa consists of medium soil (SD) and soft soil (SE) classifications. The characteristics of medium soil site (SD) class with N-SPT values 15 to 50 were found at BH-1, BH-2, BH-4, and BH-5 to BH-7 locations. Meanwhile, the soft soil site class (SE) characteristics with N-SPT value ≤ 15 were found at BH-3 and BH-6 locations.
The peak ground acceleration (PGA) on each borehole was determined based on the Indonesian National Standard or SNI 03-1726-2019. PGA was calculated by observing the specific site and considering the amplification factor, including the site coefficient (F PGA ) value. The maximum ground acceleration at the surface (PGA M ) is shown in Table 2.
The maximum PGA value was found in BH-1, which was located close to the coast, while the minimum PGA value was in the hilly area of BH-7. The areas near the coast indicated extreme ground motion due to its maximum soil density, which was stable against landslides. The gentle slope of the Balaroa hilly indicated maximum ground motion, in addition to the loose sand layer that enabled liquefaction to occur. The results of semi-empirical analyses are as shown in the liquefaction potential in Fig. 8. The analysis results prove that the liquefaction that occurred in 2018 could be predictable by using the microtremor seismic susceptibility index data (K g ).
The soil density parameters used in the analysis were the dry and wet soil unit weight of 17.0 kN/m 3 and 19.5 kN/m 3 . The groundwater level influenced the calculation of total stress and not the effective vertical stress. The depth of the groundwater level is presented in Table 2. The results of the potential liquefaction calculation in BH-6 are presented in Fig. 8. Figure 8 shows that the depth of liquefaction was less than 16 m below the ground surface. The findings from the N-SPT and microtremor seismic susceptibility index data indicate a significant factor of the deposited soil layer on the response susceptibility and manifestation of liquefaction damage. These locations experienced the liquefaction during the 2018 earthquake. The analysis of the potential liquefaction of all boreholes is shown on Table 3.

Flow liquefaction
The penetration resistance plotted in Fig. 9 are the values of (N 1 ) 60 which are limited to less than 12 based on the data in Fig. 5. The value of yield strength ratio, as shown on Table 4, is generally referred to as residual strength and the denoted, S u .
A study by Ishihara (2019) showed that the calculation failure of the liquefaction flow was the back-calculation of the residual strength ratio based on the actual case of slide flow occurring in the Palu plains after the 2018 earthquake. The resulted S u values for Palu mainland ranges within 0.05-0.25 which are in accordance with the range values of residual strength ratio of 0.240-0.272 as shown on Table 4.  The nature of loose sand soil fabric with flow behavior can be evaluated based on the penetration resistance data (Fear and Robertson 1995). Olson and Stark (2003) conducted a static stability analysis to predict the migration of the soil particle with residual strength due to an earthquake. The SPT data based on the flow failure susceptibility relationships taken from the case history data proposed by (Fear and Robertson, 1995) was used to carry out the susceptibility analysis of flow liquefaction. The susceptibility flow liquefaction in the soil layer is identified by the soil contractive state (Fig. 10). The evaluation of liquefaction based on data (N 1 ) 60 and the vertical effective stress is plotted in Fig. 10. These plots show the layers of contractive soil and identify the flow failure. The plotting results, as shown in Table 5, indicate that the region experiences flow failure, which is consistent with the 2018 liquefaction occurring in Palu City.
In addition, Robertson et al. (2019) also reported that based on the field monitoring, flow liquefaction taking place in the Balaroa, Petobo, and Jono Oge areas at the slope ground levels were 3.4%, 2%, and 1.4% respectively. Such softening soil moved to a distance up to 1000 m, 2200 m, and 3900 m for the area of Balaroa, Petobo, and Jono Oge, respectively. The soil movement was a distant   Fig. 9 The corrected SPT blow count relationship to estimate liquefied strength ratios and for liquefaction flow failures, as proposed by Olson and Stark (2003). Red arrows indicate the BH-3, BH-4, BH-5, and BH-6 estimation trajectory that buried thousands of houses and other facilities in Balaroa, Petobo, and Jono Oge.

A phenomena of flashflood like: the caused mudflows
Flow liquefaction occurrence in Palu City was similar to flash flood, which considers as rapid flooding of water above the ground caused by a sudden burst of water from an aquifer strip, or the release of retained water (e.g., the rupture of an impermeable confined aquifer). Flow liquefaction or flash flood occurs for a short period, usually within a few minutes to a few hours, distinguishing it from fluvial flooding. The phenomenon of liquefaction in Palu City showed similar behavior to the flash floods event.
The uniqueness of the geological system in the region of Palu determines its hydro-morphological state, including the groundwater basin. On these sediment layers lie Petobo, Balaroa, and Jono Oge that suffered a wide degree of liquefaction-induced-flow failure. This aquifer layer is the Palu groundwater basin, which has a very high average aquifer transmissivity value. The landslide event within the earthquake in Balaroa is illustrated in (Fig. 11).
The groundwater level in the Petobo area was shallow and was affected by the infiltration of the Gumbasa irrigation channel (Fig. 12). The soil layer was very saturated in addition to the presence of a confined aquifer layer. The water flow in the channel enters the porous soil layer through the tear gap resulted from the earthquake shocks. In this condition, the soil turned into pulp and was carried away by the flash flood waves. The soil types in Petobo and Jono Oge are shown in Fig. 13a, b, respectively. The soil layer was permeable and porous so that the sub-surface water pressed out through the gaps and raised the sediment by approximately 5 m (JICA 2019). Such rapid change of stress in the soil made the effective  (Fear and Robertson 1995) Fig. 11 The illustration earthquake induced flow liquefaction in Balaroa vertical soil stress become zero and reached the maximum pore water pressure. Thus, the soil layer induced enormous strain and liquify. The movement of soil masses and rock debris down a slope or the rapid and gradual sinking of the ground surface in a dominant vertical direction are in response to the gravitational force. Ground movement velocity significantly affects the volume of the landmass on the downslope. Soil mass movement velocity strongly depends on groundwater level fluctuation, dynamic viscosity coefficient, and sliding surface mass engineering properties Naing et al. (2018). Considering the impacts and damages caused by the mass movement, it is essential to predict its runout distance, velocity, moving volume, and coverage area Fathani et al. (2017). The mass movement model was introduced by several researchers Fathani et al. (2017), Syah et al. (2019. Small earthquakes of M w 4 can produce landslides, whereas the event M w > 7.5 usually creates a landslide of thousands of landslides on vulnerable slopes in an area of hundred square kilometers with millions of cubic volumes of loose material (Keefer, 1994). Keefer (1994) proposed a relationship between the magnitude of earthquakes and the historical number of landslides (Fig. 14). If Fig. 14 is applied for the M w 7.5 landslide in Palu City, it is found that the volume of landslide material is 2.3 × 108 m 3 , which is considered as an enormous volume of landslides that buried and swept away residential areas in Palu City. However, the volume of the landslide triggered by M w 7.5 Palu earthquake on September 28, 2018, was estimated at 9.5 × 10 6 m 3 , distributed in Balaroa, Petobo, and Jono Oge equivalent with 3 m in the height of debris.
By applying Eqs. (6), (7), and (8), the results of the flash flood discharge and velocity can be obtained, as shown in Table 6. The vertical permeability coefficient used in the calculations is based on the results of a field survey by JICA (2019). Therefore, the peak runoff of water flash floods for the Palu area could be calculated as the averaged value of the three locations from Table 6. The result is 11.31 cm 3 /s. The saturated mudflow material contained water that was equal to or higher than the liquid limit. Mudflow is a viscous liquid state with the criteria within the ratio of the average width and length of the soil mass movement, which ranges from 0.05 to 0.3 (Liu and Mason 2009). Meanwhile, the mudflow ratio in Palu City ranged from 0.2 to 0.4. However, there was some accordance between the (Liu and Mason 2009) criteria and the mudflow event in the Palu region. Meixler (2018) reported that the shaking of the earthquake generated the loose material to become wave-like mudflows and move down the slope was responsible for the massive damage in the Palu area.
The landslide that turned into mudflows in Balaroa was tremendously influenced by a 3-5 m displacement of the total length of the 140 km of Palu Koro fault Valkaniotis et al. (2018). The morphology mapping of landslides in Balaroa was impacted by secondary fault strands, which resulted in erosion. Furthermore, the major mudflow in Balaroa damaged the residential communities of the main downstream. The two landslide flows in Petobo and Jono Oge in the Palu basin were underwent an enormous land deformation. However, the most significant land displacement scale was taken place in Jono Oge village Hazarika et al. (2021).
The displacement pattern of the area after the earthquake was able to be clearly defined. The rice cultivation area experienced displacement of more than 5 m was affected by the Gumbasa irrigation channel damage and the rice planting scheme at the local location, which raised the groundwater level (Bradley et al. 2019). Based on the rational method approach, the average runoff peak at the Balaroa, Petobo, and Jono Oge locations is presented in Table 6.
The limitations of this study in determining the water runoff induced by the flash floods during the Palu earthquake required further research. Groundwater depth position data was used to predict the emergence of runoff peaks which were assumed to be the same in all areas.

Conclusions
This study presents a liquefaction analysis of the 7.5 M w earthquake in Palu, Indonesia, on September 28, 2018. The phenomenon of massive liquefaction led to mudflows in the Balaroa, Petobo, and Jono Oge areas. Microtremor and N-SPT seismic measurements were carried out to determine the soil susceptibility to the dynamic influence of the command area.
Liquefaction is controlled by the geological formation of alluvial deposits that are dominated by gravel sand and muddy sand. The Balaroa liquefaction was the only failure observed in the north to the south side of the Palu valley. The displacement of the Balaroa slope, which was most likely caused by a secondary strand fault movement, triggered mudflow. As Petobo and Jono Oge areas are located on the western slope of the Gumbasa irrigation channel, the water flow on the channel also accelerated the emergence of mudflows through the ground fractures. It can be concluded that the average peak of water runoff causing the mudflow was approximately 11.31 cm 3 /s. The shaking of the 2018 earthquake affected the softening of the soil layer and consequently produced liquefaction in the Balaroa, Petobo, and Jono Oge regions since the areas have shallow groundwater levels. Liquefaction occurence in Balaroa could be confirmed by using the semi-empirical analysis and susceptibility index approach from the microtremor survey. The result indicates that liquefaction can occur at a maximum depth of 16 m below ground level. The depth of liquefaction was found at 16 m below the ground surface. The N-SPT and microtremor seismic susceptibility index data indicate a significant factor of the deposited soil layer on the response susceptibility and manifestation of liquefaction damage.