[Lower Freshwater Creek]

Freshwater Creek, California...

Annual and Storm Analysis, Turbidity Threshold Sampling (TTS)


The equipment for this monitoring station was purchased by Salmon Forever. The station was cooperatively installed by the Forest Service, Redwood Sciences Laboratory (RSL), Salmon Forever, Humboldt State University Fisheries Professor Terry Roelofs, and members of the Freshwater community. Because the monitoring station was activated on January 13, 1999, a significant portion of hydrologic year 1999 was not included in the data and, therefore, not included in the analysis presented below. For example, the first major storm event of the hydrologic year is not included in the annual load estimates presented.  Although this storm, which peaked on November 21, had the highest peak of the 1999 hydrologic year, no flow or sediment data were collected and no attempt was made to quantify the amount of sediment delivered.

The purpose of this analysis was to estimate sediment loads using the relations between suspended sediment concentration, turbidity and discharge in Freshwater Creek.  The analysis as presented below was carried out at RSL on a Unix operating system using S-Plus and Perl.  The majority of the analysis tools and plotting routines were written in S-Plus because of its powerful ability to manipulate and graphically represent the data.

Instrumentation included a Campbell data logger, a pressure transducer for recording stage height (water depth), an ISCO sampler to extract pumped sediment samples from the stream, and an infrared turbidity probe.  The pressure transducer was calibrated in the field, and the electronic measurements were adjusted to observations made systematically in the field when appropriate.  Stage height and turbidity were recorded at 15-minute intervals.  Manual discharge measurements were made in the field to develop a relationship between stage height and discharge.  Pumped samples were collected automatically at pre-determined turbidity thresholds according to a program loaded in the data logger.  These sample bottles were analyzed for sediment concentration at RSL and the Salmon Forever Sunnybrae Sediment Lab using methods described in Standard Methods for the Examination of Water and Wastewater.  The turbidity measurements recorded by the Campbell data logger were converted to suspended sediment concentrations using regressions based on the pumped sample concentrations.  The regression methodology is described in following sections of the web page.

Depth-integrated water/sediment samples were also collected at various stages so that the point samples taken by the ISCO sampler could be adjusted to account for sediment delivered across the entire cross section.  However due to an insufficient number of depth-integrated samples, this portion of the analysis was not carried out.  Thus all load estimates presented are based on the point samples taken by the ISCO sampler.  If an adjustment based on an adequate number of samples were to be carried out, it is likely that it would result in higher sediment loads than those presented below, because sediment concentrations are expected to be greater near the stream bed.


Development of the Stage - Discharge Rating Curve at the Freshwater Station



A prerequisite for estimation of the sediment load is development of a discharge rating curve for converting stage height to water discharge. This was carried out using the standard method of discharge measurement using either a Price AA or pygmy current meter to measure water velocities at the stream cross-section (refer to Stream and Channel Reference Sites: An illustrated guide to field technique by C.C. Harrelson, C.L. Rawlins, and J.P. Potyondy, USDA Forest Service General Technical Report RM-245, April 1994).  Discharge was computed by summing up the products of velocity and partial cross-sectional area over the entire cross-section. To obtain a discharge rating curve, a power relationship, as shown in the figure below, was fit to the data presented in the following table.  The power relationship assumes a cross-section that is non-deforming in time. However, the cross section at the Freshwater station was not located on bedrock and is directly downstream from a meander bend in the channel. These characteristics of the sampling site are clearly undesirable and add substantial uncertainty to the load estimates given in the following sections of this web page.  More discharge measurements over time would quantify the uncertainty involved with this aspect of the analysis.
 

Stage - Discharge Rating Curve Data

DATE TIME QUALITY STAGE (ft) DISCHARGE (cfs)
12/17/98 1545 good 0.60 31.37
01/16/99 1635 good 1.05 67.09
01/19/99 1300 fair 1.60 160.18
01/23/99 1116 good 2.37 328.38
01/24/99 1000 good 1.41 109.80
02/06/99 921 poor 1.11 74.25
02/06/99 1009 poor 1.65 215.46
02/06/99 1054 poor 2.61 413.01
02/06/99 1150 poor 3.65 674.87
02/06/99 1306 fair 4.00 831.00
02/06/99 1524 good 3.30 423.70
02/06/99 1602 good 3.10 363.17
02/06/99 1701 fair 2.75 325.92
02/06/99 1758 fair 2.58 285.34
02/08/99 1556 good 1.42 116.50
02/15/99 905 good 1.01 68.60

Stage/Discharge Rating Curve


Annual Regression Discussion



This first plot shows the complete data set used for the annual load estimates. The numbers correspond to bottles analyzed in a given dump (batch of up to 24 bottles). For example, all number 14's correspond to bottles analyzed in dump 14. A total of 18 dumps, including approximately 154 bottles, were used for the regression shown in plot 1.

Plot 1

A linear equation seemed to do a poor job in fitting the lower portion of the data set. A quadratic, a cubic polynomial, and a smooth loess were also fit to the data set. "Loess", short for LOcal regrESSion, is flexible tool for fitting almost any shape curve. There is a good description of loess in William S. Cleveland's book "Visualizing Data" (Hobart Press, 1993). The loess fit was determined to give the best fit to the data set. The results of the linear and the loess fit are given below including the annual load predictions. Notice that the linear fit gives a fairly good R-SQUARED value and that the linear and loess predictions are very close.
 

DESCRIPTOR LINEAR LOESS
NUMBER OF OBSERVATIONS 154 154
R-SQUARED 0.971 0.980
RESIDUAL STANDARD ERROR 65.3 57.17
PREDICTED LOAD (kg) 2845365 2800470
PREDICTED LOAD (kg/ha) 826 813
PREDICTED LOAD (ton/mi2) 236 232
Note: Freshwater Watershed Area Above Study Site: 13.3 mi2 (3445 ha)

Plot 2 (shown below) summarizes the entire data set for hydrologic year 1999 at Freshwater (January 13, 1999 to August 1, 1999). The upper plot is the sedigraph and the lower plot is the hydrograph. The sedigraph is based on the linear regression model presented above where suspended sediment concentration is assumed to be a linear function of turbidity. The suspended sediment concentration is shown on the upper left axis (in mg/L) and the turbidity scale is on the right axis (in NTU). For the entire period, the load estimate (shown below the legend) matches the prediction made by the linear model presented above. The sedigraph peaks tended to occur before the hydrograph peaks.

Plot 2


Individual Storm Analysis


How well does the loess regression predict the load delivered by individual storms? Eight storms were selected for analysis (these storms are numbered 1-8 on plot 2 above). Each storm period was first designated. Then a linear regression was carried out on the samples analyzed under the dump(s) that corresponded to the designated storms. In some cases a rating curve relating discharge and suspended sediment was used where the linear regression between suspended sediment and turbidity predicted negative values. This was often at the very beginning or end of the storm. As an example, we shall select Storm 3 (start:02/06/99, 0300 hrs; end:02/08/99, 1315 hrs). The initial regression using all ISCO sample bottles taken during the storm was carried out and is shown below in plot 3.

Plot 3

However, this linear regression model predicted negative values at the very beginning and end of the storm. A power relationship between discharge and suspended sediment was then determined using the samples taken for Storm 3. This plot is shown below as plot 4. The plot shows a characteristic hysteresis loop in the suspended sediment - discharge rating curves. This can be seen by tracing the bottle numbers in plot 4 for either dump.

Plot 4

Notice the hysteresis loop is formed because a greater suspended sediment concentration was observed on the rising limb of the hydrograph, while the suspended sediment concentration was lower for a given discharge on the falling limb of the hydrograph.  Because the objective was to correct the negative predictions at the front end of the storm, another rating curve was developed using only the first three sample bottles taken. This rating curve is shown below in plot 5.

Plot 5

Finally, a regression using suspended sediment concentration and turbidity values from the last three sample bottles were used to remove the negative predictions at the end of the storm given by the original regression (refer to plot 3 above). This final regression is shown below as plot 6.

Plot 6

All three of these relationships are then pieced together to form the final estimate for the storm. Plot 7 (below) shows the resulting sedigraph and hydrograph for Storm 3. Notice the load estimate for the storm is given under the legend.

Plot 7

The table below shows each storm estimate compared to an estimate computed from the annual loess fit. Notice that all eight storms accounted for 86% of the total annual load estimate. As shown in the table below, the loess fit gave fairly good estimates for individual storms with the largest deviation observed to be 28%. The deviations are explained by observing the distribution of points in plot 1. By observing how the dumps (i.e. ISCO sample bottles) corresponding to a given storm fall in relation to the loess curve, one can predict whether the loess model will under- or over-estimate the sediment load delivered by the storm event. The bias in using an annual relationship to estimate individual storms was limited in this example because there was relatively little scatter in the annual data set.
 

STORM LOAD ESTIMATE (kg) 
STORM REGRESSION
LOAD ESTIMATE (kg) 
ANNUAL LOESS
% DIFFERENCE
1 273493 303014 10.8%
2 309508 277954 -10.2%
3 892185 852181 -4.5%
4 212461 179649 -15.4%
5 56836 41171 -27.6%
6 614202 561425 -8.6%
7 23069 20070 -13.0%
8 60676 56878 -6.3%
TOTAL 2442430 2292342 -6.2%
LOESS LOAD ESTIMATE OF ANNUAL LOAD: 2845365kg
TOTAL ESTIMATED STORM LOAD AS A PERCENTAGE OF ESTIMATED ANNUAL LOAD: 85.8%

[Coho Banner]
| Turbidity Threshold Publications | Back to Main Page |
- Salmon Forever - Redwood Community Action Agency - Redwood Sciences Lab[Coho Banner]

Freshwater // Last updated on October 19, 2000, by Ben Bray