The contrasted evolution of high and low flows and precipitation indices in the Duero basin (Spain)

Abstract

Trends in high and low flows are valuable indicators of hydrological change because they highlight changes in various parts of the frequency distribution of streamflow series. This enables improved assessment of water availability in regions with high seasonal and inter-annual variability. There has been a substantial reduction in water resources in the Duero basin (Iberian Peninsula, Spain) and other areas of the Mediterranean region during the last 50 years, and this is likely to continue because of climate change. In this study, we investigated the evolution and trends in high and low flows in the Spanish part of the Duero basin, and in equivalent or closely-related precipitation indices for the period 1961–2005. The results showed a general trend of decrease in the frequency and magnitude of high flows throughout most of the basin. Moreover, the number of days with low flows significantly increased over this period. No clear relationship was evident between the evolution of high/low flows and changes in the distribution frequencies of the precipitation series. In contrast to what was expected, the number of days with heavy precipitation and the mean annual precipitation did not show significant trends across the basin, and the number of days without rainfall decreased slightly. The divergence between precipitation and runoff evolution was more accentuated in spring and summer. In the absence of trends in precipitation, it is possible that reforestation processes in the region, and increasing temperatures in recent decades, could be related to the decreasing frequency of high flows and the increasing frequency of low flows.

Editor Z.W. Kundzewicz; Associate editor S. Grimaldi

Citation Morán-Tejeda, E., López-Moreno, J.I., Vicente-Serrano, S.M., Lorenzo-Lacruz, J. and Ceballos-Barbancho, A., 2012. The contrasted evolution of high and low flows and precipitation indices in the Duero basin (Spain). Hydrological Sciences Journal, 57 (4), 591–611.     

Validation of ERS scatterometer‐derived soil moisture data in the central part of the Duero Basin,Spain

The objective of this study was to validate the soil moisture data derived from coarse‐resolution active microwave data (50 km) from the ERS scatterometer. The retrieval technique is based on a change detection method coupled with a data‐based modelling approach to account for seasonal vegetation dynamics. The technique is able to derive information about the soil moisture content corresponding to the degree of saturation of the topmost soil layer (∼5 cm). To estimate profile soil moisture contents down to 100 cm depth from the scatterometer data, a simple two‐layer water balance model is used, which generates a red noise‐like soil moisture spectrum. The retrieval technique had been successfully applied in the Ukraine in a previous study. In this paper, the performance of the model in a semi‐arid Mediterranean environment characterized by low annual precipitation (400 mm), hot dry summers and sandy soils is investigated. To this end, field measurements from the REMEDHUS soil moisture station network in the semi‐arid parts of the Duero Basin (Spain) were used. The results reveal a significant coefficient of determination (R² = 0·75) for the averaged 0–100 cm soil moisture profile and a root mean square error (RMSE) of 2·2 vol%. The spatial arrangement of the REMEDHUS soil moisture stations also allowed us to study the influence of the small‐scale variability of soil moisture within the ERS scatterometer footprint. The results show that the small‐scale variability in the study area is modest and can be explained in terms of texture fraction distribution in the soil profiles. Copyright © 2004 John Wiley & Sons, Ltd.     

Naturally occurring arsenic in groundwater and identification of the geochemical sources in the Duero Cenozoic Basin, Spain

Arsenic concentrations surpassing potability limit of 10 μg/L in the groundwater supplies of an extensive area in the Duero Cenozoic Basin (central Spain) have been detected and the main sources of arsenic identified. Arsenic in 514 samples of groundwater, having mean values of 40.8 μg/L, is natural in origin. Geochemical analysis of 553 rock samples, assaying arsenic mean values of 23 mg/kg, was performed. Spatial coincidence between the arsenic anomaly in groundwater and the arsenic lithogeochemical distribution recorded in the Middle Miocene clayey organic-rich Zaratan facies illustrates that the rocks of this unit are the main source of arsenic in groundwater. The ferricretes associated to the Late Cretaceous–Middle Miocene siliciclastics also constitute a potential arsenic source. Mineralogical study has identified the presence of arsenic in iron oxides, authigenic pyrite, manganese oxides, inherited titanium–iron oxides, phyllosilicates and organomineral compounds. Arsenic mobilization to groundwater corresponds to arsenic desorption from iron and manganese oxides and from organic matter.     

Evaluating Duero's basin (Spain) response to the NAO phases: spatial and seasonal variability

In this study, we analysed the influence of the North Atlantic Oscillation (NAO) on the hydrological response of the Duero River (central Spain) and its tributaries. We identified the positive and negative phases of the winter NAO for the period 1961–2006 and related precipitation and river discharge anomalies. Significant differences in precipitation and river discharge were found between the positive and negative NAO phases with negative anomalies (dry conditions) during positive NAO periods, and positive anomalies (wet conditions) during negative NAO periods. Marked differences were also found in the precipitation and river discharge response times to the NAO episodes, with an immediate response for precipitation but a lag and more sustained response for river discharges. Differing patterns were also identified in the response of river discharge to the NAO throughout the Duero basin. The physical characteristics of watersheds (including area, altitude and permeability explained most of the differences in the timing and magnitude of anomalies in the river discharge in response to the NAO. The findings highlight the great variability in the hydrological response of rivers to the NAO episodes. Copyright © 2010 John Wiley & Sons, Ltd.     

Origin and evolution of playas and blowouts in the semiarid zone of Tierra de Pinares (Duero Basin, Spain)

A large number of blowouts and playas occur in the marginal sectors of the aeolian deposits located in the southern sector of the Duero Depression (Tierra de Pinares) in Spain. The blowouts are relict landforms that were developed on sand sheets by deflation during dry periods with lower vegetation cover and a deeper water table. The studied blowouts form complexes of NW–SE and NNW–SSE elongated hollows with accompanying dunes up to 4 km long in the leeward margins. Some hollows host lakes or swampy areas related to a shallow water table. The dunes formed by NE–ENE winds show steep windward slopes and gentle leeward slopes. The studied playas, with prevalent NNW–SSE orientations, result from the aeolian excavation of terrace deposits and the underlying marly bedrock. It is probable that the formation of these depressions in an initial stage was related to deflation processes affected preferentially NNW–SSE sandy channels perpendicular to the dominant wind direction. The precipitation of salts in the playas generates aggregates of clay particles (peloids) that are easily removed by the wind. Once the bottom reached the substratum, the deepening of the depressions progressed by the deflation of particles produced by weathering of the argillaceous bedrock.     

Erosion surfaces and Neogene landscape evolution in the NE Duero Basin (north-central Spain)

We present a chronological model of erosion surface development in the Iberian and Cantabrian Ranges of north-central Spain. We map four erosion surfaces and interpret these to be related to Duero basin continental sediments and tectonic activity from Upper Oligocene to Plio-Pleistocene. The oldest erosion surface, SE1, formed across Upper Oligocene–Lower Miocene synorogenic deposits; while surface SE2 was contemporaneous with the Middle Miocene alluvial systems, ending with an uplift stage in the Astaracian. The two most recent erosion surfaces, SE3 and SE4, developed during extensional tectonic episodes and are associated with the deposition of Upper Páramo sedimentary units at the end of the Miocene (Upper Turolian) and alluvial fan deposits, known as rañas (Plio-Pleistocene). With the exception of SE1, which seems to be associated with a relatively wet climate, the surfaces formed during periods of marked aridity and generally warm temperatures. Through geostatistical reconstruction of the best preserved surface (SE2), applying ordinary kriging method to the topography (DEM) of the erosion surface and its correlating sedimentary plains, we identify the deformation processes which occurred on this surface after its formation.     

A GIS method to identify flat surfaces and restore relict fluvial long-profiles from terrace remnants gives new clues on how large basins respond to endorheic–exorheic transitions (Duero basin,Iberian Peninsula)

Fluvial terraces are used as geomorphic indicators for deciphering long-term landscape evolution. Knowing the distribution of fluvial terraces is essential for establishing former river profiles and their tectonic significance, for studying climate-modulated processes of terrace development, or for defining fluvial network adjustments in response to sudden base-level changes like those produced by fluvial captures. Multiple methods for automatic map production have been proposed based on the comparison of morphometric indices with those of the modern river course. Here we propose an alternative method to identify flat surfaces and scarps separating them from digital elevation models without setting comparisons with a modern river course and thus fully applicable to study flat landforms whatever their origin. Its application to the low-relief landscape of the Cenozoic Duero basin has allowed the improvement of previous geomorphological maps and the analysis of fluvial network adjustments in response to a sudden base-level fall after the opening of the Neogene endorheic basin towards the Atlantic Ocean. Reconstructed terrace long-profiles suggest an initial episode of fast vertical incision followed by a period of repeated planation–aggradation–incision with the formation of 14 to 13 unpaired terrace levels. Changes observed in the pattern of terrace profiles are discussed with regard to changes in regional tectonics and base-level variations. © 2019 John Wiley & Sons, Ltd.     

Soil-water behaviour of sandy soils under semi-arid conditions in the Duero Basin (Spain)

This paper examines the soil-water behaviour of sandy soils representative of a broad area in the Duero river basin (Spain), with semi-arid bioclimatic characteristics. Soil-moisture measurements of the first 25 cm of soil profile were taken using time domain reflectometry (TDR) methodology. The infiltration process was studied by means of field experiments combining rainfall simulations and TDR techniques. Finally, a set of hydrodynamic variables was calculated in the laboratory (i.e. available soil-water content, soil-water-retention curves and hydraulic conductivity).The average available soil-water content is 0·07 cm³ cm⁻³, in accordance with the predominance of the sandy textural fraction (85%). The rainfall simulation experiments show that infiltration prevails over overland flow, with a steady-state infiltration rate higher than rainfall application (40mmh⁻¹) in most cases; a low water-retention capacity was also observed. The evolution of soil-moisture over time indicates a close dependence on intra-annual rainfall variability due to the low water-retention capacity of sandy soils. With a maximum water-storage capacity of around 0·18 cm³ cm⁻³, the soils studied show high sensitivity to dry spells and may remain with <50% of the available soil-water content for several months during the year.     

Nocturnal basin low-level jets: an integrated study

Low-level jets (LLJs) are a very common feature in the nocturnal stably stratified boundary layer. Many factors can intervene in their generation, linked basically to effects of baroclinity. A special kind of low-level jets is composed by the nocturnal katabatic and basin flows, generated over terrain slopes. A study of observed LLJs in the Duero Basin is shown here, combining observational data and modelling experiments. Normalized in respect to the maximum wind height, the dynamic characteristics of the jets are similar: a two-layer system, with a stably stratified layer below the jet maximum and a near neutral layer above, with a very stable layer separating them at the level of the wind maximum. There is vertical mixing above and below the jet, and the connection between these layers takes place occasionally in a very turbulent manner.     

Quantitative reconstruction of Late Cenozoic landscapes: a case study in the Sierra de Atapuerca (Burgos,Spain)

We have developed a method to reconstruct palaeorelief by means of detailed geomorphological and geological studies, geostatistical tools, GIS and a DEM. This method has been applied to the Sierra de Atapuerca (NE Duero Basin, Burgos, Spain), allowing us to model a three‐dimensional reconstruction of the relief evolution from the Middle Miocene to the present. The modelling procedure is based on geostatistical recovery of the palaeosurfaces characteristic of each geomorphological evolution stage, using polynomial regressions, trend surfaces and kriging. The modelling of morphology trends has been useful in establishing new geological and geomorphological relationships in the geodynamic evolution of this basin, such as uplift quantification, correlation of erosion surfaces and sedimentary units, and the evolution of fluvial base levels. The palaeosurface reconstruction together with an analysis of the slope retreat have allowed us to reconstruct the palaeoreliefs that define the Late Cenozoic landscape evolution of this area, where the Lower and Middle Pleistocene archaeopalaeontological sites of the Sierra de Atapuerca are located. Copyright © 2007 John Wiley & Sons, Ltd.