Incipient motion conditions for single rock fragments in simulated rill flow

Laboratory experiments using simulated tilting rill channels with different bed roughnesses have been conducted to determine hydraulic conditions for incipient motion of single rock fragments lying on rill beds. These threshold conditions can be described by a modified Shields' entrainment parameter (θ_m) and relative depth (Z), adjusted for rill conditions. The negative relation found between θ_m and Z also fits steep (stream) channel data collected by other investigators. Plots of the classical Shields' entrainment parameter (θ) versus both relative depth and relative particle size also showed clear negative trends. Special care has been taken to avoid false relations that can arise between compound parameters. This occurs when the latter are not independent from each other because of their definitions.     

Channel stability in steep gravel–cobble streams is controlled by the coarse tail of the bed material distribution

Researchers have associated channel-forming flows with reach-average shear stresses close to the entrainment threshold for the surface D₅₀ . We conducted experiments using a model of a generic steep, gravel–cobble stream to test this association. Our results suggest that channel-forming flows fully mobilize the D₅₀ , and produce shear stresses close to the entrainment threshold for the largest grains in the bed. The channel dimensions were set by flows capable of mobilizing between 85% and 90% of the bed surface, which produced a brief period of lateral instability lasting about 1 h, followed by a prolonged period of relative stability during which modest adjustments occurred, but during which the reach-average hydraulics remained about the same. The adjustments during the unstable phase of the experiments are characterized by rapid bank erosion, extensive deposits on the channel bed and a restructuring of the major morphologic elements of the stream. The adjustments during the stable phase of the experiments involved barform migration and bed surface coarsening but did not appreciably modify the physical template established by the end of the unstable phase. The behaviour we observed is not consistent with the concept of a dynamic equilibrium associated with a formative flow that is just capable of entraining the bed surface D₅₀ . Instead, it suggests that rapid adjustments occur once a stability threshold is exceeded, which creates a template that constrains channel activity until another event drives the system across the stability threshold, and re-sets the template. While we believe that it is probably too simplistic to associate a channel-forming discharge with the entrainment threshold for a single grain size, our results suggest that the D₉₅ is a more logical choice than the D₅₀ © 2020 John Wiley & Sons, Ltd.     

A rational sediment transport scaling relation based on dimensionless stream power

The concept of stream channel grade – according to which a stream channel reach will adjust its gradient, S, in order to transport the imposed sediment load having magnitude Q_b and characteristic grain size D_b, with the available discharge Q (Mackin, 1948 , Geological Society of America Bulletin 59 : 463–512; Lane, 1955 , American Society of Civil Engineers, Proceedings 81 : 1–17) is one of the most influential ideas in fluvial geomorphology. Herein, we derive a scaling relation that describes how externally imposed changes in either Q_b or Q can be accommodated by changes in the channel configuration, described by the energy gradient, mean flow depth, characteristic grain size and a parameter describing the effect of bed surface structures on grain entrainment. One version of this scaling relation is based on the dimensionless bed material transport parameter (W*) presented by Parker and Klingeman ( 1982 , Water Resources Research 18 : 1409–1423). An equivalent version is based on a new dimensionless transport parameter (E*) using dimensionless unit stream power. This version is nearly identical to the relation based on W*, except that it is independent of flow resistance. Both versions of the scaling relation are directly comparable to Lane's original relation. In order to generate this stream power‐based scaling relation, we derived an empirical transport function relation relating E* to dimensionless stream power using data from a wide range of stable, bed load‐dominated channels: the form of that transport function is based on the understanding that, while grain entrainment is related to the forces acting on the bed (described by dimensionless shear stress), sediment transport rate is related to the transfer of momentum from the fluid to the bed material (described by dimensionless stream power). Copyright © 2010 John Wiley & Sons, Ltd.     

Dimensionless critical shear stress evaluation from flume experiments using different gravel beds

Flume experiments were conducted using four different gravel beds (D₅₀ + 12–39 mm) and a range of marked particles (10–65 mm). The shear stresses were evaluated from friction velocities, when initial movement of marked particles occurred. Two kinds of equations were produced: first for the threshold of initial movement, and second for generalized movement. Equations of the type 0_c + a(D_i/D₅₀)^b, as proposed by Andrews (1983) are applicable even if the material is relatively well sorted. However, the values of a and b are lower (respectively 0·050 and -0·70) for initial movement. Generalized movement requires a higher shear stress (a + 0·068 and b + -0·80). D₉₀ of the bed material and y₀ (the bed roughness parameter) were also used as reference values in place of D₅₀. They produced lower values than in natural streams, mainly owing to the fact that the material used in the flume is better sorted: clusters are less well developed and the bed roughness is lower.     

Channel flow competence and sediment transport in upland streams in southeast Australia

Bedload transport data from planebed and step‐pool reach types are used to determine grain size transport thresholds for selected upland streams in southeast Australia. Morphological differences between the reach types allow the effects of frictional losses from bedforms, microtopography and bed packing to be incorporated into the dimensionless critical shear stress value. Local sediment transport data are also included in a regime model and applied to mountain streams, to investigate whether empirical data improve the delineation of reach types on the basis of dimensionless discharge per unit width (q*) and dimensionless bedload transport (q_b*). Instrumented planebed and step‐pool sites are not competent to transport surface median grains (D_50s) at bankfull discharge (Q_bf). Application of a locally parametrized entrainment equation to the full range of reach types in the study area indicates that the majority of cascades, cascade‐pools, step‐pools and planebeds are also not competent at Q_bf and require a 10 year recurrence interval flood to mobilize their D_50s. Consequently, the hydraulic parameters of the regime diagram, which assume equilibrium conditions at bankfull, are ill suited to these streams and provide a poor basis of channel delineation. Modifying the diagram to better reflect the dominant transported bedload size (equivalent to the D₁₆ of surface sediment) made only slight improvements to reach delineation and had greatest effect on the morphologies with smaller surface grain sizes such as forced pool‐riffles and planebeds. Likewise, the Corey shape factor was incorporated into the regime diagram as an objective method for adjusting a base dimensionless critical shear stress (τ*_c50b) to account for lithologically controlled grain shape on bed packing and entrainment. However, it too provided only minor adjustments to reach type delineation. Copyright © 2007 John Wiley & Sons, Ltd.     

The impact of stress history on bed structure

Recent research has started to focus on how prolonged periods of sub‐threshold flows may be capable of imparting structural changes that contribute to increased bed stability. To date, this effect (termed ‘stress history’) has been found to be significant in acting to increase a bed's critical shear stress at entrainment threshold. However, it is supported by only limited, qualitative and often speculative information on the mechanisms of this stabilization process in grade‐specific studies. As such, this paper uses high resolution laser scanning to quantitatively ascertain the granular mechanics underpinning the relationship between stress history and entrainment threshold for beds of a range of grain size distributions. Employing a bed slope of 1/200, three grain size distributions with median grain sizes (D₅₀) of 4·8 mm [uniform (σ_g = (D₈₄/D₁₆)^0.5 = 1·13; bimodal (σ_g = 2·08); and, unimodal (σ_g = 1·63)] were exposed to antecedent stress histories of 60 and 960 minutes duration. Antecedent shear stress magnitude was set at 50% of the critical shear stress for the D₅₀ when no stress history period was employed. Two laser displacement scans of the bed surface (approximate area 100 mm × 117 mm) were taken, one prior to the antecedent period and one after this period, so that changes to surface topography could be quantified (resolution of x = 0·10 mm, y = 0·13 mm and z = 0·24 mm). Rearrangement of bed surface structure is described using statistical analysis and two‐dimensional (2D) semi‐variograms to analyse scaling behaviour. Results reveal vertical settlement, changes to bed roughness and particle repositioning. However, the bed grain size distribution influences the relative importance of each mechanism in determining stress history induced bed stability; this is the focus of discussion in this paper. Copyright © 2012 John Wiley & Sons, Ltd.     

Aerodynamic entrainment thresholds and dislodgement rates on impervious and permeable beds

In the literature it has been suggested that on permeable, granular beds, both the threshold and rate of aerodynamic entrainment may be affected significantly by seepage flows into and out of the bed induced by fluctuating pressures in the overlying turbulent boundary layer. Using a range of grain sizes and flow conditions, the series of laboratory experiments reported here compares directly the aerodynamic entrainment of loose grains overlying fixed permeable sediment beds with that occuring over fixed impervious beds. For a given granular material, no significant differences in entrainment dynamics on the two types of bed were observed and in the range of flow conditions examined both the threshold shear velocity (U _*T) and the aerodynamic entrainment coefficient (k) were found to be independent of bed permeability.     

Relation between flow,surface‐layer armoring and sediment transport in gravel‐bed rivers

This study investigates trends in bed surface and substrate grain sizes in relation to reach‐scale hydraulics using data from more than 100 gravel‐bed stream reaches in Colorado and Utah. Collocated measurements of surface and substrate sediment, bankfull channel geometry and channel slope are used to examine relations between reach‐average shear stress and bed sediment grain size. Slopes at the study sites range from 0·0003 to 0·07; bankfull depths range from 0·2 to 5 m and bankfull widths range from 2 to 200 m. The data show that there is much less variation in the median grain size of the substrate, D_50s, than there is in the median grain size of the surface, D₅₀; the ratio of D₅₀ to D_50s thus decreases from about four in headwater reaches with high shear stress to less than two in downstream reaches with low shear stress. Similar trends are observed in an independent data set obtained from measurements in gravel‐bed streams in Idaho. A conceptual quantitative model is developed on the basis of these observations to track differences in bed load transport through an idealized stream system. The results of the transport model suggest that downstream trends in total bed load flux may vary appreciably, depending on the assumed relation between surface and substrate grain sizes. Copyright © 2007 John Wiley & Sons, Ltd.     

A multicatchment analysis of headwater and downstream temperature effects from contemporary forest harvesting

Stream temperature is a key physical water‐quality parameter, controlling many biological, chemical, and physical processes in aquatic ecosystems. Maintenance of cool stream temperatures during summer is critical for high‐quality aquatic habitat. As such, transmission of warm water from small, nonfish‐bearing headwater streams after forest harvesting could cause warming in downstream fish‐bearing stream reaches with negative consequences. In this study, we evaluate (a) the effects of contemporary forest management practices on stream temperature in small, headwater streams, (b) the transmission of thermal signals from headwater reaches after harvesting to downstream fish‐bearing reaches, and (c) the relative role of lithology and forest management practices in influencing differential thermal responses in both the headwater and downstream reaches. We measured summer stream temperatures both preharvest and postharvest at 29 sites—12 upstream sites (4 reference, 8 harvested) and 17 downstream sites (5 reference, 12 harvested)—across 3 paired watershed studies in western Oregon. The 7‐day moving average of daily maximum stream temperature (T_7DAYMAX) was greater during the postharvest period relative to the preharvest period at 7 of the 8 harvested upstream sites. Although the T_7DAYMAX was generally warmer in the downstream direction at most of the stream reaches during both the preharvest and postharvest period, there was no evidence for additional downstream warming related to the harvesting activity. Rather, the T_7DAYMAX cooled rapidly as stream water flowed into forested reaches ~370–1,420 m downstream of harvested areas. Finally, the magnitude of effects of contemporary forest management practices on stream temperature increased with the proportion of catchment underlain by more resistant lithology at both the headwater and downstream sites, reducing the potential for the cooling influence of groundwater.     

Entrainment,displacement and transport of tracer gravels

The mass and size distribution of grain entrainment per unit bed area may be measured by replacing a volume of the bed with tracer gravels and observing the mass difference before and after a transport event. This measure of spatial entrainment is relevant to any process involving size-selective exchange of sediment between transport and bed and may be directly used in calculations of sediment transport rate using an elementary relation for fractional transport components presented here. This relation provides a basis for evaluating tracer data collected by different methods and may be used to provide physical insight regarding the expected behaviour of tracer grains. The variation with grain size of total displacement length L_ti depends on the degree of mobilization of the individual fractions on the bed surface: L_ti is independent of D_i for smaller, fully mobile sizes and decreases rapidly with D_i for larger fractions in a state of partial transport (in which a portion of the surface grains remain immobile through the flow event). The boundary between fully and partially mobile grain sizes increases with flow strength. These inferences are supported by values of L_ti calculated from flume experiments and provide a physical explanation for a summary relation between L_ti and D_i based on field data. © 1997 John Wiley & Sons, Ltd.     

When does bedload transport begin in steep boulder‐bed streams?

The paper presents the results of field measurements of critical conditions for bedload motion conducted in the Rio Cordon, a steep boulder‐bed stream in the Italian Alps, under conditions of high Reynolds numbers and low relative submergence poorly explored before. Two methods have been used to determine threshold of motion: the displacement of marked clasts and the flow competence approach, which uses the largest grain size diameter transported by each flood event. The high values of confirm the great relevance of non‐bedload effective shear stresses in step–pool streams given by the additional form drag associated with this morphology. Relative submergence effects on the dimensionless critical shear stress have been quantified by considering the relative submergence ratio R_h/D₈₄, and the major effect of relative size on the mobility of each particle in steep, widely graded bed mixtures has been evaluated. Finally, the dimensionless critical unit discharge has also been adopted in the regression equations as the critical hydraulic parameter, because it may represent an easier parameter to use than the critical shear stress for steep, rough mountain rivers. Copyright © 2006 John Wiley & Sons, Ltd.     

The impact of inter-flood duration on non-cohesive sediment bed stability

Limited field and flume data suggests that both uniform and graded beds appear to progressively stabilize when subjected to inter-flood flows as characterized by the absence of active bedload transport. Previous work has shown that the degree of bed stabilization scales with duration of inter-flood flow, however, the sensitivity of this response to bed surface grain size distribution has not been explored. This article presents the first detailed comparison of the dependence of graded bed stability on inter-flood flow duration. Sixty discrete experiments, including repetitions, were undertaken using three grain size distributions of identical D₅₀ (4.8 mm); near-uniform (σ_g = 1.13), unimodal (σ_g = 1.63) and bimodal (σ_g = 2.08). Each bed was conditioned for between 0 (benchmark) and 960 minutes by an antecedent shear stress below the entrainment threshold of the bed (τ*_c50). The degree of bed stabilization was determined by measuring changes to critical entrainment thresholds and bedload flux characteristics. Results show that (i) increasing inter-flood duration from 0 to 960 minutes increases the average threshold shear stress of the D₅₀ by up to 18%; (ii) bedload transport rates were reduced by up to 90% as inter-flood duration increased from 0 to 960 minutes; (iii) the rate of response to changes in inter-flood duration in both critical shear stress and bedload transport rate is non-linear and is inversely proportional to antecedent duration; (iv) there is a grade dependent response to changes in critical shear stress where the magnitude of response in uniform beds is up to twice that of the graded beds; and (v) there is a grade dependent response to changes in bedload transport rate where the bimodal bed is most responsive in terms of the magnitude of change. These advances underpin the development of more accurate predictions of both entrainment thresholds and bedload flux timing and magnitude, as well as having implications for the management of environmental flow design. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Self-diffusion studies of pore fluids in unconsolidated sediments by PFG NMR

The self-diffusion of water and hexadecane in medium and coarse sands from glacial sand deposits in central Germany were investigated by pulsed field gradient nuclear magnetic resonance (PFG NMR). Due to the restriction of the diffusion path at the pore/grain interface, the measured apparent self-diffusion coefficients (D(Δ)) in the pore space depend on the observation time (Δ) in the PFG NMR experiment. Although the bulk self-diffusion coefficients of water and hexadecane differ by about one order of magnitude, the apparent self-diffusion coefficients in the pore space obey the same characteristic time-behaviour, which depends only on geometrical properties of the pore system. Using the “short-time diffusion” model, surface-to-volume (S/V) ratios and inherent self-diffusion coefficients (D₀) of the pore fluids were extracted from these diffusion measurements. The S/V ratios obtained are independent of the pore fluid used and agree with known geometrical properties of the sand grains. Moreover, the D₀ values are consistent with the corresponding bulk self-diffusion coefficients measured separately. In contrast to these results of PFG NMR, simultaneous investigations of longitudinal (T₁) nuclear magnetic relaxation reveal that the relaxation time of the pore fluid is a less suitable parameter for a quantitative estimation of geometrical properties of the pore/grain interface in these unconsolidated sediments since it depends on chemical properties of the fluid/grain interface.     

Woody buffer effects on water temperature: The role of spatial configuration and daily temperature fluctuations

Water temperature is a key driver for riverine biota and strongly depends on shading by woody riparian vegetation in summer. While the general effects of shading on daily maximum water temperature T_max are well understood, knowledge gaps on the role of the spatial configuration still exist. In this study, the effect of riparian buffer length, width, and canopy cover (percentage of buffer area covered by woody vegetation) on T_max was investigated during summer baseflow using data measured in seven small lowland streams in western Germany (wetted width 0.8–3.7 m). The effect of buffer length on T_max differed between downstream cooling and heating: T_max approached cooler equilibrium conditions after a distance of 0.4 km (~45 min travel-time) downstream of a sharp increase in canopy cover. In contrast, T_max continued to rise downstream of a sharp decrease in canopy cover along the whole 1.6 km stream length investigated. The effect of woody vegetation on T_max depended on buffer width, with changes in canopy cover in a 10 m wide buffer being a better predictor for changes in T_max compared to a 30 m buffer. The effect of woody vegetation on T_max was linearly related to canopy cover but also depended on daily temperature range T_range, which itself was governed by cloudiness, upstream canopy cover, and season. The derived empirical relationship indicated that T_max was reduced by −4.6°C and increased by +2.7°C downstream of a change from unshaded to fully shaded conditions and vice versa. This maximum effect was predicted for a 10 m wide buffer at sunny days in early summer, in streams with large diel fluctuations (large T_range). Therefore, even narrow woody riparian buffers may substantially reduce the increase in T_max due to climate change, especially in small shallow headwater streams with low baseflow discharge and large daily temperature fluctuations.     

Calculation of magnetic field turbulent diffusivities and α-effect

Abstract

First the exact numerical solutions of DIA system of equations describing the transportation of magnetic field in an infinite medium are presented. It is assumed that the turbulence is stationary, homogeneous, isotropic and incompressible. The spectra of turbulence of δ-type and Kolmogorov's type were used. The steady state values of magnetic field diffusivity D_T and the α-effect coefficient α _T were calculated for various values of space-scale and lifetimes of these spectra and the spectra of helicity. Also investigated is the dependence of D_T and α _T on the degree of helicity. The corrections to the α _T -coefficient due to the contribution of four-order velocity correlators are given. The results are compared with those due to the self-consistent technique.     

Fractal approach of the temporal earthquake distribution in the Hellenic arc-trench system

The time clustering of earthquakes occurring in the Hellenic arc-trench system is quantitatively analyzed by means of the fractal dimension,D, of their time distribution in the time intervals of 1950–1985 (M _s >-4.5) and 1964–1985 (M _s 4.0). The results obtained imply that scale-invariant clustering holds over very large scale lengths of time,T, with 2²–2⁸T (in min) 2²⁰–2²², depending on the seismotectonic segment considered. In all segments a common feature is the relation between theD ₁,D ₂ andD ₃-values found for shallow, intermediate-depth and all-depth shocks, respectively:D ₃>D₁>D₂. TheD-values found for shallow shocks range between 0.137 and 0.191 with the exception of the Ionian Islands and Cretan segments where anomalously high values (D=0.221–0.251) have been determined. We discuss possible seismotectonic interpretations of the results.     

Horizontal hydraulic conductivity of shallow streambed sediments and comparison with the grain‐size analysis results

Streambed horizontal hydraulic conductivity (K_h) has a substantial role in controlling exchanges between stream water and groundwater. We propose a new approach for determining K_h of the shallow streambed sediments. Undisturbed sediment samples were collected using tubes that were horizontally driven into streambeds. The sediment columns were analysed using a permeameter test (PT) on site. This new test approach minimizes uncertainties due to vertical flow in the vicinity of test tube and stream stage fluctuations in the computation of the K_h values. Ninety‐eight PTs using the new approach were conducted at eight sites in four tributaries of the Platte River, east‐central Nebraska, USA. The K_h values were compared with the nondirectional hydraulic conductivity values (K_g) determined from 12 empirical grain‐size analysis methods. The grain‐size analysis methods used the same sediment samples as K_h tests. Only two methods, the Terzaghi and Shepherd methods, yielded K_g values close to the K_h values. Although the Sauerbrei method produced a value relatively closer to K_h than other nine grain‐size analysis methods, the values from this method were not as reliable as the methods of Terzaghi and Shepherd due to the inconsistent fluctuation of the average estimates at each of the test sites. The Zunker, Zamarin, Hazen, Beyer, and Kozeny methods overestimated K_h, while the Slichter, US Bureau of Reclamation (USBR), Harleman, and Alyamani and Sen methods underestimated K_h. Any of these specific grain‐size methods might yield good estimates of streambed K_h at some sites, but give poor estimates at other sites, indicating that the relationship between K_g and K_h is significantly site dependent in our study. Copyright © 2011 John Wiley & Sons, Ltd.     

Predictive equation for longitudinal dispersion coefficient

Dye tracing field data were collected in small, steep streams in Ontario and used to calculate longitudinal dispersion coefficients for these headwater streams. A predictive equation for longitudinal dispersion coefficient is developed using combined data sets from five steeper head – water streams and 24 milder and larger rivers. The predictive equation relates the longitudinal dispersion coefficient to hydraulic and geometric parameters of the stream and has been developed using multiple regression analysis. The newly developed equation shows impressive accuracy of predictions for longitudinal dispersion coefficient (R² = 0.86, RMSE = 25, Nash–Sutcliffe coefficient E_ns = 0.86 and Index of Agreement D = 0.96) for both small, steep headwater streams as well as large, mild rivers. The Froude number has been introduced as a third key parameter to capture the effect of slope of the reach – in addition to the aspect ratio and bed material surface roughness – on the longitudinal dispersion coefficient. The pronounced improvement in the accuracy of the prediction is due to the addition of the Froude number to capture the effect of the slope of the reach on longitudinal dispersion coefficient. Copyright © 2014 John Wiley & Sons, Ltd.     

On the informative value of the largest sample element of log-Gumbel distribution

Extremes of stream flow and precipitation are commonly modeled by heavytailed distributions. While scrutinizing annual flow maxima or the peaks over threshold, the largest sample elements are quite often suspected to be low quality data, outliers or values corresponding to much longer return periods than the observation period. Since the interest is primarily in the estimation of the right tail (in the case of floods or heavy rainfalls), sensitivity of upper quantiles to largest elements of a series constitutes a problem of special concern. This study investigated the sensitivity problem using the log-Gumbel distribution by generating samples of different sizes (n) and different values of the coefficient of variation by Monte Carlo experiments. Parameters of the log-Gumbel distribution were estimated by the probability weighted moments (PWMs) method, method of moments (MOMs) and maximum likelihood method (MLM), both for complete samples and the samples deprived of their largest elements. In the latter case, the distribution censored by the non-exceedance probability threshold, F _T , was considered. Using F _T instead of the censored threshold T creates possibility of controlling estimator property. The effect of the F _T value on the performance of the quantile estimates was then examined. It is shown that right censoring of data need not reduce an accuracy of large quantile estimates if the method of PWMs or MOMs is employed. Moreover allowing bias of estimates one can get the gain in variance and in mean square error of large quantiles even if ML method is used.     

Characterizing stream temperature hysteresis in forested headwater streams

Stream temperature (T_s) is a key water quality parameter that controls several biological, ecological, and chemical processes in aquatic systems. In forested headwaters, exchanges of energy across air-water-streambed interfaces may influence T_s regimes, especially during storm events as the sources of runoff change over space and time. Analysis of the hysteretic behaviour of T_s during storm events may provide insights into rainfall-runoff responses, but such relationships have not been thoroughly investigated. As such, our objectives were to (a) quantify the variability of stream temperature hysteresis across seasons in different sub-regions and (b) investigate the relationship between the hysteretic response and catchment characteristics. T_s hysteresis during storm events was assessed based on the hysteresis index (HI), which describes the directionality of hysteresis loops, and the temperature response index (TRI), which indicates whether T_s increased or decreased during a storm event. We analysed T_s data from 10 forested headwater reaches in two sub-regions (McGarvey and West Fork Tectah) in Northern California. We also performed a clustering analysis to examine the relationship amongst HI, TRI, topographic metrics, and meteorological characteristics of the study areas. Overall, the hysteretic behaviour of T_s varied across seasons—the greatest HI occurred during spring and summer. Interestingly, in the McGarvey streams the variability in T_s hysteresis co-varied strongly with topographic metrics (i.e., upslope accumulative area, average channel slope, topographic wetness index). Comparatively, in West Fork Tectah the variability of T_s hysteresis co-varied most strongly with meteorological metrics (i.e., antecedent rainfall events, solar radiation, and air temperature). Variables such as the gradient between stream and air temperatures, slope, and wetted width were significant for both sub-regional hysteretic patterns. We posit that the drivers of T_s response during storms are likely dependent on catchment physiographic characteristics. Our study also illustrated the potential utility of stream temperature as a tracer for improving the understanding of hydrologic connectivity and shifts in the dominant runoff contributions to streamflow during storm events.