REDD+

Reducing Emission from Deforestation
and Forest Degradation-plus

応用講習a

平成25年度第４章
森林炭素変化量の推定について（一社）日本森林技術協会
七海 崇
1

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

目次
1. 森林炭素変化量の推定方法について
2. 参照排出レベルの設定について
3. まとめ4. 実習（GISによる２時点地図の合成～マトリックスによる変化量の解析）

2

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

1. 森林炭素変化量の推定方法について
2. 参照排出レベルの設定について
3. まとめ4. 実習（GISによる２時点地図の合成～マトリックスによる変化量の解析）

3

REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

土地利用に関連した活動からの温暖化ガス排出量を評価する
IPCCの基本的な計算式

排出量の推定活動度データActivity data

排出係数

Emission factor

換算値

UNFCCCのもとでモニタリングシステムは、森林面積及びその変化だけでなく、炭素蓄積量及びその変化までのデータを提供する。

4

REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

森林炭素変化量の推定フロー1990年 森林分布図

森林インベントリ情報
１巡目 1991-1995年
２巡目 1996-2000年
４巡目 2006-2010年

2000年 森林分布図
2010年 森林分布図

マトリックスを用いた解析が有効！

森林動態地図森林面積変化の把握

森林タイプ別の蓄積

Emission factor

Activity data
炭素蓄積変化量5REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

マトリックスによる森林面積変化量の推定
森林タイプ （２０００年）
単位：CtEvergreen
forest （rich）Category of Forest
Evergreen Evergreen Evergreen
Rrehabilita Deciduous
forestforestforesttion forest
forest（rich）（medium） （poor）

Year 2000

Mixed
Mixed
broadleaf
timber and Coniferous
Mangrove&bamboo
forestforestconiferous
forestforestLimestone
Plantation
forestLimestone
area (noforest)

Bare land,shub land,
Residential
Water body
Other land
fragmented
areatrees

Grand Total

14,537-34,093-4,552-61,960

0-5,660

0000000-22,157

00-2,940-116,826914

0-190,937-539,492

0-3,576-32,688000-982-613-2,305-894,207-19-6,542-156,667-1,827,115

Evergreen
forest （poor）

0

15,269Rrehabilitation
forest6,101

11,103

Evergreen
forest（medium）

Category ofForest

Forest remaining Forest
・Degradation Forest
・Conservation Forest
・Enhancement Forest
・Sustainable FM Forest-62,052-110,588

0-464

9,029

000-62-169-4,990-338,247

0-3,667-56,681-552,624

2,224-80,949

0

161

7,520

000

996

0-1,620-170,550

0-1,211-35,910-262,137

0000000000000000

431-7,56907,685

3,032

000231-92-714-76,732-1-479-15,750-88,509-647-81,499

0-3-220,984

000-1,681-255-86-208,340-3-1,105-52,055-566,294

0000000000000

00000000000000000000000000000000

000

Deciduous
forest00

Bamboo forest

01,449Mixed timberand bamboo
forest0364

Coniferous
forest00

Mixedbroadleaf &coniferous

00

Mangroveforest199 Limestone
forest0000104-95-18,455

000000-252,364-1,533-16,309-138,619

0-596-6,901-434,769

Plantation

000

0

0000

000

1,0840-840

00-202

420

579

336

1,131

0

0

0

0

0

0

2,384

1,582

0

0

0

0

0

6,012

2,068

35,232

Water body

0

1

7

Residential
area

00

Other land

13

36,838Grand Total

23,633

66,845Year

Limestone
area (noforest)
Bare land,shub land,
fragmentedtrees

Category of
Non Forest

森林タイプ（１９９０年）

BambooforestCategory of
Non Forest

Afforestation / Reforestation
・Enhancement Forest
26,213 2,099,813

Deforestation

078,228

81,328000

47,949162,31400

0002,533,146

744

03

1,117

000

487

1,172

0

0

0

0

0

3,529

0

1,002

0

10

1

0

0

0

24

4,109

0

0

0

0

0

5,147

26,704200,00304,540114,830000

39,64023,129

0

0

0

0

0

445,697-202,368 1,402,182

080,925-36,817

000-163,380

190,728-26,024 -1,849,693-24-13,602-327,107-854,701

6

REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

どの様にマトリックスを作成するか？
1990
Forest Map 2000
2010

マトリックスの作成手順
２時点地図のオーバーレイ
処理 （インターセクト）
GISによる処理
ポリゴン面積計算

クエリによるテーブル集計
森林面積変化マトリックス

表計算ソフトによる集計
ピボットテーブル作成（面積変化マトリックス）
7

REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

マトリックスでどの様な結果が得られるか？
（５時点森林分布図を利用した例）
(ha)

(ha)

350,000

150,000

300,000
100,000

Limestone forest

200,000

Mixed timber andbamboo

150,000

Rehabilitationr
forest100,000

同じデータソース

250,000

50,000
1st D
2nd D

090～95’95～00'

00～05'

05～10’

EnrichReforest(50,000)

Poor forest(100,000)

50,000Medium forestRich forest

01990

1995

2000

2005

2010

(150,000)（例）A省の森林減少面積と森林回復面積の推移（例）A省の森林面積の推移

8

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

1. 森林炭素変化量の推定方法について
2. 参照排出レベルの設定について
3. まとめ4. 実習（GISによる２時点地図の合成～マトリックスによる変化量の解析）

9平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

排出量削減の考え方
削減努力がない場合の予想排出量
排出量
実際の排出量
予想排出量

排 出 量

REDD+排出削減量

削減努力

年
10

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

参照排出レベル設定における留意事項
REL設定に必要な過去の森林分布図の時点数
REL設定に必要な過去の森林分布図の時点数⇒ Robustness（頑健性）やUncertainties（不確実性）を考慮
すると3時点（変化量では2つ）では少なすぎる？
⇒ RELの設定方法（外挿モデル）について検討の余地がある。
REL設定のスケールによるトレンドの違い⇒ Sub-Nationalアプローチを検討する余地がある。
⇒ RELとRLを分離して解析すると、より傾向がわかる。
（マトリックスを用いた解析が有効）排出量の推定における排出係数（emission factor）の乗算方法
排出量の推定における排出係数（emissionfactor）の乗算方法⇒ 各時点の排出係数を個別に用いるか、平均値を用いるか？

11

REDD+Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

REL設定に必要な過去の森林分布図の時点数について
REL
設定に必要な過去の森林分布図の時点数について

A省の森林減少面積の推移（３時点森林分布図から作成）
A省の森林減少面積の推移（５時点森林分布図から作成）⇒RELの設定方法（回帰モデル、平均値等）によって得られる
インセンティブに影響があることにも留意
12

REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

REL設定のスケールによるトレンドの違い（１）REL
設定のスケールによるトレンドの違い（１）140,000
120,000

過去5年の森林減少・劣化量を全国集計

100,000
80,000

全国値

60,00040,00020,0000National
アプローチ

90～95’ 95～00' 00～05' 05～10’

10,000

過去5年の森林減
少・劣化量を集計
（A省）

2000

5,000

省別にRELを外挿する

20100

過去5年の森林減
少・劣化量を集計（B省）

Sub-National
アプローチ

1990

25～30’

過去5年の森林減
少・劣化量を集計
（C省）

1,000
19902000
2010

2030

1990
2000

500

Σ省A+省B
＋・・

2010
0199020002010
4,000

2030
1990

2,000

2000

0

2010

13

REDD+

平成25年度 応用講習a

Reducing Emission from Deforestation and Forest Degradation-plus

REL設定のスケールによるトレンドの違い（２）REL
設定のスケールによるトレンドの違い（２）(M Ct)
60

全国(M Ct)

40150

20100

050
090～95’95～00’00～05’05～10’50100(M Ct)

2060

404060

2080

0

北部

15020(M Ct)

Deforestation

Forest degradation

4060

Carbon enrichment

Reforestation

60

408020中部⇒⇒

地域別にRELとRLを分離し解析するとトレンドがわかりやすい。（さらに省・県別ではどうか？）
Driving forceの影響が推計できるREL/RLを設定できる。

020
406080

南部
14REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

排出量の推定における排出係数（emission
排出係数（emission factor）の乗算方法factor）の乗算方法
リモートセンシング情報から得られた森林変化に関する活動量を炭素に換算する（例えば）
t1からt2において常緑林(rich)が100haだけ常緑林(medium)に減少した。
常緑林(rich)の平均材積が400m3/ha
常緑林(medium)の平均材積が200m3/ha↓
森林劣化量=100×(400-200)=20,000m3⇒

このように単純に平均材積を用いて森林減少・劣化量を計算してしまうと、排出量を過少に推定してしまう場合がある。
15

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

1. 森林炭素変化量の推定方法について
2. 参照排出レベルの設定について
3. まとめ4. 実習（GISによる２時点地図の合成～マトリックスによる変化量の解析）

16

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus平成25年度 応用講習a

3．まとめ
活動度データ（森林分布図等）は、時点数に留意が必要
排出量の推定やREL排出量の推定やREL設定では、スケール感に留意が必要
設定では、スケール感に留意が必要（国レベル、準国レベル、プロジェクトレベル）

RELとRLRELとRLを分離して解析すると地域のトレンドがわかりやすいを分離して解析すると地域のトレンドがわかりやすい⇒ マトリックスによる森林面積動態の解析が有効
森林変化マトリックスでは、面積変化より蓄積の変化に着目
排出係数を適用する際は、適用方法に留意が必要
RELの設定方法（回帰モデル、平均値・・・）は、検討が必要
RELの設定方法（回帰モデル、平均値・・・）は、検討が必要
その他、National
その他、National Circumstances（国状）を考慮する必要性
Circumstances（国状）を考慮する必要性
17

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 応用講習a

1. 森林炭素変化量の推定方法について
2. 参照排出レベルの設定について
3. まとめ4. 実習（GISによる２時点地図の合成～マトリックスによる変化量の解析）

18