The impact of changing storage area on flood magnitude and occurrence.
Lakes and swamps form an important part of tropical water resources. They contribute significantly to the water balance acting as storage in the basin. Swamps occupy much larger tropical areas than lakes. Swamps are considered as standing water series in which the water motion is not that of a continuous flow in a definite direction, although a certain amount of water flow may occur as internal currents in the vicinity of the inlets and outlets . Lakes are regarded as inland bodies of standing water. Swamps are described as lowlands flooded in the rainy season and usually watery at all times. The presence of storages such as lakes and swamps in the basin greatly reduces the total runoff, introduces intermittency to the overall catchment runoff response, due to the effect of these stores attenuating or quite often terminating runoff. Smoothed values of the runoff depend on the percentage of storages in the basin and on the annual rainfall.
Storages, especially swamps, are often located near a river. The location of a storage is very crucial toward the determination of the dominant sources of stream-flow reaching the storages and also the catchment outlet.
In addition to the relative location of the storages within a catchment, contributing areas to the storages also was found to influence storage-overflow behaviour as well as the overall catchment response. For example, lakes have been found to be efficient in smoothing the hydrographs of incoming catchment runoff responses only when they comprised more than 5% of the catchment area [2, 3].
The present study was motivated by general observations drawn from previous studies [2, 3, 4], but specifically by the flooding problems experienced in Rantau Jaya Village, located in Way Pegadungan catchment, a sub-catchment of Way Seputih in Lampung Province, Sumatera. Natural swamps exist throughout the Way Pegadungan as well as Way Seputih catchments. Way Pegadungan catchment have experienced severe flood events, causing overflow into the nearby village, Rantau Jaya. The downstream part of Way Pegadungan river is swamp areas which include Rawa Betik and Tanjung Kramat. Rantau Jaya village is located in the downstream of Rawa Betik. This village experiences flooding every year, with flood duration can be two to three months which inundate some part of the village located on the side of the river. Due to this annual flooding problems, people who live in that area ask the government to build dyke.
In addition to that, this study has also been motivated by flooding problem in the area nearby Nagara River, South Kalimantan. Swamp areas along the river have acted as flood plain. The problem arises since there was a plan a few years ago to change the land use for plantation. To prevent the plantation from flooding coming from the Nagara River, there were scenarios to built dyke at one side of the river or at both sides. Nowadays the plan has been implemented.
This paper consists of two main parts of analysis. The first part is sensitivity analysis upon landscape controls on the catchment toward the occurrence and the peak of floods. The analysis involves utilizing rainfall model for data generation and developing a conceptual rainfall-runoff model. The first part of the paper present the results of a sensitivity analysis that had the objective of exploring and exposing the critical process controls of flood frequencies in a catchment-storage system, in this way providing guidance toward developing an appropriate monitoring scheme to predict frequency of occurrence of such catastrophic floods in the future, and preventing them from ever happening again, or mitigating their effects. The second part of the paper will look at the case studies, Way Pegadungan and Nagara River. Utilizing available data, hydrology and hydraulic data, some analysis will be done to present the results of hydraulic modeling with regards to the effect of changing the area of storages on flood peaks.
Flood Frequency and Flood Occurrences; Sensitivity Analysis
The flood frequency curve, typically estimated from observed flood records and widely used in flood estimation practice, is the culmination of complex interactions between climatic inputs (rainfall intensities, evaporation demand) and those landscape properties that have a bearing on the rainfall to runoff to flood peak transformation, presented within a stochastic framework [4, 5, 6, 7, 8]. While the approach adopted in part of this study is general enough, climate conditions and catchment (and storages) characteristics typical of Nagara River and Way Pegadungan catchments (where storages such as swamp areas exist along the river) will be used to parameterize the adopted conceptual models.
For the sensitivity analysis, a model which can simulate rainfall and responding runoff for very long period needs to be developed, in order to construct a long period flood frequency and to figure out the frequency of storage overflow. In order to do this, it needs a rainfall model which can generate as long period as it needs [9, 10]. Variables on the rainfall model are determined based on the analysis of the rainfall data from the study area. Rainfall intensity, duration and interstorm period are based on rainfall data analysis from the study area. Average annual rainfall in study area is 2300 mm/year. The rainfall runoff model will use a bucket model which has three catchment-lake/swamps thresholds, i.e. the catchment field capacity storage governing subsurface stormflow, total storage capacity governing catchment surface runoff, and lake/swamp storage capacity governing lake/swamp overflow.
The water balance equation for the catchment bucket is:
dS/dt = i(t) - [Q.sub.se](t) - [Q.sub.ss](t) - [E.sub.bs](t) - [E.sub.veg] (t) (1)
where S is catchment storage (mm), i is rainfall (mm/hr), [Q.sub.se] is surface runoff (mm/hr), [Q.sub.ss] is subsurface runoff (mm/hr), [E.sub.bs] is bare soil evaporation (mm/hr), [E.sub.veg] is transpiration (mm/hr) and t is time in hour. Transpiration ([E.sub.veg]) is a function of the fraction of vegetation covering the catchment, M, potential evaporation, [E.sub.p], and storage at field capacity ([S.sub.fc]). Bare soil evaporation ([E.sub.bs]) is a function of potential evaporation and the fraction of the catchment with no vegetation as shown in Equations 2 and 3.
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)
Subsurface flow ([Q.sub.ss]) occurs when the soil moisture storage in the bucket exceeds the field capacity threshold and is therefore a function of the dynamic storage, S(t), as shown in Equation 4. The nonlinear storage-discharge relationship for subsurface flow is represented by two coefficients a and b, which are derived from recession curve.
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)
Surface runoff ([Q.sub.se]) occurs when the soil is fully saturated, which occurs when the soil moisture storage S(t) exceeds the bucket capacity (Equation 5).
[Q.sub.Se] = S(t) - [S.sub.b] if S(t) > [S.sub.b] (5)
The water balance equation for the storage (lake/swamp) bucket is given by:
[dS.sub.S]/dt = i(t) + [A.sub.C]/[A.sub.S][Q.sub.in](t) - 0.7[E.sub.p](t) - [Q.sub.ses](t) (6)
where [S.sub.S] is lake/swamp storage, [A.sub.C] and [A.sub.S] are catchment and storage (lake/swamp) area respectively, [Q.sub.in] is inflow from the contributing catchment area,
[Q.sub.in] (t) = [Q.sub.ss] (t) + [Q.sub.se] (t) (7)
and [Q.sub.ses] is outflow from the storage (lake/swamp).
[Q.sub.ses] (t) = [S.sub.s] (t) - [S.sub.bs] if [S.sub.s] (t) > [S.sub.bs] (8)
where [S.sub.bs] is lake/swamp storage capacity. The values of parameters used in the simulation are presented in Table 1.
Results of Sensitivity Analysis
The ratio of contributing catchment area to storage area ([A.sub.C]/[A.sub.S]) significantly impacts the one-dimensional rate of inflow into the storage, and therefore the incidence of storage-overflow. Figure 1 shows the relative per-storm frequency of storage-overflow for a range of sensitivity experiments involving systematic alteration of landscape parameter values. The frequency of storage-overflow events is quantified by the ratio of the number of storage-overflow triggering storms to the total number of storms. Figure 1 presents the frequency of storage-overflows as generated using average rainfall intensity of 2300 mm/year which represents the average annual rainfall in the study area. Small values of [A.sub.C]/[A.sub.S] mean that storage overflow is contributed more by direct rainfall falling in the storage. As the value of [A.sub.C]/[A.sub.S] increases, catchment contribution supplements direct rainfall, and the frequency of storage overflow increases significantly until it reaches high value of [A.sub.C]/[A.sub.S] and then the increase of [A.sub.C]/[A.sub.S] will not increase frequency of storage overflow significantly.
[FIGURE 1 OMITTED]
Runoff generation and flooding are non-linear, threshold-driven processes. In a contributing catchment of a catchment-storage system, runoff which can be produced includes saturation excess runoff and subsurface storm flow. The existence of storages such as lakes or swamps in a catchment-storage system introduces additional runoff, storage overflow, in which storage overflow occurs when the storage capacity is exceeded. Figure 2 presents the results of an examination of the impact of landscape control upon both the frequency and magnitude of runoff response expressed in terms of the flood frequency curve. Note that the total area of the landscape, [A.sub.C] + [A.sub.S], is the same for all simulations reported here. It can be seen that the flood frequency simulated using [A.sub.C]/[A.sub.S] = 100 has higher values for every return periods compared with that simulated using [A.sub.C]/[A.sub.S] = 5. It means that the smaller the storage area or the larger the contributing catchment area, the larger the floods peaks that occur.
Hydraulic and Hydrology Analysis on Case Studies
To see the impact of Ac/As on the flood peaks, some case studies have been investigated. One of the case studies is Way Pegadungan catchment. The methods involved in the analysis include hydrology and hydraulic analysis. Hydrology analysis for Way Pegadungan used Synthetic Unit Hydrograph GAMA I . Some parameters used in the calculation of synthetic unit hydrograph are shown in Table 2.
[FIGURE 2 OMITTED]
Hydraulic analysis used HEC-RAS software , which assumes that the flow is non-uniform one dimensional and steady flow. Another widely used software for hydraulic analysis is Duflow , a computer program for one dimensional hydraulic modeling of surface water, which has also been used for flood modeling and mitigation [14, 15]. The result of simulation using HEc-RAS for condition without dyke (existing condition) and with dyke is presented in Figure 3.
The results presented in Figure 3 used 25 years return period of river discharge as it is shown in the legend as Q-25yr. Furthermore the legend shows EG Q-25yr, WS Q-25yr and Crit Q-25yr which represents Energy Grade, Water Surface and Critical Depth using design discharge 25 years return period respectively, and Bank Sta which represents the main channel bank station. The peak discharge is 1200 [m.sup.3]/sec. Water surface elevations at the station observed shows the difference between the results of without and with dyke on the sides of the river. Although it is obvious that the impact of dyke is an increase of water level in the channel, the extent of water increase remains an important question to be examined. From the hydraulic simulations using HEC-RAS  it is found that water surface elevations without dyke is 2.99 m (Figure 3a) and with dyke is 3.99 m (Figure 3b). Therefore the difference of water surface elevations with and without dyke is 1 m.
[FIGURE 3 OMITTED]
Relating this finding with our concern regarding the ratio of contributing catchment area to storage area, it shows that the ratio of contributing area to storage area has a significant effect in flood peaks. On the existing condition (without dyke) with wide flood plain due to swampy area on the side of the river, the ratio of contributing catchment area to storage area is about 30. Using dyke means reducing storage area, and for those simulations the ratio of contributing catchment area to storage area is about 250.
Nagara River Catchment
The length of Nagara River is 155 km and the catchment area is 7801.48 [km.sup.2]. Based on the computation using 2 day rainfall module (i.e. maximum total two day rainfall from available rainfall data), it is found that the rainfall module was 164 mm. Discharge is computed using Rational method and the value is 1750 [m.sup.3]/sec . The computed discharge is further used in hydraulic computation. Simulation using HEC-RAS  shows that the discharge value is represented by PF 1 as it is shown in the legends of Figure 4 (EG PF 1, WS PF 1 and Crit PF1).
When flooding happened, flood plain was inundated. As the area is flat, water moves to swamp area near the river which can reach few kilometers. In modelling using HEC-RAS  during flooding the swamp area became extension of flood plain. This can be seen in Figure 4 that river geometry from the top shows that flood plain is significantly wide in condition without dyke. Figure 4 also shows that along the river, it is divided into more than 100 cross sections. The arrow displays flow direction.
[FIGURE 4 OMITTED]
Simulation on the condition with and without dyke on the sides of the river shows the difference between those results. From the hydraulic simulations using HEC-RAS  it is found that water surface elevations without dyke is 3.49 m (Figure 4a) and with dyke is 3.79 m (Figure 4b). The difference of water surface elevations with and without dyke is 0.3 m.
This study again shows that the change of storage area has an impact on flood peaks. In the existing condition (without dyke) with significantly wide flood plain due to swamp area on the sides of the river, the ratio of contributing catchment area to storage area is about four. When the dyke is considered to protect the area near the river, the ratio of contributing catchment area to storage area increases for nearly four times.
Building a dyke will protect the corresponding area on the sides of the river from flooding. In those two case studies, people who live in the area with annual flooding problem as well as the plantation will be protected from flooding. However, building a dyke may create a new problem. Water level in the river will increase and propagate so that it may cause flooding in another area elsewhere which has no dyke protection.
1. This study shows that building a dyke on the sides of a river will decrease the storage area and increase the ratio of Ac/As. Consequently water level in the river will increase. The extent of water increase in the river is affected by the change of Ac/As ratio. Despite protecting the corresponding area from flooding, the increase of water level in the river may cause flooding elsewhere.
2. This study highlights the importance of the ratio of Ac/As. When ratio of Ac/As is small, then direct rainfall on the storage (lake/swamp) is the dominant external input. However, their probability of triggering flood will be enhanced if ratio of Ac/As is larger, since there is increased chance for the storage capacity deficit to be exceeded by a combination of direct rainfall on the storage and runoff contribution from the upstream catchment area. This clearly points to the importance of Ac/As as a critical parameter governing the frequency and magnitude of storage overflow.
3. The improved understanding of the process control on catchment-storage landscape will be useful toward the estimation of flood and in assisting the management of the combined catchment-storage system for Nagara River and Way Pegadungan catchments and other similar system in the same region and regions elsewhere which contain swamps and/or lakes.
[1.] Balek, J., Hydrology and Water Resources in Tropical Regions, Elsevier, New York, 1983.
[2.] Fitzgibbon, J., and Dunne, T., Land Surface and Lake Storage During Snowmelt Runoff in A Subarctic System, Arct. Alp. Res., 13, 1981, pp. 277-285.
[3.] Mielko, c., and Woo, M.K., Snowmelt Runoff Processes in A Headwater Lake and Its Catchment, Subarctic Canadian Shield, Hydrol. Process, 20, 2006, pp. 987 - 1000.
[4.] Kusumastuti, D.I., Sivapalan, M., Struthers, I., Reynolds, D.A., Murray, K., and Turlach, B.A., Thresholds in the Storm Response of A Catchment-Lake System and the Occurrence and Magnitude of Lake Overflows: Implications for Flood Frequency. Water Resour. Res, 44, 2008, W02438, doi:10.1029/2006WR005628.
[5.] Sivapalan, M., Wood, E. F., and Beven, K.: On Hydrologic Similarity: 3. A Dimensionless Flood Frequency Model Using A Generalized Geomorphologic Unit Hydrograph and Partial Area Runoff Generation, Water Resour Res, 26(1), 1990, pp. 43-58.
[6.] Sivapalan, M., Bloschl, G., Merz, R., and Gutknecht, D., Linking Flood Frequency to Long-term Water Balance: Incorporating the Effects of Seasonality. Water Resour. Res., 41, 2005, W06012, doi:10.1029/2004WR003439.
[7.] Kusumastuti, D.I., Struthers, I., Sivapalan, M., and Reynolds, D.A., Threshold Effects in Catchment Storm Response and the Occurrence and Magnitude of Flood Events: Implications for Flood Frequency, Hydrol. Earth Syst. Sci., 11, 2007, pp. 1515-1528.
[8.] Kusumastuti, D.I., Sivapalan, M., Struthers, I., Reynolds, D.A., Thresholds in the Storm Response of A Lake Chain System and the Occurrence and Magnitude of Lake Overflows: Implications for Flood Frequency. Advances in Water Resources, 31, 2008, pp. 1651-1661.
[9.] Kusumastuti, D.I., The Impact of Within-Storm Temporal Pattern on Flood Frequency; Application on A Catchment in Lampung. International Journal of Science Engineering and Technology, 2(3), 2009, pp. 1-6.
[10.] Kusumastuti, D.I., The Impact of Rainfall Variability and Hydrological Regimes on Flood Frequency, Proceeding International Seminar on Water Related Risk Management, Jakarta-Indonesia, July 15 - 17, 2011, 114-122.
[11.] Kusumastuti, D.I., Konsep Tampungan dan Perlindungan Banjir. Prosiding Seminar Hasil Penelitian & Pengabdian Kepada Masyarakat, Lembaga Penelitian. Universitas Lampung. 2008, pp. 55-61.
[12.] Brunner, G.W., HEC-RAS River Analysis User's Manual Version 4.1. US Army Corps of Engineers, Institute for Water Resources, Hydrology Engineering Center (HEC), 2010.
[13.] Duflow Modelling Studio, User's Guide, Stowa, EDS, 2000.
[14.] Joleha, J., V. Maino, A. Adyabadam, W. Zhao, S. Ke, T.H., Le., Flood Mitigation of Nyando River Using Duflow Modelling, civil Engineering Dimension, 11(1), 2009, pp. 46-57.
[15.] Badilla, R.A., Flood Modelling in Pasiq-Mariikina River Basin. M.Sc. Thesis. 2008. International Institute for Geo-Information Science and Earth Observation, Enschede, The Netherlands.
[16.] Kusumastuti, D.I., Analisis Dampak Tanggul terhadap Elevasi Banjir Sungai Nagara Kalimantan Selatan. Rekayasa, 15(3), 2011, pp. 183-192. Kusumastuti, D.I. (1)
Kusumastuti, D.I. (1)
(1) Civil Engineering Department, Lampung University (UNILA), Jl. Sumantri Brojonegoro No. 1 Bandar Lampung 35145, INDONESIA Email: firstname.lastname@example.org
Note: Discussion is expected before November, 1st 2012, and will be published in the "Civil Engineering Dimension" volume 15, number 1, March 2013.
Received 17 November 2011; revised 01 February 2012; accepted 30 May 2012
Table 1. Rainfall-runoff Model Parameters Parameter Value Units a 70 [mm.sup.0.5] [hr.sup.0.5] b 0.5 [S.sub.b] 150 mm [S.sub.fc] 45 mm M 0.6 [S.sub.bl] 1500 mm [E.sub.p] 1700 mm/year Tabel 2. Some Parameters Describing the Catchment of Way Pegadungan Parameter Value catchment area 1200,58 [km.sup.2] River length 36,06 km River slope 0,0167% Number of first order branch 8 Number of all order branches 48 Length of first order branch 71,67 km Length of all order branches 91,11 km
|Printer friendly Cite/link Email Feedback|
|Publication:||Civil Engineering Dimension|
|Date:||Sep 1, 2012|
|Previous Article:||Application of prefabricated vertical drain in soil improvement.|
|Next Article:||Soil liquefaction in Padang due to Padang earthquake 30 September 2009.|