Modelling increased soil cohesion due to roots with EUROSEM

As organic root exudates cause soil particles to adhere firmly to root surfaces, roots significantly increase soil strength and therefore also increase the resistance of the topsoil to erosion by concentrated flow. This paper aims at contributing to a better prediction of the root effects on soil erosion rates in the EUROSEM model, as the input values accounting for roots, presented in the user manual, do not account for differences in root density or root architecture. Recent research indicates that small changes in root density or differences in root architecture considerably influence soil erosion rates during concentrated flow. The approach for incorporating the root effects into this model is based on a comparison of measured soil detachment rates for bare and for root‐permeated topsoil samples with predicted erosion rates under the same flow conditions using the erosion equation of EUROSEM. Through backwards calculation, transport capacity efficiencies and corresponding soil cohesion values can be assessed for bare and root‐permeated topsoils respectively. The results are promising and present soil cohesion values that are in accordance with reported values in the literature for the same soil type (silt loam). The results show that grass roots provide a larger increase in soil cohesion as compared with tap‐rooted species and that the increase in soil cohesion is not significantly different under wet and dry soil conditions, either for fibrous root systems or for tap root systems. Power and exponential relationships are established between measured root density values and the corresponding calculated soil cohesion values, reflecting the effects of roots on the resistance of the topsoil to concentrated flow incision. These relationships enable one to incorporate the root effect into the soil erosion model EUROSEM, through adapting the soil cohesion input value. A scenario analysis shows that the contribution of roots to soil cohesion is very important for preventing soil loss and reducing runoff volume. The increase in soil shear strength due to the binding effect of roots on soil particles is two orders of magnitude lower as compared with soil reinforcement achieved when roots mobilize their tensile strength during soil shearing and root breakage. Copyright © 2008 John Wiley & Sons, Ltd.     

Large debris avalanches and associated eruptions in the Holocene eruptive history of Shiveluch Volcano, Kamchatka, Russia

 Shiveluch Volcano, located in the Central Kamchatka Depression, has experienced multiple flank failures during its lifetime, most recently in 1964. The overlapping deposits of at least 13 large Holocene debris avalanches cover an area of approximately 200 km² of the southern sector of the volcano. Deposits of two debris avalanches associated with flank extrusive domes are, in addition, located on its western slope. The maximum travel distance of individual Holocene avalanches exceeds 20 km, and their volumes reach ∼3 km³. The deposits of most avalanches typically have a hummocky surface, are poorly sorted and graded, and contain angular heterogeneous rock fragments of various sizes surrounded by coarse to fine matrix. The deposits differ in color, indicating different sources on the edifice. Tephrochronological and radiocarbon dating of the avalanches shows that the first large Holocene avalanches were emplaced approximately 4530–4350 BC. From ∼2490 BC at least 13 avalanches occurred after intervals of 30–900 years. Six large avalanches were emplaced between 120 and 970 AD, with recurrence intervals of 30–340 years. All the debris avalanches were followed by eruptions that produced various types of pyroclastic deposits. Features of some surge deposits suggest that they might have originated as a result of directed blasts triggered by rockslides. Most avalanche deposits are composed of fresh andesitic rocks of extrusive domes, so the avalanches might have resulted from the high magma supply rate and the repetitive formation of the domes. No trace of the 1854 summit failure mentioned in historical records has been found beyond 8 km from the crater; perhaps witnesses exaggerated or misinterpreted the events. Received: 18 August 1997 / Accepted: 19 December 1997     

Biot dispersion for P‐ and S‐wave velocities in partially and fully saturated sandstones

Ultrasonic compressional‐ and shear‐wave velocities have been measured on 34 samples of sandstones from hydrocarbon reservoirs. The sandstones are all of low clay content, high porosity, and cover a wide range of permeabilities. They were measured dry and brine‐saturated under hydrostatic effective stresses of 10, 20 and 40 MPa. For eight of the sandstones, ultrasonic velocity measurements were made at different partial water saturations in the range from dry to fully saturated. The Gassmann–Biot theory is found to account for most of the changes in velocities at high effective stress levels when the dry sandstones are fully saturated with brine, provided the lower velocities resulting when the dry sandstone initially adsorbs small amounts of moisture are used to determine the elastic properties of the ‘dry’ sandstone. At lower effective stress levels, local flow phenomena due to the presence of open microcracks are assumed to be responsible for measured velocities higher than those predicted by the theory. The partial saturation results are modelled fairly closely by the Gassmann–Biot theory, assuming heterogeneous saturation for P‐waves.     

Process‐based modeling of kilometer‐scale alongshore sandbar variability

Subtidal nearshore sandbars may exhibit cyclic net offshore migration during their multi‐annual lifetime along many sandy coasts. Although this type of behavior can extend continuously for several kilometers, alongshore variations in cross‐shore bar position and bar amplitude are commonly observed. Alongshore variability is greatest when bars display km‐scale disruptions, indicative of a distinct alongshore phase shift in the bar cycle. An outer bar is then attached to an inner bar, forming a phenomenon known as a bar switch. Here, we investigate such large‐scale alongshore variability using a process‐based numerical profile model and observations at 24 transects along a 6 km section of the barred beach at Noordwijk, The Netherlands. When alongshore variability is limited, the model predicts that the bars migrate offshore at approximately the same rate (i.e. the bars remain in phase). Only under specific bar configurations with high wave‐energy levels is an increase in the alongshore variability predicted. This suggests that cross‐shore processes may trigger a switch in the case of specific antecedent morphological configurations combined with storm conditions. It is expected that three‐dimensional (3D) flow patterns augment the alongshore variability in such instances. In contrast to the observed bar behaviour, predicted bar morphologies on either side of a switch remain in different phases, even though the bars are occasionally located at a similar cross‐shore position. In short, the 1D model is not able to remove a bar switch. This data‐model mismatch suggests that 3D flow patterns are key to the dissipation of bar switches. Copyright © 2014 John Wiley & Sons, Ltd.     

Multi‐year water balance assessment of a newly constructed wetland,Fort McMurray,Alberta

Oil sands mining in Alberta transforms the boreal landscape of forests and wetlands into open pits, tailings ponds and overburden piles. Whereas reclamation efforts have primarily focused on upland forests, rebuilding wetland systems has recently become a motivation for research. Wetland creation and sustainability in this region is complicated by the sub‐humid climate and salinity of underlying mining material. In 2012, Syncrude Canada Ltd. completed the construction of the Sandhill Fen Watershed (SFW), a 52‐ha upland‐wetland system to evaluate wetland reclamation strategies on soft tailings. SFW includes an active pumping system, upland hummocks, a fen wetland and underdrains. To evaluate the influence of management practices on the hydrology of the system, this study reports the water balance from January 2013 to December 2014, the first 2 years after commissioning. A semi‐distributed approach was taken to examine the fluxes and stores of water in uplands and lowlands. Natural and artificial inputs and outputs were measured using a series of precipitation gauges and pumps, and evapotranspiration was quantified using three eddy covariance towers. A series of near surface wells recorded water table position. Both 2013 and 2014 were normal rainfall years, with 2013 having more and 2014 less snow than normal. In 2013, inflow/outflow from pumping was the predominant hydrological fluxes, resulting in considerable variability in water table position and storage changes throughout the summer. In 2014, the artificial addition of water was negligible, yet the water table remained near the surface in lowland locations, suggesting that wetland conditions could be maintained under current conditions. Evapotranspiration rates between uplands and lowlands were similar between years and sites, ranging from 2.2 ± 1.8 to 2.5 ± 1.2 mm/day and were largely controlled by climate. These rates were less than nearby older upland systems, suggesting that water balance partitioning will change as vegetation develops. Comparison between years and with natural systems provides insight on how management practices influence hydrologic dynamics and the overall water balance of the SFW. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of pressure on aqueous chemical equilibria at subzero temperatures with applications to Europa

Pressure plays a critical role in controlling aqueous geochemical processes in deep oceans and deep ice. The putative ocean of Europa could have pressures of 1200 bars or higher on the seafloor, a pressure not dissimilar to the deepest ocean basin on Earth (the Mariana Trench at 1100 bars of pressure). At such high pressures, chemical thermodynamic relations need to explicitly consider pressure. A number of papers have addressed the role of pressure on equilibrium constants, activity coefficients, and the activity of water. None of these models deal, however, with processes at subzero temperatures, which may be important in cold environments on Earth and other planetary bodies. The objectives of this work were to (1) incorporate a pressure dependence into an existing geochemical model parameterized for subzero temperatures (FREZCHEM), (2) validate the model, and (3) simulate pressure-dependent processes on Europa. As part of objective 1, we examined two models for quantifying the volumetric properties of liquid water at subzero temperatures: one model is based on the measured properties of supercooled water, and the other model is based on the properties of liquid water in equilibrium with ice.The relative effect of pressure on solution properties falls in the order: equilibrium constants(K) > activity coefficients (γ) > activity of water (a_w). The errors (%) in our model associated with these properties, however, fall in the order: γ > K > a_w. The transposition between K and γ is due to a more accurate model for estimating K than for estimating γ. Only activity coefficients are likely to be significantly in error. However, even in this case, the errors are likely to be only in the range of 2 to 5% up to 1000 bars of pressure. Evidence based on the pressure/temperature melting of ice and salt solution densities argue in favor of the equilibrium water model, which depends on extrapolations, for characterizing the properties of liquid water in electrolyte solutions at subzero temperatures, rather than the supercooled water model. Model-derived estimates of mixed salt solution densities and chemical equilibria as a function of pressure are in reasonably good agreement with experimental measurements.To demonstrate the usefulness of this low-temperature, high-pressure model, we examined two hypothetical cases for Europa. Case 1 dealt with the ice cover of Europa, where we asked the question: How far above the putative ocean in the ice layer could we expect to find thermodynamically stable brine pockets that could serve as habitats for life? For a hypothetical nonconvecting 20 km icy shell, this potential life zone only extends 2.8 km into the icy shell before the eutectic is reached. For the case of a nonconvecting icy shell, the cold surface of Europa precludes stable aqueous phases (habitats for life) anywhere near the surface. Case 2 compared chemical equilibria at 1 bar (based on previous work) with a more realistic 1460 bars of pressure at the base of a 100 km Europan ocean. A pressure of 1460 bars, compared to 1 bar, caused a 12 K decrease in the temperature at which ice first formed and a 11 K increase in the temperature at which MgSO₄·12H₂O first formed. Remarkably, there was only a 1.2 K decrease in the eutectic temperatures between 1 and 1460 bars of pressure. Chemical systems and their response to pressure depend, ultimately, on the volumetric properties of individual constituents, which makes every system response highly individualistic.     

Depths to density discontinuities beneath the Adamawa Plateau region, Central Africa, from spectral analyses of new and existing gravity data

New gravity data from the Adamawa Uplift region of Cameroon have been integrated with existing gravity data from central and western Africa to examine variations in crustal structure throughout the region. The new data reveal steep northeast-trending gradients in the Bouguer gravity anomalies that coincide with the Sanaga Fault Zone and the Foumban Shear Zone, both part of the Central African Shear Zone lying between the Adamawa Plateau and the Congo Craton. Four major density discontinuities in the lithosphere have been determined within the lithosphere beneath the Adamawa Uplift in central Cameroon using spectral analysis of gravity data: (1) 7–13 km; (2) 19–25 km; (3) 30–37 km; and (4) 75–149 km. The deepest density discontinuities determined at 75–149 km depth range agree with the presence of an anomalous low velocity upper mantle structure at these depths deduced from earlier teleseismic delay time studies and gravity forward modelling. The 30–37 km depths agree with the Moho depth of 33 km obtained from a seismic refraction experiment in the region. The intermediate depth of 20 km obtained within region D may correspond to shallower Moho depth beneath parts of the Benue and Yola Rifts where seismic refraction data indicate a crustal thickness of 23 km. The 19–20 km depths and 8–12 km depths estimated in boxes encompassing the Adamawa Plateau and Cameroon Volcanic Line may may correspond to mid-crustal density contrasts associated with volcanic intrusions, as these depths are less than depths of 25 and 13 km, respectively, in the stable Congo Craton to the south.