A network theory approach for a better understanding of overland flow connectivity

Hydrological connectivity describes the physical coupling (linkages) of different elements within a landscape regarding (sub‐) surface flows. A firm understanding of hydrological connectivity is important for catchment management applications, for example, habitat and species protection, and for flood resistance and resilience improvement. Thinking about (geomorphological) systems as networks can lead to new insights, which has also been recognized within the scientific community, seeing the recent increase in the use of network (graph) theory within the geosciences. Network theory supports the analysis and understanding of complex systems by providing data structures for modelling objects and their linkages, and a versatile toolbox to quantitatively appraise network structure and properties. The objective of this study was to characterize and quantify overland flow connectivity dynamics on hillslopes in a humid sub‐Mediterranean environment by using a combination of high‐resolution digital‐terrain models, overland flow sensors and a network approach. Results showed that there are significant differences between overland flow connectivity on agricultural areas and semi‐natural shrubs areas. Significant positive correlations between connectivity and precipitation characteristics were found. Significant negative correlations between connectivity and soil moisture were found, most likely because of soil water repellency and/or soil surface crusting. The combination of structural networks and dynamic networks for determining potential connectivity and actual connectivity proved a powerful tool for analysing overland flow connectivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Vegetation and peat accumulation steer Holocene tidal–fluvial basin filling and overbank sedimentation along the Old Rhine River,The Netherlands

In the transformation from tidal systems to freshwater coastal landscapes, plants act as eco-engineering species that reduce hydrodynamics and trap sediment, but nature and timing of the mechanisms of land creation along estuaries remains unclear. This article focuses on the Old Rhine estuary (The Netherlands) to show the importance of vegetation in coastal landscape evolution, predominantly regarding tidal basin filling and overbank morphology. This estuary hosted the main outflow channel of the river Rhine between ca 6500 to 2000 cal bp , and was constrained by peat during most of its existence. This study reconstructs its geological evolution, by correlating newly integrated geological data and new field records to varying conditions. Numerical modelling was performed to test the inferred mechanisms. It was found that floodbasin vegetation and resulting organic accumulation strongly accelerated back-barrier infill, by minimizing tidal influence. After tidal and wave transport had already sufficiently filled the back-barrier basin, reed rapidly expanded from its edges under brackish conditions, as shown by diatom analysis and datings. Reed growth provided a positive infilling feedback by reducing tidal flow and tidal prism, accelerating basin infilling. New radiocarbon dates show that large-scale crevassing along the Old Rhine River – driven by tidal backwater effect – only started as nutrient-rich river water transformed the floodbasin into an Alder carr in a next phase of estuary evolution. Such less dense vegetation promotes crevassing as sediments are more easily transported into the floodbasin. As river discharge increased and estuary mouth infilling progressed, crevasse activity diminished around 3800 to 3000 cal bp , likely due to a reduced tidal backwater effect. The insights from this data-rich Holocene study showcase the dominant role that vegetation may have in the long-term evolution of coastal wetlands. It provides clues for effective use of vegetation in vulnerable wetland landscapes to steer sedimentation patterns to strategically adapt to rising water levels.     

Stratigraphic architecture of Solander Basin records Southern Ocean currents and subduction initiation beneath southwest New Zealand

Solander Basin is characterized by subduction initiation at the Pacific‐Australia plate boundary, where high biological productivity is found at the northern edge of the Antarctic Circumpolar Current. Sedimentary architecture results from tectonic influences on accommodation space, sediment supply and ocean currents (via physiography); and climate influence on ocean currents and biological productivity. We present the first seismic‐stratigraphic analysis of Solander Basin based on high‐fold seismic‐reflection data (voyage MGL1803, SISIE). Solander Trough physiography formed by Eocene rifting, but basinal strata are mostly younger than ca. 17 Ma, when we infer Puysegur Ridge formed and sheltered Solander Basin from bottom currents, and mountain growth onshore increased sediment supply. Initial inversion on the Tauru Fault started at ca. 15 Ma, but reverse faulting from 12 to ca. 8 Ma on both the Tauru and Parara Faults was likely associated with reorganization and formation of the subduction thrust. The new seabed topography forced sediment pathways to become channelized at low points or antecedent gorges. Since 5 Ma, southern Puysegur Ridge and Fiordland mountains spread out towards the east and Solander Anticline grew in response to ongoing subduction and growth of a slab. Solander Basin had high sedimentation rates because (1) it is sheltered from bottom currents by Puysegur Ridge; and (2) it has a mountainous land area that supplies sediment to its northern end. Sedimentary architecture is asymmetric due to the Subtropical Front, which moves pelagic and hemi‐pelagic sediment, including dilute parts of gravity flows, eastward and accretes contourites to the shelf south of Stewart Island. Levees, scours, drifts and ridges of folded sediment characterize western Solander Basin, whereas hemi‐pelagic drape and secondary gravity flows are found east of the meandering axial Solander Channel. The high‐resolution record of climate and tectonics that Solander Basin contains may yield excellent sites for future scientific ocean drilling.     

The unusual and large drawdown response of buried-valley aquifers to pumping

The buried-valley aquifers that are common in the glacial deposits of the northern hemisphere are a typical case of the strip aquifers that occur in many parts of the world. Pumping from a narrow strip aquifer leads to much greater drawdown and much more distant drawdown effects then would occur in a sheet aquifer with a similar transmissivity and storage coefficient. Widely used theories for radial flow to wells, such as the Theis equation, are not appropriate for narrow strip aquifers. Previously published theory for flow to wells in semiconfined strip aquifers is reviewed and a practical format of the type curves for pumping-test analysis is described. The drawdown response of strip aquifers to pumping tests is distinctive, especially for observation wells near the pumped well. A case study is presented, based on extensive pumping test experience for the Estevan Valley Aquifer in southern Saskatchewan, Canada. Evaluation of groundwater resources in such buried-valley aquifers needs to take into account the unusually large drawdowns in response to pumping.     

Soil Organic Carbon mapping of partially vegetated agricultural fields with imaging spectroscopy

Soil Organic Carbon (SOC) is one of the key soil properties, but the large spatial variation makes continuous mapping a complex task. Imaging spectroscopy has proven to be an useful technique for mapping of soil properties, but the applicability decreases rapidly when fields are partially covered with vegetation. In this paper we show that with only a few percent fractional maize cover the accuracy of a Partial Least Square Regression (PLSR) based SOC prediction model drops dramatically. However, this problem can be solved with the use of spectral unmixing techniques. First, the fractional maize cover is determined with linear spectral unmixing, taking the illumination and observation angles into account. In a next step the influence of maize is filtered out from the spectral signal by a new procedure termed Residual Spectral Unmixing (RSU). The residual soil spectra resulting from this procedure are used for mapping of SOC using PLSR, which could be done with accuracies comparable to studies performed on bare soil surfaces (Root Mean Standard Error of Calibration = 1.34 g/kg and Root Mean Standard Error of Prediction = 1.65 g/kg). With the presented RSU approach it is possible to filter out the influence of maize from the mixed spectra, and the residual soil spectra contain enough information for mapping of the SOC distribution within agricultural fields. This can improve the applicability of airborne imaging spectroscopy for soil studies in temperate climates, since the use of the RSU approach can extend the flight-window which is often constrained by the presence of vegetation.     

Critical soil conditions for oxygen stress to plant roots: Substituting the Feddes-function by a process-based model

Effects of insufficient soil aeration on the functioning of plants form an important field of research. A well-known and frequently used utility to express oxygen stress experienced by plants is the Feddes-function. This function reduces root water uptake linearly between two constant pressure heads, representing threshold values for minimum and maximum oxygen deficiency. However, the correctness of this expression has never been evaluated and constant critical values for oxygen stress are likely to be inappropriate. On theoretical grounds it is expected that oxygen stress depends on various abiotic and biotic factors. In this paper, we propose a fundamentally different approach to assess oxygen stress: we built a plant physiological and soil physical process-based model to calculate the minimum gas filled porosity of the soil (_{gas_min}) at which oxygen stress occurs.First, we calculated the minimum oxygen concentration in the gas phase of the soil needed to sustain the roots through (micro-scale) diffusion with just enough oxygen to respire. Subsequently, _{gas_min} that corresponds to this minimum oxygen concentration was calculated from diffusion from the atmosphere through the soil (macro-scale).We analyzed the validity of constant critical values to represent oxygen stress in terms of _{gas_min}, based on model simulations in which we distinguished different soil types and in which we varied temperature, organic matter content, soil depth and plant characteristics. Furthermore, in order to compare our model results with the Feddes-function, we linked root oxygen stress to root water uptake (through the sink term variable F, which is the ratio of actual and potential uptake).The simulations showed that _{gas_min} is especially sensitive to soil temperature, plant characteristics (root dry weight and maintenance respiration coefficient) and soil depth but hardly to soil organic matter content. Moreover, _{gas_min} varied considerably between soil types and was larger in sandy soils than in clayey soils. We demonstrated that F of the Feddes-function indeed decreases approximately linearly, but that actual oxygen stress already starts at drier conditions than according to the Feddes-function. How much drier is depended on the factors indicated above. Thus, the Feddes-function might cause large errors in the prediction of transpiration reduction and growth reduction through oxygen stress.We made our method easily accessible to others by implementing it in SWAP, a user-friendly soil water model that is coupled to plant growth. Since constant values for _{gas_min} in plant and hydrological modeling appeared to be inappropriate, an integrated approach, including both physiological and physical processes, should be used instead. Therefore, we advocate using our method in all situations where oxygen stress could occur.