Investigating the Spatial and Temporal Variability of Precipitation using Entropy Theory
Abstract
This study uses entropy theory to develop a novel application of the apportionment entropy disorder index (AEDI) to capture both spatial and temporal variability in monthly precipitation for various types of hydrologic modeling. In total, 41 Environment Canada stations across Ontario with long term (1955 to 2005) records and a very low percentage of missing data were selected. It was found that the fall and summer seasons are the major contributors to annual precipitation variability. Spatial variability of annual precipitation was observed to be increasing from southern to northern Ontario. The AEDI index map of Ontario, developed in this study, has been successfully integrated into several hydrologic models.
1 Introduction
Prediction of precipitation variability and its influence on the hydrological cycle and water quality is an increasingly active area in hydrological research (Zhang et al. 2000; Zhang et al. 2001; Lemmen et al. 2008; Finney et al. 2010; Ahmed et al. 2013; Disley et al. 2015; Liu et al. 2015; Mattar et al. 2016; Atieh et al. 2017). The distribution of precipitation spatially and temporally directly affects the hydrological cycle within a watershed, altering flow regimes and complicating the accurate predictions of flow (Das et al. 2008; Boyer et al. 2010; Coulibaly 2006; Khan and Coulibaly 2010; Liu and Cui 2011; Trenouth et al. 2013; Trenouth and Gharabaghi 2015; Trenouth et al. 2015; Thompson et al. 2016). An understanding of precipitation variability is also needed for flood and drought preparedness planning. There is a need to better understand the spatial and temporal variability in precipitation to more accurately model the hydrologic system and the health of aquatic life (Perera et al. 2009; Perera et al. 2010; Perera et al. 2013; Betts et al. 2014; Betts et al. 2015; Sabouri et al. 2013; Trenouth and Gharabaghi 2015; Gazendam et al. 2016). Parameters that represent this variability may be important factors to include in modeling studies (Vasiljevic et al. 2012; Rudra et al. 2015).
Understanding precipitation variability is a necessity for water resource management as any change in pattern influences the availability of water resources. As such, studies have used different methods that assess trends in hydrological variables. These methods may be divided into bootstrap (Burn et al. 2010), regression (Cheung et al. 2008), and nonparametric methods (Tabari et al. 2015). Trend and variation analysis of hydrologic parameters including flow and precipitation has commonly been studied by applying statistical analysis methods such as linear regression (Akinremi et al. 1999; Asnaashari et al. 2015; Atieh, Gharabaghi et al. 2015; Atieh, Mehltretter et al. 2015). However, the expected accuracy of linear regression analysis is limited because the relationships between hydrologic parameters are generally nonlinear. The nonparametric Mann–Kendall (MK) test has been the most popularly used for trend analysis (Gajbhiye et al. 2016; Palizdan et al. 2014; Sayemuzzaman and Jha 2014). The major limitation to the MK test is its sensitivity to serial correlation in time series data which is available in precipitation data (Palizdan et al. 2014; Sayemuzzaman and Jha 2014). Time series analysis has been the most popular method for hydrological trend analysis as it does not assume linearity (Sang 2013; Yu et al. 2014). The Fourier transform has been widely adopted as a method of distinguishing the high frequency modes within a time series (Fotiadi et al. 1999). Major drawbacks of using the Fourier transform are that it is not well adapted to non-stationary data and does not maintain temporal information (Drago and Boxall 2002; Coulibaly 2006; Palizdan et al. 2014). More recently, wavelet analysis has been adopted for analysing hydrologic parameter trends and variations because it overcomes these limitations of the Fourier transform (Coulibaly 2006; Adamowski et al. 2009; Nalley et al. 2012; Adamowski et al. 2013). Discrete wavelet analyses have also been coupled with the MK test to determine the dominant periodic components affecting precipitation trends (Palizdan et al. 2014). Application of the concept of entropy in water resource engineering is an evolving idea that has attracted more research in the past decade. The application of entropy concepts has facilitated advancements in a range of hydrologic problems including the assessment of model performance, parameter estimation, derivation of functional relationships, development of flow duration curves, streamflow forecasting, uncertainty estimation, and assessment of the efficiency of monitoring networks (Cui and Singh 2016; Hao and Singh 2011; Moramarco and Singh 2010; Singh 1997; Singh 2013; Singh et al. 2014). Liu et al. (2008) studied the basin-scale spatial patterns of precipitation time series data using Shannon‘s entropy theory in the Yellow River basin, China. Entropy was used to detect the spatial structure of precipitation that was found to increase with longitude. Mishra et al. (2009) developed a methodology using entropy theory to investigate the intra-annual and decadal distributions of precipitation across Texas. They were able to determine the spatial variability of precipitation amounts within a year across Texas. Precipitation variability, quantified using entropy, is also an important parameter to include in various hydrologic modeling studies. Atieh, Gharabaghi et al. (2015) and Atieh, Mehltretter et al. (2015) developed an entropy parameter that captures both the temporal and spatial variability of precipitation in models for predicting flow duration curves and sediment parameters at ungauged sites, and found that including the entropy parameter improved predictions. Using entropy concepts for estimating the variability of hydrologic parameters is still a novel and appealing option that has not been fully investigated.
1.1 Entropy
Entropy can be considered as a measure of uncertainty associated with a random process, in this case precipitation, and may be described by a probability function. Reduction of uncertainty is equal to increasing the amount of information gained from a set of observations and thus entropy can be used to indirectly quantify the information content of a given data time series. The concept of Boltzmann entropy is a result of the second law of thermodynamics. Maxwell’s law for the distribution of velocities among molecules of a gas (1859) was used as a starting point for Boltzmann’s research into determining a statistical representation of entropy. Initially, Boltzmann indicated that in addition to understanding the velocity distribution of gas molecules, a time-independent energy distribution among any large scale collection of molecules at equilibrium should also exist (Dugdale 1996). The Boltzmann fundamental entropy equation is:
(1) |
where:
k | = | Boltzmann constant with a value of 1.32 × 10−23J/K, and |
r | = | multiplicity, i.e. the number of microstates that could occur in the system being considered. |
In a large scale system, the natural logarithm of r is equal to the natural logarithm of the maximum probability (Dugdale 1996).
Shannon entropy, also called informational entropy, was first introduced by Shannon in 1948. In 1997, Singh provided an extensive review of the uses of entropy, including Shannon entropy, in hydrology and water resources. The concept underlying Shannon entropy is that the more uncertain an event, the more information is needed to characterize it. Entropy concepts can be used to measure the uncertainty of a collection of events by determining the average amount of information needed to eliminate uncertainty (Singh 2013).
This study adopts entropy theory as a means to investigate the temporal and spatial variability of precipitation in Ontario, Canada. The uncertainty associated with different time series of precipitation measurements was determined and variation in precipitation was analysed relative to other stations in the area. Precipitation analysis is a common need in the water resource sector, and is necessary for determining the effect of annual, seasonal and decadal variability on the availability of water. This information directly impacts future water resource decision-making. The key objective of this study was to determine temporal and spatial variability of precipitation in Ontario on annual and decadal scales.
2 Methodologies
In this study, temporal variability is measured as the randomness of a time series (annual and decadal) over different time intervals (monthly, seasonal, and yearly). Variability is defined by a disorder index that is the difference between maximum and observed entropy within an individual series. A higher disorder index indicates greater variability.
Precipitation time series were considered individually to determine variability within each time series (temporal variability). Variability was then compared between stations (spatial variability). The time series used in this study are mean daily precipitation, annual precipitation and decadal precipitation time series. The variability of the annual time series was analysed over monthly intervals. The variability of the decadal time series was analysed for monthly, seasonal, and annual time intervals separately. The winter season was considered to include the months of December, January and February. Further, spatial variability in annual and decadal precipitation across Ontario was addressed by analysing average annual and decadal measurements respectively at each location.
2.1 Apportionment Entropy (AE)
Apportionment entropy (AE) is a measure of the temporal variability of monthly precipitation over a year.
(2) | |
(3) |
The aggregate precipitation in a year (A) is the sum of the individual aggregate precipitation during each month (ai). AE is highest (AE = log212 = 3.585) when annual precipitation is equally apportioned between each month (probability 1/12 for each month), while it is lowest (AE = 0) when all precipitation occurs in a single month.
2.2 Decadal Apportionment Entropy (DAE)
Decadal apportionment entropy (DAE) measures the randomness of precipitation over a time interval of a decade:
(4) | |
(5) |
The aggregate decadal precipitation (D) is the sum of the individual aggregate annual precipitations for each year (di). DAE quantifies the variability of the decadal annual precipitation over a 10 y period. DAE is highest when annual precipitation is equally apportioned to each year (DAE = log210). Similarly, Equations 4 and 5 may be applied to any time series (monthly, seasonal). For example, Equation 4 could be applied to the aggregate precipitation values for the month of January for a decade, which would allow the determination of the randomness that occurred for that month over the decade.
2.3 Disorder Index (DI)
A disorder index is the entropy calculated, based on the observed precipitation, divided by the maximum possible entropy (even distribution of precipitation throughout the time interval). The disorder index may be calculated using an annual or a decadal dataset, and is thus labelled the annual apportionment entropy disorder index (AEDI) and decadal apportionment entropy disorder index (DEDI) respectively. The disorder index is a measure of the variability within a time series. A mean disorder index is used to compare temporal variability (i.e. from one year to another) and spatial variability (i.e. from one station to another). A long term mean is used as a threshold value to determine high versus low variability.
3 Dataset
This project focused on the province of Ontario, located in east central Canada with a total area of 1 000 000 km2. Precipitation is lowest in northwest Ontario and increases in a southeast direction; annual mean precipitation in northern Ontario ranges between 400 mm and 600 mm whereas in the southern parts of Ontario it is between 800 mm and 1 200 mm (Atlas of Canada 2014).
Using the Environment Canada Adjusted Precipitation Dataset stations located in Ontario with total monthly precipitation data available for the 50 y period 1955 to 2005 were selected. This period is similar to that adopted by Nalley et al. (2012) for studying the variability of precipitation in Ontario and Quebec. Kumar et al. (2009) also indicated this period was sufficient for analyzing long term precipitation trends. The 41 stations with the lowest percentages of missing data were selected. Figure 1 shows the location of the selected stations within Ontario. Stations are numbered from the most southern station (station 1) to the most northern one (station 41). Of the 41 stations, 30 had <4% missing data for the 1955 to 2005 period, whereas the highest percentage of missing data among the 41 stations was 22%. This was included in order to allow a wide distribution of stations in northern Ontario. Mishra and Singh (2010) indicated that when missing precipitation data is less than 3%, the data can be considered of good quality. However, they also included in their study a few stations with missing values >3% (Mishra and Singh 2010).
Figure 1 Distribution of selected precipitation stations in Ontario.
4 Results and Discussion
4.1 Statistical Measures
Using ArcGIS 10.1 and Environment Canada Climate Normals (1971 to 2000), a total of 245 monitoring stations with 30 y of data were selected in Ontario and surrounding provinces to generate a mean annual precipitation (MAP) map using the inverse distance weighted method (Atieh, Gharabaghi et al. 2015). Figure 2 shows that the MAP varies between 500 mm and 1 200 mm across Ontario with the lowest precipitation (500 mm to 700mm) in northern Ontario. Standard deviation (SD) of the mean daily precipitation was calculated for all stations and was represented as a percentage of MAP. It ranged between a minimum of 22% (northern Ontario) to a maximum of 42% (southern Ontario).
Figure 2 Mean annual precipitation in Ontario in mm (Atieh, Gharabaghi et al. 2015).
4.2 Variability of Precipitation Distribution Within a Year
Using apportionment entropy, the distribution of annual variability of precipitation in a year was investigated during a 50 y period (1955 to 2005). The apportionment entropy disorder index (AEDI) was used as an indicator of variability of precipitation distribution for different months and different seasons within a year. The higher the AEDI value, the higher the variation within the annual time series. The AEDI was first calculated for each station individually to study temporal variation, after which the stations were compared to determine the spatial distribution of precipitation variability. The degree of variability for specific years and stations can be identified based on a threshold for temporal and spatial variation respectively.
Temporal Variability of Annual Time Series
Apportionment entropy (AE) was calculated at the 41 sites for each month over a year using Equation 2. The AEDI value (Equation 4) was then calculated using the difference between maximum possible apportionment entropy and observed apportionment entropy. The mean AEDI of the stations over 50 y was computed to examine average temporal variation in the region. The long term average of mean AEDI was also computed and used as a threshold indicator to identify the years with high precipitation variability versus those with low variability over a regional basis. Results presented in Figure 3 indicate that within the variability range 0.18 to 0.21, there appear to be a few high and low peaks. Years 1993 and 2000 had the highest annual variability, with a few smaller peaks in 1968, 1977 and 1987. The lowest variation was observed in years 1972, 1975, 1990, 1996, and 2003.
Figure 3 Mean AEDI for all stations using monthly time interval.
Spatial Variability of Annual Time Series
The means of the computed AEDI for all years across the selected 41 stations for annual and seasonal time series are presented in Figures 4 and 5 respectively to show the spatial variability of annual precipitation time series in Ontario. The long term mean is used as the threshold to determine stations with high variability versus those with low variability.
Figures 4 and 5 both show that stations 1 to 28 were associated with lower variability (both monthly and seasonal) whereas the remaining stations had higher variability. The highest variability was observed at the westernmost stations in Ontario. The trend of spatial variability is similar for both time scales over an annual time series with a few distinctions. An interesting observation is that variability of precipitation within a year increases from southern Ontario to northern Ontario (Figure 6).
Figure 6 Distribution of apportionment entropy disorder index across Ontario.
4.3 Variability of Precipitation Within a Decade
The decadal apportionment entropy disorder index (DEDI) was calculated for each station to investigate the spatial and temporal variability. DEDI was calculated for five decades (1955 to 2005). Results are presented in Figure 7. For the monthly scale, there is a large variation between decades. The highest peak was observed in the decade 1956 to 1965, for the month of October. A high variability peak is observed in February for all decades. These peaks in variability of precipitation may be due to greater ranges of snowfall in February of different years and similar patterns of wide ranges of rainfall in October in Ontario.
Figure 7 Mean DEDI for all stations over 5 decades.
4.4 Applications of AEDI in Hydrologic Modeling
The methodology used in this paper is based on entropy, which is a different perspective on analyzing hydrological trends than has been seen in other papers that used statistical and signal analysis. Using AEDI the variability of the annual and seasonal time series was analyzed. The years with higher variability are determined and any trends throughout the analysed 50 y may be identified. Through collapsing the time scale and analyzing AEDI spatially, across the 41 stations, the spatial distribution of precipitation variability was presented. Similarly, using DEDI, the decadal trends in precipitation was analyzed.
Incorporating precipitation variability in various hydrological models allows a more accurate presentation of the hydrologic system (Hwang et al. 2012; Masinde 2014; Rudra et al. 2015). Different strategies are used to represent precipitation variability within hydrologic models, such as comparing long term precipitation records and dividing precipitation time series into different periods (Masinde 2014; Hiltner et al. 2016; Li et al. 2016). The entropy based methodology presented in this paper generates a simple, easy to use parameter (AEDI) that allows incorporation of temporal and spatial variability of precipitation into any hydrologic model. An AEDI map may be developed for any region in the world using the methodology presented in this paper for generating a spatial AEDI map of Ontario. Atieh, Gharabaghi et al. (2015) and Atieh, Mehltretter et al. (2015) were the first to introduce AEDI as an input parameter to an artificial neural network (ANN) model to more accurately predict flow duration curve (FDC) and sediment rating curve (SRC) parameters at ungauged sites. The results of Atieh, Gharabaghi et al. (2015) indicated that the parameters defining the log normal distribution for the FDCs (location and scale) were highly sensitive to AEDI. Incorporating AEDI in the location and scale ANN models improved prediction performance by 7% and 21%, respectively. Atieh, Mehltretter et al. (2015) reported that removing climatic parameters from the SRC prediction model decreased the correlation coefficient by 40% and the Nash–Sutcliffe coefficient by 42%. This paper recommends using AEDI in future hydrological modeling research such as runoff estimation, flood and drought prediction, and water quality studies.
5 Conclusions
This study presents an alternative method for quantifying spatial and temporal variability of precipitation for hydrologic modeling. The apportionment entropy disorder index (AEDI) was calculated using 50 y (1955 to 2005) historic precipitation records for 41 Environment Canada Stations across Ontario. The AEDI values of the stations were compared to determine the spatial distribution of precipitation variability. Precipitation variability was observed to increase from southern Ontario to northern Ontario.
The decadal entropy disorder index (DEDI) was calculated for 5 decades (1955 to 2005). The highest monthly variability peaks were observed in February. Spatial variability increased from southern Ontario to northern Ontario.
This paper presents a novel approach for a better understanding of the spatial and temporal variability of precipitation which is useful for a variety of hydrologic modeling projects, including flood and drought preparedness planning.
Acknowledgments
The authors would like to acknowledge the generous support of the Natural Sciences and Engineering Research Council of Canada (NSERC) for their support of this research. The authors also acknowledge the support received from Dr Coulibaly and Dr Mishra, including constructive comments and the initial Matlab code.
References
- Adamowski, J., K. Adamowski and A. Prokoph. 2013. “Quantifying the Spatial Temporal Variability of Annual Streamflow and Meteorological Changes in Eastern Ontario and Southwestern Quebec Using Wavelet Analysis and GIS.” Journal of Hydrology 499:27–40.
- Adamowski, K., A. Prokoph and J. Adamowski. 2009. “Development of a New Method of Wavelet Aided Trend Detection and Estimation.” Hydrological Processes 23 (18): 2686–96.
- Ahmed, S., A. Singh, R. Rudra and B. Gharabaghi. 2013. “Comparison of CANWET and HSPF for Water Budget and Water Quality Modeling in Rural Ontario.” Water Quality Research Journal of Canada 49 (1): 53–71.
- Akinremi, O., S. Mcginn and H. Cutforth. 1999. “Precipitation Trends in the Canadian Prairies.” Journal of Climatology 12:2996–3003.
- Asnaashari, A., B. Gharabaghi, E. McBean and A. Mahboubi. 2015. “Reservoir Management under Predictable Climate Variability and Change.” Journal of Water and Climate Change 6 (3): 472–85.
- Atieh, M., B. Gharabaghi and R. Rudra. 2015. “Entropy-Based Neural Networks Model for Flow Duration Curves at Ungauged Sites.” Journal of Hydrology 529 (3): 1007–20.
- Atieh, M., S. Mehltretter, B. Gharabaghi and R. Rudra. 2015. “Integrated Neural Networks Model for Prediction of Sediment Rating Curve Parameters for Ungauged Basins.” Journal of Hydrology 531 (3): 1095–107.
- Atieh, M., G. Taylor, A. M. Sattar and B. Gharabaghi. 2017. “Prediction of Flow Duration Curves for Ungauged Basins.” Journal of Hydrology 545:383–394. https://doi.org/10.1016/j.jhydrol.2016.12.048.
- Atlas of Canada. 2014. Atlas of Canada. Ottawa: Natural Resources Canada. http://atlas.nrcan.gc.ca/site/english/index.html.
- Betts, A. R., B. Gharabaghi and E. A. McBean. 2014. “Salt Vulnerability Assessment Methodology for Urban Streams.” Journal of Hydrology 517:877–88.
- Betts, A. R., B. Gharabaghi, E. A. McBean, J. Levison and B. Parker. 2015. “Salt Vulnerability Assessment Methodology for Municipal Supply Wells.” Journal of Hydrology 531 (3): 523–33.
- Boyer, C., D. Chaumont, I. Chartier and A. G. Roy. 2010. “Impact of Climate Change on the Hydrology of St. Lawrence Tributaries.” Journal of Hydrometeorology 384 (1–2): 65–83.
- Burn, D. H., M. Sharif and K. Zhang. 2010. ”Detection of Trends in Hydrological Extremes for Canadian Watersheds.” Hydrological Processes 24:1781–90.
- Coulibaly, P. 2006. “Spatial and Temporal Variability of Canadian Seasonal Precipitation (1900–2000).” Advances in Water Resources 29:1846–65.
- Cheung ,W. H., G. B. Senay and A. Singh. 2008. “Trends and Spatial Distribution of Annual and Seasonal Rainfall in Ethiopia.” International Journal of Climatology 28:1723–34.
- Cui, H. and V. Singh. 2016. “Maximum Entropy Spectral Analysis for Streamflow Forecasting.” Physica A: Statistical Mechanics and its Applications 442:91–9.
- Das, S., R. Rudra, B. Gharabaghi, S. Gebremeskel, P. K. Goel and W. T. Dickinson. 2008. “Applicability of AnnAGNPS for Ontario Conditions.” Canadian Biosystems Engineering, 50 (1): 1–11.
- Disley, T., B. Gharabaghi, A. Mahboubi and E. McBean. 2015. “Predictive Equation for Longitudinal Dispersion Coefficient.” Hydrological Processes 29 (2): 161–72.
- Drago, A. F. and S. R. Boxall. 2002. “Use of the Wavelet Transform on Hydro-Meteorological Data.” Physics and Chemistry of the Earth, Parts A/B/C 27 (32–34): 1387–99. https://doi.org/10.1016/S1474-7065(02)00076-1.
- Dugdale, J. S. 1996. Entropy and its Physical Meaning. Abingdon: Taylor and Francis.
- Finney, K., B. Gharabaghi, E. A. McBean, R. P. Rudra and G. MacMillan. 2010. “Compost Biofilters for Highway Stormwater Runoff Treatment.” Water Quality Research Journal of Canada, 45 (4): 391–402.
- Fotiadi. A. K., D. A. Metaxas and A. Bartzokas. 1999. ”A Statistical Study of Precipitation in Northwest Greece.” International Journal of Climatology 19:1221–32.
- Gajbhiye, S., C. Meshram, S. K. Singh, P. K. Srivastava and K. Islam. 2016. “Precipitation Trend Analysis of Sindh River Basin, India from 102 Year Record (1901–2002).” Atmospheric Science Letters 17:71–7.
- Gazendam, E., B. Gharabaghi, J. Ackerman and H. Whiteley. 2016. “Integrative Neural Networks Models for Stream Assessment in Restoration Projects.” Journal of Hydrology 536:339–50.
- Hao, P. and V. P. Singh. 2011. “Single-Site Monthly Stream Flow Simulation Using Entropy Theory.” Water Resource Research 47 (9): W09528
- Hiltner, U., A. Brauninga, A. Gebrekirstos, A. Huth and R. Fischer. 2016. “Impacts of Precipitation Variability on the Dynamics of a Dry Tropical Montane Forest.” Ecological Modeling 320:92–101
- Hwang, Y., M. Clark, B. Rajagopalan and G. Leavesley. 2012. “Spatial Interpolation Schemes of Daily Precipitation for Hydrologic Modeling.” Stochastic Environmental Research Risk Assessment 26:295–320.
- Khan, M. and P. Coulibaly. 2010. “Assessing Hydrologic Impact of Climate Change with Uncertainty Estimates: Bayesian Neural Network Approach.” Journal of Hydrology 11:482–95.
- Kumar, S., V. Merwade, J. Kam and K. Thurner. 2009. ”Streamflow Trends in Indiana: Effects of Long Term Persistence, Precipitation and Subsurface Drains.” Journal of Hydrology 374 (1–2): 171–83.
- Lemmen, D., F. Warren and J. Lacroix. 2008. “Synthesis.” In From Impacts to Adaptation: Canada in a Changing Climate 2007, edited by D. Lemmen, F. Warren, J. Lacroix and E. Bush, 1–20. Ottawa: National Resources Canada.
- Li, J., G. Li, S. Zhou and F. Chen. 2016. “Quantifying the Effects of Land Surface Change on Annual Runoff Considering Precipitation Variability by SWAT.” Water Resources Management 30 (3): 1071–84
- Liu, Q. and B. Cui. 2011. “Impacts of Climate Change/Variability on the Streamflow in the Yellow River Basin, China.” Ecological Modeling 222 (2): 268–74.
- Liu, Q., Z. Yang and B. Cui. 2008.” Spatial and Temporal Variability of Annual Precipitation During 1961–2006 in Yellow River Basin, China.” Journal of Hydrology 361 (3-4): 330–8.
- Liu, Y., W. Yang, Z. Yu, I. Lung and B. Gharabaghi. 2015. “Estimating Sediment Yield from Upland and Channel Erosion at a Watershed Scale Using SWAT.” Water Resources Management 29 (5): 1399–412.
- Masinde, M. 2014. “Artificial Neural Networks Models for Predicting Effective Drought Index: Factoring Effects of Rainfall Variability.” Mitigation and Adaptation Strategies for Global Change 19 (8): 1139–62.
- Mattar, M. A., A. A. Alazba, B. Alblewi, B. Gharabaghi and M. A. Yassin. 2016. “Evaluating and Calibrating Reference Evapotranspiration Models Using Water Balance under Hyper-Arid Environment.” Water Resources Management 30:3745.
- Mishra, A. K., O. Mehmet and V. P. Singh. 2009. “An Entropy Based Investigation into the Variability of Precipitation.” Journal of Hydrology 270:139–54.
- Mishra, A. K. and V. P. Singh. 2010. “Changes in Extreme Precipitation in Texas.” Journal of Geophysical Research 115 (D14): D14106.
- Moramarco, T. and V. P. Singh. 2010. “Formulation of the Entropy Parameter Based on Hydraulic and Geometric Characteristics of River Cross Sections.” Journal of Hydrologic Engineering 15 (10): 852–8.
- Nalley, D., J. Adamowski and B. Khalil. 2012. “Using Discrete Wavelet Transforms to Analyse Trends in Streamflow and Precipitation in Quebec and Ontario (1954–2008).” Journal of Hydrology 475:204–28.
- Palizdan, N., Y. Falamarzi, F. Y. Huang, T. S. Lee and A. Ghazali. 2014. “Regional Precipitation Trend Analysis at the Langat River Basin, Selangor, Malaysia.” Theory of Applied Climatology 117:589–606.
- Perera, N., B. Gharabaghi and K. Howard. 2013. “Groundwater Chloride Response in the Highland Creek Watershed Due to Road Salt Application: A Re-Assessment After 20 Years.” Journal of Hydrology 479:159–68.
- Perera, N., B. Gharabaghi and P. Noehammer. 2009. “Stream Chloride Monitoring Program of City of Toronto: Implications of Road Salt Application.” Water Quality Research Journal of Canada 44 (2): 132–40.
- Perera, N., B. Gharabaghi, P. Noehammer and B. Kilgour. 2010. “Road Salt Application in Highland Creek Watershed, Toronto, Ontario.” Water Quality Research Journal of Canada 45 (4): 451–61.
- Rudra, R., W. T. Dickinson, S. I. Ahmed, P. Patel, J. Zhou, B. Gharabaghi and A. A. Khan. 2015. “Changes in Rainfall Extremes in Ontario.” International Journal of Environmental Research 9 (4): 1117–26.
- Sabouri, F., B. Gharabaghi, A. Mahboubi and E. McBean. 2013. “Impervious Surfaces and Sewer Pipe Effects on Stormwater Runoff Temperature.” Journal of Hydrology 502:10–7.
- Sang, Y. F. 2013. “A Review on the Applications of Wavelet Transform In Hydrology Time Series Analysis.” Atmospheric Research 122:8–15.
- Sayemuzzaman, M. and M. Jha. 2014.” Seasonal and Annual Precipitation Time Series Trend Analysis in North Carolina, United States.” Atmospheric Research 137:183–94.
- Singh, V. P. 1997. “The Use of Entropy in Hydrology and Water Resources.” Hydrological Processes 11:587–626.
- Singh V. P. 2013. Entropy Theory and its Applications in Environmental and Water Engineering. New York: Wiley.
- Singh V. P., A. Byrd and H. Cui. 2014. “Flow Duration Curve Using Entropy Theory.” Journal of Hydrological Engineering 19 (7): 1340–8.
- Tabari, H., M. T. Taye and P. Willems. 2015. ”Statistical Assessment of Precipitation Trends in the Upper Blue Nile River Basin.” Stochastic Environmental Research Risk Assessment 29:1751. https://doi.ORG/10.1007/s00477-015-1046-0.
- Thompson, J., A. Sattar, B. Gharabaghi and R. Warner. 2016. “Event-Based Total Suspended Sediment Particle Size Distribution Model.” Journal of Hydrology, 536: 236-246.
- Trenouth, W. and B. Gharabaghi. 2015. “Event-Based Design Tool for Construction Site Erosion and Sediment Controls.” Journal of Hydrology 528:790–5.
- Trenouth, W., B. Gharabaghi, A. Bradford and G. MacMillan. 2013. “Better Management of Construction Sites to Protect Inland Waters.” Inland Waters 3 (2): 167–78.
- Trenouth, W., B. Gharabaghi and N. Perera. 2015. “Road Salt Application Planning Tool for Winter De-Icing Operations.” Journal of Hydrology 524:401–10.
- Vasiljevic, B., E. McBean and B. Gharabaghi. 2012. “Trends in Rainfall Intensity for Stormwater Designs in Ontario.” Journal of Water and Climate Change 3 (1): 1–10.
- Yu, Z. G., Y. Leung, Y. D. Chen, Q. Zhang, V. Anh and Y. Zhou. 2014. ”Multifractal Analyses of Daily Rainfall Time Series in Pearl River Basin of China.” Physica A: Statistical Mechanics and its Applications 405:193–202.
- Zhang, X., K. D. Harvey, W. D. Hogg and T. Yuzyk. 2001. “Trends in Canadian Streamflow.” Water Resources Research 37:987–98.
- Zhang, X., L. A. Vincent, W. D. Hogg and A. Niitsoo. 2000. “Temperature and Precipitation Trends in Canada During the 20th Century.” Atmosphere-Ocean 38:395–429.