Stream and bed temperature variability in a coastal headwater catchment: influences of surface‐subsurface interactions and partial‐retention forest harvesting

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.     

Stream temperature,specific conductance and runoff process in mountain watersheds

Stream temperature ranged from 3 to 4°C at an experimental site during snowmelt on Hokkaido Island, Japan, which provided direct evidence of major contributions of subsurface water to stream water. In contrast, stream temperatures during rainstorms in summer decreased gradually after stream flow peaked, attaining a nearly constant temperature ranging from 9 to 11°C. During storm flow recession, stream temperatures during summer or snowmelt were similar to the soil temperature at 1·8 m below the land surface, suggesting that subsurface water contributions to stream flow are derived from this depth. The hygrographs during two rainstorms, August 1987 and September 1989, were separated using temperature. The stream temperature was assumed to depend on the mixing of surface flow, having a temperature ranging from that of rainfall to that of shallow (50 cm deep) soil water, and subsurface flow, having the temperature of the soil at 1·8 m below the land surface. Subsurface flow was estimated to contribute 85–90% of the total stream flow during each rainstorm. A two‐component hydrograph separation was also evaluated using specific conductance. Runoff contributions from the two sources for the temperature and specific conductance analysis were similar. Analysis of the temperature and conductance–discharge hysteresis loop, and of individual flow components for storm hygrographs, provide a general picture of the runoff process in the experimental basin. Copyright © 1999 John Wiley & Sons, Ltd.     

Air/water/sediment temperature contrasts in small streams to identify groundwater seepage locations

The need to identify groundwater seepage locations is of great importance for managing both stream water quality and groundwater sourced ecosystems due to their dependency on groundwater‐borne nutrients and temperatures. Although several reconnaissance methods using temperature as tracer exist, these are subjected to limitations related to mainly the spatial and temporal resolution and/or mixing of groundwater and surface water leading to dilution of the temperature differences. Further, some methods, for example, thermal imagery and fiber optic distributed temperature sensing, although relative efficient in detecting temperature differences over larger distances, these are labor‐intensive and costly. Therefore, there is a need for additional cost‐effective methods identifying substantial groundwater seepage locations. We present a method expanding the linear regression of air and stream temperatures by measuring the temperatures in dual‐depth; in the stream column and at the streambed‐water interface (SWI). By doing so, we apply metrics from linear regression analysis of temperatures between air/stream and air/SWI (linear regression slope, intercept, and coefficient of determination), and the daily water temperature cycle (daily mean temperatures, temperature variance, and the mean diel temperature fluctuation). We show that using metrics from only single‐depth stream temperature measurements are insufficient to identify substantial groundwater seepage locations in a head‐water stream. Conversely, comparing the metrics from dual‐depth temperatures show significant differences; at groundwater seepage locations, temperatures at the SWI merely explain 43–75% of the variation opposed to ? 91% at the corresponding stream column temperatures. In general, at these locations at the SWI, the slopes ( < 0.25) and intercepts ( > 6.5 °C) are substantially lower and higher, respectively, while the mean diel temperature fluctuations ( < 0.98 °C) are decreased compared to remaining locations. The dual‐depth approach was applied in a post‐glacial fluvial setting, where metrics analyses overall corroborated with field measurements of groundwater fluxes and stream flow accretions. Thus, we propose a method reliably identifying groundwater seepage locations along streambeds in such settings.     

Local‐scale variability of seepage and hydraulic conductivity in a shallow gravel‐bed river

Seepage rate and direction measured with a seepage metre modified for use in flowing water were greatly variable along a 300‐m reach of a shallow, gravel‐bed river and depended primarily on the local‐scale bed topography. The median value of seepage measured at 24 locations was 24 cm/day, but seepage measured at specific sites ranged from ?340 to +237 cm/day. Seepage also varied substantially over periods of hours to days and occasionally reversed direction in response to evolution of the sediment bed. Vertical hydraulic conductivity was related to seepage direction and was larger during upward seepage than during downward seepage; with differences ranging from 4 to 40% in areas of active sediment transport to more than an order of magnitude in areas where current was too slow to mobilize bed sediment. Seepage was poorly related to hydraulic gradient measured over vertical distances of 0·3 m and appeared to be opposite the hydraulic gradient at 18% of the locations where both parameters were measured. Results demonstrate the scale dependence of these measurements in coarse‐grained hyporheic settings and indicate that hydraulic gradients should be determined over a much shorter vertical increment if used to indicate exchange across the sediment–water interface. Published in 2009 by John Wiley & Sons, Ltd.     

Sheet flow hydrodynamics over a non-uniform sand bed channel

The current study experimentally investigates the flow characteristics and temporal variations in the sheet flow profile of a non-uniform sand bed channel. Experiments were done to explore turbulent structures in the presence of a sheet flow layer with and without seepage. The turbulent events, such as stream wise velocity, Reynolds shear stresses, and turbulence intensities were found to be increasing and vertical velocity was found decreasing with a sheet layer. The presence of a sheet layer also effects the turbulent energy production and energy dissipation. All the turbulence parameters with and without a sheet layer have also been influenced by the presence of downward seepage. The rate of sheet flow movement is increased with seepage, owing to increased turbulence with seepage. The current study used wavelet analysis on temporally lagged spatial bed elevation profiles obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuation at different spatial scales. A spatial cross correlation analysis at multiple scales, based on the wavelet coefficients, has been done on these bed elevation datasets to observe the effect of downward seepage on the dynamic behavior of sheet flow at different length scales. It is found that seepage increases average bed celerity and also increases the celerity of sheet flow of similar length scales. This increase in the celerity has been hypothesized as the increase of sheet flow movement as well as the increase in turbulent parameters with seepage, which destabilizes the bed particles resulting in a disruption in the continuous propagation pattern of the sheet flow. The increase of sheet flow celerity with seepage is confirmed from the saturation level of the wavelet power spectra of the bed elevation series. The presence of seepage also affects the non-uniformity of collective sheet material.     

Thermal structure and thickness of the lithospheric mantle underlying the Siberian Craton from the kraton and kimberlit superlong seismic profiles

A self-consistent approach is proposed for the investigation of the thermal conditions, chemical composition, and internal structure of the upper mantle of the Earth. Using this approach, the thermal state of the lithospheric mantle beneath the Siberian Craton (SC) is reconstructed from P velocities, taking into account the phase transitions, anharmonicity, and the effects of anelasticity. The velocities of seismic waves are more sensitive to temperature than to the composition of the mantle rocks, which allows the velocity models to be effectively used for reconstruction of the thermal regime of the mantle. The temperature at depths 100–300 km is reconstructed by inversion of the Kraton and Kimberlit superlong seismic profiles for compositions of the garnet harzburgite, lherzolite, and intermediate composition of garnet peridotite. The averaged temperature in the normal continental mantle is reconstructed by inversion of the IASP91 reference model for depleted and fertile substance. One-dimensional models and two-dimensional thermal fields undergo a substantial fall in temperature (～300–600°C) beneath the Siberian Craton as compared to the temperatures of the continental mantle and paleotemperatures inferred from the thermobarometry of xenoliths. Temperature profiles of the Siberian Craton deduced from seismic data lie between the conductive geotherms of 32.5–40.0 mW/m² and below the P(H)-T values obtained for low- and high-temperature xenoliths from the Mir, Udachnaya, and Obnazhennaya kimberlite pipes. The thickness of the thermal lithosphere estimated from the intersection with the potential adiabat is 300–320 km, which is consistent with the data on heat flows and seismotomographic observations. This provides grounds for the assumption that the low-temperature anomalies (thermal roots of continents) penetrate down to a depth of 300 km. The analysis of the sensitivity of seismic velocity and density to the variations in temperature, pressure, and chemical and phase composition of petrological models shows that recognition of fine differences in chemical composition of the lithospheric rocks by seismic methods is impossible.     

Macroporosity and infiltration in blanket peat: the implications of tension disc infiltrometer measurements

Little is known about the processes of infiltration and water movement in the upper layers of blanket peat. A tension infiltrometer was used to measure hydraulic conductivity in a blanket peat in the North Pennines, England. Measurements were taken from the surface down to 20 cm in depth for peat under four different vegetation covers. It was found that macropore flow is a significant pathway for water in the upper layers of this soil type. It was also found that peat depth and surface vegetation cover were associated with macroporosity and saturated hydraulic conductivity. The proportion of macropore flow was found to be greater at 5 cm depth than at 0, 10 and 20 cm depth. Peat beneath a Sphagnum cover tends to be more permeable and a greater proportion of macropore flow can occur beneath this vegetation type. Functional macroporosity and matrix flow in the near‐surface layers of bare peat appear to have been affected by weathering processes. Comparision of results with rainfall records demonstrates that infiltration‐excess overland flow is unlikely to be a common runoff‐generating mechanism on blanket peat; rather, a saturation‐excess mechanism combined with percolation‐excess above much less permeable layers dominates the runoff response. Copyright © 2001 John Wiley & Sons, Ltd.     

The effect of roughness configuration on velocity profiles in an artificial channel

In order to determine the effect of bed roughness on velocity distribution, we used seven different configurations of bed roughness, with 16 test runs of varying discharge and slope for each configuration. For each run, one-dimensional velocity profiles were measured at 1 cm vertical increments over the crest of the roughness element, and at intervals of 4·25 cm downstream. Results indicate that velocity profile shape remains fairly constant for a given slope and roughness configuration as discharge increases. As slope increases, the profiles become less linear, with a much larger near-bed velocity gradient and a more pronounced velocity peak close to 0·6 flow depth at the measurement point immediately downstream from the roughness element. The zone of large near-bed velocity gradients increases in both length and depth as roughness concentration decreases, up to a length/height ratio of about 9, at which point maximum flow resistance occurs. Longitudinal roughness elements do not create nearly as much flow resistance as do transverse elements. Rates of velocity increase suggest that roughness elements spaced at a length/height ratio of about 9 are most effective at creating flow resistance over a range of discharges in channels with steeper slopes. © 1998 John Wiley & Sons, Ltd.     

A Comparison of Tools and Methods for Estimating Groundwater-Surface Water Exchange

A comparison of tools for measuring discharge rates in a sandy streambed was conducted along a transect near the north bank of the Grindsted Å (stream), Denmark. Four tools were evaluated at six locations spaced 3 m apart in the stream: mini-piezometers, streambed point velocity probes (SBPVPs), temperature profilers, and seepage meters. Comparison of the methods showed that all identified a similar trend of low to high groundwater discharges moving westward along the transect. Furthermore, it was found that the differences between discharges estimated from Darcy calculations (using the mini-pizometers), and SBPVPs were not statistically different from zero, at the 90% confidence level. Seepage meter estimates were consistently lower than those of the other two methods, but compared more reasonably with the application of a correction factor of 1.7, taken from the literature. In contrast, discharges estimated from temperature profiling (to a depth of 40 cm) were found to be about an order of magnitude less than those determined with the other methods, possibly due to interferences from horizontal hyporheic flow. Where the various methods produced statistically different discharge estimations at the same location, it is hypothesized that the differences arose from method-specific sources of bias, including installation depths. On the basis of this work, practitioners interested in measuring flow across the groundwater-surface water interface achieve the least variability with seepage meters and the SBPVP. However the accuracy of the seepage meter depended on a calibrated correction factor while that of the SBPVP did not.     

Investigation of Seepage Meter Measurements in Steady Flow and Wave Conditions

Water exchange between surface water and groundwater can modulate or generate ecologically important fluxes of solutes across the sediment‐water interface. Seepage meters can directly measure fluid flux, but mechanical resistance and surface water dynamics may lead to inaccurate measurements. Tank experiments were conducted to determine effects of mechanical resistance on measurement efficiency and occurrence of directional asymmetry that could lead to erroneous net flux measurements. Seepage meter efficiency was high (average of 93%) and consistent for inflow and outflow under steady flow conditions. Wave effects on seepage meter measurements were investigated in a wave flume. Seepage meter net flux measurements averaged 0.08 cm/h—greater than the expected net‐zero flux, but significantly less than theoretical wave‐driven unidirectional discharge or recharge. Calculations of unidirectional flux from pressure measurements (Darcy flux) and theory matched well for a ratio of wave length to water depth less than 5, but not when this ratio was greater. Both were higher than seepage meter measurements of unidirectional flux made with one‐way valves. Discharge averaged 23% greater than recharge in both seepage meter measurements and Darcy calculations of unidirectional flux. Removal of the collection bag reduced this net discharge. The presence of a seepage meter reduced the amplitude of pressure signals at the bed and resulted in a nearly uniform pressure distribution beneath the seepage meter. These results show that seepage meters may provide accurate measurements of both discharge and recharge under steady flow conditions and illustrate the potential measurement errors associated with dynamic wave environments.     

Base flow‐driven shifts in tropical stream temperature regimes across a mean annual rainfall gradient

Climate change is expected to affect air temperature and watershed hydrology, but the degree to which these concurrent changes affect stream temperature is not well documented in the tropics. How stream temperature varies over time under changing hydrologic conditions is difficult to isolate from seasonal changes in air temperature. Groundwater and bank storage contributions to stream flow (i.e., base flow [BF]) buffer water temperatures against seasonal and daily fluctuations in solar radiation and air temperature, whereas rainfall‐driven runoff produces flooding events that also influence stream temperature. We used a space‐for‐time substitution to examine how shifts in BF and runoff alter thermal regimes in streams by analyzing hydrological and temperature data collected from similar elevations (400–510 m above sea level) across a 3,500‐mm mean annual rainfall gradient on Hawai'i Island. Sub‐daily water temperature and stream flow gathered for 3 years were analyzed for daily, monthly, and seasonal trends and compared with air temperature measured at multiple elevations. Results indicate that decreases in median BF increased mean, maximum, and minimum water temperatures as well as daily temperature range. Monthly and daily trends in stream temperature among watersheds were more pronounced than air temperature, driven by differences in groundwater inputs and runoff. Stream temperature was strongly negatively correlated to BF during the dry season but not during the wet season due to frequent wet season runoff events contributing to total flow. In addition to projected increases in global air temperature, climate driven shifts in rainfall and runoff are likely to affect stream flow and groundwater recharge, with concurrent influences on BF resulting in shifts in water temperature that are likely to affect aquatic ecosystems.     

Models of stream bed topography at the reach scale

A simple model of channel geometry is presented which focuses on the way that a single form variable may be autocorrelated with upstream values of that variable and which is amenable to analysis using time series methods. Stream bed topography at the local scale is the variable of interest. Bed height series were obtained for four stream lengths located at different positions along a mountain stream to test the hypothesis that, as control variables change spatially, the form of the model will also change. The Box-Jenkins models describing the series are dominated by autoregressive terms and indicate a non-random variation of bed height. An increase in the order and parameter magnitude of the models is associated with an increase in discharge and a decrease in bed material size, suggesting that the structure of the model reflects the adjustability of the stream bed and, specifically, the extent of riffle-pool development, although channel gradients may be too high for a well-defined riffle-pool sequence to be developed. Changes in the structure of models obtained for bed height and other form variables may indicate the variable ability of a stream to adjust its channel morphology.     

Effect of seepage on flow and bedforms dynamics

This study, using an experimental approach, focuses on the effect of downward seepage on a threshold alluvial channel morphology and corresponding turbulent flow characteristics. In all the experiments, we observed that the streamwise time‐averaged velocities and Reynolds shear stresses were increased under the influence of downward seepage. Scales of eddy length and eddy turnover time were significantly increased with the application of downward seepage, leading to sediment transport and initiation of bedforms along the channel length. As the amount of seepage discharge increased, eddy length and turnover time were further increased, causing the development of larger bedforms. It was revealed that the geometry of bedforms was linked with the size of eddies. In this work, statistics of bedform dynamics are presented in terms of multi‐scalar bedforms in the presence of seepage. These multi‐scalar ubiquitous bedforms cast a potential impact on flow turbulence as well as stream bed morphology in channels. We used wavelet to analyse temporally lagged spatial bed elevation profiles that were obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuations at different length scales. A spatial cross‐correlation analysis, based on the wavelet coefficients, was performed on these bed elevation datasets to observe the effect of downward seepage on the dynamic behaviour of bedforms at different length scales. It was found that celerity of bedforms reduced with increase in seepage percentage. Bedform celerity was best approximated by a probability density function such as Rayleigh distribution under varying downward seepage. Further, statistical analysis of physical parameters of bedforms ascertained that the reduction in bedform celerity was a result of increased bedform size. Copyright © 2017 John Wiley & Sons, Ltd.     

Two-dimensional mantle flow beneath a rigid accreting lithosphere

This study considers two-dimensional mantle flow beneath a rigid lithosphere. The lithosphere which forms the upper boundary of a convecting region moves with a prescribed uniform horizontal velocity, and thickens with distance from the accreting plate boundary as it cools. Beneath the lithosphere, the mantle deforms viscously by diffusion creep and is heated radiogenically from within. Solutions for thermal convection beneath the lithosphere are obtained by finite-difference methods. Two important conclusions have resulted from this study: (1) convective patterns of large aspect ratio are stable beneath a rigid moving lithosphere; (2) even for a lithosphere velocity as small as 3 cm/yr. and a Rayleigh number as large as 10⁶, mantle circulation with large aspect ratio is driven dominantly by the motion of the lithosphere rather than by temperature gradients within the flow. Gravity, topography and heat flow are determined and implications for convection in the upper mantle are discussed.     

Scour Downstream of Grade Control Structures under the Influence of Upward Seepage

The installation of free falling jet grade control structures has become a popular choice for river bed stabilization. However, the formation and development of scour downstream of the structure may lead to failure of the structure itself. The current approaches to scour depth prediction are generally based on studies conducted with the absence of upward seepage. In the present study, the effects of upward seepage on the scour depth were investigated. A total of 78 tests without and with the application of upward seepage were carried out using three different sediment sizes, three different tailwater depths, four different flow discharges, and four different upward seepage flow discharge rates. In some tests, the three-dimensional components of the flow velocity within the scour hole were measured for both the cases with and without upward seepage. The scour depth measured for the no-seepage results compared well with the most accurate relationship found in the literature. It was found that generally the upward seepage reduced the downward velocity components near the bed, which led to a decrease in the maximum scour depth. A maximum scour depth reduction of 49% was found for a minimum tailwater depth, small sediment size, and high flow discharge. A decay of the downward velocity vector within the jet impingement was found due to the upward seepage flow velocity. The well known equation of D’Agostino and Ferro was modified to account for the effect of upward seepage, which satisfactorily predicted the experimental scour depth, with a reasonable average error of 10.7%.     

Fine sediment dynamics in a shallow lake and implication for design of hydraulic works

Lake Markermeer is a large (680?km²), shallow body of water in the middle of the Netherlands, with a mean water depth of 3.6?m. One of the major problems in the lake is its decreasing ecological value which is, among other reasons, caused by a gradual increase of suspended sediment concentration and associated increase of light attenuation in the water column. A thorough understanding of fine sediment dynamics in the lake is a prerequisite for solving this problem. This paper addresses the 3D nature of near-bed sediment dynamics in Lake Markermeer, based on data sampled from a 1-month field experiment in autumn 2007. The campaign involved the collection of 71 bed samples across the lake. At each location, dual-frequency echo soundings were carried out to assess the thickness of the silt layer, and sediment concentration throughout the water column was measured with an Optical Backscatter Sensor (OBS). Moreover, 2-week time series of wave height, water level, current velocities, and near-bed sediment concentration were collected at a single location. The time series of sediment concentration were measured with a regular OBS and an Argus Surface Meter IV (ASM). During the measurement period, flow velocities ranged between 2 and 15?cm/s, wave heights up to 1.2?m were observed and turbidity levels varied between 40?mg/l to more than 300?mg/l. The ASM data generally showed uniform concentration profiles. However, profiles with steep concentration gradients near the bed were found for wave heights above 0.5?m. The field experiments further revealed pronounced 3D structures near the bed during discrete storms. The results are generalized for a wider range of conditions and across the full water depth through application of a 1DV point model, using a two-fraction representation of the grain size distribution. The fine and coarse fractions are found to resuspend rapidly for wind speeds above 5?m/s and 10??2?m/s, respectively, forming a uniform concentration profile if these wind conditions persists. High-concentration (???g/l) layers near the bed, containing the coarse sediment fraction, only occur at the onset and towards the end of a storm, when wind speed changes rapidly. It is under these conditions that horizontal gradients in layer density or thickness can transport considerable fine sediment. This transport provides an additional mechanism for the infill of, for instance, silt traps and navigation channels.     

塔里木盆地岩石层热结构特征

在大地热流密度分布的基础上，研究了塔里木盆地中库尔勒－若羌和阿克苏－叶城两条剖面岩石层热结构特征．由岩石层Ｐ波速度分布转换成生热率剖面，用二维数值模型获得了岩石层热结构和热状态特征．结果表明，塔里木盆地壳幔边界温度的高低与其埋深密切相关．居里等温面深度大，地幔热流密度较低．岩石层厚度变化与其新生代期间挠曲过程密切相关．在岩石层温度分布基础上，确定了深部脆－韧性过渡带深度和岩石层屈服强度，表明塔里木盆地岩石层相对较冷，且具有刚性的地球动力学特征．     

Influence of a thin veneer of low‐hydraulic‐conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells. Published in 2011 by John Wiley & Sons, Ltd.     

Response of soil detachment rate to the hydraulic parameters of concentrated flow on steep loessial slopes on the Loess Plateau of China

Using hydraulic parameters is essential for describing soil detachment and developing physically based erosion prediction models. Many hydraulic parameters have been used, but the one that performs the best for describing soil detachment on steep slopes when the lateral expansion (widening) of rills is not limited has not been identified. An indoor concentrated flow scouring experiment was performed on steep loessial slopes to investigate soil detachment rates for different flow rates and slope gradients. The experiments were conducted on a slope‐adjustable plot (5 m length, 1 m width, 0.5 m depth). Sixteen combinations of 4 flow rates (10, 15, 20, and 25 L/min) and 4 slope gradients (17.6%, 26.8%, 36.4%, and 46.6%) were investigated. The individual and combined effects of slope gradient and flow hydraulic parameters on soil detachment rate were analysed. The results indicated that soil detachment rate increased with flow rate and slope gradient. Soil detachment rate varied linearly and exponentially with flow rate and slope gradient, respectively. Multivariate, nonlinear regression analysis indicated that flow depth exerted the greatest influence on the soil detachment rate, followed by unit discharge per unit width, slope gradient, and flow rate in this study. Shear stress and stream power could efficiently describe the soil detachment rate using a power equation. However, the unit stream power and unit energy of the water‐carrying section changed linearly with soil detachment rate. Stream power was an optimal hydraulic parameter for describing soil detachment. These findings improve our understanding of concentrated flow erosion on steep loessial slopes.     

Thermal modeling of the Southern Alps,New Zealand

Finite-element modeling of the thermal regime across the Southern Alps of New Zealand has been carried out along two profiles situated near the Franz Josef and Haast valleys. The modeling involves viscous deformation beneath the Southern Alps, including both uplift and erosion, and crustal/lithospheric thickening, as a result of crustal shortening extending to 20 mm/y of a 25-km thick crust. Published uplift rates and crustal thickness variations along the two profiles are used to constrain the modeled advection of crustal material, and results are compared with the recent heat flow determinations, 190±50 mW/m² in the Franz Josef valley and 90±25 mW/m² in the Haast valley. Comparisons of the model with published K–Ar and fission track ages, show that the observed heat flow in the Franz Josef valley is consistent with observed zircon fission track ages of around 1 Ma, if the present-day uplift rate is close to 10 mm/y. Major thermal differences between the Franz Josef and Haast profiles appear to be due to different uplift and erosion rates. There is weak evidence that frictional heating close to the Alpine fault zone is not significant. The modeling provides explanations for the distribution of seismicity beneath the Southern Alps, and predicts a low surface heat flow over the eastern foothills due to the dominant thermal effect of crustal thickening beneath this region. Predicted temperatures at mid-crustal depth beneath the zone of maximum uplift rate are 50–100°C cooler than those indicated in previously published models, which implies that thermal weakening of the crust may not be the main factor causing the aseismicity of the central Southern Alps. The results of the modeling demonstrate that the different types of reset age data in the region within 25 km of the Alpine fault are critical for constraining models of the deformation and the thermal regime beneath the Southern Alps.