In the SSP_Examples.ssp study, the Bulletin 17 analysis Bushkill_ExpandingWindow illustrates how to utilize the Expanding Window computational option.  In this example, an initial period of 20 years is incrementally expanded by 5 years to explore how individual flood peaks impact a flow-frequency curve as well as the parameters of the Log Pearson Type III Distribution.

Input Data

The Bushkill Creek at Shoemakers, PA stream gage has an annual peak record consisting of 111 peaks beginning in 1909 and ending in 2019, as shown in Figure 1. The annual peak flow record is tabulated within Table 1.

Figure 1. Bushkill Creek Annual Peak Flow Record

Table 1. Bushkill Creek at Shoemakers, PA Annual Peak Flow Record

DateFlow (cfs)
20 Feb 19092270
07 Mar 19101190
14 Jun 19111120
15 Mar 19121190
28 Mar 19131350
29 Mar 19141350
24 Feb 19152070
02 Apr 19161430
27 Mar 19171510
26 Feb 19181970
22 Jul 19192380
24 Jul 19205250
10 Mar 19211780
08 Mar 19221780
24 Mar 19231270
30 Sep 19242830
12 Feb 19251430
16 Nov 19252380
16 Nov 19262830
19 Oct 19271970
15 Mar 19291190
18 Nov 19291150
29 Mar 19311080
01 Apr 19321430
16 Sep 19333070
12 Apr 19341040
10 Jul 19353330
18 Mar 19364770
06 Apr 19371240
25 Jan 19382070
06 Dec 19382540
31 Mar 19402790
06 Apr 19411220
28 Sep 19422920
30 Dec 19421710
25 Apr 19441500
19 Jul 19453900
28 May 19462080
06 Apr 19471580
21 Mar 19482080
31 Dec 19481780
29 Mar 19501540
31 Mar 19512860
08 Nov 19512980
11 Dec 19524680
07 Dec 1953928
19 Aug 195523400
16 Oct 19553680
06 Apr 19571460
21 Dec 19573050
06 Mar 19591260
05 Apr 19601800
26 Feb 19611890
08 Apr 19621060
27 Mar 19632680
10 Mar 19641800
09 Feb 1965804
06 Mar 1966880
15 Mar 19671200
30 May 19682290
28 Jul 19697300
02 Apr 19703130
23 Oct 19702070
23 Jun 19723360
03 Feb 19731730
21 Dec 19734280
08 Dec 19742910
28 Jan 19762830
25 Feb 19772120
09 Jan 19784720
25 Jan 19792830
22 Mar 19802670
12 May 19812770
04 Feb 19821300
16 Apr 19833160
05 Apr 19843340
27 Sep 19851420
15 Mar 19862730
13 Sep 19872580
30 Nov 19871260
06 May 19892120
20 Oct 19893260
04 Dec 19902470
31 May 19921780
01 Apr 19933210
28 Mar 19941750
08 Mar 19951280
27 Jan 19964990
02 Dec 19963060
12 May 19981450
24 Jan 19991360
28 Feb 2000895
17 Dec 20001930
29 May 20021210
22 Jun 20032340
18 Sep 20045330
03 Apr 20055670
28 Jun 20065380
16 Apr 20072790
09 Mar 20082750
12 Dec 20082590
25 Jan 20101940
28 Aug 20115480
08 Dec 20111480
11 Jun 20131740
01 May 20142030
01 Jul 20152790
25 Feb 20162220
07 Apr 20171970
25 Feb 20181550
24 Jan 20192070

A Bulletin 17 Analysis using Bulletin 17C procedures and the Variable Time Window | Expanding Window compute option has been developed for this example. To open the analysis, either double-click on the analysis labeled Bushkill_ExpandingWindow from the Study Explorer or from the Analysis menu select open, then select Bushkill_ExpandingWindow from the list of available analyses. When Bushkill_ExpandingWindow is selected, the Bulletin 17 analysis editor General tab will appear as shown in Figure 2.  As shown, the 17C EMA computational option was selected along with all the subsequent default options.
 
Figure 2. Bulletin 17 Analysis General Tab for Bushkill_ExpandingWindow

On the Options tab, the Compute Using Expanding Window option was selected within the Variable Time Window panel.  An Initial Window Length of 20 years and an Expansion Increment of 5 years was entered, as shown in Figure 3.

 Figure 3. Bulletin 17 Analysis Options Tab for Bushkill_ExpandingWindow

The EMA Data tab for this example is shown in Figure 4. Since this example uses an annual maximum series consisting entirely of systematic data with a complete record, a single zero – inf perception threshold is adequate. No modifications to the default flow ranges and data types are necessary.
Figure 4. Bulletin 17 Analysis EMA Data Tab for Bushkill_ExpandingWindow

Computing the Analysis

Once all of the General, Options, and EMA Data tab settings are set or selected, the user can press the Compute button to perform the analysis. A Compute Warnings message will be shown noting that the final expansion increment was less than 5 years.  Once the computations have been completed, a message window will open stating Compute Complete.

Close these windows and then select the Tabular Results tab. The analysis window should resemble Figure 5.

Figure 5. Bulletin 17 Analysis Tabular Results Tab for Bushkill_ExpandingWindow

When using the Expanding Time Window compute option, the Tabular Results tab will only contain the results from the final time window computation.  In this example, the final time window spans water years 1909 - 2019, which is 111 total years.

In addition to the tabular results, a graphical plot of the computed frequency curves for the final time window can be obtained by pressing the Plot Curve button at the bottom of the analysis window.

Computed frequency curve

Expanding Time Window Results

Upon clicking the Variable Time Window Results tab, two sub-tabs become available: Flow and Statistics.  The Flow sub-tab contains the evolution of quantile information for each time window and each of the frequency ordinates specified within the Options tab | Output Frequency Ordinates table, as shown within Figure 6.

Figure 7. Quantile information for each time window in the expanding window analysis

Within this example, 19 time windows were derived from the input data, as shown in Table 2.  

Table 2. Time Windows for Bushkill_ExpandingWindow Bulletin 17 Analysis

Begin YearEnd YearNumber of Years
1909193325
1909193830
1909194335
1909194840
1909195345
1909195850
1909196355
1909196860
1909197365
1909197870
1909198375
1909198880
1909199385
1909199890
1909200395
19092008100
19092013105
19092018110
19092019111

The Statistics sub-tab contains the evolution of Log Pearson Type III parameters for each time window, as shown within Figure 7.

Figure 7. Log Person III parameters for each 20 year window in the expanding window analysis

Notice the extreme event in 1955. The peak flow rate for the 1955 event was 23,400 cfs; whereas, the next largest flood was 7,300 cfs in 1969.  The inclusion (or exclusion) of this flood event causes dramatic changes to the computed quantiles (e.g. 1/100 annual exceedance probability flow rate) as well LPIII parameters (e.g. skew).

Report File

In addition to the tabular and graphical results, there is a report file that shows the order in which the calculations were performed. To review the report file, press the View Report button at the bottom of the analysis window. When this button is selected, a text viewer will open the file and display it on the screen. Shown in Figure 8 is the report file. The report file contains a listing of the input data, preliminary results, outlier and historical data tests, additional calculations needed, and the final frequency curve results for each of the time windows.  Additionally, a final summarization of all the time windows is included. Different types and amounts of information will show up in the report file depending on the data and the options that have been selected for the analysis. The user should review the report file to understand how the Bulletin 17C procedures were applied for each time window.

Figure 8. Report File for Bushkill_ExpandingWindow