Land-use affects the radiocarbon age, storage and depth distribution of soil organic carbon in Eastern Australia

Land-use has been shown to affect soil organic carbon (SOC) storage, with natural systems generally storing larger quantities of SOC than anthropogenically managed systems in surface soils. However, these effects are often difficult to detect deeper in the soil profile. Little is known regarding the effects of land-use on the radiocarbon age of SOC, both at the surface and deeper in the soil profile. We investigated the storage, radiocarbon content and depth distribution of soil organic carbon from across the state of NSW, Australia. A total of 100 profiles were analysed for total SOC concentration at numerous depths (up to 1 m) and a machine learning approach implementing tree ensemble methods was used to identify the key drivers of SOC depth distribution. Surface SOC storage was strongly associated with climate (predominately precipitation, to a lesser degree relative humidity and temperature), whereas SOC depth distribution was predominately influenced by land-use, soil type and to a lesser extent temperature. A subset of 12 soil profiles from a range of climate zones were analysed for radiocarbon content with a view to contrasting three land-use systems: natural, cleared/grazed and cropped. Radiocarbon content was affected strongly by land-use, with effects most pronounced at depth. Native systems appeared to have the youngest carbon throughout the profile, with cropped and grazed systems having older SOC. Radiocarbon content was also strongly associated with SOC content. Our results indicate that natural systems act as a carbon pump into the soil, injecting young, fresh organic carbon which is vertically distributed throughout the profile. In contrast, managed systems are deprived of this input and are depleted in SOC at all depths, leading to higher radiocarbon ages throughout the profile

A lattice grain model of hillslope evolution

Thispaperdescribesandexploresanewcontinuous-timestochasticcellularautomatonmodelofhill- slope evolution. The Grain Hill model provides a computational framework with which to study slope forms that arise from stochastic disturbance and rock weathering events. The model operates on a hexagonal lattice, with cell states representing fluid, rock, and grain aggregates that are either stationary or in a state of motion in one of the six cardinal lattice directions. Cells representing near-surface soil material undergo stochastic disturbance events, in which initially stationary material is put into motion. Net downslope transport emerges from the greater likelihood for disturbed material to move downhill than to move uphill. Cells representing rock undergo stochas- tic weathering events in which the rock is converted into regolith. The model can reproduce a range of common slope forms, from fully soil mantled to rocky or partially mantled, and from convex-upward to planar shapes. An optional additional state represents large blocks that cannot be displaced upward by disturbance events. With the addition of this state, the model captures the morphology of hogbacks, scarps, and similar features. In its simplest form, the model has only three process parameters, which represent disturbance frequency, characteris- tic disturbance depth, and base-level lowering rate, respectively. Incorporating physical weathering of rock adds one additional parameter, representing the characteristic rock weathering rate. These parameters are not arbitrary but rather have a direct link with corresponding parameters in continuum theory. Comparison between observed and modeled slope forms demonstrates that the model can reproduce both the shape and scale of real hillslope profiles. Model experiments highlight the importance of regolith cover fraction in governing both the downslope mass transport rate and the rate of physical weathering. Equilibrium rocky hillslope profiles are possible even when the rate of base-level lowering exceeds the nominal bare-rock weathering rate, because increases in both slope gradient and roughness can allow for rock weathering rates that are greater than the flat-surface maximum. Examples of transient relaxation of steep, rocky slopes predict the formation of a regolith-mantled pediment that migrates headward through time while maintaining a sharp slope break

Fresh shallow valleys in the Martian midlatitudes as features formed by meltwater flow beneath ice

Significant numbers of valleys have been identified in the Martian midlatitudes (30–60°N/S), spatially associated with extant or recent ice accumulations. Many of these valleys date to the Amazonian, but their formation during these cold, dry epochs is problematic. In this study, we look in detail at the form, distribution, and quantitative geomorphology of two suites of these valleys and their associated landforms in order to better constrain the processes of their formation. Since the valleys themselves are so young and thus well preserved, uniquely, we can constrain valley widths and courses and link these to the topography from the Mars Orbiter Laser Altimeter and High-Resolution Stereo Camera data. We show that the valleys are both qualitatively and quantitatively very similar, despite their being >5000 km apart in different hemispheres and around 7 km apart in elevation. Buffered crater counting indicates that the ages of these networks are statistically identical, probably forming during the Late Amazonian, ~100 Ma. In both localities, at least tens of valleys cross local drainage divides, apparently flowing uphill. We interpret these uphill reaches to be characteristic of flow occurring beneath a now absent, relatively thin (order 101–102 m), regionally extensive ice cover. Ridges and mounds occasionally found at the foot of these valley systems are analogous to eskers and aufeis-like refreezing features. On the basis of their interaction with these aufeis-like mounds, we suggest that this suite of landforms may have formed in a single, short episode (perhaps order of days), probably forced by global climate change

Dynamics of long term fluvial response in postglacial catchments of the Ladakh Batholith, Northwest Indian Himalaya

Upland rivers control the large-scale topographic form of mountain belts, allow coupling of climate and tectonics at the earth’s surface and are responsible for large scale redistribution of sediment from source areas to sinks. However, the details of how these rivers behave when perturbed by changes to their boundary conditions are not well understood. I have used a combination of fieldwork, remotely sensed data, mathematical analysis and computer modelling to investigate the response of channels to well constrained changes in the forcings upon them, focussing in particular on the effects of glacial remoulding of the catchments draining the south flank of the Ladakh batholith, northwest Indian Himalaya. The last glacial maximum for these catchments is atypically old (~100 ka), and this allows investigation of the response to glaciation on a timescale not usually available. The geomorphology of the catchments is divided into three distinct domains on the basis of the behaviour of the trunk stream – an upper domain where the channel neither aggrades above or incises into the valley form previously carved by glacial abrasion, a middle domain where the channel incises a gorge down into glacial sediments which mantle the valley floor, and a lower domain where the channel aggrades above this postglacial sediment surface. This landscape provides a framework in which to analyze the processes and timescales of fluvial response to glacial modification. The dimensions of the gorge and the known dates of glacial retreat record a time averaged peak river incision rate of approximately 0.5 mm/y; the timescale for the river long profile to recover to a smooth, concave up form must exceed 1 Ma. These values are comparable with those from similarly sized catchments that have been transiently perturbed by changing tectonics, but have never been quoted for a glacially forced basin-scale response. I have also demonstrated that lowering of the upper reaches of the Ladakh channel long profiles by glacial processes can systematically and nonlinearly perturb the slope-area (concavity) scaling of the channel downstream of the resulting profile convexities, or knickzones. The concavity values are elevated significantly above the expected equilibrium values of 0.3-0.6, with the magnitude controlled by the relative position of the knickzone within the catchment, and thus also by the degree of glacial modification of the fluvial system. This work also documents the existence of very similar trends in measured concavities downstream of long profile convexities in other transiently responding river systems in different tectonoclimatic settings, including those responding to changes in relative channel uplift. This previously unrecognised unity of response across a wide variety of different environments argues that such a trend is an intrinsic property of river response to perturbation. Importantly, it is consistent with the scaling expected from variation in incision efficiency driven by evolving sediment flux downstream of knickzones. The pervasive nature of this altered scaling, and its implications for fluvial erosion laws in perturbed settings, have significant consequences for efforts to interpret past changes in forcings acting on river systems from modern topography. I follow this by examining in detail the channel hydraulics of the Ladakh streams as they incise in response to the glacial perturbation. I present a new framework under which the style of erosion of a natural channel can be characterized as either detachment- or transport-limited based upon comparison of the downstream distribution of shear stress with the resulting magnitude of incision. This framework also allows assessment of the importance of sediment flux driven effects in studied channels. This approach is then used to demonstrate that fluvial erosion and deposition in the Ladakh catchments is best modelled as a sediment flux dependent, thresholded, detachment-limited system. The exceptional quality of the incision record in this landscape enables an unprecedented calibration of the sediment flux function within this incision law for three different trunk streams. The resulting curves are not compatible with the theoretically-derived parabolic form of this relation, instead showing nonzero erosion rates at zero sediment flux, a rapid rise and peak at relative sediment fluxes of less than 0.5 and a quasi exponential decrease in erosional efficiency beyond this. The position of the erosional efficiency peak in relative sediment flux space and the magnitude of the curve are shown to be both variable between the catchments explored and also correlated with absolute sediment flux in the streams

Modeling the shape and evolution of normal-fault facets

Facets formed along the footwalls of active normal-fault blocks display a variety of longitudinal profile forms, with variations in gradient, shape, degree of soil cover, and presence or absence of a slope break at the fault trace. We show that a two-dimensional, process-oriented cellular automaton model of facet profile evolution can account for the observed morphologic diversity. The model uses two dimensionless parameters to represent fault slip, progressive rock weathering, and downslope colluvial-soil transport driven by gravity and stochastic disturbance events. The parameters represent rock weathering and soil disturbance rates, respectively, scaled by fault slip rate; both can be derived from field-estimated rate coefficients. In the model's transport-limited regime, slope gradient depends on the ratio of disturbance to slip rate, with a maximum that represents the angle of repose for colluvium. In this regime, facet evolution is consistent with nonlinear diffusion models of soil-mantled hillslope evolution. Under the weathering-limited regime, bedrock becomes partly exposed but microtopography helps trap some colluvium even when facet gradient exceeds the threshold angle. Whereas the model predicts a continuous gradient from footwall to colluvial wedge under transport-limited behavior, fully weathering-limited facets tend to develop a slope break between footwall and basal colluvium as a result of reduced transport efficiency on the rocky footwall slope. To the extent that the model provides a reasonable analogy for natural facets, its behavior suggests that facet profile morphology can provide useful constraints on relative potential rates of rock weathering, soil disturbance, and fault slip

The role of superficial geology in controlling groundwater recharge in the weathered crystalline basement of semi-arid Tanzania

Study region: Little Kinyasungwe River Catchment, central semi-arid Tanzania. Study focus: The structure and hydraulic properties of superficial geology can play a crucial role in controlling groundwater recharge in drylands. However, the pathways by which groundwater recharge occurs and their sensitivity to environmental change remain poorly resolved. Geophysical surveys using Electrical Resistivity Tomography (ERT) were conducted in the study region and used to delineate shallow subsurface stratigraphy in conjunction with borehole logs. Based on these results, a series of local-scale conceptual hydrogeological models was produced and collated to generate a 3D conceptual model of groundwater recharge to the wellfield. New hydrological insights for the region: We propose that configurations of superficial geology control groundwater recharge in dryland settings as follows: 1) superficial sand deposits act as collectors and stores that slowly feed recharge into zones of active faulting; 2) these fault zones provide permeable pathways enabling greater recharge to occur; 3) ‘windows’ within layers of smectitic clay that underlie ephemeral streams may provide pathways for focused recharge via transmission losses; and 4) overbank flooding during high intensity precipitation events increases the probability of activating such permeable pathways. These conceptual models provide a physical basis to improve numerical models of groundwater recharge in drylands, and a conceptual framework to evaluate strategies (e.g., Managed Aquifer Recharge) to artificially enhance the availability of groundwater resources in these regions

The Landlab v1.0 OverlandFlow component: a Python tool for computing shallow-water flow across watersheds

Representation of flowing water in landscape evolution models (LEMs) is often simplified compared to hydrodynamic models, as LEMs make assumptions reducing physical complexity in favor of computational efficiency. The Landlab modeling framework can be used to bridge the divide between complex runoff models and more traditional LEMs, creating a new type of framework not commonly used in the geomorphology or hydrology communities. Landlab is a Python-language library that includes tools and process components that can be used to create models of Earth-surface dynamics over a range of temporal and spatial scales. The Landlab OverlandFlow component is based on a simplified inertial approximation of the shallow water equations, following the solution of de Almeida et al. (2012). This explicit two-dimensional hydrodynamic algorithm simulates a flood wave across a model domain, where water discharge and flow depth are calculated at all locations within a structured (raster) grid. Here we illustrate how the OverlandFlow component contained within Landlab can be applied as a simplified event-based runoff model and how to couple the runoff model with an incision model operating on decadal timescales. Examples of flow routing on both real and synthetic landscapes are shown. Hydrographs from a single storm at multiple locations in the Spring Creek watershed, Colorado, USA, are illustrated, along with a map of shear stress applied on the land surface by flowing water. The OverlandFlow component can also be coupled with the Landlab DetachmentLtdErosion component to illustrate how the nonsteady flow routing regime impacts incision across a watershed. The hydrograph and incision results are compared to simulations driven by steady-state runoff. Results from the coupled runoff and incision model indicate that runoff dynamics can impact landscape relief and channel concavity, suggesting that on landscape evolution timescales, the OverlandFlow model may lead to differences in simulated topography in comparison with traditional methods. The exploratory test cases described within demonstrate how the OverlandFlow component can be used in both hydrologic and geomorphic applications

CellLab-CTS 2015: continuous-time stochastic cellular automaton modeling using Landlab

CellLab-CTS 2015 is a Python-language software library for creating two-dimensional, continuous-time stochastic (CTS) cellular automaton models. The model domain consists of a set of grid nodes, with each node assigned an integer state code that represents its condition or composition. Adjacent pairs of nodes may undergo transitions to different states, according to a user-defined average transition rate. A model is created by writing a Python code that defines the possible states, the transitions, and the rates of those transitions. The code instantiates, initializes, and runs one of four object classes that represent different types of CTS models. CellLab-CTS provides the option of using either square or hexagonal grid cells. The software provides the ability to treat particular grid-node states as moving particles, and to track their position over time. Grid nodes may also be assigned user-defined properties, which the user can update after each transition through the use of a callback function. As a component of the Landlab modeling framework, CellLab-CTS models take advantage of a suite of Landlab's tools and capabilities, such as support for standardized input and output

Hesperian-Amazonian Transition Mid-Latitude Valleys: Markers of a Late Martian Climate Optima?

Recently the inventory of fluvial features that have been dated to the late Hesperian to early Amazonian epoch has increased dramatically, including a reassessment of the ages of the large alluvial fans and deltas (e.g., Eberswalde) to this time period. Mid-latitude Valleys (MLVs) are distinct from the older, more integrated Noachian-Hesperian Valley Networks which are deeply dissected, are generally of much larger spatial extent, and are more degraded. Although some MLVs involve rejuvenation of older Valley Networks, many MLVs are carved into smooth or rolling slopes and intercrater terrain. The MLVs range from a few meters to < 300 m in width, with nearly parallel valley walls and planforms that are locally sinuous. Although the MLVs in Newton and Gorgonum basins extend from the basin rims up to 75 km into the basin interior, most MLVs are shorter and often discontinuous. The occurrence of widespread MLVs suggest the possibility of their formation during one or perhaps more regional to global climatic episodes, possibly due to melting of seasonal to long-term accumulations of snow and ice. Temperatures warm enough to cause extensive melting may have occurred during optimal orbital and obliquity configurations, perhaps in conjunction with intensive volcanism releasing moisture and greenhouse gasses, or as a result of a brief episode of warming from a large impact. The concentration of MLVs to the northern and western basin slopes of Newton and Gorgonum basins suggests a possible aspect control to ice accumulation or melting. MLV activity occurred about at the same time as formation of the major outflow channels. A possible scenario is that delivery of water to the northern lowlands provided, through evaporation and sublimation, water that temporarily accumulated in the mid-southern latitudes as widespread ice deposits whose partial melting formed the MLVs and small, dominantly ice-covered lakes

Sedimentology and climatic environment of alluvial fans in the Martian Saheki Crater and a comparison with terrestrial fans in the Atacama Desert

The deflated surfaces of the alluvial fans in Saheki crater reveal the most detailed record of fan stratigraphy and evolution found, to date, on Mars. During deposition of at least the uppermost 100 m of fan deposits, discharges from the source basin consisted of channelized flows transporting sediment (which we infer to be primarily sand- and gravel-sized) as bedload coupled with extensive overbank mud-rich flows depositing planar beds of sand-sized or finer sediment. Flow events are inferred to have been of modest magnitude (probably less than ~60 m3/s), of short duration, and probably occupied only a few distributaries during any individual flow event. Occasional channel avulsions resulted in the distribution of sediment across the entire fan. A comparison with fine-grained alluvial fans in Chile’s Atacama Desert provides insights into the processes responsible for constructing the Saheki crater fans: sediment is deposited by channelized flows (transporting sand through boulder-sized material) and overbank mudflows (sand size and finer) and wind erosion leaves channels expressed in inverted topographic relief. The most likely source of water was snowmelt released after annual or epochal accumulation of snow in the headwater source basin on the interior crater rim during the Hesperian to Amazonian periods. We infer the Saheki fans to have been constructed by many hundreds of separate flow events, and accumulation of the necessary snow and release of meltwater may have required favorable orbital configurations or transient global warming