Bulletin of JESC No.43. 201662

大気移流拡散モデルの改良とその適用
Improvement and Application of an Air Advection Diffusion Model
池本久利 1，植田洋匡 2
Hisatoshi IKEMOTO 1 and Hiromasa UEDA 2
近年、計算機の性能向上やコミュニティシステムが普及してきたこと等により、モデル計算に関する取り組みが進んできている。本研究では、近年のモデル計算に関する状況を整理し、使用環境を整えるとともに、水銀(Hg)の全球規模の挙動に関するモデル計算の精度向上への取組み、及び冬季北京のPM 2.5 汚染を対象として大気汚染と都市気象変化の相乗作用の研究を行った。その結果、Hgの予測精度は臭素(Br)の反応を組み込むことで精度が向上すること、大気汚染と都市気象変化はブラックカーボン(BC)による日射の吸収、地表気温の低下、逆転層の強度と高度の変化を促し、さらなる高濃度の汚染状況を作り出していることが導き出された。

キーワード：air advection diffusion model, chemical transport model, mercury, PM 2.5

１．研究目的きた。また、気象関係従事者による天気予報や、近年、計算機の性能向上（並列処理などの高速一般の研究者による研究活動等に利用されており、化、大容量化）により、PCクラスタマシンを用いMM5 （ Fifth-Generation Mesoscale Model ） はたモデル計算が可能になり、環境関連分野で広く

Pennsylvania State UniversityとNational Center for

用いられはじめた。このことを受け、本研究の目

Atmospheric Research（NCAR）が開発したもので的は、モデルの種類やその現在の状況を整理し、韓国などでは現業業務、天気予報にも使用されて水銀(Hg)等の化学物質への適用を検討することにい る 。 そ の 後 継 モ デ ル 、 WRF-ARW (Weatherある。

Research and Forecasting, Advanced Research WRF )大気汚染移流拡散の予測システムは、気象モデ

モデルはNCARが中心になって開発したもので、

ル、化学輸送モデル、エミッションインベントリ、コミュニティもNCAR中心で運営され、システム

グ ラ フ ィ ック ス モ デ ルか ら な る フレ ー ム ワ ーク

も進化して、現在はバージョン3.8が公開されてい（枠組み）である。各モデルは物理的、化学的なる。

素過程（プロセス）をパラメタライズした種々の化学輸送予測システムとしては、US-EPAが開発モジュールから構成されている。例えば、化学輸

し た CMAQ （ Community Multi-scale Air Quality

送モデルには、大気拡散、化学反応、乾性、湿性Model）モデルやNCARのWRF-ARWモデルに対応沈着などの素過程が含まれ、化学反応過程には光す る WRFchem モ デ ル 、 ハ ー バ ー ド 大 学 の化学反応、エアロゾル動力学等のモジュールが含GEOS-chemなどがある。

まれる。

我々は、これまで大気乱流と乱流拡散過程、光気象予測システムの構築は、先進国を中心に各化学反応、エアロゾル動力学等の素過程の理論的、国気象庁で進められてきた。1980年代後半には、実験的研究を行い、これらのモデリング、モジュ欧米を中心に「共用」を目的とした、コミュニテール作成を通して、MM5、WRF-ARW、WRFchemィシステム（コミュニティモデルともいう）が開のシステム進化に貢献してきた。また、これらを発され、計算プログラム･コード、理論・内容説明用いて、海陸風、山岳気象、台風、集中豪雨等の書、ユーザーズガイドを公開し、コミュニティにメソ異常気象や光化学大気汚染、酸性雨や越境大参加する研究者によるシステム進化が進められて気 汚 染 研 究 を 行 っ て き た （ Uno,I., Ueda,H. and1 一般財団法人日本環境衛生センター東日本支局環境工学部計画建設技術課 Japan Environmental Sanitation
Center Environmental Engineering Dept. Municipal Solid Waste Consulting Div.
2 一般財団法人日本環境衛生センター顧問 Japan Environmental Sanitation Center adviser

日 環 セ 所 報 No.43. 平成28年

63

Wakamatsu,S. 1989: Boundary-Layer Meteorology,Data, National Centers for Environmental Prediction）

49, 77-98など）。のものを用いる。これは、6時間/ 日、水平最大解

本研究では、（１）これらのモデルシステムを像度 1 °（約100 km）、鉛直26層（最高高度10 hPa

日本環境衛生センターのPCクラスタマシンにイ高）の気象要素がリストアップされており、1960ンストールして、種々のユーザ、目的に対応でき年以降、今日までのものがインターネットからフる計算環境を構築する。（２）さらに、これら予リーでダウンロードできる。

測システムを用いて、１）Hgの全球的な挙動、２）気象モデルは、大気中での運動量、熱、水分のPM2.5 の大気汚染と都市気象変化の相乗作用汚染、輸送、地表面、水面での交換過程、大気放射、積それに伴う都市気象変化と気象変化に伴う大気汚

雲対流、雲物理過程等を計算して、風向風速、気染の更なる悪化についての研究を進める。

温、湿度、雲、雨、雪等の分布を高い時間及び空間分解能で予測する。各過程には、理論及び実験

２．大気移流拡散予測システムの概要
大気移流拡散予測システムの構成は Fig.1 のとおりである。に基づくモジュールが数種類用意されており、ユーザは目的に合わせて選択できる。本研究では、気象モデルとして、MM5とその後継、WRF-ARWをインストールして実行する。
2.2

化学輸送モデル気象場での化学物質の輸送、拡散、反応、沈着過程を計算する。拡散過程は気象モデルでも共通に使用する。拡散過程には、単純なものから高次の乱流理論に基づくモジュールまで、数種用意されている。
Fig.2 に示したように、排出された一次汚染物質は、大気中を輸送、拡散されるうちに、光化学反応などの大気反応により、二酸化窒素(NO2)や酸化Fig.1

大気移流拡散予測システム

性の二次汚染物質オゾン, 過酸化水素(H2O2)などに変換さる。これが硫黄酸化物(SOx), 窒素酸化物

2.1

気象モデル

(NOx),炭化水素類と反応して酸性物質をつくり、まず、化学物質を輸送する大気運動について気一部は粒子化して硫酸、硝酸エアロゾルや有機エ象モデルを用いて計算する。計算に用いる初期条アロゾルなどの二次粒子になる。二次的に生成さ件、境界条件には、世界各国の気象官署で観測される粒子はPM2.5 など微小粒子として分類される。

れたデータを収録し、客観解析（外挿法の一種）

これらガス状、エアロゾル状の汚染物質は、大気した気象客観解析データMeteorological re-analysis

乱流中で地表に沈着除去される（乾性沈着）。さdata（Fig.1の最上段）を用いる。

らに、雲粒の凝結核として、また衝突併合によりこの他、気象客観解析データとしては、日本気雲粒、雨滴に取り込まれ、液相反応を経て酸性雨、象庁の数値予報に基づいて気象業務支援センター

雪、霧になって直接地表に沈着する（湿性沈着）。

が提供している気象予報数値データGPVやヨーロ

大気反応過程については、通常、化学種110個、素ッパ中期予報センター (ECMWF)が提供している反応120種の気相反応、化学種30個、反応および気ERA-40 データ（1957年9月以降）が用いられる。

液平衡40種の液相反応のモジュールが用いられる。

本研究では、主に米国のNCEP-fnl（ Final Analysis

エアロゾルについては、粒子直径別に、粒子生成、

Bulletin of JESC No.43. 201664

衝突併合、成長、気相－エアロゾル相間の平衡、ミッションインベントリを示す。全球については表面反応など多数の過程を含むモジュールが数種

1 °×1 °メッシュで、東アジアなど領域規模で用意されている。

は0.5 °×0.5 °メッシュのものが入手できる。日

本研究では、WRFchemおよびGEOS-chemをイン

本国内については、汚染物質によっては500 mや1ストールして実行する。WRFchemの反応過程、エkmメッシュのものが作成されている。本研究では、

ア ロ ゾ ル の モ ジ ュ ー ル と し て は 、 RADM2 －全球にはGEIA（Global Emissions InitiAtive）、東アMADE-SORGANおよびCBM-Ⅳ－MOSAICを用い

ジ ア 、 日 本 に つ い て は REAS（ Regional Emissionる。

inventory in ASia）およびStreets（2012）のインベントリを用いる。
2.4

グラフィックスモデル

計算結果の時間－空間分布（1-5次元）および時間－空間断面の等値線図（地表濃度の水平分布など）とそれらのアニメーションを描画する。本研究ではGrADS及びVIS-5Dをインストールして用いる。
３．輸送･拡散モデルの構築
発電所や廃棄物焼却施設などの大気環境影響ア
Fig.2

輸送中に発現する様々な形態の大気汚染セスメントにおいて、大気拡散予測にはプルーム、パフモデルが用いられている。これに加えて、｢特2.3エミッションインベントリ

二酸化硫黄(SO2)、NOx、PM2.5 、ブラックカーボン（BC）等の大気汚染物質や温室効果物質など大気微量物質について排出インベントリが整備され、1 ～数年毎に更新されている。Table 1に、種々の研究機関で整備されている各種大気微量物質のエ

殊気象条件下での拡散の予測｣（内部境界層の形成
による「いぶし現象」（フュミゲーション）、逆転層の形成によるトラッピング現象、建物の影響（ダウンドラフト、ダウンウォッシュ））と施設周辺に排煙がぶつかるような高い山などがある場
合には、地形影響の予測が要求される。これらにTable.1 Emission inventories in Asia and the WorldUNFCC inventory has been developed by United Nations for Climate Changes mainly for GHG but also for air pollutants
RAINS・GAINS is developed by International Institute for Applied System Analysis (IIASA) to estimate emissions of air
pollutants including GHG.

EDGAR is developed by National Institute for Public Health and the Environment (RIVM) to estimate emissions of air
pollutants and GHG.

GEIA has been developing inventories of global gas and aerosol emissions as part of International Geosphere - Biosphere
Programme (IGBP).

LTP is a joint research program among China, Japan and Korea to improve understanding of trans- boundary air pollutants inNortheast Asia.

ACESS is developed by David Streets in Argonne Nat. Lab. to support the Aerosol Characterization Experiments (ACE-Asia),
Transport and Chemical Evolution over the Pacific Experiments and the MICS-Asia project.

REAS is developed by Frontier Research Center for Global Change and National Institute for Environmental Studies tounderstand the role of trace constituents in the atmosphere.

EA-Grid is developed by the Ministry of the Environment in Japan to understand transboundary air pollutants in NortheastAsia.

日 環 セ 所 報 No.43. 平成28年

65

対応するためのモデルとして、ISC-PRIMEモデル、子運動を定式化したものである。Fig.3の一番下のMETI-LISモデル、トラッピングモデル、Lyons-Cole

図に、本研究の結果がLuhr and Sawfordの計算結果モデルなどが用いられている。しかし、これらはと良く一致することを示す。

いずれもプルーム、パフモデルを拡張したものであり、予測精度に課題がある。また、地形影響の予測には、電中研モデルが用いられているが、大気 拡 散 に 大き な 影 響 を及 ぼ す 大 気安 定 度 の 効果（温度成層効果）が考慮されていない。
3.1

統計的ラグランジェ拡散モジュール

プルーム、パフモデルの欠点は、不安定成層境界層（混合層、Convective boundary layerという）で顕著になる。混合層中では、地表からの強いプルームによる上昇流域updraft regionと緩やかな下Fig.3 Contour of concentration when pollutant
is emitted from top of the convective boundarylayer.

降流領域downdraft regionの2つで構成されている。
上昇流域での鉛直風速は下降流域での鉛直風速よ3.2

りはるかに強く、また、それぞれの領域の面積比モデルWRFと化学輸送モデルWRF-chem

aupdreft/ adowndraft は < 1になっている。このような鉛直風速の非対称性（大きなskewness歪度）のために、煙の拡散は、上下に対称なプルーム、パフモデルの濃度分布（ガウス分布）から偏倚し、煙の主軸は地面に平行ではなく上下に大きく蛇行する。
Fig.3は、Willmarthの水槽実験とLuhr and Sawfordの数値モデルの結果を比較した。
図に示したように、混合層（厚さδ）中のトップzs=δに排出された場合、煙が徐々に下降して地面に達し、その後徐々に上昇する。そのため、プルーム、パフモデルの最大着地濃度は1/1.6－1/2.0程度の過小評価になってしまう。
本研究では統計的ラグランジェ拡散モジュール（Lagrangian stochastic turbulent diffusion module）を開発した。これは、鉛直風速の確率密度分布の非対称性を考慮し、乱流拡散をマルコフ過程と仮定し、多数の拡散粒子の運動をLengevin equationで記述して、それを数値的に追尾するものである。
混合層の内部では、粒子運動の時間スケール
Lagrangian time scaleが上昇（下降）流域の時間ス

成層乱流拡散モジュールを組み込んだ気象

乱流拡散の理論としては、従来k-εモデルやk-ℓモデルが使われてきたが、温度成層状態（鉛直方向に温度分布があり、流体運動に浮力が作用する状態）によって数千倍変化する乱流拡散を記述するには十分な精度がなかった。我々は、乱流に及ぼす成層効果を室内実験と野外観測で調べ、「成層乱流理論」（Algebraic Stress Model, ASM：Uno,
Wakamatsu & Ueda 1989, Ueda et al. 2012）を提唱してきた。また、このASMをプログラミングして成層乱流拡散モジュールを作成した。
このASMは、これまで特に予測精度が低いとされていた「安定」成層乱流に対しても、高い精度を有する。即ち、安定成層中では、乱流拡散は急激に減衰するが、その影響は鉛直方向と流れ方向、スパン方向とで大きく異なり、運動量拡散と物質の拡散とで異なり、接地気層とその上空（自由大気）とでも異なること、さらに、強い安定成層では逆勾配拡散（低温側から高温側への熱拡散）さえ起こることを、高精度で予測できる（Komori,Ueda, Ogino & Mizushina, 1978, 1983; Ueda,

ケールよりずっと長いことから、上昇（下降）流Mitsumoto

域内の粒子は、そのまま同領域に止まって上昇（下Hanazaki & Hunt, 2004）。現時点では、このASM降）すると仮定でき、それぞれの流域内部での粒&Komori,

1981;

Hanazaki,

1993;が成層乱流理論として最も再現性の高いものの一Bulletin of JESC No.43. 2016

66
つと考える。

これまで、この成層乱流拡散モジュールを組み4.2

WRF-chemモデル

込んだ熱流体力学モデルSIMPLER-3Dを、熱的にWRF（Weather Research and Forcasting）モデルは誘起される局地風（海陸風や山谷風）とその中で大気科学・気象の研究と天気の予報業務の両方のの 光 化 学 大 気 汚 染 の 研 究 に 用 い て き た （ Sha,

ニーズに対応するために設計された次世代数値予Ueda,1990,1991; Chang, et al. 1989,1990, 1991）。

報システムである。データ同化システム、ネステ

本研究では、気象モデルWRFと化学輸送モデル

ィングおよび並列計算が容易に行える。モデルはWRFchemに、この成層乱流拡散モジュールを組み数千kmから数10m四方の領域を対象とした計算に込んだ予測システムを整備した。

対応できる。地球規模の計算も可能である。

このシステムを用いれば、複雑地形上に建設さWRF-chemはWRFモデルに大気化学をカップルれる廃棄物焼却施設や発電所などの排煙拡散の合したモデルである。モデルは、気象現象と同時に理的な予測、特に、これまで予測が困難であった微量ガスとエアロゾルの排出量、輸送、混合、お「地形影響の予測」と｢特殊気象条件下での拡散のよび化学変換をシミュレートする。モデルは、地予測｣に利用できる。特殊気象条件のうち「逆転層域規模の大気質、野外観測の解析、および雲と化の形成によるトラッピング現象」にはWRF-chem学の間の雲スケールの相互作用の研究に使用できがそのまま用いられ、ダウンドラフト、ダウンウる。

ォッシュの予測には、熱流体力学モデル

本研究では、気象モデルWRFと化学輸送モデル

SIMPLER-3Dと統計的ラグランジェ拡散モジュー

WRF-chemに上記の「成層乱流拡散モジュール」をルを組み合わせて、またフュミゲーションの予測組み込んだ。また、東アジア全域について、エミには、WRFモデルと統計的ラグランジェ拡散モジョンインベントリをWRF-chemにリンクすれば、東ュールとを組み合わせて用いれば、これまでになアジア全域の気象及び種々の大気汚染が計算できい高精度予測が可能になる。

る。また、WRF-chemのネスティング機能を用いて、東アジアの一部（例えば、日本、関東地方、川崎４．地球規模、領域規模の大気化学輸送モデル市など）の気象、大気汚染の計算ができる。さらに、大気汚染による気象変化、気象変化に伴う大4.1

GEOS-chemモデル

GEOS-chemは気象データとゴダード地球観測システム（GEOS）をもとに実行するグローバルな３次元大気化学輸送モデル（CTM）であり、NASAのグローバルモデリングチームと関係するオフィ気汚染のさらなる悪化など、気象と大気汚染の相乗作用（フィードバック）が計算できる。
WRFモデル、WRF-chemモデルについての解説書、ユーザーズマニュアルおよびモデル結果の可視化グラフィックスモデルの説明書をTable 2に示す。

スが提供しているものである。本モデルは、世界中の研究グループによって大気組成の問題の広い範囲に適用されている。モデルの科学的な方向性は、国際GEOS-CHEMの運営委員会とユーザワーキンググループによって提供されている。モデル
は、米国NASAの地球科学部門とカナダ国立工学研究評議会からの支援を受けて、ハーバード大学、ダルハウジー大学内に設置されるGEOS-CHEMサポートチームによって管理されている。

Table.2 References of WRF model, WRF-chemmodel and graphics model GrADS.

日 環 セ 所 報 No.43. 平成28年

５．研究結果

3.0

3.0

II

2.5

2.5

THg (ng m -3)

GEOS-chemを用いた、Hgの全球規模の挙動

3.5

THg (ng m -3)

5.1

67

2.0
1.5

1.5
1.0

1.0の研究

2.0

0.5

0.5

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

5.1.1

目的4.5

3.5

4.0

2.5
2.0
1.5
1.0

GEOS-chemにも Hgモデルが 公表されて いること

0.5

THg (ng m -3)のフラックスを文献によって見直し、臭素(Br)及び0オゾン(O3)による無機水銀(Hg )の酸化反応を考慮した。
Fig.4の領域ⅡからⅦに対して、環境省「有害大気汚染物質モニタリング調査」によるHgの測定結1.50.5Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

4.0

5.0

3.5

4.5
4.0

3.0
2.5
2.0
1.5

0.5

V0.0

3.5
3.0
2.5
2.0
1.51.0

産業革命以降蓄積されたHgの土壌からの放出

2.0

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

ライズの修正を検討する。
結果

2.5

1.0IV0.0から、本研究では精度向上を目的としてパラメタ

5.1.2

THg (ng m -3)

3.0

THg (ng m -3)

おり、モデルの検討など取り組みが進んでいる。

3.03.5THg (ng m -3)

欧米を中心に水銀問題に対する関心が高まってIII

0.00.01.0

VI0.50.0

VIIObs.
M1
M2-2

BASE
M2-1

0.0
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

BASE：土壌からのHgﾌﾗｯｸｽ→土壌Hg濃度及び日射量を考慮 Hg0の酸化反応：Br
M1：土壌からHgﾌﾗｯｸｽ→土壌Hg濃度、日射量、及び気温を考慮 Hg0の酸化反応：Br
M2-1：土壌からのHgﾌﾗｯｸｽ→土壌Hg濃度、日射量、及び気温を考慮 Hg0の酸化反応：Br+O3
M2-2：土壌からのHgﾌﾗｯｸｽ→土壌Hg濃度、日射量、及び気温を考慮 Hg0の酸化反応：Br+O3（極域はBrのみ）

果と比較を行った結果を示す（Fig.5）。各グリッドの観測値は、グリッド内に存在する全ての一般環境測定局における測定結果の平均値及び標準偏差（エラーバー）である。夏季のHg濃度について着目すると、Br及び日射量のみを考慮したBASE

Fig.5 日本列島上の各陸面グリッド(図4-2のII～VII)
における、2007年の地表層大気中総水銀(THg)濃度の観測値(Obs.)及びモデル計算値。観測値は、当該グリッド内に存在する全ての一般環境測定局における、｢有害大気汚染物質モニタリング調査｣によるHgの測定結果の平均値(Cave)及び標準偏差(σ; エラーバー)。ただし、Cave ± 3σの範囲から外れる測定結果は外れ値として除外したでは過小評価であったのに対し、パラメータに修正を加えた他の3ケースでは再現性が向上した。

5.2

WRFchemを用いた、大気汚染と都市気象変化の相乗作用の研究－冬季北京のPM 2.5 汚染－
5.2.1

目的

発展途上国の巨大都市では、冬季晴天日に環境基準の数10倍の高濃度エアロゾル汚染が発生する。
本研究は、大気汚染と都市気象変化との相互作用の解明と排出削減効果の評価を目的とする。
5.2.2

結果

エアロゾルは日射の散乱により地表での直達日射量を減少させ、地上気温の低下を招く。また、
Fig.4 GEOS-Chem-Hgによる日本周辺の計算領域。
グリッド解像度は2°(緯度方向) × 2.5°(経度方向)

ブラックカーボン（BC、すす）は日射を吸収する
ため、日出後、午前中にはエアロゾル層（地表から200-300 mの気層）の上面を加熱して逆転層を強68

Bulletin of JESC No.43. 2016

化する。夜間には、通常、地面との長波放射の授６．まとめ

受により地表気温を上昇させるが、冬季北京の場

１）気象モデル、化学輸送モデルとして、WRF、合、夜間の気温低下と高濃度大気汚染のために大WRFchem、GEOS-chemを日本環境衛生センター

量の凝結核から放射霧が生成し、霧粒から天空へのPCクラスタマシンにインストールして予測の長波放射により地表気温は１～２℃低下する。

システムを構築した。多くのユーザが種々の目このように、鉛直方向の熱的構造、特に、逆転層的に対して計算実行できるように、操作マニュの強度と高度を大きく変化させる。これがエアロアルを整備した。

ゾル濃度のさらなる上昇（正のフィードバック）をもたらす。Fig.6には、一例として、2013年1月26２）グラフィックモデルとして、GrADSをインストールして描画システムを構築した。

日から2月1日にかけての北京市の高濃度大気汚染

３）Hgの全球挙動について土壌からの再飛散過程、とそれに伴って生じた都市気象の変化を示す。

金属水銀の酸化過程のモジュールを精緻化し、上図には汚染物質（SO2+NOx）排出量が現状のわが国の有害大気汚染物質モニタリングのデー0.001倍、下図には1.0倍時の地上気温、日射量とBCタおよび海外での観測データがよく再現できるの地上濃度の日変化を示す。

ことを示した。

比 較 す る と 、 BC の 濃 度 上 昇 後 に 汚 染 物 質（SO2+NOx）排出量の濃度上昇が認められ、汚染物質（SO2+NOx）排出量が1.0倍時のものは、1月29日頃から気温及び日射量の上昇が抑えられ、その後に汚染物質（SO2+NOx）排出量の増加が認められた。

【参考文献】
1) Uno,I., Ueda,H. and Wakamatsu,S. 1989:Boundary-Layer Meteorology, 49, 77-98.
2) Ueda, H., Fukui, T., Kajino, M., Horiguchi, M.,
Hashiguchi, H. and Fukao, S. 2012: J. Atmos. Sci.,
69, 323-337.
3) Komori,S., Ueda,H., Ogino,F. & Mizushina,T.
1978: Phys. Fluids, 25, 1539-1546.
4) Komori,S., Ueda,H., Ogino,F. & Mizushina,T.
1983: J. Fluid Mech., 130, 13-26.
5) Ueda,H., Mitsumoto,S. & Komori,S. 1981: Quart.
J.Roy. Met. Soc. 107,561-571.
6) Hanazaki,H, & Hunt,J.C.R. 2004: J. Fluid Mech.,507, 1-42.
7) Hanazaki,H. 1993: J. Fluid Mech., 249, 415-44

Summary
Model calculations of meteorology and air
pollution have become popular in recent years. The
calculations have shown great progress, due toadvances in computer performance and community
modeling systems. In this study, we integrated ourprevious research achievements on physical andchemical processes into the frame work of thecommunity model. Moreover, we constructed an airadvection diffusion model system and set the computerenvironment. Subsequently, first, we tried to simulate
global-scale mercury (Hg) behavior in atmospherewith high accuracy. Second, we conducted
Fig.6 大気汚染と都市気象変化の相乗作用．汚染
calculations of the synergy effects between urban
物質排出量が現状の0.001倍と1.0倍時の地上気温、
climate change and air pollution. The PM2.5 winter
日射量とBC、SO2+NOxの地上濃度の日変化

日 環 セ 所 報 No.43. 平成28年
pollution in Beijing was used as a target. We improved
the prediction accuracy of the behavior of atmosphericHg, by refining parameterizations of re-entrained Hgfrom soil and incorporating the photo-oxidation
reaction of Hg with bromine(Br). In other cases, wederived high concentrations of air pollutant. Moreover,
we showed that black carbon promotes the absorption
ofsolarradiationanddecreases

thesurfacetemperature. A change in the strength and height of theinversion layer causes a higher concentration of air
pollutants.

69