Simulating the effects of geologic heterogeneity and transient boundary conditions on streambed temperatures — Implications for temperature-based water flux calculations

Analytical solutions to the one-dimensional heat transport equation for steady-state conditions can provide simple means to quantify groundwater surface water exchange. The errors in exchange flux calculations that are introduced when the underlying assumptions of homogeneous sediments and constant temperature boundary conditions are violated were systematically evaluated in a simulation study. Temperatures in heterogeneous sediments were simulated using a numerical model. Heterogeneity in the sediments was represented by discrete, binary geologic units. High contrasts between the hydraulic conductivities (K) of the geologic units were found to lead to large errors, while the influence of the structural arrangement of the units was smaller. The effects of transient temperature boundary conditions were investigated using an analytical equation. Errors introduced by transient boundary conditions were small for Darcy-velocities > 0.1 m d^− 1 in the period near maximum and minimum annual surface water temperatures. For smaller fluxes, however, errors can be large. Assuming steady-state conditions and vertical flow in homogeneous sediments is acceptable at certain times of the year and for medium to high exchange fluxes, but pronounced geologic heterogeneity can lead to large errors.     

Rangeland hydrology and erosion model (RHEM) enhancements for applications on disturbed rangelands

The rangeland hydrology and erosion model (RHEM) is a new process‐based model developed by the USDA Agricultural Research Service. RHEM was initially developed for functionally intact rangelands where concentrated flow erosion is minimal and most soil loss occurs by rain splash and sheet flow erosion processes. Disturbance such as fire or woody plant encroachment can amplify overland flow erosion by increasing the likelihood of concentrated flow formation. In this study, we enhanced RHEM applications on disturbed rangelands by using a new approach for the prediction and parameterization of concentrated flow erosion. The new approach was conceptualized based on observations and results of experimental studies on rangelands disturbed by fire and/or by tree encroachment. The sediment detachment rate for concentrated flow was calculated using soil erodibility and hydraulic (flow width and stream power) parameters. Concentrated flow width was calculated based on flow discharge and slope using an equation developed specifically for disturbed rangelands. Soil detachment was assumed to begin with concentrated flow initiation. A dynamic erodibility concept was applied where concentrated flow erodibility was set to decrease exponentially during a run‐off event because of declining sediment availability. Erodibility was estimated using an empirical parameterization equation as a function of vegetation cover and surface soil texture. A dynamic partial differential sediment continuity equation was used to model the total detachment rate of concentrated flow and rain splash and sheet flow. The enhanced version of the model was evaluated against rainfall simulation data for three different sites that exhibit some degree of disturbance by fire and/or by tree encroachment. The coefficient of determination (R²) and Nash–Sutcliffe efficiency were 0.78 and 0.71, respectively, which indicates the capability of the model using the new approach for predicting soil loss on disturbed rangeland. By using the new concentrated flow modelling approach, the model was enhanced to be a practical tool that utilizes readily available vegetation and soil data for quantifying erosion and assessing erosion risk following rangeland disturbance. Copyright © 2014 John Wiley & Sons, Ltd.     

Gully cut‐and‐fill cycles as related to agro‐management: a historical curve number simulation in the Tigray Highlands

Gully cut‐and‐fill dynamics are often thought to be driven by climate and/or deforestation related to population pressure. However, in this case‐study of nine representative catchments in the Northern Ethiopian Highlands, we find that neither climate changes nor deforestation can explain gully morphology changes over the twentieth century. Firstly, by using a Monte Carlo simulation to estimate historical catchment‐wide curve numbers, we show that the landscape was already heavily degraded in the nineteenth and early twentieth century – a period with low population density. The mean catchment‐wide curve number (> 80) one century ago was, under the regional climatic conditions, already resulting in considerable simulated historical runoff responses. Secondly, twentieth century land‐cover and runoff coefficient changes were confronted with twentieth century changing gully morphologies. As the results show, large‐scale land‐cover changes and deforestation cannot explain the observed processes. The study therefore invokes interactions between authigenic factors, small‐scale plot boundary changes, cropland management and sociopolitical forces to explain the gully cut processes. Finally, semi‐structured interviews and sedistratigraphic analysis of three filled gullies confirm the dominant impact of (crop)land management (tillage, check dams in gullies and channel diversions) on gully cut‐and‐fill processes. Since agricultural land management – including land tenure and land distribution – has been commonly neglected in earlier related research, we argue therefore that it can be a very strong driver of twentieth century gully morphodynamics. Copyright © 2014 John Wiley & Sons, Ltd.     

Application of Multiple Point Geostatistics to Non-stationary Images

Simulation of flow and solute transport through aquifers or oil reservoirs requires a precise representation of subsurface heterogeneity that can be achieved by stochastic simulation approaches. Traditional geostatistical methods based on variograms, such as truncated Gaussian simulation or sequential indicator simulation, may fail to generate the complex, curvilinear, continuous and interconnected facies distributions that are often encountered in real geological media, due to their reliance on two-point statistics. Multiple Point Geostatistics (MPG) overcomes this constraint by using more complex point configurations whose statistics are retrieved from training images. Obtaining representative statistics requires stationary training images, but geological understanding often suggests a priori facies variability patterns. This research aims at extending MPG to non-stationary facies distributions. The proposed method subdivides the training images into different areas. The statistics for each area are stored in separate frequency search trees. Several training images are used to ensure that the obtained statistics are representative. The facies probability distribution for each cell during simulation is calculated by weighting the probabilities from the frequency trees. The method is tested on two different object-based training image sets. Results show that non-stationary training images can be used to generate suitable non-stationary facies distributions.     

Detrital Cr-spinel in the Šambron–Kamenica Zone (Slovakia): evidence for an ocean-spreading zone in the Northern Vardar suture?

The Šambron–Kamenica Zone is situated on the northern margin of the Levočské vrchy mountains and Šarišskà vrchovina Highland, where the Central Carpathian Paleogene joins the Pieniny Klippen Belt. Sandstone outcrops in this area. From Cretaceous to Late Oligocene in age, these sediments suggest transport directions from S and SE. The heavy mineral assemblages of this sandstone include Cr-spinel grains, mainly displaying types II and III alpine-peridotite affinities, and are representative of Ocean Island Basalt volcanism. A sample from Upper Eocene sediments at Vit’az shows a clear change in Cr-spinel composition, which turns out to have types I and II peridotite affinities, and to derive from arc and Middle Ocean Ridge Basalt volcanism, with sediment transport directions from SW and WSW. These data indicate major variations in the Upper Eocene tectonic setting, giving constraints to paleogeographic reconstruction of the Slovak Central Carpathians.     

Accuracy of Laser Ablation U-Pb Zircon Dating: Results from a Test Using Five Different Reference Zircons

Using a state‐of‐the‐art 193 nm‐LA‐MC‐ICP‐MS system and with careful control of analytical procedures, the long term external reproducibility and accuracy of the ages Phanerozoic zircons measured over a period of months using calibrator bracketing for the ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²⁰⁶Pb ages were ca. 1% (2 RSD) if a single reference zircon was used for the matrix‐matched calibration. When different reference zircons were used for the calibration, suspicious systematic shifts in the obtained ages were observed and thus a reduction in the overall accuracy of the dating method became obvious. Such shifts were within a few percent range of the U‐Pb and Pb/Pb ages and seemed to vary independently of zircon age and composition. A “test of accuracy” experiment was conducted reducing instrumental effects as far as possible by analysing five different reference zircons mounted on a single mount eight times during the same session. An identical protocol was used for all analyses, with unchanged instrument parameters and with ion beam intensities kept as identical as possible. For data reduction, every zircon served consecutively as the reference zircon for calibration, with the others in the batch treated as unknowns. The known reference age and the four calculated ages obtained using the four other RMs for calibration were then compared. Even using such a strict analytical protocol, shifts in ²⁰⁶Pb/²³⁸U, ²⁰⁷Pb/²³⁵U and ²⁰⁷Pb/²⁰⁶Pb ratios were still present. They varied non‐systematically and ranged from ?4.35% to 3.08% for the investigated age range (1065 Ma to 226 Ma). Assuming the absence of instrumental effects (i.e., memory, dead‐time correction, non‐linearity of ion counters and interdetector calibration, crystallographic orientation, ablation cell geometry and setup, gas flows), the observed shifts were attributed to matrix and/or ablation related effects. It is proposed that non‐spectral matrix effects in the Ar plasma torch resulted in non‐uniform signal enhancement (or depression?) leading to shifts both in elemental and Pb isotopic ratios. Additionally, the ablated particle size distribution could be an important factor controlling plasma conditions and thus mass bias and fractionation. Until such effects are well understood and controlled, it would seem that any LA‐ICP‐MS zircon U‐Pb and ²⁰⁷Pb/²⁰⁶Pb age determination cannot be meaningfully interpreted at below a ca. 3% to 4% (2 RSD) confidence level.     

Biomarker lipids of the freshwater fern Azolla and its fossil counterpart from the Eocene Arctic Ocean

Eocene sediments recovered from the Lomonosov Ridge in the central Arctic Ocean during Integrated Ocean Drilling Program Expedition 302 contain high amounts of fossil remains of the free floating freshwater fern Azolla. Both extant Azolla and the sediments from the Arctic Azolla interval were found to contain relatively high quantities of compounds identified as 1,ω20 C₃₀–C₃₆ diols. Furthermore, structurally related mid-chain hydroxy fatty acids, long chain n-alkanols and keto-ols were discovered. The different series have a common feature: there is always a hydroxy group at the ω20 position. In addition, structurally related C₂₉ ω20,ω21 diols, C₂₉ 1,20,21 triols, C₂₉ dihydroxy fatty acids as well as a series of wax esters containing these mono- and dihydroxy lipids are reported. Selective extraction of Azolla surface lipids revealed that these compounds are most likely present in the leaf waxes of the fern. The suite of long-chain, mid-chain ω20 hydroxy wax constituents is described for the first time. As they are well preserved in the Eocene sediments, these lipids may serve as palaeo-environmental indicators of the Arctic Azolla interval and as markers for the past occurrence of Azolla in general.