In-Well Degassing of Monitoring Wells Completed in Gas-Charged Aquifers

Total dissolved gas pressure (P_TDG) measurements are useful to measure accurate in situ dissolved gas concentrations in groundwater, but challenged by in-well degassing. Although in-well degassing has been widely observed, its cause(s) are not clear. We investigated the mechanism(s) by which gas-charged groundwater in a recently pumped well becomes degassed. Vertical P_TDG and dissolved gas concentration profiles were monitored in the standing water column (SWC) of a groundwater well screened in a gas-charged aquifer for 7 days before and 15 days after pumping. Prior to pumping, P_TDG values remained relatively constant and below calculated bubbling pressure (P_BUB) at all depths. In contrast, significant increases in P_TDG were observed at all depths after pumping was initiated, as fresh groundwater with elevated in situ P_TDG values was pumped through the well screen. After pumping ceased, P_TDG values decreased to below P_BUB at all depths over the 15-day post-pumping period, indicating well degassing was active over this time frame. Vertical profiles of estimated dissolved gas concentrations before and after pumping provided insight into the mechanism(s) by which in-well degassing occurred in the SWC. During both monitoring periods, downward mixing of dominant atmospheric and/or tracer gases, and upwards mixing of dominant groundwater gases were observed in the SWC. The key mechanisms responsible for in-well degassing were (i) bubble exsolution when P_TDG exceeded P_BUB as gas-charged well water moves upwards in the SWC during recovery (i.e., hydraulic gradient driven convection), (ii) microadvection caused by the upward migration of bubbles under buoyancy, and (iii) long-term, thermally driven vertical convection.     

Total Dissolved Gas Pressure Measurements for Assessing Contaminant Degradation Processes in Groundwater

The total dissolved gas pressure (P_TDG ) probe has been used in groundwater studies for over a decade, but rarely in assessing contaminant degradation, despite the many degradation reactions that produce or consume dissolved gases. Here we present three studies to demonstrate the application of P_TDG measurements to groundwater experiencing contaminant degradation, with discussion of its benefits and limitations. The first study is a pilot‐scale laboratory experiment simulating dissolved ethanol contamination of an anaerobic sand aquifer. Continuous monitoring of P_TDG showed the rapid onset of microbial hydrocarbon degradation via denitrification and fermentation. The subsequent formation of a gas phase was revealed when P_TDG began mimicking the bubbling pressure (P_G *; sum of hydrostatic and atmospheric pressure), fluctuating with atmospheric pressure. Some deviations of P_TDG above P_G * occurred also, which may hold promise for signalling substantial changes in the rate or type of degradation process (here, the onset of methanogenesis). In the second study, synoptic field measurements at a petroleum plume site demonstrated how elevated P_TDG could identify wells with evidence of hydrocarbon degradation (denitrification and/or methanogenesis). And finally, combined field measurements of dissolved oxygen (DO) and P_TDG in monitoring wells of a nitrate‐contaminated aquifer (Abbottsford‐Sumas) revealed areas where denitrification was likely occurring. Limitations to P_TDG use identified in these studies included the masking of degradation processes by the presence of a gas phase, as when trapped following water table fluctuations or formed from rigorous degradation reactions, and confounded assessment of P_TDG patterns from other natural or anthropogenic processes that can also influence groundwater P_TDG .     

Changes in Entrapped Gas Content and Hydraulic Conductivity with Pressure

Water table fluctuations continuously introduce entrapped air bubbles into the otherwise saturated capillary fringe and groundwater zone, which reduces the effective (quasi‐saturated) hydraulic conductivity, K_quasi, thus impacting groundwater flow, aquifer recharge and solute and contaminant transport. These entrapped gases will be susceptible to compression or expansion with changes in water pressure, as would be expected with water table (and barometric pressure) fluctuations. Here we undertake laboratory experiments using sand‐packed columns to quantify the effect of water table changes of up to 250 cm on the entrapped gas content and the quasi‐saturated hydraulic conductivity, and discuss our ability to account for these mechanisms in ground water models. Initial entrapped air contents ranged between 0.080 and 0.158, with a corresponding K_quasi ranging between 2 and 6 times lower compared to the K_s value. The application of 250 cm of water pressure caused an 18% to 26% reduction in the entrapped air content, resulting in an increase in K_quasi by 1.16 to 1.57 times compared to its initial (0 cm water pressure) value. The change in entrapped air content measured at pressure step intervals of 50 cm, was essentially linear, and could be modeled according to the ideal gas law. Meanwhile, the changes in K_quasi with compression–expansion of the bubbles because of pressure changes could be adequately captured with several current hydraulic conductivity models.     

Evidence that piezometers vent gas from peat soils and implications for pore‐water pressure and hydraulic conductivity measurements

Entrapped gas bubbles in peat can alter the buoyancy, storativity, void ratio and expansion/contraction properties of the peat. Moreover, when gas bubbles block water‐conducting pores they can significantly reduce saturated hydraulic conductivity and create zones of over‐pressuring, perhaps leading to an alteration in the magnitude and direction of groundwater flow and solute transport. Some previous researches have demonstrated that these zones of over‐pressuring are not observed by the measurements of pore‐water pressures using open‐pipe piezometers in peat; rather, they are only observed with pressure transducers sealed in the peat. In has been hypothesized that open‐pipe piezometers vent entrapped CH₄ to the atmosphere and thereby do not permit the natural development of zones of entrapped gas. Here we present findings of the study to investigate whether piezometers vent subsurface CH₄ to the atmosphere and whether the presence of piezometers alters the subsurface concentration of dissolved CH₄. We measured the flux of methane venting from the piezometers and also determined changes in pore‐water CH₄ concentration at a rich fen in southern Ontario and a poor fen in southern Quebec, in the summer of 2004. Seasonally averaged CH₄ flux from piezometers was 1450 and 37·8‐mg CH₄ m^?2 d^?1 at the southern Ontario site and Quebec site, respectively. The flux at the Ontario site was two orders of magnitude greater than the diffusive flux at the site. CH₄ pore‐water concentrations were significantly lower in open piezometers than in water taken from sealed samplers at both the Ontario and Quebec sites. The flux of CH₄ from piezometers decreased throughout the season suggesting that CH₄ venting through the piezometer exceeded the rate of methanogenesis in the peat. Consequently we conclude that piezometers may alter the gas dynamics of some peatlands. We suggest that less‐invasive techniques (e.g. buried pressure transducers, tracer experiments) are needed for the accurate measurement of pore‐water pressures and hydraulic conductivity in peatlands with a large entrapped gas component. Furthermore, we argue that caution must be made in interpreting results from previous peatland hydrology studies that use these traditional methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

AMSR‐E algorithm for snowmelt onset detection in sub‐arctic heterogeneous terrain

The onset of snowmelt in the upper Yukon River basin, Canada, can be derived from brightness temperatures (T_b) obtained by the Advanced Microwave Scanning Radiometer for EOS (AMSR‐E) on NASA's Aqua satellite. This sensor, with a resolution of 14 × 8 km² for the 36·5 GHz frequency, and two to four observations per day, improves upon the twice‐daily coverage and 37 × 28 km² spatial resolution of the Special Sensor Microwave Imager (SSM/I). The onset of melt within a snowpack causes an increase in the average daily 36·5 GHz vertically polarized T_b as well as a shift to high diurnal amplitude variations (DAV) as the snow melts during the day and re‐freezes at night. The higher temporal and spatial resolution makes AMSR‐E more sensitive to sub‐daily T_b oscillations, resulting in DAV that often show a greater daily range compared to SSM/I. Therefore, thresholds of T_b > 246 K and DAV > ± 10 K developed for use with SSM/I have been adjusted for detecting the onset of snowmelt with AMSR‐E using ground‐based surface temperature and snowpack wetness relationships. Using newly developed thresholds of T_b > 252 K and DAV > ± 18 K, AMSR‐E derived snowmelt onset correlates well with SSM/I observations in the small subarctic Wheaton River basin through the 2004 and 2005 winter/spring transition. In addition, the onset of snowmelt derived from AMSR‐E data gridded at a higher resolution than the SSM/I data indicates that finer‐scale differences in elevation and land cover affect the onset of snowmelt and are detectable with the AMSR‐E sensor. On the basis of these observations, the enhanced resolution of AMSR‐E is more effective than SSM/I at delineating spatial and temporal snowmelt dynamics in the heterogeneous terrain of the upper Yukon River basin. Copyright © 2007 John Wiley & Sons, Ltd.     

Flow resistance in gravel‐bed channels with large‐scale roughness

A previously published mixing length (ML) model for evaluating the Darcy–Weisbach friction factor for a large‐scale roughness condition (depth to sediment height ratio ranging from 1 to 4) is brie?y reviewed and modi?ed (MML). Then the MML model and a modi?ed drag (MD) model are experimentally tested using laboratory measurements carried out for gravel‐bed channels and large‐scale roughness condition. This analysis showed that the MML gives accurate estimates of the Darcy–Weisbach coef?cient and for Froude number values greater than 0·5 the MML model coincides with the ML one. Testing of the MD model shows limited accuracy in estimating ?ow resistance. Finally, the MML and MD models are compared with the performance of a quasi‐theoretical (QT) model deduced applying the P‐theorem of the dimensional analysis and the incomplete self‐similarity condition for the depth/sediment ratio and the Froude number. Using the experimental gravel‐bed data to calibrate the QT model, a constant value of the exponent of the Froude number is determined while two relationships are proposed for estimating the scale factor and the exponent of the depth/sediment ratio. This indirect estimate procedure of the coef?cients (b₀, b₁ and b₂) of the QT model can produce a negligible overestimation or underestimation of the friction factor. Copyright © 2003 John Wiley & Sons, Ltd.     

Laboratory measurements of ‘porosity‐free’ intrinsic Vp and Vs in an olivine gabbro of the Oman ophiolite: Implication for interpretation of the seismic structure of lower oceanic crust

In order to determine ‘porosity‐free’ intrinsic ultrasonic compressional (Vp) and shear wave (Vs) velocities and Vp/Vs of an olivine gabbro from the Oman ophiolite, we developed a new experimental system using a piston‐cylinder type high‐pressure apparatus. The new system allowed us to measure velocities at pressures ranging from 0.20 to 1.00 GPa and at temperatures up to 300°C for Vp and 400°C for Vs. At room temperature, the Vp and Vp/Vs increase rapidly with pressure up to 0.40 GPa, while between 0.45 and 1.00 GPa the increase is more gradual. The change in increasing rate is attributed to closure of porosity at pressures above 0.45 GPa. Based on the linear regression of data obtained at higher pressures (0.45–1.00 GPa) and extrapolation to the lower pressures, combined with temperature derivatives of velocities of the sample measured at 1.00 GPa, we determined the intrinsic Vp and Vs of the sample as a function of pressure (P, in GPa) and temperature (T, in °C). The intrinsic velocities can be expressed as Vp (km/s) = 7.004 + 0.096 × P ? 0.00015 × T, and Vs (km/s) = 3.827 + 0.007 × P ? 0.00008 × T. We evaluated the intrinsic Vp and Vs of the olivine gabbro at oceanic crustal conditions and compared them with a velocity depth‐profile of the borehole seismic observatory WP‐2 area in the northwestern Pacific Basin. Although the intrinsic Vp (～7.0 km/s) and Vs (～3.8 km/s) for the olivine gabbro studied are comparable to those of seismic layer 3 in the WP‐2 area, the estimated vertical gradients of intrinsic velocities are significantly smaller than those reported from layer 3. These results suggest that velocity profiles of layer 3 in the WP‐2 area may reflect the presence of a minor porosity in lower oceanic crust, which closes with increasing depth and/or continuous changes in mineralogy of layer 3 rocks.     

Diurnal and seasonal variations of sap flow of Caragana korshinskii in the arid desert region of north‐west China

The aim of this study was to obtain the diurnal and seasonal changes of trunk sap flow in desert‐living Caragana korshinskii so as to understand its water requirement and ecological significance. The experiment was carried out with 15‐year old Caragana korshinskii grown in north‐west China under natural conditions. Heat pulse sensors based on the heat compensation theory were applied to measure the trunk sap flow, and soil moisture content at 0–300 cm layer, using tube‐type time domain reflectometry (Tube‐TDR). The solar radiation, the maximum and minimum air temperatures, relative humidity, wind speed, wind direction and precipitation were measured at a standard automatic weather station. The diurnal and seasonal variations of sap flow rate, the sap velocity at different positions in the trunk and the sap flow rate under different weather conditions were analysed. And the correlation between the sap flow rate and the meteorological factors was also analysed. Results showed that the trunk sap flow varied regularly in the diurnal term and the sap flow velocity decreased with the probe‐inserted depth into the sapwood. Magnitude of sap flow changed considerably between sunny and rainy days. The order of the main meteorological factors affecting the sap flow rate of Caragana korshinskii shrubs were: vapour pressure deficit > solar radiation > air temperature > wind speed. The close correlation between daily sap flow rate and meteorological factors in the whole growing season can be used to estimate the transpiration of Caragana korshinskii. Copyright © 2007 John Wiley & Sons, Ltd.     

Multi‐decadal water table manipulation alters peatland hydraulic structure and moisture retention

A peatland complex disturbed by berm construction in the 1950s was used to examine the long‐term impact of water table (WT) manipulation on peatland hydraulic properties and moisture retention at three adjacent sites with increasing depth to WT (WET, INTermediate reference and DRY). Saturated hydraulic conductivity (K_s) was found to decrease with depth by several orders of magnitude over a depth of 1–1.5 m at all sites. The depth dependence of WT response to rainfall was similar across sites: WT response increased from 1 : 1 at the surface, to 5 : 1 at 50 cm depth. While surface specific yield (S_y) values were similar across all sites, it decreased with depth at a rate of 0.014 cm^?1 in hollows and 0.007 cm^?1 in hummocks. Bulk density (ρ_b) exhibited similar depth‐dependent trends as S_y and explains a high amount of variance (r² > 0.69) in moisture retention across a range of pore water pressures (?15 to ?500 cm H₂O). Because of higher ρ_b, hollow peat had greater moisture retention, where site effects were minimal. However, the estimated residual water content for surface Sphagnum samples, while on average lower in hummocks (0.082 m³ m^?3) versus hollows (0.087 m³ m^?3), increased from WET (0.058 m³ m^?3) to INT (0.088 m³ m^?3) to DRY (0.108 m³ m^?3) which has important implications for moisture stress under conditions of persistent WT drawdown. Given the potential importance of microtopographic succession for altering peatland hydraulic structure, our findings point to the need for a better understanding of what controls the relative height and proportional coverage of hummocks in relation to long‐term disturbance‐response dynamics. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimation of tree water stress from stem and soil water monitoring with time‐domain reflectometry in two small forested basins in Spain

Soil‐tree water relationships were studied using time domain reflectometry (TDR) in two small forested basins in Spain. The stem water content of two Mediterranean Quercus species (Quercus pyrenaica and Quercus rotundifolia) was measured using previously constructed species‐specific equations. To monitor soil moisture, a TDR station network was used in both cases. Sixteen Q. pyrenaica and six Q. rotundifolia individuals were selected to install two TDR probes in their trunks (at 20 and 120 cm above the ground) to monitor stem water content. Stem and soil water contents were measured fortnightly. The stem water content of both species showed a similar temporal trend for the period studied. A spring maximum (0·654 cm³ cm^?3 for Q. pyrenaica and 0·568 cm³ cm^?3 for Q. rotundifolia) was found to be associated with high transpiration and no soil moisture deficit, and a late‐summer minimum (0·520 cm³ cm^?3 for Q. pyrenaica and 0·426 cm³ cm^?3 for Q. rotundifolia) was associated with the end of the dry season. This drop in stem water content occurs when the available water in the soil decreases. This seasonal difference presumably reflects water withdrawn from stem storage to support the transpirational demands of the tree. Since plant water stress results in reduced stem water content and since this drop can be measured by TDR, it may be concluded that this technology offers a suitable tool for detecting plant water stress. Copyright © 2007 John Wiley & Sons, Ltd.     

Determining thermal diffusivity using near‐surface periodic temperature variations and its implications for tracing groundwater movement at the eastern margin of the Tibetan Plateau

The thermal diffusivity is the key parameter that controls near‐surface temperature where periodic temperature variation is progressively attenuated and delayed with depth. This article presents the results of apparent thermal diffusivity using temperatures recorded by a bedrock temperature measurement network in the fault zones of western Sichuan. High sensitivity temperature sensors (10^?4 K) were installed at a maximum depth reaching 30 m. The apparent thermal diffusivities were deduced from both amplitude damping and phase shifting of annual temperature variations between two different depths. Under pure conduction, the thermal diffusivity determined through the phase method (α_Φ) should be equivalent to that determined through the amplitude method (α_A), whereas effects of the upward (downward) water flow are evidently reflected in the amplitude decay to make α_Φ larger (lesser) than α_A. The discrepancy between α_Φ and α_A can thus be a tracer of water movement or convective heat transfer. The calculated α_Φ of the measurement stations varies from 1.22 × 10^?6 to 3.00 × 10^?6 m²/s, and the estimated α_A ranges from 0.93 × 10^?6 to 2.41 × 10^?6 m²/s. Two regimes of heat transfer underground were suggested from the results. Conductive heat transport prevails over the nonconductive processes at five stations, which is characterized by α_Φ coincident with α_A for the same depth pair. On the contrary, the values of α_Φ differ from α_A at six stations in the intersection area of the Y‐shaped fault system, implying that convective heat transfer also plays a comparably important role. This finding is consistent with the hot springs distribution of the area. The results also indicate that water moves upward with an average Darcy velocity of approximately ?1 × 10^?7 m/s in this region. Our research provides new evidence for the hydrothermal activity in the fault zones at the eastern margin of the Tibetan Plateau.     

Application of an S‐shaped curve model to the temporal development of tafoni of salt‐weathering origin

To describe temporal change in tafone development, an S‐shaped curve equation is proposed: Z = Z_c [1 ? (n + 1) exp (? β t ) + n exp (? (1 + 1/n) β t )] , where Z is observed tafone depth, Z_c is ultimate tafone depth, t is time, and n and β are constants. The applicability of this model is examined using tafone data selected from seven sites, which are categorized into three different salt‐weathering environments: a spray/splash‐dominant (occasionally wave‐affected) supra‐tidal zone, aerosol‐affected coastal regions, and inland desert areas. The results indicate that the equation can well describe tafone development in each of these environments. An investigation based on the values of n and β, determined through a best fit of the equation to the data, suggests that n characterizes site‐specific environmental conditions and β reflects the magnitude of factors controlling the recession mechanism of tafone surfaces. It is found that (1) the maximum rate of tafone growth dramatically decreases from supra‐tidal, through coastal, to desert environments, and (2) the growing mode of tafoni is different depending on the environmental settings. The erosional force to facilitate the development of tafoni at supra‐tidal sites is estimated to be about 400 times greater than that in the general coastal area. Copyright © 2011 John Wiley & Sons, Ltd.     

Transfer of carbon dioxide and methane through the soil‐water‐atmosphere system at Mer Bleue peatland,Canada

Surface waters associated with peatlands, supersaturated with CO₂ and CH₄ with respect to the atmosphere, act as important pathways linking a large and potentially unstable global repository of C to the atmosphere. Understanding the drivers and mechanisms which control C release from peatland systems to the atmosphere will contribute to better management and modelling of terrestrial C pools. We used non‐dispersive infra‐red (NDIR) CO₂ sensors to continuously measure gas concentrations in a beaver pond at Mer Bleue peatland (Canada); measurements were made between July and August 2007. Concentrations of CO₂ in the surface water (10 cm) reached 13 mg C l^?1 (epCO₂ 72), and 26 mg C l^?1 (epCO₂ 133) at depth (60 cm). The study also showed large diurnal fluctuations in dissolved CO₂ which ranged in amplitude from ～1·6 mg C l^?1 at 10 cm to ～0·2 mg C l^?1 at 60 cm depth. CH₄ concentration and supersaturation (epCH₄) measured using headspace analysis averaged 1·47 mg C l^?1 and 3252, respectively; diurnal cycling was also evident in CH₄ concentrations. Mean estimated evasion rates of CO₂ and CH₄ over the summer period were 44·92 ± 7·86 and 0·44 ± 0·25 µg C m^?2s^?1, respectively. Open water at Mer Bleue is a significant summer hotspot for greenhouse gas emissions within the catchment. Our results suggest that CO₂ concentrations during the summer in beaver ponds at Mer Bleue are strongly influenced by biological processes within the water column involving aquatic plants and algae (in situ photosynthesis and respiration). In terms of carbon cycling, soil‐stream connectivity at this time of year is therefore relatively weak. Copyright © 2008 John Wiley & Sons, Ltd.     

Atmospheric Concentrations and Phase Partitioning of Polycyclic Aromatic Hydrocarbons in Izmir,Turkey

Ambient air polycyclic aromatic hydrocarbon (PAH) samples were collected at a suburban (n = 63) and at an urban site (n = 14) in Izmir, Turkey. Average gas‐phase total PAH (∑₁₄PAH) concentrations were 23.5 ng m^?3 for suburban and 109.7 ng m^?3 for urban sites while average particle‐phase total PAH concentrations were 12.3 and 34.5 ng m^?3 for suburban and urban sites, respectively. Higher ambient PAH concentrations were measured in the gas‐phase and ∑₁₄PAH concentrations were dominated by lower molecular weight PAHs. Multiple linear regression analysis indicated that the meteorological parameters were effective on ambient PAH concentrations. Emission sources of particle‐phase PAHs were investigated using a diagnostic plot of fluorene (FLN)/(fluorine + pyrene; PY) versus indeno[1,2,3‐cd]PY/(indeno[1,2,3‐cd]PY + benzo[g,h,i]perylene) and several diagnostic ratios. These approaches have indicated that traffic emissions (petroleum combustion) were the dominant PAH sources at both sites for summer and winter seasons. Experimental gas–particle partition coefficients (K_P) were compared to the predictions of octanol–air (K_OA) and soot–air (K_SA) partition coefficient models. The correlations between experimental and modeled K_P values were significant (r² = 0.79 and 0.94 for suburban and urban sites, respectively, p < 0.01). Octanol‐based absorptive partitioning model predicted lower partition coefficients especially for relatively volatile PAHs. However, overall there was a relatively good agreement between the measured K_P and soot‐based model predictions.     

Effect of SO on 1,1,1‐Trichloroethane Degradation by Fe0 in Aqueous Solution

Sulfate in groundwater has been previously shown to change the reactivity of Fe⁰ in permeable reactive barriers for reducing chlorinated organics. To better understand the effect and mechanism of SO, the degradation of 1,1,1‐trichloroethane (TCA) by Fe⁰ in unbuffered aqueous solutions with and without SO was investigated. In a Fe⁰‐TCA‐H₂O system with initial pH of 2.0 to 10.0, the maximum removal rate of TCA was achieved at the initial pH 6.0 with pseudo‐first‐order constant K_obs 9.0 × 10^?3/min. But in a Fe⁰‐TCA‐Na₂SO₄‐H₂O system, the removal rate of TCA decreased remarkably with a reduction in K_obs to 1.0 × 10^?3/min, and the pH varied from 6.0 to 9.6, indicating an inhibition of TCA dehydrochlorination by SO. Sulfate remarkably inhibited TCA degradation via changing the route of Fe⁰ dissolution. It accelerated the dissolution of Fe⁰ and transformed the intermediate form Fe(OH)_ads to Fe₂(SO₄)_ads, which weakened the affinity between Fe and TCA, and thus depressed the degradation of TCA by Fe⁰.     

Magnetic enhancement in wildfire‐affected soil and its potential for sediment‐source ascription

Intense rainfall following wildfire can cause substantial soil and sediment redistribution. With concern for the increasing magnitude and frequency of wildfire events, research needs to focus on hydrogeomorphological impacts of fire, particularly downstream fluxes of sediment and nutrients. Here, we investigate variation in magnetic enhancement of soil by fire in burnt eucalypt forest slopes to explore its potential as a post‐fire sediment tracer. Low‐frequency magnetic susceptibility values (χ_lf) of <10 µm material sourced from burnt slopes (c. 8·0–10·4 × 10^?6 m³ kg^?1) are an order of magnitude greater than those of <10 µm material derived from long‐unburnt areas (0·8 × 10^?6 m³ kg^?1). Susceptibility of anhysteretic remanent magnetization (χARM) and saturation isothermal remanent magnetization (SIRM) values are similarly enhanced. Signatures are strongly influenced by soil and sediment particle size and storage of previously burnt material in footslope areas. Whilst observations indicate that signatures based on magnetic enhancement show promise for post‐fire sediment tracing, problems arise with the lack of dimensionality in such data. Magnetic grain size indicators χ_fd%, χARM/SIRM and χ_fd/χARM offer further discrimination of source material but cannot be included in numerical unmixing models owing to non‐linear additivity. This leads to complications in quantitatively ascribing downstream sediment to source areas of contrasting burn severity since sources represent numerical multiples of each other, indicating the need to involve additional indicators, such as geochemical evidence, to allow a more robust discrimination. Copyright © 2005 John Wiley & Sons, Ltd.     

Time‐lapse electrical resistivity imaging of solute transport in a karst conduit

The use of electrical resistivity surveys to locate karst conduits has shown mixed success. However, time‐lapse electrical resistivity imaging combined with salt injection improves conduit detection and can yield valuable insight into solute transport behaviour. We present a proof of concept above a known karst conduit in the Kentucky Horse Park (Lexington, Kentucky). A salt tracer solution was injected into a karst window over a 45‐min interval, and repeat resistivity surveys were collected every 20 min along a 125‐m transect near a monitoring well approximately 750 m downgradient from the injection site. In situ fluid conductivity measurements in the well peaked at approximately 25% of the initial value about 3 h after salt injection. Time‐lapse electrical resistivity inversions show two broad zones at the approximate conduit depth where resistivity decreased and then recovered in general agreement with in situ measurements. Combined salt injection and electrical resistivity imaging are a promising tool for locating karst conduits. The method is also useful for gaining insight into conduit geometry and could be expanded to include multiple electrical resistivity transects. Copyright © 2015 John Wiley & Sons, Ltd.     

Accuracy of kinematic wave and diffusion wave for spatial‐varying rainfall excess over a plane

The kinematic‐wave and diffusive‐wave approximations were investigated for unsteady overland flow resulting from spatially varying rainfall excess. Three types of boundary conditions were adopted: zero flow at the upstream end, and critical flow and zero depth‐gradient at the downstream end. Errors were derived by comparing the dimensionless profiles of the flow depth over the plane with those computed from the dynamic‐wave solution. It was found that the mean errors for both the approximations were independent of the type of rainfall excess distribution for KF₀² > 5, where K is the kinematic‐wave number and F₀ is the Froude number. Therefore, the regions (KF₀², F₀) where the kinematic‐wave and diffusive‐wave solutions would be fairly accurate and for any distribution of spatially varying rainfall, were characterized. The kinematic‐wave approximation was reasonably accurate, with a mean error of less than 5% and for the critical depth at the downstream end, for KF₀² ≥ 20 with F₀ ≤ 1; if the rainfall excess was concentrated in a portion of the plane, the field where the kinematic‐wave solution was found accurate, it was more limited and characterized for KF₀² > 35 with F₀ ≤ 1. The diffusive‐wave solution was in good agreement with the dynamic‐wave solution with a mean error of less than 5%, in the flow depth, for KF₀² ≥ 15 with F₀ ≤ 1; for rainfall excess concentrated in a portion of the plane, the accuracy of the diffusion wave solution was in a region more restricted and defined for KF₀² ≥ 30 with F₀ ≤ 1. For zero‐depth gradient at the downstream end, the accuracy field of the kinematic‐wave was found to be greater and characterized for KF₀² > 10 with F₀ ≤ 1; for rainfall excess concentrated in a portion of the plane, the region was smaller and defined for KF₀² > 15 with F₀ ≤ 1. The diffusive‐wave solution was found accurate in the region defined for KF₀² > 7·5, whereas for rainfall excess concentrated in a portion of the plane, the field of accuracy was for KF₀² > 12·5 with F₀ ≤ 1. The lower limits of the regions, defined on KF₀², can be considered generally valid for both approximations, but for F₀ < 1 smaller lower limits were also characterized. Finally, the accuracy of these approximations was influenced significantly by the downstream boundary condition. Copyright © 2002 John Wiley & Sons, Ltd.