福岡水海技セ研報 第 32 号 2022 年３月
Bull.Fukuoka.Fisheries.Mar.Technol.Res.Cent.No32.March 2022

豊前海の長期環境変動に関する研究
後川 龍男・黒川 皓平・田中 慎也・鹿島 祥平（豊前海研究所）

地球温暖化について世界的に議論される中，日本近海でも長期的な水温上昇傾向が報告されていることから，豊前海における長期水温変動について 48 年分の浅海定線調査データを元に解析した。また同データを用いて長期栄養塩濃度の変動についても解析した。その結果，豊前海の水温は長期上昇傾向，DIN 濃度は長期減少傾向にあることが確認され，48 年間の年平均全層平均値では，水温が 0.89℃上昇，DIN 濃度が 2.38µM 減少したと試算された。水温に関しては夏季以外の上昇傾向が顕著であった。DIN 濃度に関しては夏季を除いて減少傾向にあり，特に秋季と冬季の減少傾向が顕著であった。
キーワード：豊前海，水温，DIN

地球温暖化の議論は近年世界的に盛んとなっており，海が行われている。初期の調査定点は 3～5 点で位置も度々

水温の上昇についても同様である。日本近海でも長期的な変更されてきたものの，1971（昭和 46）年 4 月以降は，現海水温の上昇傾向が報告されており，瀬戸内海の周辺海域

在と同一定点で浅海定線調査が行われている 4）（図 1）。

でも，直近の報告では四国・東海沖で＋1.24℃／100 年，

これを踏まえて使用可能なデータを精査した結果，1972 年

東シナ海で＋1.25℃／100 年

1）の長期的な水温上昇傾向が1 月から 2019 年 12 月までの 48 年間の水温および DIN，
DIP

報告されている。瀬戸内海西部の周防灘の一部である豊前

濃度について，観測データの抽出，取りまとめおよび解析

海でも過去に佐藤，神薗 2）が長期的な水温上昇傾向を報告を行った。

しており，その要因としてレジームシフトが示唆されているが，それ以降の観測データによる検討は未実施である。
また，瀬戸内海の水質汚濁防止対策として 1973 年に瀬戸内海環境保全特別措置法（以下，特措法）が制定され，海域への窒素やリンの排出が強力に規制された結果，栄養塩不足による養殖ノリの色落ちや基礎生産力の低下による漁獲量の減少など，瀬戸内海の栄養塩濃度が低下し水質が大幅に改善した一方で当初想定されていなかった負の影響も指摘されるようになった 3）。これに対応するため 2021 年に特措法が改正され，これまでの栄養塩削減のみを進める対策だけではなく，関係府県知事が策定した計画に基づき，特定の海域への栄養塩供給も可能となった。
こうした状況の中，今後の水産施策に活用する基礎デー
タとして，温暖化に関連する水温と，特措法改正に関係する栄養塩（DIN および DIP）濃度について，豊前海における長期変動をとりまとめたので報告する。

方

法
図1

豊前海では昭和 20 年代から月 1 回の定期海洋観測調査- 43 -

調査定点（数字は Stn 番号，○囲みは代表点）

後川・黒川・田中・鹿島

浅海定線調査は月 1 回，原則として月上旬に行われてった。すなわち豊前海全体の月別平均値については，観いる。今回使用したデータの調査点は図 1 に示した 12測自体の中止による全点欠測のみ欠測値とし，一部定定点，観測層は 0m 層（表層），5m 層，B-1m 層の 3 層点の欠測や，一部定点の一部観測層における欠測の場とし，抽出した観測項目は水温及び DIN，DIP（リン酸合は，これを無視して層別平均値および全層平均値を態リン）濃度とし，DIN については，アンモニア態窒素算出した。代表点別の月平均値については，調査水深 3および硝酸態窒素，亜硝酸態窒素の合計とした。観測項

層 の中 に欠測 値があ る定 点で は全 層平均 値を算 出 せ

目ごとの観測方法を表 1 に示した 4）。使用期間等につず，欠測値とした。また豊前海全体，代表点別とも，欠いては事業報告書等に記録のない場合もあったため，測月がある年の年平均値は欠測値とした。

一部は機器の購入履歴等から推定し記載した。
解析にあたっては，観測値の平均値を算出してデー

結

果タを取りまとめた。豊前海全体の層別月平均値は同月の 12 調査点の同一観測層における観測値の単純平均，

豊前海の水温の解析結果を表 2，3 に，DIN，DIP 濃度の

全層月平均値は層別月平均値の単純平均，層別年平均解析結果を表 4 にそれぞれ示した。表には回帰直線の傾き

値および全層年平均値はそれぞれ層別月平均値およびを「変化率」，回帰直線の R2 乗値から算出された「P 値」，

全層月平均値の単純平均とした。これらを 48 年分取りP＜0.05 となり回帰直線の有意性が 5％水準で確認された

まとめて解析した。

項目について変化率から試算した 48 年間の上昇（下降）幅

水温については更に，海域別の特徴を検討するため図 1を「変化量」として数値を記載した。なお P≧0.05 となりに示した沖合 1 点（Stn.3：沖合）と沿岸 4 点（Stn.4

5％水準で有意性が確認されなかった項目の変化量は「－」：北部，Stn.12：人工島，Stn.10：中部，Stn.1：南で記載した。

部）の計 5 点を代表点として解析した。代表点別の層

豊前海全体の全層平均水温は，年平均値および 2～4 月，

別月観測値はそれぞれの生観測値を使用し，全層月平10～12 月の月平均値で上昇傾向を示し，全層年平均水温は

均値は層別月観測値（0m，5m，B-1m の 3 層）の単純

48 年間で 0.89℃上昇したと試算された（図 2）。層別月平平均，層別年平均値および全層年平均値はそれぞれ層

均水温は 2 月の 0m 層を除き 2～4 月，10～12 月に各層で上別月平均値および全層月平均値の単純平均とした。こ

昇傾向を示した。

れらを 48 年分取りまとめて解析した。

代表点別の全層平均水温は，年平均値及び 3～4 月，10～

データの解析は，各データを 48 年分並べて最小二乗

12 月の全点，1～2 月の北部，2 月の中部，南部， 5 月の人法による線形回帰分析を行い，P 値から有意水準 95％

工島で上昇傾向を示した。層別月平均水温は 3 月，11 月，

で回帰直線の有意性を判断した。有意と判定されたデ12 月に全点の各層で上昇傾向を示した。また，1 月に北部ータ項目のみ上昇（あるいは下降）傾向があると評価の各層，2 月に北部の 5m，B-1m 層および中部，南部の各し，回帰直線の傾きから 48 年間での上昇（下降）幅を層，4 月に沖合と中部の 0m 層を除く全点の各層，5 月に北試算した。

部の B-1m と人工島の 5m 層，10 月に人工島の 0m 層を除く

なお欠測値についてはこれを含めずに回帰分析を行

表1

全点の各層で上昇傾向を示した。

観測項目別の主な観測機器，分析方法，使用期間（一部重複期間あり）
栄養塩

観測項目

水温

DIN
NH4-N

使用期間
観測機器・手法（メーカー）

1972-1976
電気温度計（電気化学計器 K.K.製）
1976-1982
棒状水銀温度計
1982-1985
電気温度計（東邦電探 ET-5D）
1985-1989
CSTD（アレック電子 ADR-500)
1989-2001
STD（アレック電子 AST-1000M）
2000-2015
クロロテック（アレック電子 ACL-208）
2015クロロテック（JFEアドバンテック ASTD-102）

DIP

NO2-N

NO3-N

PO4-P

1972-1983

1972-1974

1972-1980

1972-1974

インドフェノール法

グリース・ロミン法

銅還元法

ドンゼ法

1974-1983

1980-1983

1974-1983

ストリクランド法

Cu-Cdカラム
還元法

ストリクランド法

1984-

テクニコンオートアナライザーⅡ（1984-1998）
オートアナライザートラックス800（1998-）

インドフェノール青
吸光光度法

ナフチルエチレンジアミン
吸光光度法- 44 -

Cu-Cdカラム
還元法

モリブデン青ｰアスコルビン酸還元
吸光光度法

豊前海の長期環境変動

表2

水深

豊前海における水温の長期変動解析結果（豊前海全体（全点平均），沖合（Stn.3），北部（Stn.4））

時期

0m
5m
1月
B-1ｍ
全層
0m
5m
2月
B-1ｍ
全層
0m
5m
3月
B-1ｍ
全層
0m
5m
4月
B-1ｍ
全層
0m
5m
5月
B-1ｍ
全層
0m
5m
6月
B-1ｍ
全層
0m
5m
7月
B-1ｍ
全層
0m
5m
8月
B-1ｍ
全層
0m
5m
9月
B-1ｍ
全層
0m
5m
10月
B-1ｍ
全層
0m
5m
11月
B-1ｍ
全層
0m
5m
12月
B-1ｍ
全層
0m
5m
年平均
B-1ｍ
全層

豊前海全体
（全点平均）
変化率
（℃/年）

P値

0.02130.02370.02430.02310.0201
0.02120.02230.02120.02580.02980.02900.02820.02440.02850.02680.02660.00410.01230.01610.0108-0.0106
0.0022-0.0006-0.0030-0.00490.00420.00180.00040.01390.00770.01090.0108-0.0018
0.00700.00810.00440.02660.03250.03360.03090.02930.03120.03190.03080.04340.04440.04450.04410.01490.01970.02090.01850.14460.1051
0.09500.11290.0581
0.04460.03530.04480.02520.00830.00880.01210.01400.00120.00210.00290.73680.28020.12010.32080.35190.83250.95430.76450.72740.70690.88130.97440.30680.54500.46040.36480.87910.54620.50440.69870.01260.00170.00080.00270.00740.00390.00160.00360.00090.00060.00060.00070.00400.00020.00030.0005沿岸
北部（Stn.4）

沖合（Stn.3）（℃/48年）

変化率
（℃/年）

P値

1.016
1.070
1.017
1.238
1.430
1.392
1.354
1.170
1.370
1.286
1.275
1.278
1.559
1.615
1.484
1.407
1.496
1.530
1.478
2.083
2.130
2.136
2.115
0.717
0.944
1.001
0.887

0.02280.02300.01640.02070.01400.01490.01540.01480.02380.03000.02700.02690.01570.02650.02030.0209-0.01160.00200.0039-0.0019-0.01210.0010-0.0176-0.0096-0.0006
0.0104-0.01140.00130.01410.01220.01140.01250.00410.0047-0.0108-0.00070.02920.0346
0.03390.0326
0.03650.03740.03820.03740.04770.04770.05280.04940.01490.01970.01630.01720.14130.13680.32880.18860.13900.11860.13960.12560.01900.00240.00390.00450.07840.00170.01280.00720.38010.87480.67340.84300.75790.82880.3061
0.39800.96520.41950.50530.90620.32800.39300.53260.22520.74240.69780.47970.95190.00240.00030.00030.00050.00050.00040.00080.00050.00010.00010.00030.00020.00940.00080.02340.0042

変化量- 45 -

変化量（℃/48年）

1.141
1.438
1.297
1.292
1.274
0.977
1.001
1.403
1.663
1.625
1.564
1.752
1.796
1.834
1.794
2.288
2.288
2.535
2.371
0.714
0.947
0.783
0.824

変化率
（℃/年）

P値

0.03730.04400.03940.04030.02310.02700.02900.02640.03300.03680.03550.03520.02950.03540.03840.03440.00850.01370.02420.0154-0.00480.0026-0.0068-0.0030-0.00670.00110.01270.00350.01490.01290.02370.0172-0.00190.00530.00840.00390.02820.03180.03400.03140.03770.03910.0396
0.03880.04340.04230.03980.04180.01870.02300.02550.02230.01510.00420.00750.00700.09460.04700.02500.04720.00150.00040.00070.00070.00110.00010.00010.00020.44800.18770.01410.12680.68030.81640.5626
0.78340.62670.82460.47250.77170.28840.33400.27570.23530.86570.62520.56410.73590.00710.00210.00110.00250.00110.00060.00060.00070.00020.00030.00030.00030.00050.00000.00010.0001

変化量（℃/48年）

1.791
2.114
1.892
1.933
1.296
1.393
1.265
1.584
1.767
1.702
1.690
1.415
1.702
1.842
1.653
1.162
1.355
1.528
1.631
1.505
1.812
1.876
1.898
1.862
2.082
2.030
1.911
2.008
0.896
1.105
1.223
1.069

後川・黒川・田中・鹿島

表3

水深

時期

0m
5m
1月
B-1ｍ
全層
0m
5m
2月
B-1ｍ
全層
0m
5m
3月
B-1ｍ
全層
0m
5m
4月
B-1ｍ
全層
0m
5m
5月
B-1ｍ
全層
0m
5m
6月
B-1ｍ
全層
0m
5m
7月
B-1ｍ
全層
0m
5m
8月
B-1ｍ
全層
0m
5m
9月
B-1ｍ
全層
0m
5m
10月
B-1ｍ
全層
0m
5m
11月
B-1ｍ
全層
0m
5m
12月
B-1ｍ
全層
0m
5m
年平均
B-1ｍ
全層

豊前海における水温の長期変動解析結果（人工島（Stn.12），中部（Stn.10），南部（Stn.12））

人工島（Stn.12）
変化率
（℃/年）

P値

0.01720.01970.01800.01960.02100.02540.02340.02330.02740.03190.03080.03230.02680.03090.02710.03250.00640.02700.02250.0261-0.00480.02190.00810.0188-0.0034
0.01420.02220.02120.01400.01320.01760.0242-0.00310.00920.01140.00870.02080.02940.03000.02850.02480.02680.02890.0256
0.03610.0336
0.04340.03410.01450.02820.02290.01870.32800.2646
0.28120.85600.11140.05600.08880.08120.04670.04920.01770.01730.00750.00450.00540.00300.60570.02470.05810.02750.68030.07400.45520.11900.80040.29730.2693
0.09600.31510.38830.26700.08590.81080.46230.40160.51220.06200.00590.00440.00730.04350.03550.01380.04010.01010.01300.00210.01790.01200.00010.00020.0035沿岸
中部（Stn.10）（℃/48年）

変化率
（℃/年）

P値

1.316
1.532
1.479
1.548
1.287
1.484
1.300
1.558
1.295
1.253
1.413
1.442
1.367
1.192
1.288
1.389
1.231
1.734
1.613
2.082
1.636
0.697
1.355
1.099
0.898

0.01560.02060.02050.01890.02630.02670.02780.02690.02730.03400.03020.03170.02220.02870.02600.02600.01100.01740.01340.0207-0.0048-0.0012-0.0034-0.0001-0.0058
0.0021-0.01190.00320.01710.00170.00250.0070-0.0057
0.01090.01190.00800.02430.03690.03510.03500.02940.02950.03310.02930.04510.04740.0471
0.04650.01570.01810.01940.02100.32630.19570.19050.23010.03650.02790.01910.02630.03040.00970.01370.01280.08370.00370.00450.01120.40080.16820.26520.07870.68030.91800.77880.98870.71580.88780.93100.80570.24160.91560.88410.62210.6547
0.39640.3776
0.51630.03750.00080.00090.00130.01200.01230.00520.01380.00150.00090.00080.00090.00430.01480.00110.0015

変化量- 46 -

南部（Stn.1）（℃/48年）

変化率
（℃/年）

P値

1.264
1.280
1.334
1.292
1.310
1.630
1.449
1.522
1.378
1.247
1.250
1.168
1.771
1.685
1.679
1.410
1.414
1.588
1.407
2.167
2.274
2.259
2.233
0.754
0.867
0.931
1.007

0.02220.0216
0.02650.02140.02530.02330.02490.02380.03270.03210.03350.03280.02900.02960.02700.02740.01680.01830.01980.0183-0.0048-0.0110-0.0058-0.0108-0.00930.0020-0.0016-0.0009
0.02500.00540.00970.01290.00020.00200.00970.00020.02950.03910.03680.03540.03510.03520.03310.03480.0466
0.0486
0.04130.04920.01870.02100.02040.02190.14790.13440.05880.13950.02620.04380.02230.03500.01280.01340.00700.01100.01200.00260.00410.00500.22340.16200.13590.15390.68030.39020.62890.36380.60380.89870.89650.94890.13200.76270.55090.42440.98800.86420.42020.98650.02540.00180.00100.00510.00180.00310.00430.00300.00350.00040.00200.00040.00050.00190.00040.0008

変化量

変化量（℃/48年）

1.216
1.119
1.197
1.140
1.569
1.542
1.610
1.573
1.393
1.420
1.295
1.317
1.418
1.877
1.768
1.700
1.686
1.692
1.587
1.669
2.237
2.332
1.983
2.363
0.896
1.006
0.980
1.053

豊前海の長期環境変動

表4

豊前海における栄養塩類の長期変動解析結果

DIN
水深

時期

0m
5m
1月
B-1ｍ
全層
0m
5m
2月
B-1ｍ
全層
0m
5m
3月
B-1ｍ
全層
0m
5m
4月
B-1ｍ
全層
0m
5m
5月
B-1ｍ
全層
0m
5m
6月
B-1ｍ
全層
0m
5m
7月
B-1ｍ
全層
0m
5m
8月
B-1ｍ
全層
0m
5m
9月
B-1ｍ
全層
0m
5m
10月
B-1ｍ
全層
0m
5m
11月
B-1ｍ
全層
0m
5m
12月
B-1ｍ
全層
0m
5m
年平均
B-1ｍ
全層

変化率
（µM/年）

P値-0.0708-0.0699-0.0743-0.0717-0.0313-0.0323-0.0350-0.0328-0.0651-0.0578-0.0618-0.0546-0.0476-0.0326-0.0380-0.0394-0.0250-0.0282-0.0393-0.0308-0.0048-0.0114-0.0227-0.0151-0.0055-0.0091-0.0545-0.0230-0.0032-0.0081-0.0514-0.0209-0.0177-0.0297-0.0610-0.0361-0.0593-0.0661-0.0959-0.0738-0.1005-0.1006-0.1311-0.1107-0.0718-0.0740-0.0983-0.0814-0.0422-0.0431-0.0635-0.04960.00130.00150.00080.00110.01640.01570.01080.01340.00060.00120.00120.00330.00050.00770.00290.00190.00410.00400.00010.00110.68030.15890.02280.0712
0.79290.49380.01540.15710.72530.35840.0006
0.03140.38810.12580.0267
0.09140.03960.02490.00680.01710.00740.00530.00200.00400.03830.03170.00670.01990.00010.00000.00000.0000DIP（µM/48年）

変化率（µM/年）

P値-3.399-3.357-3.567-3.441-1.500-1.549-1.682-1.577-3.125-2.774-2.967-2.623-2.285-1.565-1.822-1.891-1.1982-1.3516-1.8847-1.4782-1.0893-2.6166-2.4654-1.0035-2.927-2.845-3.171-4.605-3.541-4.824-4.827-6.292-5.315-3.449-3.551-4.716-3.905-2.025-2.069-3.050-2.382-0.0028-0.0006-0.0015-0.0016-0.0034-0.0015-0.0051-0.0033-0.0007-0.0005-0.0008-0.00070.00030.00010.00020.0002-0.0006-0.0006-0.0013-0.0008-0.0048-0.0001-0.0004-0.00030.0000-0.0001-0.0003-0.0001-0.0021-0.0020-0.0049-0.00300.00180.0011-0.00100.00060.00180.00180.00070.00150.00190.00120.00030.00110.00090.00060.00030.0006-0.00030.0000-0.0011-0.0005

0.12520.60960.23290.23830.23000.44910.16780.23310.53700.59870.49580.53570.71650.85130.76450.76600.44180.55220.14940.32840.68030.92080.52320.62540.98020.96980.9226
0.94620.06920.06890.00100.0076
0.19440.47290.66060.70250.19010.20480.65000.31410.2356
0.57710.84680.51000.52100.70290.83130.67020.71390.94590.21940.5323

変化量- 47 -

変化量（µM/48年）-0.2331-0.1442-

後川・黒川・田中・鹿島

豊前海の全層平均 DIN 濃度は，年平均値および 6 月，7

り大きい。福岡の気温の上昇幅は 2.45℃／100 年 1）で先述

月，9 月を除く月平均値で減少傾向を示し，全層年平均 DINの水温上昇幅より大きいことが報告されており，水深が浅濃度は 48 年間で 2.38µM 減少したと試算された（図 3）。

い豊前海は気温の影響を受けやすいため，東シナ海北部や特に 10～1 月の減少幅は 3.44～5.32µM と試算され，年平

四国・東海沖より水温上昇幅が大きくなったと考えられる。

均値の減少幅よりも大きくなった。

逆に，4 月の Stn.3 や Stn.10 等，表層水温のみ上昇傾向を層別でも，全層平均 DIN 濃度が減少傾向を示した月のほ

示さない事例は，直前の気象の影響を受けたために生じた

とんどで各層とも減少傾向を示したが，全層平均値では減可能性があり，豊前海のような浅海域において表層水温の少傾向が示されなかった 6 月，7 月，9 月でも，B-1m 層のみで変動傾向を評価する際には注意が必要である。

み減少傾向を示した。一方全層平均値で減少傾向を示した

次に月別の水温上昇傾向を見ると，2～4 月および 10～

8 月は，0m，5m 層で減少傾向が示されず，B-1m 層のみ減少

12 月の水温が上昇しており，特に 12 月の水温上昇が顕著傾向を示した。

であった。福岡の気温は 3～5 月（春季）と 9～11 月（秋豊前海の DIP 濃度については，年平均値では有意な変動

季）に，東シナ海北部や四国・東海沖の水温は 10～3 月に，傾向を示さず，増加あるいは減少の傾向は認められなかっそれぞれ上昇幅が他の時季より大きいと報告されている

た（図 4）。全層月平均値では，8 月のみ減少傾向を示し，

1）

48 年間で 0.14µM 減少したと試算された。層別月平均値でとしては，もともと東シナ海北部や四国・東海沖の水温上もほとんどの層，月で有意な変動傾向は示されなかったも昇幅が他の時季より大きいことに加え，9～11 月の気温上のの，8 月の B-1m 層のみ減少傾向を示した。

昇幅が大きいため水温が以前より下がりにくくなっている。すなわち豊前海の 12 月の水温上昇幅が特に大きい理由

可能性がある。

考

周辺海域の水塊が豊前海に及ぼす影響としては，1 月に察

北部のみで水温上昇傾向が確認されていることが挙げられ
本研究により確認された豊前海の全層年平均水温の上昇る。北部（Stn. 4），沖合（Stn.3），豊前海全体および幅は 48 年間で 0.89℃であった。1973 年 1 月から 2003 年

福岡県の筑前海側で行われる浅海定線調査のうち豊前海に

12 月までの豊前海における浅海定線調査データを解析し最も近い定点（図 1 の筑前海 Stn.9）6）における 1999 年

た佐藤，神薗 2）によれば，解析期間の 28 年間で豊前海のから 2019 年までの 1 月上旬の水温（3 層平均水温）を図 5

水温は約 0.57℃上昇していると報告されており，本研究にに示した。これらを比較すると，1 月上旬の水温は概ね

よればその後も現在まで上昇傾向が続いていることが確認された。また豊前海に隣接する山口県の周防灘海域でも1972 年 4 月から 2003 年 3 月までの 30 年間で約 0.71℃上昇していることが報告されている 5）。いずれも解析期間が異なるため今回得られた上昇幅との差については論じないが，豊前海を含む周防灘およびその周辺海域で水温が長期にわたって上昇傾向にあることは確実である。
この豊前海の水温上昇幅 0.89℃／48 年は，同期間に換算して比較すると気象台から報告されている東シナ海北部
や四国・東海沖の水温上昇幅 1.24～1.25℃／100 年 1）よ図2

豊前海全体の全層年平均水温（1972～2019）- 48 -

図3

豊前海全体の全層年平均 DIN 濃度（1972～2019）

図4

豊前海全体の全層年平均 DIP 濃度（1972～2019）

豊前海の長期環境変動

豊前海平均＜北部＜筑前海 Stn.9 となっており，筑前海は特措法により水質改善が着実に進んできたことを示す一

Stn.9 の水温は豊前海平均よりも約 3℃高く，北部の水温方，同じ瀬戸内海の播磨灘等で報告される DIN の減少によは豊前海平均よりも約 1℃高かった。また伊予灘や豊後水るノリ色落ち被害 3）など，豊前海でもノリ養殖への影響が道など南方由来の海水の影響を受けやすいと思われる沖合懸念されることを示している。豊前海区のノリ養殖経営体

よりも北部のほうが高水温となる年が多く，1 月の水温が数は，1974 年には 700 経営体を上回っていたが 7），現在は北部＞沖合となった年では両定点の水温差は 2℃を超える

数経営体まで減少している。豊前海では特にノリ養殖の最こともあった。すなわち，豊前海全体の長期変動傾向とは盛期である 10～1 月に DIN 濃度の減少幅が大きいことか異なり北部のみで 1 月に水温の上昇傾向が見られたのは，ら，DIN の減少が海区内のノリ養殖業の衰退につながった

関門海峡から流入する筑前海由来の高水温の水塊による影要因の一つである可能性がある。また DIN が減少傾向にあ響が北部で特に大きくなることが要因と示唆された。

ることで，海域の基礎生産を担う珪藻類などの植物プラン

月別の水温上昇傾向から，豊前海の海況を過去と比較しクトンも減少傾向にある可能性がある。カキの身入りとクた場合，「秋の水温低下時期が遅れてきている」「冬が早ロロフィル濃度との相関が中川ら 8）によって示唆されていく終わり春の水温上昇時期が早まっている」という傾向がるが，本研究で DIN 濃度の減少幅が特に大きかった 10～１認められた。こうした海況の変化は漁業にも影響を与えて月はカキの出荷時期と一致することから，DIN の減少が植いる可能性がある。例えばカキ養殖では高水温化に伴う秋物プランクトンの減少につながれば，出荷時期のカキの身

7），実際過去には入りに悪影響を及ぼす可能性もある。今後この時期の植物

10 月下旬頃であったカキの販売開始時期が，近年では 11

プランクトン量の増減，あるいはクロロフィル濃度の長期

月中旬以降となっている。また豊前海におけるノリ養殖の変動傾向についても解析する必要がある。

季のカキの身入り遅れが示唆されており採苗時期は，1970 年代から 2000 年ごろまでは 10 月上旬だ

一例として 8 月の底層 DIN 濃度を図 6 に示したとおり，ったが，2000 年台半ばには 10 月中旬となり，近年では 11

本研究において 6～9 月の DIN 濃度が底層のみで減少傾向

月 1 日頃の採苗となっている

6）。採苗方法の違いによる差を示した要因は次のとおりと考えられる。すなわち，夏季

はあるものの，ノリの採苗時期の変化は，秋の水温低下時

は成層が発達し海水の上下混合が起こりにくいため，夏季期が遅れていることによる影響と考えられる。これらの他の底層 DIN の供給源は陸域由来よりも底質からの溶出が主にも，水温が有意に上昇した季節に関連する代表的な魚種体と考えられる。底質中の有機物量は過去と比較して減少として，春のコウイカや，秋のガザミ，サワラ，冬のカレしており底質が改善したことが報告されている 3）ことかイ類などの来遊や漁獲にも影響が及んでいる可能性が考えら，夏季の底層 DIN の主供給源である底質が改善されたこられる。水温上昇と水産業との関連については今後，市場とにより底質からの溶出が減少傾向となり，6～9 月の底層への魚種別水揚げ状況や，各漁業種の操業実態，特に操業

DIN 濃度が減少したと考えられる。底質改善による栄養塩

位置が固定されている枡網に着目した魚種別の来遊状況を

溶出量の減少は， 8 月の底層で DIP 濃度も減少傾向を示し解析する等の方法で，より詳細に影響を検討する必要があたことからも示唆される。一方で，覆砂によって窒素溶出るだろう。

速度が大きく低下することが実験的に確かめられているこ

栄養塩のうち，DIN 濃度については減少傾向が認められとから 9－11），底質改善を目的として豊前海で行われてきた

た一方，DIP 濃度については減少傾向が認められなかった。

覆砂の効果により夏季の底質からの栄養塩溶出量が減少し

DIN 濃度については夏季を除いて減少傾向が認められ，特た可能性もある。
以上述べてきたように，豊前海の長期環境変動としては，に底層では一年を通じて減少傾向が認められた。このこと

図5

1 月の水温比較（豊前海，筑前海 Stn.9）

図6- 49 -

8 月の平均底層 DIN 濃度

後川・黒川・田中・鹿島

水温の長期的な上昇傾向と，DIN 濃度の長期的な減少傾向

239-244．

が挙げられる。長期的な水温の上昇傾向については地球規５）和西昭仁．山口県周防灘海域における最近 30 年間

模での気候変動と関連すると思われるため人為的な制御はの水温変動．山口県水産研究センター研究報告

難しいものの，昨今全国各地で導入され豊前海でも試験を2004；2：1-6．

開始した ICT ブイのようなリアルタイムの海況観測システ

６）平成11年度～令和元年度福岡県水産海洋技術セン

ムを整備し，漁業関係者が水温などの海況情報にアクセス

ター事業報告．福岡県水産海洋技術センター，福しやすい環境を整備することで，実際の海況に合わせて効岡．2001-2021．

率的な操業や養殖管理ができる体制をとることが有効と考

７）第21次福岡農林水産統計年報水産編．九州農政局福えられる。また DIN 濃度の長期的な減少傾向については，

岡統計情報事務所，福岡．1974．

現状を関係機関で共有するとともに，水産部局からは特措

８）中川浩一，俵積田貴彦，中村優太，近年の「豊前海法の改正に伴い漁業生産の維持増大につながる提案を行っ一粒かき」の成育状況と漁場環境との関係．福岡県ていく必要があると考えられる。提案の一例としては，赤水産海洋技術センター研究報告

潮発生のリスクが少なく DIN 濃度の減少幅が大きい秋～冬

114．
９）神薗真人，江藤拓也，上妻智行．覆砂による豊前海

季に，沿海下水処理場での栄養塩管理運転を行い海域に栄養塩を供給する試みなど

12）が挙げられるだろう。

の底質改善効果．福岡県水産海洋技術センター研究報告文2009；19：109-

1994；2：129-134．

10）江藤拓也，中川浩一，佐藤博之．豊前浅海域におけ献る覆砂による底質改善効果．福岡県水産海洋技術センター研究報告

１）福岡管区気象台．九州・山口県の気候変動監視レポー

1999；9：61-65．

11）江藤拓也，佐藤利幸，長本篤，上妻智行．豊前海沖

ト2020．
２）佐藤利幸，神薗真人．豊前海の水温・塩分変動と漁合域における覆砂による底質改善効果．福岡県水獲量の変動．福岡県水産海洋技術センター研究報

産海洋技術センター研究報告2006；16：115-120．

告12）水産研究・教育機構瀬戸内海区水産研究所他．ノリ，

2006；16：121-127．

ワカメ養殖場における栄養塩供給技術実証試験事

３）瀬戸内海における今後の環境保全の方策の在り方について（答申）．中央環境審議会，東京．2021．

内海区水産研究所，広島，2018．

４）神薗真人．周防灘西部の海洋観測について．昭和58
年度福岡県豊前水産試験場研究業務報告

例集．国立研究開発法人水産研究・教育機構

1985：- 50 -

瀬戸