Open Access
EUSEDcollab: a network of data from European catchments to monitor net soil erosion by water
(Nature Research, 2023-08-04) Matthews, Francis; Verstraeten, Gert; Borrelli, Pasquale; Vanmaercke, Matthias; Poesen, J.; Steegen, An; Degré, Aurore; Cárceles Rodríguez, Belén; Bielders, Charles; Franke, Christine; Alary, Claire; Zumr, David; Patault, Edouard; Nadal-Romero, Estela; Smolska, Ewa; Licciardello, Feliciana; Swerts, Gilles; Thodsen, Hans; Casalí Sarasíbar, Javier; Eslava, Javier; Richet, Jean-Baptiste; Ouvry, Jean-François; Farguell, Joaquim; Święchowicz, Jolanta; Nunes, João Pedro; Pak, Lai Ting; Liakos, Leonidas ; Campo-Bescós, Miguel; Żelazny, Mirosław; Delaporte, Morgan; Pineux, Nathalie; Henin, Nathan; Bezak, Nejc; Lana Renault, Noemí ; Tzoraki, Ourania; Giménez Díaz, Rafael; Li, Tailin; Durán Zuazo, Víctor Hugo; Bagarello, Vincenzo; Pampalone, Vincenzo; Ferro, Vito; Úbeda, Xavier; Panagos, Panos; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
As a network of researchers we release an open-access database (EUSEDcollab) of water discharge and suspended sediment yield time series records collected in small to medium sized catchments in Europe. EUSEDcollab is compiled to overcome the scarcity of openaccess data at relevant spatial scales for studies on runoff, soil loss by water erosion and sediment delivery. Multi-source measurement data from numerous researchers and institutions were harmonised into a common time series and metadata structure. Data reuse is facilitated through accompanying metadata descriptors providing background technical information for each monitoring station setup. Across ten European countries, EUSEDcollab covers over 1600 catchment years of data from 245 catchments at event (11 catchments), daily (22 catchments) and monthly (212 catchments) temporal resolution, and is unique in its focus on small to medium catchment drainage areas (median = 43 km2, min = 0.04 km2, max = 817 km2) with applicability for soil erosion research. We release this database with the aim of uniting people, knowledge and data through the European Union Soil Observatory (EUSO).
Open Access
Evaluation of 2D models for the prediction of surface depression storage using realistic reference values
(Wiley, 2016) Giménez Díaz, Rafael; Mezkiritz Barberena, Irantzu; Campo-Bescós, Miguel; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Martínez de Aguirre Escobar, Alejandro; Casalí Sarasíbar, Javier; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Depression storage (DS) is the maximum storage of precipitation and runoff in the soil surface at a given slope. The DS is determined by soil roughness that in agricultural soils is largely affected by tillage. The direct measurement of DS is not straightforward because of the natural permeability of the soil. Therefore, DS has generally been estimated from 2D/3D empirical relationships and numerical algorithms based on roughness indexes and height measurements of the soil surface, respectively. The objective of this work was to evaluate the performance of some 2D models for DS, using direct and reliable measurements of DS in an agricultural soil as reference values. The study was carried out in experimental microplots where DS was measured in six situations resulting from the combination of three types of tillage carried out parallel and perpendicular to the main slope. Those data were used as reference to evaluate four empirical models and a numerical method. Longitudinal altitudinal profiles of the relief were obtained by a laser profilometer. Infiltration measurements were carried out before and after tillage. The DS was largely affected by tillage and its direction. Highest values of DS are found on rougher surfaces mainly when macroforms cut off the dominant slope. The empirical models had a limited performance while the numerical method was the most effective, even so, with an important variability. In addition, a correct hydrological management should take into account that each type of soil tillage affects infiltration rate differently.
Open Access
Evaluation of nitrate soil probes for a more sustainable agriculture
(MDPI, 2022) Bellosta Diest, Amelia; Campo-Bescós, Miguel; Zapatería Miranda, Jesús; Casalí Sarasíbar, Javier; Arregui Odériz, Luis Miguel; Agronomía, Biotecnología y Alimentación; Agronomia, Bioteknologia eta Elikadura; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Gobierno de Navarra / Nafarroako Gobernua
Synthetic nitrogen (N) fertilizers and their increased production and utilization have played a great role in increasing crop yield and in meeting the food demands resulting from population growth. Nitrate (NO3−) is the common form of nitrogen absorbed by plants. It has high water solubility and low retention by soil particles, making it prone to leaching and mobilization by surface water, which can seriously contaminate biological environments and affect human health. Few methods exist to measure nitrate in the soil. The development of ion selective sensors provides knowledge about the dynamics of nitrate in the soil in real time, which can be very useful for nitrate management. The objective of this study is to analyze the performance of three commercial probes (Nutrisens, RIKA and JXCT) under the same conditions. The performance was analyzed with respect to electrical conductivity (EC) (0–50 mS/cm) and nitrate concentration in aqueous solution and in sand (0–180 ppm NO3−) at 35% volumetric soil moisture. Differences were shown among probes when studying their response to variations of the EC and, notably, only the Nutrisens probe provided coherent accurate measurements. In the evaluation of nitrate concentration in liquid solution, all probes proved to be highly sensitive. Finally, in the evaluation of all probes’ response to modifications in nitrate concentration in sand, the sensitivity decreased for all probes, with the Nutrisens probe the most sensitive and the other two probes almost insensitive.
Embargo
Introducing QAnnAGNPS - A QGIS plugin to facilitate the use of AnnAGNPS (Annualized Agricultural Nonpoint source model)
(Elsevier, 2024) Barberena Ruiz, Íñigo; Campo-Bescós, Miguel; Casalí Sarasíbar, Javier; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
AnnAGNPS is a watershed-scale hydrological model designed to analyze the impact of nonpoint source pollution in agricultural environments. Its unique capabilities have defined it as an essential model that is used globally to evaluate agricultural management scenarios. However, it does not currently have a user-friendly graphical interface that provides a simple way for users to perform simulations. This article presents QAnnAGNPS, a plugin developed in QGIS to facilitate access to the simulation capabilities of AnnAGNPS through a user-friendly interface and the addition of extra features, including data visualization. QAnnAGNPS, in addition to fulfilling this valuable task, opens the door to the incorporation of additional functions already included in other similar hydrological models. The plugin has been tested in the Latxaga cereal basin in Navarra, Spain, and has demonstrated that it provides a simpler way to perform simulations and visualize results compared to AnnAGNPS.
Open Access
Relationship of weather types on the seasonal and spatial variability of rainfall, runoff, and sediment yield in the western Mediterranean basin
(MDPI, 2020) Peña-Angulo, Dahis; Nadal-Romero, Estela; Campo-Bescós, Miguel; Casalí Sarasíbar, Javier; Giménez Díaz, Rafael; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Ingeniería
Rainfall is the key factor to understand soil erosion processes, mechanisms, and rates. Most research was conducted to determine rainfall characteristics and their relationship with soil erosion (erosivity) but there is little information about how atmospheric patterns control soil losses, and this is important to enable sustainable environmental planning and risk prevention. We investigated the temporal and spatial variability of the relationships of rainfall, runoff, and sediment yield with atmospheric patterns (weather types, WTs) in the western Mediterranean basin. For this purpose, we analyzed a large database of rainfall events collected between 1985 and 2015 in 46 experimental plots and catchments with the aim to: (i) evaluate seasonal differences in the contribution of rainfall, runoff, and sediment yield produced by the WTs; and (ii) to analyze the seasonal efficiency of the different WTs (relation frequency and magnitude) related to rainfall, runoff, and sediment yield. The results indicate two different temporal patterns: the first weather type exhibits (during the cold period: autumn and winter) westerly flows that produce the highest rainfall, runoff, and sediment yield values throughout the territory; the second weather type exhibits easterly flows that predominate during the warm period (spring and summer) and it is located on the Mediterranean coast of the Iberian Peninsula. However, the cyclonic situations present high frequency throughout the whole year with a large influence extended around the western Mediterranean basin. Contrary, the anticyclonic situations, despite of its high frequency, do not contribute significantly to the total rainfall, runoff, and sediment (showing the lowest efficiency) because of atmospheric stability that currently characterize this atmospheric pattern. Our approach helps to better understand the relationship of WTs on the seasonal and spatial variability of rainfall, runoff and sediment yield with a regional scale based on the large dataset and number of soil erosion experimental stations.
Open Access
Extended assessment of sprinkler irrigation uniformity in greenhouses using GIS and hydraulic modeling
(MDPI, 2022) Barberena Ruiz, Íñigo; Campo-Bescós, Miguel; Casalí Sarasíbar, Javier; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Ingeniería; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Traditionally, distribution uniformity has been obtained by using rain gauges, which makes it a very expensive process. This paper sought to create a simulation strategy using QGIS and EPANET, both free software, that allowed the simulation of the water application results of all the emitters of an irrigation installation. In this way, it was possible to obtain the geospatial representation of the applied water and finally to know the distribution uniformity in the whole installation. The simulation finally fulfilled its objective and was compared with a study of distribution uniformity with rain gauges. The biggest difference between the measured and simulated data was a difference of 5.76% among the sectors. The simulated uniformity was very similar to the measured uniformity, which allowed us to affirm that the proposed simulation methodology was adequate. We believe that the methodology proposed in this article could be very useful in improving the management of sprinkler irrigation systems, particularly those in which distribution uniformity is of special importance. These improvements in management can also result in savings in water and other inputs, which are becoming increasingly important in the current context of climate change and the reduction in the impact of agriculture on the environment. Finally, similar studies could be carried out with the same tools for other pressurized irrigation systems, such as sprinkler irrigation outside greenhouses and drip irrigation.
Open Access
Accuracy of methods for field assessment of rill and ephemeral gully erosion
(Elsevier, 2006) Casalí Sarasíbar, Javier; Loizu Maeztu, Javier; Campo-Bescós, Miguel; Santisteban Comino, Luisa María de; Álvarez-Mozos, Jesús; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
To properly assess soil erosion in agricultural areas, it is necessary to determine precisely the volume of ephemeral gullies and rills in the field by using direct measurement procedures. However, little information is available on the accuracy of the different methods used. The main purpose of this paper is to provide information for a suitable assessment of rill and ephemeral gully erosion with such direct measurement methods. To achieve this objective: a) the measurement errors associated to three methods used for field assessment of channel cross sectional areas are explored; b) the influence of the number of cross sections used per unit channel length on the assessment accuracy, is analysed and; c) the effect of the channel size and shape on measurement errors is examined. The three methods considered to determine the cross sectional areas were: micro-topographic profile meter (1); detailed measurement of section characteristic lengths with a tape (2); measurement of cross section width and depth with a tape (3). Five reaches of different ephemeral gully types 14.0 or 30.0 m long and a set of six 20.4 to 29.4 m long rill reaches were selected. On each gully reach, the cross sectional areas were measured using the three above mentioned methods, with a separation (s) between cross sections of 1 m. For rills, the cross sectional areas were measured with methods 1 and 3, with s= 2 m. Then, the corresponding total erosion volumes were computed. The volume calculated with method 1 with s= 1 m for gullies and s= 2 m for rills was taken as the reference method. For each channel, and for each one of the possible combinations of s and measurement method (m), the relative measurement error and the absolute value of the relative measurement error (Ersm and |Ersm| ), defined with respect to the reference one, was calculated. |Ersm| much higher than 10% were obtained very easily, even for small s values and for apparently quasi prismatic channels. Channel size and shape had a great influence on measurement errors. In fact, the selection of the more suitable method for a certain gully shape and size seemed to be much more important than s, at least when s< 10 m. Method 1 always provided the most precise measurements, and its results were the less dependent on s. However, s must be <5 m to guarantee an error smaller than 10%. Method 2 is not recommended, because it is difficult, time consuming and can lead to large errors. Method 3 seems to be enough for small, wide and shallow gullies, and for small rills, but only if s is shorter than 5 m. Results obtained after the analysis of rill measurement errors were similar to those of gullies. The analysis of Ersm and |Ersm| when calculating channel volumes using a unique representative cross section highlighted the importance of correctly selecting the adequate cross section. Due to the high error values that this method can entail, it is not considered as advisable whenever accurate erosion measurements are pursued.
Open Access
Dissolved solids and suspended sediment dynamics from five small agricultural watersheds in Navarre, Spain: a 10-year study
(Elsevier, 2019) Merchán Elena, Daniel; Luquin Oroz, Eduardo Adrián; Hernández García, Iker; Campo-Bescós, Miguel; Giménez Díaz, Rafael; Casalí Sarasíbar, Javier; Valle de Lersundi, Jokin del; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Ingeniería; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Dissolved solids (DS) and suspended sediment (SS) loads are considered relevant environmental problems. They are related to a wide range of on-site and off-site impacts, such as soil erosion or salinization of water bodies. In this study, the dynamics of DS and SS concentrations and loads were assessed in five small watersheds covering representative agricultural land uses in Navarre (Spain). To this end, discharge, DS and SS concentration data were collected during ten hydrological years at each watershed outlet, and loads were computed from discharge and concentration values. DS concentration followed a seasonal pattern imposed by the availability of water, with higher concentrations recorded in low-flow periods and lower concentration in the high-flow period. SS concentration was extremely variable, with a range of 2–4 orders of magnitude in concentration for any specific discharge. Temporal variations (both intra- and inter-annual) in DS loads were explained by differences in runoff, whereas those of SS were not, being the SS loads associated mainly with specific high flow events. These temporal patterns were observed for both agricultural (this study) and non-agricultural (literature) watersheds. From the data in the Navarrese watersheds and those available in the literature, we inferred that agricultural land use, in general, tends to increase the concentration of both DS and SS. Regarding DS and SS yields, the effects of agricultural land use on DS yields are controlled by the changes in runoff rather than the (small) changes in DS concentration. In this sense, land uses changes expected to increase runoff (i.e., a shift from forested to arable or from rainfed to irrigated agriculture) would increase DS yields. On the other hand, agricultural land use tends to increase SS yields, although the effect is highly variable depending on site-specific factors, both natural (e.g., watershed shape) and anthropogenic (e.g., degree of soil conservation practices). In the Navarrese watersheds, DS yields ranged from 1.1 to 2.2 Mg ha−1 year−1 whereas SS yields ranged from 0.3 to 4.3 Mg ha−1 year−1. DS yields seem to dominate under non-agricultural conditions and in most agricultural land uses at the small watershed scale. On the other hand, SS yields dominate in watersheds with increased soil erosion as a consequence of arable land use over erosion-prone watersheds.
Open Access
Hydrological records can be used to reconstruct the resilience of watersheds to climatic extremes
(Nature Research, 2024) Huffaker, Ray; Campo-Bescós, Miguel; Luquin Oroz, Eduardo Adrián; Casalí Sarasíbar, Javier; Muñoz Carpena, Rafael; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
Hydrologic resilience modeling is used in public watershed management to assess watershed ability to supply life-supporting ecoservices under extreme climatic and environmental conditions. Literature surveys criticize resilience models for failing to capture watershed dynamics and undergo adequate testing. Both shortcomings compromise their ability to provide management options reliably protecting water security under real-world conditions. We formulate an empirical protocol to establish real-world correspondence. The protocol applies empirical nonlinear dynamics to reconstruct hydrologic dynamics from watershed records, and analyze the response of reconstructed dynamics to extreme regional climatic conditions. We devise an AI-based early-warning system to forecast (out-of-sample) reconstructed hydrologic resilience dynamics. Application to the La Tejería (Spain) experimental watershed finds it to be a low dimensional nonlinear deterministic dynamic system responding to internal stressors by irregularly oscillating along a watershed attractor. Reconstructed and forecasted hydrologic resilience behavior faithfully captures monthly wet-cold/dry-warm weather patterns characterizing the Mediterranean region.
Open Access
Effect of topography on retreat rate of different gully headcuts in Bardenas Reales area (Navarre, Spain)
(Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa, 2007) Campo-Bescós, Miguel; Álvarez-Mozos, Jesús; Casalí Sarasíbar, Javier; Giménez Díaz, Rafael; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
In Northeast Spain, gullying is a widespread phenomenon. This type of erosion is especially intense in Bardenas Reales (Navarre) where at least two major typical kinds of gully headcut are present. A first group developed in soil material (named, conventional gully headcut), and second group of gully headcut with a sandstone layer as a top horizon (named, sandstone gully headcut). In addition, within the former group, we can distinguish a subgroup of gully headcuts developed in soils particularly prone to piping and tunnelling due to the dispersive condition of the materials (named piping associated gully headcut). In this situation, a question arises: to what extent simple topographic parameters account for the retreat rate of the different kind of gully headcuts observed in the region of Bardenas Reales? The aim of this study was to investigate and gain insight in this issue.