Next Article in Journal
A Unique Conditions Model for Landslide Susceptibility Mapping
Previous Article in Journal
Geological Assessment of Faults in Digitally Processed Aerial Images within Karst Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and Verification of the Movement of the Hidden Active Fault Using Electrical Resistivity Tomography and Excavation

by
Rungroj Arjwech
1,*,
Sutatcha Hongsresawat
2,
Suriyachai Chaisuriya
3,
Jetsadarat Rattanawannee
1,
Pitsanupong Kanjanapayont
4 and
Winit Youngme
1
1
Department of Geotechnology, Faculty of Technology, Khon Kaen University, Khon Kaen 40002, Thailand
2
Division of Geoscience, School of Interdisciplinary Studies, Mahidol University, Kanchanaburi Campus, Kanchanaburi 71150, Thailand
3
Office of Mineral Resources Region 3, Department of Mineral Resources, Ministry of Natural Resources and Environment, Pathum Thani 10120, Thailand
4
Department of Geology, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
*
Author to whom correspondence should be addressed.
Geosciences 2024, 14(8), 196; https://doi.org/10.3390/geosciences14080196
Submission received: 5 June 2024 / Revised: 18 July 2024 / Accepted: 23 July 2024 / Published: 24 July 2024

Abstract

:
Identifying the movement of the branches of the hidden Thakhek fault in Thailand is challenging due to the absence of evident landforms indicating an active fault. In this study, we analyzed a digital elevation model (DEM) to identify potential landforms. A 2D Electrical Resistivity Tomography (ERT) survey was conducted to locate the hidden Thakhek fault. The results reveal vivid images of resistivity contrast, interpreted as two reverse faults, with mudstone exhibiting low resistivity in the middle, flanked by thick sediment layers with higher resistivity. Three trenches were excavated perpendicular to the two interpreted reverse faults. The displacement of reverse faulting appears to have shifted mudstone over Quaternary sediments, with vertical offsets revealed in trenches NWY-1, NWY-2, and NWY-3. This movement could be identified as a positive flower structure. Additionally, lakes are identified as a negative flower structure along the traces. These features result from strike-slip strains under a locally appropriate compressional and extensional environment in a shearing strike-slip fault.

1. Introduction

The active fault zone along the eastern edge of the Khorat Plateau demonstrates strike-slip movement. This zone reportedly comprises three significant faults oriented NW–SE, showing evidence of neotectonic displacements [1,2]. The fault zone is clearly characterized by a linear or curvilinear principal lateral displacement zone in satellite image investigations [3]. Landforms of strike-slip faults can be identified in satellite images due to significant lateral displacement, such as wrench-fault troughs, pull-apart basins, or a combination of both [4,5]. Wrench-fault troughs can form relatively long and narrow basins [6]. Moreover, they can create a depressed area, such as a sag pond, a rhomb graben, or a pull-apart basin, identified on a larger scale [7].
The Thakhek-Sepon-Danang active fault zone runs predominantly NW–SE [8,9,10] and extends from Laos into Vietnam, spanning over 550 km. This fault zone comprises three major faults, as illustrated in Figure 1 [9,11,12]. During the neotectonic stage, the fault zone experienced two main phases: the early phase (Late Miocene—early Pliocene) with left-lateral slip, and the late phase (Late Pliocene-Quaternary) with right-lateral slip [9]. The fault zone partly runs along the Mekong River at the Thailand–Laos border, characterized by left-lateral slip, and is mapped as a steeply dipping NW-striking fault [13,14]. Fault motion formed the plateau edge on the Lao side, caused by compressional forces resulting in older rocks moving upwards [14]. Many prominent active fault morphologies are visible in satellite images of Laos [1].
Flower structure is one of the important features that can be used to identify strike-slip faults and shear zones [5,15]. The structural features are mainly characterized by their internal fault and fold architectures in cross-sections that resemble the petals of a flower [15,16]. According to their internal fault structural architecture, flower structures can be classified into two main types: positive and negative structures [15,16,17]. A positive flower structure, also termed palm-tree structures, consists of a shallow antiform displaced by upward diverging strands of a wrench fault with mostly reverse separations. A negative flower structure, also called tulip structures, consists of a shallow synform bounded by upward spreading strands of a wrench fault with mostly normal separations [15]. These structures are bounded by strike-slip faults on the two parallel sides [7].
The active faults on the northeast edge of the Khorat Plateau are hidden beneath loose sediments, with fault morphology eroded and quickly covered by new deposits, making fault traces hard to detect. Arjwech et al. [2] studied an active fault by excavating two trenches across a fault line identified using 2D Electrical Resistivity Tomography (ERT). Several studies [1,18,19,20] have successfully used the 2D ERT method to detect subsurface resistivity contrasts in unexposed sediment outcrops. However, a reverse fault, bounded by strike-slip faults on both sides of the compressive zone, was not verified. Thus, this study aims to collect additional data from 2D ERT and trenching to reassess the reverse faults, with implications for understanding the morphology affected by compression and dilation along the strike-slip zone.

2. Geological Setting

Thailand primarily consists of the Indochina and Sibumasu tectonic blocks [21]. The Khorat Plateau and the surrounding region have been significantly impacted by tectonic events linked to the prolonged Himalayan orogeny during the Miocene–Pleistocene, characterized by NW–SE compression resulting in fold and fault structures [22]. This plateau stretches from northeast Thailand, extending northward and eastward across the Mekong River into Laos, within the Indochina tectonic block [23,24]. Major active fault zones in this region predominantly exhibit strike-slip movements with NE–SW or NW–SE orientations [25]. Most active faults in Thailand are located in the North and West, while those in the Northeast are mainly found along the plateau’s edges [1,26]. Clear surface expressions of active faults are primarily visible along the plateau rim in Laos [27].
The Khorat Plateau basin is relatively flat with low, small undulating hills [28]. Most of the basin’s current margins are formed by mountain ranges. The Khorat Plateau is characterized by thick sequences of gently folded Mesozoic terrestrial non-marine sedimentary rocks known as the Khorat Group [29,30]. During the mid-Cretaceous, a deformation phase led to the E–W-trending Phu Phan uplift, which traverses the plateau and divides the Khorat Plateau basin into the Sakon Nakhon basin to the north and the larger Khorat basin to the south [31]. Fractures aligned with the major structural trend of the isolated Phu Sing Mountain range are observable in satellite imagery [32]. Isolated mountain ranges, mainly following a NW trend and bounded by regional faults, exhibit synclinal structures with gently dipping limbs at angles of 5–10° [32,33,34]. The fracture system aligns with the orientation of the mountain ranges, where escarpments likely represent fault-generated structures. However, clear evidence of active faulting morphology is not apparent in the field in Thailand.
The geology of the study area is predominantly covered by two members of the Phu Tok Formation of the Khorat Group. The lowest member, the Nawa member, consists of reddish-brown or reddish-orange claystone and mudstone with low resistance [35]. The Phu Tok Noi is the upper member of the late Cretaceous–early-Tertiary Phu Tok Formation. This upper member primarily features large-scale cross-bedded yellowish-red fine- to medium-grained sandstones, interbedded with reddish-yellow very-fine-grained sandstones and siltstones [26,36,37]. The paleo Mekong River carved the sandstone mountains of the Phu Wua, Phu Sing, and Phu Tok ranges, leaving them lying parallel or subparallel to the flow of the powerful river [38]. Additionally, Quaternary gravel sediments are dispersed as hilly areas along the paleo Mekong River [1].

3. Study Area and Methods

The study area is situated in Bung Kan and Nakhon Phanom Provinces in the far northeast of Thailand, near the northeastern rim of the Khorat Plateau, as illustrated in Figure 2. The Nonwangyiam and Nano sites are characterized by low-lying paddy fields flanked by rolling hills with rubber tree plantations. These rubber plantations are primarily concentrated on gravel hills. Partially exposed weathered reddish-brown mudstone can be observed in excavated ponds within the paddy fields. Swamps or lakes are found along the fault trace to the northwest and southeast of the excavation site.
In Thailand, morphological features of active faults are difficult to identify. The ERT survey is ideal for this study because it provides high resolution and wide cross-sections. The data can be processed and interpreted immediately on-site, enabling quick decision-making. In the previous study, high resolution 2D ERT surveys were successful in determining the location of hidden faults in the Thakhek-Bung Kan 1 (TB1) trace. Two reverse faults (F1 and F2) were interpreted in the BK01 tomogram, but only one fault was verified by trench Nonwangyiam 1 (NWY-1), a parallel bounded fault (F2) in the BK01 tomogram was not verified yet. One reverse fault was interpreted in the BK02 tomogram. The fault was verified by trench NWY-2. Two trenches (NWY-1 and NWY-2) were excavated across the interpreted faults F1.
In this study, we analyzed a digital elevation model (DEM) and the most recent high-resolution satellite images available on Google Earth to identify morphotectonic landforms. QGIS software (open-source geographic information system (GIS) application version 3.28) was employed for the analysis of active fault traces using a 30 m resolution DEM. We revisited the Nonwangyiam site and extended the Thakhek fault analysis with one more excavation (Trench NWY-3) and interpretations. Subsequently, trench-log stratigraphy was conducted, and sediments were sampled for dating to determine the ages of sedimentary layers in the trench. Finally, paleo-seismology was studied (though not mentioned in this paper). More than 40 ERT profiles, each consisting of an electrode array with length ranging from 200 to 750 m, were conducted across the interpreted fault, the TB2 trace. We present only the 2D ERT profile BK03 from the Nano site, which displayed similar morphological features to the Nonwangyiam site. This profile revealed the most significant fault-related anomaly.

4. Results

Field investigations, including geomorphological and geological reconnaissance and ERT surveys, were conducted after DEM interpretation. Distinct fault traces with a NW–SE orientation are easily seen as sharp lineaments in satellite images and DEM on the Laos side (Figure 1). These features are evident along escarpments of Mesozoic rock-resistant beds and at the boundaries between Paleozoic and Mesozoic units (Figure 2). However, the branching faults that are clearly defined and prominent on the Laos side quickly fade into unconsolidated Quaternary sediments in Thailand. The surface traces continuing southeast from the Laos side likely represent a distinct fault. These lineaments align with natural creeks, lakes, and swamps on the Thailand side (Figure 1, Figure 2 and Figure 5c). Two parallel or subparallel lineaments, TB1 and TB2, intersect the Late-Cretaceous–Early-Tertiary Phu Tok mudstone and unconsolidated Quaternary sediments.

4.1. ERT Result

The ERT survey was performed using a dipole–dipole configuration with 0.75 and 2.5 m electrode spacing. The acquisition dataset was optimized to take full advantage of the available 10 channels in the SYSCAL PRO resistivity meter from IRIS Instruments (www.iris-instruments.com). The data were processed using 2D tomographic reconstruction methods employing the finite-element modeling inversion algorithms of Res2DInv software (version 4.08). The L1-norm (blocky optimization) inversion was applied, minimizing the sum of absolute deviations (residuals) and resulting in blocky profiles with vertical gradients. This geometric approach provides clear boundaries that separate low- and high-resistivity layers, thereby highlighting significant resistivity contrasts in subsurface geology. For dip-slip faults, a distinct vertical boundary was detected, emphasizing vertical flatness over horizontal flatness. As a result, robust model inversion constraints were found to be more appropriate. The ERT data for profiles BK01 and BK02 were acquired from previous studies [1,2], and the data were reprocessed. The results offer clear details for the excavation of trench NWY-3 and its geometrical features. Additionally, the BK03 profile was conducted, and its results supported these findings.
At the Nonwangyiam site in Bung Kan Province (Figure 5a–c), two profiles, BK01 and BK02, were conducted across the fault trace TB1. Each survey profile measured 235.0 m and 70.5 m, respectively, traversing low-lying paddy fields flanked by gravel hills. The 2D ERT tomograms converged with absolute errors of 2.6 and 1.9 after six iterations (Figure 3 and Figure 4), indicating a good fit for the model inversions. The maximum penetration depths are approximately 25 m and 16 m.
The BK01 tomogram shows two sharp vertical sections (Figure 3a,b). A heterogeneous zone of relatively low resistivity (<50 Ωm) is evident in the middle of the line between 65 and 142 m, interpreted as mudstone. This zone is flanked on both sides by sharply and nearly vertical higher resistivity (>50 Ωm), interpreted as Quaternary sediments throughout the entire depth of the tomogram. The ERT tomogram clearly reveals the dip fault plane, and the contact between rock and sediments is identified as putative faults F1 at 65 m and F2 at 142 m along the survey line. The trench sites in this study are determined by anomalies F1 and F2.
The shorter profile BK02 was deployed across the interpreted F1 fault in profile BK01, approximately 150 m away from BK01. The tomogram clearly shows a sharp vertical resistivity contrast at 22 m, interpreted as the fault. The low resistivity zone is identified as mudstone, while the higher resistivity is attributed to Quaternary sands, silts, and gravel clasts, consistent with the interpretation of profile BK01.
At the Nano site in Nakhon Phanom Province, profile BK03 was conducted across the fault trace TB2 (Figure 5b). The survey profile measured 355.0 m, traversing a low-lying paddy field flanked by hills. Six iterations yielded a final absolute error of 2.5 (Figure 4). The maximum penetration depth is approximately 40 m. The tomogram shows sharp vertical sections. A relatively high resistivity zone (>50 Ωm) on the northeast side is flanked by sharply and nearly vertical low resistivity (<50 Ωm), interpreted as Quaternary sediments. The contact between rock and sediments is identified as the putative fault F2 at 270 m along the survey line. Consequently, the locations of faults are indicated with dashed lines.

4.2. Trench Excavations and Stratigraphy

At the Nonwangyiam site (Figure 5c), Trenches NWY-1 (Figure 5d and Figure 6a,b) and NWY-2 (Figure 5e and Figure 6c,d) were previously excavated perpendicular to anomaly F2. In this study, Trench NWY-3 (Figure 5f and Figure 7a,b) was excavated perpendicular to anomaly F1. The trench measures approximately 20 m long, 4 m wide, and 5 m deep. A 1 m2 string grid was used for the entire trench, following standard paleo-seismic methods, and a lithologic log was carried out. Three depositional units (A–C) were exposed in both sidewalls. Unit A, the oldest, consists of clast-supported gravel with clay, indicative of a fluvial deposit. The gravels are typically subrounded to rounded cobbles of sandstone and quartz, with laterite accretions. Unit B is a fluvio-lacustrine deposit containing a few gravels within a clayey sand matrix. Unit C is a fluvial deposit of clast-supported gravel with clay, with laterite formed at the bottom. The succession is capped with a layer of dark topsoil composed of sand and clay.
These units are displaced by dip-slip faults F11–F14. Reverse faulting appears to have shifted reddish, pale, blown-weathered mudstone and cut through the upper sediment layers. The vertical offset along the faults observed in the trench exceeds 2 m. The fault cuts through the strata, with the hanging wall of mudstone pushed up, cutting through the upper sediment layers. We interpret the movement within the trenches as compressional deformation of a positive flower structure produced by sinistral strike-slip movement.

5. Discussion

Two main subparallel lineaments, TB1 and TB2, trending NW–SE, were previously identified (Figure 1). Fault trace TB1 cuts through low-relief terrain to the floodplain, ending at the Song Kham River, while TB2 runs through varying terrain, passing through Bung Kong Long Lake and joining the Mekong River in the southeast. The Thakhek fault system extends into Thailand, influencing the geomorphological evolution of the Mekong River, where faulting controls flow diversions. Gravel hills, elongated remnants aligned with the interpreted faults, serve as geomorphological indicators of sediment transport processes of the river. Many boulder-strewn areas are currently used as rubber tree farms.
The combination of shear and compression due to transgression could explain the lineaments relative to the regional horizontal stress field observed in satellite images. Evidence includes R-shears along the Phu Sing Mountain ridge and the crescent-shaped Kud Ting Lake. The stress field causes compressional strain on the lower mudstone, accumulating vertical displacement over time relative to the overlying sediment layers. This structure likely corresponds to splays in the main Thakhek fault, developing local reverse faults connected at depth to form a flower structure. ERT geophysical surveys successfully identified these flower structures (Figure 3a). The trench log shows good spatial correlation with the ERT images, displaying large vertical offsets in the trench wall (Figure 6a,b and Figure 7a,b). Vertical deformation associated with the sinistral strike-slip motion caused mudstone on either side of the fault to laterally override younger sediment layers, as seen in Trench NWY-2 (Figure 6c,d). However, the 2D ERT tomogram only images the upper parts of the complex flower structure.
In wrench fault zones, both positive and negative flower structures can form. In this study, a positive flower structure was identified by the 2D ERT method and verified by excavation along the TB1 trace of the Thakhek fault. The general structure includes the predominant crescent-shaped Bung Kong Long Lake located in the TB2 trace, indicating a large pull-apart basin, interpreted as a negative flower structure in a divergent wrench fault zone under transtensional deformation. Divergent-wrench fault landforms are widespread, exemplified by Kud Ting Lake, Nong Bang Bat Lake, Nong Chaiya Wan Lake, and Nong Bung Khra Lake along lineament TB1. These consist of a shallow synform bounded by upward-spreading strands of a wrench fault with mostly normal separations. However, positive flower structures’ landforms could be eroded quickly in tropical rainforest regions and rapidly covered by new deposits, making them less evident at the Earth’s surface. Therefore, sinistral or dextral movements of these branch faults cannot be identified in this study.

6. Conclusions

High-resolution 2D ERT tomograms were used to identify hidden faults along the transpressional faults TB1 and TB2. The comparison of ERT profiles with the excavated trench logs confirmed the success of the ERT survey in detecting these faults. From our integrated study, we interpret the steeply inclined faults in the 2D ERT tomogram as a positive flower structure. The middle section of the tomogram shows a low-resistivity zone corresponding to a mudstone layer with sharp contacts flanked by higher-resistivity sediment layers. Excavations of trenches NWY-1 and NWY-3 confirmed faults F1 and F2, further indicating a positive flower structure. This fault was caused by compression from strike-slip fault movement, as observed in Trench NWY-2. We extended the analysis of strike-slip fault movement using DEM and satellite images, which revealed the sagging of the large crescent-shaped Kud Ting Lake along the TB2 trace. Numerous large lakes, extending over approximately 100 km, are found along the lengths of the TB1 and TB2 traces.

Author Contributions

Conceptualization, R.A. and S.H.; methodology, R.A., S.H., J.R., S.C. and W.Y.; software, R.A.; validation, R.A., S.H. and P.K.; formal analysis, R.A.; investigation, R.A., S.H., J.R., S.C. and W.Y.; resources, R.A. and S.H.; data curation, R.A.; writing—original draft preparation, R.A. and S.H.; writing—review and editing, R.A. and P.K.; visualization, R.A.; supervision, R.A.; project administration, R.A.; funding acquisition, R.A. and S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by the Earth Science Research Cluster, Mahidol University, Thailand Science research and Innovation Fund Chulalongkorn University, and the Electricity Generating Authority of Thailand (EGAT) contract number 60-E303000-11-IO.

Data Availability Statement

The datasets presented here are available from the corresponding author upon request.

Acknowledgments

We would like to thank the landlord, Nay Noo, for allowing us to excavate trenches on his property. We are also grateful to Sriboy, K-win, Tithiwarada, Sakhon, Romyupa, and Pakawat, who assisted with the ERT data acquisition and trench works. Fieldwork was conducted using equipment and vehicles provided by the Department of Geotechnology, Khon Kaen University.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Arjwech, R.; Everett, M.E.; Chaisuriya, S.; Youngme, W.; Rattanawannee, J.; Saengchomphu, S.; Thitimakorn, T.; Somchat, K. Electrical resistivity tomographic detection of the hidden thakek fault, Northeast Thailand. Near Surf. Geophys. 2021, 19, 489–501. [Google Scholar] [CrossRef]
  2. Arjwech, R.; Chaisuriya, S.; Rattanawannee, J.; Youngme, W.; Thitimakorn, T.; Srikraiwest, R. Recent paleoseismic investigations at the hidden Thakhek fault in Thailand. J. Asian Earth Sci. 2022, 236, 105315. [Google Scholar] [CrossRef]
  3. DMR. Active Fault Map of Thailand. 2006. Available online: https://www.researchgate.net/figure/Map-of-active-faults-in-Thailand-modified-DMR-2006_fig1_334469653 (accessed on 1 October 2023).
  4. Noda, A. Strike-slip basin–its configuration and sedimentary facies. In Mechanism of Sedimentary Basin Formation; Intech: Houston, TX, USA, 2013; pp. 27–57. [Google Scholar]
  5. Huang, L.; Liu, C.Y. Three types of flower structures in a divergent-wrench fault zone. J. Geophys. Res. Solid Earth 2017, 122, 10–478. [Google Scholar] [CrossRef]
  6. Christie-Blick, N.; Biddle, K.T. Deformation and basin formation along strike-slip faults. Soc. Econ. Paleontol. Mineral. Spec. Publ. 1985, 37, 1–34. [Google Scholar]
  7. Azizi, R.; Ghannem, N.; Mahmoudi, N.; Chihi, L.; Regaya, K. Active negative flower structure at the western edge of the Tunisian Atlas: Morphotectonic evidence from the case study of Borj Edouane-El Gara Quaternary basin. Quat. Int. 2021, 604, 93–112. [Google Scholar] [CrossRef]
  8. Raksaskulwong, L.; Assavapatchara, S.; Bamroongsong, P.; Khaowiset, K.; Krawchan, W.; Sumart, O.; Keohavong, C. Continuation of the Mesozoic continental deposits across Thailand and Lao PDR. In Proceedings of the Thai-Lao Technical Conference on Geology and Mineral Resources, Department of Mineral Resources Bangkok, Bangkok, Thailand, 7–8 September 2010; pp. 7–8. [Google Scholar]
  9. Van Thom, B.; Van Hung, N.; Sunlinthone, O.; Duangpaseuth, S.; Markvilay, B. Characteristics of the Thakhet-Sepon active fault zone. Vietnam. J. Earth Sci. 2015, 37, 36–47. [Google Scholar]
  10. Zhang, J.; Feng, Q.; Zhang, Z. Tracing escaping structure in the Northern Indo-China Peninsula by Openness and remote sensing. J. Earth Sci. 2017, 28, 147–160. [Google Scholar] [CrossRef]
  11. Duong, N.A.; Van Chinh, V. Some geomorphic indices in the North Central Vietnam. Geosci. J. 2021, 25, 813–829. [Google Scholar] [CrossRef]
  12. Xuan, P.T.; Duong, N.A.; Van Chinh, V.; Dang, P.T.; Qua, N.X.; Van Pho, N. Soil gas radon measurement for identifying active faults in Thua Thien Hue (Vietnam). J. Geosci. Environ. Prot. 2020, 8, 44. [Google Scholar] [CrossRef]
  13. Morley, C.K.; Charusiri, P.; Watkinson, I.M.; Ridd, M.F. Structural geology of Thailand during the Cenozoic. In The Geology of Thailand; Geological Society of London: London, UK, 2011; pp. 273–334. [Google Scholar]
  14. Sattayarak, N. Phu Tok Formation; before and after Tok. Geosocial KKU 2020, 1, 21–25. [Google Scholar]
  15. Harding, T.P. Seismic characteristics and identification of negative flower structures, positive flower structures, and positive structural inversion. AAPG Bull. 1985, 69, 582–600. [Google Scholar]
  16. Harding, T.P. Identification of wrench faults using subsurface structural data: Criteria and pitfalls. AAPG Bull. 1990, 74, 1590–1609. [Google Scholar]
  17. Holdsworth, R.E.; Strachan, R.A.; Dewey, J.F. Continental Transpressional and Transtensional Tectonics; Special Publication: London, UK, 1998; Volume 135, pp. 1–14. [Google Scholar]
  18. Arjwech, R.; Phothaworn, T.; Chaisuriya, S.; Thitimakorn, T.; Pondthai, P. Evaluation of Slope Susceptibility Using 2D Electrical Resistivity Tomography Supplemented with Spatial Resistivity Change. Geotech. Geol. Eng. 2023, 41, 4023–4039. [Google Scholar] [CrossRef]
  19. Arjwech, R.; Everett, M.E.; Saengchomphu, S.; Somchat, K.; Pondthai, P. Geophysical mapping of gypsum for exploration of reserves in the Nong Bua area of Thailand. Q. J. Eng. Geol. Hydrogeol. 2021, 54, qjegh2019-180. [Google Scholar] [CrossRef]
  20. Sarntima, T.; Arjwech, R.; Everett, M.E. Geophysical mapping of shallow rock salt at Borabue, Northeast Thailand. Near Surf. Geophys. 2019, 17, 403–416. [Google Scholar] [CrossRef]
  21. Ridd, M.F.; Barber, A.J.; Crow, M.J. The Geology of Thailand; Geological Society of London: London, UK, 2011. [Google Scholar]
  22. Harnpattanapanich, T.; Luddakul, A. Seismic Hazard of the Khorat Plateau: Preliminary Review. In International Conference on Geology, Geotechnology and Mineral Resources of Indochina (GEOINDO 2011), Thailand. 2011. Available online: http://home.kku.ac.th/geoindo2011/A6-457-502.pdf (accessed on 1 October 2023).
  23. Hite, R.J.; Japakasetr, T. Potash deposits of the Khorat plateau, Thailand and Laos. Econ. Geol. 1979, 74, 448–458. [Google Scholar] [CrossRef]
  24. El Tabakh, M.; Utha-Aroon, C.; Schreiber, B.C. Sedimentology of the Cretaceous Maha Sarakham evaporites in the Khorat Plateau of northeastern Thailand. Sediment. Geol. 1999, 123, 31–62. [Google Scholar] [CrossRef]
  25. DMR. Atlas Book, Active Faults in Thailand; Department of Mineral Resources, Ministry of Natural Resources and Environment: Bangkok, Thailand, 2018. [Google Scholar]
  26. Imsamut, S. Magnetostratigraphy of the Phu Thok Formation at Phu Thok and Phu Wua areas, Changwat Nong Khai. Master’s Thesis, Chulalongkorn University, Bangkok, Thailand, 1996; 306p. [Google Scholar]
  27. Pailoplee, S.; Charusiri, P. Analyses of seismic activities and hazards in Laos: A seismicity approach. TAO Terr. Atmos. Ocean. Sci. 2017, 28, 8. [Google Scholar] [CrossRef]
  28. Suwanich, P. Geological map of Phutok formation improvement explored from Potash and Rock Salt drilled holes, topography and outcrops on the Khorat Plateau. Asia-Pac. J. Sci. Technol. 2012, 17, 58–70. [Google Scholar]
  29. Meesook, A. Cretaceous. In The Geology of Thailand; Ridd, M.F., Barber, A.J., Crow, M.J., Eds.; Geological Society of London: London, UK, 2011; pp. 169–184. [Google Scholar]
  30. Mouret, C. Geological history of northeastern Thailand since the Carboniferous. In Proceedings of the Relations with Indochina and Carboniferous to Early Cenozoic evolution model. In Proceedings of the International Symposium on Stratigraphic Correlation of Southeast Asia, Bangkok, Thailand, 15–20 November 1994; Department of Mineral Resources: Bangkok, Thailand, 1994; pp. 132–158. [Google Scholar]
  31. Utha-aroon, C.; Coshell, L.; Warren, J.K. Early and late dissolution in the Maha Sarakham Formation: Implications for basin stratigraphy. In Proceedings of the International Conference on Geology, Geotechnology and Mineral Resources of Indochina (GEOINDO′95), Khon Kaen, Thailand, 22–25 November 1995; pp. 275–286. [Google Scholar]
  32. Raksaskulwong, L. Mapping of the Stratigraphic Mess of the Phu Thok Formation in the Vicinity of Phu Thok Noi Area, Siwilai District, Nong Kai Province; Geology Division, Department of Mineral Resources: Bangkok, Thailand, 2002; 57p. [Google Scholar]
  33. Imsamut, S. Stratigraphy and Petrology of Phu Thok Formation at Khao Phu Thok and Phu Wue, Nong Khai Province; Geology Division, Department of Mineral Resources: Bangkok, Thailand, 1997. [Google Scholar]
  34. Hasegawa, H.; Imsamut, S.; Charusiri, P.; Tada, R.; Horiuchi, Y.; Hisada, K.I. ‘Thailand was a desert’ during the mid-Cretaceous: Equatorward shift of the subtropical high-pressure belt indicated by eolian deposits (Phu Thok Formation) in the Khorat Basin, northeastern Thailand. Isl. Arc 2010, 19, 605–621. [Google Scholar] [CrossRef]
  35. Assavapatchara, S. Lithostratigraphy of the Nawa Member (Preliminary Investigation). In Proceedings of the 12th Regional Congress on Geology, Mineral and Energy Resources of Southeast Asia GEOSEA 2012, Geoscience in Response to the Changing Earth, Bangkok, Thailand, 7–8 March 2012; p. 57. [Google Scholar]
  36. Sattayarak, N. Review of the continental Mesozoic stratigraphy of Thailand. In Workshop on Stratigraphic Correlation of Thailand and Malaysia; Geological Society of Thailand: Bangkok, Thailand; Geological Society of Malaysia: Kuala Lumpur, Malaysia, 1983; pp. 127–148. [Google Scholar]
  37. DMR. Geologic Map of Nong Kai Province; Department of Mineral Resources (DMR), Ministry of Natural Resources and Environment: Bangkok, Thailand, 2009. [Google Scholar]
  38. DMR. The Invaluable Natural Heritage: Thai Geological Site; Department of Mineral Resources (DMR), Ministry of Natural Resources and Environment: Bangkok, Thailand, 2015. [Google Scholar]
Figure 1. Interpretation of the Thakhek fault using DEM, showing major faults extending from Laos into Thailand. The TB1 and TB2 traces are inferred to branch from the main strike-slip fault, extending into Northeast Thailand with a general NW–SE orientation. The locations of Figures 3 and 5 are indicated on TB1 trace [9].
Figure 1. Interpretation of the Thakhek fault using DEM, showing major faults extending from Laos into Thailand. The TB1 and TB2 traces are inferred to branch from the main strike-slip fault, extending into Northeast Thailand with a general NW–SE orientation. The locations of Figures 3 and 5 are indicated on TB1 trace [9].
Geosciences 14 00196 g001
Figure 2. Geologic map of study area in Northeast Thailand and western Laos with the indicated locations of Figure 5a and a cross-section presenting interpreted faults (modified from [2]).
Figure 2. Geologic map of study area in Northeast Thailand and western Laos with the indicated locations of Figure 5a and a cross-section presenting interpreted faults (modified from [2]).
Geosciences 14 00196 g002
Figure 3. Two-dimensional ERT tomogram of profile BK01 (a) displaying resistivity tomography based on the L1–norm (robust) inversion method. The images show mudstone (low resistivity) flanked by coarse-grained sediment layers (high resistivity). Profile BK02 (b) displays a sharp contact between mudstone and sediment layers.
Figure 3. Two-dimensional ERT tomogram of profile BK01 (a) displaying resistivity tomography based on the L1–norm (robust) inversion method. The images show mudstone (low resistivity) flanked by coarse-grained sediment layers (high resistivity). Profile BK02 (b) displays a sharp contact between mudstone and sediment layers.
Geosciences 14 00196 g003
Figure 4. Two-dimensional ERT tomogram of profile BK03 displaying resistivity tomography based on the L1–norm (robust) inversion method. The image shows coarse-grained sediment layers (high resistivity) flanked by mudstone (low resistivity).
Figure 4. Two-dimensional ERT tomogram of profile BK03 displaying resistivity tomography based on the L1–norm (robust) inversion method. The image shows coarse-grained sediment layers (high resistivity) flanked by mudstone (low resistivity).
Geosciences 14 00196 g004
Figure 5. Photographs displaying the deployment of the 2D ERT profiles BK02 (a) and BK03 (b). An overview of the trench locations (c) and exposed walls of the NWY-1, with an inset showing the size of gravel compared to a water bottle (d), and NWY-2 (e) trenches. The inset photo in (e) shows a close-up view of the contact between mudstone and sediments, depicting the grain size distribution of unconsolidated sediments, while (f) shows the exposed walls of the NWY-3.
Figure 5. Photographs displaying the deployment of the 2D ERT profiles BK02 (a) and BK03 (b). An overview of the trench locations (c) and exposed walls of the NWY-1, with an inset showing the size of gravel compared to a water bottle (d), and NWY-2 (e) trenches. The inset photo in (e) shows a close-up view of the contact between mudstone and sediments, depicting the grain size distribution of unconsolidated sediments, while (f) shows the exposed walls of the NWY-3.
Geosciences 14 00196 g005
Figure 6. Photographs of the trench walls in NWY-1 (a) and NWY-2 (c), along with the stratigraphy of NWY-1 (b) and NWY-2 (d) trenches. Thrust faulting results in a main strike slip movement characterized by as a positive flower structure and antiform (e).
Figure 6. Photographs of the trench walls in NWY-1 (a) and NWY-2 (c), along with the stratigraphy of NWY-1 (b) and NWY-2 (d) trenches. Thrust faulting results in a main strike slip movement characterized by as a positive flower structure and antiform (e).
Geosciences 14 00196 g006
Figure 7. Photograph of the trench wall in NWY-3 (a) along with the stratigraphy (b).
Figure 7. Photograph of the trench wall in NWY-3 (a) along with the stratigraphy (b).
Geosciences 14 00196 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arjwech, R.; Hongsresawat, S.; Chaisuriya, S.; Rattanawannee, J.; Kanjanapayont, P.; Youngme, W. Identification and Verification of the Movement of the Hidden Active Fault Using Electrical Resistivity Tomography and Excavation. Geosciences 2024, 14, 196. https://doi.org/10.3390/geosciences14080196

AMA Style

Arjwech R, Hongsresawat S, Chaisuriya S, Rattanawannee J, Kanjanapayont P, Youngme W. Identification and Verification of the Movement of the Hidden Active Fault Using Electrical Resistivity Tomography and Excavation. Geosciences. 2024; 14(8):196. https://doi.org/10.3390/geosciences14080196

Chicago/Turabian Style

Arjwech, Rungroj, Sutatcha Hongsresawat, Suriyachai Chaisuriya, Jetsadarat Rattanawannee, Pitsanupong Kanjanapayont, and Winit Youngme. 2024. "Identification and Verification of the Movement of the Hidden Active Fault Using Electrical Resistivity Tomography and Excavation" Geosciences 14, no. 8: 196. https://doi.org/10.3390/geosciences14080196

APA Style

Arjwech, R., Hongsresawat, S., Chaisuriya, S., Rattanawannee, J., Kanjanapayont, P., & Youngme, W. (2024). Identification and Verification of the Movement of the Hidden Active Fault Using Electrical Resistivity Tomography and Excavation. Geosciences, 14(8), 196. https://doi.org/10.3390/geosciences14080196

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop