Open Access
Assessment of the main factors affecting the dynamics of nutrients in two rainfed cereal watersheds
(Elsevier, 2020) Hernández García, Iker; Merchán Elena, Daniel; Aranguren Erice, Itxaso; Casalí Sarasíbar, Javier; Giménez Díaz, Rafael; Campo-Bescós, Miguel; Valle de Lersundi, Jokin del; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Ingeniería
Nutrient dynamics and factors that control nutrient exports were observed in two watersheds, namely Latxaga and La Tejería, with similar climatic and management characteristics throughout 10 years (2007–2016). Similar patterns were observed in intra-annual and inter-annual dynamics with higher NO3 − concentration and NO3 −-N yield during the humid seasons (i.e., winters and hydrological year 2013). Regarding concentration, Latxaga showed a higher decrease of nitrate due to a higher development of vegetated areas. High discharge events produced nitrate dilution due to the presence of tile-drainage at La Tejeria. At Latxaga, where tile-drainage was not observed, an increase in concentration occurred as a response to high discharge events. Comparing both watersheds, La Tejería presented ca. 73 ± 25 mg NO3 − L−1 while at Latxaga, the concentration observed was almost three times lower, with ca. 21 ± 15 mg NO3 − L−1 throughout the study period. Similar patterns were observed for the NO3 −-N yield, with 32 kg NO3 −-N ha−1 year−1 and 17 kg NO3 −-N ha−1 year−1 at La Tejería and Latxaga, respectively. Regarding phosphorous, the observed concentrations were 0.20 ± 0.72 mg PO4 3− L−1 and 0.06 ± 0.38 mg PO4 3− L−1 at La Tejería and Latxaga, respectively, with PO4 3−-P yields being 71 kg PO4 3−-P ha−1 year−1 and 33 kg PO4 3−-P ha−1 year−1. Annual phosphate-P yield distribution in both watersheds followed similar patterns to those observed for the nitrate-N yield, with higher yields in the humid season. Regarding concentration, highly erosive rainfall that occurred in summer, mobilizing sediments and probably generating desorption of phosphorous in the stream channel, increased phosphate concentration. This research adds to the knowledge base regarding the dynamics of nutrients and the controlling factors in complex agricultural systems with Mediterranean characteristics.
Open Access
Photogrammetrical and field measurement of gullies with contrasting morphology
(Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa, 2007) Marzolff, I.; Giménez Díaz, Rafael; Seeger, M.; Campo-Bescós, Miguel; Ries, J. B.; Casalí Sarasíbar, Javier; Álvarez-Mozos, Jesús; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Despite a wealth of studies on monitoring different types of gullies by using remote-sensing technique such as photogrammetry, relatively few efforts have been made to test their accuracy. Therefore the question arises as to what extent the accuracy of gully monitoring using photogrammetric technique depends on gully morphology. The objective of this work is to investigate this issue. To do that, we confront field measurements of cross-sectional areas of gullies with contrasting morphology with a similar dataset obtained using photogrammetry. Below, we present the first findings of this investigation.
Open Access
Spatial variability of the relationships of runoff and sediment yield with weather types throughout the Mediterranean basin
(Elsevier, 2019) Peña-Angulo, Dahis; Nadal-Romero, Estela; González-Hidalgo, J. Carlos; Albaladejo, J.; Andreu, V.; Bagarello, Vincenzo; Barhi, H.; Batalla, R.J.; Bernal, S.; Bienes, Ramón; Campo, J.; Campo-Bescós, Miguel; Canatário-Duarte, Antonio; Cantón, Yolanda; Casalí Sarasíbar, Javier; Castillo, Víctor M.; Cerdà, Artemi; Cheggour, A.; Cid, Patricio; Cortesi, N.; Desir, G; Díaz-Pereira, E.; Espigares, T.; Estrany, J.; Fernández-Raga, M.; Ferreira, Carla S.S.; Ferro, Vito; Gallart, F.; Giménez Díaz, Rafael; Gimeno, E.; Gómez, Jose A.; Gómez Gutiérrez, Á.; Gómez-Macpherson, Helena; González-Pelayo, O.; Hueso-González, P.; Kairis, O.; Karatzas, G.P.; Klotz, S.; Kosmas, C.; Lana Renault, Noemí; Lasanta, T.; Latron, J.; Lázaro, Roberto; Le Bissonnais, Y.; Le Bouteiller, Caroline; Licciardello, Feliciana; López-Tarazón, J.A.; Lucía, Ana; Marín, C.; Marqués, M.J.; Martínez-Fernández, J.; Martínez-Mena, M.; Martínez-Murillo, J.F.; Mateos, Luciano; Mathys, N.; Merino-Martin, L.; Moreno de las Heras, Mariano; Moustakas, N.; Nicolau, J.M.; Novara, A.; Pampalone, Vincenzo; Raclot, D.; Rodríguez-Blanco, M. Luz; Rodrigo-Comino, J.; Romero Díaz, A.; Roose, E.; Rubio, J.L.; Ruiz-Sinoga, J.D.; Schnabel, Susanne; Senciales-González, J.M.; Simonneaux, V.; Solé-Benet, A.; Taguas, Encarnación V.; Taboada-Castro, M. Mercedes; Taboada-Castro, M.T.; Todisco, F.; Úbeda, Xavier; Varouchakis, E.A.; Vericat, D.; Wittenberg, L.; Zabaleta, Ane; Zorn, M.; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
Open Access
Evaluation of the impact of changing from rainfed to irrigated agriculture in a mediterranean watershed in Spain
(MDPI, 2023) Oduor, Brian Omondi; Campo-Bescós, Miguel; Lana Renault, Noemí; Alfaro Echarri, Alberto; Casalí Sarasíbar, Javier; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD
The conversion of cultivated areas from rainfed to irrigated agriculture alters the watershed’s hydrology and could affect the water quality and quantity. This study examined how streamflow, nitrate load, and nitrate concentration changed after irrigation implementation in a Mediterranean watershed in Navarre, Spain. The Soil Water Assessment Tool (SWAT) model was applied in the Cidacos River watershed to simulate streamflow and nitrate load under rainfed conditions. The simulated outputs were then compared with the post-irrigation observed values from mid-2017 to 2020 at the watershed outlet in Traibuenas to determine the irrigation impact. The model calibration (2000–2010) and validation (2011–2020) results for streamflow (NSE = 0.82/0.83) and nitrate load (NSE = 0.71/0.68) were satisfactory, indicating the model’s suitability for use in the watershed. A comparison of the rainfed and post-irrigation periods showed an average annual increase in streamflow (952.33 m3 ha−1, +18.8%), nitrate load (68.17 kg ha−1, +62.3%), and nitrate concentration (0.89 mg L−1 ha−1, +79%) at the watershed outlet. Irrigation also caused seasonal changes by altering the cropping cycle and increasing the streamflow and nitrate export during the summer and autumn when irrigation was at its peak. The increases in the post-irrigation period were attributed to the added irrigation water for streamflow and increased nitrogen fertilizer application due to changes in cropping for nitrate concentration and export. These findings are useful to farmers and managers in deciding the best nitrate pollution control and management measures to implement. Furthermore, these results could guide future development and expansion of irrigated lands to improve agricultural sustainability.
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 spatial variability and temporal dynamics on the retrieval of soil moisture from SAR observations
(MDPI, 2009) Álvarez-Mozos, Jesús; Verhoest, Niko E. C.; Larrañaga Urien, Arantzazu; Casalí Sarasíbar, Javier; González de Audícana Amenábar, María; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Radar-based surface soil moisture retrieval has been subject of intense research during the last decades. However, several difficulties hamper the operational estimation of soil moisture based on currently available spaceborne sensors. The main difficulty experienced so far results from the strong influence of other surface characteristics, mainly roughness, on the backscattering coefficient, which hinders the soil moisture inversion. This is especially true for single configuration observations where the solution to the surface backscattering problem is ill-posed. Over agricultural areas cultivated with winter cereal crops, roughness can be assumed to remain constant along the growing cycle allowing the use of simplified approaches that facilitate the estimation of the moisture content of soils. However, the field scale spatial variability and temporal variations of roughness can introduce errors in the estimation of soil moisture that are difficult to evaluate. The objective of this study is to assess the impact of roughness spatial variability and roughness temporal variations on the retrieval of soil moisture from radar observations. A series of laser profilometer measurements were performed over several fields in an experimental watershed from September 2004 to March 2005. The influence of the observed roughness variability and its temporal variations on the retrieval of soil moisture is studied using simulations performed with the Integral Equation Model, considering different sensor configurations. Results show that both field scale roughness spatial variability and its temporal variations are aspects that need to be taken into account, since they can introduce large errors on the retrieved soil moisture values.
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
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
Progress in gully erosion research: IV International Symposium on Gully Erosion
(Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa, 2007) Casalí Sarasíbar, Javier; Giménez Díaz, Rafael; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
The classic forms of water erosion of the soil comprise sheet, rill, and gully erosion. According to the concept most generalized, in sheet erosion, thin layers of material are uniformly removed from the soil surface due to the action of an overland flow, in a homogeneous manner, over the area affected. This results in a normally very gradual and inappreciable loss of soil. However, in rill and gully erosion, the soil loss is caused by the intense action of a concentrated flow, which thus triggers the formation of small or large channels, i.e. rills or gullies, although the mechanisms implicated in either form of concentrated flow erosion are not identical. One of the main differences probably lies in the (much) greater interrelation between the roughness of the channel bed and the hydraulics of the water flow observed in a eroded rill, with respect to what occurs in a typical gully. So, the latter, as a physical process, deserves a special and specific study. Gully erosion, of world-wide importance, is catalogued by some of the principal centres devoted to soil resource conservation as being the foremost problem to be solved. As it is one of the most serious forms of water erosion, this phenomenon is capable of generating major soil losses even though it covers limited land surfaces. Additionally, the damage caused by this type of erosion frequently spreads beyond the area directly affected, i.e. through the siltation of lakes and reservoirs due to the large amounts of sediments it originates. Nevertheless, gully erosion has not received the attention that it warrants from the scientific community. For instance, a rapid search through any important virtual library shows that only less than 10% of soil erosion studies published up to now in international scientific journals deal directly and specifically with gully erosion. More research and surveys are required in order to obtain a better understanding of the physical mechanisms involved in this type of erosion, with the ultimate aim of developing accurate prediction algorithms and efficient control and damage prevention systems. In fact, there are so many unanswered questions on this important environment topic that scientists all over the world have been holding periodic meetings, in which the latest knowledge and advances in the study of gully erosion have been expounded. The first of these meetings was held in Leuven (Belgium) in 2000, the second in Sicuani (China) in 2002 and the third in Oxford, Mississippi (U.S.A.) in 2004. On that last occasion, the participants proposed that Pamplona (Spain) should be the seat of the following meeting, to be held in September 2007. This book contains the abstracts both from the key speeches and from the contributions presented in the IV International Symposium on Gully Erosion held in Pamplona, in September, 2007.