1. About SepHydro

SepHydro is a hydrograph separation (or baseflow separation) web tool developed through a collaborative research effort between Canadian Rivers Institute (CRI), University of New Brunswick (UNB), Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC). SepHydro is a result of the research effort aimed at evaluating the effects of agricultural production systems on groundwater and surface water quality and on the quality of downgradient aquatic ecosystems.
SepHydro is part of the Hydrology Tool Set (HTS; https://www.hydrotools.tech), a suite of tools that can be used for advancing the understanding of various local and watershed scale hydrological processes. HTS also includes ETCalc (potential, reference and actual evapotranspiration), SWIB (crop water deficit and excess, irrigation requirement and its impact on aquifer storage, water budget), SNOSWAB (water balance - snowfall, snowmelt, snowpack, soil water content, evapotranspiration, drainage, infiltration, surface runoff, etc.), RECHARGE BUDDY (groundwater recharge), and SNOWFALL BUDDY (snowfall and rainfall amounts).
The tool offers several customizable filtering algorithms for performing hydrograph separation and assessing surface runoff and groundwater contributions to streamflow, as well as various output and data visualization options. The web-based tool has been designed with a user-friendly interface and provides flexibility in the format of user-data input and output files through a streamlined process.
For detailed instructions on how to use this tool please refer to the User Guide section.
If you used SepHydro, please use the following citation in your publication(s):
Danielescu S, MacQuarrie KTB, Popa A (2018) SEPHYDRO: A Customizable Online Tool for Hydrograph Separation. Groundwater 56: 589-593. DOI: https://doi.org/10.1111/gwat.12792

2. Background

Knowledge of baseflow contribution to streamflow is important for understanding watershed scale hydrology, including groundwater-surface water interactions, impact of geology and landforms on baseflow, estimation of groundwater recharge rates, etc.
Hydrograph separation (also known as baseflow separation) methods can be used as supporting tools in many areas of environmental research, such as the assessment of the impact of agricultural practices, urbanization, climate change etc. on surface water and groundwater.
Conceptually, baseflow is the portion of the streamflow that has a different source than surface runoff. Often times, baseflow is considered to represent the sum of both deep and shallow subsurface contributions to streamflow.
Recursive digital filters are one of the most popular methods used for hydrograph separation. Digital filters, routinely used for signal analysis, allow for separation of “high” (i.e. surface runoff) and “low” (i.e. baseflow) frequency signals using a mathematical formulation. Each digital filter method is dependent on one or more filter parameters, which control the shape and magnitude of baseflow.
One advantage of the digital filters is that they can be applied in multiple passes, and therefore the streamflow can be separated in more than two components. The components can then be attributed to various sources of streamflow such as surface runoff, delayed shallow flow, tile drain flow, deep groundwater flow, etc.

3. Method Overview

Currently, the tool offers a total of 11 methods for the purpose of hydrograph separation. The methods are based on either digital filters or graphical techniques. Brief descriptions of each of the methods are included below.

3.1. Lyne and Hollick Method

This hydrograph separation method was first introduced by Lyne and Hollick in 1979. It has been designed as a means of analyzing the time-dependent and slow stream responses to precipitation. The method assumes high frequency signals to represent direct runoff, while low frequency signals are associated with the baseflow component. The method calculates the surface runoff component which, together with the streamflow, is then used to derive the baseflow component.
Lyne-Hollick Baseflow Separation Method
q – surface or quick runoff [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the surface runoff is calculated;
α - catchment constant [values between 0 and 1];

Lyne V., Hollick M. (1979) - Stochastic time-variable rainfall-runoff modelling. Institute of Engineers Australia National Conference. Publ. 79/10, pp. 89-93.


3.2. Chapman Method

The Chapman method was introduced in 1991 as a response to Lyne-Hollick algorithm's results, which incorrectly provided a constant streamflow or baseflow even when direct runoff has ceased.
Chapman Baseflow Separation Method
b - baseflow [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the baseflow is calculated;
α - hydrological recession constant [values between 0 and 1];

Chapman T. (1991) - Comment on evaluation of automated techniques for base flow and recession analyses, by RJ Nathan and TA McMahon. Water Resources Research, 27(7), pp. 1783-1784.


3.3. Eckhardt Method

In 2004, Eckhardt presented a new two-parameter algorithm for hydrograph separation that was inspired by the previous one-parameter filters.
Eckhardt Baseflow Separation Method
b - baseflow [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the baseflow is calculated;
α – groundwater recession constant [values between 0 and 1];
BFImax - long-term ratio of baseflow to total streamflow [values between 0 and 1];

Eckhardt K. (2005) - How to construct recursive digital filters for baseflow separation. Hydrol. Process. 19, pp. 507–515.


3.4. Pettyjohn and Henning - Fixed Interval Method

This hydrograph separation method is one of the three methods developed by Pettyjohn and Henning in 1979. The method systematically draws connecting lines between the low points of the streamflow hydrograph to determine the baseflow hydrograph. The low points of the streamflow hydrograph are determined using a fixed window of specified width (i.e. equal to a set number of streamflow readings in the source dataset). All baseflow values in a given interval or window are set to the minimum streamflow value in that respective interval. The window width is the nearest integer between 3 and 11 that is equal to 2N, where N is empirically determined as follows:
Pettyjohn and Henning Fixed Interval Baseflow Separation Method
N - number of days after which runoff ceases;
A - drainage area [km2]*;
* the original equation uses [sq mi], however the tool requires [km2] and converts the drainage area to the proper units internally;

Pettyjohn W.A., Henning R. (1979) - Preliminary estimate of ground-water recharge rates, related streamflow and water quality in Ohio: Ohio State University Water Resources Center Project Completion Report Number 552, p. 323.


3.5. Pettyjohn and Henning - Sliding Interval Method

The low points of the streamflow hydrograph are determined using a moving interval or sliding window of specified width (i.e. equal to a set number of streamflow readings in the source dataset). The baseflow value in the middle of each given interval or window is set to the minimum streamflow value in that respective interval. The window width is the nearest integer between 3 and 11 that is equal to 2N, where N is empirically determined as follows:
Pettyjohn and Henning Sliding Interval Baseflow Separation Method
N - number of days after which runoff ceases;
A - drainage area [km2]*;
* the original equation uses [sq mi], however the tool requires [km2] and converts the drainage area to the proper units internally;

Pettyjohn W.A., Henning R. (1979) - Preliminary estimate of ground-water recharge rates, related streamflow and water quality in Ohio: Ohio State University Water Resources Center Project Completion Report Number 552, p. 323.


3.6. Pettyjohn and Henning - Local Minimum Method

The low points of the streamflow hydrograph are determined using a moving interval or sliding window of specified width (i.e. equal to a set number of streamflow readings in the source dataset). Each streamflow value is checked whether it is the lowest in a given interval or not. If it is, then it is considered a local minimum and is joined with adjacent local minimums using a straight line. The resulting line, across the whole dataset, represents the baseflow hydrograph. The window width is the nearest integer between 3 and 11 that is equal to 2N, where N is empirically determined as follows:
Pettyjohn and Henning Local Minimum Baseflow Separation Method
N - number of days after which runoff ceases;
A - drainage area [km2]*;
* the original equation uses [sq mi], however the tool requires [km2] and converts the drainage area to the proper units internally;

Pettyjohn W.A., Henning R. (1979) - Preliminary estimate of ground-water recharge rates, related streamflow and water quality in Ohio: Ohio State University Water Resources Center Project Completion Report Number 552, p. 323.


3.7. TR-55 Method

TR-55 is a hydrograph separation method developed by USDA starting with 1975 and provides a number of techniques used for small watersheds, especially urbanized watersheds. The presented method predicts the peak rate of runoff as well as the total volume.
TR-55 Baseflow Separation Method
TR-55 Baseflow Separation Method
TR-55 Baseflow Separation Method
q - surface runoff [m3/s];
P - precipitation [mm/d]*;
t - the time (e.g. day) for which the surface runoff is calculated;
CN - runoff curve number;
S - potential maximum retention after runoff begins [in]; calculated based on CN;
Ia - initial abstraction (i.e. sum of all losses before runoff begins) [in]**; calculated based on CN;
* the original equation uses [in/d], however the tool requires [mm/d] and converts the precipitation to the proper units internally;
** if this value is not calculated based on CN but is overridden by the user, the value must be provided in [in]

Cronshey R. et. al (1986) - Urban Hydrology for Small Watersheds, Technical Release 55, June 1986, United States Department of Agriculture.


3.8. Szilagyi Method

Introduced in 1997, this hydrograph separation method estimates the value of the baseflow maximum as well as the baseflow recession hydrograph by fitting the analytical solutions of the Boussinesq equation to the observed discharge values for individual flood events. The method provides both a short-term and a long-term solution, as follows:
Szilagyi Baseflow Separation Method
b - baseflow [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the baseflow is calculated;
d - average interval length between adjacent discharge measurements [h];
a – catchment parameter (see below);
b – catchment parameter (see below);

The short-time solution defines constants a and b as follows:
Szilagyi Baseflow Separation Method
Szilagyi Baseflow Separation Method

For the long-time solution, constants a and b are defined as follows:
Szilagyi Baseflow Separation Method
Szilagyi Baseflow Separation Method

A - catchment area [km2];
D - acquifer depth [m];
L - length of contributing channels [km];
K – saturated hydraulic conductivity [m/s];
Φ - drainable porosity or specific yield [values between 0 and 1];

Szilagyi J., Parlange M. B. (1997) - Baseflow separation based on analytical solutions of the Boussinesq equation, Journal of Hydrology 204 (1998), pp. 251-260.


3.9. Boughton (AWBM 1993) Method

Developed by Walter Boughton in 1993, together with the widespread AWBM rainfall-runoff model. The method is using groundwater levels for hydrograph separation and essentially works as a one-pass filter.
Boughton (AWBM 1993) Baseflow Separation Method
b - baseflow [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the baseflow is calculated;
k – groundwater recession constant [values between 0 and 1];
C - shape constant [positive values];

Boughton W.C. (1993) - A hydrograph-based model for estimating water yield of ungauged catchments. Institute of Engineers Australia National Conference. Publ. 93/14, pp. 317-324.


3.10. Furey and Gupta Method

This hydrograph separation filter, introduced in 2001, is based on a mass balance equation for baseflow through a hillside, and its construction is founded on a physical-statistical theory of low streamflows developed by Furey and Gupta.
Furey and Gupta Baseflow Separation Method
b - baseflow [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the baseflow is calculated;
γ – recession constant [values between 0 and 1];
c1 – ratio of overland flow to precipitation;
c3 – ratio of groundwater recharge to precipitation;
d - time-delay between precipitation and groundwater recharge [d]*;
* d is assumed to be equal to zero, however the user can change this value during method selection;

Gupta V. K., Furey P. R. (2001) - A physically based filter for separating base flow from streamflow time series, Water Resources Research, Vol. 27, No. 11, November 2001, pp. 2709-2722.


3.11. Chapman and Maxwell (1996) Method

The Chapman and Maxwell hydrograph separation digital filter, introduced in 1996, can be viewed as a simplified version of the Boughton filter from 1993.
Chapman and Maxwell Baseflow Separation Method
b - baseflow [m3/s];
Q - streamflow [m3/s];
t - the time (e.g. day) for which the baseflow is calculated;
k - recession constant [values between 0 and 1];

Chapman, T. G., Maxwell, A. I. (1996) - Baseflow separation - comparison of numerical methods with tracer experiments, in Hydrol. and Water Resour. Symp., Institution of Engineers Australia, Hobart. pp. 539-545.


4. User Guide

SepHydro is a web application that performs hydrograph separation (also known as baseflow separation) of streamflow using a multitude of digital filter methods. The tool provides tabular and graphical representations of the input and output data as well as associated representative statistics and exporting functions.
On the left, a contextual menu provides access to the various components of the tool. At the top of the page, a progress bar displays the current status of the analysis.
The three steps required for the use of this tool (1. Load Streamflow Data; 2. Perform Hydrograph Separation; 3. Investigate Results and Export Data) are presented in detail below.

4.1. Load Streamflow Data
The first step in any analysis is to choose a streamflow time series to be used for hydrograph separation. This can be done using the Upload User Data menu entry under the SOURCE DATA menu section. SepHydro allows uploading of files with maximum 7500 rows (~20 years of daily data). It is recommended to split the input data set in blocks of 20 years daily timeseries when the intent is to analyze longer time periods.
Alternatively, the user can use the sample data set provided by using the Load Sample Data menu entry. The data set contains one year of daily streamflow and precipitation records.
A streamflow dataset is a collection of flow records (i.e. discharge) over a period of time. The dataset should contain at least two columns: a time series (e.g. daily, monthly, etc.) and a flow data series. The tool can also accommodate datasets with multiple columns (e.g. with precipitation data series). In this case, the user has to select the columns representing the time series and the flow and precipitation records, respectively.
The tool accepts source data sets in Excel (xlsx) and Comma Separated File (csv) formats. The tool can accommodate missing data or erroneous data entries, by removing the respective lines, but we recommend the user to verify the integrity of the source data before uploading.
Once a streamflow dataset is loaded to the tool, a series of options become available in the SOURCE DATA menu section. These allow for investigation of the streamflow data using a tabular format (Table View) as well as for visualization of the stream hydrograph through a graphical representation (Graphical View). The Graphical View allows for inspection of the full data set or of a reduced subset by changing the beginning and the ending values of the time interval that is displayed. If the dataset contains precipitation, the user can choose to view the precipitation graph by toggling the "Show Precip." box (precipitation is displayed above the streamflow graph). Another toggle, called "Normalize Y", allows the user to redraw the graph with the X axis crossing the Y axis at the minimum Y value, instead of crossing Y axis at Y=0. Statistics for both the entire dataset and the selected subset are shown in both the Table View and the Graphical View.

4.2. Perform Hydrograph Separation
Hydrograph Separation can be performed once the streamflow data is loaded.
The Analyze menu entry becomes available once the streamflow data is loaded. By following this link, the user is presented with a set of methods that can be used for hydrograph separation.
Multiple methods can be used at the same time by clicking the checkboxes next to their names. For each selected method, parameters can be configured using the input boxes.
A preview is also available. The preview is generated only for a limited interval around a local maximum in the data. More in depth visualization is available once the analysis is executed and completed.
The algorithm is forcing baseflow to be equal to streamflow (b = Q) for the instances when the calculated baseflow is larger than the streamflow. The number of such values in the resulting dataset is recorded and can be used for evaluating the parameter selection and overall method reliability.
Once the method selection and parameter configuration is final, the user can proceed to analyze / perform hydrograph separation by clicking the Run Analysis button at the bottom of the page.

4.3. Investigate Results and Export Data
Following completion of the data analysis / hydrograph separation (or baseflow separation), the tool redirects to the Output Data menu section, where the user can investigate the results in tabular (Table View) and graphical form (Graphical View). The OUTPUT DATA menu section becomes visible only when data analysis is complete.
The Graphical View allows the user to choose the desired results for review, as well as a custom viewing interval, by using drop-down lists. If the dataset contains precipitation, the user can choose to view the precipitation graph by toggling the "Show Precip." box (precipitation is displayed above the baseflow graph). Another toggle, called "Normalize Y", allows the user to redraw the graph with the X axis crossing the Y axis at the minimum Y value, instead of crossing Y axis at Y=0. Several statistics, including the baseflow contribution to streamflow are shown on this page. A more comprehensive list of statistics is available under the Table View section.
When multiple methods were used concurrently for baseflow separation, a menu entry named Compare Methods also becomes available. By following this link, the user can configure and view comparison hydrographs between the methods, as well as representative statistics and graphical representations related to the results.
The data set resulted from the analysis (i.e. output dataset), as well as the graphical representations and statistics relevant for the analysis, can be downloaded by following the specific links available in the left-hand menu as well as in the Table View, Graphical View and Compare Methods pages.

5. Limitations

a) SepHydro allows uploading of files with maximum 7500 rows (~20 years of daily data). It is recommended to split the input data set in blocks of 20 years daily timeseries when the intent is to analyze longer time periods;
b) Large data sets (i.e. over 3000 data points) take longer to load. Please be patient until loading is completed. Once the loading is complete the user is redirected to the Source Data Table View page;
c) For large data sets (i.e. over 3000 data points) using all filters simultaneously is not recommended. Users are encouraged to use only up to three hydrograph separation methods at a time;
d) The optimum number of filters that can be used simultaneously can be determined by trial.

6. Terms of Use

SepHydro can be used freely.
The authors do not assume any responsibility for the tool's operation, output, interpretation, or use of results. Please contact the author for additional information and feedback.

7. References

Lyne V., Hollick M. (1979) - Stochastic time-variable rainfall-runoff modelling. Institute of Engineers Australia National Conference. Publ. 79/10, pp. 89-93.
Chapman T. (1991) - Comment on evaluation of automated techniques for base flow and recession analyses, by RJ Nathan and TA McMahon. Water Resources Research, 27(7), pp. 1783-1784.
Eckhardt K. (2005) - How to construct recursive digital filters for baseflow separation. Hydrol. Process. 19, pp. 507–515.
Pettyjohn W.A., Henning R. (1979) - Preliminary estimate of ground-water recharge rates, related streamflow and water quality in Ohio: Ohio State University Water Resources Center Project Completion Report Number 552, p. 323.
Cronshey R. et. al (1986) - Urban Hydrology for Small Watersheds, Technical Release 55, June 1986, United States Department of Agriculture.
Szilagyi J., Parlange M. B. (1997) - Baseflow separation based on analytical solutions of the Boussinesq equation, Journal of Hydrology 204 (1998), pp. 251-260.
Boughton W.C. (1993) - A hydrograph-based model for estimating water yield of ungauged catchments. Institute of Engineers Australia National Conference. Publ. 93/14, pp. 317-324.
Gupta V. K., Furey P. R. (2001) - A physically based filter for separating base flow from streamflow time series, Water Resources Research, Vol. 27, No. 11, November 2001, pp. 2709-2722.
Chapman, T. G., Maxwell, A. I. (1996) - Baseflow separation - comparison of numerical methods with tracer experiments, in Hydrol. and Water Resour. Symp., Institution of Engineers Australia, Hobart. pp. 539-545.

8. Contact

Serban Danielescu, Ph.D.
Research Scientist | Chercheur scientifique
Environment and Climate Change Canada | Environnement et Changements Climatiques Canada
Agriculture and Agri-Food Canada | Agriculture et Agroalimentaire Canada
Fredericton Research and Development Centre | Centre de recherche et développement de Fredericton
95 Innovation Rd., Fredericton, NB, E3B 4Z7
Telephone/Téléphone: 506-460-4468
Facsimile/Télécopieur: 506-460-4377