Dissolved solid (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 highflow period. SS concentration was extremely variable, with a range of 2-4 orders of magnitude in concentration for any specific discharge. Temporal variations (both intraand 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

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.

<p>The identification of areas susceptible to gully formation is an objective that has important consequences for erosion control. It allows for the optimization of resources by focusing on prevention and control efforts on the most susceptible areas, avoiding the frequent evolution of ephemeral to permanent gullies. The issue is of great interest in Spanish olive groves, many of which are affected by serious problems of gully erosion.</p><p>Gullies are formed in the swales, which allows the use of topography-based tools to predict their location.</p><p>The Compound Topographic Index (CTI) proposed by Thorne et al. (1986) is calculated for each pixel as an estimate of the flow capacity to cause erosion, as it includes the product of the pixel draining area and its slope. Its application requires the identification of a critical value of the CTI (CTIc), above which the potential areas of gully occurrence will be located. Using historical orthophotos, the gullies observed were digitized for 2011 in the experimental areas called Morente (11 km<sup>2</sup> of traditional olive groves on degraded and poor vertisols) and Matasanos (6 km<sup>2</sup> of intensive olive groves also on vertisols) and nearby area, with cereal crops.</p><p>The objectives of this work are: to identify CTIc values corresponding to cultivated areas in Cordoba, mainly olive groves; to develop and evaluate an application that allows a user without great technical skills to obtain the CTI; to evaluate the capacity of this CTIc to reproduce gullies observed in nearby areas or in different time periods (2005) to establish cause-effect relationships between changes in landuse in this type of phenomenon, using the aforementioned tool.</p><p>Part of the digitized gullies, representative of olive grove areas, were used to obtain the CTIc of each gully, by modifying it until the best reproduction of the gullies observed was achieved, then their average value was taken as CTIc. To calculate the CTI, a 5m resolution DEM was used, obtained from LiDAR PNOA 2014.</p><p>In the framework of the Innolivar project, a desktop GIS application has been developed in a free software environment such as QGIS, which allows the calculation of the CTI. The APET tool (AGNPS Potential Ephemeral Gully Evaluation Tool) recently implemented has helped in the development of this application.</p><p>The CTI calculation by the application, after the determination of the CTIc threshold, serves to identify critical areas from a DEM, which is free and available in many countries. A first qualitative evaluation by visual verification indicates a very good characterization of the gullies. Subsequently, the goodness of fit of the gully position between the digitized gullies and the app-calculated gullies according to the CTIc is evaluated quantitatively by obtaining a binary confusion matrix by lengths. In the Morente area, an error of omission of 29% and of commission of 16% was obtained.</p><p>It can be concluded that the application generated that allows the application of the CTI methodology makes identification of areas susceptible to gully formation possible in an efficient and relatively simple manner, helping to achieve a more sustainable agriculture.</p>

<p>Agricultural activity can have a significant effect on the environment. Often, the lack of experimental data leaves simulation models as the only alternative for understanding and assessing such effects and they can be useful for exploring the response of agricultural systems to different scenarios, in order, for example, to minimize soil erosion or the pollution of watercourses by agrochemicals.</p><p>In this work we present a simulation exercise of the runoff and erosion in two typical olive groves of the Cordoba countryside with contrasting characteristics during the 2009-19 period. The model used is AnnAGNPS, widely tested and very well suited for use in agricultural environments. The specific objectives are: to analyze the applicability of the model confronting its results with data from other nearby areas; to determine the controlling factors of runoff and erosion, such as seasonality; to quantify the importance of the main types of erosion; to explore the response to two different management scenarios.</p><p>The study areas were two, Matasanos (189.4 ha of intensive olive groves on vertisols) and Morente (4.2 ha of traditional olive groves on degraded and poor vertisols). The first scenario (TC) consists of maintaining the soil bare by means of continuous conventional tillage. The second (CC) considers a temporary vegetation cover (around 70 %) on the lanes. All the possible types of erosion in those areas are considered: sheet and rill, ephemeral gullies (EG) and permanent gullies (PG). For the purposes of the simulations, the EGs are tilled while the PGs are not. The latter show more constant characteristics over time (although they also evolve), and are larger in size (i.e., they were assigned a greater depth).</p><p>The results show a significant decrease in average annual runoff in CC with respect to TC (38% in Matasanos and 55% in Morente), which is concentrated in the late autumn and winter months. Thus, according to our simulations, still preliminary, the implementation of covers would have achieved one of its objectives, which is to reduce the runoff generated in the watersheds.</p><p>The sediment yields in both watershed outlets also suffered a significant decrease in CC with respect to TC, going from 4.75 to 1.66 Mg/ha/year and from 16.2 to 6.9 Mg/ha/year in Matasanos and Morente respectively. The simulated erosion rates are consistent with observations made in the area and with other previous simulation exercises. Both sediment export and runoff show a marked seasonality, although erosion occurs somewhat more distributed throughout the year. The different types of erosion take on different importance in each watershed. For example, permanent gullies play a very important role in Morente (46% in TC and 44% in CC), despite they are active at very specific times, probably with extreme events, which is reasonable according to the observations made in the area. The results show that the model is apparently useful with respect to the proposed objectives, allowing the effect of different uses and management on the environment to be contrasted in the medium and long term.</p>

<p>Ephemeral gullies (EG) are linear erosion features located in swales where surface and/or subsurface runoff concentrate during or immediately after rainfall events. As its name states, EGs are temporary because they are easily filled by conventional machinery, but when filled they reform if the area is not appropriately managed. Downstream water quality issues and decreased soil productivity are the main environmental impacts. EGs are frequently identified as (the most) relevant sediment sources in agricultural areas but their dimensions and particular contribution to the total erosion under different temporal, spatial, climate and land use condition is still unknown. Therefore, the objective of this study is to obtain ephemeral gully erosion rates and estimate the main morphological characteristics of the ephemeral gullies (width, length and depth) and their evolution both in relation to time and position on the landscape.</p><p>The studied EGs, B6 with a 0.94 ha watershed and I3 with a 0.95 ha watershed formed in two fields located in the Walnut Creek watershed, Iowa (US). The field-sized watersheds are less than 1.5 Km apart and have similar topography and soils. The cropping system consists of a two-year corn-soybean rotation managed by one farmer using no-till and other standard management practices. EG were measured using close range photogrammetry techniques. In order to achieve a suitable characterization of the EG evolution over time and space, EGs were divided in three sections (bottom, middle and top) of equal length. Photographs were taken at least once in 2013, 2014 and 2018 (a total of five in I3 and three in B6). Cross section profiles along the EG perpendicular to the flow path direction were selected and their width, area and depth were determined from a graphical representation of the cross sections. EG volumes were estimated by the sum of interpolating sequential cross-section areas and multiplying by the distance between them.</p><p>Average EG erosion rates during 2013-2014 were 3.19 Mg ha<sup>-1</sup> year<sup>-1</sup> for B6 and 3.63 Mg ha<sup>-1</sup> year<sup>-1</sup> for I3. Values in agreement with rates estimated by the United States Department of Agriculture (USDA) Natural Resources Conservation Service (NRCS) of 0.49 to 5.18 Mg ha<sup>-1</sup> year<sup>-1 </sup>across the USA and other simulated values of 4.00 ± 1.76 Mg ha<sup>-1</sup> year<sup>-1</sup> for no till systems in the state of Iowa. The current study shows evidences that EG in no till systems may not stabilize after their formation. EG dimensions (depth, width and length, thus volume) varied over time and space during the continuously monitored period. In general, volumes tend to increase in the middle position while depths decrease in the bottom position. When the EG was filled, it reformed again in approximately the same location showing similar dimensions to that which existed prior to filling.</p>