KURVATEK Vol. No. April 2025, pp. e-ISSN: 2477-7870 p-ISSN: 2528-2670 ANALYSIS OF NATURAL COMPACTION IN MRT TUNNEL CP201. CENTRAL JAKARTA Muhammad Adi Naufal1*. Imam Achmad Sadisun1. Rendy Dwi Kartiko1 Department of Geological Engineering. Institut Teknologi Bandung Jalan Ganesa No. Lebak Siliwangi. Coblong. Bandung. Jawa Barat 40132. Indonesia *Email corresponding: adinaufalmuhammad. How to cite: M. Naufal. Sadisun, and R. Kartiko, "Analysis of Natural Compaction in MRT Tunnel CP201. Central Jakarta". Kurvatek, vol. 10, no. 1, pp. April 2025. 33579/krvtk. 5539 [Onlin. Abstract Ai This study analyzes the land subsidence caused by natural consolidation in the Thamrin Area. Central Jakarta, along the MRT CP201 tunnel route. JakartaAos soft soils, particularly clay layers, contribute to land subsidence, with 0. 2 to 3. 6 meters and rate of subsidence is 1,24 cm/years. Rate of subsidence refers to the speed at which the ground surface lowers over a specific period due to many factors . Calculating rate of subsidence is essential for assessing its impact on infrastructure and predicting subsidence predictions. Using both analytical . ne-dimensional consolidatio. and numerical . inite elemen. methods, subsidence predictions were made based on geotechnical data from three boreholes which is TS-01,TS-02 and TS-03. Results show positive correlation between the two methods, with numerical approach yielding slightly smaller values due to its consideration of both vertical and horizontal deformations. This analysis highlights unavoidable nature of subsidence due to natural consolidation and emphasizes the need for mitigation and adaptation measures, such as spatial planning and monitoring systems, to minimize impacts on infrastructure. Keywords: Land subsidence. Natural consolidation. One-Dimensional Terzaghi. Finite element INTRODUCTION MRT Jakarta has been designed as one of the integrated modes of transportation to overcome the problem of congestion. The current Jakarta MRT development project work program in phase 2 stretches along about 11. 8 kilometers from Bundaran HI Station to West Ancol Station. Phase 2 development is a national strategic project based on Presidential Regulation Number 56 of 2018 about Second Amendment to Presidential Regulation Number 3 of 2016 on Accelerating the Implementation of National Strategic Projects . Jakarta MRT Tunnel CP201 construction is an important part of the development of the public transportation system in Central Jakarta with total length of 2. 6 km and depth of 17 to 36 meters underground . Jakarta, which is located in land subsidence prone area, has geological conditions dominated by soft soils, such as clay and silt, which can affect the stability of the tunnel. Research shows that the conditions of the soil below the surface are various, with soft soil layers having high plasticity, which has the potential to cause deformation and surface subsidence during and after the tunnel excavation process . Land subsidence is a geological phenomenon that occurs when the soil surface decreases due to various factors, including human activities and geological conditions. Land subsidence are problem of several cities in Indonesia, especially in the coastal area of the North Coast of Java (Pantur. , such as Jakarta. Semarang. Pekalongan, and Semarang showed in Figure 1. In addition. Bandung is also an area with the potential for land subsidence. Land subsidence in Jakarta has been going on since 1970 . This subsidence can be caused by various factors, both natural and anthropogenic. In Jakarta, natural factors include natural compaction and tectonic activity, while anthropogenic factors primarily consist of groundwater extraction and structural loads from buildings. Land subsidence in Jakarta occurs spatially at a rate of 3Ae10 cm/year . Recent studies analyzing land subsidence in Central Jakarta using InSAR technology have indicated a slower rate, ranging from 0 to 2. 5 cm/year . One of the causes of land subsidence in Central Jakarta is natural compaction. Natural sediment compaction is a process where the volume of sediment layers decreases due to the weight of overlying sediments . Natural compaction referred to the study is natural consolidation. The assessment of land subsidence due to natural consolidation in Central Jakarta is essential for evaluating necessary mitigation and adaptation measures. This study aims to analyze the rate of land subsidence caused by natural compaction factors using analytical and numerical methods. Geotechnical data are collected Received November 30, 2024. Revised March 21, 2025. Accepted March 27, 2025 DOI : https://doi. org/10. 33579/krvtk. through secondary sources and laboratory tests to determine the parameters required for the subsidence Figure 1. Land subsidence of several city in Indonesia . Geology of Study Area Physiographically. West Java is divided into four zones: Batavia Plain Zone (Jakart. Bogor Zone. Bandung Zone, and Southern Mountain Zone. The study area is located in Batavia Plain Zone, which extends along the northern coast from Serang to Cirebon. This zone spans approximately 40 km in width and features relatively flat morphology. It is composed primarily of alluvial deposits transported by rivers. The Jakarta Basin consists of a thick sequence of Quaternary deposits measuring approximately 200250 meters, which overlay Tertiary basement rocks. The basal contact of these Quaternary layers serves as the lower boundary of the groundwater aquifer system. This sedimentary sequence can be divided into three main stratigraphic units based on their relative ages, arranged from oldest to youngest: . Pleistocene marine and non-marine sediments, . volcanic fan deposits from the late Pleistocene, and . Holocene deposits, including marine and floodplain sediments . Based on the regional geological map. Turkandi et al. classified the rocks in the study area, ranging from the Early Miocene to the Holocene, into several formations. Regional stratigraphy of the study area, from youngest to oldest, as follows: aluvium (Q. Coastal Ridge Deposits (Qb. Alluvial Fan Deposits (Qa. Coral Limestone (Q. Andesite of Sundamanik Mount (Qva. Young Volcanic Rock (Q. Banten Tuff (QTv. Serpong Formation (Tpa. Genteng Formation (Tp. Bojongmanik (Tm. Klapanunggal Formation (Tm. Jatiluhur Formation (Tm. and Rengganis Formation (Tmr. The research location is in Central Jakarta, with geographic coordinates ranging from 106A48Ao44Au106A49'58. 61Au East longitude and 6A9'35. 72Au-6A11'27. 60Au South latitude. Figure 2 shows that the geology at the research area alluvium fan deposits . This deposit consists of fine-grained layered tuff, sandy tuff, interspersed with conglomeratic tuff. The formation process of this deposit originates from the young volcanic rocks in the highlands of Bogor, which were deposited in a terrestrial environment, forming a fan These deposits are estimated to be of Pleistocene age. The soil conditions in Central Jakarta, particularly in the Menteng subdistrict, are dominated by soft to very soft clay layers. KURVATEK Vol. No. April 2025: 59 Ae 66 KURVATEK e-ISSN: 2477-7870 p-ISSN: 2528-2670 Figure 2. Modification geological map of study area . II. METHODS This research used analytical methods and numerical approach based on soil properties. We use boreholes on the track of the CP201 MRT tunnel to describe soil and soil parameter for analysis subsidence. Total of three boreholes core samples were described each to determine the value of soil properties. The model's parameters were derived from laboratory measurements conducted on undisturbed soil samples retrieved from boreholes TS-1. TS-2, and TS-3 in Thamrin-Monas Area. TS-1 and TS-2 are located in Menteng Subdistrict, whereas TS-3 is located is Gambir Subdistrict. Clay at the study area assumed homogeneous at the same layer. Table 1 shows details measurements for land subsidence analysis. Table 1: Engineering properties of Thamrin-Monas Area Section. Section ThamrinMonas Saturated unit weight (Kn/mA) Unsaturated unit weight (Kn/mA) Void Coefficient of . eEyeE) Coefficient of . eEye ) . A/yea. Swelling Index Cohesion . Internal friction . Silty Clay Sand Clay Litology Modulus . 5x10A to 6x10A to 4x10A 1x10A to Permea . 1x10OA 1x10OA 1x10OA Figure 3 shows research flow diagram. The research stage starts from the preliminary stage as literature study, then continue data collection stage, the analysis stage, and the conclusion. This analysis carried out was to calculate the amount of land subsidence that occurred at the research area. Idenfitication of consolidation parameters used for land subsidence calculation to produce a comparison of analytical and numerical methods. Analysis of Natural Compacton in MRT Tunnel CP201. Central Jakarta (Muhammad Adi Naufal. Imam Achmad Sadisun and Rendy Dwi Kartik. DOI : https://doi. org/10. 33579/krvtk. Figure 3. Research flow diagram. LAND SUBSIDENCE ANALYSIS Land subsidence analysis using one-dimensional consolidation and numerical method using finite Borehole used for analysis are TS-1. TS-2, and TS-3. The three boreholes were analyzed at each When stress increases in saturated soil, it generates pore water pressure. In highly permeable sand, water drains rapidly, resulting in both immediate settlement and consolidation occurring simultaneously. Soils with lower permeability, such as clay, water disipated gradually under the applied load, causing a slow transfer of excess pore water pressure . to effective stress (EA. , as represented in Eq. OIyua = OIyuaAo OI yc Where are: OIyuaAo Oe increase in the effective stress OI yc Oe increase in the pore water pressure Rate of consolidation is governed by an equation derived by Terzaghi, based on assumptions that soil compression and pore water flow occur solely in the vertical direction, the soil is homogeneous, and both total and effective stress are uniform on any horizontal plane . Predictions of one-dimensional consolidation settlement are typically made using results from oedometer tests on representative clay soil KURVATEK Vol. No. April 2025: 59 Ae 66 KURVATEK e-ISSN: 2477-7870 p-ISSN: 2528-2670 Primary consolidation can be calculated using Equation 2, while secondary consolidation, which occurs after the pore water pressure has fully dissipated, can be estimated with Equation 3 . ya ya ycy OIycy Sp = 1 yca yce ycoycuyci ycoycuyci ( 0ycy Ss = yaAyu ya ycoycuyci ( ycy0 OIycy ycy0 Where are: = settlement due to primary consolidation = settlement due to secondary consolidation yaAyu = yayu /. yceycy ) ycy0 = effective stress OIycy = vertical stress yayca = compression index ya = thickness of layer yce0 = void ratio The finite element method for calculating soil subsidence is based on Biot's consolidation theory . which employs coupled consolidation equations. These equations enable the simultaneous calculation of excess pore pressure and deformation in a porous medium over time. In this model, several assumptions are made: the soil is homogeneous, uniform, and saturated. fluid flow adheres to Darcy's law. and a small strain poroelastic model is applied. s = s' m . steady pexces. Where are: = total stress vector = effective stress vector = unit vector for stress components psteady = pressure determined as the end of consolidation pexcess = excess pore water pressure i. RESULTS AND DISCUSSION Analytical Method Land subsidence is assumed to have started in 2018, with a subsidence value of 0. 058 meters. The analysis shows that the maximum subsidence will occur at point TS-02, reaching 3. 60 meters, with a maximum time of 242 years. Rate of subsidennce at this point is 1. 24 cm/year. The vertical and effective stress, representing the overburden are calculated for each clay layer, with an effective pressure of 486. kN and a vertical stress of 398. 9 kN. At point TS-02, the depth to the maximum groundwater level is 3. meters, and the clay layer thickness is 36 meters (Table . Higher the compression index value, the higher coefficient of consolidation value. This will affect the consolidation process which can occur faster. Result of the primary decrease at this point is 3. 51 m, while the secondary decrease value is 0. 09 m. Relationship between the degree of consolidation (U) and the time factor (T. is used to determine the maximum descent Table 2: Primary and secondary subsidence due to natural consolidation. Code TS-01 TS-02 TS-03 H . yyOAo yyOya yaeyau . 3,51 2,19 yaeyae . 0,07 0,09 0,07 yaeyayaya . 1,77 2,26 yaeyeayeCyeo . The slowest land subsidence in Thamrin-Monas Area is the TS-01 borehole with a sinking value of 77 m until 4001. (Table . The sinking rate on the TS-01 reaches 0. 09 cm/year. This indicates that the thicker the clay, the slower the time required for consolidation . This happens because clay is impermeable that can hold water, while the permeable layer at all drill depths is not described. Previous studies, consolidation process in clay is time-dependent, with the dissipation of excess pore water pressure occurring gradually, especially in soils with low permeability like clay. These studies highlight how water retention and slow dissipation contribute to prolonged settlement and consolidation time in clays . Analysis of Natural Compacton in MRT Tunnel CP201. Central Jakarta (Muhammad Adi Naufal. Imam Achmad Sadisun and Rendy Dwi Kartik. DOI : https://doi. org/10. 33579/krvtk. Figure 4. Land subsidence graph TS-02. Maximum subsidence time during the consolidation process can be very slow, making it impractical for construction purposes. The subsidence after 90% consolidation is insignificant, thus resulting in minimal impact. Table 3 shows the land subsidence until the 90% consolidation process. The slowest maximum subsidence occurs at TS-01, with a rate of 0. 17 cm/year and a maximum subsidence time of up to 944 years. On the other hand. TS-02 experiences the largest maximum subsidence of 3. 25 meters in 138 years, with a subsidence rate of 2. 82 cm/year. The consolidation process at 90% for point TS-02 will continue until the year 2156. Table 2: Primary and secondary subsidence due to natural consolidation. Code Subsidence . TS-01 TS-02 TS-03 3,24 2,04 Rate of Subsidence . m/tahu. 0,17 2,35 1,94 Time Subsidence . Years Subsidence Numerical Method Numerical method analysis using finite element with Plaxis 2D Software. Figure 5 shows the model and element mesh used for analysis by the finite element method at point TS-02, as the largest subsidence. The loads used for the analysis are based on the respective loads on the soil layer. Based on the model, land subsidence that occurred due to natural consolidation was 3. 25 m. The value of the land subsidence occurs until the consolidation process is completed. Figure 5. Land subsidence graph TS-02. KURVATEK Vol. No. April 2025: 59 Ae 66 KURVATEK e-ISSN: 2477-7870 p-ISSN: 2528-2670 Comparison of Analytical and Numerical Methods Results of land subsidence analysis using both analytical and numerical approaches show varying These discrepancies arise from the parameters used in the analysis as well as the stress considered in each method. The deformation assumed in the Terzaghi method approach occurs in only one direction, whereas the finite element method assumes deformation in two directions, vertical and horizontal. The percentage difference is calculated to assess the similarity of the results between the two methods. this point, there is a difference in the results between the analytical and numerical methods. The comparison of subsidence calculations using these two methods show positive correlation across all borehole points, with subsidence difference due to natural consolidation at TS-01 being 0. 45 meters. The comparison of subsidence calculations due to natural consolidation using the Terzaghi one dimentional analytical method also yields a positive correlation, with a difference of less than 15%. Figure 5. Land subsidence graph TS-02. IV. CONCLUSION Soil condition in MRT CP201 Tunnel especially in Thamrin Area is predominantly composed of soft to very soft clay layers. Loose to very loose sand lenses are trapped in the upper layers, while sand lenses are also found in the lower layers with the soil consistency becomes stiff to very stiff, dominated by clay and silt materials. Discrepancies in interpretation based on the regional geological map and borehole sample descriptions need further examination. These interpretations influence the formation present at the research The presence of clay layers beneath the surface contributes to natural consolidation. Overall, natural consolidation has impact on land subsidence ranging from 0. 2 to 3. 6 meters with rate of subsidence 35/year. Land subsidence in Central Jakarta due to natural consolidation is unavoidable. Adaptation measures are necessary to minimize the potential impacts, such as spatial planning, early warning systems, and monitoring rate of land subsidence. ACKNOWLEDGEMENTS The author would like to express my sincere gratitude to Pak Dr. Eng. Imam A. Sadisun. and Dr. Rendy Dwi Kartiko. , for their invaluable guidance and support throughout my research. Their expertise, patience, and insightful feedback have been instrumental in shaping this work. The author wishes to extend heartfelt gratitude to PT MRT Jakarta for their invaluable support and cooperation in providing the data crucial for this land subsidence study. The successful completion of this research would not have been possible without their generous contribution of data and permission to use it. Analysis of Natural Compacton in MRT Tunnel CP201. Central Jakarta (Muhammad Adi Naufal. Imam Achmad Sadisun and Rendy Dwi Kartik. DOI : https://doi. org/10. 33579/krvtk. REFERENCES