A hydropedological approach to simulate streamflow and soil water contents with SWAT+

Reflecting internal catchment hydrological processes in hydrological models is important for accurate predictions of the impact of climate and land-use change on water resources. Characterizing these processes is however difficult and expensive due to their dynamic nature and spatio-temporal variability. Hydropedology is a relatively new discipline focusing on the synergistic integration of hydrology, soil physics and pedology. Hydropedological interpretations of soils and soil distribution can be used to characterize key hydrological processes, especially in areas with no or limited hydrometric measurements. Here we applied a hydropedological approach to reflect flowpaths through detailed routing in SWAT+ for a 157 ha catchment (Weatherley) in South Africa. We compared the hydropedological approach and a standard (no routing) approach against measured streamflow (two weirs) and soil water contents (13 locations). The catchment was treated as ‘ungauged’ and the model was not calibrated against hydrometric measurements in order to determine the direct contribution of hydropedology on modelling efficiency. Streamflow was predicted well without calibration (NSE > 0.8; R² > 0.82) for both approaches at both weirs. The standard approach yielded slightly better streamflow predictions. The hydropedological approach resulted in considerable improvements in the simulation of soil water contents (R² increased from 0.40 to 0.49 and PBIAS decreased from 40% to 20%). The routing capacity of SWAT+ as employed in the hydropedological approach reduced the underestimation of wetland water regimes drastically and resulted in a more accurate representation of the dominant hydrological processes in this catchment. We concluded that hydropedology can be a valuable source of ‘soft data’ to reflect internal catchment structure and processes and, potentially, for realistic calibrations in other studies, especially those conducted in areas with limited hydrometric measurements.     

Spatiotemporal variability in hydrochemistry of shallow groundwater in a small pre‐alpine catchment: The importance of landscape elements

Topography and landscape characteristics affect the storage and release of water and, thus, groundwater dynamics and chemistry. Quantification of catchment scale variability in groundwater chemistry and groundwater dynamics may therefore help to delineate different groundwater types and improve our understanding of which parts of the catchment contribute to streamflow. We sampled shallow groundwater from 34 to 47 wells and streamflow at seven locations in a 20‐ha steep mountainous catchment in the Swiss pre‐Alps, during nine baseflow snapshot campaigns. The spatial variability in electrical conductivity, stable water isotopic composition, and major and trace ion concentrations was large and for almost all parameters larger than the temporal variability. Concentrations of copper, zinc, and lead were highest at sites that were relatively dry, whereas concentrations of manganese and iron were highest at sites that had persistent shallow groundwater levels. The major cation and anion concentrations were only weakly correlated to individual topographic or hydrodynamic characteristics. However, we could distinguish four shallow groundwater types based on differences from the catchment average concentrations: riparian zone‐like groundwater, hillslopes and areas with small upslope contributing areas, deeper groundwater, and sites characterized by high magnesium and sulfate concentrations that likely reflect different bedrock material. Baseflow was not an equal mixture of the different groundwater types. For the majority of the campaigns, baseflow chemistry most strongly resembled riparian‐like groundwater for all but one subcatchment. However, the similarity to the hillslope‐type groundwater was larger shortly after snowmelt, reflecting differences in hydrologic connectivity. We expect that similar groundwater types can be found in other catchments with steep hillslopes and wet areas with shallow groundwater levels and recommend sampling of groundwater from all landscape elements to understand groundwater chemistry and groundwater contributions to streamflow.     

Preservation of Re‐Os isotope signatures in pyrite throughout low‐T,high‐P eclogite facies metamorphism

Pyrite in LT–HP eclogites from the western Tianshan orogenic belt yields a Re‐Os age of 378.1 ± 8.9 Ma, which is 30–70 Ma older than ages previously obtained for the same rocks using the Rb–Sr, Sm–Nd, Ar–Ar, U–Pb, and Lu–Hf isotope systems. The Tianshan LT–HP eclogite experienced temperatures of up to ~570 °C combined with pressures of up to 2.1 GPa during metamorphism. These conditions are below the transition of pyrite to pyrrhotite, which defines both pyrite stability and possibly its closure temperature for Re‐Os. Pyrite can preserve Re‐Os signatures through eclogite facies peak metamorphic conditions, and thus allow determination of the formation age of pyrite in the protolith.     

Mantle resistance against Gibraltar slab dragging as a key cause of the Messinian Salinity Crisis

The Messinian Salinity Crisis (5.97–5.33 Ma) was caused by the closure of the Atlantic‐Mediterranean gateways that cut through the Gibraltar orogenic system. The geodynamic drivers underlying gateway closure and re‐opening are still debated. Here, we interrogate the gateway successions to find the imprints of surface deformation, infer the timing and nature of associated geodynamic drivers, and test such inferences against numerical simulations of slab dynamics. We find that since the latest Miocene, a tectonic framework was established in the gateway region dominated simultaneously by (a) relative plate convergence, (b) slab tearing under the eastern Betic Cordillera and (c) mantle resistance against north‐northeastward dragging of the Gibraltar slab by the African plate's absolute motion. We propose that mantle‐resisted slab dragging and slab tearing operated in concert closing the gateways that caused the Messinian Salinity Crisis, whereas sinking of heavy oceanic lithosphere located between buoyant continental plates re‐opened the Strait of Gibraltar at 5.33 Ma.     

Exploring hydroclimatic change disparity via the Budyko framework

The Budyko framework characterizes landscape water cycles as a function of climate. We used this framework to identify regions with contrasting hydroclimatic change during the past 50 years in Sweden. This analysis revealed three distinct regions: the mountains, the forests, and the areas with agriculture. Each region responded markedly different to recent climate and anthropogenic changes, and within each region, we identified the most sensitive subregions. These results highlight the need for regional differentiation in climate change adaptation strategies to protect vulnerable ecosystems and freshwater resources. Further, the Budyko curve moved systematically towards its water and energy limits, indicating augmentation of the water cycle driven by changing vegetation, climate and human interactions. This finding challenges the steady state assumption of the Budyko curve and therefore its ability to predict future water cycles. Copyright © 2013 John Wiley & Sons, Ltd.     

The Late Cretaceous Klepa basalts in Macedonia (FYROM)—Constraints on the final stage of Tethys closure in the Balkans

The waning stage(s) of the Tethyan ocean(s) in the Balkans are not well understood. Controversy centres on the origin and life‐span of the Cretaceous Sava Zone, which is allegedly a remnant of the last oceanic domain in the Balkan Peninsula, defining the youngest suture between Eurasia‐ and Adria‐derived plates. In order to investigate to what extent Late‐Cretaceous volcanism within the Sava Zone is consistent with this model we present new age data together with trace‐element and Sr–Nd–Pb isotope data for the Klepa basaltic lavas from the central Balkan Peninsula. Our new geochemical data show marked differences between the Cretaceous Klepa basalts (Sava Zone) and the rocks of other volcanic sequences from the Jurassic ophiolites of the Balkans. The Klepa basalts mostly have Sr–Nd–Pb isotopic and trace‐element signatures that resemble enriched within‐plate basalts substantially different from Jurassic ophiolite basalts with MORB, BAB and IAV affinities. Trace‐element modelling of the Klepa rocks indicates 2%–20% polybaric melting of a relatively homogeneously metasomatised mantle source that ranges in composition from garnet lherzolite to ilmenite+apatite bearing spinel–amphibole lherzolite. Thus, the residual mineralogy is characteristic of a continental rather than oceanic lithospheric mantle source, suggesting an intracontinental within‐plate origin for the Klepa basalts. Two alternative geodynamic models are internally consistent with our new findings: (1) if the Sava Zone represents remnants of the youngest Neotethyan Ocean, magmatism along this zone would be situated within the forearc region and triggered by ridge subduction; (2) if the Sava Zone delimits a diffuse tectonic boundary between Adria and Europe which had already collided in the Late Jurassic, the Klepa basalts together with a number of other magmatic centres represent volcanism related to transtensional tectonics.     

The Dynamics of Small Postglacial Lake Basins and the Nature of Their Archaeological Record: A Case Study of the Middle Palaeolithic Site Neumark‐Nord 2, Germany

Due to favorable conditions of preservation, sedimentary basins provide rich records of human behavior and its environmental context. The conditions for the preservation of archaeological material however vary between basin types (large, river‐fed or small, closed basins), while conditions also differ within a particular basin environment. The goal of this paper is to understand how the dynamics of a small postglacial basin such as Neumark‐Nord 2, a context that dominates the Eemian archaeological record, affected the archaeology situated at its basin margin. The approach used here is to correlate the archaeological record with reconstructions of patterns of deposition and the water conditions within the basin, using lithology, micromorphology, pollen, and macroremains from a transect running from the basin center to the margin. The results show that (1) find levels were exposed to overland flow‐induced winnowing, which vertically concentrated finds but did not cause significant transport, (2) find levels correspond to phases of increased water presence in the basin, and (3) lateral shifts in hominin activity areas may reflect adjustments in the water level. The research shows the importance of large‐scale archaeological excavations and a multidisciplinary sampling strategy that covers both the basin center and the margins, when studying postglacial basin localities like Neumark‐Nord 2.     

Autogenic avulsion,channelization and backfilling dynamics of debris‐flow fans

Alluvial fans develop their semi‐conical shape by quasi‐cyclic avulsions of their geomorphologically active sector from a fixed fan apex. On debris‐flow fans, these quasi‐cyclic avulsions are poorly understood, partly because physical scale experiments on the formation of fans have been limited largely to turbidite and fluvial fans and deltas. In this study, debris‐flow fans were experimentally created under constant extrinsic forcing, and autogenic sequences of backfilling, avulsion and channelization were observed. Backfilling, avulsion and channelization were gradual processes that required multiple successive debris‐flow events. Debris flows avulsed along preferential flow paths given by the balance between steepest descent and flow inertia. In the channelization phase, debris flows became progressively longer and narrower because momentum increasingly focused on the flow front as flow narrowed, resulting in longer run‐out and deeper channels. Backfilling commenced when debris flows reached their maximum possible length and channel depth, as defined by channel slope and debris‐flow volume and composition, after which they progressively shortened and widened until the entire channel was filled and avulsion was initiated. The terminus of deposition moved upstream because the frontal lobe deposits of previous debris flows created a low‐gradient zone forcing deposition. Consequently, the next debris flow was shorter which led to more in‐channel sedimentation, causing more overbank flow in the next debris flow and resulting in reduced momentum to the flow front and shorter runout. This topographic feedback is similar to the interaction between flow and mouth bars forcing backfilling and transitions from channelized to sheet flow in turbidite and fluvial fans and deltas. Debris‐flow avulsion cycles are governed by the same large‐scale topographic compensation that drives avulsion cycles on fluvial and turbidite fans, although the detailed processes are unique to debris‐flow fans. This novel result provides a basis for modelling of debris‐flow fans with applications in hazards and stratigraphy.