national institute of environmental research atmospheric physics division national institute of...

National Institute of Environmental Research

Atmospheric Physics Division National Institute of Environmental ResearchFax:82-32-568-2040, E-mail:[email protected]

Park, Il-SooPark, Il-Soo

The 6th Workshop on the Transport of Air Pollutants in Asia(Model Intercomparison Study-MICS-Asia)

February 9-10, 2004Vienna, Austria

Modeling of Long-range Transboundary Air PollutantsModeling of Long-range Transboundary Air Pollutants

In Northeast Asia during the Spring of 2002In Northeast Asia during the Spring of 2002


Emission rates for SOEmission rates for SO22(left) and NO(left) and NOxx(right)(right)

Emission data of LTP projectEmission data of LTP project


Emission rates for SOEmission rates for SO22(left) and NO(left) and NOxx(right)(right)

Made by WhiskyMade by Whisky Made by WhiskyMade by Whisky

Emission data of CGREREmission data of CGRER(Center for Global and Regional Environmental Research)(Center for Global and Regional Environmental Research)


Latitude Longitude NH3 [ton]

37 126 20,904.0037 127 30,112.7037 128 8,448.6037 129 977.836 126 16,893.3036 127 23,679.6036 128 16,190.7036 129 3,768.9035 126 15,674.5035 127 27,830.8035 128 15,831.1035 129 21,679.9034 126 15,740.0034 127 7,870.0034 128 1,803.3033 126 4,801.00

232,205.80Total

Total ammonia emission in 2001 (1°×1°)Total ammonia emission in 2001 (1°×1°)

Distribution of grid emission in 2001Distribution of grid emission in 2001

(Unit : Ton)


Latitude Longitude VOC[ton]

37 126 126,26237 127 154,35037 128 85,29037 129 10,63436 126 32,67536 127 84,22236 128 76,26036 129 17,68135 126 49,94135 127 83,80235 128 92,09735 129 128,37034 126 43,09634 127 21,54834 128 9,42433 126 7,993

1,020,645Total

Total VOC emission in 2001(1°×1°)Total VOC emission in 2001(1°×1°)

Distribution of grid emission in 2001Distribution of grid emission in 2001

National Institute of Environmental Research SO2

(2002. 3. 5 - 3. 15)

Sites

Dalian Xiamen Gangwha Geoje Rishiri Oki

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

60

NOx

(2002. 3. 5 - 3. 15)

Sites

Dalian Xiamen Ganghwa Geoje Rishiri

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

60

O3

(2002. 3. 5 - 3. 15)

Sites

Dalian Xiamen Gangwha Geoje Rishiri Oki

Co

nce

ntr

atio

n(p

pb)

0

20

40

60

80

100

Comparison with measurement stations for 5-15 March 2002.Comparison with measurement stations for 5-15 March 2002.

National Institute National Institute of Environmental Researchof Environmental Research

National Institute of Environmental ResearchSO2

0.0

0.5

1.0

1.5

2.0

2.5

3.0

3.5

4.0

4.5

02/3/5 02/3/6 02/3/7 02/3/8 02/3/9 02/3/10 02/3/11 02/3/12 02/3/13 02/3/14 02/3/15

Days

Log(c

on.*

100

)

Dalian Xiamen GanghwaGeoje Rishiri Oki

NO2

0.0

0.5

1.0

1.5

2.0

2.5

3.0

3.5

02/3/5 02/3/6 02/3/7 02/3/8 02/3/9 02/3/10 02/3/11 02/3/12 02/3/13 02/3/14 02/3/15

Days

Log(C

on.*

10)

Dalian Xiamen GanghwaGeoje Rishiri

O3

0.0

10.0

20.0

30.0

40.0

50.0

60.0

70.0

80.0

02/3/5 02/3/6 02/3/7 02/3/8 02/3/9 02/3/10 02/3/11 02/3/12 02/3/13 02/3/14 02/3/15

Days

Conc.(

ppb)

Dalian Xiamen

Ganghwa Geoje

Rishiri Oki

The variations of SOThe variations of SO22. NO. NO22, O, O33 concentration for 5 ~ 15 March 2002. concentration for 5 ~ 15 March 2002.



Dalia

n

Gangw

ha

Taea

n

Gosa

n

Geo

je

Ris

hiri

Oki

SO

42-

NO

3- Cl-

NH

4+ K+

Ca2+

Mg2+

Na+ F-

EC

arb

on

OC

arb

on

Mass

Conc.

0.00

10.00

20.00

30.00

40.00

50.00

60.00

70.00

Conc.

(μg

/m3)

Dalia

nG

angw

haTaea

n

Gosa

n

Geo

je

Ris

hiri

Oki

SO

42-

NO

3- Cl-

NH

4+ K+

Ca2+

Mg2+

Na+ F-

EC

arb

on

OC

arb

on

Mass

Conc.

0.00

20.00

40.00

60.00

80.00

100.00

Conc.

(μg

/m3)

Ratio of ion composition Ratio of ion composition in PMin PM2.52.5 Mass(2002.3.5- Mass(2002.3.5-

3.15)3.15)

Ratio of ion composition Ratio of ion composition in PMin PM1010 Mass(2002.3.5-3.15) Mass(2002.3.5-3.15)


CADMCADM (Comprehensive Acid (Comprehensive Acid

Deposition Model)Deposition Model)


Surface weather map for 00UTC 7, 8, 9 and 11 March 2002 Surface weather map for 00UTC 7, 8, 9 and 11 March 2002


Weather map Weather map of 850hPa level of 850hPa level forfor 00UTC 7, 8, 9 and 11 March 2002 00UTC 7, 8, 9 and 11 March 2002

National Institute of Environmental ResearchSimulated sea-level Simulated sea-level

pressure pressure

and wind fields and wind fields

for 9 March 2002.for 9 March 2002.

National Institute of Environmental ResearchSimulated geopotential Simulated geopotential height and wind fields height and wind fields

of 850hPa level of 850hPa level for 9 March 2002.for 9 March 2002.

National Institute of Environmental Research The variations of The variations of

SOSO22 concentration concentration

on 9 March 2002.on 9 March 2002.


Concentrations in the x-z cross section(y= 34) on 9 March 2002. Concentrations in the x-z cross section(y= 34) on 9 March 2002.


Concentrations in the y-z cross section(x= 20) on 9 March 2002. Concentrations in the y-z cross section(x= 20) on 9 March 2002.

National Institute of Environmental ResearchA: SeoulB: TatenoC: Tangshan

A: SeoulB: TatenoC: Tangshan



Vertical distribution of SOVertical distribution of SO22 concentration on 9 March 2002. concentration on 9 March 2002.


SOSO22 dry deposition on 9 March 2002. dry deposition on 9 March 2002.


SOSO22 wet deposition on 9 March 2002. wet deposition on 9 March 2002.


Flux (ton/hr)

0 50 100 150 200 250 300

Ho

ur

(Ma

rch

9,

20

02

)

3

6

9

12

15

18

21

24

NOx

Flux (ton/hr)

0 10 20 30 40 50 60 70

Ho

ur

(Ma

rch

9,

20

02

)

3

6

9

12

15

18

21

24

O3

Flux (ton/hr)

0 200 400 600 800 1000 1200 1400 1600

Ho

ur

(Ma

rch

9,

20

02

)

3

6

9

12

15

18

21

24

SO2

Flux (ton/hr)

0 100 200 300 400

Ho

ur

(Ma

rch

9,

20

02

)

3

6

9

12

15

18

21

24

NOx

Flux (ton/hr)

0 2 4 6 8 10 12

Ho

ur

(Ma

rch

9,

200

2)

3

6

9

12

15

18

21

24

O3

Flux (ton/hr)

0 200 400 600 800 1000 1200 1400

Ho

ur

(Ma

rch

9,

20

02

)

3

6

9

12

15

18

21

24

The variations of The variations of in/out fluxin/out flux on 9 March 2002. on 9 March 2002.

124o 130o

32o

38o

Dry depositio

n(ug/m2)

Suspended concentration

ug/m3 ppb

SO2 0.01439934.4096

41.20434E-05

NOx 0.000956 7.74908 4.15546E-06

O3 0.257263169.792

07.92363E-05

National Institute of Environmental Research SO2

(2002. 3. 9 ~ 3. 10)

Sites

Dalian Xiamen Gangwha Gosan Rishiri Oki

Co

nce

ntr

atio

n(p

pb

)

0

50

100

150

200

250ObservationPrediction

NOx

(2002. 3. 9 ~ 3. 10)

Sites

Dalian Xiamen Ganghwa Gosan Rishiri

Co

nce

ntr

atio

n(p

pb

)

0

20

40

60


O3

(2002. 3. 9 ~ 3. 10)

Sites

Dalian Xiamen Gangwha Gosan Rishiri Oki

Co

nce

ntr

atio

n(p

pb

)

0

20

40

60

80


Comparison with monitoring data for 9-10 March 2002. Comparison with monitoring data for 9-10 March 2002.

Model results are represented by dots, while measurements are presented as whiskersModel results are represented by dots, while measurements are presented as whiskers


National Institute of Environmental ResearchDalian

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

50

100

150

200

250

ObservationPrediction

9/March/02 10/March/02

Gangwha

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

10

20

30

40

50



Gosan

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)0

1

2

3

4

5

6

7



Rishiri

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

1

2

3

4

5



Oki

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

2

4

6

8

10



Comparison with monitoring data and LTP model results Comparison with monitoring data and LTP model results of SOof SO22 concentration for 9~10 March 2002. concentration for 9~10 March 2002.



Comparison with monitoring data and LTP model results Comparison with monitoring data and LTP model results of NOof NOxx concentration for 9~10 March 2002. concentration for 9~10 March 2002.

Dalian

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

60

70

80



Gosan

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Co

nce

ntr

atio

n(p

pb

)

0

1

2

3

4

5



Rishiri

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Co

nce

ntr

atio

n(p

pb

)

0

1

2

3

4

5



Gangwha

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Co

nce

ntr

atio

n(p

pb

)

0

20

40

60

80

100

120

140


National Institute of Environmental ResearchDalian

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

20

40

60

80

100



Gangwha

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

20

40

60

80



Gosan

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

20

40

60

80

100

120



Rishiri

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

10

20

30

40

50

60



Oki

Days

3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48

Con

cent

ratio

n(pp

b)

0

20

40

60

80

100



Comparison with monitoring data and LTP model results Comparison with monitoring data and LTP model results of Oof O33 concentration for 9~10 March 2002. concentration for 9~10 March 2002.


National Institute of Environmental ResearchAircraft ObservationAircraft ObservationAircraft ObservationAircraft Observation

Variations of SOVariations of SO22, NOx, O, NOx, O33 measured in 7, 8 March 2002 aircraft flight. measured in 7, 8 March 2002 aircraft flight.


SO2

0

1

2

3

4

5

6

7

0 500 1000 1500 2000 2500 3000 3500

Altitude(m)

ppb

3/7 PM 3/8 PM

3/9 PM 3/11 PM

NOx

0

2

4

6

8

10

12

0 500 1000 1500 2000 2500 3000 3500

Altitude(m)

ppb

3/7 PM 3/8 PM

3/9 PM 3/11 PM

O3

0

10

20

30

40

50

60

70

0 500 1000 1500 2000 2500 3000 3500

Altitude(m)

ppb

3/7 PM 3/8 PM

3/9 PM 3/11 PM

Distribution of SODistribution of SO22, NO, NOxx and O and O33 with height for 7, 8, 9 and 11 March 2002. with height for 7, 8, 9 and 11 March 2002. Distribution of SODistribution of SO22, NO, NOxx and O and O33 with height for 7, 8, 9 and 11 March 2002. with height for 7, 8, 9 and 11 March 2002.



SO2

(2002. 3. 7., 124.30o, 37.10og 35.08o)

Height(m)

0 500 1000 1500 2000 2500 3000 3500

Concentr

ation(p

pb)

0

1

2

3

4

AircraftLTPMICS

SO2

(2002. 3. 7., 124.30o, 35.11og 37o)

Height(m)

0 500 1000 1500 2000 2500 3000 3500

Concentr

ation(p

pb)

0

1

2

3

4

AircraftLTPMICS

SO2

(2002. 3. 8., 124.30o, 36.56og 34.59o)

Height(m)

0 500 1000 1500 2000 2500 3000

Co

nce

ntr

atio

n(p

pb

)

0

1

2

3

AircraftLTPMICS

SO2

(2002. 3. 8., 124.30o, 35.07og 36.45o)

Height(m)

200 400 600 800 1000 1200 1400 1600 1800 2000 2200

Co

nce

ntr

atio

n(p

pb

)

0

1

2

3

4

AircraftLTPMICS

SO2

(2002. 3. 9.,125.30o~126.34o, 36.48o~36.11o)

Height(m)

550 600 650 700 750 800 850 900 950

Con

cen

tratio

n(p

pb)

0

1

2

3

4

AircraftLTPMICS

SO2

(2002. 3. 11., 124.40o, 37.14o~36.46o)

Height(m)

0 500 1000 1500 2000 2500 3000

Con

centr

ation(p

pb)

0

1

2

3

4

5

6

7

AircraftLTPMICS

Comparison with Comparison with

aircraft observation aircraft observation

and model results and model results

of LTP and MICS of LTP and MICS

for 7, 8, 9 and for 7, 8, 9 and

11 March 2002.11 March 2002.

Comparison with Comparison with

aircraft observation aircraft observation

and model results and model results

of LTP and MICS of LTP and MICS

for 7, 8, 9 and for 7, 8, 9 and

11 March 2002.11 March 2002.

SOSO22 Concentration ConcentrationSOSO22 Concentration Concentration


NOx

(2002. 3. 7., 124.30o, 37.10og 35.08o)

Height(m)

0 500 1000 1500 2000 2500 3000 3500

Co

nce

ntr

atio

n(p

pb

)

0

1

2

3

4

AircraftLTPMICS

NOx

(2002. 3. 7., 124.30o, 35.11og 37o)

Height(m)

0 500 1000 1500 2000 2500 3000 3500

Co

nce

ntr

atio

n(p

pb)

0

1

2

3

4

AircraftLTPMICS

NOx

(2002. 3. 8., 124.30o, 36.56og 34.59o)

Height(m)

0 500 1000 1500 2000 2500 3000

Concentr

atio

n(p

pb

)

0

2

4

6

8

10

AircraftLTPMICS

NOx

(2002. 3. 8., 124.30o, 35.07og 36.45o)

Height(m)

200 400 600 800 1000 1200 1400 1600 1800 2000 2200

Co

nce

ntr

atio

n(p

pb

)

0

2

4

6

8

10

AircraftLTPMICS

NOx

(2002. 3. 9.,125.30o~126.34o, 36.48o~36.11o)

Height(m)

550 600 650 700 750 800 850 900 950

Concentr

ation(p

pb)

0

1

2

3

4

5

AircraftLTPMICS

NOx

(2002. 3. 11., 124.40o, 37.14o~36.46o)

Height(m)

0 500 1000 1500 2000 2500 3000

Concentr

ation(p

pb)

0

1

2

3

4

AircraftLTPMICS

NONOxx Concentration ConcentrationNONOxx Concentration Concentration


O3

(2002. 3. 7., 124.30o, 37.10og 35.08o)

Height(m)

0 500 1000 1500 2000 2500 3000 3500

Co

nce

ntr

atio

n(p

pb

)

0

20

40

60

80

AircraftLTPMICS

O3

(2002. 3. 7., 124.30o, 35.11og 37o)

Height(m)

0 500 1000 1500 2000 2500 3000 3500

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

AircraftLTPMICS

O3

(2002. 3. 8., 124.30o, 36.56og 34.59o)

Height(m)

0 500 1000 1500 2000 2500 3000

Co

nce

ntr

atio

n(p

pb

)

0

20

40

60

80

AircraftLTPMICS

O3

(2002. 3. 8., 124.30o, 35.07og 36.45o)

Height(m)

200 400 600 800 1000 1200 1400 1600 1800 2000 2200

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

AircraftLTPMICS

O3

(2002. 3. 9.,125.30o~126.34o, 36.48o~36.11o)

Height(m)

550 600 650 700 750 800 850 900 950

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

60

70

AircraftLTPMICS

O3

(2002. 3. 11., 124.40o, 37.14o~36.46o)

Height(m)

0 500 1000 1500 2000 2500 3000

Co

nce

ntr

atio

n(p

pb

)

0

10

20

30

40

50

60

AircraftLTPMICS

OO33 Concentration ConcentrationOO33 Concentration Concentration


Days

02/3/7 02/3/8 02/3/9 02/3/11

Co

rre

latio

n c

oe

ffic

ien

ts(R

)

-1.0

-0.5

0.0

0.5

1.0

LTPMICS

NOx

Days

02/3/7 02/3/8 02/3/9 02/3/11

Co

rre

latio

n c

oe

ffic

ien

ts(R

)

-1.0

-0.5

0.0

0.5

1.0

LTPMICS

O3

Days

02/3/7 02/3/8 02/3/9 02/3/11

Co

rre

latio

n c

oe

ffic

ien

ts(R

)

-1.0

-0.5

0.0

0.5

1.0

LTPMICS

Correlation between aircraft observed and predicted data Correlation between aircraft observed and predicted data of LTP and MICSof LTP and MICS for 7, 8, 9 and 11 March 2002.for 7, 8, 9 and 11 March 2002. Correlation between aircraft observed and predicted data Correlation between aircraft observed and predicted data of LTP and MICSof LTP and MICS for 7, 8, 9 and 11 March 2002.for 7, 8, 9 and 11 March 2002.



Emissions

Duration of modeling

Source-Receptor

Work PlanWork Plan

PM10& CO in 1°×1° In China, Emission will depend on funding, but Ammonia and VOC emission data will be reported to the Secretariat by the next meeting

March 1 to 31, 2002 July 1 to 31, 2002

Sources 1) N.E. China 2) C.E. China 3) S. China 4) S. Korea 5) Japan Receptors 1) The same as sources


Spatial distribution of concentration (SO2, SO42-, NO2, O3)

Half-month average surface concentrations(first grid level above the surface) for 00, 06, 12, 18UTC and 24hours

Surface concentrations and wet deposition(SO2, SO42-, NO2, O3) at 6 stations

(Dalian, Xiamen, Gosan, Kangwha, Oki, Rishiri)

Vertical profiles at inland stations(Tateno in Japan, Seoul in Korea, Tangshan in China) to see vertical diffusion(SO2, NO2, O3) at 00, 06, 12, 18UTC for

10th, 20th & 30th days of each month

Model intercomparisons

Daily average


Dry and wet deposition spatial distribution

Source - receptor relationship

Components: SO2, SO42-, NOx, ,NO3

-, for dry deposition

Monthly deposition amount

※ Reference for Identical Condition for Model Intercomparison

Components : sulfur

Analyzing

Contribution ratio and deposition by source


40, 114

45, 114

52, 132

42, 132

30, 114 30, 123

21, 114 21, 119

25, 122

SO2 : 6,206,755 ton/yr

NOx : 2,968,761 ton/yr

SO2 : 1,603,167 ton/yr

NOx : 678,800 ton/yr

Source - Receptor L

atit

ud

e

Longitude

42, 124

SO2 : 1,618,422 ton/yr

NOx : 1,192,899 ton/yr

45, 120

40, 124

40, 123

52, 120

25, 119

¥ °

¥ ±

¥ ²

38, 126 38, 130

33, 126 33, 130

SO2 : 1,146,000 ton/yrNOx : 1,083,700 ton/yr

¥ ³

¥ ´SO2 : 916,782 ton/yr

NOx : 1,992,535 ton/yr

115 120 125 130 135 140 14520

25

30

35

40

45

50

55

(%)100

1

n

iij

H

ijH

ijR

Rij: Contribution rate to jth

receptor by ith source

Hij : Concentration or

deposition of jth receptor by only ith source

Where,

national institute of environmental research atmospheric physics division national institute of...

Documents