Temporal and spatial variation of soil erosion in wooded rangelands of southwest Spain

Changes in land use are common in Mediterranean areas and are reported as having produced changes in the intensity of soil erosion. Dehesas are rangelands with a disperse tree cover, widespread in the south-western part of the Iberian Peninsula and similar ecosystems are also common in other areas with a Mediterranean climate. The aim of the present study is to analyse temporal and spatial variations of soil erosion rates estimated along three hillsides, located in two farms (Buitrera and Parapuños) in southwest Spain. To understand the temporal variation, soil erosion rates were studied in light of land use-management changes that took place during the last few centuries. Results indicate very low erosion rates prior to the 18^th century in both farms. In Buitrera, a first increase of soil loss rates was identified during the period 1831-1897, amounting to 7.4 t ha^-1 y^-1. A further increase took place during the 20^th century, reaching a mean erosion rate of 29.1 t ha^-1 y^-1. In Parapuños, data points to a significant increase from 1881 onwards, with an estimated mean erosion rate of 18.5 t ha^-1 y^-1. Those increases were presumably connected with an intensification of land use, such as cultivation and excessive livestock populations. Regarding spatial variation, the bare surface and the erosive power of run-off along the hillsides accounts for 76% of the soil erosion rates dispersion. At a local scale, the variability of erosion rates could not be explained, because of (i) uncertainty related to the micromorphology of the past soil surface and (ii) the role of tillage erosion in the past. However, the results obtained offer valuable data on the temporal and spatial variation of erosion rates in dehesas at the hillslope scale and a similar approach could be used for other rangelands with a disperse tree cover. © 2019 John Wiley & Sons, Ltd.     

Biocrust landscape-scale spatial distribution is strongly controlled by terrain attributes: Topographic thresholds for colonization in a semiarid badland system

Biological soil crust, or biocrust communities, are the dominating life form in many extreme habitats, such as arid and semiarid badlands, where water scarcity and highly erodible substrates limit vegetation cover. While climate, soil and biotic factors have been described as environmental filters influencing biocrust distribution in such biomes, little is known about the effect of terrain attributes on creating specific microhabitats that promote or restrict biocrust colonization. This study aimed to identify the main terrain attributes controlling biocrust distribution in the driest badland system in Europe, the Tabernas Badlands (SE Spain). To do this, we analysed the influence of different terrain attributes related to landscape stability and microclimate formation on the spatial distribution of lichen and cyanobacteria, using field measurements and topographical information from a LiDAR survey. Our results showed that the spatial distribution of cyanobacteria-dominated biocrusts, which are physiologically and morphologically adapted to extreme drought and high UVA radiation, was mostly associated with areas of high potential incoming solar radiation. The exception was bare south-aspect hillslopes with very high sediment transport potential, where bare physically crusted soils were the dominant ground cover. Lichen-dominated biocrusts, in contrast, colonized near the top of north-aspect hillslopes, characterized by low potential incoming solar radiation and potential evapotranspiration, and their cover decreased downstream, as conditions became good enough for vascular plants. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Hydrogeophysical modelling of Hisarcik (Kütahya) geothermal field,western Turkey

Western Anatolia hosts many low-to-moderate and high-temperature geothermal sources in which active faults play a dominant role to control the recharge and the discharge of geothermal fluid. In this study, we used the two-dimensional geoelectric structure of Kütahya Hisarcık geothermal field, and created a conceptual hydrogeophysical model that includes faults, real topographical variations and geological units. The temperature distribution and fluid flow pattern are also investigated. The depth extension of Hisarcık Fault, electrical basement and low resistivity anomalies related to the presence of geothermal fluid are determined by using resistivity studies in the area. Numerical simulations suggest that Hisarcık fault functioning as a fluid conduit primarily enables hot fluid to be transported from depth to the surface. It is shown that the locations of predicted outflow vents coincide with those of hot springs in the area.     

Changes in phytoplankton production and upwelling intensity off A Coruña (NW Spain) for the last 28 years

Ocean Dynamics - The analysis of a 28-year-long (1989–2016) series of monthly measurements of chlorophyll concentrations and primary production rates at a shelf station off A Coruña (NW...     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 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.     

Paleoflood events recorded by speleothems in caves

Speleothems are usually composed of thin layers of calcite (or aragonite). However, cemented detrital materials interlayered between laminae of speleothemic carbonate have been also observed in many caves. Flowstones comprising discontinuous carbonate layers form due to flowing water films, while flood events introduce fluviokarstic sediments in caves that, on occasion, are recorded as clayey layers inside flowstones and stalagmites. This record provides a potential means of understanding the frequency of palaeofloods using cave records. In this work, we investigate the origin of this type of detrital deposit in El Soplao Cave (Northern Spain). The age of the lowest aragonite layer of a flowstone reveals that the earliest flood period occurred before 500 ka, though most of the flowstone formed between 422 +69/‐43 ka and 400 +66/‐42 ka. This suggests that the cave was periodically affected by palaeoflood events that introduced detrital sediments from the surface as a result of occasional extreme rainfall events, especially at around 400 ka. The mineralogical data enable an evolutionary model for this flowstone to be generated based on the alternation of flood events with laminar flows and carbonate layers precipitation that can be extrapolated to other caves in which detrital sediments inside speleothems have been found. Copyright © 2014 John Wiley & Sons, Ltd.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

Piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events on the Loess Plateau,China

Establishing a universal watershed‐scale erosion and sediment yield prediction model represents a frontier field in erosion and soil/water conservation. The research presented here was conducted on the Chabagou watershed, which is located in the first sub‐region of the hill‐gully area of the Loess Plateau, China. A back‐propagation artificial neural model for watershed‐scale erosion and sediment yield was established, with the accuracy of the model, then compared with that of multiple linear regression. The sensitivity degree of various factors to erosion and sediment yield was quantitatively analysed using the default factor test. On the basis of the sensitive factors and the fractal information dimension, the piecewise prediction model for erosion and sediment yield of individual rainfall events was established and further verified. The results revealed the back‐propagation artificial neural network model to perform better than the multiple linear regression model in terms of predicting the erosion modulus, with the former able to effectively characterize dynamic changes in sediment yield under comprehensive factor conditions. The sensitivity of runoff erosion power and runoff depth to the erosion and sediment yield associated with individual rainfall events was found to be related to the complexity of surface topography. The characteristics of such a hydrological response are thus closely related to topography. When the fractal information dimension is greater than the topographic threshold, the accuracy of prediction using runoff erosion power is higher than that of using runoff depth. In contrast, when the fractal information dimension is smaller than the topographic threshold, the accuracy of prediction using runoff depth is higher than that of using runoff erosion power. The developed piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events, which introduces runoff erosion power and runoff depth using the fractal information dimension as a boundary, can be considered feasible and reliable and has a high prediction accuracy. Copyright © 2013 John Wiley & Sons, Ltd.