POLDER I-II and POLDER/PARASOL instruments
The POLDER instruments consist of a digital camera with a 274×242-pixel CCD detector array, wide-field telecentric optics and a rotating filter wheel enabling measurements in 9 spectral channels with bandwidths between 20nm and 40nm. Because it acquires a sequence of images every 20 seconds, the instrument can observe ground targets from different view directions. The POLDER I-II instruments onboard ADEOS 1 and 2 are identical while the instrument on the PARASOL platform has been turned 90 degrees to favor multidirectional viewing (maximum of 16 directions compared to 14) over daily global coverage (swath of 2400 km compared to 1600 km). Depending on the altitude of the platforms, the size of the images varies from 2400 x 1800 km2 to 1600 x 2100 km2 with a corresponding ground resolution of 7×6 km2 and 5.3×6.2 km2 at nadir. The PARASOL platform is part of the A-Train and takes advantage of the other instruments in the constellation. The 3 instrument spectral coverage ranges from blue (443µm) through near-infrared (0.91µm) with 3 polarized spectral bands. For POLDER/PARASOL, the bluest polarized channel has been moved from 0.443µm to 0.490µm and a 1.02 µm waveband has been added.
GRASP (Generalized Retrieval of Aerosol and Surface Properties) algorithm is designed to retrieve complete aerosol and surface properties globally. In order to achieve reliable retrieval from satellites observations even over very reflective desert surfaces, the algorithm was designed as simultaneous inversion of a large group of pixels within one or several images. Such, multi-pixel retrieval regime takes advantage from known limitations on spatial and temporal variability in both aerosol and surface properties. Specifically the variations of the retrieved parameters horizontally from pixel-to-pixel and/or temporary from day-to-day are enforced to be smooth by additional appropriately set a priori constraints. This concept provides satellite retrieval of higher consistency. The details can be found in the articles by Dubovik et al., 2011, 2014 and in the technical documentation.
POLDER/GRASP aerosol and surface date product
NOTE: the same PARASOL/GRASP data but in somewhat in different format can be obtained from AERIS/ICARE Data and Services Center (http://www.icare.univ-lille1.fr/) that is official distributor of POLDER Level-1 and Level-2 data. Also, ICARE center provided detailed visualization of retrievals.
Bellow there are two archives of POLDER/PARASOL data processed by GRASP:
(1) PARASOL/GRASP «optimized» (in sense that radiative transfer calculations were optimize to best tread-off between speed of processing and accuracy of results);
(2) PARASOL/GRASP «high-precision» (the accurate radiative transfer calculations were used).
The results from «optimized» and «high-precision» processing are quite similar, however the analysis and validation against AERONET showed that «high-precision» overall provide more accurate retrievals (especially for coarse mode dominated aerosols as dust) and, therefore, these data can be considered as recommended for the users.
This data sets includes the aerosol and surface reflectance retrieval product obtained from the POLDER/PARASOL data processed by the GRASP code.
The POLDER/PARASOL data including measurements of all angular measurements intensity at 0.443, 0.490, 0.565 , 0,670, 0, 765, 0.865 and a 1.02 µm and Q and U polarization component at 0.490, 0,670 and 0.865 at the instrument native resolution were inverted in cloud-free conditions as determined by the original cloud-mask algorithm developed by the Laboratoire d’Optique Atmosphérique. The data at polar regions at latitudes higher than 70 degrees north and south were not processed. The observations with less than 10 viewing directions were ignored.
The retrievals were performed using one aerosol component model with 5 size bins size distribution and spectrally dependent complex refractive index. The aerosol vertical distribution was modelled using exponential profile and scale height retrieved. The “optimized” settings (inversion strategy used by GRASP) were used in RT modeling of multiple-scattering in the atmosphere. The POLDER data were processed using GRASP v0.6.5. The full list of the retrieved parameters, as well as additional characteristics derived from retrieved parameters, can be found in the product list provided below. The further details of the retrieval setting can be provided by request.
The raw results files from GRASP contains three months of data of small regions (162×162 pixels). For the convenience of the users, the daily global files were created as well as their aggregates in time and space. Only the most reliable and demanded data were exported for final public distribution. Data format of this public archives is NetCDF.
In order to assure higher quality of the data, some post-processing was applied in order to eliminate the low quality points resulted from cloud contamination, bad surface description near the coast, etc. The post-processing screening follows these steps:
- First, from raw GRASP output we create global daily files. We call them level 1.
- Pixels with AOD443 > 10 are removed.
- The coast is removed so all pixels with land percent between 1 and 99 are removed. Also, to guarantee a proper coast elimination, the first pixel into ocean and land is removed.
- We remove unphysical values like water surface model over land and the other way around.
- We apply a criteria over to screen bad pixels (retrieval error is high so we remove completely the pixel):
- ocean → (residual <= 0.13)
- If precondition ndvi < 0.1
- if (dhr670 >= 0.3) then
- if (AOD670<1)
then (residual <= 0.04)
- else (residual <= 0.06)
- if (AOD670<1)
- (dhr670 < 0.3) -> (residual <= 0.07)
- if (dhr670 >= 0.3) then
- # precondition 0.1 <= ndvi < 0.4
- (dhr670 >= 0.25) -> (residual <= 0.075)
- (dhr670 < 0.25) -> (residual <= 0.085)
- # precondition 0.4 <= ndvi < 0.6
- (residual <= 0.1)
- # precondition 0.6 <= ndvi < 1
- (residual <= 0.12)
- If precondition ndvi < 0.1
The date passed all above filter files are stored under level 1.5. These data are not publicly available, but can be provided by a special request. This archive ensure good surface retrieval but quality of aerosol products is not guaranteed.
- the outliers are screened analyzing groups of 20×20 pixels and iteratively we remove the worst pixel (farther to the mean of AOD870) if the group does not fulfil the following condition “std of data <= 1.5 AND std / mean of data <= 0.5”. If the result group has less than 60 pixels it is completely removed (noisy area probably due to cloud contamination).
- Some aerosol products can only be calculated when there is some aerosol loading (otherwise it is difficult to measure them). So, we apply extra filters. Angstrom Exponent is only provided if AOD560 is higher than 0.02 over ocean or 0.2 over land. Other more complex products (SSA, Re(m), Im(m), AAOD,Size distribution,SphereFraction) are strongly filtered:
- land: AOD443 >= 0.3 and 0.65 <= SSA <= 1.
- ocean: AOD443 >= 0.02 and 0.65 <= SSA <= 1.
- These results are exported as level 2. Level 2 is also temporally aggregated. The temporal aggregation pixels are removed if AOD443>4.
Level 3 is created as regridding products at 0.1 and 1 degrees resolution in WGS84 projection of data available in level 2. In the regridding process we applied median filter instead of average or any other sophisticated filter.
In summary, the following data sets are prepared:
Level 0: raw results from grasp
Level 1: daily files from the output
Level 1.5: data softly screened. All surface pixels are good but we cannot guarantee the quality of aerosol information
Level 2: Full resolution data filtered and aggregations (daily, monthly, yearly, seasonal, and climatologically monthly, seasonal and full archive).
Level 3: Regrid at 0.1 and 1 degree of level 2 (including temporal aggregations).
Only levels 2 and 3 are publicly available since levels – 0 to 1.5 are considered internal files results of intermediate processings. Climatological monthly aggregation is the aggregation of all Januarys, Februaries, … in the archive, as seasonal climatology is all winters, all springs, …. .
The access to the internal archive could be provided under request. Please, contact us for further information.
Any publication, presentation, or other derivative work based on these results shall acknowledge its use. Example of acknowledgment:
For PARASOL/GRASP «optimized»: The authors would like to acknowledge the use of POLDER data "POLDER/PARASOL Level-1 data originally provided by CNES (http://www.icare.univ-lille1.fr/) processed at Laboratoire d'Optique Atmosphérique with GRASP software (https://www.grasp-open.com) developed by Dubovik et al. (2011, 2014). Short credit: "Credit: CNES/GRASP/LOA"
For PARASOL/GRASP “high-precision”: The authors would like to acknowledge the use of POLDER data "POLDER/PARASOL Level-1 data originally provided by CNES (http://www.icare.univ-lille1.fr/) processed at AERIS/ICARE Data and Services Center with GRASP software (https://www.grasp-open.com) developed by Dubovik et al. (2011, 20114). Short credit: "Credit: CNES/GRASP/LOA/ICARE"
Dubovik, O., T. Lapyonok, P. Litvinov, M. Herman, D. Fuertes, F. Ducos, A. Lopatin, A. Chaikovsky, B. Torres, Y. Derimian, X. Huang, M. Aspetsberger, and C. Federspiel: GRASP: a versatile algorithm for characterizing the atmosphere, SPIE: Newsroom, Published Online: September 19, 2014. doi:10.1117/2.1201408.005558 Dubovik, O., Herman, M., Holdak, A., Lapyonok, T., Tanré, D., Deuzé, J. L., Ducos, F., Sinyuk, A., and Lopatin, A.: Statistically optimized inversion algorithm for enhanced retrieval of aerosol properties from spectral multi-angle polarimetric satellite observations, Atmos. Meas. Tech., 4, 975-1018, 2011. doi:10.5194/amt-4-975-2011