Sea water turbidity analysis from Sentinel-2 images: atmospheric correction and bands correlation
This is a modal window.
The media could not be loaded, either because the server or network failed or because the format is not supported.
Formal Metadata
Title |
| |
Title of Series | ||
Number of Parts | 351 | |
Author | ||
Contributors | ||
License | CC Attribution 3.0 Unported: You are free to use, adapt and copy, distribute and transmit the work or content in adapted or unchanged form for any legal purpose as long as the work is attributed to the author in the manner specified by the author or licensor. | |
Identifiers | 10.5446/68964 (DOI) | |
Publisher | ||
Release Date | ||
Language | ||
Production Year | 2022 |
Content Metadata
Subject Area | ||
Genre | ||
Abstract |
| |
Keywords |
00:00
AreaMathematical analysisCross-correlationPresentation of a groupComputer animation
00:19
SatelliteMathematical analysisObservational studyMusical ensembleCross-correlationComputer fontMathematical analysisSatelliteObservational studyMappingInformationCASE <Informatik>Computer animation
00:56
MathematicsWater vaporNegative numberMeasurementParameter (computer programming)Sound effectParticle systemIntegrated development environmentComputer animation
01:45
InformationMatter waveAngleExecution unitScatteringDegree (graph theory)SpacetimeMeasurementSurfaceElectronic signatureSpectrum (functional analysis)Function (mathematics)VotingSatelliteInformation retrievalOpen sourceSoftware frameworkStrategy gameMaß <Mathematik>Observational studyAreaMultiplicationCompilation albumProduct (business)Musical ensembleVaporAngular resolutionMusical ensembleExecution unitSet (mathematics)InformationView (database)AreaSatelliteSurfacePosition operatorCondition numberMedical imagingDifferenz <Mathematik>Water vaporEvent horizonComputer programProduct (business)HistogramObservational studyIntegrated development environmentRemote procedure callMultiplication signMeasurementImage resolutionObject (grammar)Dependent and independent variablesElectronic signatureParameter (computer programming)TesselationComplex systemMetreAngleMatter wavePoint (geometry)Information retrievalReflection (mathematics)Component-based software engineeringVotingSpacetimeFood energyDegree (graph theory)Greatest elementParticle systemMusical ensembleState observerType theoryOpen sourceMereologyIncidence algebraGreen's functionSpectrum (functional analysis)ScatteringComputer animation
06:49
Process (computing)ScatteringProduct (business)Level (video gaming)GeometryGreatest elementObservational studyOpen sourceGreatest elementParticle systemElectronic signatureEnergy levelSound effectReflection (mathematics)Computer animation
07:42
SoftwareGrass (card game)Density of statesInformationChainProcess (computing)DataflowProduct (business)ImplementationMusical ensembleSystem identificationPrice indexPoint cloudGroup theoryDifferenz <Mathematik>Computer-generated imageryLogical constantPixelSatelliteSpectrum (functional analysis)SimulationRayleigh-ZahlSound effectObject (grammar)Characteristic polynomialLevel (video gaming)Set (mathematics)Numbering schemePresentation of a groupInformationMedical imagingGrass (card game)Musical ensembleUniform resource locatorGreatest elementReflection (mathematics)PixelSatelliteObject (grammar)Particle systemSoftwareResultantPoint cloudProduct (business)Water vaporLibrary (computing)Energy levelSound effectDifferenz <Mathematik>WordMappingMultiplication signOpen sourceComputer animation
09:59
Standard deviationService (economics)Population densityMetadataDiscrete element methodProduct (business)Digital signaloutputParameter (computer programming)Parameter (computer programming)Arithmetic meanSoftwarePixelComputer fileOpticsProduct (business)Multiplication signMedical imagingPosition operatorGeometryService (economics)Level (video gaming)InformationGrass (card game)MetadataComputer animation
11:15
Workstation <Musikinstrument>WebsiteDistanceParameter (computer programming)Characteristic polynomialComputer networkMeasurementRoboticsOpen sourceWebsiteRoboticsCharacteristic polynomialGrass (card game)SoftwareMeasurementParameter (computer programming)AreaCASE <Informatik>Workstation <Musikinstrument>Observational studyComputer animation
11:59
Cellular automatonSound effectParameter (computer programming)Workstation <Musikinstrument>Statistical dispersionComputer-generated imageryLogical constantMusical ensembleProcess (computing)PixelMeasurementPairwise comparisonCross-correlationMusical ensembleSingle-precision floating-point formatMatter waveMedical imagingCross-correlationParameter (computer programming)Greatest elementWorkstation <Musikinstrument>CASE <Informatik>Sound effectSet (mathematics)Square numberGrass (card game)Reflection (mathematics)SoftwareError messageSubject indexingInformationCoefficient of determinationArithmetic meanEnergy levelFehlerschrankeMusical ensembleMusical ensembleValidity (statistics)Right angleRootVariable (mathematics)BendingRange (statistics)Resultant1 (number)Computer animationDiagram
14:46
Compilation albumEnergy levelMusical ensembleLevel (video gaming)Set (mathematics)Medical imagingSound effectMappingRight angleMaxima and minimaWater vaporCartesian coordinate systemShape (magazine)Computer animation
16:34
Product (business)Latent heatComputer-generated imageryOpticsInformationGreatest elementLimit (category theory)Sound effectObservational studyCross-correlationMaxima and minimaImage resolutionSpacetimeAreaVirtual machinePredictionImage resolutionPoint cloudCondition numberSoftware developerInformationGreatest elementSpacetimeMeasurementSet (mathematics)State observerGoodness of fitProduct (business)Observational studySound effectCore dumpComplex analysisOpticsPredictabilityVirtual machineSystem callCross-correlationReflection (mathematics)CASE <Informatik>Musical ensembleComputer animation
Transcript: English(auto-generated)
00:00
Thank you Bianca for your kind presentation and thank you everyone for attending. So yes, today I will present our work that was done to obtain, to analyze sea water turbidity from Sentinel-2 images.
00:20
So first of all I will say something about turbidity, what turbidity is and why we study for water quality purposes. I will recall briefly some satellite remote sensing basics information that make it useful for water quality analysis and then I will go deeper into our case study.
00:42
I will present all the steps for the study that goes from atmospheric correction to band analysis up to the realization of turbidity maps. I will show some results and get some conclusions. So what is turbidity? I think we're all pretty familiar about turbidity but what is it exactly?
01:03
Turbidity is a measure of the amount of light scattered by the particles that are in water. For sea water, changes in turbidity are due to, for example, algal blooms but the most relevant contribution to changes in turbidity is due to sediment.
01:22
For example, sediment coming from river flows. Especially when related to anthropogenic activities like beach nourishment or dredging activities, the increase of turbidity can be associated to negative environmental effects. That's why turbidity is one of the parameters that is considered for European legislation
01:42
for water quality monitoring. In situ, turbidity monitoring is done by a probe. As said, turbidity is not the only measure that is taken so usually it's measured with a multi-parametric probe that also measures other parameters.
02:01
The measured turbidity depends of course by the particles that are present in water but also on the instrument, on the wavelength of the light and the angle at which the detector is positioned. For example, nephilometric turbidity unit, NTU, measures the scattered light at 90 degrees from the incident light beam measured with a white light.
02:25
In situ monitoring offers high accuracy data but it is of course time consuming and offers information that is sparse both in time and space. On the other hand, remote sensing and Earth observation can provide fast and economic
02:44
when we consider free available data and most of all a synoptic view of the situation. Remote sensing and in particular optical remote sensing measures the solar energy that is reflected by the Earth's surface.
03:02
The different objects and features of the Earth's surface have a unique spectral response called spectral signature. And this is the fact that we can use to study seawater surface. In fact, the seawater surface has a spectral signature that depends on the optical significant
03:27
constituents. There are phytoplankton pigments, color-disturbed organic matter and suspended matter. We can see in the right image the spectral signature of the different components and we can notice that water absorbs energy of most of the wavelength while the presence
03:47
of chlorophyll causes a small peak of reflectance in the green area. But sediments is the particle that causes the higher values of reflectance not only
04:04
in the visible spectrum but also in the near infrared region. So the aim of the study was to retrieve turbidity information in a study area. The study area is here reported. It's in the North Iranian Sea that includes Liguria region and Tuscany region.
04:25
You can see here the study area and it's a very complex area from the turbidity point of view because it comprehends many different sources of turbidity. We have in the northern part a marine protected area. So it's supposed to be a very clean type of water.
04:42
But then we have a harbor. We have two river mouth, Magra river mouth and Arno river mouth. And again we have another harbor that is Livorno harbor.
05:03
So it's a very complex system. We use in situ turbidity measure collected from 2015 up to 2021 by the regional agency for the environmental protection of Liguria where I work and Tuscany. According to the official monitoring programs.
05:21
Turbidity is measured in a philometric turbidity units. You can see on the bottom left image the histogram of the turbidity data. This data set is not just the data that were collected in these six years but the data corresponding to the availability also of satellite image
05:44
because satellite image is not always available. We can immediately notice how the data set is very imbalanced towards low data and especially data lower than three NTU which means clear water. This is also due to the fact that the high turbidity events are associated
06:03
to intense meteorological events during which the sea condition does not allow in situ monitoring or the sky is cloudy so the satellite information is not available.
06:24
We use the Sentinel-2 multispectral instrument products because of its time and special resolution. Especially in the visible bands we have a special resolution of 10 meters. The products are available in areas of fixed size called tiles
06:41
and our study area is totally included in one tile. The revisit time in this area is about three days. Sentinel-2 provides data since 2015 but only since March 2018 the data are provided as so-called level 2 that is bottom of atmosphere reflectance.
07:02
We know in fact that the atmosphere interfere with the sunlight and changes the spectral signature. So the effects of the atmosphere are basically absorption. For example UV rays are blocked by ozone while carbon dioxide stops the thermal infrared
07:21
and by scattering due mostly by aerosol particles. So the so-called top of atmosphere reflectance which is the one that is measured by the sensor is composed by the bottom of atmosphere or surface reflectance plus the contribution of atmosphere.
07:44
So this is the scheme of our work. So first of all we collected of course the in-situ data and we collected the satellite data for the same day. As said we had to remove all the data with the presence of cloud
08:01
and then because not all the Sentinel data were available as a bottom of atmosphere information we had to do some atmospheric correction. And three open source softwares were tested and compared the QGIS, the GRASS GIS and the SNAP. And the results of the bottom of atmosphere obtained reflectance were compared
08:21
to the level 2 Sentinel-2 data because we wanted to have a homogeneous data set to use. So once the best method was identified and the whole data set was used to identify which bands were more correlated to turbidity and we only considered the bands that are known to be correlated to particles in water
08:42
that we saw before. So we found calibrated and validated a turbidity model that was used finally to obtain some turbidity maps. About atmospheric correction, just a few words. There are two main groups of atmospheric correction methods. The first one and is the easiest one are the image-based approaches.
09:02
They consist of the subtraction of a constant value from each pixel of the processed image. So this is a faster way because it doesn't require any information about the atmosphere. Basically it is based on the idea that the darkest pixel in the image
09:20
have a reflectance that is the lowest reflectance of the image that is due not because of the reflectance of the object but the effect of the atmosphere. So that is the value that is considered at the atmospheric effect and is removed by the other pixel. The physically based approach instead, they rely on the use of a radiative transfer model
09:42
to accurately compute the atmospheric effect. Examples are the 6S method which is used by GRASS GIS and the LibRATRAN method which is used by SNAP. Of course they need some atmospheric information at the time and location of the image acquisition.
10:01
So what parameters are required for the two methods? Of course some sensor and solar geometries information, the date and time of acquisition and the position of the image. But these data are easy to find because they are provided by the metadata file. We need the mean target elevation above sea level which is usually derived by
10:22
digital elevation model. And then we need the most difficult one are the data related to the atmosphere. So atmospheric model and aerosol model. A very important parameter is the aerosol optical depth or thickness at 550 nanometers. This is a dimensionless parameter that is related to the amount
10:43
of aerosol that is present in the vertical column of the atmosphere above the target. The higher the IOD, the less sunlight passes through the atmosphere. So for SNAP, this is a parameter that is automatically calculated by the
11:00
software from visibility using dark reference areas. If no dark pixels are present in the image, a forecast data obtained from Copernicus atmosphere monitoring service is used. While from GRASS GIS, we used an IOD parameter estimated from the aerosol robotic network.
11:23
It is possible to search for the nearest measurement site or the site that has the atmospheric characteristic closest to the study area. In our case, we identified five stations, Ispra and Modena because they are the closest to the study area,
11:42
but they're not coastal sites so they have high altitude. So we also analyzed three other station, Rome, Naples, and Palma de Mallorca in Spain, because they are at the coast, so they should have more similar characteristics. From the image on the left, we have an example of the AOD daily variability
12:05
for the different stations. What we can immediately notice is that the values are not constantly available. For example, for that date, the data from Naples stations were not available. On the right image, we can see the effect of the AOD parameter
12:21
on the atmospheric correction in GRASS GIS. Again, for just one day, an example, and one band. What we can notice is that changing the AOD parameter changes a lot the bottom of atmosphere reflectance obtained. And most of all, we can see that the AOD parameter of 0.15,
12:42
which is the one that we would use by looking at the stations, is not actually the one to give us the best bottom of atmosphere values compared to the level 2 Sentinel-2. So results obtained with the different softwares were compared to the level 2 data.
13:06
We can say that the QGIS is very fast, but is not very accurate, at least for this case. While the GRASS GIS software, which is a physically based approach, it can be very accurate, but in our case,
13:21
because we didn't have information about the AOD parameter, it's very difficult to use. On the other hand, SNAP, because it evaluates automatically, the AOD parameter was very efficient and actually accurate when compared to level 2 data. So the bottom of atmosphere reflectance obtained either directly from level 2
13:46
or after atmospheric correction was then used to identify which bands are more correlated to turbidity. Again, I said just the visible and near infrared bands were considered. And we can see here that the bands in the visible range
14:01
are the ones that are most correlated to turbidity. They were used to identify an index, were combined together to identify an index to relate it to turbidity. The correlation is reported here. To validate the correlation, a validation data set was used,
14:23
and you can see here the indexes R square, root mean square error, and mean absolute error were calculated to evaluate the accuracy. And you can see that the correlation that was found is linear with respect to the index,
14:41
but is actually non-linear with respect to the single bands. So the relationship was finally used to obtain turbidity maps. You can see in the left image the RGB composite or true color image for one specific day.
15:00
And on the right image, we have the turbidity map corresponding to the same day. Well, of course, we can notice that the land is evaluated as very high turbidity levels. Yes, of course. But we can clearly see that the shapes of the turbidity plumes are well represented.
15:21
The values are good values. I mean, they are physically good values. But even though they are very correlated to the calibration data set that we had, I mean, we cannot think to estimate values that are lower than the minimum value that we had,
15:41
or higher than the maximum value that we had for the calibration data set. This is another example of a turbidity map. Again, we can see how the turbidity plumes are well represented. I just wanted to recall that this is the Arno River, and we have here. Another river.
16:01
What we can see, though, in this image is that there are some effects due to the image, the Sentinel-2 image. We can also notice, not so well, but we can also notice them in the RGB composite. This is probably due to the fact that the Sentinel-2 image is actually not thought for application on water.
16:21
For example, they don't consider sun-glint correction, which is done instead with the Sentinel-3 ocean land color instrument. To conclude, some conclusion about the atmospheric correction. So as we said, the semi-empirical physically based model to obtain
16:43
water quality information from optical data required the use of bottom of atmosphere reflectance, which has to be obtained from atmospheric correction. From our experience, the physically based methods are performed better. However, the IOD, or in general,
17:03
the atmospheric information that are required can make it difficult to apply. Sentinel core from ESA was found to be the most suitable for our case study. To identify correlation, we think we perform good. We can actually say that earth observation can be an important tool
17:24
to support water quality monitoring. It provides especially wide information at significant resolution in space and time. The relationship shows a good agreement between measures and prediction. However, as said, it is highly dependent on the data set available.
17:43
Of course, it could benefit from new data, but they are very hard to obtain because they are limited to the sea condition and from cloud condition. To limit the effects of the atmosphere correction and improve the complexity of
18:01
the relationship, for example, including more bands, machine learning techniques are now under development. And before concluding and thanking you all, I would like to make an open call saying that if any of you works with water quality and turbidity, it's very welcome to cooperate and test our model with your data.
18:22
Thank you so much.