Differential Interferometry over Sentinel-1 TopSAR Images as a Tool for Water and Tillage Soil Erosion Analysis

: Erosion is a major problem on agricultural lands in Europe. Erosion measurement tools have traditionally been focused on delocalized quantiﬁcation but without mapping the real places inside the basin where the erosion took effect. In this article, we use new space missions, such as Sentinel-1, and the opportunity they offer to obtain SAR (Synthetic Aperture Radar) images with high frequency, resolution, range, and, above all, availability to enable the application of techniques, like differential interferometry, in new ﬁelds. We propose to measure ground deformation caused by water and tillage erosion in small agricultural basins using TopSAR (Terrain Observation with Progressive Scans SAR, Synthetic Aperture Radar) images acquired by the Sentinel-1 mission, after previously verifying the accuracy of these measurements through comparison with measurements from a LIDAR (Light Detection and Ranging) system. The results of this work conﬁrm the potential of monitoring erosion in agricultural basins with differential interferometry over Sentinel-1 TopSAR images. Its capabilities have been successfully tested in different conditions related to agricultural tasks without precipitation or storm events. This technique makes it possible to study both water and tillage erosion and sedimentation processes and even to test the efﬁciency of anti-erosion measures in the ﬁeld or to verify the results of different management practices over time.


Introduction
Europe's agricultural land, both active and in disuse, is substantially affected by erosion [1][2][3][4][5] from two main sources: water and tillage [6,7]. Water erosion is usually observed in lower, concave slope positions following the drainage network, while tillage erosion tends to be found on higher, convex slope positions following the direction of tillage. The rate of tillage erosion may be much greater than water erosion [6,7].
Initial research into erosion measurement tools focused on delocalized quantification, such as total material lost by basins [8], but without mapping the real places inside the basin where the erosion took effect. Subsequent mathematical models based on geographical data were developed to make quantitative assessments of local erosion [9][10][11], but these methods still require field measurements [12][13][14] for calibration and commonly use data from other models to simulate rainfall [15,16]. Chemical soil-particle tracing of soil movement due to erosion using multispectral remote sensing data was also proposed [17] but had availability and price drawbacks. In the search for a methodology without these constraints, airborne, full-waveform LIDAR systems offered remote sensing data to model erosion processes but was expensive and had limited coverage and low visit frequency [18].
The Sentinel-1 mission (a joint venture between the European Space Agency and the European Commission) was an opportunity to study localized erosion by using differential interferometry through its TopSAR images, which are frequently produced and freely distributed, and have a high resolution and broad coverage [19,20]. Interferometry uses SAR images of the same area to generate interferograms in which the information is closely related to the topography of the ground [21] and can be used to represent deformation in the form of maps [22]. Differential interferometry has been used in the past to study the localized effect on the Earth's surface of natural phenomena, such as volcanic eruptions [23,24], earthquakes [25], subsidence [26], slope movements [27], landslides [28,29], and even a "jökulhlaup" (a type of glacial outburst flood) [30], as well as human-induced effects, such as coal mining [31], fracking [32], and engineering works [33,34].
The present work describes the potential of differential interferometry for water and tillage surface erosion analysis. The methodology previously developed by the authors [35] was tested in a small plot in El Molar (north of Madrid, Spain) comparing ground deformation measurements made by differential interferometry and the gold standard laser-based terrestrial LIDAR (Light Detection and Ranging) system. This work demonstrated that differential interferometry measures are sufficiently similar to LIDAR measures (for 38 out of a total 40 observations, the difference between both methods was within the range of the maximum possible error associated with the LIDAR system, and even in the remaining two, the range was only exceeded by 0.0003 m and 0.0049 m) to make them useful in larger scale studies [36]. The methodology proposed in this work has been shown to be able to monitor spatially distributed erosion, thus enhancing the understanding of the water and tillage erosion processes involved. The results would increase the capacity for undertaking corrective or protective measures [37] since erosion and sedimentation processes could be modelled on a grid perspective at the scale of just a few square metres or days, regardless of weather or sunlight conditions [19,38].

Study Area
The study area is a small (7.05 km 2 ), monitored agricultural basin in Daganzo de Arriba, located slightly to the north of the city of Madrid, Spain, with centroid coordinates EPSG25830: 461.659; 4.494.032 ( Figure 1). proposed [17] but had availability and price drawbacks. In the search for a methodology without these constraints, airborne, full-waveform LIDAR systems offered remote sensing data to model erosion processes but was expensive and had limited coverage and low visit frequency [18]. The Sentinel-1 mission (a joint venture between the European Space Agency and the European Commission) was an opportunity to study localized erosion by using differential interferometry through its TopSAR images, which are frequently produced and freely distributed, and have a high resolution and broad coverage [19,20]. Interferometry uses SAR images of the same area to generate interferograms in which the information is closely related to the topography of the ground [21] and can be used to represent deformation in the form of maps [22]. Differential interferometry has been used in the past to study the localized effect on the Earth's surface of natural phenomena, such as volcanic eruptions [23,24], earthquakes [25], subsidence [26], slope movements [27], landslides [28,29], and even a "jökulhlaup" (a type of glacial outburst flood) [30], as well as human-induced effects, such as coal mining [31], fracking [32], and engineering works [33,34].
The present work describes the potential of differential interferometry for water and tillage surface erosion analysis. The methodology previously developed by the authors [35] was tested in a small plot in El Molar (north of Madrid, Spain) comparing ground deformation measurements made by differential interferometry and the gold standard laser-based terrestrial LIDAR (Light Detection and Ranging) system. This work demonstrated that differential interferometry measures are sufficiently similar to LIDAR measures (for 38 out of a total 40 observations, the difference between both methods was within the range of the maximum possible error associated with the LIDAR system, and even in the remaining two, the range was only exceeded by 0.0003 m and 0.0049 m) to make them useful in larger scale studies [36]. The methodology proposed in this work has been shown to be able to monitor spatially distributed erosion, thus enhancing the understanding of the water and tillage erosion processes involved. The results would increase the capacity for undertaking corrective or protective measures [37] since erosion and sedimentation processes could be modelled on a grid perspective at the scale of just a few square metres or days, regardless of weather or sunlight conditions [19,38].

Study Area
The study area is a small (7.05 km 2 ), monitored agricultural basin in Daganzo de Arriba, located slightly to the north of the city of Madrid, Spain, with centroid coordinates EPSG25830: 461.659; 4.494.032 ( Figure 1).  It is an elongated basin with topographic heights ranging between 689.23 and 789.70 m above sea level. The basin has slopes of less than 10%, except near the border of the watershed, and a main drainage network developed in the shape of a fishbone, with a secondary network heavily anthropized by the effect of agriculture, with a predominant arrangement of parallel channels following the direction of the maximum slope, EW, and WE ( Figure 2). It is an elongated basin with topographic heights ranging between 689.23 and 789.70 metres above sea level. The basin has slopes of less than 10%, except near the border of the watershed, and a main drainage network developed in the shape of a fishbone, with a secondary network heavily anthropized by the effect of agriculture, with a predominant arrangement of parallel channels following the direction of the maximum slope, EW, and WE ( Figure 2). The soil coverage in the basin alternates between cereal crops and fallow land due to its agricultural use, so no interference from vegetation is expected from C-band microwave measurements [39].

Scenarios Studied
This study includes five scenarios: three to analyse tillage erosion and two to analyse water erosion. It should be noted that the NDVI (Normalized Difference Vegetation Index) is obtained from Sentinel-2 images whose acquisition dates are different from the Sentinel-1 images, so dates for broader periods were selected.

Agricultural Tasks without Storm Events
The influence of the agricultural tasks of tillage and harvesting on the vertical deformation measures calculated by differential interferometry over TopSAR images was studied in three different scenarios, each consisting of 12 days: Scenario AW-1: During this period (12 days, from 15 June 2017 to 27 June 2017), there was almost no rainfall (3.6 mm in 12 days), but harvesting was carried out in multiple plots. The NDVI acquisition period was from 11 June 2017 to 1 July 2017. The soil coverage in the basin alternates between cereal crops and fallow land due to its agricultural use, so no interference from vegetation is expected from C-band microwave measurements [39].

Scenarios Studied
This study includes five scenarios: three to analyse tillage erosion and two to analyse water erosion. It should be noted that the NDVI (Normalized Difference Vegetation Index) is obtained from Sentinel-2 images whose acquisition dates are different from the Sentinel-1 images, so dates for broader periods were selected.

Agricultural Tasks without Storm Events
The influence of the agricultural tasks of tillage and harvesting on the vertical deformation measures calculated by differential interferometry over TopSAR images was studied in three different scenarios, each consisting of 12 days: Scenario AW-1: During this period (12 days, from 15 June 2017 to 27 June 2017), there was almost no rainfall (3.6 mm in 12 days), but harvesting was carried out in multiple plots. The NDVI acquisition period was from 11 June 2017 to 1 July 2017.
Scenario AW-2: During this period (12 days, from 17 June 2019 to 29 June 2019), there was no rain, but harvesting was carried out in multiple plots. The NDVI acquisition period was from 16 June 2019 to 1 July 2019.
Scenario AW-3: During this period (12 days, from 29 June 2019 to 11 July 2019), there was almost no rainfall (less than 0.3 mm), but agricultural tasks, such as tillage, were The moisture conditions were similar for each pair of acquisitions in the three scenarios, considering the hours of sunlight and the lack of precipitation the days before each acquisition.

Significant Storm Events without Agricultural Tasks
Free from the influence of agricultural tasks, vertical deformations caused by storms were measured by differential interferometry over TopSAR images in two different scenarios. Both scenarios include large storm events and the absence of agricultural tasks but differ in the land coverage: one is well covered with cereal crops, whereas the other has large extensions of naked soil in fallow lands.
Scenario STM-1: During this period (6 days, from 28 February 2018 to 6 March 2018), there was significant rainfall with a total of 95.1 mm, and most plots had the protection of vegetation cover (mainly wheat and barley). Moisture conditions were quite different for both SAR acquisitions. The first acquisition came after a long period with numerous hours of sun (about 10 h/day) without precipitation, except the day before the acquisition, which was totally overcast and rainy, while the second acquisition occurred after seven days of precipitation and almost no insolation, and the soil was almost saturated and even waterlogged in some areas. This meant that the soil moisture at the second acquisition was much higher than at the first one, and soil moisture affects the soil-penetrating capacity of the C-band microwaves; so in this case a slight increase can be expected in elevation values in the second acquisition. The NDVI acquisition dates were 26 February 2018 and 18 March 2018.
Scenario STM-2: During this period (12 days, from 16 August 2019 to 28 August 2019), there was significant rainfall with a total of 44.4 mm, coinciding with a large number of fallow plots, some of them recently tilled and therefore lacking the protection of vegetation cover. Moisture conditions were quite similar for each acquisition: on the first acquisition day, the soil was dry after a long period without precipitation and high insolation values, and on the second acquisition day, the soil was slightly wetter due to rain two days earlier but followed by two days of high insolation that dried the soil. The NDVI acquisition dates were 15 and 30 August 2019.

The Sentinel-1 Mission
The Sentinel-1 mission involves two satellites with almost polar orbits with a phase shift of 180 • on the same plane. These satellites are equipped with microwave C-band SAR devices, allowing image acquisition both day and night regardless of the cloud conditions [19,38]. Each satellite has a re-visiting frequency of 12 days, which is reduced to six days when both satellites are considered.
The SM mode has a high resolution (5 m × 5 m) and good coverage (80 km swath) but very limited availability (currently only available for support for catastrophic events and for studies already under way). The EWS mode has very good coverage (400 km swath) and availability but only medium resolution (20 m × 40 m). The W mode is a SM variation with better availability but poorer coverage. The IW mode has a medium-high resolution (5 m × 20 m), good coverage (250 km swath), and very good availability. The IW mode information stream was therefore adopted for the present work.
The IW mode uses the TOPSAR (Terrain Observation with Progressive Scans SAR) technique to capture images in the form of three sub-swaths, each composed of several radar bursts [40]; this improves signal-to-noise (SNR) and distributed target ambiguity (DTAR) ratios thanks to the electronic control and synchronization of bursts and passes and produces images with less noise and fewer deformations [41].

Data Acquisition
The data used in this research consist of downloaded TopSAR images for conducting the differential interferometric analysis and data gathered on inspection visits during on-site field visits to the small basin studied in this research.
TopSAR images were downloaded from the Open Hub of the ESA (https://scihub. copernicus.eu/, accessed on 16 October 2021) just a few hours after they were obtained. In order to minimize quality loss during simultaneous registration, it is important that the images in the interferogram should be taken in the same ascending or descending orbit and overlapped as much as possible.
The images used were captured by the Sentinel 1A satellite with IW mode single-look captures. The preferred polarization considered was VV, which minimizes the effect of low-growing vegetation [42,43]. Table 1 contains the TopSAR images used in the present work. Table 1. The ascending orbit S-1A satellite TopSAR images used in the present work.
TopSAR Image S1A_IW_SLC__1SDV_20170615T181057_20170615T181125_017048_01C665_D0FE S1A_IW_SLC__1SDV_20170627T181058_20170627T181126_017223_01CBB8_D37D S1A_IW_SLC__1SDV_20190617T181110_20190617T181138_027723_032115_1922 S1A_IW_SLC__1SDV_20190629T181111_20190629T181139_027898_03264F_2B36 S1A_IW_SLC__1SDV_20190711T181111_20190711T181139_028073_032B9D_2FF8 S1A_IW_SLC__1SDV_20180228T061751_20180228T061819_020803_023AB2_8227 S1A_IW_SLC__1SDV_20180306T181059_20180306T181127_020898_023D9D_F125 S1A_IW_SLC__1SDV_20190816T181114_20190816T181142_028598_033C18_9BA7 S1A_IW_SLC__1SDV_20190828T181114_20190828T181142_028773_034240_9BA9 Before and after each significant event studied (storm, harvesting, or tillage), a field visit was made to verify the following: • Photography: each plot was photographed to record its condition.

•
Erosion indicators, such as the presence of gullies or formerly buried stones.

•
Sedimentation indicators, such as the existence of sediment accumulations. • Vegetation coverage: on each visit, the plots were marked as fallow land, fallow land with remnants of crop harvesting, cereal crops, or forest. • Agricultural tasks: evidence of tillage or harvesting, such as the presence of working agricultural machinery, recently tilled plots, or recently harvested crops.

•
Other interferences, such as the presence of flooded areas, stacks of straw bales, grazing animals, etc.
On-site field visits were supplemented with vegetation index NDVI obtained from Sentinel 2 satellites in order to evaluate the vegetation status of each plot and any changes caused by agricultural activities [44][45][46]. Table 2 contains the Sentinel-2 images used in the present work.

NDVI Image Satellite
S-2A

Differential Interferometry
The interferometry technique is based on accurately measuring the distance travelled by electromagnetic radiation between a SAR radar antenna and the Earth's surface and back, and the intensity and phase of the returning radiation is recorded in SAR images. By comparing two images of the same area of the Earth's surface, it is possible to make an interferogram, in which the information on the phase difference is closely related to the terrain topography, and its deformations can be graphically represented through maps.
SAR sensors accurately measure the intensity and phase of electromagnetic radiation during their journey (Figure 3), which is recorded in a SAR image by means of a twodimensional mosaic of elements named pixels, each one referring to a small area of the Earth's surface. Each pixel contains information on the intensity and phase of the reflected radiation for each cell of surface area. S-2A

Differential Interferometry
The interferometry technique is based on accurately measuring the distance travelled by electromagnetic radiation between a SAR radar antenna and the Earth's surface and back, and the intensity and phase of the returning radiation is recorded in SAR images. By comparing two images of the same area of the Earth's surface, it is possible to make an interferogram, in which the information on the phase difference is closely related to the terrain topography, and its deformations can be graphically represented through maps.
SAR sensors accurately measure the intensity and phase of electromagnetic radiation during their journey (Figure 3), which is recorded in a SAR image by means of a twodimensional mosaic of elements named pixels, each one referring to a small area of the Earth's surface. Each pixel contains information on the intensity and phase of the reflected radiation for each cell of surface area.
According to Equation (1), the phase measurement depends on the wavelength (λ) of the electromagnetic signal and the distance travelled (R).
SAR interferometry can determine the distance to the ground ( Figure 4) by using the phase difference between two SAR observations of the same cell taken from slightly different positions [47].
According to Equation (1), the phase measurement depends on the wavelength (λ) of the electromagnetic signal and the distance travelled (R).
SAR interferometry can determine the distance to the ground ( Figure 4) by using the phase difference between two SAR observations of the same cell taken from slightly different positions [47].
According to Equations (2)-(4), the phase difference between two acquisitions depends on the wavelength (λ) of the electromagnetic signal used and the difference between the distances that the signal travelled in each acquisition (∆R).
SAR interferometry requires the acquisition of both images and the subsequent construction of an interferogram by combining the phase differences, whose information is closely related to the topography of the terrain [21] and allows the mapping of surface deformations.
According to previous research [47], the phase difference obtained can be caused by five factors (Equation (5)). where: flat = the curvature of the Earth elevation = topography displacement = deformations of the terrain between the two information acquisitions atmosphere = differences in relative humidity, pressure, and temperature between information acquisitions noise = changes over time in volume scattering and acquisition angles, etc.
The interferometric process removes all these factors except those of interest: elevation (to generate digital models of the terrain) and displacement (to analyse surface changes). Therefore, if the difference in phase caused by the topography is removed from the interferogram, the resulting difference corresponds only to surface deformations between the acquisition dates.
A coherence band is obtained along with the difference phase band and is calculated within the interferogram to estimate the similarity between both images at a pixel scale between 0 (no similarity) and 1 (total similarity). According to previous researches [48,49], the coherence loss between a pair of images can be explained through Equation (6). where: γ T = Temporal factor: cannot be avoided and is due to differences on the ground occurring between two information acquisitions; this is the object of the present work. γ G = Geometric factor: caused by errors in the orbits followed by satellites; can be partially removed from calculations. γ V = Volumetric factor: the presence of vegetation causes the appearance of the volumetric factor, which cannot be removed. γ P = Processing factor: caused by calculation errors, which must be minimized and avoided at all costs.

Software
The ESA's Sentinel Application Platform (SNAP) was used for the differential interferometric analysis of pairs of TopSAR images [50,51], applying the SNAPHU algorithm (Copyright 2002-2020 Board of Trustees, Leland Stanford Jr. University) produced at Stanford University for the unwrapping phase [52][53][54] inside a CYGWIN (Copyright © Cygwin authors) environment (a collection of tools which provide functionality like a Linux dis- tribution on Windows). The ESA's Sentinel application platform (SNAP) was also used to calculate the NDVI index, and the differences in the NDVI index between days were computed with QGIS software.

Interferometry Calculation Steps
To build the interferograms, interferometric pairs of SAR images were coregistered, aligning the images at the pixel level and using the older one as the image reference. This coregistration was done using 21-point Bisinc interpolation, enhanced spectral display, and orbital correction (with the Sentinel Precise polynomial degree 3 option) [55].
Once the images were aligned, their interferograms (intensity bands, phase, and coherence) were calculated and corrected by considering the Earth's curvature (using a fifth-degree polynomial) and performing an interpolation of the capturing satellite's orbit (using a third-degree polynomial).
One peculiarity of TopSAR images is that they are acquired by radar bursts [41], so the spaces between the bursts were eliminated to advance the interferometric analysis [40]. The image was then cropped to the study area, which drastically reduces the time required for the subsequent operations.
After obtaining the interferogram, the dephasing due to the topography was eliminated using an SRTM 1sec HGT digital terrain model. Then, after reducing the noise effect by Goldstein phase filtering [56], the resulting phase difference corresponds to the deformation and atmosphere differences occurring between data acquisitions [21]. The coherence band was obtained along with the difference phase band, and its information is used in the final step to make a vertical adjustments to the deformation raster to compensate, at least partially, for the atmosphere factor.
The interferograms were unwrapped with the SNAPHU algorithm in order to calculate the vertical displacements [52][53][54].
The vertical displacement in each pixel was calculated using Equation (7) [57].
where: d = the deformation according to the vertical axis λ = the wavelength used ∆φ d = the phase displacement for each cell between data acquisitions θ inc = incident angle Differential interferometry can obtain deformation only in the direction of line of sight (LOS), so vertical deformation is estimated as a vertical projection of the LOS in the vertical direction [58].
The terrain was geometrically corrected to compensate for the distortions in the SAR image due to variations in the terrain with respect to the inclination of the sensor. The aim was to produce a two-dimensional representation that faithfully reflected the realworld situation. This was done using a Range-Doppler orthorectification method [59] to geodecodify the SAR images, considering information for the capturing satellite's orbit, time records, conversion factors for the oblique distance to the Earth's surface, and a reference digital terrain model to determine the exact location represented by each pixel.
This geometric correction was done using the ETRS89 projection system, which produced two raster products: the terrain deformation between data acquisitions and the coherence between the acquisitions.
Finally, the coherence band was used to make a vertical adjustment required to compensate vertical error in the estimation of local deformations. As the flatness, noise, and elevation factors were previously removed, the remaining errors were due to atmosphere factor [47]. Representative cells with very high coherence, hence with little likelihood of being deformed, were chosen to provide the zero-deformation reference. The value of its twin cells in the deformation raster were then used to create an adjustment raster by interpolating its values [30,36,60]. Finally, the adjustment raster was subtracted from the deformation raster to obtain the adjusted deformation raster.

Normalized Difference Vegetation Index
The data for the calculation of the NDVI were obtained from the ESA's Sentinel-2 satellite platform. It is important to note that the satellite overpass calendar for Sentinel-2 is different from Sentinel-1, which is the reason the NDVI and interferometry dates do not match exactly. Sentinel-2 is also weather dependent, so images collected on cloudy days were rejected.
If the exact date of the Sentinel-1 overpass was not available, another date was selected in order to cover a longer period, seeking an earlier date for the first acquisition in the interferogram pair and a subsequent date for the second one.
The NDVI algorithm shows the strength and vitality of the vegetation on the Earth's surface [44][45][46]61,62] and is obtained through Equation (8) by using the reflection in the near-infrared spectrum (NIR) obtained by Sentinel-2 band B8 (842 nm) and the reflection in the red range of the spectrum (RED) obtained by Sentinel-2 band B4 (665 nm) The difference in the NDVI index between two days is calculated by simply subtracting the raster values of the second day from the raster values of the first day using QGIS software. The NDVI index has values between −1.0 and 1.0 [44][45][46]63], with negative values corresponding to water or snow areas. The

Agricultural Tasks without Storm Events
In NDVI acquisition (Figure 5a,b), areas shaded yellow and light green represent ripe cereal crops ready to harvest, while red and orange represent fallow land. Plots that change from yellow-green ( Figure 5a) to red-orange ( Figure 5b) indicate plots that were harvested between both acquisition dates. Figure 5c shows the variation in the NDVI between the acquisition dates, where green denotes crop growth, orange-red represents harvesting tasks, and yellow and yellow-orange represent no change or crop ripening.
As can be seen from the NDVI variations ( Figure 5), harvesting was carried out in multiple plots in scenarios AW-1 and AW-2 (verified during field inspection).
Tillage was observed in scenario AW-3 but no harvesting. As tillage has no effect on the NDVI index, it shows no changes.
Coherence values in scenario AW-3 ( Figure 6) were somewhat poorer than in previous scenarios owing to the agricultural tasks (mainly tillage) that affected the surface of the land (plots with a dominance of orange and red). high coherence values correspond to non-tilled plots, and low coherence values correspond to tilled plots. This divergence produces a greater standard deviation and a wider range of values, from 0.2 to 0.7 ( Figure 6).   High coherence (value > 0.95) points were selected for each scenario both inside and around the basin, and their matching deformation values were used to calculate the vertical adjustment raster by means of interpolation (Figure 7) to compensate for noise and atmospheric artifacts. The correction raster was subtracted from the deformation raster to obtain the corrected deformation raster for each scenario (Figure 8).  Table 3 shows the total deformation range inside the basin-that is, the difference between the maximum positive and negative deformations-and the mean deformations inside the basin (with the standard deviation) for scenarios AW-1 to AW-3. There was no rainfall during the period in scenario AW-3 (29 June 2019 to 11 July 2019), but there were tillage works in several plots. As expected, there was no basin level erosion, as could be verified from the mean value of the basin deformation that tended to zero (0.00021 m), although high local deformation values were detected inside the plots with tillage works performed during scenario AW-3, ranging from −6.3 cm to +6.7 cm. These local deformation values correspond to tillage erosion and sedimentation processes limited to the interior of each plot (Figure 8 AW-3).

Significant Storm Events without Agricultural Tasks
NDVI STM-1: As can be seen in Figure 9c, NDVI differences were positive (green) according to the vegetation growth and absence of harvesting, except for two plots to the south and one plot to the north. Figure 9a,b both show a large number of plots with a predominance of green shades, implying that they had a cereal crop coverage. Some plots The corrected deformation raster for scenarios AW-1, AW-2, and AW3 and its histogram (Figure 8) show the following:

•
In harvesting scenarios (AW-1 and AW-2), no deformation differences are found between harvested and non-harvested plots.

•
In the tillage scenario (AW-3), no deformation is observed (mainly in yellow) in nontilled plots, while a mix of positive (green) and negative (red) deformations of the Earth's surface is observed in tilled plots. Table 3 shows the total deformation range inside the basin-that is, the difference between the maximum positive and negative deformations-and the mean deformations inside the basin (with the standard deviation) for scenarios AW-1 to AW-3. However, during this period, there were major changes in the height of cereal crops (harvesting tasks) that were not detected as deformation by differential interferometry over C-Band microwaves of TopSAR images and thus were not reflected in the deformation raster.
There was no rainfall during the period in scenario AW-3 (29 June 2019 to 11 July 2019), but there were tillage works in several plots. As expected, there was no basin level erosion, as could be verified from the mean value of the basin deformation that tended to zero (0.00021 m), although high local deformation values were detected inside the plots with tillage works performed during scenario AW-3, ranging from −6.3 cm to +6.7 cm. These local deformation values correspond to tillage erosion and sedimentation processes limited to the interior of each plot (Figure 8 AW-3).

Significant Storm Events without Agricultural Tasks
NDVI STM-1: As can be seen in Figure 9c, NDVI differences were positive (green) according to the vegetation growth and absence of harvesting, except for two plots to the south and one plot to the north. Figure 9a,b both show a large number of plots with a predominance of green shades, implying that they had a cereal crop coverage. Some plots show a predominance of orange or yellow shades, so they lack vegetation coverage, as these are fallow lands with bare soil or with the remnants of crop harvesting.
NDVI STM-2: Figure 9a,b both show most plots with predominance of orange, indicating a lack of vegetation coverage, and are fallow lands with bare soil or with the remnants of crop harvesting. Only a few plots are green, implying that they had some vegetation coverage. These plots are also green in Figure 9c, where crops are still growing. The field visit confirmed that the coverage of these green NDVI plots in Figure 9c were weeds growing between the remnants of the harvested cereal crops. The rest of the plots are coloured yellow in the image on the right, so they remained fallow during the study period. Owing to the erosion expected due to the heavy rainfall, the coherence values of both scenarios STM-1 and STM-2 ( Figure 10) were low, as erosion corresponds to the temporal factor, an unavoidable component of coherence loss. Figure 9. NDVI at the beginning and end of AW-1, 2, and 3 (a,b) and NDVI variation (c). Zebra frame segments of 1 km.
Owing to the erosion expected due to the heavy rainfall, the coherence values of both scenarios STM-1 and STM-2 ( Figure 10) were low, as erosion corresponds to the temporal factor, an unavoidable component of coherence loss.
A different coherence behaviour is observed between plots with and without vegetal coverage in scenario STM-1. As can be seen, the greener areas in Figure 9a,b (higher NDVI index and more developed vegetation coverage) correspond to the highest coherence values in Figure 10, while the red-orange areas in Figure 9a A different coherence behaviour is observed between plots with and without vegetal coverage in scenario STM-1. As can be seen, the greener areas in Figure 9a,b (higher NDVI index and more developed vegetation coverage) correspond to the highest coherence values in Figure 10, while the red-orange areas in Figure 9a,b (lower NDVI index, fallow land) correspond to the lowest coherence values in Figure 10. The basin in scenario STM-1 had a mean coherence value of 0.437, a standard deviation of 0.174, and a maximum of 0.950.
A different coherence behaviour ( Figure 10) is observed in scenario STM-2 between plots with a similar low NDVI value corresponding to fallow land (Figure 9a,b), so these plots should have a similar lack of vegetation coverage. However, two kinds of fallow land were registered on field visits: recently tilled fallow land with bare soil (see the agricultural tillage tasks detected in scenario AW-3, agricultural tasks without storm events, from 29 June 2019 to 11 July 2019), with low coherence values (predominance of red-orange), and fallow land without tillage and with remnants of crop harvesting, with higher coherence values (predominance of yellow-green). The basin in scenario STM-2 had a mean coherence value of 0.401, a standard deviation of 0.178, and a maximum of 0.889. A different coherence behaviour ( Figure 10) is observed in scenario STM-2 between plots with a similar low NDVI value corresponding to fallow land (Figure 9a,b), so these plots should have a similar lack of vegetation coverage. However, two kinds of fallow land were registered on field visits: recently tilled fallow land with bare soil (see the agricultural tillage tasks detected in scenario AW-3, agricultural tasks without storm events, from 29 June 2019 to 11 July 2019), with low coherence values (predominance of red-orange), and fallow land without tillage and with remnants of crop harvesting, with higher coherence values (predominance of yellow-green). The basin in scenario STM-2 had a mean coherence value of 0.401, a standard deviation of 0.178, and a maximum of 0.889.
Both scenarios showed fewer high-coherence (value > 0.95) points for selection than in the periods without heavy rain but enough to interpolate the vertical adjustment raster (Figure 11). Agronomy 2021, 11, 2075 18 Both scenarios showed fewer high-coherence (value > 0.95) points for selection in the periods without heavy rain but enough to interpolate the vertical adjustment r ( Figure 11).

STM-1
STM-2 Figure 11. Matching deformation values to high-coherence points and its interpolation. zebra frame segments of 1 km.
The corrected deformation raster is obtained for each scenario after subtracting correction raster from the deformation raster ( Figure 12). Table 4 shows the deformation range and the mean deformations inside the basin (with the stan deviation) for scenarios STM-1 and STM-2. The corrected deformation raster is obtained for each scenario after subtracting the correction raster from the deformation raster ( Figure 12). Table 4 shows the total deformation range and the mean deformations inside the basin (with the standard deviation) for scenarios STM-1 and STM-2. From the comparison between the deformation raster for scenario STM-1 (Figure 12, STM-1) and the information collected in field visits: • Dark green areas: visual inspection of the basin showed that in this case, dark green areas corresponded to sediment (max. 4.4 cm) and flooded areas (Images 2-3 in Figure 13). • Light green areas: this area requires interpretation. There are extensive light green areas that cannot correspond to sedimentation areas (Image 4 in Figure 13). It is likely that the slight elevation recorded is due to clay expansion phenomena or changes in the dielectric constant of the soil that cause a more superficial rebound of the C-band microwaves rather than sediment phenomena. • Yellow areas: no erosion, mainly plots with a coverage of cereal crops.

•
Orange to red areas: areas with moderate to high erosion (max. 4.03 cm), forming even rills and gullies in the reddest areas, as observed in the field. Mainly fallow lands without vegetation coverage (Image 1 in Figure 13). From the comparison between the deformation raster for scenario STM-1 (Figure 12, STM-1) and the information collected in field visits: • Dark green areas: visual inspection of the basin showed that in this case, dark green areas corresponded to sediment (max. 4.4 cm) and flooded areas (Images 2-3 in Figure 13). • Light green areas: this area requires interpretation. There are extensive light green areas that cannot correspond to sedimentation areas (Image 4 in Figure 13). It is likely that the slight elevation recorded is due to clay expansion phenomena or changes in the dielectric constant of the soil that cause a more superficial rebound of the C-band microwaves rather than sediment phenomena. From the comparison between the STM-2 corrected deformation raster ( Figure 12) and the information collected on field visits:

•
Dark green areas: visual inspection of the basin showed that in this case, the dark green areas correspond to sediment areas with a maximum of 4.87 cm (Image 2 in Figure 14).
• Yellow areas: no erosion, mainly fallow land with remnants of harvested cereal crops (Image 3 and the right part of Image 4 in Figure 14). • Red and orange areas: areas with moderate to high erosion. Orange areas (from 0.25 to 1 cm) predominate over the fewer and scattered reddish areas (from 1 to 5 cm) ( Figure 15). Mainly bare soil on recently tilled fallow lands without vegetation coverage (Image 1 and the left part of Image 4 in Figure 14), forming rills and gullies.

•
Purple areas: areas with very high erosion (from 5 to 6.21 cm). There are only a few scattered points with no relevance in the global erosion values (Figure 15).
Some plots were tilled about three weeks before the storms in STM-2 and did not show the usual water erosion patterns ( Figure 15) but a combined effect of tillage and water erosion. In addition to the usual pattern of water erosion following the drainage network, another pattern can be observed showing high erosion values in convex areas on low slopes (sub-watershed border areas), characteristic of tillage erosion.
Hence, even if no tillage works were done in the STM-2 period, the effects of the previous tillage became evident with the first storm and increased the erosion rates in zones where no water erosion was expected. From the comparison between the STM-2 corrected deformation raster ( Figure 12) and the information collected on field visits: • Dark green areas: visual inspection of the basin showed that in this case, the dark green areas correspond to sediment areas with a maximum of 4.87 cm (Image 2 in Figure 14).

•
Yellow areas: no erosion, mainly fallow land with remnants of harvested cereal crops (Image 3 and the right part of Image 4 in Figure 14).

•
Red and orange areas: areas with moderate to high erosion. Orange areas (from 0.25 Some plots were tilled about three weeks before the storms in STM-2 and did not show the usual water erosion patterns ( Figure 15) but a combined effect of tillage and water erosion. In addition to the usual pattern of water erosion following the drainage network, another pattern can be observed showing high erosion values in convex areas on low slopes (sub-watershed border areas), characteristic of tillage erosion.
Hence, even if no tillage works were done in the STM-2 period, the effects of the previous tillage became evident with the first storm and increased the erosion rates in zones where no water erosion was expected.

Discussion
Each of the scenarios studied shows the different capabilities and limitations of the present methodology: Scenarios AW-1 and AW-2: Scenarios without soil works and without rainfall but

Discussion
Each of the scenarios studied shows the different capabilities and limitations of the present methodology: Scenarios AW-1 and AW-2: Scenarios without soil works and without rainfall but with a major change in the height of the cereal crops due to harvesting. No deformations were recorded in the terrain surface, as was expected from the lack of tillage or water erosion, so differential interferometry is not affected by vegetation, such as winter cereal crops, because of the wavelength of the C-Band microwaves used.
Scenario AW-3: Scenario with almost no rain, without harvesting, but with tillage tasks. Deformations were recorded only in the tilled plots, while the global deformation of the basin was near zero, so differential interferometry was able to detect tillage erosion. In this case, coherence loss was considerable, so a clear temporal limitation was detected, as coherence loss in longer periods could be a problem.
Scenario STM-1: Scenario without agricultural tasks, with high rainfall, and with soil coverage in most of the plots. Differential interferometry recorded deformations as expected from water erosion, which was greater in plots without soil coverage.
Scenario STM-2: Scenario without agricultural tasks but with recently tilled plots, high rainfall, and without soil coverage in most of the plots. Deformations recorded with differential interferometry showed that even with lower total rainfall values than STM-1, the erosion values in STM-2 were greater, as the rainfall in STM-1 was spread over several days, and the vegetation coverage exerted a defensive effect. In contrast, most of the rainfall in STM-2 was concentrated in a single day, there was no vegetation cover to serve as a defence against erosion, and the terrain in some plots had been recently tilled, increasing the effect of water erosion. It is worth noting the synergistic effect between water erosion and previously tilled plots, with a much greater erosion value recorded than in nearby non-tilled plots, and the effect of crop remnants as soil protection against erosion.
Scenarios STM-1 and STM-2 show that differential interferometry is capable of detecting water erosion and even the synergy between water erosion and recent tillage and can be used to measure the effects of different factors, such as soil coverage. As in scenario AW-3, coherence loss was considerable, highlighting the temporal limitation of this methodology, so coherence loss over a longer period could be problematic.
In any case, the temporal limitation due to coherence loss could be solved by studying longer periods as a sum of shorter periods as long as coherence loss is limited in each of the shorter periods.
Another factor to take into account is the differences in soil moisture between acquisitions. In scenario STM-1, the moisture conditions varied widely between acquisitions and were quite dry at first and saturated and even waterlogged in some areas in the second. Moisture conditions have an effect on acquisition (the wetter the soil, the lower the penetrating capabilities of the C-Band microwaves), so in this case, erosion may have been underestimated in some areas. This is the most likely the reason for the light green areas in Figure 12 with a slight elevation recorded that does not correspond to sedimentation processes and may explain why the mean value of basin deformation is almost zero and not a negative value. To eradicate this effect, a longer time series should be studied to obtain the effect of moisture on measurements in this basin, or-even better-high-coherence points should be provided and used in the affected zones to measure their distortion in the deformation band.
The use of the NDVI combined with the coherence band allowed the researchers (verified with fieldwork) to identify the events related to the erosion/sedimentation in each scenario (Table 5).

Conclusions
The results of this work show that erosion monitoring in agricultural basins is possible with differential interferometry over Sentinel-1 TopSAR images. Its capabilities have been successfully tested in different conditions related to agricultural tasks without precipitation and storm events.
Two different scenarios were analysed in the case of agricultural tasks without precipitation: harvesting and tillage. In harvesting scenarios without precipitation, no erosion or sedimentation was detected, so changes in cereal crop height do not significantly modify vertical deformation measurements, and there is no evidence of erosion processes in this scenario. For the tillage scenario without rain, local deformations were detected inside plots, but there was no global basin deformation, which is consistent with the tillage erosion expected.
Two different time periods, including storm events with heavy rain, were analysed. The analyses conducted using satellite images during both periods showed erosion and sedimentation processes. These results agree with the information gathered during the field visits, which revealed evidence, such as rills and gullies (erosion), sedimentation areas, and vegetation coverage, in areas with less erosion measurement (even the remnants of cereal crops after harvesting demonstrated protection capabilities against erosion). Plots previously tilled before STM-2 deserve a special mention. Without a direct tillage erosion measurement (no tilling was done during the period but some weeks before), an increase was observed in the expected water erosion due to the combined effect of the previously tilled terrain with the first storm after tilling.
The output erosion raster has a cell size of 13.93 × 13.93 m, a good resolution for basin scale erosion studies, and sentinel satellites have a frequency that provides images of a given zone every six days (12 if using the same satellite to avoid system differences). This, combined with the free availability of Sentinel-1 TopSAR images, points to differential interferometry as a powerful tool for the analysis and management of erosion in agricultural basins.
It not only makes it possible to study erosion and sedimentation processes but even to test the efficiency of anti-erosion measures in the field and verify the results of different management practices over time. However, some issues still require further implementation and consideration, such as the influence of the soil moisture condition between both acquisitions.
The capacity of differential interferometry for monitoring basin erosion over Sentinel-1 TopSAR images, combined with the interpretation of its coherence band and the NDVI index obtained from Sentinel-2, drastically minimizes the field visits required. It also allows large areas to be monitored with minimum effort and only one qualified professional to interpret the combined information from the deformation/coherence/NDVI raster and distinguish between water and tillage erosion/sedimentation processes from other processes, such as subsidence, landslides, or any other event capable of deforming the Earth's surface.