Next Article in Journal
Handling Missing Data in Large-Scale MODIS AOD Products Using a Two-Step Model
Next Article in Special Issue
A Review of Remote Sensing Challenges for Food Security with Respect to Salinity and Drought Threats
Previous Article in Journal
First Evidences of Ionospheric Plasma Depletions Observations Using GNSS-R Data from CYGNSS
Previous Article in Special Issue
Improving Soil Thickness Estimations Based on Multiple Environmental Variables with Stacking Ensemble Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Letter

Separability of Mowing and Ploughing Events on Short Temporal Baseline Sentinel-1 Coherence Time Series

1
KappaZeta Ltd., 51007 Tartu, Estonia
2
Tartu Observatory, University of Tartu, Tõravere, 61602 Tartumaa, Estonia
3
University Medical Center Groningen, Department of Epidemiology, University of Groningen, 9713GZ Groningen, The Netherlands
4
Institute of Ecology and Earth Sciences, University of Tartu, 51003 Tartu, Estonia
5
Software Technology and Applications Competence Center, 51009 Tartu, Estonia
6
Flow Traders, 1018LL Amsterdam, The Netherlands
7
Institute for Environmental Solutions, Lidlauks, LV-4101 Priekuļu parish, Latvia
8
School of Electrical Engineering, Aalto University, 02150 Espoo, Finland
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(22), 3784; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223784
Submission received: 21 October 2020 / Revised: 13 November 2020 / Accepted: 15 November 2020 / Published: 18 November 2020

Abstract

:
Short temporal baseline regular Synthetic Aperture Radar (SAR) interferometry is a tool well suited for wide area monitoring of agricultural activities, urgently needed in European Union Common Agricultural Policy (CAP) enforcement. In this study, we demonstrate and describe in detail, how mowing and ploughing events can be identified from Sentinel-1 6-day interferometric coherence time series. The study is based on a large dataset of 386 dual polarimetric Sentinel-1 VV/VH SAR and 351 Sentinel-2 optical images, and nearly 2000 documented mowing and ploughing events on more than 1000 parcels (average 10.6 ha, smallest 0.6 ha, largest 108.5 ha). Statistical analysis revealed that mowing and ploughing cause coherence to increase when compared to values before an event. In the case of mowing, the coherence increased from 0.18 to 0.35, while Sentinel-2 NDVI (indicating the amount of green chlorophyll containing biomass) at the same time decreased from 0.75 to 0.5. For mowing, there was virtually no difference between the polarisations. After ploughing, VV-coherence grew up to 0.65 and VH-coherence to 0.45, while NDVI was around 0.2 at the same time. Before ploughing, both coherence and NDVI values were very variable, determined by the agricultural management practices of the parcel. Results presented here can be used for planning further studies and developing mowing and ploughing detection algorithms based on Sentinel-1 data. Besides CAP enforcement, the results are also useful for food security and land use change detection applications.

Graphical Abstract

1. Introduction

Grasslands constitute at least one third of the utilised agricultural area in the European Union (EU) [1]. Besides its agricultural value, grasslands support biodiversity and tourism, and provide protection from soil erosion [2]. Protection of grasslands is one of the goals of the EU Common Agricultural Policy (CAP), by encouraging regular mowing and prohibiting ploughing of permanent grasslands, but allowing it for agricultural grassland. Given the agenda to introduce satellite monitoring for subsidy checks, remote sensing will play an increasing role in the context of CAP [3].
Sentinel-1 (S1) is well suited for a potential grassland monitoring service due to its frequent revisit times, close to all-weather operations, and free data policy. Sentinel-1 constellation consists of two satellites—S1A and S1B, and is tailored for frequent interferometric capability. This allows to produce coherence maps on a wide scale, several times per week. Attempts to relate coherence to mowing and ploughing events on grasslands, are relatively recent, but based on well established studies. Changes in scatterer placement in soil and canopy lead to loss of coherence over time [4,5]. Low coherence can be an indication of presence of vegetation or farming activity, and high coherence can indicate sparse vegetation or lack of it due to, for example, mowing [6]. However, changes due to wind and precipitation, cause coherence to decrease over time, [5,7,8,9] making the interpretation ambiguous.
As earlier studies show, in the case of COSMO-Skymed interferometric pairs with a one-day temporal baseline, short vegetation after mowing results in high coherence [10]. A similar effect is present on TerraSAR-X Staring Spotlight coherence images with a temporal baseline of 11 days [11]. A comprehensive effort by [12] in relating Sentinel-1 12-day coherence to mowing events demonstrated the value of Sentinel-1 coherence for agricultural applications. Since then, grassland and other agricultural monitoring with Sentinel-1 SAR data has been a topic with increasing research interest. Applying artificial neural networks (ANN) on second order statistics of Sentinel-1 backscatter time series has shown reasonable accuracy for grassland mowing detection [13]. In another research, grassland mowing events were studied by modelling NDVI on a parcel level with recurrent neural networks, using Sentinel-1 backscatter and coherence time series as an input [14]. Sentinel-1 coherence time series have also been demonstrated to be useful for crop harvesting detection [15,16] and land cover mapping [17].
With the recent availability of S1B data, the experiment conducted in [12] can be extended by reducing the temporal baseline to six days, which would reduce the temporal decorrelation between the images [4,18]. This should lead to higher coherence values for post-mowing short vegetation conditions, thus enhancing the signatures of the farming events in the coherence time series. In addition, any potential operational agriculture monitoring system that uses coherence would benefit from the more frequent updates in the case of 6-day interferometric pairs.
This paper characterises the dynamics of Sentinel-1 6-day coherence with respect to mowing and ploughing events on agricultural grasslands in Estonia. Time series of coherence from more than a thousand grassland parcels (ranging from 0.6 to 108.5 ha, average 10.6 ha) are aggregated and presented in respect to the recorded agricultural management events. We show the increase of coherence after an event, and how it differs for mowing and ploughing. Known to the authors, it is the first study including such a large set of ground reference data, and has thus a significant value for conclusion and generalisation.

2. Data and Methods

2.1. Test Site and Farming Events

The study was carried out in Estonia, Northern Europe (between 57.5°N, 21.5°E and 59.8°N, 28.2°E). The study area is relatively flat with no steep slopes and altitudes ranging between 0 and 200 m above sea level (asl). Data about events were collected directly from the farmers and consisted of information about the start and end date of the activity and the area covered. In total, data about 1695 mowing and 195 ploughing events on 1159 and 177 unique agricultural parcels, respectively, were collected. The average area of mowing and ploughing events, respectively, is 10.4 and 12.2 ha, the minimum 0.6 and 0.9 ha, and the maximum 108.5 and 92.0 ha. Considering the main agricultural areas of the country (excluding forests and marshlands), events are geographically evenly distributed as seen in Figure 1. The earliest reported mowing event occurred on 20 May, while the latest on 12 October 2017. In the case of ploughing, the earliest event was recorded on 1 July and the latest on 8 September 2017. The dataset covers majority of the agricultural season in Estonia, where green vegetation emerges in the beginning of May and lasts until October. The dataset includes also information about the main soil type (peat or mineral soil).
Agricultural grasslands on mineral soils are major ecosystems, but some less intensively managed grassland parcels are also in wetlands on peat-soils. Main grass species are from Poaceae (e.g., reed canary grass and timothy) and Fabaceae (e.g., clover and alfalfa) families. There are parcels with pure species and also the mixtures of them. The fields are cultivated in spring or late autumn and mown one to three times per season. In most of the fields the grass is collected and used as animal fodder. Some parcels are also grazed, which is less common than mowing. Figure 2 presents typical Estonian Poaceae and Fabaceae agricultural grassland parcels before and after mowing. Parcels with Poaceae reed-like grasses have dominantly vertical structures. Parcels with Fabaceae clover-type grasses have also vertical structures from the stems, but the orientation of the canopy is more random.

2.2. Sentinel-1 Data

In total, 386 S1A/B SLC IW products acquired between 1 May and 30 October 2017, were processed. Local time of overpass was around 18:00 (15:00 GMT). A total of 87 products were from relative orbit number (RON)160, 62 from RON131, 84 from RON87, 93 from RON58, and 60 from RON29. These were organised into S1A/S1B 6-day pairs. Only nine products in the given time period were missing from the Copernicus Open Access Hub, slightly decreasing the data density for RON160 and RON87.
Besides different viewing geometries, the value of using multiple RONs is a shifted grid with independent temporal sampling. It should increase robustness to disturbances that are corrupting a single RON. The value of using multi-RON data is discussed in Section 3.3 below.
Incidence angles range between 29° and 46°. Spatial resolution ranges between 2.7 and 3.5 m in slant range, but is fixed at 22 m in azimuth.

2.3. Sentinel-2 Data

In total, 351 Sentinel-2A and -2B L2A products acquired between 1 May and 30 October 2017, were processed. Local time of overpass was around 13:00 (10:00 GMT). Each Sentinel-2 image is a maximum of three days off from the closest Sentinel-1 image. A total of 14 products were from RON22, 74 from RON36, 67 from RON79, 50 from RON93, 78 from RON122 and 68 from RON136. As the focus was on Sentinel-1 coherence analysis, only the NDVI [19], computed according to Equation (1), was derived from Sentinel-2 to have a comparison with a well-known optical vegetation remote sensing parameter.
NDVI = ( band 8 band 4 ) ( band 8 + band 4 )
where band8 is the 8th band (near infrared with central wavelength of 842 nm) and band4 is the 4th band (red with central wavelength of 665 nm) of Sentinel-2.

2.4. Ancillary Data

Weather and precipitation data were derived from Estonian Weather Service database. The weather information is collected from a network of weather stations over Estiona, which form a sufficiently dense grid, so that to every location in our field database, a weather station not farther than 40 km can be found. The database includes observation history for many years and is publicly available at Estonian Weather Service homepage [20].

2.5. Methods

The goal of data analysis was to characterize temporal changes in Sentinel-1 interferometric coherence caused by mowing and ploughing events on the fields and relate it also to the changes in optical images. For this, coherence time series were calculated about every field in the database in relation to the event. Average coherence of a field, imaging geometry parameters, imaging time and average NDVI was stored in a database. Analysis of the database of coherence values with different temporal distance from a farming event allowed to extract the typical behaviour by means of statistical analysis. The behaviour of Sentinel-2 NDVI was analysed in a similar way. The database formation process involved preprocessing a large amount of satellite images where average coherence and NDVI value was calculated for every parcel for every available date (constrained by image availability and cloud cover).
The Sentinel Application Platform (SNAP) toolbox was used for processing. The following steps were applied: apply orbit file, back-geocoding (using Shuttle Radar Topography Mission (SRTM) data), coherence calculation, deburst, terrain correction, and reprojection to the local projection (EPSG:3301). In the last step, data were resampled to 4 m resolution to maintain the maximum spatial resolution, and having square-shaped pixels. The terrain around the study areas is relatively flat, resulting in negligible topographic distortions of the SAR data. Coherence was calculated for each swath separately. Only pixels that were completely inside the parcel boundaries (also considering the averaging window used for coherence calculation) were used for calculation of results and to discard any interference beyond parcel borders. Coherence was computed pair-wise with 6-day time steps. For every coherence computation there was a new master and slave. The results were stored in a database following a forward-looking convention—coherence about date X means the coherence between Sentinel-1 images from date X and X + 6 days.
To understand how coherence ( γ ) varies in relation to mowing and ploughing events, γ from all parcels with recorded events were aggregated. All time series were shifted to the same grid with a step of one day by placing the event at time 0. The grid extends from 30 (i.e., 30 days before an event) to + 30 (Figure 3). The median value and third and first quartiles were calculated for each grid point. Ploughing signatures are presented in limited −10 to +30 day range, as the status of a grassland before a ploughing event was too variable for a meaningful generalisation.

2.6. Interferometric Coherence and Temporal Decorrelation

The amplitude of the complex correlation coefficient also known as coherence is calculated for two complex SAR acquisitions s 1 and s 2 , as
γ = | s 1 s 2 * | s 1 s 1 * s 2 s 2 * , 0 γ 1 .
Estimated coherence γ is a product of several terms specific to the measurement system, data processing steps and the targets that were imaged. It can be modelled as follows:
γ = γ t γ n γ o t h e r s ,
where γ t refers to temporal decorrelation, γ n to decorrelation due to signal-to-noise ratio, and γ o t h e r s to other terms that may be present but are negligible in our case [21].
According to [4,22,23], γ t can be expressed, as
γ t = e t / τ , with
τ = 2 σ r 2 λ 4 π 2 ,
where t is the temporal baseline between acquisitions, λ is the radar wavelength, and σ r is the mean motion of scatters in the line of sight (LOS). γ t is a measure of change over time, it approaches 1 in the case as σ r approaches 0.
In the case of management of agricultural grasslands, there are three distinct phases that we need to evaluate in relation to γ : growth, mowing or ploughing event, and developed state. Growth occurs in the early vegetation season and after events. Developed state occurs between the previous states when vegetation is steady and relatively tall.
Movement and growth of vegetation in the LOS causes temporal decorrelation. Coherence is also decreased due to variations in the moisture content of soil and vegetation, which affects scattering centre location, as seen by SAR. During the growth phase, physical and moisture changes are gradual. Changes accumulate as more vegetation is visible in the LOS, causing decorrelation. Coherence would remain low throughout the developed phase.
Removal of vegetation during a mowing event acts in two ways in relation to γ t . The actual event drastically changes the scatterers. γ t calculated for an image pair where the event occurred between the acquisitions would be very low. Depending on the canopy cover, it can be even lower than the coherence of the developed state. For image pairs after the actual event, γ t would be higher when compared to pairs where vegetation was present—as the vegetation is removed, returns from more coherent soil increase.
A similar behaviour applies to ploughing. The event itself causes decorrelation, and conditions after it would cause higher γ t . Total lack of vegetation would suggest even higher γ t values than in the case of mowing.
Precipitation is a source of decorrelation. Variations of relative permittivity in the resolution cell between acquisitions cause coherence to decrease, as scattering centre changes according to permittivity change. The variations are hard to predict as one must consider a multitude of factors, including air temperature, wind, soil and vegetation conditions.
The estimated γ is also a function of Equivalent Number of Looks (ENL) defined by the size of the averaging window in (2) [24]. Contrast between areas of low coherence will be lost if a small averaging window (resulting in low ENL) is used in (2) [25]. We used an averaging window ensuring a square footprint on the ground (approximately 70 × 70 m) and an ENL value between 40 and 50, depending on the swath. According to [24,26], ENL of 45 would ensure contrast for γ values larger than 0.13. Correction for γ n term in (2) was applied according to [27].

3. Results

3.1. Aggregated Time Series

γ v v and γ v h signatures of mowing and ploughing events together with the respective NDVI curves, smoothed with a 3-day averaging window (to reduce the noise), are presented in Figure 3. Understanding the signatures of the farming events are important for adapting satellite monitoring for CAP subsidy checks [28], where the signatures are called ‘markers’.
In the case of mowing, γ decreases slightly but gradually before the event reaching a minimum value just before it. After the event, γ reaches a maximum after which it starts to decrease gradually again. Signatures for both γ v v and γ v h follow each other very closely suggesting invariance to polarisation. Mowing signature is smoother due to more data points (1695 mowing vs 195 ploughing events). Additional variability in ploughing signature is due to other management practices taking place around the time of a ploughing event. Ploughing is often performed closely after mowing and other events—such as seeding or cultivation—might follow. This proximity introduces the large variability evident in Figure 3.
In general, the coherence signature of ploughing is similar to that of mowing. A minimum value is reached just before the event followed by a maximum after it. The maximum, however, occurs much later, between 10 and 20 days after the event.
The NDVI signature of mowing follows the coherence signature inversely. While the mowing event causes coherence to increase from 0.18 to around 0.35, NDVI decreases from 0.75 to 0.5. The NDVI signature of ploughing is less clear. From 0 to 12 days after ploughing NDVI is relatively stable at 0.2, but after that the growth of it has larger variability.
To better understand the dynamics of γ around mowing and ploughing events, coherence values were divided into five time groups. The first group contains measurements where the first acquisition in the interferometric pair was acquired seven to one days before an event and it is labelled as −1 W. The second group labelled as +1 W corresponds to values from up to one week after an event, i.e., acquisition occurred one to seven days after the event, etc. Measurements that were taken on the same day as the event were discarded as the exact time of the activity was not known. Boxplots of these time groups are shown in Figure 4 and the median values are presented in Table 1.
In the case of mowing, the highest median γ v v value is for +1 W group (0.34) with the increase from −1 W being 0.17. γ v h , on the other hand, increases by 0.14, from 0.18 at −1 W to 0.32 at +1 W. On average, γ v v exhibits the largest increase between −1 W and the maximum post-event median. In the following weeks, the difference between γ v v and γ v h continue to be in the range of 0.02 units.
As discussed before, differences between polarisations are larger in the case of ploughing. γ v v increases between −1 W and +1 W by 0.33 and γ v h by 0.10. For both polarisations maximum γ is reached at +3 W being 0.63 for VV and 0.39 for VH polarisation.
Growth of vegetation varies depending on the time of the year. To see whether it influences results of this study we further divided the −1 W and +1 W groups into subgroups depending on the month when the event occurred. Figure 5 and Table 2 contain the results. Results are presented only for months when there were enough events.
For both event types and polarisations, −1 W values stay stable for most of the season, being between 0.16 and 0.18. +1 W values vary more. In the case of mowing, γ v v median values increase from 0.31 in June to 0.35 in July, decrease temporarily to 0.33 in August and reach the maximum at 0.42 in September. γ v h medians follow the same pattern with slightly different absolute values.
In the case of ploughing, the difference between +1 W γ v v in June and August is 0.16. The biggest difference in the case of γ v h is also between June and August—0.10. Higher variability in the case of ploughing is likely due to fewer events. Also, vegetation generally allows to maintain more steady levels of moisture on the parcel, and total lack of it after ploughing would lead to higher sensitivity to moisture changes.

3.2. Relation of Sentinel-1 Coherence to Sentinel-2 NDVI

Figure 6 presents a series of Sentinel-1 backscatter and coherence and Sentinel-2 natural colours and NDVI subsets about a grassland area of Western Estonia. It provides a visual comparison of the traces of mowing events on Copernicus satellite imagery derivatives.
Due to cloud cover and data availability it was not possible to get the Sentinel-1 and -2 images about the exact same dates. Yet the time series gives a good example about the dynamic range of Sentinel-1 coherence and Sentinel-2 NDVI of different mown and unmown grassland parcels. Before mowing the NDVI of tall grass can be above 0.8, whereas after mowing it will fall around 0.4 (Figure 6g,h). Freshly mown areas seem to be well visible also on Sentinel-2 natural colours RGB (Figure 6e,f).
In Sentinel-1 coherence (Figure 6c,d) the pre-event tall grass values are close to the noise floor at 0.15, whereas after mowing it can reach to 0.5 and sometimes even higher. Figure 6 illustrates very well also the value of interferometric processing of Sentinel-1 data. An area mown between 19 and 25 July is indicated with a red oval. Notice that on backscatter image there is almost no difference—grassland backscattering remains low no matter the grass height and the mowing status (Figure 6a,b). While the coherence image is not visually appealing, it contains the information and freshly mown areas are well visible as the bright high coherence areas (Figure 6d).
To explore better the relation between Sentinel-1 coherence and Sentinel-2 NDVI a scatterplot was compiled. Only such X-Y pairs were included where the Sentinel-1 VV-coherence and Sentinel-2 NDVI was available about the very same date of a given grassland parcel. The result is shown in Figure 7.
It is well visible that there is no strong direct correlation (Figure 7), because they measure different physical parameters. Exact Sentinel-1 coherence and Sentinel-2 NDVI values depend on multiple factors. NDVI is mainly sensitive to the green, chlorophyll containing biomass, whereas interferometric coherence is a normalised measure of change between two complex radar images, which in turn in natural conditions are mainly sensitive to water and objects containing water. Yet, the tendency is clear—high coherence corresponds to low NDVI and vice versa. The small cluster in the bottom left (low coherence and low NDVI) corresponds most probably to the parcels, which had low and sparse vegetation, but the coherence is low due to rainfall happening right before either of the InSAR pair images.

3.3. Instructive Individual Time Series

Whereas the general behaviour of Sentinel-1 repeat pass coherence is intuitive to understand—short and sparse vegetation corresponds to high coherence, tall and dense vegetation to low coherence—there are still some details to be considered for correct coherence time series interpretation.
Below, the time series of three instructive parcels are given to illustrate the dynamics of coherence on typical Estonian grasslands. Figure 8 presents an unmown grassland parcel affected by heavy rain. Figure 9 corresponds to a grassland growing on a wet marshy soil and Figure 10 a parcel, which was first mown and then ploughed. The rain, marshland and ploughing effects illustrated below were common in the 2017 dataset used for the study. Understanding them will help to design better monitoring systems for grasslands as well as other agricultural applications exploiting Sentinel-1 coherence time series.
Figure 8 shows the S1 coherence and S2 NDVI time series of a typical unmown grassland. NDVI has a relatively smooth arc like curve being maximum at mid summer and decreasing slowly in the second half of the summer as the senescence of vegetation occurs. The behaviour of coherence is inverse to NDVI—being high in the early and late season and low in the middle.
It is important to notice the effect of rainfall on the time series. The average accumulated past 6 hours precipitation prior to the Sentinel-1 data take on DOY 262 at 15:57 UTC in the nearby Türi, Pajusi and Viljandi weather stations was 22 mm [20]. By DOY 262 due to senescence the vegetation was already drier and sparser resulting in coherence increase as one can see from RON 58 data points. Still, two consecutive RON 160 coherence data points (DOY 256 and 262) are pushed down close to the noise floor leaving the RON 58 time series unaffected. This is a major difference with the effect of farming events as explained in Figure 10 below.
Daily 24 h precipitation sums (accumulated at 18:00 UTC) are also given in Figure 8. Again they are the averages of the closest nearby weather stations from Türi, Pajusi and Viljandi [20], which are not further than 40 km from the particular grassland parcel. One can see that high precipitation values on the same day with the Sentinel-1 acquisition tend to result in low coherence of two consecutive coherence values of a single RON. Precipitation between the InSAR pair acquisitions often do not alter the coherence, as the field can quickly dry out depending on the weather.
Figure 9 presents the S1 VV-coherence and S2 NDVI time series of a typical grassland on a wet peat-soil (e.g., close to peat extraction areas or river valleys). The coherence for the unmown tall grass state is higher-around 0.4 to 0.5, whereas for drier mineral-soils it is below 0.2. NDVI is also slightly lower—never exceeding 0.7, whereas for mineral-soil grasslands it can be over 0.8 in the case of unmown grass. Another property typical for the peat-soil grassland is a general instability of the coherence; even without any farming events, the coherence can fluctuate by 0.2 or more.
The signature of the mowing events on peat-soils is different. Whereas on mineral-soils the mowing event is indicated by the transition of coherence from low to high, it is different on peat-soils. The general coherence levels before and after the mowing event seem virtually unchanged, the only characteristic is a temporary decrease of coherence just before the mowing event to 0.15–0.25 (taking into account that here the coherence at date X presents the coherence between images of date X and date X + 6 days). It is important to notice that all of the RONs are affected and only the interferometric pairs, where the mowing event took place in-between, as it totally changed the vegetation structure.
Figure 10 shows S1 VV-coherence and S2 NDVI time series of a grassland on mineral-soil, which was first mown on DOY 183 and a second time on DOY 222. After the second mowing the parcel was also ploughed twice (on DOYs 227 and 240). The first mowing increased the coherence from below 0.2 close to 0.3, where NDVI decreased from 0.73 to 0.57 at the same time. After the second mowing the coherence increase and the NDVI decrease is almost undetected as it is followed by ploughing just five days after the mowing event. Yet the two ploughing events have very distinct signatures of the coherence time series. The coherence is very low (close to the noise floor) just before the ploughing event and much higher afterwards. It is important to notice how the farming event affects a single data point in all of the RONs—a remarkable difference compared to how rain is affecting coherence, which is decreasing the coherence for two consecutive coherence data points of just a single RON.

4. Discussion

In the following section, possible reasons behind the observed behaviour in the Results part is discussed. In case, available references to the previous studies are given.
Observed coherence behaviour respect to mowing and ploughing events (Figure 3) is in general well in line with the theory in Section 2.6. There is still a detail to be noticed - maximum coherence occurs with a delay. If some vegetation was left on the parcel after mowing it would take some time for it to dry out and stop interacting with the SAR signal. Similar behaviour is reported in a study with Sentinel-2 optical data. When the grass was mown but not removed, the NDVI was still relatively high, it decreased to even lower levels in time when the vegetation was removed or it dried out [29].
The coherence maximum is even more delayed for ploughing events. This might be due to several reasons. Overturned soil with no vegetation cover tends to dry out. Moisture levels change between acquisitions causing decorrelation. After they reach a stable level, coherence can increase. Also, possible follow-up activities such as seeding and fertilisation cause decorrelation, and coherence increases only after all activities cease. Maximum γ value is reached when the soil has dried to a stable moisture level and hence also the scattering patterns between the six-day interferometric pair change less.
Another dissimilarity between mowing and ploughing is the difference due to polarisation. γ v v values are remarkably higher than γ v h after the ploughing event. This is most probably due to the structure of the overturned soil that favours VV backscatter. This should be considered when designing algorithms for classification of management practices. For both event types coherence values reach a minimum just before the actual event. When an event happens between acquisitions, coherence for the affected pair is the lowest due to profound changes in the soil and vegetation caused by mowing or ploughing.
Out of the signatures presented in Figure 3, the NDVI behaviour respect to ploughing event was the least clear. The variability evident 12-30 days after a ploughing event could be caused by different grass species different growth tempo and possible follow-up farming activities such as seeding. Even more variable is the NDVI just before the ploughing event. While for mowing it is very stable high NDVI (as only tall well-developed grassland is usually mown), the parcel status before a ploughing event can be very different: spring-time ploughing of an almost bare field, ploughing after mowing, second ploughing event right after the first one, etc.
VV- and VH-coherence follow generally the same trend, but the behaviour is not exactly one to one. On the level of the individual parcel, the differences can be even larger than on the aggregated plots (Figure 3, Figure 4 and Figure 5) and are related to structural differences of vegetation (primarily vertically oriented stems and leaves or more randomly oriented canopy) [30]. In case of soil scattering, γ is not polarisation dependent considering the moisture component [31]. We must consider that mowing does not completely remove vegetation from a parcel. Unfortunately we do not have information about vegetation height in parcels used in this study. Yet from the previous studies about Estonian grasslands we know that unmown tall grass ranged typically between 50–75 cm and mown short grass between 5–30 cm. Wet biomass ranged from 0.1 (mown grasslands) to 4.5 (unmown tall grass) kg/m 2 [32,33]. It is very likely that similar characteristics were also present during 2017.
Pre-event values of γ v v and γ v h are very stable and are largely determined by ENL used in coherence computation [34], as described in Section 2.6. For the −1 W group, events occurring between the acquisitions cause total decorrelation. The ’noise floor’ due to ENL in our case is 0.13, i.e., coherence cannot be lower than that.
We mentioned in Section 2.2 that there were nine S1 products missing for RON160 and RON87 for a period between 1 May and 30 October. It has negligible influence on the results presented in Figure 3 and Figure 4. However, on a parcel level, lack of S1 acquisitions may be problematic. If a single acquisition is missing, it influences two pairs of acquisitions and a period of 12 days, creating a gap of 18 days between usable coherence values. The severity of this varies depending on the latitude and regional data acquisition plan that influences the data density [35]. It is possible that a single missing acquisition may cause an event to have no visible effect on γ time series. Using multi-RON data can help to solve this challenge as described in the paragraph below.
Whereas a farming event affects a single coherence data point of all the RONs, causing it to be very low for interferometric pairs, where the farming event took place in-between, rainfall caused coherence loss is image and RON specific. It affects just a single RON and leaves other RON coherences unaffected. A significant rainfall just before an SAR image acquisition totally alters the scattering properties of the respective image [36]. As the image is part of two consecutive interferometric pairs (first time when the one is the after image and the second time when the one is the before image)—two consecutive coherences of a certain RON will be temporarily lowered. This underlines the value of using multiple RONs for operational agricultural applications as it helps to solve the ambiguity between rainfall and farming events on coherence time series. The advantage of using multiple RONs is mainly in the shifted temporal sampling and not so much in the different viewing geometries. Rainfall right before a Sentinel-1 data take can mask the traces of a farming event in the time series of a single RON, but (if available) on the other RONs the traces are still visible.
Figure 8 presented also the 24 h accumulated precipitation data, allowing to compare the rainfall info with coherence time series. Care must still be taken when interpreting it. During mid-season (DOY 150–230) due to tall well-developed vegetation the coherence was already very low and the precipitation effects are therefore not so pronounced for this period. The resolution of the given precipitation data from the Estonian Weather Service is 24 h, most of the rainfall could have been many hours before the Sentinel-1 data take (around 15:00 UTC) and by that time the fields could have already been dried out. Studying precipitation effects on Sentinel-1 coherence was not the main focus of our research. For the readers interested in this topic it is recommended to select short grass low biomass grassland parcels as close as possible to the weather stations. There, the dynamics of coherence due to precipitation should come out in the most pure and clear way.
The fluctuating coherence levels of the grassland parcel on peat-soil are likely connected with the moisture changes in the soil and vegetation. This is altering the scattering properties and as a result also the coherence. High Sentinel-1 coherence on peatland vegetation is also reported in [37], which discusses in detail the relation between the coherence and water regime.

5. Conclusions

According to the study, mowing and ploughing can be detected from Sentinel-1 coherence time series with high confidence. We demonstrate it based on a large ground reference data set covering 1695 mowing and 195 ploughing events on agricultural grassland in Estonia. While the vegetation is present, the values are close to the noise floor. For both event types, coherence tends to increase after an event as scattering from the soil becomes dominant. As the vegetation regrows, coherence slowly decreases to pre-event values. On average, ploughing causes a larger increase, at least for VV polarisation. In the case of mowing, both polarisations exhibit very similar behaviour to each other. The magnitude of increase varies depending not only on the event type but also on vegetation species and specific management practices within one event type.
We discussed the effect of rainfall and farming events on Sentinel-1 coherence time series and explained the means to distinguish them using the time series of alternative viewing geometries and independent temporal sampling (different RONs). We also compared how the coherence time series for grassland on mineral soil differ from the one on peat-soil.
Precipitation right before a Sentinel-1 image acquisition can suppress coherence, masking the increase due to event and hiding the signature of mowing or ploughing. For reliable operational applications, a dense multi-RON time series should be used for the analysis to reduce the effect of precipitation caused disturbances.
The results presented here show the traces of mowing and ploughing events on Sentinel-1 coherence time series convincingly. However, the detection limits are not fully understood. Future studies should compare the coherence values with in situ grass height and wet biomass measurements. This could help to find how much remainder biomass could be left on ground to leave the farming event still visible on Sentinel-1 coherence time series and possibly contribute towards a grassland biomass retrieval algorithm based on Sentinel-1 and other data sources.

Author Contributions

Conceptualization, K.V. and K.Z.; methodology, K.Z. and K.V.; software, I.S., T.V. and T.T.; validation, J.P., A.B. and D.J.; investigation, K.Z. and K.V.; resources, K.V. and K.K.; data curation, I.S. and T.T.; writing—original draft preparation, K.Z. and K.V.; writing—review and editing, K.V. and J.P.; visualization, I.S., K.V. and K.Z.; supervision, J.P.; project administration, K.V. and K.K.; funding acquisition, K.V. and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded through the following projects: “Data Analytics for Optimizing Agricultural Monitoring” funded in part by the Software Technology and Applications Competence Centre, “Assessment of Grassland Quality and Quantity Parameters and Management Activities Using Sentinel-1&2 data (SentiGrass)” financed by the Plan for European Cooperating States programme of the European Space Agency (ESA), and European Regional Development Fund within National Programme for Addressing Socio-Economic Challenges through R&D (RITA) and under the grant No EU48684.

Acknowledgments

We acknowledge the support from the Estonian Agricultural Registers and Information Board as well as farmers for providing the data. Special thanks to Dr Aire Olesk for careful language review of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Eurostat. Farm Structure Survey 2013—Main Results. Available online: https://ec.europa.eu/eurostat/statistics-explained/index.php?title=Archive:Farm_structure_survey_2013_-_main_results (accessed on 14 December 2015).
  2. Hopkins, A.; Holz, B. Grassland for agriculture and nature conservation: Production, quality and multi-functionality. Agron. Res. 2006, 4, 3–20. [Google Scholar]
  3. d’Andrimont, R.; Lemoine, G.; van der Velde, M. Targeted Grassland Monitoring at Parcel Level Using Sentinels, Street-Level Images and Field Observations. Remote Sens. 2018, 10, 1300. [Google Scholar] [CrossRef] [Green Version]
  4. Zebker, H.A.; Villasenor, J. Decorrelation in interferometric radar echoes. IEEE Trans. Geosci. Remote Sens. 1992, 30, 950–959. [Google Scholar] [CrossRef] [Green Version]
  5. Hagberg, J.O.; Ulander, L.M.H.; Askne, J. Repeat-pass SAR interferometry over forested terrain. IEEE Trans. Geosci. Remote Sens. 1995, 33, 331–340. [Google Scholar] [CrossRef]
  6. Wegmuller, U.; Werner, C. Retrieval of vegetation parameters with SAR interferometry. IEEE Trans. Geosci. Remote Sens. 1997, 35, 18–24. [Google Scholar] [CrossRef]
  7. Askne, J.; Santoro, M.; Smith, G.; Fransson, J.E.S. Multitemporal repeat-pass SAR interferometry of boreal forests. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1540–1550. [Google Scholar] [CrossRef]
  8. Santoro, M.; Wegmuller, U.; Askne, J.I.H. Signatures of ERS–Envisat interferometric SAR coherence and phase of short vegetation: An analysis in the case of maize fields. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1702–1713. [Google Scholar] [CrossRef]
  9. Lavalle, M.; Simard, M.; Hensley, S. A temporal decorrelation model for polarimetric radar interferometers. IEEE Trans. Geosci. Remote Sens. 2012, 50, 2880–2888. [Google Scholar] [CrossRef]
  10. Zalite, K.; Antropov, O.; Praks, J.; Voormansik, K.; Noorma, M. Monitoring of agricultural grasslands with time series of X-band repeat-pass interferometric SAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 3687–3697. [Google Scholar] [CrossRef]
  11. Ali, I.; Barrett, B.; Cawkwell, F.; Green, S.; Dwyer, E.; Neumann, M. Application of Repeat-Pass TerraSAR-X Staring Spotlight Interferometric Coherence to Monitor Pasture Biophysical Parameters: Limitations and Sensitivity Analysis. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 3225–3231. [Google Scholar] [CrossRef] [Green Version]
  12. Tamm, T.; Zalite, K.; Voormansik, K.; Talgre, L. Relating Sentinel-1 interferometric coherence to mowing events on grasslands. Remote Sens. 2016, 8, 802. [Google Scholar] [CrossRef] [Green Version]
  13. Taravat, A.; Wagner, M.; Oppelt, N. Automatic Grassland Cutting Status Detection in the Context of Spatiotemporal Sentinel-1 Imagery Analysis and Artificial Neural Networks. Remote Sens. 2019, 11, 711. [Google Scholar] [CrossRef] [Green Version]
  14. Garioud, A.; Valero, S.; Giordano, S.; Mallet, C. On the Joint Exploitation of Optical and SAR Imagery for Grassland Monitoring. Available online: https://hal.archives-ouvertes.fr/hal-02571832/document (accessed on 21 August 2020).
  15. Kavats, O.; Khramov, D.; Sergieieva, K.; Vasyliev, V. Monitoring Harvesting by Time Series of Sentinel-1 SAR Data. Remote Sens. 2019, 11, 2496. [Google Scholar] [CrossRef] [Green Version]
  16. Shang, J.; Liu, J.; Poncos, V.; Geng, X.; Qian, B.; Chen, Q.; Dong, T.; Macdonald, D.; Martin, T.; Kovacs, J.; et al. Detection of Crop Seeding and Harvest through Analysis of Time-Series Sentinel-1 Interferometric SAR Data. Remote Sens. 2020, 12, 1551. [Google Scholar] [CrossRef]
  17. Jacob, A.; Vicente-Guijalba, F.; Lopez-Martinez, C.; Lopez-Sanchez, J.; Litzinger, M.; Kristen, H.; Mestre-Quereda, A.; Ziółkowski, D.; Lavalle, M.; Notarnicola, C.; et al. Sentinel-1 InSAR Coherence for Land Cover Mapping: A Comparison of Multiple Feature-Based Classifiers. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 535–552. [Google Scholar] [CrossRef] [Green Version]
  18. Monti-Guarnieri, A.; Manzoni, M.; Giudici, D.; Recchia, A.; Tebaldini, S. Vegetated Target Decorrelation in SAR and Interferometry: Models, Simulation, and Performance Evaluation. Remote Sens. 2020, 16, 2545. [Google Scholar] [CrossRef]
  19. Rouse, J.W.; Haas, R.H.; Scheel, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS. In Proceedings of the 3rd Earth Resource Technology Satellite (ERTS) Symposium, Washington, DC, USA, 10–14 December 1973; Volume 1, pp. 48–62. [Google Scholar]
  20. EWS, E.W.S. Estonian Weather Observations Archive. Available online: http://www.ilmateenistus.ee/ilm/ilmavaatlused/vaatlusandmed/?lang=en (accessed on 10 November 2020).
  21. Hanssen, R.F. Radar Interferometry: Data Interpretation and Error Analysis; Springer Science & Business Media: Berlin, Germany, 2001; Volume 2. [Google Scholar]
  22. Rocca, F. Modeling interferogram stacks. IEEE Trans. Geosci. Remote Sens. 2007, 45, 3289–3299. [Google Scholar] [CrossRef]
  23. Morishita, Y.; Hanssen, R.F. Temporal decorrelation in L-, C-, and X-band satellite radar interferometry for pasture on drained peat soils. IEEE Trans. Geosci. Remote Sens. 2015, 53, 1096–1104. [Google Scholar] [CrossRef]
  24. Lee, J.; Hoppel, K.W.; Mango, S.A.; Miller, A.R. Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery. IEEE Trans. Geosci. Remote Sens. 1994, 32, 1017–1028. [Google Scholar]
  25. Touzi, R.; Lopes, A.; Bruniquel, J.; Vachon, P.W. Coherence estimation for SAR imagery. IEEE Trans. Geosci. Remote Sens. 1999, 37, 135–149. [Google Scholar] [CrossRef] [Green Version]
  26. Cumming, I.G.; Wong, F.H. Digital processing of synthetic aperture radar data. Artech House 2005, 1, 3. [Google Scholar]
  27. Just, D.; Bamler, R. Phase statistics of interferograms with applications to synthetic aperture radar. Appl. Opt. 1994, 33, 4361–4368. [Google Scholar] [CrossRef] [PubMed]
  28. Devos, W.; Lemoine, G.; Milenov, P.; Fasbender, D.; Loudjani, P.; Wirnhardt, C.; Sima, A.; Griffiths, P. Second Discussion Document on the Introduction of Monitoring To Substitute OTSC: Rules for Processing Applications in 2018–2019; JRC Science Hub: Ispra, Italy, 2018. [Google Scholar]
  29. Kolecka, N.; Ginzler, C.; Pazur, R.; Price, B.; Verburg, P.H. Regional scale mapping of grassland mowing frequency with Sentinel-2 time series. Remote Sens. 2018, 10, 1221. [Google Scholar] [CrossRef] [Green Version]
  30. Treuhaft, R.N.; Siqueira, P.R. Vertical structure of vegetated land surfaces from interferometric and polarimetric radar. Radio Sci. 2000, 35, 141–177. [Google Scholar] [CrossRef] [Green Version]
  31. Zwieback, S.; Hensley, S.; Hajnsek, I. Assessment of soil moisture effects on L-band radar interferometry. Remote Sens. Environ. 2015, 164, 77–89. [Google Scholar] [CrossRef] [Green Version]
  32. Voormansik, K.; Jagdhuber, T.; Olesk, A.; Hajnsek, I.; Papathanassiou, K. Towards a detection of grassland cutting practices with dual polarimetric TerraSAR-X data. Int. J. Remote Sens. 2013, 34, 8081–8103. [Google Scholar] [CrossRef]
  33. Voormansik, K.; Jagdhuber, T.; Zalite, K.; Noorma, M.; Hajnsek, I. Observations of cutting practices in agricultural grasslands using polarimetric SAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 9, 1382–1396. [Google Scholar] [CrossRef]
  34. Lee, J.-S.; Pottier, E. Polarimetric Radar Imaging—From Basics to Applications; CRC Press: Boca Raton, FL, USA, 2010; pp. 120–121. [Google Scholar]
  35. Torres, R.; Snoeij, P.; Geudtner, D.; Bibby, D.; Davidson, M.; Attema, E.; Potin, P.; Rommen, B.; Floury, N.; Brown, M.; et al. GMES Sentinel-1 mission. Remote Sens. Environ. 2012, 120, 9–24. [Google Scholar] [CrossRef]
  36. Ulaby, F.T.; Brisco, B.; Dobson, C. Improved spatial mapping of rainfall events with spaceborne SAR imagery. IEEE Trans. Geosci. Remote Sens. 1983, 1, 118–121. [Google Scholar] [CrossRef]
  37. Tampuu, T.; Praks, J.; Uiboupin, R.; Kull, A. Long Term Interferometric Temporal Coherence and DInSAR Phase in Northern Peatlands. Remote Sens. 2020, 12, 1566. [Google Scholar] [CrossRef]
Figure 1. Geographic distribution of events used in this study. Size of markers not to scale.
Figure 1. Geographic distribution of events used in this study. Size of markers not to scale.
Remotesensing 12 03784 g001
Figure 2. Typical grassland in Estonia before (left column, (a,c)) and after mowing (right column, (b,d)) with dominant Poaceae (top row, (a,b)) and dominant Fabaceae (bottom row, (c,d)) vegetation.
Figure 2. Typical grassland in Estonia before (left column, (a,c)) and after mowing (right column, (b,d)) with dominant Poaceae (top row, (a,b)) and dominant Fabaceae (bottom row, (c,d)) vegetation.
Remotesensing 12 03784 g002
Figure 3. Median values of S1 coherence (top) and S2 NDVI (bottom) respect to mowing (left) and ploughing (right) events. Coloured regions represent the interquartile range. Events are from 2017, spanning between May and October (mowing), and June and September (ploughing).
Figure 3. Median values of S1 coherence (top) and S2 NDVI (bottom) respect to mowing (left) and ploughing (right) events. Coloured regions represent the interquartile range. Events are from 2017, spanning between May and October (mowing), and June and September (ploughing).
Remotesensing 12 03784 g003
Figure 4. Evolution of γ over time around mowing and ploughing events represented by boxplots. Top row contains boxplots for mowing while the bottom one for ploughing. Left and right columns correspond to VV and VH polarisation, respectively.
Figure 4. Evolution of γ over time around mowing and ploughing events represented by boxplots. Top row contains boxplots for mowing while the bottom one for ploughing. Left and right columns correspond to VV and VH polarisation, respectively.
Remotesensing 12 03784 g004
Figure 5. Boxplots of −1 W and +1 W subgroups defined by months when events were recorded. Top row contains boxplots for mowing while the bottom one for ploughing. Left and right columns correspond to VV and VH polarisation, respectively.
Figure 5. Boxplots of −1 W and +1 W subgroups defined by months when events were recorded. Top row contains boxplots for mowing while the bottom one for ploughing. Left and right columns correspond to VV and VH polarisation, respectively.
Remotesensing 12 03784 g005
Figure 6. A subset of Western Estonian grassland area in 2017 on Sentinel-1 backscatter image (Red = VV, Green = VH and Blue = VV + VH) from 19 July (a), 25 July (b), coherence between 19 and 25 July (c), 25 and 31 July (d). Red oval denotes a freshly mown area. The same area on Sentinel-2 natural colours and NDVI image on 14 July (e and g respectively) and on 3 August (f and h).
Figure 6. A subset of Western Estonian grassland area in 2017 on Sentinel-1 backscatter image (Red = VV, Green = VH and Blue = VV + VH) from 19 July (a), 25 July (b), coherence between 19 and 25 July (c), 25 and 31 July (d). Red oval denotes a freshly mown area. The same area on Sentinel-2 natural colours and NDVI image on 14 July (e and g respectively) and on 3 August (f and h).
Remotesensing 12 03784 g006
Figure 7. Scatterplot between Sentinel-1 VV-coherence (X-axis) and Sentinel-2 NDVI (Y-axis), selected are only those parcel averages, where a Sentinel-1 and -2 image was available about the same day.
Figure 7. Scatterplot between Sentinel-1 VV-coherence (X-axis) and Sentinel-2 NDVI (Y-axis), selected are only those parcel averages, where a Sentinel-1 and -2 image was available about the same day.
Remotesensing 12 03784 g007
Figure 8. S1 VV-coherence (relative orbit number (RON) 160—solid black line, RON 58—solid grey line) and S2 NDVI (grey dotted line) time series of an unmown grassland parcel. Daily 24 h precipitation sums (accumulated at 18:00 UTC) are shown on the secondary y-axis with orange. Dashed oval notes the coherence data points affected by the heavy rainfall on Day of Year (DOY) 262.
Figure 8. S1 VV-coherence (relative orbit number (RON) 160—solid black line, RON 58—solid grey line) and S2 NDVI (grey dotted line) time series of an unmown grassland parcel. Daily 24 h precipitation sums (accumulated at 18:00 UTC) are shown on the secondary y-axis with orange. Dashed oval notes the coherence data points affected by the heavy rainfall on Day of Year (DOY) 262.
Remotesensing 12 03784 g008
Figure 9. S1 VV-coherence (RON 160—solid black line, RON 7—dashed black line, RON 87—dashed grey line) and S2 NDVI (grey dotted line) time series of a grassland growing on marshy soil. Thin vertical lines indicate the mowing events.
Figure 9. S1 VV-coherence (RON 160—solid black line, RON 7—dashed black line, RON 87—dashed grey line) and S2 NDVI (grey dotted line) time series of a grassland growing on marshy soil. Thin vertical lines indicate the mowing events.
Remotesensing 12 03784 g009
Figure 10. S1 VV-coherence (RON 160—solid black line, RON 7—dashed black line, RON 58—solid grey line) and S2 NDVI (grey dotted line) time series of a grassland parcel, which was mown and ploughed twice. Thin vertical lines indicate the mowing and ploughing events.
Figure 10. S1 VV-coherence (RON 160—solid black line, RON 7—dashed black line, RON 58—solid grey line) and S2 NDVI (grey dotted line) time series of a grassland parcel, which was mown and ploughed twice. Thin vertical lines indicate the mowing and ploughing events.
Remotesensing 12 03784 g010
Table 1. Median values of time groups in Figure 4 for combinations of event type and polarisation.
Table 1. Median values of time groups in Figure 4 for combinations of event type and polarisation.
−1 W+1 W+2 W+3 W+4 W
M γ v v 0.170.340.300.260.25
γ v h 0.180.320.300.270.26
P γ v v 0.160.490.570.630.57
γ v h 0.190.290.360.390.36
Table 2. Median values of −1 W and +1 W subgroups defined by month when events were recorded. ‘−’ and ‘+’ correspond to −1 W and +1 W, respectively. M corresponds to Mowing and P to Ploughing.
Table 2. Median values of −1 W and +1 W subgroups defined by month when events were recorded. ‘−’ and ‘+’ correspond to −1 W and +1 W, respectively. M corresponds to Mowing and P to Ploughing.
JunJulAugSep
++++
M γ v v 0.160.310.170.350.170.330.180.42
γ v h 0.170.290.180.340.180.310.180.37
P γ v v 0.170.370.160.380.160.53
γ v h 0.180.230.180.240.170.33
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Voormansik, K.; Zalite, K.; Sünter, I.; Tamm, T.; Koppel, K.; Verro, T.; Brauns, A.; Jakovels, D.; Praks, J. Separability of Mowing and Ploughing Events on Short Temporal Baseline Sentinel-1 Coherence Time Series. Remote Sens. 2020, 12, 3784. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223784

AMA Style

Voormansik K, Zalite K, Sünter I, Tamm T, Koppel K, Verro T, Brauns A, Jakovels D, Praks J. Separability of Mowing and Ploughing Events on Short Temporal Baseline Sentinel-1 Coherence Time Series. Remote Sensing. 2020; 12(22):3784. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223784

Chicago/Turabian Style

Voormansik, Kaupo, Karlis Zalite, Indrek Sünter, Tanel Tamm, Kalev Koppel, Tarvi Verro, Agris Brauns, Dainis Jakovels, and Jaan Praks. 2020. "Separability of Mowing and Ploughing Events on Short Temporal Baseline Sentinel-1 Coherence Time Series" Remote Sensing 12, no. 22: 3784. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12223784

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop