The influence of tree root water uptake on the long term hydrology of a clay fill railway embankment

This paper uses a numerical model to investigate the influence of tree root water uptake and tree removal on pore water pressures and the vertical movement of a clay fill railway embankment. Simulated results of soil wetting and drying are compared with field measurements from an instrumented railway embankment before and after tree removal. A parametric study compares the influence of vegetation on the seasonal movement of the embankment slope. The simulations and field measurements show that while trees cause significant seasonal variations in pore water pressure and water content near the soil surface, they can maintain persistent soil suctions at depth within the tree rooting zone. Demonstration of this result using a numerical model requires a root water uptake function that separates spatially the processes of water infiltration, evaporation and transpiration. When all of the trees are removed, the persistent soil suctions established by the trees are lost as water infiltrates from the soil surface. Leaving the trees in place over the bottom third of the slope can maintain persistent suctions at the slope toe, while potentially also reducing seasonal ground movements at the crest that may adversely affect railway track geometry

The Kolkata Paise Restaurant Problem and Resource Utilization

We study the dynamics of the "Kolkata Paise Restaurant problem". The problem is the following: In each period, N agents have to choose between N restaurants. Agents have a common ranking of the restaurants. Restaurants can only serve one customer. When more than one customer arrives at the same restaurant, one customer is chosen at random and is served; the others do not get the service. We first introduce the one-shot versions of the Kolkata Paise Restaurant problem which we call one-shot KPR games. We then study the dynamics of the Kolkata Paise Restaurant problem (which is a repeated game version of any given one shot KPR game) for large N. For statistical analysis, we explore the long time steady state behavior. In many such models with myopic agents we get under-utilization of resources, that is, we get a lower aggregate payoff compared to the social optimum. We study a number of myopic strategies, focusing on the average occupation fraction of restaurants.Comment: revtex4, 8 pages, 3 figs, accepted in Physica

Accurate crop yield predictions from modelling tree-crop interactions in gliricidia-maize agroforestry

Agroforestry systems, containing mixtures of trees and crops, are often promoted because the net effect of interactions between woody and herbaceous components is thought to be positive if evaluated over the long term. From a modelling perspective, agroforestry has received much less attention than monocultures. However, for the potential of agroforestry to impact food security in Africa to be fully evaluated, models are required that accurately predict crop yields in the presence of trees. The positive effects of the fertiliser tree gliricidia (Gliricidia sepium) on maize (Zea mays) are well documented and use of this tree-crop combination to increase crop production is expanding in several African countries. Simulation of gliricidia-maize interactions can complement field trials by predicting crop response across a broader range of contexts than can be achieved by experimentation alone. We tested a model developed within the APSIM framework. APSIM models are widely used for one dimensional (1D), process-based simulation of crops such as maize and wheat in monoculture. The Next Generation version of APSIM was used here to test a 2D agroforestry model where maize growth and yield varied spatially in response to interactions with gliricidia. The simulations were done using data for gliricidia-maize interactions over two years (short-term) in Kenya and 11 years (long-term) in Malawi, with differing proportions of trees and crops and contrasting management. Predictions were compared with observations for maize grain yield, and soil water content. Simulations in Kenya were in agreement with observed yields reflecting lower observed maize germination in rows close to gliricidia. Soil water content was also adequately simulated, except for a tendency for slower simulated drying of the soil profile each season. Simulated maize yields in Malawi were also in agreement with observations. Trends in soil carbon over a decade were similar to those measured, but could not be statistically evaluated. These results show that the agroforestry model in APSIM Next Generation adequately represented tree-crop interactions in these two contrasting agro-ecological conditions and agroforestry practices. Further testing of the model is warranted to explore tree-crop interactions under a wider range of environmental conditions

Hillslope and stream connectivity: simulation of concentration-discharge patterns using the HYDRUS model

Abstract: Nutrient concentrations and loads in streamflow are sensitive to rapidly changing stream chemistry and discharge during storms. Mechanistic models that can simulate water and solute movement at hillslope scales could be useful for predicting concentration-discharge (C-Q) patterns and thereby improve our quantitative understanding of terrestrial-aquatic linkages for targeted catchment management. Our objective was to use the HYDRUS model to represent hydro-biogeochemical processes in soils that drive seepage of water and solutes from soil profiles into streams. Specifically we compared measurements in the literature with HYDRUS outputs using two methods for simulating runoff. This model predicts runoff (R) as rainfall that is instantaneously in excess of infiltration, but it is not designed to route runoff as overland flow. Post-HYDRUS addition of seepage to runoff was used to simulate the delivery of dissolved or particulate constituents to a stream (method A). Alternatively, we demonstrated how simulations using HYDRUS could include a hypothetical layer at the top of the soil profile with extremely high porosity and hydraulic conductivity that enabled overland flow and down-slope infiltration, but in this case only dissolved constituents could be considered (method B). These methods were evaluated by comparing the simulated temporal patterns of discharge and concentration with observed patterns. The catchments considered were in Slovenia (4210 ha) and in Australia (11.9 ha). Methods A and B were shown to adequately simulate some aspects of published discharge-concentration patterns, e.g. runoff dilution or concentration effects, but the temporal patterns of discharge for both methods did not precisely match those measured at small time-steps (e.g. 15 minutes). This limitation was due mainly to inadequate simulation of the down-slope movement of runoff and down-slope infiltration of a portion of this runoff. Method A was generally more useful than method B. Despite this limitation, both methods, if used carefully, should be adequate for many purposes, especially when simulating longer time-steps. Additional hypothetical simulations illustrated the significance of soil hydraulic conductivity, soil water content, and vertical gradients in solute concentrations in soil. Two temporal types of dischargeconcentration patterns were observed; short-term hysteresis caused by runoff during and shortly after a rainfall event, and longer-term trends associated with infiltration and seepage. Clockwise and anti-clockwise hysteresis was demonstrated to be potentially due to the temporal asynchrony of peak discharge and peak concentration in runoff. Simulations also demonstrated advantages over using the more common approach of a 2-or 3-component mixing model. Our results suggest that the HYDRUS model will be useful for the mechanistic simulation of within-soil processes that are needed to predict discharge-concentration patterns at hillslope scales

Distinct responses of neurons and astrocytes to TDP-43 proteinopathy in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal and incurable neurodegenerative disease caused by motor neuron loss, resulting in muscle wasting, paralysis and eventual death. A key pathological feature of ALS is cytoplasmically mislocalized and aggregated TDP-43 protein in >95% of cases, which is considered to have prion-like properties. Historical studies have predominantly focused on genetic forms of ALS, which represent ∼10% of cases, leaving the remaining 90% of sporadic ALS relatively understudied. Additionally, the role of astrocytes in ALS and their relationship with TDP-43 pathology is also not currently well understood. We have therefore used highly enriched human induced pluripotent stem cell (iPSC)-derived motor neurons and astrocytes to model early cell type-specific features of sporadic ALS. We first demonstrate seeded aggregation of TDP-43 by exposing human iPSC-derived motor neurons to serially passaged sporadic ALS post-mortem tissue (spALS) extracts. Next, we show that human iPSC-derived motor neurons are more vulnerable to TDP-43 aggregation and toxicity compared with their astrocyte counterparts. We demonstrate that these TDP-43 aggregates can more readily propagate from motor neurons into astrocytes in co-culture paradigms. We next found that astrocytes are neuroprotective to seeded aggregation within motor neurons by reducing (mislocalized) cytoplasmic TDP-43, TDP-43 aggregation and cell toxicity. Furthermore, we detected TDP-43 oligomers in these spALS spinal cord extracts, and as such demonstrated that highly purified recombinant TDP-43 oligomers can reproduce this observed cell-type specific toxicity, providing further support to a protein oligomer-mediated toxicity hypothesis in ALS. In summary, we have developed a human, clinically relevant, and cell-type specific modelling platform that recapitulates key aspects of sporadic ALS and uncovers both an initial neuroprotective role for astrocytes and the cell type-specific toxic effect of TDP-43 oligomers

Challenges in monitoring and managing engineered slopes in a changing climate

Geotechnical asset owners need to know which parts of their asset network are vulnerable to climate change induced failure in order to optimise future investment. Protecting these vulnerable slopes requires monitoring systems capable of identifying and alerting to asset operators changes in the internal conditions that precede failure. Current monitoring systems are heavily reliant on point sensors which can be difficult to interpret across slope scale. This paper presents challenges to producing such a system and research being carried out to address some of these using electrical resistance tomography (ERT). Experimental results show that whilst it is possible to measure soil water content indirectly via resistivity the relationship between resistivity and water content will change over time for a given slope. If geotechnical parameters such as pore water pressure are to be estimated using this method then ERT systems will require integrating with more conventional geotechnical instrumentation to ensure correct representative information is provided. The paper also presents examples of how such data can be processed and communicated to asset owners for the purposes of asset management

Galaxy Zoo: evidence for diverse star formation histories through the green valley

Does galaxy evolution proceed through the green valley via multiple pathways or as a single population? Motivated by recent results highlighting radically different evolutionary pathways between early- and late-type galaxies, we present results from a simple Bayesian approach to this problem wherein we model the star formation history (SFH) of a galaxy with two parameters, [t, τ] and compare the predicted and observed optical and near-ultraviolet colours. We use a novel method to investigate the morphological differences between the most probable SFHs for both disc-like and smooth-like populations of galaxies, by using a sample of 126 316 galaxies (0.01 < z < 0.25) with probabilistic estimates of morphology from Galaxy Zoo. We find a clear difference between the quenching time-scales preferred by smooth- and disc-like galaxies, with three possible routes through the green valley dominated by smooth- (rapid time-scales, attributed to major mergers), intermediate- (intermediate time-scales, attributed to minor mergers and galaxy interactions) and disc-like (slow time-scales, attributed to secular evolution) galaxies. We hypothesize that morphological changes occur in systems which have undergone quenching with an exponential time-scale τ < 1.5 Gyr, in order for the evolution of galaxies in the green valley to match the ratio of smooth to disc galaxies observed in the red sequence. These rapid time-scales are instrumental in the formation of the red sequence at earlier times; however, we find that galaxies currently passing through the green valley typically do so at intermediate time-scales.

In situ measurements of near-surface hydraulic conductivity in engineered clay slopes

In situ measurements of near-saturated hydraulic conductivity in fine grained soils have been made at six exemplar UK transport earthwork sites: three embankment and three cutting slopes. This paper reports 143 individual measurements and considers the factors that influence the spatial and temporal variability obtained. The test methods employed produce near-saturated conditions and flow under constant head. Full saturation is probably not achieved due to preferential and by-pass flow occurring in these desiccated soils. For an embankment, hydraulic conductivity was found to vary by five orders of magnitude in the slope near-surface (0 to 0.3 metres depth), decreasing by four orders of magnitude between 0.3 and 1.2 metres depth. This extremely high variability is in part due to seasonal temporal changes controlled by soil moisture content, which can account for up to 1.5 orders of magnitude of this variability. Measurements of hydraulic conductivity at a cutting also indicated a four orders of magnitude range of hydraulic conductivity for the near-surface, with strong depth dependency of a two orders of magnitude decrease from 0.2 to 0.6 metres depth. The main factor controlling the large range is found to be spatial variability in the soil macro structure generated by wetting/drying cycle driven desiccation and roots. The measurements of hydraulic conductivity reported in this paper were undertaken to inform and provide a benchmark for the hydraulic parameters used in numerical models of groundwater flow. This is an influential parameter in simulations incorporating the combined weather/vegetation/infiltration/soil interaction mechanisms that are required to assess the performance and deterioration of earthwork slopes in a changing climate