Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.     

Watershed structural influences on the distributions of stream network water and solute travel times under baseflow conditions

Watershed structure influences the timing, magnitude, and spatial location of water and solute entry to stream networks. In turn, stream reach transport velocities and stream network geometry (travel distances) further influence the timing of export from watersheds. Here, we examine how watershed and stream network organization can affect travel times of water from delivery to the stream network to arrival at the watershed outlet. We analysed watershed structure and network geometry and quantified the relationship between stream discharge and solute velocity across six study watersheds (11.4 to 62.8 km²) located in the Sawtooth Mountains of central Idaho, USA. Based on these analyses, we developed stream network travel time functions for each watershed. We found that watershed structure, stream network geometry, and the variable magnitude of inputs across the network can have a pronounced affect on water travel distances and velocities within a stream network. Accordingly, a sample taken at the watershed outlet is composed of water and solutes sourced from across the watershed that experienced a range of travel times in the stream network. We suggest that understanding and quantifying stream network travel time distributions are valuable for deconvolving signals observed at watershed outlets into their spatial and temporal sources, and separating terrestrial and in‐channel hydrological, biogeochemical, and ecological influences on in‐stream observations. Copyright © 2016 John Wiley & Sons, Ltd.     

Interpreting zirconium‐in‐rutile thermometric results

At first sight, experimental results and observations on rocks suggest that the Zr content in rutile, where equilibrated with quartz and zircon, should be a useful thermometer for metamorphic rocks. However, diffusion data for Zr in rutile imply that thermometry should not, for plausible rates of cooling, give the high temperatures commonly observed in high‐grade metamorphic rocks. It is suggested here that such observations can be accounted for by high‐T diffusive closure of Si in rutile, causing the interior of rutile grains to become insensitive to the thermometer equilibrium well above the temperature of Zr diffusive closure. Paired with comparatively slow grain boundary diffusion and problematic zircon nucleation, this allows for cases of Zr retention in rutile through temperatures where Zr is still diffusively mobile within rutile grains. Other observations that may be accounted for in this context are large inter‐grain ranges of rutile Zr contents uncorrelated with rutile grain size, and flat Zr profiles across individual rutile grains, counter to what would be expected from diffusive closure. A consequence is that it is unlikely that Zr‐in‐rutile thermometry will be useful for estimating rock cooling rates.     

Diatom-inferred depth models in 8 Canadian boreal lakes: inferred changes in the benthic:planktonic depth boundary and implications for assessment of past droughts

Assessment of past drought in boreal regions, a region predicted to be strongly affected by climate warming, can provide insights into future availability of water. However, limited instrumental data and paleoclimatic data are available for this assessment. To address this lack of data in the boreal region of northwest Ontario, a regional study of lakes in the Winnipeg River Drainage Basin was initiated. Diatom-inferred (D-I) depth models were developed based on surface samples collected along a depth gradient within 8 small boreal lakes. Weighted-averaging and modern analog approaches provided robust within-lake depth models for each of the study lakes, with bootstrapped r² values ranging from 0.90 to 0.98, and root-mean-squared-errors of prediction (RMSEP) between 1.1 and 2.5 m. Large differences in the estimated depth optima for three representative, but common diatom species across our 8 study lakes suggested that within-lake calibration datasets are more appropriate for inferring past drought based in depth models than a regional multi-lake calibration dataset, and that light and related variables are controlling factors governing the maximum depth of benthic taxa. A down-core application of the D-I depth models on a near-shore core from Meekin Lake, retrieved near the present-day ecotone between the benthic and planktonic diatom assemblages indicated highly similar trends in inferred depth (r = 0.96). The models have significant correlations with other metrics of changes in depth including diatom species richness (r = 0.74–0.78) and evenness (r = 0.76–0.8), thereby allowing a check on the strength and direction of the depth inferences down-core. Near-shore cores located near the benthic:planktonic transition is a sensitive region that can provide estimates of past droughts in lakes where such inferences have been difficult to estimate.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.     

Stromatoporoid growth forms and Devonian reef fabrics in the Upper Devonian Alexandra Reef System,Canada – Insight on the challenges of applying Devonian reef facies models

Existing facies models for Devonian reef systems can be divided into high‐energy and low‐energy types. A number of assumptions have been made in the development of these models and, in some cases, criteria that distinguish important aspects of the models are poorly defined. The Upper Devonian Alexandra Reef System contains a variety of reef fabrics from different depositional environments and is ideal for studying the range of environments in which stromatoporoids thrived and the facies from these different environments. A wide variety of stromatoporoid growth forms including laminar, tabular, anastamosing laminar and tabular, domal, bulbous, dendroid, expanding conical, concave‐up whorled‐laminar, concave‐up massive tabular and platy‐multicolumnar are present in the Alexandra Reef System. The whorled‐laminar and massive tabular concave‐up growth forms are virtually undocumented from other Devonian reefs but were common in the reef front of the Alexandra, where they thrived in a low‐energy environment around and below fair‐weather wave base. In contrast, high‐energy parts of the reef margin were dominated by bioclastic rubble deposits with narrow ribbon‐like discontinuous bodies of laminar stromatoporoid framestone. In the lagoon, laminar stromatoporoids formed steep‐sided sediment‐dominated bioherms in response to sea‐level rise and flooding. Relying mostly on the different reef facies in the Alexandra system, a new classification scheme for Devonian reef fabrics has been developed. Devonian reef fabrics can be classified as being: (i) sediment‐laden metazoan dominated; (ii) metazoan–microbial dominated (boundstone); (iii) metazoan dominated (framestone); or (iv) metazoan–marine cement dominated. Distinction of these fabrics carries important sedimentary and palaeoecological implications for reconstructing the depositional environment. With examples from the Alexandra Formation, it is demonstrated that reef facies accumulated in a range of depositional environments and that the simple observation of massive stromatoporoids with or without microbial deposits does not automatically imply a high‐energy reef margin, as otherwise portrayed in a number of the existing facies models for these systems.