Abstract. Gully erosion is an important erosive process in Mediterranean basins. However, the long-term dynamics of gully networks and the variations in sediment production in gullies are not well known. Available studies are often conducted only over a few years, while many gully networks form, grow, and change in response to environmental and land use or management changes over a long period. In order to clarify the effect of these changes, it is important to analyse the evolution of the gully network with a high temporal resolution. This study aims at analysing gully morphodynamics over a long timescale (1956–2013) in a large Mediterranean area in order to quantify gully erosion processes and their contribution to overall sediment dynamics.

A gully network of 20 km² located in southwestern Spain has been analysed using a sequence of 10 aerial photographs in the period 1956–2013. The extension of the gully network both increased and decreased in the study period. Gully drainage density varied between 1.93 km km⁻² in 1956, a minimum of 1.37 km km⁻² in 1980, and a maximum of 5.40 km km⁻² in 2013. The main controlling factor of gully activity appeared to be rainfall. Land use changes were found to have only a secondary effect. A new Monte Carlo-based approach was proposed to reconstruct gully erosion rates from orthophotos. Gully erosion rates were found to be relatively stable between 1956 and 2009, with a mean value of 11.2 t ha⁻¹ yr⁻¹. In the period 2009–2011, characterized by severe winter rainfalls, this value increased significantly to 591 t ha⁻¹ yr⁻¹. These results show that gully erosion rates are highly variable and that a simple interpolation between the starting and ending dates greatly underestimates gully contribution during certain years, such as, for example, between 2009 and 2011. This illustrates the importance of the methodology applied using a high temporal resolution of orthophotos.