SPECTRA BARRAX CAMPAIGN SPARC

Transcription

1 SPECTRA BARRAX CAMPAIGN SPARC Overview and first results from CHRIS data Jose F. Moreno Depart. Thermodynamics Faculty of Physics University of Valencia, Spain Second CHRIS / PROBA Workshop ESRIN, April 2004

2 STUDY AREA vegetation map Barrax Core Site soil map Barrax Barrax

3 Major on going activities in 2003 DEMETER (DEMonstration of Earth observation Technologies in Routine irrigation advisory services) European Commission contract EVG1-CT , 12/ /2005 (Energy, Environment and Sustainable Development Programme) CHRIS / PROBA time series acquisition (AO Core Site) ASAR / MERIS / AATSR ENVISAT AO Projects ENVISAT / MERIS Biophysical Products Validation MSG / SEVIRI Biophysical Products Validation (VALERI extension) Canopy reflectance modelling studies (FLUORMOD) SPECTRA Phase-A Preparatory Studies...

4 ITAP, Albacete European Spatial Agency (ESA) University of Valencia - Remote Sensing Unit University of Valencia - Global Change Unit University of Valencia - Solar Radiation Unit University of Valencia GPDS University of Castilla-La Mancha Institute of Regional Development (IDR), Albacete Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas, Madrid Universitá degli Studi di Napoli Federico II, Italy University of Thessaly, Greece INRA-CSE, Avignon Laboratoire du Télédétection et SIRS, Tunisia Meteo-France, Toulouse, France

5 SPECTRA Reference Dataset - CHRIS / PROBA multitemporal dataset - Ground measurements: soil/vegetation/atmosphere - Complementary satellite data: scaling issues

6 SPARC Campaign CHRIS/PROBA Barrax

7 1 m / 115 bands 5 m / 126 bands 34 m / 62 bands Field measurements ROSIS HYMAP CHRIS/PROBA Scaling satellite / airborne data available in SPARC SPOT (3 July 2003) LANDSAT (15 July 2003) MODIS ASTER 300 m / 15 bands MERIS-FR MERIS-RR SEVIRI 1200 m / 15 bands

8 CHRIS/PROBA 12th July Acquisition Min Zen Ang: 20º Fly-by Time: 11:07 UT Solar Zen Ang: 22º 3598 [0º] 3599 [+36º] 359A [-36º]

9 CHRIS/PROBA 13th July Acquisition 35A0 35A1 Fly-by Zen Ang: -4º 359D Fly-by Time: 11:20 UT 359E 359F

10 CHRIS/PROBA 14th July Acquisition Min Zen Ang: -27º Fly-by Time: 11:32 UT Solar Zen Ang: 20º 35A2 [0º] 35A3 [+36º] 35A4 [-36º]

11 CHRIS/PROBA 12th July Acquisition 3598 [0º] 3599 [+36º] 359A [-36º] Min Zen Ang: 20º Fly-by Time: 11:07 UT Solar Zen Ang: 22º Pixel Resolution: 34m/pixel 359B [+55º] 359C [-55º]

12 CHRIS/PROBA 780 nm, MZA= nm, MZA= nm, MZA= nm, MZA= A5 [+55º] 35A3 [+36º] 35A2 [0º] 35A4 [-36º] 35A6 [-55º] Min Zen Ang: -27º Fly-by Time: 11:32 UT Solar Zen Ang: 20º

13 Sampling Points for Radiometric Calibration S6 S5 A S3 W S0

14 ground reflectance 1, CHRIS radiance ground radiance

15 CHRIS/PROBA DATA PROCESSING geometric processing retrieved atmospheric information angular atmospheric corrections spectral atmospheric corrections

16

17 IMU & GPS ROSIS HyMap 4 lines in total: 1-4 with ROSIS 1&4 with HyMap

18

19

20

21

22 R: Line C G: Line B B: Line A

23

24 MERIS Full Resolution coverage 14 July :34:09 ENVISAT Orbit 7162 Track 94 / Frame m

25

26 TOA reflectance Surface reflectance Geometric and atmospheric correction

27 Atmospheric characterisation Solar irradiance

28 LIDAR mobile laboratory (CIEMAT) 39.03º N, 2.07º W ( 703 m asl ) Las Tiesas TECHNICAL CHARACTERISTICS Laser source Frequency-doubled Nd:YAG (Continumm) Wavelength 532 nm Energy per pulse 20 mj Frequency 20 Hz Pulse width < 10 ns Beam diameter/divergence ~ 5 cm / ~0.2 mrad Telescope type Newtonian Telescope diameter 30 cm Telescope FOV ~ 1.3 mrad PMT Hamamatsu R928 Maximum range ~ 6 Km Spatial resolution 3-30 m

29 12/07/2003 Sounding & LIDAR profiles 11:07 UTC free troposphere mixed layer

30 Sounding & LIDAR profiles 11:32 UTC 14/07/2003 free troposphere mixed layer

31 LIDAR VERTICAL PROFILES

32 Use of standard atmospheric profiles for atmospheric corrections 2 Sounding (day/hour(utc)) Water integrated value (g/cm ) 12 th July 09: (1.52) 12 th july 11: ( 1.3 ) 13 th july 06: (1.42) 13 th july 11: (1.47) 14 th july 06: (1.94) 14 th july 10: (1.36) 1.72 (1.50) Standard atmospheres 2 Water integrated value (g/cm ) Tropical 3.25 Mid-latitude summer 2.28 Subartic summer 1.66 Mid-latitude winter 0.69 Subartic winter 0.36 US

33 Landuse Map and Sampling Points

34 CHLOROPHYLL SPECTROPHOTOMETER LABORATORY MEASUREMENTS Between 4-12 samples colected at the same time that fresh matter content 3 CCM measurements Circular sample:7.5 mm of diameter For each leaf sample NN-Dimetilforamide disolution (DFM) during hours Absorptivity for wavelengths 647 nm λ= nm measured by means of spectrophotometer ECUATIONS TO QUANTIFY: Chl a = A(664.5) 2.79 A(647) (mg/l) Chl b = A(647) 4.62 A(664.5) (mg/l) Chl total = A(647) A(664.5) (mg/l) CARY-UV-visible Spectrophotometer

35 Vegetation measurements: Leaf chlorophyll Absolute chlorophyll values were measured by means of leaf pigments extraction and absorption spectrometry. Technique: N,N-Dimethylformamide (DMF) for 20h Extinction coefficients for Chla and Chlb in and 647 (A664.5 and A647) were calculated by means of absorption measurements with an spectrometer. Finally, chlorophyll contents were determined as: Chl = 17.90A A leaf chlorophyll (mg/cm 2 ) Alfalfa Corn Sugarbeet Wheat stubble Garlic Onions Potatoes Alfalfa Corn Sugarbeet Wheat Garlic Onions Chl Chl Potatoes a b = Chl ab 0.053± ± ± ± ± ± ± = 20.70A A Chl a Chl b (mg/cm 2 ) (mg/cm 2 ) (mg/cm 2 ) 0.043± ± ± ± ± ± ± A ± ± ± ± ± ± ±0.001

36 CHLOROPHYLL IN-SITU MEASUREMENTS Chlorophyll content meter CCM-200 SUGAR-BEET GARLIC ONION B3 ON1 G1 A9 P1 ON2 C1 Transmisivity ALFALFA POTATO CORN

37 Chlorophyll (µ g/cm 2 ) Chlorophyll (µg/cm 2 ) CALIBRATION FUNCTION Digital Counts CCM Digital Counts CCM-200 m1 m2 R Alfalfa Corn Sugar Beet Wheat Garlic Onion Potato y = m1+m2*log(x) Value Error

38 Crop ESU Samples Chlorophyll Mean Value (µg cm -2 ) Crop Mean Values Corn C1-A1 C1-A2 C1-A ± 3 52 ± 4 50 ± 3 51 ± 3 B1-T ± 7 B1-T ± 6 Sugar Beet B1-T3 B1-T ± 8 43 ± 6 44 ± 7 B1-T ± 7 B1-T ± 6 On1-B ± 20 On1-A ± 20 On1-B ± 20 Onion On1-B4 On1-B ± ± ± 20 On2-A ± 20 On2-A ± 20 G1-A ± 15 Garlic G1-A8 G1-A ± ± ± 13 G1-A ± 11 P1-T ± 5 Potato P1-T11 P1-T ± 4 35 ± 5 36 ± 5 P1-T ± 5 Alfalfa A9-T ± 9 48 ± 9

39 LAI measurements 113 Elementary Sampling Units (24 data samples each ESU) b b b bbb b bb b b b b b b b b b bb b bb bbb bb b b b b b b bbbb b bb b b b b b b b b b b b b b b b b bbb b bbb bbb b b bb b b b b b b b b b b b b b bb b b b b

40 LAI 6 Alfalfa Corn Sugarbeet LAI measurements types of crop: alfalfa corn sugarbeet onions garlic potato papaver 12% 10% 8% 6% 4% 2% 0% avg = 3.07; std = LAI

41 BIOMASS Methodology applied for Biomass measurements 8 Samples 3 Leaves Dry matter content Area Fresh matter content into an oven Sampling name B3-1 1 Leaf SUGAR-BEET 8 Samples 1 Leaf 8 Samples CORN A1-1 ALFALFA C1-1 ON1 B3 G1 A9 P1 A1 C1 ON1-1 8 Samples 1 Leaf ONION Ruled paper G1-1 8 Samples 1 Leaf GARLIC

42 FRACTIONAL VEGETATION COVER Hemispherical Photographs (+) for LAI and FVC estimation ESUs description CROP FIELDS ESUs Corn C3,C2,C1,C 9 Alfalfa A9,A10,A1 15 Sugarbeet B3,B 6 Onions ON1 2 Garlic G1 3 Potato P1 6 Vineyard V 3 Papaver PA 1 Grass Grass 4 Fruit trees 2 4 Sparse Vegetation SV 2 TOTAL Mean values per crops CROP Chlorophyll mean value per crop (mg*cm -2 ) FVC mean value per crop C B, B On G P ,8 0,6 0,4 0,2 0 Fraction of Vegetation Cover (FVC) FVC_HP_A9 FVC_HP_A1 FVC_HP_A10 FVC_HP_P1 FVC_HP_C2 FVC_HP_C1 FVC_HP_C3 FVC_HP_C FVC_HP_G1 FVC_HP_B3 FVC_HP_B FVC_HP_ON FVC_HP_SV FVC_HP_VINE FVC_HP_PA FVC_HP_GRASS A9 A1 A10 P1 C2 C1 C3 C G1 B3 B ON1SV Vine GS

43 RETRIEVALS: - Small reference database - only sampled points, easy to handle - useful to define / validate algorithms - Maps of biophysical variables - maps of zenith / azimuth angles - optimal exploitation of multiangular information (CHRIS and HYMAP) - Validation - ROSIS data for spatial reference - statistical robustness of reference data

44 DATABASE FORMAT C1 corn C1-A AM

45 CLOUD SCREENING CHRIS / PROBA - 14 July 2003

46 CHRIS ENVISAT ASAR

47 12/05/03 20/05/03 12/06/03 12/07/03 14/07/03 CHRIS multitemporal series

48 MAJOR PENDING ISSUES: - Multitemporal aspects - CHRIS / PROBA as unique opportunity - data assimilation (crop growth models) - Thermal information - decoupling soil / leaves temperature - link of VIS/NIR/SWIR and TIR information

49 Multitemporal aspects: Key issue for 2004 activities

50 Optimised multitemporal acquisitions NDVI NDVI SPARC July 12th - 14th 2003 fcover Fc Fcv Wheat Corn 0.8 1) 20/4-30/4 NDVI fcover 2) 23/5-3/ ) 22/6-18/7 (most critical period) ) around 28/ ) 13/8 2/9 DOY

51 WORKPLAN FOR 2004: - Exploit multitemporal aspects in CHRIS / PROBA acquisitions (together with Landsat, MERIS and other data) - SPARC 2004 in preparation for July 2004