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
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
Implications of recent experimental findings for rill erosion modeling
(Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa, 2007) Govers, G.; Giménez Díaz, Rafael; Van Oost, K.; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
The major aim of this paper is to critically review the theoretical concepts that are underpinning current models of rill flow and sediment detachment in the light of recent experimental results and, when necessary, to propose modifications to the theoretical formulations so that they are in agreement with experimental evidence. We also investigate to what extent a detachment model of reduced complexity, which is based on experimental observations, is consistent with field observation on the effect of topography on rill erosion.
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
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
Gully geometry: what are we measuring?
(European Geosciences Union, 2015) Casalí Sarasíbar, Javier; Giménez Díaz, Rafael; Campo-Bescós, Miguel; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Much of the research on (ephemeral) gully erosion comprises the determination of the geometry of these eroded channels, especially their width and depth. This is not a simple task due to uncertainty generated by the wide range of variability in gully cross section shapes found in the field. However, in the literature, this uncertainty is not recognized so that no criteria for their measurement are indicated. The aim of this work is to make researchers aware of the ambiguity that arises when characterizing the geometry of an ephemeral gully and similar eroded channels. In addition, a measurement protocol is proposed with the ultimate goal of pooling criteria in future works. It is suggested that the geometry of a gully could be characterized through its mean equivalent width and mean equivalent depth, which, together with its length, define an “equivalent prismatic gully” (EPG). The latter would facilitate the comparison between different gullies.
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
Evaluación de la producción de sedimentos y calidad de las aguas en cuencas agrarias de Navarra
(Universidad de La Rioja, 2012) Giménez Díaz, Rafael; Casalí Sarasíbar, Javier; Díez Beregaña, Javier; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Gobierno de Navarra / Nafarroako Gobernua
The Government of Navarre (Spain) established a series of experimental watersheds on different sites of its territory in order to assess the impact of the agricultural activities on the environment. Hydrological/meteorological and water quality data are continuously monitored. In this paper, data recorded –precipitation, runoff, sediment, nitrate and phosphate– in two grain-sown watersheds (La Tejería and Latxaga) and in a third one (Oskotz) covered with forest and pasture (cattle-breeding) with a sub-watershed almost entirely under forest, are analyzed; the information covers a period of ca. 15 years. The inter/intra annual variability of the involved processes are studied as well as their impact on the water quality and on the ecosystem. The majority of the sediments yielded in the watersheds are registered during winter, though most of the erosive rainfalls occurred in summer. This is because of the easy-prone runoff that prevailed during winter due to the (much) higher water content of the soils in this period. The grain-sown watersheds –roughly similar at first view– present however important differences regarding sediment yield: 1800 kg/ha year at La Tejería and only 450 kg/ha year at Latxaga. This could be mainly due to differences in morphology, topography, and amount of stream channel vegetation between both sites. At Oskotz this figures are as a whole lower than those registered in the cultivated watersheds but eventually larger than those at Latxaga; this depending on the forest management (logging). Water quality is seriously affected at La Tejería with nitrate concentration constantly over critical threshold. By contrast, nitrate values are much lower in Latxaga watershed. However, phosphate concentrations were similar in both grain-sown watersheds corresponding to water with a significant risk of eutrophication. The cattle-breeding activity at Oskotz generated unexpected levels of phosphate much higher than those registered in the cultivated wathersheds.
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
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.