Climate driven coevolution of weathering profiles and hillslope topography generates dramatic differences in critical zone architecture

Considerable debate revolves around the relative importance of rock type, tectonics, and climate in creating the architecture of the critical zone. We demonstrate the importance of climate and in particular the rate of water recharge to the subsurface, using numerical models that incorporate hydrologic flowpaths, chemical weathering, and geomorphic rules for soil production and transport. We track alterations in both solid phase (plagioclase to clay) and water chemistry along hydrologic flowpaths that include lateral flow beneath the water table. To isolate the role of recharge, we simulate dry and wet cases and prescribe identical landscape evolution rules. The weathering patterns that develop differ dramatically beneath the resulting parabolic interfluves. In the dry case, incomplete weathering is shallow and surface parallel, whereas in the wet case, intense weathering occurs to depths approximating the base of the bounding channels, well below the water table. Exploration of intermediate cases reveals that the weathering state of the subsurface is strongly governed by the ratio of the rate of advance of the weathering front itself controlled by the water input rate, and the rate of erosion of the landscape. The system transitions between these end‐member behaviours rather abruptly at a weathering front speed ‐ erosion rate ratio of approximately 1. Although there are undoubtedly direct roles for tectonics and rock type in critical zone architecture, and yet more likely feedbacks between these and climate, we show here that differences in hillslope‐scale weathering patterns can be strongly controlled by climate.     

Ecosystem engineering by hummock-building earthworms in seasonal wetlands of eastern South Africa: Insights into the mechanics of biomorphodynamic feedbacks in wetland ecosystems

This paper resolves the origin of clay hummock micro-topography in seasonal wetlands of the Drakensberg Foothills, providing a review and appraisal of previously-suggested mechanisms of hummock formation in the context of new field and laboratory data. Field surveys revealed neo-formation of clay hummocks in a river channel that had been abandoned in c.1984. Fresh earthworm castings were located atop hummocks protruding from inundated abandoned channel margins. Earthworm castings, and sediment cores taken in hummocks and adjacent hollows, were analysed for soil-adsorbed carbon and nitrogen using an HCN analyser, and for ²¹⁰Pb activity using alpha-geochronology. ²¹⁰Pb activity profiles suggest relative enrichment of the isotope in hummocks, and relative depletion in adjacent hollows. Earthworm castings are characterised by very high ²¹⁰Pb activity, as well as high C and N contents. Hummocks have significantly higher C and N contents than adjacent hollows. Results suggest that it is the foraging activity of earthworms in litter-rich seasonal wetland hollows, and repeated excretion of castings atop adjacent hummocks, that is responsible for the elemental enrichment observed. The paper presents a conceptual model of hummock formation in wetlands through interactions between hydrogeomorphology and earthworm activity, and illustrates a mechanism of biogeomorphic inheritance through which ordered patterns of preferential flow can emerge in ecosystems. Further implications of hummock formation and nodal accumulation of nutrients are considered in relation to wetland resilience and regulatory ecosystem service provision.© 2018 John Wiley & Sons, Ltd.     

Understanding river dune splitting through flume experiments and analysis of a dune evolution model

Forecasts of water level during river floods require accurate predictions of the evolution of river dune dimensions, because the hydraulic roughness of the main channel is largely determined by the bed morphology. River dune dimensions are controlled by processes like merging and splitting of dunes. Particularly the process of dune splitting is still poorly understood and – as a result – not yet included in operational dune evolution models. In the current paper, the process of dune splitting is investigated by carrying out laboratory experiments and by means of a sensitivity analysis using a numerical dune evolution model. In the numerical model, we introduced superimposed TRIAS ripples (i.e. triangular asymmetric stoss side‐ripples) on the stoss sides of underlying dunes as soon as these stoss sides exceed a certain critical length. Simulations with the model including dune splitting showed that predictions of equilibrium dune characteristics were significantly improved compared to the model without dune splitting. As dune splitting is implemented in a parameterized way, the computational cost remains low which means that dune evolution can be calculated on the timescale of a flood wave. Subsequently, we used this model to study the mechanism of dune splitting. Literature showed that the initiation of a strong flow separation zone behind a superimposed bedform is one of the main mechanisms behind dune splitting. The flume experiments indicated that besides its height also the lee side slope of the superimposed bedform is an important factor to determine the strength of the flow separation zone and therefore is an important aspect in dune splitting. The sensitivity analysis of the dune evolution model showed that a minimum stoss side length was required to develop a strong flow separation zone. Copyright © 2014 John Wiley & Sons, Ltd.     

Shenandoah Watershed Study-Virginia Trout Stream Sensitivity Study (SWAS-VTSSS): Stream water quality and hydrologic monitoring data for mid-Appalachian headwater streams

The Shenandoah Watershed Study (established in 1979) and the Virginia Trout Stream Sensitivity Study (established in 1987) serve to increase understanding of hydrological and biogeochemical changes in western Virginia mountain streams that occur in response to acidic deposition and other ecosystem stressors. The SWAS-VTSSS program has evolved over its 40+ year history to consist of a temporally robust and spatially stratified monitoring framework. Currently stream water is sampled for water quality bi-hourly during high-flow events at three sites and weekly at four sites within Shenandoah National Park (SHEN), and quarterly at 72 sites and on an approximately decadal frequency at ~450 sites within the wider western Virginia Appalachian region. Stream water is evaluated for pH, acid neutralizing capacity (ANC), base cations (calcium, magnesium, sodium and potassium ion), acid anions (sulphate, nitrate and chloride), silica, ammonium, and conductivity with a subset of samples evaluated for monomeric aluminium and dissolved organic carbon. Hourly stream discharge (four sites) and in-situ measurements of conductivity, water and air temperature (three sites) are also measured within SHEN. Here we provide an overview and timeline of the SWAS-VTSSS stream water monitoring program, summarize the field and laboratory methods, describe the water chemistry and hydrologic data sets, and document major watershed disturbances that have occurred during the program history. Website links and instructions are provided to access the stream chemistry and time-series monitoring data in open-access federal databases. The purpose of this publication is to promote awareness of these unique, long-term data sets for wider use in catchment studies. The water chemistry and hydrologic data can be used to investigate a wide range of biogeochemical research questions and provide key inputs for models of these headwater stream ecosystems. SWAS-VTSSS is an ongoing program and quality assured data sets are uploaded to the databases annually.     

Setting nutrient thresholds to support an ecological assessment based on nutrient enrichment, potential primary production and undesirable disturbance

The EU Water Framework Directive recognises that ecological status is supported by the prevailing physico-chemical conditions in each water body. This paper describes an approach to providing guidance on setting thresholds for nutrients taking account of the biological response to nutrient enrichment evident in different types of water. Indices of pressure, state and impact are used to achieve a robust nutrient (nitrogen) threshold by considering each individual index relative to a defined standard, scale or threshold. These indices include winter nitrogen concentrations relative to a predetermined reference value; the potential of the waterbody to support phytoplankton growth (estimated as primary production); and detection of an undesirable disturbance (measured as dissolved oxygen). Proposed reference values are based on a combination of historical records, offshore (limited human influence) nutrient concentrations, literature values and modelled data. Statistical confidence is based on a number of attributes, including distance of confidence limits away from a reference threshold and how well the model is populated with real data. This evidence based approach ensures that nutrient thresholds are based on knowledge of real and measurable biological responses in transitional and coastal waters.     

Topographic monitoring of a middle estuary mudflat, Humber estuary, UK – Anthropogenic impacts and natural variation

Annual topographic surveys were carried out at the Saltend mudflat (Humber estuary, UK) between 1998 and 2006. These surveys formed part of an ongoing monitoring programme to examine the potential effects on the mudflat topography of the construction and operation of a waste water treatment works (WwTW) development by Yorkshire Water. Of particular concern was the potential disruption to the sedimentological regime within the special protection area (SPA) and candidate special area of conservation (cSAC) which could affect the invertebrate communities and ornithological functioning of the site. In addition to the development of the WwTW located to the extreme north-west of the site, a port extension removing 10 ha of the Saltend intertidal mudflat (outside the SPA but immediately south east of the WwTW) also occurred between 1999 and 2006. Minimal change was noted across the site following the construction and operation of the WwTW between 1998 and 2000. However, the construction of the bund in closer proximity to the SPA and cSAC masked any potential impact the WwTW could have had across the site after 2000. Profiles and contour mapping indicate that significant mudflat accretion occurred in the immediate area of the bund, with a general increase recorded across the western section of the site since 2000. In contrast the alternations to channel planform and subsequent rapid accretion of the mudflat to the east of the jetty, being a significant distance from the developments, are attributed to natural cyclical changes.     

Quantifying preferential flow in soils: A review of different techniques

Preferential flow (PF) in soil has both environmental and human health implications since it favours contaminant transport to groundwater without interaction with the chemically and biologically reactive upper layer of soil. PF is, however, difficult to measure and quantify. This paper reviews laboratory and field techniques, such as breakthrough curves, dye tracing, and scanning techniques, for evaluating PF in soil at different scales. Advanced technologies, such as scanning techniques, have increased our capability to quantify transport processes within the soil with minimal soil disturbance. Important issues with respect to quantifying PF concern large-scale studies, frozen soil conditions, tracing techniques for particles and gases, a lack of simple mathematical tools for interpreting field data, and the lack of a systematic approach for comparing PF data resulting from different measurement techniques. Also, more research is required to quantify the relative importance of the various PF processes that occur in soil rather than the integrated result of all PF processes in soils.     

The Low Metallicity ISM of Dwarf Galaxies

The properties of the low metallicity environments of dwarf galaxies are studied through dust observations in conjunction with the FIR fine-structure cooling lines. There is a striking enhancement of the I[CII]/I(CO) in dwarf galaxies that is explained by the decreased attenuation of the UV light in molecular clouds. An important consequence is that a significant mass of the molecular gas mass inventory can be missed through CO observations alone. Modeling the infrared spectral energy distribution into submillimeter wavelengths in dwarf galaxies reveals the presence of very cold (∼ 8K) dust,which accounts for a large fraction of the dust mass, until now missed by models using IRAS observations alone. In spite of the striking defficiency of the mid-infrared aromatic band carriers, cooling in the photodissociation regions, via [CII] line emission is a very efficient process. This revised version was published online in September 2006 with corrections to the Cover Date.     

Seasonal inundation patterns in two large savanna floodplains of South America: the Llanos de Moxos (Bolivia) and the Llanos del Orinoco (Venezuela and Colombia)

Inundation patterns in two of the largest savanna floodplains of South America were studied by analysis of the 37‐GHz polarization difference observed by the Scanning Multichannel Microwave Radiometer (Nimbus‐7 satellite). Flooded area was estimated at monthly intervals for January 1979 through to August 1987 using mixing models that account for the major landscape units with distinctive microwave emission characteristics. Results are presented separately for five subregions in each of the two floodplain regions to show the spatial as well as temporal variability in inundation patterns. The total area inundated during the 9 years varied between 2069 and 78 460 km² in the Llanos de Moxos (also spelled as Mojos; median area, 23 383 km²) and 1278 and 105 454 km² in the Llanos del Orinoco (median, 25 374 km²), not including the open‐water area of permanent lakes and river channels. The correlation between flooded area and river stage was used to extend the inundation records over a 30‐year period in the Moxos (1967–97) and a 58‐year period (1927–85) in the Orinoco. Interannual variability in inundation is greater in the Moxos than the Orinoco. Comparison of these data, however, with a previously published analysis of the Pantanal wetland shows that inundation patterns in these two floodplain regions are not as variable across years as they are in the Pantanal. Copyright © 2004 John Wiley & Sons, Ltd.     

Recent and fossil resins from New Zealand and Australia

The carbon-13 nuclear magnetic resonance spectra of fossil resins from New Zealand and Australia have been compared with those of modern and semifossilized materials. The great majority of the fossilized samples have strong spectral similarities to modern Agathis resins and to North American fossil resins, which have been attributed to Agathis. The Agathis-related spectra are different from those of modern Hymenaea and Araucaria. A small subgroup of Late Cretaceous resins from Australia and Papua New Guinea appears to derive from a different botanical source and shows strong resemblances to Claiborne amber from Arkansas. The spectral resonances of the exomethylene carbons degrade over time and on average provide an approximate measure of the geological age of Agathis-related fossil resins. © 1993 John Wiley & Sons, Inc.