Stochastic Modeling of Flow and Sediment of the Rivers at Delta Head, East Coast of India ()
1. Introduction
The Central Water Commission (CWC) of India has classified rivers in India based on their Catchment area as major (>20,000 sqkm), medium (2000 - 20,000 sqkm) and minor (<2000 sqkm). India has 12 major, 46 medium and 14 minor and desert flowing river basins comprising of areas 252.8 MH, 24.6 MH and 110 MH respectively (CWC Hand Book 2005) [1] . All peninsular rivers are either west or east flowing in India. The major east flowing rivers are the Subarnarekha, the Mahanadi, the Godavari, the Krishna, the Cauvery and the Vaigai. Large, medium and small prominent rivers debouching in Bay of Bengal are about 88 numbers (Table 1). The length of east coast (EC) is 2729 km from the river Hooghly (north) to Tambaraparani (south) (Figure 1).
Figure 1. Major basins in east coast of India. Source: NRSC, Bhuvan, Hyd.
Table 1. The rivers and the lagoons that transport sediment to Bay of Bengal.
The Subarnarekha, the northern most east flowing peninsular river, is originating from rifts of the Chota Nagpur Hills (lat 23˚18'E and long 85˚11'N). The small river, Thamarparani in extreme south starts from the Western Ghats (lat 8˚46'E and long 77˚15'N) near Vana-theertham waterfalls near Agastyarkoodam peak of Western Ghats in Ambasamudram district.
The east flowing rivers originate from eastern side of the Western Ghats range and flow through the basalts of the Deccan plateau. They take course along the Deccan traps and crosses through the gorges of the Eastern Ghats belt to join the Bay of Bengal.
The lengths of east flowing ephemeral rivers are long and catchment areas of these basins are large with wide deltas than that of West Coast Rivers. The east coast is emergent whereas the west coast submerging. They carry huge sediment to form large deltas. The longest river is the Godavari (1465 km) of catchment area 312,812 km2 (Figure 2) which carries average annual flow of 83.94 Bcum at deltahead (till 1986-2011) and average sediment influx from the basin (45.5 MMT till 1993-2012) (Table 2).
East coast of India is extending from Kanya Kumari in south to Haladia port in north with a continental shelf of area 22,411 sqkm up to 50 m., Kulkarni et al., (1985) [2] . It has five coastal states having coastal length i.e. West Bengal (158 km), Orissa (480 km), Andhra Pradesh (974 km) and Tamil Nadu (1076 km) and Pondicherry Union Territory (41 km). Rivers joining the Bay of Bengal in West Bengal are either from NW or north. The mega cities along East Coast are Calcutta, Vishakhapatnam and Chennai etc. Indian rivers carry 1645 × 109 km3
Table 2. The major 12 rivers along EC of India and their length, area, time series for sediment and flow (source: CWC Integrated Hand book 2005, 2006, 2009, 2012 and 2015) [30] .
Krishna Basin Godavari Basin Mahanadi Basin
Figure 2. Three large basins of the rivers Mahanadi, the Godavari and Krishna in east coast of India.
of water with huge quantities of sediment in average 1.2 × 1012 kg to sea which is one tenth of sediment flux of the globe Chandra Mohan et al. (2001) [3] .
The climate of the region is subtropical to tropical. The temperature of the basin area is in average 18˚C to 45˚C with rainfall ranging from 1400 - 1500 mm in northern part and 900 - 1000 mm in the southern region of EC of India. The South West monsoon (JJAS months) and North East monsoon (Nov. and Dec.) influence the quantity of rainfall of the area. The rivers in south of the Godavari basin receives rainfall from both the south west monsoon and north east monsoon. The east flowing rivers to the north of the Godavari basin receives rainfall from the south west monsoon. The Bay of Bengal disturbances through the Odisha coast (48.2%), West Bengal coast (25.5%, Andhra coast (19.3%) and Tamilnadu coast (6.7%) Mohapatra et al., 2010 [4] .
The East coast of India and the adjoining Bay of Bengal was rifted from Antarctica during the early Cretaceous period (65.5 Mya BP). The conjugate oceanic crust of the Bay of Bengal was assumed to be between Ender Bay basin and East Antarctica. Fuloria et al. 2014 [5] , described a major fault separate the Ganga basin from the Mahanadi basin. The Visakhapatnam high separates the Mahanadi from the Godavari and the highs of Naidupeta and Chingle put between the Krishna and the Cauvery basins in the peninsular India (Figure 2). The rivers took thousands of years to manifest by deposition in layers of clastic sediments till extinction of prolonged volcanic activities.
Tandon et al., 2007 [6] reported the high sediment flux in the peninsular rivers in India are not the outcrop of recent tectonic setting like Himalayan rivers but due to the result of tropical weathering of rocks in their basin. Lal et al. 2009 [7] , identified four stages of formation of the east coast of India which was Cretaceous origin but large modifications occurred during Quaternary period due to global cooling and warming. Dune dynamics and reactivation along coast have changed the river dynamics from late Holocene to recent. Deltaic Changes of river course were identified geospatially in Odisha, Andhra and Tamilnadu coasts (Mahalik et al. (2006) [8] , Kumaran et al. (2012) [9] , Kakani et al. (2003) [10] , Prabakaran et al. 2010 [11] and Sundaresh et al., 2014 [12] .
Review of Literature
The suspended sediment (coarse, medium and silt loads) are 90% of the total annual sediment load which is carried by the flood flow d/s to Bay retaining only 10% behind the water resources structures. Dams (Major, Medium and minor) and BWA’s (barrages, weirs and anicuts) in Odisha constructed along these rivers are 180 and 43 respectively (INDIA-WRIS WIKI-2014) [13] . Increased anthropogenic activities with natural processes have accelerated the sediment fluxes in the sea Milliman et al. 1992 [14] , Syvitsky et al. (1999) [15] . The net effect of sediment load changes on the global land-ocean sediment flux suggests that the “natural” or “pre-human” flux has been almost doubled by anthropogenic impacts. But 50% of the gross flux has been trapped in reservoirs, Zhou Zhide et al. 2015 [16] .
The deltas of rivers are sinking, shrinking and under subsidence due to anthropologic interventions like dams and BWA’s. Sediment studies have been done by various authors like Gamage et al. (2009) [17] , Rao et al. (2010) [18] Syvitski et al. (2009) [19] , (2011) [20] , Gupta et al. (2012) [21] and Dandekar (2014) [22] for only east flowing large rivers in India, like the Ganges, Mahanadi, the Godavari, the Krishna and the Cauvery but not for rivers of small catchment. The coastal plains of E.C. of India grow deltas by the sediments of the major rivers, but some patches are of Aeolian and marine origin. There are seven important deltas in the east coast, namely the Subarnarekha, the Mahanadi, the Godavari, the Krishna, the Penner and the Cauvery host estuaries with dense mangroves reported by Kumaran et al., 2012 [8] .
Sanil Ku. V. et al., 2006 [23] reported that shoreline along east coast oscillate and seasonal. The places where the cycle of one oscillation is not complete within a year the beach suffer from net erosion and vice versa. Ports and harbours suffer from erosion to its north as a result of the construction of breakwaters as sedimentation is not allowed.
Nnaji et al. 2014 [24] , reported Gumbel and Log Pearson Type III (LPT III) and Weibull distributions can reasonably predict of flow frequency. LPT III distribution can be provided and more appropriate. Vogel R. A, (1993) [24] reported LPT III method of flood flow modeling was recommended by United states as flood frequency guide lines. The Gumbel method of prob. distribution model accepted by chi-square goodness of fit tests was considered as the best model in the Seyhan river Teheran, Topaloglu (2002) [25] , Luo P. et al., 2013 [26] classified rivers into three parts basing on the quantum of suspended sediment and the found the R2 values of the time series. Solomon et al., 2013 [27] , considered Gumbel method of distribution was suitable for predicting expected flow of the river OSSE, Argentina. Ghosh K. G. et al. 2015 [28] reported that after construction of dam and barrages flood frequency have raised and analysis of the flood frequency by LPT III method is acceptable. Mishra et al. 2015 [29] , reported Gumbel and LPT III methods can be used for probabilistic prediction of monsoon rainfall in the Mahanadi basin.
Less research work have been done on the sediment flow of small river basin. This paper is an attempt to study the water and the sediment flow in 12 major rivers flowing into the Bay of Bengal along east coast of India at their delta head.
2. Methods
Spatiotemporal trends of annual discharge of (monsoon + non-monsoon) flow and sediment were collected from Integrated Hydrological data book (non-classified basins), 2005, 2006, 2009, 2012 and 2015 of Central Water Commission, Government of India [30] . The annual flow and corresponding sediment data was considered of the stations at the delta head of the twelve rivers. The data for the major 12 major and medium rivers debouching Bay of Bengal were compiled as a time series and analyzed (Figure 3). The rate of suspended
Figure 3. Rivers and their Gauge sites along east coast of India as per Table 1.
sediment discharge against annual flow and their trend with time were studied. The linear and nonlinear regression equations of the twelve rivers were developed using different curves such as polynomial, sigmoid, peak, power, logarithm, rational and wave form by the statistical tools and EXCEL methods and others. The equations that had highest R2 value (>60%) were considered as the fitted curve.
Flood and sediment modeling and prediction were done by deterministic modeling by using a long data series. Different models used are deterministic, stochastic, statistical or soft computing. Attempt has been made in this paper to predict the annual flow and corresponding sediment at various recurrence intervals though the time series length is short. The common statistical probability density functions used were Gumble, Log Pearson Type III (LPT III) and compared the results with their present trend.
3. Methodology
The annual time series data were collected to find a trend in the annual flow and the sediment influx and used to predict future values of the time series. Regression equations both linear and nonlinear equations used were polynomials, peak, sigmoid, wave form, power, exponential and rational etc. Correlation is done to test how closely the experimental data (Variables) are interdependent. The correlation coefficient and coefficient of determination were found to verify the reliability of the data.
Correlation Coefficient (R)
It is the data in a time series are linearly related the correlation coefficient (r) is calculated by formulae
It is the ratio of standard deviation to covariance and how a time series close to the mean value of data series,
where
and
where “d” represents as the deviation from the mean. R values lies between −1 < R < 1. When R is zero or near zero there is almost zero correlation, +1 is a perfect fit with same slope and the time series is positively and linearly related. For R-values to be negative it indicates perfect negative fit with opposite slope and the series is not related. OR, the other form can be: The
where SSE is the sum of squares due to error and
where,
= mean value of the series y.
Coefficient of determination (R2)
In a regression model the ratio of the variance of the fitted values to the observed values are called coefficient of determination. If
is the observed value,
= mean value of the series y and
is the fitted value then the coefficient of
determination =
where, R2 has values lies between 0 and 1.
Higher the R2 values better the result of the model and hence the prediction and the curve fitting. It indicates about % of data lies close to best fit line. Negative values of R2 indicate the regression equation does not fit the data.
The dimensions like shape, size, and flow parameters differ. The regression equations that fit the curve with highest R2 values also differ from river to river. The common used statistical projections for annual flow and sediment were made by using Gumbel II or Log Pearson Type III methods for recurrence intervals 1.11, 2, 5, 10, 20, 25, 50, 100 ∙∙∙ upto 10,000 years for 12 major rivers in Odisha, Andhra Pradesh, Puducherry and Tamil Nadu.
4. Discussion
The rivers, their length, the location of the gauge and discharge site, the average annual flow and the average annual sediment of the twelve rivers are given in Table 2.
The total catchment area for the peninsular east flowing rivers debouching the Bay of Bengal upto their delta head found to be 0.97 million sqkm. The average annual flow, the corresponding average annual sediment flowing at the delta head of the rivers found to be 227.3 Billon Cum and 84.6 MMT respectively (About 10% allowances provided for the small rivers and lagoons which are not considered).
5. Regression Equations
5.1. Annual Flow Series (Curve Fitting)
The annual flow for the largest and the longest 12small and large rivers in East coast of India were taken for the years 1986 to 2011. Various stochastic models tried and the curve fitting equations having highest R2 values are given in Table 3.
Using methods of both linear and nonlinear regression techniques, fitted curves were tried. The equation of the curve considered to be the best fit whose coefficient of determination values (R2) found to be the highest. From time series of annual flow for the northern east flowing rivers in India (the Subarnarekha to the Godavari) supports rational function whereas the small southern rivers (Krishna to Thambarparani) follows Gaussian or Lognormal distribution except the river Vaigai, The fittest distribution have high R2 values for the southern rivers than the northern rivers Table 3.
5.2. Annual Sediment Series (Curve Fitting)
The annual sediment flow for the largest and the longest 12small and large rivers in East coast of India were taken for the years 1993 to 2012. Various stochastic models tried and the curve fitting equations having highest R2 values for the quantum of sediment flow in the rivers are given in Table 4.
On linear and nonlinear curve fitting, the annual sediment flow exhibit nonlinear correlation of Gaussian type and the coefficient of determination values have shown higher values 50% - 99% except rivers Brahmani and Godavari. Maximum numbers of large dams in the sub- basins of the river Godavari might not be giving R2-value <50% (Table 5). Both the rivers have surplus water potential considered by NWDA, the National Water Development agency (an agency to interconnect rivers of India).
6. Estimation of Flow and Sediment
Return period or recurrence interval (Tr) is the average interval of any event in the years equaling or exceeding a given magnitude. Probability exceedence (p) is
Table 3. The model fit to the time series for annual flow of the rivers (monsoon + non-monsoon) to Bay of Bengal (BoB) via East coast of India.
Table 4. The model fit to the time series for sediment influx to BoB via the rivers in East Coast.
Table 5. Estimation of annual total flow in east flowing rivers in EC of India (Gumbel II method).
the chance that the annual event of any year will equal or exceed some given value. The probability (P) of an event with recurrence interval T is P = 1/T.
6.1. Gumbel Method
Gumbel method is one of the oldest methods for prediction of hydrologic data Al-Mashidani et al. 1978 [31] . This method is given by:
Where
= the discharge at recurrence interval of T-years,
and
are the mean of the observed time series discharge and standard deviation respectively. K is the frequency factor found from Gumbel’s table, where
where
is Euler’s constant = 0.5772. To test the goodness of fit χ2 test is conducted with hypothesis that the data fits Gumbel distribution and vice versa. (Never Mujere 2011) [32] . The estimations are done by using statistical Hydraulic package http://onlinecalc.sdsu.edu/textbookhydrologyp217.html.
From the future estimation of quantity of flow by Gumble II method in 12 major rivers in east coast of India at various recurrence intervals are observed. It is found that total annual discharge of the major rivers shall increase proportionately at higher return periods. There shall be years ahead when the monsoon shall be erratic and the rainfall/runoff shall have increasing tendency as flow prediction indicate irregularity Table 6.
Prediction of the annual sediment against various recurrence interval of the east flowing rivers in the east coast, it is ascertained that the sediment concentration/Bcum of inland inflow shall increase gradually and the delta building process shall be emerging Table 7.
Table 6. Estimation of sediment in east flowing rivers in East coast of India (Gumbel II method).
Table 7. Estimation of annual flow to BoB in major East Coast Rivers of India (LPT III method).
6.2. LPT III Method
LPT III method is one of the commonly used method hydrologic frequency analysis. The method is extremely flexible Nnaji G. A. et al., 2014 [33] . When skewness coefficient become zero when the LPT III distribution reduced to lognormal distribution. LPT III distribution is calculated using equation
where x is taken as annual flow or yearly sediment flow.
arithmetic mean of logarithm values of the variables x and K is the frequency factor found from the Log Pearson frequency factor table, and is the standard deviation of the
logarithm values given by
and the results are shown in Table 7.
On analysis of annual flow series of major rivers in the east coast falling Bay of Bengal by LPT III it is observed that at 50 years return period Godavari shall increase in flow and shall have 108.84 Bcum whereas present annual average discharge was 83.942 Bcum. The river Brahmani shall show a constant average annual flow for the coming 10,000 years Figure 4 and Table 8.
Table 8. Estimation of sediment influx to BoBin MMT different east flowing rivers in East coast of India (LPT III method).
Figure 4. Total annual flow from major rivers in east coast of India (Gumble II and LPT III).
On comparison of predicted values at various recurrence intervals of total annual flow through the major rivers along east coast falling in the Bay of Bengal it is observed that Gumbel II method had lower values than LPT III method towards tail end. The forecasted values are almost same upto 4000 years of forecast Figure 5.
On comparison of the projected values of total annual sediment influx from the inland rivers along east coast, it is found Gumbel II method predict lower values of sediment at lower probabilities whereas LPT III method very high values. For design purposes of hydraulic structures Gumbel method may be accepted as projected sediment Figure 5.
Figure 5. Total annual sediment influx from major rivers of EC of India to BoB (Gumble II and LPT III) [Estimation].
7. Conclusions
On study of major east flowing rivers along east coast of India and falling in Bay of Bengal, the following conclusions may be drawn:
1) All the east flowing rivers are emerging below the rift of Narmada valley and falling in Bay of Bengal forming large coastal plains, number of distributaries and encouraging huge sediment flow to the Bay of Bengal. The average annual flow and sediment at the delta head of the rivers are found to be 227.3 Billon Cum and 84.6 MMT respectively.
2) Erratic monsoon and hydrologic interventions have reduced sediment concentration of the rivers reaching delta head.
3) Regression analysis results of annual flow and sediment series had made a distinct division of peninsula Indian rivers above and below the Godavari River. The rivers flowing south of the river Godavari followed Gaussian model and the rivers above towards North East followed rational equations. The flow and the sediment influx are monsoon dependent.
4) Regression analysis of sediment influx time series follows Gaussian pattern towards south of Godavari river, but to the north, the sediment pattern does not obey any curve pattern, and 60% follow peak type equations except sediment flow of the Brahmani river.
5) The flow and sediment were projected by Gumbel II and Log Pearson III methods and observed both the methods gave almost equal predicted values for flow but for sediment forecast Gumbel II method gave lower values at various return periods.
From the above results, it can be inferred that for any structural design purposes, the Gumbel method projected values can be used for hydrologic study.