QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV WP4 T4.2 - Development of harmonised retrieval for surface BRDF/albedo, D4.2, D4.3
EUMETSAT, UCL, BC, NPL, JRC, RF
D4.2, D4.3 EUMETSAT, UCL, BC, NPL, JRC, RF QA4ECV Mid-Term Review - - PowerPoint PPT Presentation
QA4ECV WP4 T4.2 - Development of harmonised retrieval for surface BRDF/albedo, D4.2, D4.3 EUMETSAT, UCL, BC, NPL, JRC, RF QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016 BRDF/Albedo Plan A Overview Development of GlobAlbedo approach
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
EUMETSAT, UCL, BC, NPL, JRC, RF
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Approximate optimal solution using weighting on Cobs in time
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Formal linear optimisation framework: g smoothness control
Dx at day and year period
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Low pass filter
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
uncertainty
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Linear mapping to common spectral basis Hainich
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
MERIS l1b VGT l1p
L1 → SDR Orbit → Tile
SDR
by /sensor/tile/year/day/ti me/
SDR → BBDR/MODIS (Spectral to BB)
BBDR → BB BRDF SDR → Spectral BRDF
MOD/MYD09 (MODIS SDR) SDR → SDR/MODIS
(Spectral mapping)
BRDF-Prior
(Spectral & BB) by /tile/day/
BRDF → Albedo
Albedo (Spectral & BB)
by /tile/year/day/ & /mosaic/year/day/
Proba-V l1b
AVHRR/GEO SDR BB
SDR BB → BBDR/MODIS
(Broadband mapping )
SDR BB
by /sensor/tile/year/da y/
Orbit/disk → Tile
BRDF (Spectral & BB)
by /tile/year/day/
Tiles → Mosaics
(Upscaling by energy conservation)
BRDF (Spectral & BB)
by /mosaic/year/day/ Weights by /sensor/modis- band/time/ SDR/MODIS→ BBDR/MODIS (Spectral to BB) (Weights creation) prior creation
Uniformity
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Ingestion of AVHRR-LTDR surface reflectance products
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Ingestion of AVHRR-LTDR surface reflectance products
which have an intersection with the subset
chain and has been used for MERIS and VGT in the BBDR retrieval. Adaptation to other sensors is straightforward.
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
From global product to latitude subset:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Reprojection of latitude subset onto MODIS SIN tiles:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Ingestion of AVHRR-LTDR surface reflectance products QA quality flag coding:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Ingestion of Meteosat GEO MVIRI+SEVIRI surface reflectance products
as seen from the satellites’ Equatorial position at different longitudes
special geolocation algorithm to handle off-planet pixels in the disk product.
products , their computation was implemented separately
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
From ‘disk’ to lat/lon grid:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
MODIS SIN tiles within lat/lon area:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
GEO BRF on 4 MODIS SIN tiles:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
AVHRR + Meteosat GEO MVIRI+SEVIRI: From BRF to BBDR to BRDF/albedo
reflectances: BB_VIS, BB_NIR, BB_SW
broadbands using coefficients provided by Liang et al. (2000). Corresponding uncertainties are computed from simple error propagation arithmetics
using the method described by Loew and Govaerts (2010)) is taken as BB_SW
separately (using AVHRR or MVIRI/SEVIRI samples only) or by combining the samples
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Albedo from AVHRR and/or MVIRI: examples
tiles for the 8-day periods in year 2005
desert) and h19v08 (tropical rain forest) for 2005, DoY 121. Comparisons with GlobAlbedo results retrieved from MERIS+VGT
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h18v06. Left: AVHRR; centre: AVHRR+MVIRI; right: GlobAlbedo
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h19v08. Left: AVHRR; centre: AVHRR+MVIRI; right: GlobAlbedo
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h18v06 (desert): AVHRR and AVHRR+MVIRI vs. GlobAlbedo MERIS+VGT
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h19v08 (rain forest): AVHRR and AVHRR+MVIRI vs. GlobAlbedo MERIS+VGT
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h18v06: MSA_Albedo standard product (left) vs AVHRR+MVIRI (right)
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h19v08: MSA_Albedo standard product (left) vs AVHRR+MVIRI (right)
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
BHR/DHR h18v06/h19v08: AVHRR+MVIRI vs. MSA_Albedo standard product
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Next steps
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
MERIS l1b VGT l1p
L1 → SDR Orbit → Tile
SDR
by /sensor/tile/year/day/ti me/
SDR → BBDR/MODIS (Spectral to BB)
BBDR → BB BRDF SDR → Spectral BRDF
MOD/MYD09 (MODIS SDR) SDR → SDR/MODIS
(Spectral mapping)
BRDF-Prior
(Spectral & BB) by /tile/day/
BRDF → Albedo
Albedo (Spectral & BB)
by /tile/year/day/ & /mosaic/year/day/
Proba-V l1b
AVHRR/GEO SDR BB
SDR BB → BBDR/MODIS
(Broadband mapping )
SDR BB
by /sensor/tile/year/da y/
Orbit/disk → Tile
BRDF (Spectral & BB)
by /tile/year/day/
Tiles → Mosaics
(Upscaling by energy conservation)
BRDF (Spectral & BB)
by /mosaic/year/day/ Weights by /sensor/modis- band/time/ SDR/MODIS→ BBDR/MODIS (Spectral to BB) (Weights creation) prior creation
Uniformity
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Surface reflectance from sat X Surface reflectance from sat MODIS
Scan for samples close in angle and time and over homogenous and flat areas
MODIS Snow MODIS Landcover DEMs BRDF prior
When CoV is low use Weighted least squares regression: 1st and 2nd order Stage-1 Stage-2
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Linear least-squares Quadratic least-squares
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Tile: h20v06 (Egypt) Date_Modis: 11 Aug 2011 Date_Meris: 08 Aug 2011 Difference in angles: <5deg MERIS is transformed into the equivalent MODIS Spectral bands and compared. Example: desert MODIS MERIS MERIS vs MODIS
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
MERIS l1b VGT l1p
L1 → SDR Orbit → Tile
SDR
by /sensor/tile/year/day/ti me/
SDR → BBDR/MODIS (Spectral to BB)
BBDR → BB BRDF SDR → Spectral BRDF
MOD/MYD09 (MODIS SDR) SDR → SDR/MODIS
(Spectral mapping)
BRDF-Prior
(Spectral & BB) by /tile/day/
BRDF → Albedo
Albedo (Spectral & BB)
by /tile/year/day/ & /mosaic/year/day/
Proba-V l1b
AVHRR/GEO SDR BB
SDR BB → BBDR/MODIS
(Broadband mapping )
SDR BB
by /sensor/tile/year/da y/
Orbit/disk → Tile
BRDF (Spectral & BB)
by /tile/year/day/
Tiles → Mosaics
(Upscaling by energy conservation)
BRDF (Spectral & BB)
by /mosaic/year/day/ Weights by /sensor/modis- band/time/ SDR/MODIS→ BBDR/MODIS (Spectral to BB) (Weights creation) prior creation
Uniformity
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Time Range Dataset Satellite DOI 1982-2006 Prime at 0° Meteosat-2,7 10.15770/EUM_SEC_CLM_0001 1998-2006 Indian Ocean Data Coverage at 63°E Meteosat-5 10.15770/EUM_SEC_CLM_0002 2006-ongoing Indian Ocean Data Coverage at 57°E Meteosat-7 10.15770/EUM_SEC_CLM_0003 1991-1983 Atlantic Data Coverage at 50°W Meteosat-3 10.15770/EUM_SEC_CLM_0004 1993-1995 Extended Atlantic Data Coverage at 75°W Meteosat-3 10.15770/EUM_SEC_CLM_0005
Meteosat First Generation have different SSR. For comparing them a spectral conversion is necessary. The method by Loew and Govaerts, 2010 has been applied
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Processing of MVIRI and SEVIRI HRVIS data for generating BRF after 2006 (QA4ECV)
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
ASSUMPTION: The RPV model parameters (K, THETA, R0) retrieved for one day within a 10 days period are assumed to be applicable to the all 10
changing surface properties. The local noon Sun Zenith Angle is calculated and used for estimating the BRF for each day.
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Nadine Gobron
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
52
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
53
use actual sensor measurements but algorithm depend 1) on past sensors data and results or 2) relationship are made with in-situ measurements.
are based
statistical methods
nor ‘generic’ by nature … Hybrid
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Surface level
Reflectance Factors (BRFs) from multiple single (angular) spectral bands in visible to near- infrared (NIR) bands
reflectance in visible to near-infrared (NIR) bands
NIR broadband. Top of atmosphere level
bands in visible to near- infrared (NIR) bands
54
Time scale
scale inputs) Spatial Scale
Definition
acquisition, fixed sun zenith angle) and/or diffuse.
Assumption
‘Green’, ‘Land Cover dependency’, ‘Polychrome’
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
FAPAR comes from energy balance: it is always derived from a canopy RT model: Surface albedo products and absorption should be physically consistent QA4ECV LAI/FAPAR algorithms are generic:
discrete models)
benchmark (RAMI and RAMI4PIPLS)
55
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
56
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
– Fast computations – Applying to each single acquisition image (generic algorithm)
– Take into account sensors spectral properties – Take into account geometries of illumination and view:
effects
Continuity of ESA MERIS/OLCI products – Instantaneous ‘Green’ FAPAR – Rectified values over both vegetated and soil surfaces. JRC retrieval method assumes that the leaves are alive and photosynthesizing, hence the name “green” FAPAR. It also means that the single scattering albedo of leaves is “fixed” to only one value representing such ‘green’ leaves. JRC-FAPAR also refers to the instantaneous and green value definition. The theoretical FAPAR values are computed with RT model using the closure of the energy balance inside the plant canopy in the spectral range 400 to 700 nm.
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016 58
15 January 2015
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016
Optimization of RPV parameters and polynomial coefficients are achieved with RT canopy model simulations
VZA BRF
Polynomial FAPAR
VZA BRF
59
10 spectral responses + One virtual VIS for FAPAR simulations
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016 60
# canopies within TOA version, i.e. MERIS, OLCI, etc. More LAI values
After RPV
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016 61
where σx
2 and σy 2 are the standard
deviation error of the rectified channels in Band 1 and Band 2, respectively.
We use here the specifications of the AVHRR surface reflectance accuracies (Claverie and Vermote, 2015, personal communication) which correspond to the double of the MODIS surface reflectance specifications ones:
QA4ECV Mid-Term Review Meeting, MSSL, 14-16 June 2016 62