Gully erosion is one of the most prominent natural denudation processes of the Mediterranean. It causes significant soil degradation and sediment yield. Most traditional field methods for measurement of erosion-induced spatio-temporal changes are time and labor consuming, while their accuracy and precision are highly influenced by various factors. The main research question of this study was how the measurement approach of traditional field sampling methods can be automated and upgraded, while satisfying the required measurement accuracy. The VERTICAL method was developed as a fully automated raster-based method for detection and quantification of vertical spatio-temporal changes within a large number of gully cross-sections (GCs). The developed method was tested on the example of gully Santiš, located at Pag Island, Croatia. Repeat unmanned aerial vehicle (UAV) photogrammetry was used, as a cost-effective and practical method for the creation of very-high-resolution (VHR) digital surface models (DSMs) of the chosen gully site. A repeat aerophotogrammetric system (RAPS) was successfully assembled and integrated into one functional operating system. RAPS was successfully applied for derivation of interval (the two-year research period) DSMs (1.9 cm/pix) of gully Santiš with the accuracy of ±5 cm. VERTICAL generated and measured 2379 GCs, along the 110 m long thalweg of gully Santiš, within which 749 052 height points were sampled in total. VERTICAL proved to be a fast and reliable method for automated detection and calculation of spatio-temporal changes in a large number of GCs, which solved some significant shortcomings of traditional field methods. The versatility and adaptability of VERTICAL allow its application for other, similar scientific purposes, where multitemporal accurate measurement of spatio-temporal changes in GCs is required (e.g., river material dynamics, ice mass dynamics, tufa sedimentation and erosion).