**Journal of Water Resource and Protection**

Vol.10 No.02(2018), Article ID:82717,15 pages

10.4236/jwarp.2018.102010

Evaluating Discharge Capacity of Major Chara’s of Sylhet City Using GIS

Gulam Md. Munna^{*}, Jahir Bin Alam, Jobayer Ahmed, Abdullah Al Mamun^{ }

Department of Civil and Environmental Engineering, Shahjalal University of Science and Technology (SUST), Kumargaon, Sylhet-3114, Bangladesh

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 7, 2017; Accepted: February 25, 2018; Published: February 28, 2018

ABSTRACT

Heavy rainfall is one of the most frequent and widespread severe weather hazards that affect Sylhet city. In Sylhet, nine natural channels locally called “chara” are mainly responsible for the drainage of heavy rainfall to the Surma River. The present condition of this natural drainage is not feasible due to the unplanned land development and manipulation by people. So water logging is a common scenario of Sylhet city. Many parts of Sylhet city are experiencing the severe inundation problem due to heavy rainfall. Time series analysis of rainfall data (1974-2015) is important for knowing and predicting of rainfall variation. Mann-Kendall analysis showed no trend for monthly and yearly rainfall data. The rainfall intensity of Sylhet is higher than other districts of Bangladesh. In this study, IDF (Intensity Duration Frequency) curve has been developed that is commonly used in engineering planning and design. By using IDF curve, ArcGIS software, DEM map, and normal discharge, maximum discharge for a different point of Mongoli chara and Bolramer chara has been calculated (for 25 years return period). Discharge through Mongoli chara and Bolramer chara has been represented in ArcMap using ordinary kriging. To keep the catchment of those chara free from inundation and water logging problem, this calculated discharge needed to be managed.

**Keywords:**

ArcGIS, Canal, Chi-Square, Gamble Distribution, IDF Curve, Mann-Kendall, Rainfall

1. Introduction

Bangladesh is affected by various types of natural disasters almost every year because of the global warming as well as climate change impacts. Bangladesh is a low-lying country with more than 230 waterways, is one of the most disaster-prone nations in the world [1] . Fifteen percent of its land floods annually on average [2] . As the lowest riparian in a huge trans-boundary river basin, Bangladesh faces increasing threat of massive flood exposure [3] .

Floods are a natural and certain event. Natural flood management includes techniques and process that aim to work with natural hydrological and morphological processes, features and characteristics to deal with the sources and pathways of flood waters.

The outcome of urbanization is broadly known as hydraulic studies since there have been increasing problems in urban storm water drainage management. One of the most important facilities in preserving and improving the urban environment is an adequate and properly functioning storm water drainage system which includes storm water movement and storage capacity.

Numerous catchments in Bangladesh are currently under serious pressure from urban, industrial and infrastructure development where downstream receiving water bodies such as river, lakes, channels, pond and reservoirs have turned out to be touchy to build rates and volumes of overflow and contamination release.

Inadequate drainage causes needless death and disease and loss of home, property and live hoods. Poor storm water management also pollutes the environment and squanders limited freshwater resource. A sound investment in storm water can reduce those losses, but only after setting priorities and making different choices [4] [5] .

Integrated storm water is still a moderately new idea in many countries including Bangladesh. Present experience indicates that end-of-pipe, rapid disposal, localized reactive and mono-functional drainage concepts have been generally practiced in Bangladesh. The local design makes and professionals are beginning to recognize the need for a new and broader approach to urban drainage, flood control and storm water management in the light of development in the country processing at a tremendous stress.

2. Study Area

Our study areas are Mongoli Chara and Bolramer Chara in Sylhet City (Figure 1) of Bangladesh. The local design makes and professionals are beginning to recognize the need for a new and broader approach to urban drainage, flood control and storm water management in the light of development in the country processing at a tremendous stress.

3. Materials and Methods

The study area contains many rainfall stations of Bangladesh Meteorological Department (BMD). The daily rainfall data was collected from BMD for the years 1973 up to 2015 (Total 41 years). For the continuity of the data, the entire record of gauging station was verified.

3.1. Trend Analysis

Mann-Kendall trend test is a nonparametric test used to identify a trend in a

Figure 1. Mongoli Chara and Bolramer Chara in Sylhet, Bangladesh. Mongoli Chara: It flows from Mirer Moidan road to Old Medical, cross the Chouhatta-Rikabi bazaar road and flows southward to Soros pur, Jallar par, Kazir bazaar, Sheik Ghat and Taltota. Finally the chara falls into Surma River crossing the Kazir bazaar-Tookhana road. Bolramer Chara: Bolramer Chara is flowing from Jamtala (south of Jallar par road), Bondar bazaar and Taltota. One part of it meets with Mongoli Chara near the Mohan Market and another part falls on Surma River near Kazir bazaar.

series. The purpose of the Mann-Kendall (MK) test [6] is to statistically assess if there is a monotonic increase or decrease trend of the rainfall of 1974 to 2015 [7] .

3.2. Development of Short Duration Rainfall Data

Short duration data are rare in development countries like Bangladesh. The data of 24-hour rainfall records are available in rain gauge stations.

Indian Meteorological Department (IMD) uses an empirical reduction formula [8] .

${P}_{t}={P}_{24}{\left[t/24\right]}^{1/3}$ (1)

where, P_{t} is the required rainfall depth in mm at t hr. duration, P_{24} is the daily rainfall in mm and t is the duration of rainfall for which the rainfall depth is required in hr.

For estimation of various duration like 1-hr rainfall values from annual maximum daily rainfall data, Chowdhury et al. [8] , used Indian Meteorological Department (IMD) empirical reduction formula to estimate the short duration rainfall from daily rainfall data in Sylhet city and found that this formula gives the best estimation of short duration rainfall [9] . In this study, this empirical formula was used to estimate the short duration rainfall in SCC.

3.3. Probability Distribution and Chi-Square Test

To identify a specific theoretical distribution for the available data it is important to find an acceptable method. The aim of the test is to find how good a fit is the observed and the predicted data. Chi-square is one of the most widely used tests to find the best fit theoretical distribution of any specific dataset which is represented by this formula [10] .

${X}^{2}={\displaystyle {\sum}_{i=1}^{n}\frac{{\left({O}_{i}-{E}_{i}\right)}^{2}}{{E}_{i}}}$ (2)

where, ${O}_{i}$ and ${E}_{i}$ represent the observed and expected frequencies respectively. If the observed frequencies are close to the corresponding expected frequencies, the ${X}^{2}$ value will be small, indicating a good fit; otherwise it will be a poor fit. Chi-Square test was performed by sigma magic software.

3.4. Short Duration for Gumbel Distribution

Short duration rainfall data such as 5 min, 10 min, 15 min, 30 min, 45 min, 120 min have been developed using the depth-duration formula [10] .

$\frac{{P}_{t}^{T}}{{P}_{60}^{T}}=0.54{t}^{0.25}-0.50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(5\le t\le 120\mathrm{min}\right)$ (3)

$\frac{{P}_{t}^{T}}{{P}_{t}^{10}}=0.21\mathrm{ln}T-0.52\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(2\le T\le 100\text{\hspace{0.17em}}\text{years}\right)$ (4)

where, ${P}_{t}^{T}$ is the depth of t-minute, T-year return period rainfall, ${P}_{t}^{10}$ is the depth of a t-minute, 10-year return period rainfall and ${P}_{60}^{T}$ is the depth of a 60-minute, T-year return period rainfall.

3.5. Development of Rainfall Intensity Duration Frequency (IDF) Curve

IDF curve is used for estimating the maximum rainfall intensity for the different duration and return period. A specified return period T and duration t is calculated for the rainfall depth. Its mean intensity ${I}_{m}$ (mm/hour) has been obtained dividing it by the duration t (hour). Then the IDF curve is obtain by plotting, on a graph, the main intensity ${I}_{m}$ (mm/hour) against the duration D (min).

Mathematically, this curve can be represented in different forms as follows [9] :

$i=\frac{a}{{\left(b+D\right)}^{n}}$ (5)

$i=x\ast {\left(D\right)}^{-y}$ (6)

where, i is rainfall intensity in mm/hr, D is duration of rainfall in minutes, the parameters a, b, n, x, y define the shape and appropriate units for curve fitting to IDF data.

3.6. Catchment Area of Chara

The catchment area of Mongoli Chara and Bolramer Chara has been identified through field survey and contour map. The direction of runoff flow has been identified during field survey and consulting local people. And then the catchment area was divided into some regions.

3.7. Runoff Coefficient

The runoff coefficient C is the variable of the rational method. The runoff coefficient (C) is a dimensionless coefficient relating the amount of runoff to the amount of precipitation received. It is a larger value for areas with low infiltration and high runoff, and lower for permeable, well-vegetated areas [11] . Where watershed is not homogeneous, a weighted runoff coefficient should be determined. A weighted runoff coefficient is computed using the following equation [9] :

${C}_{w}=\frac{{\displaystyle {\sum}_{j=1}^{n}{C}_{j}{A}_{j}}}{{\displaystyle {\sum}_{j=1}^{n}{A}_{j}}}$ (7)

where, ${A}_{j}$ is the area for land cover j, ${C}_{j}$ is the runoff coefficient for area j, N is the number of distinct land covers within the watershed, and ${C}_{w}$ is the weighted runoff coefficient.

3.8. Time of Concentration

Time of concentration ( ${T}_{c}$ ) is the time required for runoff to travel from the hydraulically most distant point in the watershed to the outlet. The hydraulically most distant point is the point with the longest travel time to the watershed outlet, and not necessarily the point with the longest flow distance to the outlet. Time of concentration is generally applied only to surface runoff and may be computed using many different methods. Time of concentration will vary depending upon slope and character of the watershed and the flow path characteristics [12] [13] .

The U.S. Federal Aviation Agency (FAA) developed a simple estimation of ${T}_{c}$ to be used with the Rational Method [14] . Time of concentration ( ${T}_{c}$ ) is an important variable in many hydrologic methods, including the Rational and Natural Resources Conservation Service procedures. For the same size watershed, the shorter the ${T}_{c}$ , the larger the peak discharge [14] .

${T}_{c}=\frac{1.8\ast \left(1.1-C\right)\ast \sqrt{L}}{\sqrt[3]{S}}$ (8)

where, C is the cover factor, L is the hydraulic length in ft, and S is the slope in Percent.

3.9. Discharge Due to Rainfall

The rational formula has developed the relation between rainfall and discharge. Discharge due to rainfall calculated by this formula [15] .

${Q}_{r}=\frac{C\ast i\ast A}{360}$ (9)

where,
${Q}_{r}$ : Peak discharge (m^{3}/s) C: Runoff Coefficient.

i: Rainfall intensity during time of concentration (mm/hr).

A: Drainage Area (ha).

3.10 Discharge Due to Population

There are many formulas for population prediction. The geometric mean increase is used to find out the future increment in population. Assuming that, the future population is increasing at a constant rate. Geometric increase formula [16] is:

${p}_{n}={p}_{o}{\left(1+\frac{r}{100}\right)}^{n}$ (10)

Here, ${p}_{n}$ = population after n year, ${p}_{o}$ = Present population, r = rate of growth and n = number of year.

Discharge Due to Population [14] :

${Q}_{P}=X\ast P$ (11)

where, ${Q}_{P}$ = Discharge due to population, X = Per capita discharge and P = Total population.

4. Result and Discussion

4.1. Monthly Rainfall Trend Analysis

There are several statistical tests available for testing stationeries of time series. In the present study, the Man-Kendell test (Table 1) for linear trend has been carried out for the monthly rainfall series from 1974 to 2015 of Sylhet.

Table 1. Trend analysis.

There is no trend observed in monthly rainfall data series to find the best fit that the observed frequencies are near to the corresponding expected frequencies.

4.2. Chi-Square Test

Chi-Square test has been performed by 1-hr rainfall data that shown in Table 2.

Table 2. Chi-square test.

Beta distribution, Cauchy distribution, Gamma distribution, Rayleigh distribution, Logistic distribution, Gumbel distribution fits the data. Chi-Sq. value is low for Gumbel distribution, so the Gumbel distribution best fits the selected data.

4.3 Gumbel Distribution

1-hour 2, 5, 10, 25, 50, 100-year rainfall depth for Sylhet is estimated using Gumbel distribution formula which has been show in Table 3.

Table 3. 1-hour rainfall depths at different rainfall periods.

Rainfall depth depends on return period. While return period is 2 year rainfall depth is low. But with the increase of return period rainfall depth increase simultaneously. In 2 year return period rainfall depth is 66.847 mm and in 100 year return period rainfall depth is 123.82 mm.

4.4. Development of Rainfall Intensity Duration

Within same return period rainfall depth depends on rainfall duration. Rainfall depth is increased with the increase of rainfall duration. In 2 year return period for 5 min and 120 min rainfall depth is 20.55474 mm and 86.05012 mm respectively. Average intensities of rainfall for different return periods have been calculated by dividing t minute rainfall depths by the corresponding duration of t (hour) (Table 4).

Table 4. Rainfall depths of Sylhet for t-minute.

4.5. Intensity-Duration-Frequency (IDF) Curve

Short duration rainfall data such as 5 min, 10 min, 15 min, 30 min, 45 min, 120 min for 2, 5, 10, 25, 50, 100-year have been developed using the depth-duration formula. Within same return period, rainfall intensity depends on rainfall duration. Rainfall Intensity decrease with the increase of rainfall duration. In 2 year return period for 5 min and 120 min rainfall Intensity are 246.6569 mm and 43.02506 mm respectively. For the same duration, rainfall intensity depends on rainfall duration. Rainfall intensity is increasing with the increase of return period. For 30 min duration 2 years and 100 years rainfall intensity is 102.1142 mm and 189.1556 mm. Intensity-Duration-Frequency (IDF) curve has been developed by plotting average intensities against the duration of rainfall. IDF curve for Sylhet has been shown in Figure 2.

For same year return period, rainfall intensity decrease with the increase of time duration. For same time duration rainfall intensity is higher for higher return period. It was observed in the study that the IDF empirical equation form

$i=x\ast {\left(D\right)}^{-y}$ had the best fit rather than the equation form $=\frac{a}{{\left(b+D\right)}^{n}}$ .

4.6. Derivation of IDF Equation

To derive an equation for calculating rainfall intensity for the regions of interest

Figure 2. IDF curve for Sylhet city.

or the station, there are some required steps for establishing an equation suit the calculation of rainfall intensity for a certain recurrence interval and specific rainfall period which depends mainly on the results obtained from the intensity duration frequency (IDF) curves and the corresponding logarithmic conversion, where it is possible to convert the equation into a power equation (Table 5), and thus to calculate all the parameters related to the equation.

Table 5. Rainfall IDF curve equation and regression value for Sylhet city.

100 year rainfall intensity is higher than 2 year rainfall intensity. Rainfall intensity equation format is
$i=x\ast {\left(D\right)}^{-y}$ . R^{2} value is same for every rerun period.

4.7. Catchment Area Selection & Dividing into Sub-Areas

Catchment Area of was selected by using “Polygon” tools and then converted it to “Shape file” by using “Feature to Polygon” tools of ArcGIS. Separation into sub-area was done by the same process (Figure 3).

4.8. Development of Digital Elevation Model (DEM)

From Elevation map of Sylhet city DEM was created by using “geo-statistical analysis” tools of ArcGIS (Figure 4).

Figure 3. Catchment area of Chara.

4.9. Time of Concentration

Time of concentration ( ${T}_{c}$ ) of six zone was estimated by Federal Aviation Agency (FAA) formula (Table 6). Length of flow and average slope was found from ArcMap. Length of flow and average slope is calculated from Shape file and DEM map. Time of concentration is highest at zone 3 and lowest at zone 6.

Table 6. Estimation of time of concentration using FAA formula.

Figure 4. DEM of Sylhet city.

4.10. Discharge Due to Rainfall

Discharge was calculated by the rational formula. Rainfall intensity has been calculated from IDF curve for the duration equal to the time of concentration. 25 year return period has been selected for calculating rainfall calculated discharge has been given in below. Drainage area is calculated from Shape file. At zone 1 discharge is maximum as it has largest area, at zone 5 discharge is minimum as its catchment area is lowest. Peak discharge is the sum of discharge of this zone and peak discharge of previous zone. Peak discharge at zone 4 is the sum of discharge of zone 4, zone 6 and peak discharge of zone 3.

4.11. Discharge Due to Population

Averagely per capita use of water 100 - 120 liter/day (PHE) has been considered. Meantime, discharge due to population was measured. Table 7 shows at zone 1 discharge is maximum, at zone 5 discharges is minimum. Peak discharge is the sum of discharge of this zone and peak discharge of previous zone. Peak discharge at zone 4 is the sum of discharge of zone 4, zone 6 and peak discharge of zone 3 (Table 8).

Table 7. Discharge from different zone.

Table 8. Calculated peak discharge.

4.12. Total Discharge

Total discharge is the sum of discharge due to rainfall and discharge due to the population. Table 9 shows the peak discharge for different zones, which is maximum (60.87 m^{3}/s) at zone 5 and minimum (8.58 m^{3}/s) at zone 6.

Table 9. Design discharge.

4.13. Discharge in Chara

From the equation the corresponding discharge at different point of Chara in m^{3}/s was obtained and plot them in the ArcMap software by using “kriging” tools of ArcGIS. After that, Geostatic analysis was carried out with Kriging which gave the following map (Figure 5).

Figure 5. Discharge in Chara.

4.14. Existing and Required X-Sectional Area

The required cross-sectional area has been calculated which is shown in Table 10. Before meeting with Bolramer chara the existing x-section of Mongoli chara satisfy the required cross section. But after meeting with Bolramer chara the existing cross section is not sufficient.

Table 10. Existing & required x-sectional area.

5. Conclusions

Based on the result of the study, following conclusions were made:

・ There is no trend in rainfall that could be detected at the 5% significance level in Sylhet.

・ Average intensities of rainfall for 25 years return periods have been calculated. Average intensities for 5 min, 10 min, 15 min, 30 min, 60 min, 120 min are 387.2996, 289.8691, 236.2573, 160.3392, 104.9632, 67.55776 (mm/hr) respectively.

・ Existing cross section area of chara at Kazir Bazar (before meeting with Bolramer Chara) & Kazir Bazar (after meeting with Bolramer Chara) has been found 25.5 m^{2}, 34.96 m^{2} and required X-sectional area 21 m^{2}, 37.5 m^{2} respectively. Peak discharges of Mongoli Chara and Bolramer Chara have been calculated to 60.87 m^{3}/s and 8.58 m^{3}/s respectively for 25 years return period.

・ The required cross-sectional area has been calculated. Before meeting with Bolramer Chara, the existing x-section of Mongoli chara satisfies the required cross section. But after meeting with Bolramer Chara, the existing cross section is not sufficient.

Cite this paper

Munna, G.M., Alam, J.B., Ahmed, J. and Al Mamun, A. (2018) Evaluating Discharge Capacity of Major Chara’s of Sylhet City Using GIS. Journal of Water Resource and Protection, 10, 167-181. https://doi.org/10.4236/jwarp.2018.102010

References

- 1. Development, GBPHOCI. DFID’s Programme in Bangladesh. 2010.
- 2. Memon, N. (2012) Disasters in South Asia: A Regional Perspective. Pakistan Institute of Labour Education & Research.
- 3. Mozumder, P. (2010) Exploring Flood Mitigation Strategies in Bangladesh.
- 4. Parkinson, J. and Mark, O. (2005) Urban Stormwater Management in Developing Countries. 1st Edition, IWA Publishing, London, UK.
- 5. Kar, K.K., Yang, S.K., Lee, J.H. and Khadim, F.K. (2017) Regional Frequency Analysis for Consecutive Hour Rainfall Using L-Moments Approach in Jeju Island, Korea. Geoenvironmental Disasters, 4, 1-13. https://doi.org/10.1186/s40677-017-0082-0
- 6. Mann, H.B. (1945) Nonparametric Test against Trend. Econometrica, 13, 245-259. https://doi.org/10.2307/1907187
- 7. Mishra, N. and Gupta, S. (2015) Precipitation Trend Analysis by MK Test in Garhwal Region Districts of Uttarakhand (1901-2010). Journal of Space Science & Technology, 4, 25-35.
- 8. Chowdhury, R.K., Alam, M.J.B., Das, P. and Alam, M.A. (2007) Short Duration Rainfall Estimation of Sylhet: IMD and USWB Method. Journal-Indian Waterworks Association, 39, 285.
- 9. Rashid, M.M., Faruque, S.B. and Alam, J.B. (2012) Modelling of Short Duration Rainfall Intensity Duration Frequency (SDRIDF) Equation for Sylhet City in Bangladesh. ARPN Journal of Science and Technology, 2, 92-95.
- 10. Nhat, L.M., Tachikawa, Y. and Takara, K. (2006) Establishment of Intensity-Duration-Frequency Curves for Precipitation in the Monsoon Area of Vietnam. Annuals of Disas. Prev. Res. Inst., Kyoto Univ., No. 49 B, 2006. http://www.dpri.kyoto-u.ac.jp/nenpo/no49/49b0/a49b0p09.pdf
- 11. Thompson, D.B. (2006) The Rational Method.
- 12. Solomon, K.R., Baker, D.B., Richards, R.P., et al. (1996) Ecological Risk Assessment of Atrazine in North American Surface Waters. Environmental Toxicology and Chemistry, 15, 31-76. https://doi.org/10.1002/etc.5620150105
- 13. Kar, K.K., Yang, S.K. and Lee, J.H. (2015) Assessing Unit Hydrograph Parameters and Peak Runoff Responses from Storm Rainfall Events: A Case Study in Hancheon Basin of Jeju Island. Journal of Environmental Science International, 24, 437-447. https://doi.org/10.5322/JESI.2015.24.4.437
- 14. Matthew, A.P. (2009) Engineering 101: The FAA Equation-Time of Concentration. http://www.civil4d.com/2009/03/engineering-101-the-faa-equation-time-of-concentration/
- 15. Wong, T.S. (2002) Call for Awakenings in Storm Drainage Design. Journal of Hydrologic Engineering, 7, 1-2. https://doi.org/10.1061/(ASCE)1084-0699(2002)7:1(1)
- 16. Garg, S.K. (2005) Hydrology and Water Resources Engineering. 13th Revised Edition, Khanna Publisher, Delhi, India.