Land subsidence resulting from groundwater extractions is a global phenomenon adversely affecting many regions worldwide. Understanding the governing processes and mitigating associated hazards require knowing the spatial distribution of the implicated factors (piezometric levels, lithology, ground deformation), usually only known at discrete locations. Here, we propose a methodology based on the Kriging with External Drift (KED) approach to interpolate sparse point measurements of variables influencing land subsidence using high density InSAR measurements. In our study, located in the Alto Guadalentín basin, SE Spain, these variables are GPS vertical velocities and the thickness of compressible soils. First, we estimate InSAR and GPS rates of subsidence covering the periods 2003–2010 and 2004–2013, respectively. Then, we apply the KED method to the discrete variables. The resulting continuous GPS velocity map shows maximum subsidence rates of 13 cm/year in the center of the basin, in agreement with previous studies. The compressible deposits thickness map is significantly improved. We also test the coherence of Sentinel-1 data in the study region and evaluate the applicability of this methodology with the new satellite, which will improve the monitoring of aquifer-related subsidence and the mapping of variables governing this phenomenon.

Ephemeral streams induce flash-flood events, which cause dramatic morphological changes and impacts on population, mainly because they are intermittent and less predictable. Human pressures on the basin modify load and discharge relationships, inducing dormant instability on the fluvial system that will manifest abruptly during flood events. The flash-flood response of two ephemeral streams affected by load supply modification due to land use changes is discussed in a combination of geomorphic and hydraulic approaches. During the Rivillas flash flood, intensive clearing on the basin led to high rates of sediment flowing into an artificially straightened and inefficient channel. The stream evolved from a sinuous single channel into a shallow braiding occupying the entire width of the valley floor. Misfits and unsteady channel conditions increased velocity, stream power and sediment entrainment capacity and considerably magnified flood damage. Resulting morphosedimentary features revealed a close relationship with the valley floor post-flood hydraulic model, and pre-event awareness would have made it possible to predict risk-sensitive areas. In the second case, the Azohía stream, modelling of current pre-flood channel conditions make it possible to determine channel narrowing and entrenchment in the lower alluvial fan stretch. Abandonment of intensive agriculture, basin reforestation and urbanization diminish load contribution and trigger channel incision. This induces an increase in slope and velocity in the bankfull channel, producing renewed erosive energy and thus activating upstream propagation of incision and bank undermining. The absence of water-spreading dynamics on the alluvial fan in favour of confinement in a single channel produces an unstable dynamic in the system, also offering a false sense of stability, as long as no large magnitude floods occur. When modelling flood-prone areas and analysing hydraulic variables, it is important to detect possible anthropic disturbances that may affect basin load budgets in order to anticipate catastrophic consequences resulting from inappropriate fluvial management before the occurrence of an extraordinary event.

The Guadiana River has a significant record of historical floods, but the systematic data record is only 59 years. From layers left by ancient floods we know about we can add new data to the record, and we can estimate maximum discharges of other floods only known by the moment of occurrence and by the damages caused. A hydraulic model has been performed in the area of Pulo de Lobo and calibrated by means of the rating curve of Pulo do Lobo Station. The palaeofloods have been dated by means of ¹⁴C y ¹³⁷Cs. As non-systematic information has been used in order to calculate distribution functions, the quantiles have changed with respect to the same function when using systematic information. The results show a variation in the curves that can be blamed on the human transformations responsible for changing the hydrologic conditions as well as on the latest climate changes. High magnitude floods are related to cold periods, especially at transitional moments of change from cold to warm periods. This tendency has changed from the last medium-high magnitude flood, which took place in a systematic period. Both reasons seem to justify a change in the frequency curves indicating a recent decrease in the return period of big floods over 8000 m³ s⁻¹. The palaeofloods indicate a bigger return period for the same water level discharge thus showing the river basin reference values in its natural condition previous to the transformation of the basin caused by anthropic action.