Konstantinos KourtidisKonstantinos Kourtidis

Aristeidis GeorgouliasGeorgoulias

Antonios GasteratosAntonios Gasteratos

Konstantinos KonstaninidisKonstantinos Konstaninidis

AMFIC Web Data BaseAMFIC Web Data Base

AMFIC Final Meeting - AMFIC Final Meeting - Beijing 23 October 2009Beijing 23 October 2009

Democritus University of ThraceDemocritus University of ThraceLaboratory of Atmospheric Pollution and Pollution Laboratory of Atmospheric Pollution and Pollution Control Engineering of Atmospheric PollutantsControl Engineering of Atmospheric Pollutants

下午好下午好good afternoongood afternoon

1.1. Current status Current status::• Data brought in common gridded format and Data brought in common gridded format and

inserted in the databaseinserted in the database

• AMFIC web data base is ready and functional, AMFIC web data base is ready and functional, http://www.amfic.satellite-earth-simulator.com/

Fully available productsFully available products::SCIAMACHY WFM-DOAS v1.0 XCHSCIAMACHY WFM-DOAS v1.0 XCH44 dry air column dry air column

SCIAMACHY WFM-DOAS v0.6 CO column SCIAMACHY WFM-DOAS v0.6 CO column

SCIAMACHY SO2 column SCIAMACHY SO2 column

SCIAMACHY NO2 tropospheric and total column SCIAMACHY NO2 tropospheric and total column

The example of Methane dataThe example of Methane data::

The example of Methane dataThe example of Methane data::

FRONTFRONT

BACKBACK

Methane Methane productproduct

The example of Methane dataThe example of Methane data::

BACKBACK

The example of Methane dataThe example of Methane data::

The example of Methane dataThe example of Methane data::

The example of Methane dataThe example of Methane data::

The example of Methane dataThe example of Methane data::

The example of Methane dataThe example of Methane data::

XianXian

The example of Methane dataThe example of Methane data::

Some advantages of AMFIC web data baseSome advantages of AMFIC web data base::

• By breaking huge global files to many gridded ascii files we By breaking huge global files to many gridded ascii files we make the process of those data easier for users interested in make the process of those data easier for users interested in specific spots. specific spots.

• Format and accessibility enables easy validation/comparison of Format and accessibility enables easy validation/comparison of several products with GND-based, modelled and other data for several products with GND-based, modelled and other data for selected regionsselected regions

• The users can request either plots or ascii files for a quick look The users can request either plots or ascii files for a quick look in secondsin seconds

• Users interested in global data sets can just download the whole Users interested in global data sets can just download the whole dataset using Wget or a web ripper software dataset using Wget or a web ripper software

• Suitable for educational purposes Suitable for educational purposes

A case study on the use of Methane dataA case study on the use of Methane data3. Analysis3. Analysis

**The XCHThe XCH44 SCIAMACHY data SCIAMACHY data appearing on the 4 figures appearing on the 4 figures have been scaled by 1.02 as have been scaled by 1.02 as discussed in Schneising et al., discussed in Schneising et al., 20092009

1 2 3 4 5 6 7 8 9 10 11 121650

1700

1750

1800

1850

1900

1950 SCIAMACHY AIRS 206.25mb AIRS 306.75mb AIRS 407.25mb Flask-NOAA

Met

han

e vm

r (p

pb

)

Met

han

e vm

r (p

pb

)

Month1 2 3 4 5 6 7 8 9 10 11 12

1650

1700

1750

1800

1850

1900

1950

UUM (914m)WLG (3812m)

Month

1 2 3 4 5 6 7 8 9 10 11 121650

1700

1750

1800

1850

1900

1950

Met

han

e vm

r (p

pb

)

KZD (412m)

Met

han

e vm

r (p

pb

)

Month1 2 3 4 5 6 7 8 9 10 11 12

1650

1700

1750

1800

1850

1900

1950

KZM (2519m)

Month

A case study on the use of Methane dataA case study on the use of Methane data3. Analysis3. Analysis

1.The good agreement between the satellite 1.The good agreement between the satellite and ground observations is probably due to and ground observations is probably due to the high altitude of the station.the high altitude of the station.

2. According to Xiong et al., 2009 the 2. According to Xiong et al., 2009 the (August-September) AIRS peak is connected (August-September) AIRS peak is connected to methane transport to the upper to methane transport to the upper troposphere from the monsoon-Tibetan troposphere from the monsoon-Tibetan anticyclone system.anticyclone system.

3. According to Xiong et al., 2009 the decrease of AIRS vmr thereafter is connected to 3. According to Xiong et al., 2009 the decrease of AIRS vmr thereafter is connected to the withdrawal of the monsoon and the dissipation of the Tibetan anticyclone.the withdrawal of the monsoon and the dissipation of the Tibetan anticyclone.

4.The good agreement between SCIAMACHY and AIRS observations could be 4.The good agreement between SCIAMACHY and AIRS observations could be indicative of the transport of methane to the upper troposphere.indicative of the transport of methane to the upper troposphere.

A case study on the use of Methane dataA case study on the use of Methane data3. Analysis3. Analysis

The good agreement between The good agreement between SCIAMACHY and AIRS observations SCIAMACHY and AIRS observations could be indicative of the transport of could be indicative of the transport of methane to the upper troposphere in the methane to the upper troposphere in the region around the Tibetan Plateau.region around the Tibetan Plateau.

DiscussionDiscussion::

• No comments received so far from the partners providing data about No comments received so far from the partners providing data about the terms of use/copyright/ availability of full datasets to the users.the terms of use/copyright/ availability of full datasets to the users.

• We hence assume that partners agree with what is currently stated in We hence assume that partners agree with what is currently stated in the web data basethe web data base

• HCHO data not available to us in a form dowloadable with HCHO data not available to us in a form dowloadable with reasonable effort, hence not included (each day in separate folders, reasonable effort, hence not included (each day in separate folders, hence not downloadable with DownThemAll or WGet scripts; hence not downloadable with DownThemAll or WGet scripts; currently not available in specified folders). currently not available in specified folders).

3.3. Automatic Plume Detection Automatic Plume Detection::Within the AMFIC data base demonstration of image analysis Within the AMFIC data base demonstration of image analysis techniques for automatic plume detectiontechniques for automatic plume detection

Three-step procedure:Three-step procedure:

1.1. Reconstruct the missing parts of maps (images), which lack Reconstruct the missing parts of maps (images), which lack information due to the daily coverage stripes of several satellite information due to the daily coverage stripes of several satellite instruments or due to cloudiness, albedo characteristics etc, instruments or due to cloudiness, albedo characteristics etc, with the use of with the use of Cellular Automata (CA)Cellular Automata (CA)..

2.2. Extract only the regions on the map which are covered by Extract only the regions on the map which are covered by significant plumes with the use of a cut-off filter.significant plumes with the use of a cut-off filter.

3.3. Calculation of the area covered from each plume.Calculation of the area covered from each plume.

• Matlab image analysis algorithm runs off-line, integration not Matlab image analysis algorithm runs off-line, integration not feasiblefeasible

Cellular Automaton approachCellular Automaton approach::

A A cellular automaton requires:cellular automaton requires:1.1. a regular lattice of cells covering a portion of a four dims space a regular lattice of cells covering a portion of a four dims space

2.2. a set of variables a set of variables CC attached to each cell giving its local state at the time attached to each cell giving its local state at the time t=0, 1, 2, …t=0, 1, 2, …

3.3. a rule a rule R={RR={R11,R,R22,…,R,…,Rmm}} which specifies the time evolution of the which specifies the time evolution of the

states in the following way: states in the following way:

where designate the cells belonging to a given neighbor- hood of cell where designate the cells belonging to a given neighbor- hood of cell

1 2( , ) ( , ), ( , ),..., ( , )mC r t C r t C r t C r t

r

( , )C r t

1 2( , 1) ( ( , ), ( , ), ( , ),..., ( , ))j j qC r t R C r t C r t C r t C r t

r

r

In our case each cell is represented by an image pixelIn our case each cell is represented by an image pixel

The rule R applied here is the extraction of the The rule R applied here is the extraction of the average from pixels in the neighborhood (Moore average from pixels in the neighborhood (Moore neighborhood ) that contain informationneighborhood ) that contain information

An example for GOME-2An example for GOME-2::GOME-2 on MetOp-A (October 2006)GOME-2 on MetOp-A (October 2006)

• 44 channels cover the full spectral range from 0.240 to 0.790 µm channels cover the full spectral range from 0.240 to 0.790 µm • Resolution 0.2-0.4 nmResolution 0.2-0.4 nm• Pixel size 80x40 kmPixel size 80x40 km22

• Scan width 1920 kmScan width 1920 km• Global coverage within 3 dayGlobal coverage within 3 day

Input Input NO2 daily mapsTEMIS website

Process 1Process 1Reconstruct the

missing parts

Process 2Process 2Define areas withsignificant plumes

Process 3Process 3Calculate the area each plume covers

INPUT: INPUT: NONO22 daily maps from TEMIS web data base daily maps from TEMIS web data base www.temis.nlTime period 28/3/2008-6/4/2008Time period 28/3/2008-6/4/2008

28/3 29/3 30/3 31/3

1/4 2/4 3/4 4/4

5/4 6/4

PROCESS 1:PROCESS 1: Reconstruct full maps (images) Reconstruct full maps (images)Time period 28/3/2008-6/4/2008Time period 28/3/2008-6/4/2008

28/3 29/3 30/3 31/3

1/4 2/4 3/4 4/4

5/4 6/4

PROCESS 2:PROCESS 2: Define areas with significant plumes using color Define areas with significant plumes using color filters, time period 28/3/2008-6/4/2008filters, time period 28/3/2008-6/4/2008

28/3 29/3 30/3 31/3

1/4 2/4 3/4 4/4

5/4 6/4

PROCESS 3:PROCESS 3: Calculate the areas covered from significant plumes Calculate the areas covered from significant plumes Time period 28/3/2008-6/4/2008Time period 28/3/2008-6/4/2008

28/3 29/3 30/3 31/3

1/4 2/4 3/4 4/4

5/4 6/4

If we have larger stripes (e.g. SCIAMACHY daily maps) … If we have larger stripes (e.g. SCIAMACHY daily maps) … Example for 28/3/2008Example for 28/3/2008

vsvs vsvs

3.3. Work status Work statusObjectives

The objective of this WP is to exploit the outcome of the WP2, WP3, WP4 and WP5, to determine the applicability of satellite observations for studying and monitoring transport of air pollutants. It is to demonstrate how to use data from satellites to explore processes relevant to the transport of air pollution and to monitor these processes. It is also to integrate this ability with the GEOSS, INSPIRE and GMES requirements. Finally, to provide an extensive system assessment report, suitable for feeding GEOSS and GMES requirements, as well as an assessment of options and operational requirements for long-term operation of such a satellite system.

For the achievement of these goals, data will be homogenized and made available in a web database, and a workshop will be organized where all potential end-users will be invited and will interact with the project partners.Accomplished

3.3. Work status Work status

Description of work Task 6.1: The different data fields, observed, modelled and/or

assimilated (NO2, ozone, SO2, HCHO, CO, CH4, aerosol) will be homogenised and brought into a common spatial and temporal format, complying with INSPIRE requirements. At the project website, a web database will host these data. Underlying the database structure, a set of tools will allow easy manipulation and analysis of the data, more specifically:

Automatic detection of features such as plumes and monitoring of their evolution, Cluster analysis

Data co-variance analysis, enabling the detection of emission signatures of different sources

Export of data in different formats Accomplished

Aristeidis K. GeorgouliasAristeidis K. Georgoulias

Konstantinos KourtidisKonstantinos Kourtidis

Antonios GasteratosAntonios Gasteratos

Konstantinos KonstaninidisKonstantinos Konstaninidis

AMFIC Web Data BaseAMFIC Web Data Base

AMFIC Final Meeting - AMFIC Final Meeting - Beijing 23 October 2009Beijing 23 October 2009

Democritus University of ThraceDemocritus University of ThraceLaboratory of Atmospheric Pollution and Pollution Laboratory of Atmospheric Pollution and Pollution Control Engineering of Atmospheric PollutantsControl Engineering of Atmospheric Pollutants

谢谢谢谢Thank youThank you

THERMOPOLIS2009 ESA campaign, Athens, Greece, summer 2009 Study of urban heat island Use of satellite data Thermal Stress Indices Forecasting

THERMOPOLIS2009 Ground Stations.

ASTER imagery acquired during over Athens on July 16th (09:22 UTC), during the THERMOPOLIS campaign. This is a false colour composite image produced from the raw visible and near infrared channels of the sensor (R3G2B1). The sensor records information in 15 spectral channels whereas its spatial resolution varies from 15 m (visible to near infrared channels) to 30 m (shortwave infrared channels) and 90 m (thermal infrared channels). National Observatory of Athens.

LANDSAT TM imagery acquired on July 24th (08:53 UTC) over the Athens during the THERMOPOLIS campaign. The sensor records information in 7 spectral bands and it has a spatial resolution 30 m in the reflective bands (visible to shortwave infrared) and 120 m in the thermal infrared band. National Observatory of Athens.

Land Surface Temperature (LST) map for July 18th, one of the days of the THERMOPOLIS expedition, computed from the METEOSAT-2 SEVIRI radiometer. National Observatory of Athens.

Map with daily TERRA+AQUA (FLASHFLUX v2F) Shortwave TOA Flux (W/m2) (a), and Longwave TOA Flux (W/m2) (b) data for 17/7/2009. DUTH.

CALIPSO Lidar browse images from expedited release dataset (v2.02). Attenuated backscatter coefficient at 532nm profile google earth picture for 17/7/2009. Overpass time ~11:45UTC. DUTH.

Maps with Level 3 MODIS (C005) AOD550 data (a) and Ångström exponent (b) data for 17/7/2009. DUTH.

DUTH Sonic anemometer on 9-m tower for sensible/latent heat fluxes measurement.

Selected spectra on 17 of July 2009 at 08:00, 09:00, 10:00 and 12:00 LT, measured with the PGS-100 spectrophotometer. DUTH/NOA.

Athens transect, 21 July 2009. UVEG.

0

5

10

15

20

25

30

35

40

0:00 0:30 1:00 1:30 2:00 2:30

Time (UTC)

Te

mp

era

ture

(ºC

)

0

50

100

150

200

250

300

350

400

He

igh

t (m

)

T air

Height

Lake

UVEG

Surface and Air temperature in Bare Soil

0

10

20

30

40

50

14-7 15-7 16-7 17-7 18-7 19-7 20-7 21-7 22-7 23-7 24-7 25-7 26-7

Date (UTC)

Te

mp

era

ture

(ºC

)

Air T (ºC)

Bare Soil Ts (ºC)

Surface and Air temperature in Green Grass

0

10

20

30

40

50

14-7 15-7 16-7 17-7 18-7 19-7 20-7 21-7 22-7 23-7

Date (UTC)

Te

mp

era

ture

(ºC

)

Air T (ºC)

Green Grass Ts (ºC)

Thermography. UVEG.

Land surface temperature map (left panel) and urban heat island map (right panel) for Athens heat wave.

AMFIC continuation LST maps for Chinese cities LST + AirT Thermal Stress Indices Maps of Cooling Degree Days for energy

consumption TOA SW/LW TOA SW/LW + Aerosols for radiative

forcing of aerosol pollution