E-ISSN: 2528-388X
P-ISSN: 0213-762X
INERSIA
Vol.
No.
May 2023 Evaluation of Flood Control Performance Watershed.
Samarinda.
East Kalimantan Talangsari Joko Sugeng Hariyadi*.
Istiarto, and Adam Pamudji Rahardjo Department of Civil and Environmental Engineering.
Universitas Gadjah Mada.
Yogyakarta 55281.
Indonesia
ABSTRACT
Keywords:
River flood Flood control Normalization HEC-RAS Talangsari River The flood issues in the Talangsari River in Samarinda City were not only caused by high rainfall but also by the significant sediment accumulation and waste along the river, the encroachment of residential areas onto the river section, and the lack of public awareness regarding environmental cleanliness.
As a result, flooding occurs throughout the Talangsari River, from upstream to downstream, indicating significant damage to the river.
Considering these issues, proactive measures are necessary to address the flood problems.
One such measure is the normalization of the Talangsari River to ensure smooth land transportation from and to Apt Pranoto Airport in the Sei Siring area.
Bontang City.
Sangatta, and Kutai Kartanegara.
The flood discharge plan was estimated using Nakayasu and SCS synthetic unit hydrograph method.
The analysis revealed a peak discharge of 6.
39 m3/s for the existing conditions .
-year return perio.
38 m3/s for the normalization conditions .
-year return perio.
Subsequently, the simulations were performed using the HEC-RAS program by two scenarios: the existing river condition and the normalization scenarios.
Following the simulation with the 2-year return period discharge for the existing conditions, the Talangsari River experienced upstream to downstream flooding, underscoring the need for flood mitigation measures.
The river normalization simulation using the 25-year return period discharge considered dimensions such as a channel width of 7m and a depth of 3.
As a result of the normalization simulation, the Talangsari River successfully accommodated a flow of 13.
38 m3/s without experiencing flooding, confirming that the decision to pursue normalization was the best choice.
This is an open-access article under the CCAeBY license.
Introduction A river is a natural or artificial watercourse and container consisting of a flowing water network, starting from the source, and ending at the mouth, bounded on the right and left by its banks .
Over time, rivers have become more than just a water source.
they have become essential to communities in meeting their daily needs.
Like the Talangsari River in Samarinda City.
Samarinda Utara District.
East Kalimantan Province, it has transformed from a natural water reservoir into a commercial and residential area.
The consequence of this change is the frequent and regular flooding in the area.
Samarinda City, as a supporting city for the country's future capital, has experienced rapid population growth.
This population explosion is a contributing factor that triggers the flood issues in Samarinda.
The reduction of *Corresponding author:
E-mail: jokosugeng95@mail.
https://dx.
org/10.
21831/inersia.
Received October 28th, 2022.
Revised May 27th, 2023.
Accepted May 28th, 2023 Available online May 31th, 2023 green areas resulting from land use changes affects the hydrological characteristics and the river's potential as a natural reservoir.
In general, the current condition of the Talangsari River can no longer accommodate the flood discharge, leading to widespread flooding in the Talangsari River Basin from upstream to downstream.
This overflow is caused by high rainfall intensity, sedimentation accumulation, and the relatively small cross-sectional area of the Talangsari River, ranging from approximately 1 m to 3.
5 m.
Based on the topography of Samarinda City, about 30% of the areas or regions are situated below the river level, contributing to the high occurrence of flooding in Samarinda (Figure .
The topographic map can be seen in Figure 3.
Considering the factors contributing to the flooding in the Talangsari River, river normalization is the best countermeasure to mitigate the flood issue due to the limited conditions.
River normalization on the Talangsari Joko Sugeng Hariyadi, et al.
INERSIA.
Vol.
No.
May 2023 River can be done by dredging and widening river dimensions to increase capacity.
Along the Talangsari River are several main drainage points, including those from Jl.
Talangsari.
Jl.
Mugirejo, the Citraland housing area, and the Alaya housing area.
The Talangsari River Basin is a drainage system outlet for residential areas, commercial areas, or economic centres along Jl.
DI Panjaitan, starting from Gunung Tangga/Talangsari to Simpang 3 Jl.
DI Panjaitan in Samarinda.
Figure 2 illustrates the Talangsari River Basin Map.
In this study, the river length is 1.
7 km, with a watershed area of 2.
81 km2.
The Talangsari River mouth joins the Karangmumus River approximately 50 m from the PM Noor Bridge.
The Talangsari River Basin exhibits a subdendritic flow pattern.
This sub-dendritic flow pattern occurs when multiple tributaries or branches contribute, resembling branches of a tree merging into the main river.
Figure 1 The location map of Talangsari River in Samarinda City.
North Samarinda District.
East Kalimantan.
(Source: Google Eart.
Figure 2 Map of the Talangsari Watershed in Samarinda City.
North Samarinda District.
East Kalimantan Province.
(Map compiled based on DEMNAS dat.
INERSIA.
Vol.
No.
May 2023 Joko Sugeng Hariyadi, et al.
Figure 3 Topographic map of the Talangsari Watershed in Samarinda City.
East Kalimantan Province.
goodness-of-fit tests determine whether the selected distribution type is appropriate for the given data.
Method The research methodology includes hydrological analysis and hydraulic modelling using the HEC-RAS software.
Data preparation is required in the hydrological analysis step, including rainfall data to calculate annual rainfall and design rainfall.
The design rainfall is used to calculate the design flood discharge.
Forty-three events of extreme rainfall data .
reater than 50 mm/da.
were collected to obtain the design rainfall by frequency analysis.
Table 1 Statistic parameters requirement for probability Distributi Requirement Cs = 1.
Gumbel Ck = 5.
Cs OO 0 Normal Ck OO 3 Cs = C vA 3 Cv Log Ck = C v 8 6 C v 6 15 C v 4 Normal 16 Cv2 3 Log Pearson Aside from the above values i Among these events, the dominant rainfall duration is 11 hours, occurring in 7 rainfall events.
The data for extreme rainfall events with depths exceeding 50 mm can be obtained using WRPlot or Excel to determine the frequency of each depth and select the most frequently occurring duration as the dominant duration.
The distribution pattern calculation will utilise the data for rainfall events exceeding 50 mm.
Frequency analysis is conducted to determine the magnitude of design rainfall based on specific design Several distributions are commonly used in frequency analysis, including Normal.
Log Normal.
Gumbel, and Log Pearson i.
Each distribution requires different data parameters, namely the coefficient of skewness (C.
and coefficient of kurtosis (C.
, as specified in Table 1.
Frequency analysis calculations can be performed using the following calculations:
xT = K T E where xT is estimated value in a certain period.
A is standard deviation.
KT is frequency factor, and A is the average value of occurance.
The frequency factor equation was developed by Ven Te Chow in 1951 .
and is applicable to many probability distributions used in hydrological frequency analysis.
For specific distributions, the K-T relationship can be established between the frequency factor and the corresponding return period.
This analysis is performed using an Excel frequency analysis program .
After performing the frequency analysis calculations, the next step involves conducting goodness-of-fit tests using the Chi-square and Smirnov-Kolmogorov tests.
These Joko Sugeng Hariyadi, et al.
INERSIA.
Vol.
No.
May 2023 After calculating the design rainfall using frequency analysis and conducting goodness-of-fit tests, the next step involves computation using the synthetic unit hydrograph method.
The synthetic unit hydrograph methods used are Nakayasu and SCS.
Synthetic unit hydrographs can be utilized when historical data is data for 20 years .
obtained from the Temindung Rain Gauge Station (BMKG) .
The Temindung Rain Gauge Station is the closest station to the study area.
The summary of annual maximum daily rainfall data can be seen in Table 2.
Table 2 Daily maximum rainfall data of Temindung Station The hydraulic analysis is performed using the HEC-RAS The HEC-RAS program is used to calculate and simulate steady flow and unsteady flow conditions.
The HEC-RAS simulation is conducted using a 1D flood modelling approach.
The program also calculates water surface profiles along the river reach.
Input data for the program include river cross-sections, longitudinal profiles, hydraulic parameters (Manning's roughness and slop.
, river structures, flow discharges, and water levels.
The program outputs can consist of tables and graphs, including river schematics, cross-section plots, profiles, rating curves, and stage and flow hydrographs .
Year Rainfall .
Year Rainfall .
To determine the design rainfall, it is necessary to perform frequency analysis, which aims to establish the relationship between the magnitude of extreme events and their frequency using probability distributions, based on the annual maximum daily rainfall data as shown in Table Commonly used probability distributions include Normal.
Log Normal.
Gumbel, and Log Pearson i.
There are two methods to test whether the selected distribution is appropriate for the available data, namely the ChiSquare test and the Smirnov-Kolmogorov test .
There are two approaches to mitigate the flood issue:
structural and non-structural methods.
Non-structural methods involve sediment removal to restore the river's optimum function, while structural methods involve the construction of embankments to prevent overflow.
Structural methods are still the primary choice for flood management efforts .
The construction of embankments can only be carried out after non-structural methods have been implemented, but flooding persists .
As each phase of flood control works is completed, the capacity to handle flood discharges increases.
Therefore, when the final phase of work is completed, the flood control system can function as planned.
These efforts also significantly impact human society's economic, social, institutional, and environmental aspects .
Table 3 Chi-squared test results Probability Distribution Normal Log Normal Gumbel Log Pearson i N2 Calculated N2 Cr Information Accepted Accepted Accepted Table 4 Smirnov-Kolmogorov test results Results 1 Hydrological Analysis In this study, hydrological analysis is used to calculate the flow occurring in the Talangsari Watershed, specifically the design flood discharge.
Probability Distribution Normal Log Normal Gumbel Log Pearson i i max i criticism Information Accepted Accepted Accepted Accepted From Table 3 and Table 4, it is evident that the distributions yield different results.
The Chi-Square test was conducted using the Normal distribution, while the Smirnov-Kolmogorov test was performed using the Log Pearson i distribution.
Both distributions also differ in determining the magnitude of design rainfall.
The design rainfall for a specific return period from both distributions can be observed in Table 5.
In the rainfall-runoff analysis to estimate the design flood discharge, input of the design rainfall is required into a watershed system.
The rainfall design can be in the form of rainfall depth at a specific point or a rainfall hyetograph, representing rainfall distribution as a function of time during extreme rainfall events .
In this study, the calculation of design rainfall uses maximum daily rainfall INERSIA.
Vol.
No.
May 2023 Joko Sugeng Hariyadi, et al.
Table 5 Design rainfall with specific return period Return P Pearson i Log P Gumbel Period Distribution Distribution .
accurate with high resolution .
The analysis of rainfall distribution patterns was conducted using data on dominant rainfall events with depths exceeding 50 mm, which are considered representative of the studied area.
There were 42 events with rainfall depths exceeding 50 mm, and among them, 7 events had a dominant duration of 11 hours.
To determine the dominant duration, software like WRPlot or Excel can be utilized to identify the most frequently occurring depth values for each hour, and the duration with the highest frequency is selected as the dominant duration.
After conducting the rainfall distribution analysis, the rainfall intensity for a 2-year return period was determined to be 94 mm, while for a 25year return period, it was found to be 180 mm.
The analysis data for rainfall exceeding 50 mm is summarized in Table 6, and the rainfall distribution pattern is presented in Table 7, with graphical representations shown in Figure 4 and Figure 5.
The selected design flood return period should function effectively in terms of timing, structural considerations, and functionality .
In this study, the researcher opted for a 25-year return period, which would be used during the normalization simulation using the HEC-RAS As shown in Table 5, the design rainfall for a 25year return period using the Log Pearson i distribution is 180 mm, while for the Gumbel distribution, it is 178 mm.
Therefore, in this study, the design rainfall value to be used for further analysis is 180 mm, based on the Log Pearson i distribution.
The Log Pearson i distribution was chosen because it yielded a higher value compared to the Gumbel distribution.
Oc Extreme rainfall events at a rain gauge station can be calculated through observation using rainfall data from automatic stations.
In this study, the rainfall data used consisted of extreme rainfall events over a 5-year period .
obtained from JAXA, with the coordinates of the Temindung Rain Gauge Station, which is the nearest station to the research location .
The use of JAXA's extreme rainfall data can serve as a reference for extreme rainfall data as their global rainfall maps are highly t (%) Table 6 Rainfall distribution patterns Intensity P25 Intensity P2 180 mm 94 mm Table 7 Summary of rainfall distribution pattern > 50 mm P (%) % Average P Joko Sugeng Hariyadi, et al.
Rainfall .
INERSIA.
Vol.
No.
May 2023 and river slope in the watershed.
Based on the obtained data, the Talangsari watershed has an area of 2.
8 kmA, a river length of 1.
7 km, and a river slope of 0.
With these parameters, the analysis results for the unit hydrograph SCS method are presented in Table 8 and Figure 6, while the results for the Nakayasu method are shown in Table 9 and Figure 7.
Table 8 Nakayasu unit hydrograph of Talangsari Watershed Data Miles Analysis T0.
Tp H0.
Tp H 0.
5T0.
5 6 7
Time .
Figure 3 The 2-Year rainfall distribution pattern Rainfall .
Table 9 SCS unit hydrograph of DAS Talangsari Data Watershed Area (A) Main River Length (L) Miles Average River Slope (S) Analysis Concentration Time (T.
Time Lag .
Duration of Effective Rain .
Peak Time (T.
Base Time (T.
Peak Discharge (Q.
5 6 7 8 9 10 11 Time .
Figure 4 The 25-Year rainfall distribution pattern The parameters used in the unit hydrograph SCS and Nakayasu methods are the watershed area, river length.
Unit hydrograph Corrected Unit Figure 5 Nakayasu's unit hydrograph Graph INERSIA.
Vol.
No.
May 2023 Joko Sugeng Hariyadi, et al.
Unit hydrograph Corrected Unit Figure 6 SCS unit hydrograph Effective rainfall is a portion of excess rainfall that results in direct runoff.
The magnitude of effective rainfall is influenced by land conditions, watershed characteristics, and the amount of rainfall occurring in the river basin.
calculate losses using the SCS-CN method, several calculation parameters are required, including the Curve Number (CN) derived from land use data and soil In this study, soil type C is used, which represents soil with moderately high runoff potential and slow infiltration rate when the soil is fully saturated.
In this study, the simulation is conducted with two scenarios: the existing condition simulation and the simulation with normalization.
From the loss analysis, the following values were obtained: CN II = 81.
57, initial abstraction (I.
= 4.
66, and retention parameter (S) = 23.
28 mm.
These loss parameter calculations are used to calculate the effective rainfall (Pef.
at the study location.
The calculation of effective rainfall is based on the results of previous analyses, including the distribution pattern of rainfall, design rainfall calculation, and loss calculation.
The results of the effective rainfall calculation can be seen in Table 10 and Table 11.
Table 10 Two-years effective rainfall OcP Oc Peff 94 mm Peff Table 11 Twenty five-years effective rainfall P25 OcP Oc Peff Peff The flood hydrograph calculation is used to determine the design flood discharge.
The calculation is performed by multiplying the unit hydrograph with the flood After conducting the flood hydrograph calculation, peak discharges for the 2-year and 25-year return periods are obtained, which will be used for simulation in the HEC-RAS program.
The summary of peak discharges can be seen in Table 12.
Table 12 Peak discharge Peak Discharge .
Unit Hydrograph Nakayasu
Unit Hydrograph SCS
Q25
Q25
2 Hydraulic Modeling with HEC-RAS The hydraulic modeling simulation aims to analyze the channel's capacity to accommodate the given flow Joko Sugeng Hariyadi, et al.
INERSIA.
Vol.
No.
May 2023 Figure 7 Talangsari river geometry input The existing flow simulation aims to determine the capacity of the existing channel cross-section to receive and convey the flood discharge, allowing the identification of overflowing sections.
The first step in the flow simulation using HEC-RAS is inputting the geometry data.
The cross-sectional and longitudinal profiles of the river are created and inputted based on field measurement data, while the Manning's roughness coefficient .
is adjusted based on direct field investigations .
A higher Manning's roughness coefficient indicates a rougher surface.
The Manning's n values for each section of the channel are shown in Table 13, and the geometry of the Talangsari Watershed can be seen in Figure 8.
resembling real occurrences at the research site.
Unsteady flow simulation requires two boundary conditions:
upstream boundary conditions and downstream boundary For upstream boundary conditions, the simulation utilizes a hydrograph with a 2-year return period discharge.
It is assumed that the 2-year return period discharge can represent the existing flood events.
As for downstream boundary conditions, tidal data is used for the modeling, specifically the tidal stage hydrograph at the downstream of the Talangsari River.
The tidal data is incorporated into HEC-RAS as a stage hydrograph.
The tidal stage utilized corresponds to the tidal range between the confluence of the downstream of the Talangsari River and the Karang Mumus River, with a magnitude of 5.
Table 13 The Manning's roughness coefficient .
of the Talangsari River
Value n manning Point
LOB
Channel
ROB
Sta.
0 000 - sta.
Sta.
0 060 - sta.
Sta.
0 082 - sta.
Sta.
0 257 - sta.
Sta.
0 282 - sta.
Sta.
0 332 - sta.
Sta.
0 357 - sta.
Based on the available data, the simulation of the existing conditions is conducted using the HEC-RAS software.
the simulation of the existing conditions, the river is subjected to a discharge of 1.
28 m3/s.
Additionally, the Froude number at the maximum simulated discharge indicates a value below 1 .
r < .
, indicating that the flow type in the Talangsari River is subcritical.
The longitudinal profile .
ong sectio.
of the river in Figure 12 demonstrates the fluctuating flow, varying channel shapes that result in flow constriction, and potentially influence the upstream flow.
The cross-sectional profiles of the Talangsari River can be observed in Figures 9 to Figure In this study, the flow modeling using HEC-RAS is simulated with unsteady flow conditions.
Unsteady flow simulation is chosen to understand the changes in flow over time, aiming to depict flood events closely INERSIA.
Vol.
No.
May 2023 Joko Sugeng Hariyadi, et al.
Figure 8 Cross-sectional profile of the Talangsari River at Sta.
Figure 9 Cross-sectional profile of the Talangsari River at Sta.
Figure 10 Cross-sectional profile of the Talangsari River at Sta.
Figure 11 Longitudinal profile of the Talangsari River during the existing simulation.
Based on the simulation of the existing conditions, it is evident that the entire Talangsari River experiences Therefore, it is necessary to conduct river normalization along the entire stretch of the Talangsari River due to significant morphological changes, such as channel narrowing and bed sedimentation, as depicted in Figure 12.
The normalization simulation will cover the reach from station 0 000 to station 1 700, utilizing a 25year return period discharge of 13.
38 m3/s.
The proposed channel shape for the normalization of the Talangsari River is rectangular, with a bottom width of 7 m, channel depth of 2.
5 m, and a guard height of 1 m.
The Joko Sugeng Hariyadi, et al.
INERSIA.
Vol.
No.
May 2023 normalization will utilize concrete as the construction material, and the Manning's roughness coefficient .
used for the simulation is 0.
The crosectional profiles and longitudinal profile of the Talangsari River for the normalization can be observed in Figure 13 to Figure 16.
This indicates that the performed normalization was successful in addressing the flood issue in the river.
Normalization of the river is crucial in increasing the channel capacity and reducing flood risk in the area.
optimizing the dimensions of the channel, river structures, and minimizing flow obstructions, normalization can help control water flow and mitigate the risk of flooding.
After the normalization process with a 25-year return period discharge.
Sungai Talangsari no longer experiences Figure 12 Cross-sectional profile of the Talangsari River at Sta.
0 000 after normalization Figure 13 Cross-sectional profile of the Talangsari River at Sta.
0 762 after normalization Figure 14 Cross-sectional profile of the Talangsari River at Sta.
1 700 after normalization Figure 15 Longitudinal profile of the Talangsari River after normalization.
INERSIA.
Vol.
No.
May 2023 Joko Sugeng Hariyadi, et al.
Istiarto.
Modul Pelatihan Pemakaian HECRAS.
Yogyakarta, 2014.
Benny Sukmara and R.
Shyan Wu.
AuThe 20 th Cross-strait Symposium on Environment.
Resources and Ecological Conservation s uOAioEaIouAnAi Flood Management in Samarinda.
Indonesia:
Recent Progress,Ay 2021, doi:
20944/preprints202102.
Kodoatie.
Rekayasa dan manajemen banjir kota.
Junaidi.
Nurhamidah, and D.
Daoed.
AuFuture flood management strategies in Indonesia,Ay in MATEC Web of Conferences.
EDP Sciences.
Nov.
1051/matecconf/201822901014.
Triatmodjo.
Hidrologi Terapan.
Yogyakarta: Beta Offset, 2019.
BMKG.
AuCurah Hujan,Ay BMKG, 2022.
https://dataonline.
id/home .
ccessed Jan.
15, 2.
Harto.
Hidrologi Teori.
Masalah, dan Penyelesaian.
Yogyakarta: Nafiri Offset, 2009.
JAXA Global Rainfall Watch.
AuJAXA Global Rainfall Watch (GSMaP) ,Ay JAXA, 2022.
https://sharaku.
jp/GSMaP/index.
ccessed Jan.
15, 2.
TRMM / GPM / GSMaP.
AuHow about using rain data?Ay [Onlin.
Available:
http://sharaku.
jp/GSMaP/index_e.
Learn Culvert Studio.
AuManningAos n-values,Ay Learn Culvert Studio, 2022.
ccessed Apr.
01, 2.
Syahrin.
Rukmana.
Rachman.
Harakan, and F.
Tadjine.
AuSister city bottom up collaborative model plan as a major flooding mitigation in Samarinda,Ay European Alliance for Innovation n.
Jun.
4108/eai.
8-12-2018.
Conclusion Based on the research findings, the peak discharge for the existing condition (Q.
was determined to be 6.
39 m3/s, while the peak discharge for the normalization scenario (Q.
was found to be 13.
38 m3/s.
However, during the simulation using these discharges, it was observed that Sungai Talangsari could not convey and accommodate the increased flow, highlighting the need for further simulation to ensure Sungai Talangsari's capacity to handle larger discharges.
A comprehensive normalization approach was implemented along the entire stretch of the river, employing a rectangular channel shape with a bottom width of 7 m, a channel depth of 2.
5 m, and a guard height of 1 m.
Concrete materials were selected to normalize Sungai Talangsari, with a Manning's roughness coefficient .
After conducting the normalization modeling.
Sungai Talangsari effectively handled a discharge of 13.
38 m3/s (Q.
Additionally, proactive measures can be taken by the community to prevent flooding, such as providing training to enhance preparedness and resilience during emergency Active participation and engagement of the community are crucial aspects in effectively mitigating the impacts of floods in civil engineering projects .
References