Satellite Ocean Color Data - Wimsoft · 2019. 1. 10. · 1/10/2019 C Wimsoft SeaWiFS sensor on OrbView 2 •SeaWiFS, launched in 1997, still the ^gold standard in ocean color due

ADEOS-II/GLI Jun.-Aug. 2003

Satellite Ocean Color Data

Mati KahruScripps Institution of Oceanographymkahru@ucsd.edu

and

WimSoft, http://www.wimsoft.comEmail: wim@wimsoft.com

1/10/2019 C Wimsoft

1/10/2019 C Wimsoft

CZCS on NIMBUS-7, 1978/11 – 1986/06, 1 km, reprocessed 2007, 2011?

OCTS on ADEOS, 1996/11 – 1997/06, (1 km) 4 km, reprocessed 2006

SeaWiFS, 1997/09 – 2010, 1 km, after 2005 limited LAC (MLAC), GAC,

interrupted data flow from 2008 on

gap, 1986/07 – 1996/10, ~ 10 years (local and interrupted coverage by MOS in 1996)

gap, 1997/07 – 1997/09, ~ 2 months

MODIS on Terra, 1999-, 1 km, 250 m, bad calibrationMERIS on ENVISAT, 2002–2012, 1 km (RR), 300 m (FR)

MODIS on Aqua, 2002/11-present, 1 km, 250 m

VIIRS-SNPP, 2012-present; OLI on Landsat

OLCI-S3A, also MSI on Sentinel-2A & 2B

OLCI-S3B, VIIRS-JSPP1, SGLIPACE, GEO-CAPE, Hyspiri

GLI on ADEOS-2, failed after ~6 months, 2003

1978-

1996-

2000-

2010-

2013

2015

2019

1/10/2019 C Wimsoft C Mati Kahru, 2003-2006

http://www.ioccg.org/sensors/historical.html

1/10/2019 C Wimsoft

http://www.ioccg.org/sensors/current.htmlhttp://ioccg.org/

1/10/2019 C Wimsoft

http://www.ioccg.org/sensors/scheduled.html

1/10/2019

1/10/2019 C Wimsoft

SeaWiFS sensor on OrbView 2

•SeaWiFS, launched in 1997, still the “gold standard” in ocean color due to best orbit (noon), tilting away from the Sun glint, and best calibration

•SeaWiFS/OrbView 2 - the only sensor on the platform while MODIS/Aqua, MERIS/Envisat, etc have many instruments with conflicting requirements

•SeaWiFS calibration is a mixture of in-flight calibration against the Moon plus vicarious calibration versus MOBY

Orbit:Type: Sun-synchronousAltitude: 705 kmEquatorial Crossing: Noon ± 20 minInclination: 98.2°Period: 99 minutes

<- before launch

1/10/2019 C Wimsoft

SeaWiFS, 1997/09 – 2010, ~4 km (GAC) globally, 1 km (merged LAC) available

globally 1997-2010, ~9 km resolution

1/10/2019 C Wimsoft

Processing Level:•Level 0 - raw radiance counts from sensor, spacecraft and instrument telemetry data

•Level 1A -Level-0 data with appended calibration and navigation data, instrument and spacecraft telemetry

•Level 1B – calibrated radiances at sensor level, i.e. TOA

•Level 2 - geophysical values, derived from Level-1A by applying sensor calibration, atmospheric corrections, and bio-optical algorithms, e.g. Chl-a, Rrs412, Rrs443,..

•Level 3 – Level-2 data mapped to a known projection or a regular grid, also “Binned” (in space and time). SMI(Standard Mapped Image) at ~ 4.6 km, 9 km spatial res.

•Level 4 – From a combination of Level-3 satellite data and ancillary information, such as models, e.g. Net Primary Production, Export Flux, etc

1/10/2019 C Wimsoft

Examples:•Level 2 - geophysical values, derived by applying calibrations, atmospheric corrections, and bio-optical algorithms, e.g. chlor_a, Rrs412, Rrs443,..

•Level 3 - Binned (time and space), SMI (Standard Mapped Image) mapped to Equidistant Cylindrical projection, at ~ 4.6 km, 9 km spatial resolution

•Level 4 – After application of models (e.g. Net Primary Production, Export Flux of Carbon) to L3 data

1/10/2019 C Wimsoft

Current Sources of OC data:NASA (& NOAA & USGS) • MODIST, MODISA, VIIRS-SNPP, VIIRS-JSPP1: ~1 km

• OLI on Landsat-8: ~30 m

ESA (& EUMETSAT)• OLCI-3A, OLCI-3B (Sentinels 3A and 3B): 300 m & 1 km

• MSI-2A, MSI-2B (Sentinels 2A and 2B), ~20 m

JAXA • SGLI on GCOM-C: Global 250 m & 1 km

For 250m-1 km resolution: 4 + 2 + (1)

For 15-30 m resolution: 2 + 2

1/10/2019 C Wimsoft

NASA OC: https://oceancolor.gsfc.nasa.gov/

•Level 1&2 browser

• - geophysical values, derived by applying calibrations, atmospheric corrections, and bio-optical algorithms, e.g. chlor_a, Rrs412, Rrs443,..

1/10/2019

https://oceancolor.gsfc.nasa.gov/

•Use Level-1&2 browser to download individual L2 datasets (swaths), or order multiple sets of L2 data

•Use wget to download many global mapped L3 (SMI) data files or browser to download individual files

C Wimsoft

C Wimsoft

Level-1&2 browser

1/10/2019

Paracas2:Lat range: -13.99551 -13.11078Lon range: -76.77642 -76.17332

1/10/2019 C Wimsoft

C Wimsoft1/10/2019

1/10/2019 C Wimsoft

User Services / Order Manageror respond to email

Use Control-R to refresh the screen;Click on and save http_manifest.txt to your tmp folder 1/10/2019 C Wimsoft

• Click on and save http_manifest.txt to your tmp folder• Open cmd.exe, cd tmp• Run dload.batwget --load-cookies ~/.urs_cookies --save-cookies ~/.urs_cookies --content-disposition -i http_manifest.txt

• Note!!! Need to install 7z before!!!Run 7z1805-x64.exe in tools; check that 7z is accessible in tmp path !!!• Run extract.bat7z e *.tarREM del requested_files_*.targzip -d *.gz

1/10/2019 C Wimsoft

1/10/2019 C Wimsoft

•longitude, latitude

•aot_869, angstrom

•Rrs_412, 443, 469, 488, 531, 547, 555, 645, 667, 678

•chlor_a

•Kd_490

•pic

•poc

•cdom_index

•ipar

•nflh

•par

•l2_flags

•In earlier versions: nLw_412, nLw_443, nLw_488, nLw_531, nLw_551, nLw_667, sst

Contents of a MODIS L2 file (R2018):

1/10/2019 C Wimsoft

Sentinel-3: ESA & EUMETSAT. Copernicus Online Data Access (CODA):http://www.eumetsat.int/website/home/Data/DataDelivery/CopernicusOnlineDataAccess/index.html.

https://eoportal.eumetsat.int/userMgmt/newPasswordGenerated.faces

1/10/2019 C Wimsoft

1/10/2019

1/10/2019

1/10/2019

1/10/2019 C Wimsoft

Tools…SeaDASFree, developed by NASA, https://oceancolor.gsfc.nasa.gov/seadas/, essential for L1->L2 processing (calibration, atmospheric correction), source code in C & Python•Cons: no native Windows version (through virtual machine)

BEAN SNAPFree, developed by Brockmann Consult for ESA, https://step.esa.int/main/toolboxes/snap/ , modern user interface written in Java; Cons: primarily for MERIS (Envisat), Sentinel-3 sensors

WIM/WAM• Many high-level applications, e.g. anomalies, fronts, EOF, trend detection, native Windows interface• Cons: not free, https://www.wimsoft.com

ACOLITEFree, developed by RBINS, https://odnature.naturalsciences.be/remsem/software-and-data/acolite; For high-res sensors on Landsat and Sentinel-2A & B, atmospheric correction

1/10/2019 C Wimsoft

Level-2 data:primarily available from data originators (NASA, ESA, JAXA, NOAA), NASA also provides some JAXA, NOAAL2-data data

Level-3 data:Available from data originators (NASA, ESA, JAXA, NOAA) but also from many other data providers and data consolidators:• http://globcolour.info/ - global merged SeaWiFS, MERIS, MODISA,VIIRS-

SNNP• http://hermes.acri.fr/ - ACRI, global merged data• https://www.wimsoft.com/CAL/ - California Current 1 km, 1978-present• https://www.wimsoft.com/CC4km.htm - California Current optimally

merged 4-km

1/10/2019 C Wimsoft

Processing Level:•Level 0 - raw radiance counts from sensor, spacecraft and instrument telemetry data

•Level 1A -Level-0 data with appended calibration and navigation data, instrument and spacecraft telemetry

•Level 1B – calibrated radiances at sensor level, i.e. TOA

•Level 2 - geophysical values, derived from Level-1A by applying sensor calibration, atmospheric corrections, and bio-optical algorithms, e.g. Chl-a, Rrs412, Rrs443,..

•Level 3 – Level-2 data mapped to a known projection or a regular grid, also “Binned” (in space and time). SMI(Standard Mapped Image) at ~ 4.6 km, 9 km spatial res.

•Level 4 – From a combination of Level-3 satellite data and ancillary information, such as models, e.g. Net Primary Production, Export Flux, etc

1/10/2019 C Wimsoft

NASA Ocean Color Level-3 browserhttps://oceancolor.gsfc.nasa.gov/l3/

Better to use batch files and wget from command line (e.g. dload2019.bat); 2019, MODISA, from 1/1/2019, 9 km:

set YEAR=2019set FROM=01-01mkdir ..\%YEAR%wget -q --post-data="sensor=modisa&sdate=%YEAR%-%FROM%&edate=%YEAR%-12-31&dtype=L3m&addurl=1&results_as_file=1&search=A*DAY_CHL_chlor_a_9km.nc" -O - https://oceandata.sci.gsfc.nasa.gov/search/file_search.cgi |wget -i -move A%YEAR%*.nc ..\%YEAR%

1/10/2019 C Wimsoft

NASA Ocean Color Level-3 dataUsing wget from command line (e.g. dload2019.bat); 2019, VIIRS-SNPP, from 1/1/2019, 4 km:set YEAR=2019set START=01-01mkdir ..\%YEAR%wget -q --post-data="sensor=viirs&sdate=%YEAR%-%START%&edate=%YEAR%-12-31&dtype=L3m&addurl=1&results_as_file=1&search=V%YEAR%*DAY_SNPP_CHL_chlor_a_4km.nc" -O - https://oceandata.sci.gsfc.nasa.gov/search/file_search.cgi |wget -i -move V%YEAR%*.nc ..\%YEAR%

Using wget from command line (e.g. dload2019.bat); 2019, VIIRS-JSPP1, from 1/1/2019, 4 km:set YEAR=2019set FROM=01-01mkdir ..\%YEAR%wget -q --post-data="sensor=viirs-jpss1&sdate=%YEAR%-%FROM%&edate=%YEAR%-12-31&dtype=L3m&addurl=1&results_as_file=1&search=V*DAY_JPSS1_CHL_chlor_a_4km.nc" -O - https://oceandata.sci.gsfc.nasa.gov/search/file_search.cgi |wget -i -move V*_chlor_a_4km.nc ..\%YEAR%

1/10/2019 C Wimsoft

NASA Ocean Color Level-3 dataCheck product version:wam_attribute 2019\A2019*_9km.nc SHOW processing_version

Julian 2018 330Julian Day=330 => Year=2018 Month=11 Day=26set YEAR=2019set FROM=11-26mkdir ..\%YEAR%wget -q --post-data="sensor=modisa&sdate=%YEAR%-%FROM%&edate=%YEAR%-12-31&……..

1/10/2019 C Wimsoft

https://www.wimsoft.com/CAL/

1/10/2019 C Wimsoft

http://www.globcolour.info/FTP service: ftp.acri.fr , login globcolour_data , password fgh678

Directories: GLOB_4KM, EURO_4KM, EURO_2KM

1/10/2019 C Wimsoft

FTP service: ftp.acri.fr , login globcolour_data , password fgh678 Directories: GLOB_4KM, EURO_4KM, EURO_2KMDirectory structure under FPS:

Contents of a GLobColour 4 km CHL1 file:

CHL1_mean ↓ in ISIN projection (4.7 km)