Search

BackgroundData from the COVID-19 pandemic describes increases in drug use and related harms, especially fatal overdose. However, evidence is needed to better understand the pathways from pandemic-related factors to substance use behaviours. Thus, we investigated stockpiling drugs among people who use drugs (PWUD) in five cities in the United States and Canada.MethodsWe used data from two waves of interviews among participants in nine prospective cohorts to estimate the prevalence and correlates of stockpiling drugs in the previous month. Longitudinal correlates were identified using bivariate and multivariate generalized linear mixed-effects modeling analyses.ResultsFrom May 2020 to February 2021, we recruited 1873 individuals who completed 2242 interviews, of whom 217 (11.6%) reported stockpiling drugs in the last month at baseline. In the multivariate model, stockpiling drugs was significantly and positively associated with reporting being greatly impacted by COVID-19 (Adjusted Odds Ratio [AOR]= 1.21, 95% CI: 1.09-1.45), and at least daily use of methamphetamine (AOR = 4.67, 95% CI: 2.75-7.94) in the past month.ConclusionsWe observed that approximately one-in-ten participants reported stocking up on drugs during the COVID-19 pandemic. This behaviour was associated with important drug-related risk factors including high-intensity methamphetamine use. While these correlations need further inquiry, it is possible that addressing the impact of COVID-19 on vulnerable PWUD could help limit drug stockpiling, which may lower rates of high-intensity stimulant use.

Supraglacial debris is significant in many regions and complicates modeling of glacier melt, which is required for predicting glacier change and its influences on hydrology and sea-level rise. Temperature-index models are a popular alternative to energy-balance models when forcing data are limited, but their transferability among glaciers and inherent uncertainty have not been documented in application to debris-covered glaciers. Here, melt factors were compiled directly from published studies or computed from reported melt and MERRA-2 air temperature for 27 debris-covered glaciers around the world. Linear mixed-effects models were fit to predict melt factors from debris thickness and variables including debris lithology and MERRA-2 radiative exchange. The models were tested by leave-one-site-out cross-validation based on predicted melt rates. The best model included debris thickness (fixed effect) and glacier and year (random effects). Predictions were more accurate using MERRA-2 than on-site air temperature data, and pooling MERRA-2-derived and reported melt factors improved cross-validation accuracy more than including additional predictors such as shortwave or longwave radiation. At one glacier where monthly ablation was measured over 4 years, seasonal variation of melt factors suggested that heat storage significantly affected the relation between melt and energy exchange at the debris surface.

ObjectivesTo define the incidence of clinically-detected COVID-19 in people with HIV (PWH) in the US and evaluate how racial and ethnic disparities, comorbidities, and HIV-related factors contribute to risk of COVID-19.DesignObservational study within the CFAR Network of Integrated Clinical Systems cohort in 7 cities during 2020.MethodsWe calculated cumulative incidence rates of COVID-19 diagnosis among PWH in routine care by key characteristics including race/ethnicity, current and lowest CD4 count, and geographic area. We evaluated risk factors for COVID-19 among PWH using relative risk regression models adjusted with disease risk scores.ResultsAmong 16,056 PWH in care, of whom 44.5% were Black, 12.5% were Hispanic, with a median age of 52 years (IQR 40-59), 18% had a current CD4 count < 350, including 7% < 200; 95.5% were on antiretroviral therapy, and 85.6% were virologically suppressed. Overall in 2020, 649 PWH were diagnosed with COVID-19 for a rate of 4.94 cases per 100 person-years. The cumulative incidence of COVID-19 was 2.4-fold and 1.7-fold higher in Hispanic and Black PWH respectively, than non-Hispanic White PWH. In adjusted analyses, factors associated with COVID-19 included female sex, Hispanic or Black identity, lowest historical CD4 count ConclusionsOur results suggest that the presence of structural racial inequities above and beyond medical comorbidities increased the risk of COVID-19 among PWHPWH with immune exhaustion as evidenced by lowest historical CD4 or current low CD4:CD8 ratio had greater risk of COVID-19.

The North American Regional Reanalysis (NARR) precipitation product was evaluated using station observations and catchment water yield in British Columbia (BC), Canada, at inter-annual, monthly, an...

Stream temperature is a critical parameter for understanding hydrological and biological processes in stream ecosystems. Although a large body of research has addressed the effects of forest harvesting on stream temperature, less is known about the responses of stream temperature to the practice of strip-thinning, which produces more coherent patches of shade and sunlight areas. In this study, we examined stream temperature response to 50% strip-thinning in a 17 ha headwater catchment. The thinning lines extended through the riparian zone. Paired-catchment analysis was applied to estimate changes in daily maximum, mean, and minimum stream temperatures for the first year following treatment. Significant effects on daily maximum stream temperature were found for April to August, ranging from 0.6 °C to 3.9 °C, similar to the magnitude of effect found in previous studies involving 50% random thinning. We conducted further analysis to identify the thermal response variability in relation to hydrometeorological drivers. Multiple regression analysis revealed that treatment effects for maximum daily stream temperature were positively related to solar radiation and negatively related to discharge. Frequent precipitation during the summer monsoon season produced moderate increases in discharge (from 1 to 5 mm day−1), mitigating stream temperature increases associated with solar radiation. Catchment hydrologic response to rain events can play an important role in controlling stream thermal response to forest management practices.

Understanding seasonal snow cover dynamics is critical for management of hydrological regimes, habitat availability for wildlife species, forest fire risk assessment and recreational demands. Although data products provided at 500 m spatial resolution by the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor provide important and readily available information on snow cover, capturing snow dynamics at finer spatial resolutions remains problematic due to the lack of high temporal and spatial resolution data, which limits the number of available observations each year. In this paper we present a new approach to create a daily time-series of 30-m snow observations (called SNOWARP), derived from daily MODIS Normalised Difference Snow Index (NDSI) snow cover data to capture the temporal dynamics of snow cover and Dynamic Time Warping (DTW) to re-order historical Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and Operational Land Imager (OLI) observations to account for inter-annual variability. The SNOWARP product was produced for 2000–2018 for an area of Western Alberta (approximately 30,000 km2) and was calibrated against a network of time-lapse cameras and snow pillows. Results indicate the RMSE of the SNOWARP fractional product ranges from 31.3%–68.3%, while F score of the SNOWARP binary product ranges from 87.7% - 98.6% when compared to ground truth data. Capturing fractional snow cover at a fine spatial and temporal scale is important due to the spatial heterogeneity of snow cover, particularly in mountainous regions with implications for biodiversity assessment and monitoring. SNOWARP demonstrates a novel method to increasing the temporal resolution of Landsat-derived snow cover data, providing valuable insights on regional snow cover dynamics for use in a range of applications.

Stream surface albedo plays a key role in the energy balance of rivers and streams that are exposed to direct solar radiation. Most physically based analyses and models have incorporated a constant stream albedo between 0.03 and 0.10, based primarily on measurements from low-gradient streams with low suspended sediment concentrations. However, albedo should vary with solar elevation angle, suspended sediment concentration, aeration, and fraction of direct vs diffuse radiation. The objective of this study was to quantify the dependence of albedo of mountain streams on the controlling factors and to develop a predictive model for use in physically based analysis and modelling of stream temperature, especially for future climate and land use scenarios. Stream surface albedo was measured at nine sites with a variety of gradients and suspended sediment characteristics in the southern Coast Mountains of British Columbia, Canada. As expected, albedo of low-gradient, non-whitewater (flatwater) streams increased with solar zenith angle, suspended sediment concentration, and proportion of diffuse to direct solar radiation, ranging between 0.025 during cloudy periods over clear water to 0.25 for turbid water at zenith angles of less than 20 degrees. Albedo was enhanced in steep reaches or at channel steps and cascades where flow was visibly aerated, with a range of 0.09 to 0.33. In clear weather, albedo exhibited notable diurnal variability at flatwater sampling sites. For example, during late summer, surface albedotypically fluctuated between 0.08 and 0.15 on a daily basis at a flatwater site on the highly turbid, glacier-fed Lillooet River. Multiple regression models explained approximately 60% and 40% of the variance under cross validation for flatwater and whitewater data subsets, respectively, with corresponding root-mean-square errors of approximately 0.02 and 0.06.

We describe a model for predicting the moisture content of standardized fuel sticks. Compared to the previously published “Nelson model,” which is the basis for the model currently used by agencies involved with fire management, our model's treatment of moisture diffusion, heat conduction, precipitation interception, and evaporation/sorption processes is simplified while its treatment of radiation transfer is more sophisticated. Two independent datasets were used to calibrate and evaluate the model. The optimized models were able to predict 78–94% of the variance in the calibration observations with absolute mean biases ranging from 0.04% to 0.38% fuel moisture content. When evaluated using both an independent time period and an independent dataset, the explained variance decreased slightly, ranging from 72 to 94%, and the model bias increased with absolute mean biases ranging from 0.29% to 2.48%. The best results were achieved when the model was forced with observations taken at the same height and location as the fuel moisture sticks. Model skill was reduced for higher moisture levels. Despite the simplified approach, our model improved on the skill achieved by the Nelson model and has additional features that allow for a more realistic treatment of canopy coverage and changes in sky conditions.

Tracer dilution methods are commonly used to measure discharge in steep mountain streams. This research addressed knowledge gaps associated with dilution methods using original field data collected on nine streams in southwest British Columbia and discharge measurements conducted by Northwest Hydraulic Consultants Ltd. Minimum mixing lengths ranged between 2.4 and 24.5 stream wetted widths, but determining the mixing length can be confounded by surface-subsurface water fluxes. Probes need to be placed on opposite sides of the stream to verify adequate mixing, because probes located at different locations on the same of the stream sometimes suggested complete mixing had occurred when it in fact had not. For probes located downstream of complete mixing, breakthrough curves (BTCs) for probes located in the main current differed significantly from probes in zones with recirculating flow, even though they yielded discharge values within ± 10%. The peak of the BTC is a function of the mass of tracer injected, r...

There has been increasing attention over the last decade to the potential effects of glacier retreat on downstream discharge and aquatic habitat. This study focused on streamflow variability downstream of Bridge Glacier in the southern Coast Mountains of BC between 1979 and 2014, prior to and during a period in which the glacier experienced enhanced calving and rapid retreat across a lake-filled basin. Here we combined empirical trend detection and a conceptual-parametric hydrological model to address the following hypotheses: (1) streamflow trends in late summer and early autumn should reflect the opposing influences of climatic warming (which would tend to increase unit-area meltwater production) and the reduction in glacier area (which would tend to reduce the total volume of meltwater generated), and (2) winter streamflow should increase because of displacement of lake water as ice flows past the grounding line and calves into the lake basin. In relation to the first hypothesis, we found no significant trends in monthly discharge during summer. However, applying regression analysis to account for air temperature and precipitation variations, weak but statistically significant negative trends were detected for August and melt season discharge. The HBV-EC model was applied using time-varying glacier cover, as derived from Landsat imagery. Relative to simulations based on constant glacier extent, model results indicated that glacier recession caused a decline in mean monthly streamflow of 9% in August and 11% in September. These declines in late-summer streamflow are consistent with the results from our empirical analysis. The second hypothesis is supported by the finding of positive trends for December, January, and February discharge. Despite the modelled declines in late-summer mean monthly streamflow, recorded discharge data exhibited neither positive nor negative trends during the melt season, suggesting that Bridge Glacier may currently be at or close to the point of peak water. Further analysis of the impact of lake-terminating glaciers on downstream discharge is needed to refine the peak water model. Copyright © 2016 John Wiley & Sons, Ltd.

Meteorological observations at high elevations in mountainous regions are often lacking. One opportunity to fill this data gap is through the use of downscaled output from weather reanalysis models. In this study, we tested the accuracy of downscaled output from the North American Regional Reanalysis (NARR) against high-elevation surface observations at four ridgetop locations in the southern Coast Mountains of British Columbia, Canada. NARR model output was downscaled to the surface observation locations through three-dimensional interpolation for air temperature, vapour pressure and wind speed and two-dimensional interpolation for radiation variables. Accuracy was tested at both the 3-hourly and daily time scales. Air temperature displayed a high level of agreement, especially at the daily scale, with root mean square error (RMSE) values ranging from 0.98 to 1.21 °C across all sites. Vapour pressure downscaling accuracy was also quite high (RMSE of 0.06 to 0.11 hPa) but displayed some site specific bias. Although NARR overestimated wind speed, there were moderate to strong linear relations (r2 from 0.38 to 0.84 for daily means), suggesting that the NARR output could be used as an index and bias-corrected. NARR output reproduced the seasonal cycle for incoming short-wave radiation, with Nash–Sutcliffe model efficiencies ranging from 0.78 to 0.87, but accuracy suffered on days with cloud cover, resulting in a positive bias and RMSE ranged from 42 to 46 Wm− 2. Although fewer data were available, incoming long-wave radiation from NARR had an RMSE of 19 Wm− 2 and outperformed common methods for estimating incoming long-wave radiation. NARR air temperature showed potential to assist in hydrologic analysis and modelling during an atmospheric river storm event, which are characterized by warm and wet air masses with atypical vertical temperature gradients. The incorporation of a synthetic NARR air temperature station to better represent the higher freezing levels resulted in increased predicted peak flows, which better match the observed run-off during the event. Copyright © 2016 John Wiley & Sons, Ltd.

Paired-catchment studies conducted on small (< 10 km2) rain-dominated catchments revealed that forest harvesting resulted in a period of increased warm-season low flows ranging from less than five years to more than two decades, consistent with the results of stand-level studies and process considerations. Of the five paired-catchment studies in snow-dominated regions, none revealed a statistically significant change in warm-season low flows in the first decade following harvest, although two exhibited non-significant higher flows in August and September and one had lower flows. Two studies, one of rain-dominated catchments and one of snow-dominated catchments, found that summer low flows became more severe (i.e., lower) about two decades or so following harvest. These longer-term results indicate that indices such as equivalent clearcut area, as currently calculated using monotonic recovery curves, may not accurately reflect the nature of post-harvest changes in low flows. Studies focussed on medium to large catchments (tens to thousands of km2 in area) found either no statistically significant relations between warm-season low flows and forest disturbance, or inconsistent responses. Attempts to synthesize existing studies are hampered by the lack of a common low-flow metric among studies, as well as detailed information on post-harvest vegetation changes. Further fieldresearch and process-based modelling is required to help elucidate the underlying processes leading to the results from these paired-catchment studies and to enhance the ability to predict streamflow responses to forest harvesting, especially in the context of a changing climate. KEYWORDS: streamflow; forestry; low flows; fish habitat; hydrologic recovery

We describe a method to calculate regional snow line elevations and annual equilibrium line altitudes (ELAs) from daily MODIS imagery (MOD02QKM) on large glaciers and icefields in western North America. An automated cluster analysis of the cloud-masked visible and near-infrared bands at 250 m resolution is used to delineate glacier facies (snow and ice) for ten glacierized regions between 2000–2011. For each region and season, the maximum observed value of the 20th percentile of snow-covered pixels (ZS(20)) is used to define a regional ELA proxy (ELAest). Our results indicate significant increases in the regional ELA proxy at two continental sites (Peyto Glacier and Gulkana Glacier) over the period of observation, though no statistically significant trends are identified at other sites. To evaluate the utility of regional ELA proxies derived from MOD02QKM imagery, we compare standard geodetic estimates of glacier mass change with estimates derived from historical mass balance gradients and observations of ZS(20) at three large icefields. Our approach yields estimates of mass change that more negative than traditional geodetic approaches, though MODIS-derived estimates are within the margins of error at all three sites. Both estimates of glacier mass change corroborate the continued mass loss of glaciers in western North America. Between 2000 and 2009, the geodetic change approach yields mean annual rates of surface elevation change for the Columbia, Lillooet, and Sittakanay icefields of −0.29 ± 0.05, −0.26 ± 0.05, and −0.63 ± 0.17 m a−1, respectively. This study provides a new technique for glacier facies detection at daily timescales, and contributes to the development of regional estimates of glacier mass change, both of which are critical for studies of glacier contributions to streamflow and global sea level rise.

Glacier melt provides important contributions to streamflow in many mountainous regions. Hydrologic model calibration in glacier-fed catchments is difficult because errors in modelling snow accumulation can be offset by compensating errors in glacier melt. This problem is particularly severe in catchments with modest glacier cover, where goodness-of-fit statistics such as the Nash-Sutcliffe model efficiency may not be highly sensitive to the streamflow variance associated with glacier melt. While glacier mass balance measurements can be used to aid model calibration, they are absent for most catchments. We introduce the use of glacier volume change determined from repeated glacier mapping in a guided GLUE (generalized likelihood uncertainty estimation) procedure to calibrate a hydrologic model. This approach is applied to the Mica basin in the Canadian portion of the Columbia River Basin using the HBV-EC hydrologic model. Use of glacier volume change in the calibration procedure effectively reduced parameter uncertainty and helped to ensure that the model was accurately predicting glacier mass balance as well as streamflow. The seasonal and interannual variations in glacier melt contributions were assessed by running the calibrated model with historic glacier cover and also after converting all glacierized areas to alpine land cover in the model setup. Sensitivity of modelled streamflow to historic changes in glacier cover and to projected glacier changes for a climate warming scenario was assessed by comparing simulations using static glacier cover to simulations that accommodated dynamic changes in glacier area. Although glaciers in the Mica basin only cover 5% of the watershed, glacier ice melt contributes up to 25% and 35% of streamflow in August and September, respectively. The mean annual contribution of ice melt to total streamflow varied between 3 and 9% and averaged 6%. Glacier ice melt is particularly important during warm, dry summers following winters with low snow accumulation and early snowpack depletion. Although the sensitivity of streamflow to historic glacier area changes is small and within parameter uncertainties, our results suggest that glacier area changes have to be accounted for in future projections of late summer streamflow. Our approach provides an effective and widely applicable method to calibrate hydrologic models in glacier fed catchments, as well as to quantify the magnitude and timing of glacier melt contributions to streamflow.

The role of spatial variability in water inputs on runoff dynamics has generally not received as much research attention as topography and soils; however, the influence of topography and forest cover on snow surface energy exchanges can result in asynchronous snowmelt throughout a catchment, complicating the space–time patterns of runoff generation. This study investigates temporal variation in the relative importance of spatial controls on the occurrence, duration, and timing of shallow groundwater response, utilizing a highly distributed monitoring network in a snowmelt-dominated montane catchment in western Canada. The study findings indicate that deep-soil hydraulic conductivity is a first-order control on the spatial distribution of sites that generate shallow groundwater response versus sites that experience only deep percolation. Upslope contributing area and slope gradient are first-order controls on the duration of groundwater response during peak-flow, recession-flow, and low-flow periods. Shallow runoff response areas expand and contract throughout these periods and follow the general spatial patterns of topographic convergence. However, spatial controls on the timing, intensity, and quantity of snowmelt and controls on vertical versus lateral flux partitioning in the soil overwhelm the influence of topographic convergence on runoff patterns during early spring freshet periods. The study findings suggest that various topographic indices and topography-based rainfall runoff models would not likely be good predictors of runoff patterns in snowmelt-dominated montane catchments during early phases of the spring freshet, but would increase in importance as the freshet and post-freshet periods proceed.

Stream temperature dynamics during winter are less well studied than summer thermal regimes, but the winter season thermal regime can be critical for fish growth and development in coastal catchments. The winter thermal regimes of Pacific Northwest headwater streams, which provide vital winter habitat for salmonids and their food sources, may be particularly sensitive to changes in climate because they can remain ice-free throughout the year and are often located in rain-on-snow zones. This study examined winter stream temperature patterns and controls in small headwater catchments within the rain-on-snow zone at the Malcolm Knapp Research Forest, near Vancouver, British Columbia, Canada. Two hypotheses were addressed by this study: (1) winter stream temperatures are primarily controlled by advective fluxes associated with runoff processes and (2) stream temperatures should be depressed during rain-on-snow events, compared to rain-on-bare-ground events, due to the cooling effect of rain passing through the snowpack prior to infiltrating the soil or being delivered to the stream as saturation-excess overland flow. A reach-scale energy budget analysis of two winter seasons revealed that the advective energy input associated with hillslope runoff overwhelms vertical energy exchanges (net radiation, sensible and latent heat fluxes, bed heat conduction, and stream friction) and hyporheic energy fluxes during rain and rain-on-snow events. Historical stream temperature data and modelled snowpack dynamics were used to explore the influence of transient snow cover on stream temperature over 13 winters. When snow was not present, daily stream temperature during winter rain events tended to increase with increasing air temperature. However, when snow was present, stream temperature was capped at about 5 °C, regardless of air temperature. The stream energy budget modelling and historical analysis support both of our hypotheses. A key implication is that climatic warming may generate higher winter stream temperatures in the rain-on-snow zone due to both increased rain temperature and reduced cooling effect of snow cover.

Bridge Glacier is a lake-calving glacier in the Coast Mountains of British Columbia and has retreated over 3.55 km since 1972. The majority of this retreat has occurred since 1991. This retreat is substantially greater than what has been inferred from regional climate indices, suggesting that it has been driven primarily by calving as the glacier retreated across an overdeepened basin. In order to better understand the primary drivers of ablation, surface melt (below the equilibrium line altitude, ELA) and calving were quantified during the 2013 melt season using a distributed energy balance model (DEBM) and time-lapse imagery. Calving, estimated using areal change, velocity measurements, and assuming flotation were responsible for 23 % of the glacier's ablation below the ELA during the 2013 melt season and were limited by modest flow speeds and a small terminus cross-section. Calving and surface melt estimates from 1984 to 2013 suggest that calving was consistently a smaller contributor of ablation. Although calving was estimated to be responsible for up to 49 % of the glacier's ablation for individual seasons, averaged over multiple summers it accounted between 10 and 25 %. Calving was enhanced primarily by buoyancy and water depths, and fluxes were greatest between 2005 and 2010 as the glacier retreated over the deepest part of Bridge Lake. The recent rapid rate of calving is part of a transient stage in the glacier's retreat and is expected to diminish within 10 years as the terminus recedes into shallower water at the proximal end of the lake. These findings are in line with observations from other lake-calving glacier studies across the globe and suggest a common large-scale pattern in calving-induced retreat in lake-terminating alpine glaciers. Despite enhancing glacial retreat, calving remains a relatively small component of ablation and is expected to decrease in importance in the future. Hence, surface melt remains the primary driver of ablation at Bridge Glacier and thus projections of future retreat should be more closely tied to climate.

A growing body of research on stream thermal regimes has highlighted the importance of heat advection associated with surface water and groundwater interactions, such as hyporheic exchange, groundwater discharge, and hillslope throughflow inputs. Existing catchment models that predict stream temperature use a variety of approaches to estimate throughflow temperatures, but none has been evaluated against field measurements of throughflow temperature. In this study, throughflow temperatures were monitored over two winters at 50 locations adjacent to a headwater stream (11 ha catchment area) located in the rain-on-snow zone of the Pacific Northwest. Existing approaches to estimate throughflow temperature under or overpredicted throughflow temperatures by up to 5°C, or were unable to represent the influence of transient snow cover. Therefore, a conceptual-parametric model that is computationally efficient was developed that simulates hillslope hydrology and throughflow temperatures. The model structure includes an upslope reservoir that drains into a downslope reservoir that, in turn, drains into the stream. Vertical and lateral energy and water fluxes are simulated using simplified process representations. The model successfully predicts throughflow temperatures and highlights the dominant role of throughflow advection and the influence of snow cover on stream thermal regimes during high flow periods and rain-on-snow events.

We collected fish samples and measured physical habitat characteristics, including summer stream temperatures, at 156 sites in 50 tributary streams in two sampling areas (Upper Fraser and Thompson Rivers) in British Columbia, Canada. Additional watershed characteristics were derived from GIS coverages of watershed, hydrological and climatic variables. Maximum weekly average temperature (MWAT), computed as an index of summer thermal regime, ranged from 10 to 23 °C. High values of MWAT were associated with large, warm, low relief watersheds with a high lake influence. Measures of community similarity suggested that the fish community changed most rapidly through a lower transition zone at an MWAT of about 12 °C and an upper transition zone at an MWAT of about 19 °C. These results were confirmed using existing fisheries inventory data combined with predictions of MWAT from a landscape-scale regression model for the Thompson River watershed. For headwater sites in the Chilcotin River watershed (which drains into the middle Fraser River), the relative dominance of bull trout versus rainbow trout (based on inventory data) decreased with increasing predicted MWAT although the distinction was not as clear as for the Thompson River sites. The fish communities in these watersheds can be characterized in terms of very cold water (bull trout and some cold water species), cold water (salmonids and sculpins) and cool water (minnows and some cold water salmonids). The two transition zones (ca 12 and 19 °C) can be used to identify thresholds where small changes in stream temperature can be expected to lead to large changes in fish communities. Such clear, quantifiable thresholds are critical components of a management strategy designed to identify and protect vulnerable fish communities in streams where poor land use practices, alone or in combination with climatic change, can lead to changes in stream temperatures. Copyright © 2015 John Wiley & Sons, Ltd.

Salt dilution is a popular approach used for discharge measurement. This research focused on the procedure for determining the calibration factor (CFT) that is used to convert measured temperature corrected electrical conductivity to salt concentration for injection using dry salt. It is important to document the uncertainty in CFT because it contributes directly to uncertainty in the calculated discharge. Based on laboratory trials, it was found that the calibration should be performed as close to in situ stream temperature as possible to minimize errors. The discharge measurement and calibration procedure should be performed with the same probe to minimize uncertainty. Distilled water can be used instead of stream water for a calibration solution if an analytical correction is applied to account for differences in ionic composition of the water. The calibration factor can be determined with an uncertainty of less than ± 1% under “best-case” conditions, and the uncertainty may be as high as ± 4% under less favourable conditions. If calibration is not performed, CFT can be estimated from the relation between CFT and background temperature-corrected electrical conductivity (ECBG) with an uncertainty of about ± 2%, or estimated as a set value of 0.486 mg·cm·μS-1·L-1 with an uncertainty of about ± 2.8% for a properly calibrated probe. More testing should focus on streams with ECBG > 500 μS·cm-1, which were not well represented in this study.

Both the science and practice associated with watershed management involve the processing, presentation and analysis of quantitative information. In this article, the use of open source programming languages by watershed analysts is advocated. The R language, in particular, provides a rich set of tools for the types of data that are commonly encountered in watershed analysis. The utility of R is illustrated through three examples: intensity-duration-frequency analysis of rainfall data, baseflow separation, and watershed delineation and mapping.

The relations between groundwater level fluctuations and past climatic variations are explored using available data from the provincial observation well network, climate data and hydrometric data from the two major hydro-climatic zones of British Columbia. The majority of these well records start in the late 1970s or 1980s, providing only 20 to 30 years of record at most, and only a few of these wells monitor aquifers in “pristine” areas that reflect natural variability; the others have been influenced by human activity. Mean monthly groundwater levels were used to manually classify the wells into rainfall- and snowmelt-dominated seasonal cycles. Temporal trends of groundwater levels and a simple recharge measure were calculated using the non-parametric Spearman’s rank correlation coefficient over a common period from 1976–1999. This period was selected to coincide with the positive phase of the Pacific Decadal Oscillation (PDO) so as to eliminate the potential impact a shift in the PDO might have on grou...

In August 2003 a severe wildfire burnt 62% of Fishtrap Creek, a 158 km 2 watershed in central British Columbia, Canada. Streamflows were obtained for the period 1980–2010 and suspended sediment fluxes were determined for the period 2004–2010 for Fishtrap Creek and these were compared to data for nearby Jamieson Creek, which was not affected by the wildfire. Peak streamflows in Fishtrap Creek after the wildfire were not significantly higher than before the wildfire, although total annual runoff had increased. Perhaps the most important change in streamflows following the wildfire was that peak flows associated with the annual freshet occurred earlier in the year (by ca. 2 weeks). Following the wildfire, monthly total suspended sediment fluxes peaked in April in Fishtrap Creek and May in Jamieson Creek, which reflects the change in timing of peak streamflows in Fishtrap. Specific suspended sediment yields were low in the first year following the wildfire (2004), and peak values for the 2004–2010 monitoring period occurred in 2006. Average specific suspended sediment yields over the monitoring period were similar for both watersheds at 2.8 and 2.9 t km − 2 year − 1 for Fishtrap and Jamieson watersheds, respectively. The muted responses of streamflows and suspended sediment fluxes following this severe wildfire are due to the lack of winter precipitation and the low intensities of summer rainfall events in the first year following the wildfire. Greater winter precipitation and associated snowmelt in subsequent years coincided with vegetation recovery. The major changes in the wildfire-affected watershed were increased bank erosion and channel migration due to a loss of root strength and cohesion, which occurred 3–5 years after the fire. This work demonstrates that the hydrological and geomorphological responses of watersheds to wildfires are a function of the severity of the wildfire and the timing and nature of driving forces (i.e. rainfall intensity, winter precipitation and snowmelt) during the progression of vegetation recovery.

The objective of this study was to characterize the spatial variability of stream thermal regimes in British Columbia, Canada, with the specific goal of developing a predictive model to assist in provincial-scale assessment of fish habitat. It is part of a broader study to develop an approach to support the designation of “Temperature Sensitive Streams”, particularly in relation to the potential effects of forest harvesting and climate change. Stream temperature data were collected from researchers, consultants and government agencies. After checking for data quality, the annual maximum of a seven-day running average of mean daily water temperature (MWAT) was extracted for each station-year. A multiple regression model for the mean MWAT for each station was fitted for stations having basin areas between 1 and 104 km2. Predictor variables included the logarithm of catchment area, normal July–August air temperature for the location, the square root of the percentage of glacier cover in the catchment, the sq...

Hydrologic classification is useful for data organization, transfer of model parameters and estimation of hydrologic sensitivity to disturbance and climatic change. Stream-flow regime has frequently been used as a basis for classification, typically by mapping regimes defined by stream-flow data from a gauging network. As an alternative, we hypothesized that ecological classification systems can predict stream-flow regime because they are based on the same characteristics that control run-off generation (soils, climate and topography). A multivariate regression tree (MRT) was used to relate stream-flow regime to the fractional coverages of the Biogeoclimatic Ecological Classification (BEC) zones within the catchment for gauged streams in British Columbia, Canada. Although the MRT identified a realistic set of regimes, only a small number of BEC zones were used as predictors, reflecting bias in the gauging network. To avoid this bias, we used a water balance model to compute mean monthly stream flow for 932 ungauged basins in British Columbia that were generated with areas between 10 and 1000 km2; these monthly stream flows were used to train an MRT model based on BEC zone coverages. This model predicted the regime at gauged basins nearly as accurately as the water balance model for pluvial, nival and glacier-supported nival regimes. Difficulties occurred in smaller basins and in specific regions where the local BEC zones were not included as predictors. Coastal hybrid nivo–pluvial regimes were poorly predicted. With further development, ecological classification systems could have great value as a tool for hydrologic classification for both research and operational applications. Copyright © 2013 John Wiley & Sons, Ltd.

Stream temperature was recorded between 2002 and 2005 at four sites in a coastal headwater catchment in British Columbia, Canada. Shallow groundwater temperatures, along with bed temperature profiles at depths of 1 to 30 cm, were recorded at 10-min intervals in two hydrologically distinct reaches beginning in 2003 or 2004, depending on the site. The lower reach had smaller discharge contributions via lateral inflow from the hillslopes and fewer areas with upwelling (UW) and/or neutral flow across the stream bed compared to the middle reach. Bed temperatures were greater than those of shallow groundwater during summer, with higher temperatures in areas of downwelling (DW) flow compared to areas of neutral and UW flow. A paired-catchment analysis revealed that partial-retention forest harvesting in autumn 2004 resulted in higher daily maximum stream and bed temperatures but smaller changes in daily minima. Changes in daily maximum stream temperature, averaged over July and August of the post-harvest year, ranged from 1.6 to 3 °C at different locations within the cut block. Post-harvest changes in bed temperature in the lower reach were smaller than the changes in stream temperature, greater at sites with DW flow, and decreased with depth at both UW and DW sites, dropping to about 1 °C at a depth of 30 cm. In the middle reach, changes in daily maximum bed temperature, averaged over July and August, were generally about 1 °C and did not vary significantly with depth. The pre-harvest regression models for shallow groundwater were not suitable for applying the paired-catchment analysis to estimate the effects of harvesting. However, shallow groundwater was warmer at the lower reach following harvesting, despite generally cooler weather compared to the pre-harvest year. Copyright © 2012 John Wiley & Sons, Ltd.

This study focused on a headwater forest stream in coastal British Columbia, Canada. Air temperature, humidity and wind speed were measured at a height of 1.5 m above the stream and at a control site within a clearcut located approximately 1 km from the study stream, both before and after partial-retention forest harvesting along the stream. A specially designed evaporimeter was used to measure evaporation. Laboratory trials confirmed the reliability of evaporimeter measurements. Prior to harvesting, wind speeds were low and vapour pressure gradients above the stream were weak, leading to low rates of evaporation. Following harvesting, the decreased shading and increased ventilation over the stream led to higher wind speeds, lower vapour pressures, higher daily maximum air temperatures, higher stream temperatures and surface vapour pressures, and higher rates of evaporation. A wind function that has been used to estimate stream evaporation in a number of previous studies, but which had not been compared to measured stream evaporation, was found to overestimate evaporation. An empirical wind function fitted to the evaporimeter data differed from three others that have been derived for predicting stream evaporation in that our data did not support the inclusion of an intercept. It is hypothesized that the lack of an intercept for our data reflects the strongly stable conditions over the stream. Further research should measure stream evaporation over a broad range of streamflow and meteorological conditions, with particular attention to the role of atmospheric stability.

Uncertainty in hydrometric data is a fact of life. Basic assumptions about the nature of this uncertainty are necessary in every analysis of hydrometric data, and an understanding of the variability of uncertainty can facilitate the effective use of hydrologic information. For most of the twentieth century there has been little change in hydrometric methods and many analysts explicitly or implicitly assume that the uncertainty has not changed over the period of record. We argue that there is substantial variability in the magnitude of uncertainty in published streamflow records that is not transparent to data users. Quantifying uncertainty is particularly important in the context of the current changes in hydrometric technology and in the increasing integration of data sets from multiple providers. We recommend best practices for identifying uncertainty in field notes and propagating that observational uncertainty through the data production process. We suggest both field and reanalysis studies that could...

Stream surface albedo was measured at a location downstream of Place Glacier, Canada, in a steep bouldery channel. Portions of the water surface were visibly aerated as a result of the cascading flow even at lower discharges; at high flows, the stream was near-continuous whitewater. Albedo generally increased with discharge, from around 0.1 at the lower flows to 0.4 at the highest flows. This increase is consistent with the known effect of aeration on the reflectance of water. This discharge dependence of albedo needs to be accounted for in physically based models for predicting stream temperature to avoid biased predictions of net radiation. For steep proglacial streams that experience decreasing late-summer flows as a result of ongoing and future glacier recession, the associated decrease in albedo could promote higher stream temperatures, in addition to the effects of reduced flow depth and velocity. Copyright © 2011 John Wiley & Sons, Ltd.

Data from a paired-catchment study in south coastal British Columbia, Canada, were analyzed to assess the thermal effects of clearcut harvesting with no riparian buffer on a fish-bearing headwater stream. The approach used time series of daily mean water temperatures for East Creek (control) and A Creek (treatment), both before and after harvest. Statistical models were developed to predict (a) what the temperatures would have been in the post-harvest period had harvesting not occurred, and (b) what temperatures would have been in the pre-harvest period had harvesting already occurred. The Wisconsin Bioenergetics Model was used to simulate growth of coastal cutthroat trout (Oncorhynchus clarki clarki) for the first year following fry emergence using the predicted and observed stream temperatures to generate scenarios representing with-harvest and no-harvest thermal regimes. A Monte Carlo approach was used to quantify the effects of uncertainty associated with the regression models on predicted stream temperature and trout growth. Summer daily mean temperatures in the with-harvest scenario were up to $$5^{\circ}\hbox{C}$$ higher than those for the no-harvest scenario. Harvesting-induced warming reduced growth rates during summer, but increased growth rates during autumn and spring. In the with-harvest scenario, trout were 0.2–2.0 g (absolute weight) smaller throughout the winter period than in the no-harvest scenario. However, the bioenergetic simulations suggest that trout growth may be more sensitive to potential changes in food supply following harvesting than to direct impacts of stream temperature changes.

This study was motivated by an interest in understanding the potential effects of climate change and glacier retreat on late summer water temperatures in alpine areas. Fieldwork was carried out between July and September 2007 at Place Lake, located below Place Glacier in the southern Coast Mountains of British Columbia. Place Lake has an area of 72 000m2, a single inlet and outlet channel, and an approximate residence time of 4 days. Warming between the inlet and outlet of the lake ranged up to 3 C and averaged 1.8 C, which exceeds the amount of warming that occurred over the 1 km reach of Place Creek between the lake outlet and tree line. Over a 23-day period, net radiation totalled about 210 MJ m–2, with sensible heat flux adding another 56 MJm-2. The latent heat flux consumed about 8% of the surface heat input. The dominant heat sink was the net horizontal advection associated with lake inflow and outflow. Early in the study period, temperatures between the surface and 6-m depth were dominantly at or above 4 C and were generally neutral to thermally stable, whereas temperatures decreased with depth below 6m and exhibited irregular sub-diurnal variations. The maximum outflow temperature of almost 7 C occurred in this period. We hypothesize that turbidity currents associated with cold, sediment-laden glacier discharge formed an underflow and influenced temperatures in the deeper portion of the lake but did not mix with the upper layers. Later in the study period, the lake was dominantly well mixed with some near-surface stability associated with nocturnal cooling. Further research is required to examine the combined effects of sediment concentrations and thermal processes on mixing in small proglacial lakes to make projections of the consequences of glacier retreat on alpine lake and stream temperatures.

Interactions between mobile stream water and transient storage zones have been the subject of careful attention for decades. However, few studies have considered explicitly the influence of water exchange between the channel and neighbouring hydrological units when modelling transient storage processes, especially the lateral inflow coming from hillslope contributions and outflow to a deep aquifer or to hyporheic flow paths extending beyond the study reach. The objective of this study was to explore the influence of different conceptualizations of these hydrologic exchanges on the estimation of transient storage parameters. Slug injections of sodium chloride (NaCl) were carried out in eight contiguous reaches in the Cotton Creek Experimental Watershed (CCEW), located in south-east British Columbia. Resulting breakthrough curves were subsequently analysed using a Transient Storage Model (TSM) in an inverse modelling framework. We estimated solute transport parameters using three distinct, hypothetical spatial patterns of lateral inflow and outflow, all based on variations of the same five-parameter model structure. We compared optimized parameter values to those resulting from a distinct four-parameter model structure meant to represent the standard application of the TSM, in which only lateral inflow was implemented for net gaining reaches or only lateral outflow for net losing reaches. In the five-parameter model, solute mass was stored predominantly in the transient storage zone and slowly released back to the stream. Conversely, solute mass was predominantly removed from the stream via flow losses in the four-parameter model structure. This led to contrasting estimates of solute transport parameters and subsequent interpretation of solute transport dynamics. Differences in parameter estimates across variations of the five-parameter model structure were small yet statistically significant, except for the transient storage exchange rate coefficient α , for which unique determination was problematic. We also based our analysis on F med 200 , the fraction of median transport time due to transient storage. Differences across configurations in F med 200 estimates were consistent but small when compared to the variability of F med 200 among reaches. Optimized parameter values were influenced dominantly by the model structure (four versus five parameters) and then by the conceptualization of spatial arrangement of lateral fluxes along the reach for a set model structure. When boundary conditions are poorly defined, the information contained in the stream tracer breakthrough curve is insufficient to identify a single, unambiguous model structure representing solute transport simulations. Investigating lateral fluxes prior to conducting a study on transient storage processes is necessary, as assuming a certain spatial organization of these fluxes might set ill-defined bases for inter-reach comparisons. Given the difficulty in quantifying the spatial patterns and magnitudes of lateral inputs and outputs, we recommend small-scale laboratory tracer experiments with well-defined and variable boundary conditions as a complement to field studies to provide new insights into stream solute dynamics.

The objective of this study was to analyse stream temperature variability during summer in relation to both surface heat exchanges and reach-scale hydrology for two hydrogeomorphically distinct reaches. The study focused on a 1·5-km wildfire-disturbed reach of Fishtrap Creek located north of Kamloops, British Columbia. Streamflow measurements and longitudinal surveys of electrical conductivity and water chemistry indicated that the upper 750 m of the study reach was dominated by flow losses. A spring discharged into the stream at 750 m below the upper reach boundary. Below the spring, the stream was neutral to losing on three measurement days, but gained flow on a fourth day that followed a rain event. Continuous stream temperature measurements typically revealed a downstream warming along the upper 750 m of the study reach on summer days, followed by a pronounced cooling associated with the spring, with little downstream change below the spring. Modelled surface energy exchanges were similar over the upper and lower sub-reaches, and thus cannot explain the differences in longitudinal temperature patterns. Application of a Lagrangian stream temperature model provided reasonably accurate predictions for the upper sub-reach. For the lower sub-reach, accurate prediction required specification of concurrent flow losses and gains as a hydrological boundary condition. These findings are consistent with differences in the hydrogeomorphology of the upper and lower sub-reaches. The modelling exercise indicated that substantial errors in predicted stream temperature can occur by representing stream-surface exchange as a reach-averaged one-directional flux computed from differences in streamflow between the upper and lower reach boundaries. Further research should focus on reliable methods for quantifying spatial variations in reach-scale hydrology. Copyright © 2010 John Wiley & Sons, Ltd.

The Pacific Decadal Oscillation (PDO) is a large-scale climate system feature that influences the surface climate and hydrology of western North America. In this paper, we review the literature describing the PDO and demonstrate its effects on temperature, precipitation, snowfall, glacier mass balance, and streamflow with a focus on western Canada, and particularly British Columbia. We review how the PDO index was developed and discuss other North Pacific climate patterns that resemble the PDO. The impacts of PDO on glacier mass balance and streamflow from retrospective studies are also reviewed and illustrated with specific examples from BC. We assess the current state of knowledge regarding the PDO and provide a critical assessment of its use in hydroclimatology. This information should provide insight on the sensitivity of projects to climatic variability.