Open Access
Runoff, nutrients, sediment and salt yields in an irrigated watershed in southern Navarre (Spain)
(Elsevier, 2018) Merchán Elena, Daniel; Casalí Sarasíbar, Javier; Valle de Lersundi, Jokin del; Campo-Bescós, Miguel; Giménez Díaz, Rafael; Preciado, Beatriz; Lafarga, Alberto; Landa Ingeniaritza eta Proiektuak; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Proyectos e Ingeniería Rural
The environmental impact of irrigated agriculture on water quality was assessed in Landazuria watershed (Navarre, northeast Spain), a 479.5 ha watershed with 53% of irrigated agricultural land. In the framework of a long-term monitoring program, precipitation and discharge were measured at 10-min intervals and compound daily water samples were collected during the agricultural years (September to August) 2007–2016, and analysed for nitrate (NO3−), phosphate (PO43−), sediment and total dissolved solids (TDS) concentrations. Typical agricultural management (including crop surfaces, irrigation and fertilization rates) was obtained from inquiries to farmers. Concentration and yield of the studied variables presented a high degree of variation, both intra- and inter-annual. Median concentration for the entire study period were 185, <0.05, 31 and 2284 mg L−1 for NO3−, PO43−, sediment and TDS, respectively. NO3−-N and PO43−-P yields averaged 74 and 0.04 kg ha−1 year−1, respectively. NO3 −-N yield was higher than in other agricultural land uses in Navarre and in the order of magnitude of other irrigated areas in the Middle Ebro Valley. PO43−-P yield was in the same order of magnitude than in rainfed watersheds in Navarre but lower than in intensively grazed watersheds. Sediment yield was extremely variable, averaging 360 kg ha−1 year−1, with 44% of the total measured load recorded in a few days. It was in the lower range of those measured in Navarre for rainfed agriculture and similar to those estimated in other irrigated areas of the Middle Ebro River. TDS concentration presented a significant decreasing trend since available salts were being washed out, while TDS yield averaged 1.8 Mg ha−1 year−1. Long-term monitoring of irrigated areas is required to understand pollution processes in these agroecosystems and to adequately characterize the environmental impact of current agricultural practices on water quality, in order to implement, and adequately assess, measures to reduce agricultural pollution.
Open Access
Challenges and progresses in the detailed estimation of sediment export in agricultural watersheds in Navarra (Spain) after two decades of experience
(Elsevier, 2023) Barberena Ruiz, Íñigo; Luquin Oroz, Eduardo Adrián; Campo-Bescós, Miguel; Eslava, Javier; Giménez Díaz, Rafael; Casalí Sarasíbar, Javier; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
Soil erosion is a very serious environmental problem worldwide, with agriculture considered the main source of sediment in inland waters. In order to determine the extent and importance of soil erosion in the Spanish region of Navarra, in 1995 the Government of Navarra established the Network of Experimental Agricultural Watersheds (NEAWGN), which consists of five small watersheds representative of local conditions. In each watershed, key hydrometeorological variables, including turbidity, were recorded every 10 min, and daily samples were taken to determine suspended sediment concentration. In 2006, the frequency of suspended sediment sampling was increased during hydrologically relevant events. The main objective of this study is to explore the possibility of obtaining long and accurate time series of suspended sediment concentration in the NEAWGN. To this end, simple linear regressions between sediment concentration and turbidity are proposed. In addition, supervised learning models incorporating a larger number of predictive variables are used for the same purpose. A series of indicators are proposed to objectively characterize the intensity and timing of sampling. It was not possible to obtain a satisfactory model for estimating the concentration of suspended sediment. This would be mainly due to the large temporal variability found of the physical and mineralogical characteristics of the sediment, which would be affecting the turbidity value, independently of the sediment concentration, per se. This fact would be particularly important in small river watersheds such as those of this study, and especially if their physical conditions are spatially and temporally radically disturbed by agricultural tillage and by a constant modification of the vegetation cover, as is the case in cereal basins. Our findings suggest that better results could be obtained by including in the analysis variables such as soil texture and exported sediment texture, rainfall erosivity, and the state of vegetation cover and riparian vegetation.
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
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
Hydrological signatures based on event runoff coefficients in rural catchments of the Iberian Peninsula
(Wolters Kluwer, 2017) Taguas, Encarnación V.; Nadal-Romero, Estela; Ayuso, José L.; Casalí Sarasíbar, Javier; Cid, Patricio; Dafonte, Jorge; Canatário-Duarte, Antonio; Ferreira, Carla S.S.; Giménez Díaz, Rafael; Giráldez Cervera, Juan Vicente; Gómez-Macpherson, Helena; Gómez, Jose A.; González-Hidalgo, J. Carlos; Lana Renault, Noemí; Lucía, Ana; Mateos, Luciano; Pérez, Rafael; Rodríguez-Blanco, M. Luz; Schnabel, Susanne; Serrano-Muela, M. Pilar; Taboada-Castro, M. Mercedes; Zabaleta, Ane; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Hydrological signatures are indices that help to describe the behavior of catchments. These indices can also be used to transfer information from gauged to ungauged catchments. In this study, different approaches were evaluated to determine volumetric runoff coefficients in 18 small/ medium experimental gauged catchments of the Iberian Peninsula and to fit runoff calculations based on precipitation data for gauged and ungauged catchments. Using data derived from 1962 events, rainfall-runoff relationships were characterized and compared in order to evaluate the various hydrological response patterns. Volumetric runoff coefficients and cumulative runoff and precipitation ratios of the events that generated runoff (Rcum) minimized the root mean square error. A linear fit for the estimation of Rcum in ungauged atchments was based on mean annual precipitation, rates of infiltration, the fraction of forest-land use, and the catchment channel length. Despite high catchment heterogeneity, Rcum resulted in a suitable parameter to evaluate hydrological variability in rural gauged and ungauged catchments. In 50% of the catchments, the precipitation accounted for less than 50% of the runoff variation. Annual precipitation, antecedent rainfall, and base flow did not have a high significance in rainfall-runoff relationships, which illustrates the heterogeneity of hydrological responses. Our results highlight the need for signature characterizations of small/medium rural catchments because they are the sources of runoff and sediment discharge into rivers, and it is more economical and efficient to take action to mitigate runoff in rural locations.
Open Access
In-situ assessment of the spatial arrangement of step-pool units on eroded rills
(wiley, 2019) Govers, G.; Campo-Bescós, Miguel; Zubieta Laseca, Elena; Giménez Díaz, Rafael; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Ingeniería
Open Access
Assessment of soil factors controlling ephemeral gully erosion on agricultural fields
(Wiley, 2017) Ollobarren del Barrio, Paul; Campo-Bescós, Miguel; Giménez Díaz, Rafael; Casalí Sarasíbar, Javier; Landa Ingeniaritza eta Proiektuak; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Proyectos e Ingeniería Rural; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion.
Open Access
Influence of surface roughness sample size for C-band SAR backscatter applications on agricultural soils
(IEEE, 2017) Martínez de Aguirre Escobar, Alejandro; Álvarez-Mozos, Jesús; Lievens, Hans; Verhoest, Niko E. C.; Giménez Díaz, Rafael; Landa Ingeniaritza eta Proiektuak; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Proyectos e Ingeniería Rural
Soil surface roughness determines the backscatter coefficient observed by radar sensors. The objective of this letter was to determine the surface roughness sample size required in synthetic aperture radar applications and to provide some guidelines on roughness characterization in agricultural soils for these applications. With this aim, a data set consisting of ten ENVISAT/ASAR observations acquired coinciding with soil moisture and surface roughness surveys has been processed. The analysis consisted of: 1) assessing the accuracies of roughness parameters s and l depending on the number of 1-m-long profiles measured per field; 2) computing the correlation of field average roughness parameters with backscatter observations; and 3) evaluating the goodness of fit of three widely used backscatter models, i.e., integral equation model (IEM), geometrical optics model (GOM), and Oh model. The results obtained illustrate a different behavior of the two roughness parameters. A minimum of 10-15 profiles can be considered sufficient for an accurate determination of s, while 20 profiles might still be not enough for accurately estimating l. The correlation analysis revealed a clear sensitivity of backscatter to surface roughness. For sample sizes >15 profiles, R values were as high as 0.6 for s and ~0.35 for l, while for smaller sample sizes R values dropped significantly. Similar results were obtained when applying the backscatter models, with enhanced model precision for larger sample sizes. However, IEM and GOM results were poorer than those obtained with the Oh model and more affected by lower sample sizes, probably due to larger uncertainly of l.
Open Access
Assessing hillslope-channel connectivity in an agricultural catchment using rare-earth oxide tracers and random forests models
(Universidad de la Rioja, 2017) Masselink, Rens Hein; Temme, A.J.A.M.; Giménez Díaz, Rafael; Casalí Sarasíbar, Javier; Keesstra, Saskia D.; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
Soil erosion from agricultural areas is a large problem because of off-site effects like the rapid filling of reservoirs. To mitigate the problem of sediments from agricultural areas reaching the channel, reservoirs and other surface areas it is important to understand hillslope-channel connectivity and catchment connectivity. To determine the functioning of hillslope-channel connectivity and the continuation of transport of these sediments in the channel, it is necessary to obtain data on sediment transport from the hillslopes to the channels. Simultaneously, the factors that influence sediment export out of the catchment need to be studied. For measuring hillslope-channel sediment connectivity, rare-earth oxide (REO) tracers were applied to a hillslope in an agricultural catchment in Navarre, Spain, preceding the winter of 2014-2015. The results showed that during the winter no sediment transport from the hillslope to the channel was detected. To test the implication of the REO results at the catchment scale, two contrasting conceptual models for sediment connectivity were assessed using a random forest (RF) machine learning method. The RF method was applied using a 15-year period of measured sediment output at the catchment scale. One model proposes that small events provide sediment for large events, while the other proposes that only large events cause sediment detachment and small events subsequently remove these sediments from near and in the channel. For sediment yield prediction of small events, variables related to large preceding events were the most important. The model for large events underperformed and, therefore, we could not draw any immediate conclusions whether small events influence the amount of sediment exported during large events. Both REO tracers and RF method showed that low intensity events do not contribute any sediments from the hillslopes to the channel in Latxaga catchment. Sediment dynamics are dominated by sediment mobilization during large (high intensity) events. Sediments are for a large part exported during those events, but the system shows a memory of the occurrence of these large events, suggesting that large amounts of sediments are deposited in and near the channel after these events. These sediments are gradually removed by small events. To better understand the delivery if sediments to the channel and how large and small events influence each other, more field data on hillslope-channel connectivity and within-channel sediment dynamics is necessary.
Open Access
On the influence of spatial resolution in soil surface roughness characterization using Tls and Sfm techniques
(IEEE, 2018) Martínez de Aguirre Escobar, Alejandro; Álvarez-Mozos, Jesús; Giménez Díaz, Rafael; Milenković, Milutin; Pfeifer, Norbert; Ingeniería; Ingeniaritza
Soil surface roughness strongly affects the scattering of microwaves and determines the backscattering coefficient observed by SAR (Synthetic Aperture Radar) sensors. The aim of this study is to analyze the influence of the spatial resolution of Terrestrial Laser Scanner (TLS) and Structure from Motion (SfM) techniques to parameterize surface roughness over agricultural soils. Three experimental plots (5 x 5 meters) representing different roughness conditions were measured by TLS and SfM techniques. Roughness parameters (s and l) were calculated from profiles obtained at different spatial resolutions in parallel and in perpendicular to the tillage direction on each plot. The results showed minor differences in the parameters values between both techniques and, in general, a decreasing trend and an increasing trend for lower spatial resolutions for parameter s and l, respectively.