
Sentinel-1 Interferometric Coherence Image Services
Enhancing access to the Global Seasonal Sentinel-1 Interferometric Coherence and Backscatter Dataset with Image Services
Global Seasonal Sentinel-1 Interferometric Coherence and Backscatter Dataset (GSSICB)
This globally consistent set of Sentinel-1 coherence data was developed by Earth Big Data LLC and Gamma Remote Sensing AG , under contract for NASA's Jet Propulsion Laboratory , and published in 2021 by Kellndorfer et al. This dataset provides median 6-, 12-, 18-, 24-, 36- and 48-day seasonal coherence values and mean seasonal backscatter for Sentinel-1 acquisitions from December 1, 2019 through November 30, 2020.
Sentinel-1 is a C-band Synthetic Aperture Radar (SAR) mission developed by ESA. This dataset contains modified Copernicus Sentinel data 2019-2020, processed by ESA. The GSSICB effort used Single Look Complex (SLC) data collected in Interferometric Wide-Swath Mode (IW) mode to generate the data products.
For more information about the dataset, refer to the Global seasonal Sentinel-1 interferometric coherence and backscatter data set article published in Scientific Data .
Coherence and Backscatter
Phase
The phase indicates the location of the signal along its wave cycle when it returns to the sensor, and can be used to determine the relative distance between the sensor to the target. A common use for phase measurements is SAR Interferometry (InSAR), where phase measurements from two different acquisitions are differenced to identify and quantify deformation of the earth's surface.
The phase measurements in the two images used in InSAR must be coherent in order to detect change. Random changes in phase from one acquisition to the next can mask actual surface deformation. Vegetation is a common driver of decorrelation, as changes can easily take place in the interval between two acquisitions due to growth, seasonal changes, or wind effects.
Coherence
Coherence is a measure of the correlation between the phase measurements of two acquisitions, and is an indicator of the quality of the InSAR product that would be generated by those two acquisitions. Coherence values can also be used as an indicator of change on the landscape. We will explore this more in the Coherence section.
Amplitude
The amplitude measurements recorded by the SAR sensor indicate how much of the signal that was sent out returns to sensor to be measured. It is also referred to as backscatter or intensity. The amplitude is generally used to generate SAR imagery, and we will explore this more in the Backscatter section.
Accessing the Data
The GSSICB data is publicly available as tiled Cloud-Optimized GeoTIFFs in the Registry of Open Data on AWS , and archived in NASA’s Earthdata Cloud by the Alaska Satellite Facility Distributed Active Archive Center (ASF DAAC). As part of ASF’s mission to make SAR data more accessible, we have published image services for some of the products in this dataset.
The image services are available as REST endpoints in the GSSICB directory on NASA's Earthdata GIS Image Server, and are enabled as WMS services . They can be easily added to new or existing web maps, accessed from GIS software programs, or included in other applications that support ArcGIS Image Services or WMS layers.
A separate service is published for each season and polarization available for a data type. The services reference the source rasters hosted in the Registry of Open Data on AWS, and currently include the following:
- Median 12-day Coherence (COH12)
- both available polarizations: VV and HH
- four seasons: winter (DJF), spring (MAM), summer (JJA), fall (SON)
- Median 6-day Coherence (COH06)
- both available polarizations: VV and HH
- four seasons: winter (DJF), spring (MAM), summer (JJA), fall (SON)
- Mean Backscatter (AMP)
- all four available polarizations: VV, VH, HV, HH
- four seasons: winter (DJF), spring (MAM), summer (JJA), fall (SON)
Overview Map
The web map below displays all of the image services currently published for this dataset. For more targeted maps and tutorials, refer to the other sections of this StoryMap.
The services are grouped by data type and polarization. Within each grouping are the services for each of the four seasons (nomenclature is based on the northern hemisphere):
- Winter (DJF): December, January, February
- Spring (MAM): March, April, May
- Summer (JJA): June, July, August
- Fall (SON): September, October, November
The December data is from 2019, and data for all other months is from 2020.
Using the Web Map
Click through the slideshow for an orientation to the web map interface.
Note that some of the layer and group names in the GSSICB Web Map are different from the maps shown in this orientation, but the layer content and functionality is unchanged.
Explore the Web Map
Explore the web map embedded below, or click the Open live content in a new tab floating icon to launch it in its own browser window.
Web map displaying all of the image services for the Global Seasonal Sentinel-1 Interferometric Coherence and Backscatter Data Set (GSSICB), published by the Alaska Satellite Facility (ASF).
Polarizations
The Sentinel-1 acquisition plan generally collects data with a vertical primary polarization (VV and VH) over land, but with a horizontal primary polarization (HH and HV) over the polar regions.
The interferometric coherence was calculated for the co-polarized data only (VV and HH), and these two different polarizations are published in separate image services.
The backscatter datasets are generated for all available polarizations, with VV and VH covering most global landmasses and HH and HV covering the polar regions. Again, each polarization is published as a separate image service.
Use the swipe tool below to compare the spatial extent of the COH12 VV service on the left and HH service on the right. Both display the winter (December, January, February) data. This same geospatial separation based on primary polarization is present for all of the datasets.
Compare the spatial extent of the COH12 VV image service on the left and HH on the right.
Coherence
Interferometric SAR (InSAR) compares the phase measurements of two SAR acquisitions to look for differences, which indicate that the ground has moved relative to the sensor in the time between the two acquisitions. The spectra of two images used for InSAR must overlap well enough to generate interferometric fringes, which are used to identify and quantify surface deformation. Correlation indicates the comparability of the phase information between the two acquisitions. The magnitude of the correlation is commonly referred to as coherence.
Coherence values are also useful for change detection, or to determine the variability of surface characteristics in an area. Low coherence values serve as an indicator of change, and high coherence values identify areas with stable or persistent scatterers. When coherence values are aggregated over time, we can see patterns of variability across the landscape.
Each satellite in the Sentinel-1 constellation has a 12-day repeat cycle, so most of the globe has acquisitions over the same place on earth from the same place in space every 12 days. The Sentinel-1 mission is designed to be a two-satellite constellation, with each satellite having the same orbit, but offset by 180 degrees. If both satellites acquired imagery all the time, we would expect an acquisition over the same place on earth from the same place in space every six days. In reality, the mission plan is such that only Europe and select other regions can expect 6-day repeat imagery, while the rest of the world generally has a 12-day repeat cycle when both satellites are operational.
The GSSICB Dataset includes coherence calculations for a range of intervals, including 6, 12, 18, 24, 36 and 48 days. The first image services ASF published were for median 12-day coherence, but we are working on publishing services for the remaining intervals, and they will be added to the overview map as they become available.
The Sentinel-1A satellite was launched in 2014 and is still actively acquiring data. The Sentinel-1B satellite was launched in 2016, and stopped acquiring data at the end of December 2021. Both satellites were active during the time period used for the calculation of the GSSICB datasets. Sentinel-1C is ready for launch, and will replace Sentinel-1B once a launch vehicle is secured.
COH Data Format
Example of pop-up displaying both unscaled and scaled coherence values. Source raster values are scaled by multiplying the coherence by 100, but by default we divide those values by 100 to display unscaled coherence values between 0 and 1.
Coherence values typically range from 0 to 1. The greater the value, the higher the coherence (and the better the correlation). For the GSSICB dataset, the coherence source rasters have been converted to an integer for more efficient storage, so values range from 0 to 100.
ASF publishes raster function templates for the image services that include both a direct representation of the source raster values, which have been multiplied by 100 to give a range of values between 0 and 100, and a conversion to the standard coherence values, ranging from 0 to 1. The latter is the default template, but if you click on a pixel in the web map, both values will be displayed in the pop-up.
Dark pixels in the coherence image services indicate poor correlation, while the bright pixels indicate good correlation. Scroll through the text next to the map below and click the buttons to explore coherence values.
COH12 Data
12-day Median Interferometric Coherence
The COH12 data was generated by calculating coherence values for pairs of Sentinel-1 images, each acquired 12 days apart, then calculating the median coherence values for each season from these individual coherence rasters.
Because most of the world had 12-day repeat coverage during the time period represented in the GSSICB datasets (December 2019 - November 2020), this is the coherence dataset with the most complete coverage. We will use the COH12 dataset to explore the concept of coherence
Exploring COH12 VV
This map displays the COH12 VV data for the summer season (June/July/August) in France.
Follow the prompts to explore the map. You can also pan and zoom to take a closer look if you see something of interest. Click on the map to view pixel values and more information about the source rasters.
Urban Areas
Note that the cities show up clearly, dotted around the landscape. Urban areas typically have high coherence, as the built structures comprising cities tend to be stable, effective scatterers.
Paris is a very large bright spot on the map.
Agricultural Regions
The agricultural regions outside the city display a field-by-field patchwork of coherence values, but generally have much lower coherence than the urban areas.
During the growing season, agricultural lands undergo constant and significant change that results in decorrelation of the phase from one acquisition to the next. Changes in crop growth stages and activities around tilling and harvesting can all contribute to this decorrelation.
Seasonal Differences
This urban-rural difference is much less pronounced during the winter months in areas with strong seasonality.
The agricultural fields are still patchy, but have much higher coherence in the winter than they did in the summer. The settlements amidst the fields are much less obvious in the winter season.
Toggle between the winter (December/January/February) and summer (June/July/August) seasons to compare the coherence values across the country.
The coherence values of the rural areas are not as different from the urban areas during the winter, when many of the fields are either planted with slow-growing winter grain or cover crops, or left fallow.
Forested Areas
There is a persistent dark spot in the middle of the country. This dark area becomes even more obvious during the winter season.
This patch with consistently low coherence values is located just south of Orleans.
The Loire Valley is densely vegetated, with extensive forests, orchards and vineyards.
Areas with dense and complex vegetation typically have low coherence. This can be problematic if you're interested in generating InSAR products in densely vegetated regions, particularly when using SAR sensors with relatively short wavelengths.
Notice how low the Sentinel-1 coherence values are in the Amazon rainforest. Click on the map to view coherence values for the darkest areas in northern Brazil.
Sentinel-1 has a wavelength of about 5.6 cm, which allows some penetration of dense forest canopy, but the signal will rarely make it through to the forest floor. Forest canopy tends to have poor correlation due to changes in vegetative structure and wind-driven movement over time. This makes it very difficult to use Seninel-1 for InSAR in these areas.
SAR datasets collected with L-band sensors (~25 cm wavelength), such as the upcoming NISAR mission , may provide better results in densely forested regions.
In South America, the coherence values are even lower in the December/January/February season (southern hemisphere summer) than in June/July/August. The Amazon, however, has consistently low coherence values in all seasons.
Mountains
Mountain ranges are not particularly easy to pick out based on summer coherence values.
Some ranges may have higher coherence than surrounding areas, especially if they are rocky or sparsely vegetated. Others may exhibit similar or lower coherence compared to surrounding areas.
It's a different story in the winter, however.
Mountain ranges, especially those with very high elevation, generally have more variable snow and ice conditions from one image to the next than the surrounding areas during the winter. This results in low coherence values, clearly illustrated by very stark contrasts in the winter coherence maps.
COH12 Web Map
Explore all of the seasons of the COH12 VV dataset using the web map below. Try changing the basemap to compare the coherence data to optical imagery or terrain layers.
Map displaying all four seasons of the COH12 VV dataset.
Backscatter
The backscatter datasets, sometimes referred to as intensity or amplitude (AMP), have undergone Radiometric Terrain Correction (RTC) to produce "terrain flattened" gamma-nought (γ 0 ) backscatter images. The mean backscatter from these images was calculated for each available polarization for each 3-month season.
RTC imagery is a common way to visualize the earth's surface with SAR, and different polarizations (VV, VH, HV, HH) reveal different characteristics. The GSSICB backscatter datasets provide a way to view these characteristics through the seasons, and to examine the relationship between the AMP and COH datasets.
- VV polarization is sensitive to surface roughness
- rougher surfaces will have higher backscatter (brighter in the image), smooth surfaces tend to have low backscatter values
- VH and HV polarizations are sensitive to volume scatterers (i.e. vegetation canopies)
- Dense, complex vegetation, such as a deciduous forest, has higher cross-polarized backscatter than sparsely vegetated areas or vegetation with simple structure, such as alpine areas or grasslands
- HH polarization is sensitive to double-bounce scattering, such as what you see from buildings, tree trunks, or other abrupt angles
AMP Data Format
The AMP source rasters for the GSSICB dataset use digital numbers (DN) for the pixel values. Most RTC imagery is scaled to power, amplitude (square root of the power values), or dB (10 times the log10 of the power), and if you're accustomed to those scaled values, the DN values will not be particularly meaningful.
We include conversion raster function templates in the image services that provide the conversion to power and dB scale from the source DN values. When you click on a pixel in the map, the pop-up will display all three values, along with the conversion formulas used to calculate the different scales:
- Power values are calculated using the formula (DN) 2 / 199526231
- dB values are calculated using the formula 20 x log10(DN) - 83
Exploring AMP VV/VH
This map displays the AMP data for VV and VH polarizations.
Follow the prompts to explore the map. You can also pan and zoom to take a closer look if you see something of interest. Click on the map to view pixel values and more information about the source rasters.
Surface Water
Calm surface water generally has very low backscatter. Water effectively reflects the SAR signal, but because the signal arrives at an oblique angle to the surface, the signal bounces off in the other direction. Very little of that sent signal makes it back to the sensor to be measured.
Lac Saint-Jean is a fairly large lake in Québec. It covers more than 400 square miles (1,000 km 2 ), and is ice-free for the months included in the summer season (June/July/August) of the GSSICB dataset. The open water of this lake and other smaller lakes, ponds, streams and rivers in the area, stands out clearly as very dark features on the landscape.
Wind Impacts
On windy days, the surface of the water becomes more rough, and VV backscatter is much higher. In some conditions, it can be difficult to tell the difference between lakes and the surrounding land when looking at a single VV image. Because the GSSICB dataset averages together three months of data, however, even sizeable lakes such as Lac Saint-Jean appear to have uniformly low backscatter.
Very large lakes, such as the Great Lakes in North America, are more likely to have wave action throughout the open water season, so surface conditions are much less uniform than on smaller lakes. Even averaging the seasonal VV backscatter values together will not completely smooth out the signal caused by variable wind conditions, as seen here for Lake Ontario.
The VH polarization is much less sensitive to surface roughness, so does not exhibit the wave-driven variability seen in the VV data. The VH backscatter returns are generally much lower than VV returns overall, however, which makes them more susceptible to noise, such as the burst-level artifacts visible here. Even with the noise, the VH values over the water are still generally much lower than what is seen on land in this area, making it an effective way to identify water pixels.
Click on some of the pixels to view the dB values over the lake in the VH dataset, and compare them to the VV values in the previous view. The VH values (even the brightest ones) are all less than -27 dB, while the VV values range from about -15 to -20 dB.
Ice
Great Bear Lake is located in northern Canada, and is generally free of ice from mid-July through September.
From a distance, the summer VV values (which includes data from June, July and August) appear to be uniformly low over the lake.
If you zoom in further, however, and the dynamic range adjustment in the web map stretches the values of the dark pixels along the grayscale, you can see the impact of including acquisitions from when the lake was not yet ice-free.
Zoom out one level and pan around the lake to see the different patterns imprinted on the June/July/August dataset by the inclusion of acquisitions with ice cover.
Again, even though the lake exhibits some variability in the backscatter values, the open water pixels are still much darker than those of the surrounding uplands.
In winter, Great Bear Lake freezes over. The ice is much more effective at scattering the signal back towards the sensor, especially when the ice surface is rough. A lake this large takes some time to freeze over, and the impacts of wind, currents, and human activities can leave their mark. Even though this dataset displays mean values, you can see ice features that are persistent enough through the season to be visible here.
Use the zoom buttons and pan around the lake to explore the ice.
There is one community on Great Bear Lake: Délı̨nę.
In the winter, this community is linked to others along the Mackenzie River by winter roads, including an ice road across the southwestern tip of Great Bear Lake.
While there is an official ice road, other corridors can also be established by activity such as snowmobile travel.
Compare the AMP VV datasets from winter (December/January/February) to spring (March/April/May).
Some of the travel routes across the ice remain the same, while some new routes are established and some old routes are abandoned or adjusted.
Zoom out one level and pan around the lake to see common travel corridors.
Urban-Rural Contrasts
Let's return to Paris and take a quick look at the optical imagery northeast of the city.
Notice the urban, agricultural and forested areas in this region. Now add the VV backscatter imagery.
As with coherence values, amplitude values tend to be relatively high over urban areas, particularly in the VV polarization. Agricultural fields show up as a patchwork driven by different conditions from one field to the next.
Airports
There are some very dark pixels in the image, including those over the Seine River winding through Paris, and a smattering of lakes and other waterbodies dotting the landscape. There are also two distinct linear dark patches just outside the city.
Compare the backscatter imagery to the optical basemap.
The two dark patches are the runway areas of the Charles de Gaulle Airport (CDG). The tarmac is flat and smooth, and the grounds are flat and lack trees or shrubby vegetation, resulting in very low backscatter. The dark area to the southwest of CDG is another airport.
From a distance, these dark areas may appear to be water, but if you compare the backscatter from the airports to those of lakes, the values are often much lower for the lakes.
In this view, click on the dark patch identified as the Aérodrome des Mureaux (a small airfield) and compare the backscatter values to those from the lakes Étang de la Grosse Pierre and Étang du Gallardon.
Turn off the backscatter layer to compare the features to the optical basemap.
Agricultural Regions
The backscatter patchwork pattern of agricultural landscapes changes from one season to the next.
Adjacent fields can have very different backscatter signatures in the summer, depending on the type of crop being grown and irrigation practices. Some fields may be left fallow or be planted with a cover crop, which also impacts their backscatter signature.
The pattern in the spring can be different, especially if some fields are planted earlier than others, but are similar to the summer patchwork in terms of variability.
Fall and winter conditions are a different story. There is much less variability from one field to the next than there was in the spring and summer.
VV vs. VH
For backscatter in VV polarization, the changes in the appearance of the fields are driven largely by soil moisture and surface roughness. The fields are much more uniform in the fall and winter seasons than they are in the spring and summer.
The patterns in VH are a bit different. Click through the seasons to look at the differences.
The patchwork of fields is most obvious in the summer season for the VH backscatter. While the VV polarization had similar degrees of patchiness in the spring and summer, the VH backscatter from the fields is much lower overall in the spring than in any other season.
For the VH polarization, the main drivers of difference tend to be complex vegetation, though soil moisture may also play a role. Many of the fields in the spring VH image may have been prepared and even planted, but might not yet have much complex vegetation, resulting in the lower VH backscatter.
In general, the winter season has the highest VH backscatter values in the agricultural areas of this region. This may be due to a combination of slow-growing winter grain crops and the widespread use of cover crops after harvest throughout France.
Deserts
To wrap up our exploration of backscatter characteristics, let's take a look at some desert areas.
The Sahara and Libyan deserts cover large areas in northeastern Africa. Deserts, like surface water, have extremely low backscatter, but for completely different reasons.
Water has such high dielectric properties that the SAR signal can't penetrate at all, and just bounces off and away from the sensor. Dry sand in the desert has extremely low dielectric properties, and can absorb the incoming signal so well that very little of it is scattered back to its source.
This can be problematic when using backscatter to identify water in arid regions.
In this image, we can see Great Bitter Lake, which is connected to the Red and Mediterranean Seas by the Suez canal. East of the lake is an extremely arid part of the Sinai peninsula. In the mean backscatter products, this sandy expanse actually has even lower backscatter than the lake.
The mean backscatter products are impacted by the presence of ships on the lake, which exhibit high backscatter against the dark background of the lake. This results in a clear indication of the route used by ships navigating the Suez Canal.
To the east, the wavy structures of the sand dunes are visible, as are roads and industrial developments.
Zoom out and pan to explore the dune structure.
AMP Web Map
Explore all of the seasons and polarizations of the backscatter (AMP) dataset using the web map below. Services are grouped by polarization; expand a group to view the different seasons for that polarization and turn specific layers off or on.
Map displaying all of the backscatter (AMP) services for the GSSICB dataset.
STAC Catalog
In addition to the image services, ASF has published a SpatioTemporal Asset Catalog (STAC) catalog for the GSSICB dataset. STAC provides a simple, standardized way to explore geospatial data and metadata. A rich ecosystem of STAC tools and clients is available, including the pystac client library for Python and the STAC Browser UI.
Explore the GSSICB STAC catalog in the STAC Browser , or see sample code for searching the catalog in the Search the GSSICB STAC Catalog Jupyter Notebook .
ASF Services
The GSSICB dataset is a convenient way to explore seasonal dynamics around the globe, and provides ready access to a consistent dataset presenting both coherence and amplitude values. In some cases, however, seasonal averages do not provide the level of temporal resolution necessary for analysis.
If you are interested in generating coherence or amplitude products from specific Sentinel-1 acquisitions, ASF offers On Demand processing services for both RTC and InSAR products. This service is provided at no cost to users, and is a quick and easy way to access analysis-ready Sentinel-1 products.
On Demand processing is available directly in ASF's Data Search Vertex web portal. To access coherence values for specific Sentinel-1 pairs, use the On Demand InSAR option. The output InSAR product package includes a coherence map . To access RTC products, use the On Demand RTC option.