Thresholds

Temperature thresholds were calculated from the historical Coastal Monitoring Program data. Preliminary quality control was applied to this data, and obvious outliers and data from freshwater stations were omitted.

Separate thresholds were calculated for each county due to expected and observed spatial differences in temperature (Figure 1; Figure 2).

In the future, it may be useful to resolve temperature thresholds by depth as well as county (or smaller spatial scales). However, the data become very patchy when grouped by depth, which could result in less useful thresholds values. Additionally, the same depth can represent a different part of the water column for different stations. For example, at the Barren Island station in Guysborough County, the 15 m sensor is near the bottom. In contrast, 15 m is in the top 20 % of the water column at Tickle Island, another station in Guysborough County.

Gross Range Test

Sensor Thresholds

The sensor thresholds were determined based on the associated manual (Table 1).

Table 1: Temperature sensor thresholds for the Gross Range Test.

User Thresholds

The mean and standard deviation were used to determine \(user_{min}\) and \(user_{max}\). This provided useful thresholds, although different statistics may be more representative based on the distribution of the observations for some counties (Figure 1). The statistics and threshold values are shown in Table 2.

The \(user_{min}\) threshold is typically << 0 °C (Table 2), and is therefore not expected to flag any observations. For most counties (all except Annapolis, Queens, Shelburne, and Digby), the \(user_{min}\) is less than the \(sensor_{min}\) for the aquameasure and vr2ar sensors. In this case, any observations less than the \(sensor_{min}\) would fail the Gross Range Test (i.e., the \(user_{min}\) would be ignored).

It may be useful for other users to apply their own \(user_{min}\) threshold to highlight cold observations that are Suspect/Of Interest. For example, those interested in salmonid aquaculture may wish to flag observations at or near the superchill threshold (-0.7 °C).

Figure 1: Distribution of temperature observations (binwidth = 1 °C). Dotted orange lines indicate the user thresholds.

Table 2: Gross Range Test statistics and user thresholds for temperature.

Climatological Test

Figure 2 shows the monthly mean and standard deviation of temperature data for each county. Observations from all depths were pooled together. This resulted in high standard deviations for some counties, particularly in the summer months when temperature stratification is the strongest (?@fig-temp-climatology; Table 3).

The observations are normally distributed within each month (?@fig-do-hist-clim), and so the mean and standard deviation were used to calculate the seasonal thresholds (Table 3).

Figure 2: Mean +/- 3 standard deviations of the monthly depth observations.

Figure 3

Seasonal distribution of depth observations (binwidth = 2 °C). Dotted orange lines indicate the user thresholds.

Table 3: Seasonal statistics and thresholds for temperature.

Spike Test

For all counties, the distribution of the spike value is skewed right, and so several upper quartile values (90th, 95th, and 99.7th quartile) were evaluated to use as the \(spike_{low}\). There were relatively few large single-value spikes in the temeprature data, and so the 99.7th quartile was selected to avoid false positives. \(spike_{high}\) was defined as 3 * \(spike_{low}\) to identify especially egregious spike values.

Figure 4 shows the distribution of the spike value of temperature observations (binwidth = 0.1 °C). There is a separate tab for each county. Dotted orange line indicates \(spike_{low}\); dotted red line indicates \(spike_{high}\). Note difference in y-axis scale.

Figure 4

Table 4: Spike thresholds for temperature.

Rolling Standard Deviation Test

The distribution of rolling standard deviation is skewed right (?@fig-sal-roll-sd-dist), and so several upper quartile values were evaluated to use as the \(rolling\_sd\_max\).

The 90th, 95th, and 99.7th quartile values were each applied to the the raw data (no preliminary QC) and the results inspected. There were minimal suspect changes in rolling standard deviation for temperature, and so the 99.7th quartile was selected to avoid false positives.

Table 5: Rolling standard deviation threshold for temperature.