(57)【要約】南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法及びその応用を開示し、過去Ｎ年間にアルゼンチンマツイカが分布していた南西太平洋の南部海域の月別の南極振動指数ＡＡＯ値を取得するステップと、時系列分析方法を利用して、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去Ｎ年間の月別のＡＡＯ値に対して相関性分析を行い、統計的に有意な相関を示す月のＡＡＯ値を選択し、選択したこれらの月のＡＡＯ値を順に番号付けし、順にｘ１、ｘ２、ｘ３…ｘｍとするステップと、任意１〜ｍ個の気候因子に対して、多変数線形方程式を利用して複数のアルゼンチンマツイカ資源豊度予測モデルを作成し、各予測モデルのＰ値を計算するステップと、Ｐ値が最も小さいモデルを最適モデルとして選択するステップと、を含む。本発明の方法によれば、南極振動ＡＡＯを利用してアルゼンチンマツイカ資源豊度の予測を実現することで、漁業生産では指導の役割を果たし、漁獲効率を向上させ、漁獲コストを低減させ、将来性が期待できる。【選択図】図１
2021520575.tif 000002
【特許請求の範囲】
【請求項１】
電子機器に適用される南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法であって、
（１）過去Ｎ年間にアルゼンチンマツイカが分布していた南西太平洋の南部海域の月別の南極振動指数ＡＡＯ値を取得するステップと、
（２）時系列分析方法を利用して、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去Ｎ年間の月別のＡＡＯ値に対して相関性分析を行い、そのうち統計的に有意な相関を示す月のＡＡＯ値を選択し、選択したこれらの月のＡＡＯ値を、アルゼンチンマツイカ資源豊度に影響を与える気候因子とし、選択したこれらの月のＡＡＯ値を１、２、３…ｚ…ｍの順で番号付けし、これらの月のＡＡＯ値を順にｘ１、ｘ２、ｘ３…ｘｚ…ｘｍ（ｍは選択した月のＡＡＯ値の数）とするステップと、
（３）ｘ１、ｘ２、ｘ３…ｘｍのうちの任意１〜ｍ個の気候因子に対して、多変数線形方程式を利用して計２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルを作成し、各予測モデルの統計的Ｐ値を計算し、アルゼンチンマツイカ資源豊度予測モデルの式は、
ｌｎ（ＣＰＵＥ）＝ａ＋ｂ１＊ｘ１＋ｂ２＊ｘ２＋ｂ３＊ｘ３＋…＋ｂｚ＊ｘｚ＋…＋ｂｍ＊ｘｍ（ＣＰＵＥは１隻あたりの年間漁獲量、ａは定数、ｂ１、ｂ２、ｂ３、…ｂｚ…、ｂｎはそれぞれｘ１、ｘ２、ｘ３…ｘｚ…ｘｍに対応する係数）であるステップと、
（４）上記２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルのうち、統計的Ｐ値が最も小さいモデルを最適モデルとして選択するステップと、を含む
ことを特徴とする南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項２】
最適モデルを得た後、該最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得して最適モデルに入力し、最適モデルがアルゼンチンマツイカ資源豊度を出力すると、アルゼンチンマツイカ資源量の予測が完了する
請求項１に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項３】
アルゼンチンマツイカ資源豊度予測モデルが気候因子ｘｚに対して作成される場合、ｂ１、ｂ２、ｂ３、…ｂｚ−１、ｂｚ＋１、…、ｂｎはいずれも０である
請求項１に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項４】
前記統計的に有意な相関とは、算出されたＰ値＜０．０５である
請求項１に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項５】
ステップ（２）では、統計的に有意な相関を示す２つの月のＡＡＯ値、すなわち過去２年目の５月のＡＡＯ値及び過去１年目の８月のＡＡＯ値を選択し、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去２年目の５月のＡＡＯ値との相関性は有意であり、正の相関を示し、相関係数がそれぞれ０．５４９であり、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去１年目の８月のＡＡＯ値との相関性は有意であり、負の相関を示し、相関係数が−０．４９７である
請求項１に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項６】
ステップ（３）では、以下の計３個の予測モデルを作成し、
１）予測モデルＩ
気候因子の過去２年目の５月のＡＡＯ値ｘ１に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７３０５＋０．４５５４＊ｘ１であり、
Ｆ値は６．０３４７であり、Ｐ＝０．０２７７＜０．０５であり、
２）予測モデルＩＩ
気候因子の過去１年目の８月のＡＡＯ値ｘ２に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７２８１−０．４６９４＊ｘ２であり、
Ｆ値は４．５４８０であり、Ｐ＝０．０４９＜０．０５であり、
３）予測モデルＩＩＩ
気候因子の過去２年目の５月のＡＡＯ値ｘ１及び過去１年目の８月のＡＡＯ値ｘ２に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７０４６＋０．３９０８＊ｘ１−０．３８３６＊ｘ２であり、
Ｆ値は５．５１３５であり、Ｐ＝０．０１８＜０．０５である
請求項５に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項７】
ステップ（４）では、予測モデルＩＩＩを最適モデルとして選択し、最適モデルは、
ｌｎ  （ＣＰＵＥ）＝６．７０４６＋０．３９０８＊ｘ１−０．３８３６＊ｘ２である
請求項６に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法。
【請求項８】
電子機器であって、１つ又は複数のプロセッサ、１つ又は複数のメモリ、１つ又は複数のプログラム、及びデータ収集装置を備え、
前記データ収集装置は最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得するために使用され、前記１つ又は複数のプログラムは前記メモリに記憶され、前記１つ又は複数のプログラムが前記プロセッサにより実行されると、前記電子機器に請求項１ないし７のいずれかに記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法を実行させる
ことを特徴とする電子機器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はイカ資源量の予測の技術分野に属し、南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法及びその応用に関する。
【背景技術】
【０００２】
アルゼンチンマツイカは世界で重要な経済的頭足類資源であり、史上最高の漁獲量は１００万トンを超え、中国等の国や地域のイカ釣り漁船、トロール船等の主な漁獲対象となっている。アルゼンチンマツイカは南西大西洋海域に分布しており、ブラジル海流及びフォークランド海流の影響を受け、南北方向に回遊する。既存の研究からわかるように、アルゼンチンイレックスは一年生の種であるため、その資源補充量は産卵場の環境の影響を受けており、汪金濤ら（２０１４）は２００３−２０１１年の中国のイカ釣り漁船隊の南西大西洋での生産統計データ、及び産卵場の海面水温（ＳＳＴ）、海面水温偏差（ＳＳＴＡ）に基づいて、アルゼンチンイレックスの産卵期における月別の産卵場の最適海面水温範囲の、総面積に対する比率（ＰＳで示される）、及び海流強度を表すＳＳＴ、ＳＳＴＡ等の様々な環境変数因子と単位努力量当たりの漁獲量（ＣＰＵＥ）との相関性を計算分析し、主要環境因子に基づく様々な資源補充量予測モデルを作成した。伍玉梅ら（２０１１）は２０００−２００８年の中国のイカ釣り漁船の南西大西洋でのアルゼンチンイレックスの生産データ及び海洋衛星インバージョンデータを利用して、過去９年間の南西大西洋のアルゼンチンイレックス資源豊度の変化、及びそれと主要な生態学的因子（海面水温やクロロフィルａ）との関係を分析し、２００４−２００８年では年平均海面水温と資源豊度の変動が有意な負の相関を示し、２００５−２００８年では年平均クロロフィルａと資源豊度が高い正の相関を示すことを見出した。Ｗａｌｕｄａ及びＲｏｄｈｏｕｓｅ（２００１）は、産卵場の環境条件がアルゼンチンマツイカの資源補充量に大きな影響を与え、産卵場の水温がやや高く、産卵場の適切な海面水温の面積が高いと、アルゼンチンマツイカの発生及び成長に有利であるという見解を示した。上記研究からわかるように、現在、中国国内外の各学者はアルゼンチンマツイカの産卵場の環境による資源補充量への影響を鋭意研究して、対応する資源量予測モデルを作成したが、どのようにして気候因子を利用してその資源量を事前に予測するかは未だに空白のままである。
【０００３】
従って、気候因子に基づいてアルゼンチンマツイカ資源量を予測する方法を開発することは、非常に現実的である。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明は、従来技術に存在している問題及び欠陥を克服し、南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法及びその応用を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記目的を実現するために、本発明は以下の技術的解決手段を提供し、
電子機器に適用される南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法であって、
（１）過去Ｎ年間にアルゼンチンマツイカが分布していた南西太平洋の南部海域の月別の南極振動指数ＡＡＯ値を取得するステップと、
（２）時系列分析方法を利用して、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去Ｎ年間の月別のＡＡＯ値に対して相関性分析を行い、そのうち統計的に有意な相関を示す月のＡＡＯ値を選択し、選択したこれらの月のＡＡＯ値を、アルゼンチンマツイカ資源豊度に影響を与える気候因子とし、選択したこれらの月のＡＡＯ値を１、２、３…ｚ…ｍの順で番号付けし、これらの月のＡＡＯ値を順にｘ１、ｘ２、ｘ３…ｘｚ…ｘｍ（ｍは選択した月のＡＡＯ値の数である）とするステップと、
（３）ｘ１、ｘ２、ｘ３…ｘｍのうちの任意１〜ｍ個の気候因子に対して、多変数線形方程式を利用して計２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルを作成し、各予測モデルの統計的Ｐ値を計算し、アルゼンチンマツイカ資源豊度予測モデルの式は、
ｌｎ（ＣＰＵＥ）＝ａ＋ｂ１＊ｘ１＋ｂ２＊ｘ２＋ｂ３＊ｘ３＋…＋ｂｚ＊ｘｚ＋…＋ｂｍ＊ｘｍであり、
式中、ＣＰＵＥは１隻あたりの年間漁獲量、ａは定数、ｂ１、ｂ２、ｂ３、…ｂｚ…、ｂｎはそれぞれｘ１、ｘ２、ｘ３…ｘｚ…ｘｍに対応する係数であり、
2021520575.tif 000003
（４）上記２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルのうち、統計的Ｐ値が最も小さいモデルを最適モデルとして選択するステップと、を含む。
【０００６】
本発明に係る南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法によれば、気候因子（南極振動指数）とアルゼンチンマツイカ資源量の予測との関係を初めて作成し、気候因子（南極振動指数）によってアルゼンチンマツイカ資源量の迅速かつ正確な予測を実現でき、海洋漁業生産（アルゼンチンマツイカの漁獲）では優れた指導の役割を果たし、漁獲効率を大幅に向上させ、漁獲コストを低減させ、将来性が期待でき、また、本発明によって得られる最適モデルは一律ではなく、リアルタイムに取得される最新データに応じて最適モデルを再取得することができ、本発明の方法は適応性が高く、将来性が期待できる。
【０００７】
好適な技術的解決手段としては、
前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法において、最適モデルを得た後、該最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得して最適モデルに入力し、最適モデルがアルゼンチンマツイカ資源豊度を出力すると、アルゼンチンマツイカ資源量の予測が完了する。
【０００８】
前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法において、アルゼンチンマツイカ資源豊度予測モデルが気候因子ｘｚに対して作成される場合、ｂ１、ｂ２、ｂ３、…ｂｚ−１、ｂｚ＋１、…、ｂｎはいずれも０であり、以下同様に、アルゼンチンマツイカ資源豊度予測モデルが気候因子ｘ１及びｘｚに対して作成される場合、ｂ２、ｂ３、…ｂｚ−１、ｂｚ＋１、…、ｂｎはいずれも０である。
【０００９】
前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法において、前記統計的に有意な相関とは、算出されたＰ値＜０．０５である。Ｐ値は、統計では仮説検定の結果を判定するためのパラメータであり、Ｐ値（Ｐ  ｖａｌｕｅ）は、帰無仮説を真とするときに得られる標本観察結果又はより極端な結果が発生する確率である。Ｐ値が非常に小さいと、帰無仮説状況の発生する確率が非常に小さいことを示し、発生したら、小確率原理に従って、帰無仮説を拒否する理由があり、Ｐ値が小さいほど、帰無仮説を拒否する理由が十分になる。要するに、Ｐ値が小さいほど、結果が有意であることを示す。
【００１０】
前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法において、ステップ（２）では、統計的に有意な相関を示す２つの月のＡＡＯ値、すなわち過去２年目の５月のＡＡＯ値及び過去１年目の８月のＡＡＯ値を選択し、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去２年目の５月のＡＡＯ値との相関性は有意であり、正の相関を示し、相関係数がそれぞれ０．５４９であり、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去１年目の８月のＡＡＯ値との相関性は有意であり、負の相関を示し、相関係数が−０．４９７である。ここ及び以下で言及される３個の予測モデルは単にデータの一部で本発明の予測方法の動作ロジックを示すためのものであり、本発明の保護範囲はこれに限定されず、当業者は実際の必要に応じて適切なデータを選択してアルゼンチンマツイカ資源量を予測することができ、選択される統計的に有意な相関を示す月のＡＡＯ値の数は２個に限定されず、予測モデルの数は選択される統計的に有意な相関を示す月のＡＡＯ値の数に応じて変化する。
【００１１】
前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法において、ステップ（３）では、以下の計３個の予測モデルを作成し、
１）予測モデルＩ
気候因子の過去２年目の５月のＡＡＯ値ｘ１に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７３０５＋０．４５５４＊ｘ１であり、
Ｆ値は６．０３４７であり、Ｐ＝０．０２７７＜０．０５であり、Ｆ値はＦ検定の統計値であり、Ｆ検定は帰無仮説（ｎｕｌｌ  ｈｙｐｏｔｈｅｓｉｓ、Ｈ０）の下で統計値がＦ−分布に従う検定であり、通常、２個以上のパラメータを用いた統計モデルを分析して、該モデル中の全部又は一部のパラメータが母集団の推定に適用できるか否かを判断するために使用され、線形関係を検定するものであり、
２）予測モデルＩＩ
気候因子の過去１年目の８月のＡＡＯ値ｘ２に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７２８１−０．４６９４＊ｘ２であり、
Ｆ値は４．５４８０であり、Ｐ＝０．０４９＜０．０５であり、
３）予測モデルＩＩＩ
気候因子の過去２年目の５月のＡＡＯ値ｘ１及び過去１年目の８月のＡＡＯ値ｘ２に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７０４６＋０．３９０８＊ｘ１−０．３８３６＊ｘ２であり、
Ｆ値は５．５１３５であり、Ｐ＝０．０１８＜０．０５である。
【００１２】
前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法において、ステップ（４）では、予測モデルＩＩＩを最適モデルとして選択し、最適モデルは、
ｌｎ  （ＣＰＵＥ）＝６．７０４６＋０．３９０８＊ｘ１−０．３８３６＊ｘ２である。
【００１３】
本分野の常識をもとに上記各好適な条件を任意に組み合わせることができ、それにより本発明の各好適な例が得られる。
【００１４】
本発明は電子機器をさらに提供し、１つ又は複数のプロセッサ、１つ又は複数のメモリ、１つ又は複数のプログラム、及びデータ収集装置を備え、
前記データ収集装置は最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得するために使用され、前記１つ又は複数のプログラムは前記メモリに記憶され、前記１つ又は複数のプログラムが前記プロセッサにより実行されると、前記電子機器に前記南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法を実行させる。
【発明の効果】
【００１５】
有益な効果については、
（１）本発明に係る南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法によれば、南極振動ＡＡＯを利用してアルゼンチンマツイカ資源豊度の予測を実現することで、海洋漁業生産（アルゼンチンマツイカの漁獲）では優れた指導の役割を果たし、漁獲効率を大幅に向上させ、漁獲コストを低減させ、将来性が期待できる。
【００１６】
（２）本発明に係る南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法によれば、最適モデルは一律ではなく、リアルタイムに取得される最新データに応じて最適モデルを再取得することができ、本発明の方法は適応性が高く、将来性が期待できる。
【００１７】
（３）本発明に係る電子機器は、構成が簡単で、コストが低く、南極振動ＡＡＯに基づいてアルゼンチンマツイカ資源豊度を迅速に予測でき、将来性が期待できる。
【図面の簡単な説明】
【００１８】
【図１】図１は本発明に係る南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法のフローチャートである。
【図２】図２は２００２−２０１７年のアルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）の年間変化図である。
【図３】図３及び図４はアルゼンチンマツイカ資源豊度予測モデルの分布図である。
【図４】図３及び図４はアルゼンチンマツイカ資源豊度予測モデルの分布図である。
【図５】図５及び図６はアルゼンチンマツイカ資源豊度予測モデルの分布図である。
【図６】図５及び図６はアルゼンチンマツイカ資源豊度予測モデルの分布図である。
【図７】図７は２００２−２０１７年のアルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）の実際値及び予測値の変化分布図である。
【図８】図８は本発明に係る電子機器の構造模式図である。
【発明を実施するための形態】
【００１９】
本発明の実施例の目的、技術的解決手段及び利点をより明確にするために、以下、本発明の実施例の図面を参照しながら本発明の実施例の技術的解決手段を明確かつ完全に説明し、明らかなように、説明される実施例は本発明の一部の実施例であり、全部の実施例ではない。本発明の実施例に基づいて、当業者が創造的な努力をせずに想到し得るほかの実施例はすべて本発明の保護範囲に属する。
【００２０】
実施例１
【００２１】
図１に示すように、本実施例は南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法を提供し、
過去Ｎ年間にアルゼンチンマツイカが分布していた南西太平洋の南部海域の月別の南極振動指数ＡＡＯ値を取得するステップ１０１と、
時系列分析方法を利用して、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と過去Ｎ年間の月別のＡＡＯ値に対して相関性分析を行い、そのうち統計的に有意な相関（Ｐ値＜０．０５）を示す月のＡＡＯ値を選択し、選択したこれらの月のＡＡＯ値を、アルゼンチンマツイカ資源豊度に影響を与える気候因子とし、選択したこれらの月のＡＡＯ値を１、２、３…ｚ…ｍの順で番号付けし、これらの月のＡＡＯ値を順にｘ１、ｘ２、ｘ３…ｘｚ…ｘｍ（ｍは選択した月のＡＡＯ値の数である）とするステップ１０２と、
ｘ１、ｘ２、ｘ３…ｘｍのうちの任意１〜ｍ個の気候因子に対して、多変数線形方程式を利用して計２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルを作成し、各予測モデルの統計的Ｐ値を計算し、アルゼンチンマツイカ資源豊度予測モデルの式は、
ｌｎ（ＣＰＵＥ）＝ａ＋ｂ１＊ｘ１＋ｂ２＊ｘ２＋ｂ３＊ｘ３＋…＋ｂｚ＊ｘｚ＋…＋ｂｍ＊ｘｍであり、
式中、ＣＰＵＥは１隻あたりの年間漁獲量、ａは定数、ｂ１、ｂ２、ｂ３、…ｂｚ…、ｂｎはそれぞれｘ１、ｘ２、ｘ３…ｘｚ…ｘｍに対応する係数であるステップ１０３と、
上記２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルのうち、統計的Ｐ値が最も小さいモデルを最適モデルとして選択するステップ１０４と、
該最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得して最適モデルに入力し、最適モデルがアルゼンチンマツイカ資源豊度を出力すると、アルゼンチンマツイカ資源量の予測が完了するステップ１０５と、を含む。
【００２２】
当業者が本発明の技術的解決手段をよりよく理解できるために、以下、具体例を挙げて本発明を説明する。
【００２３】
１、資料及び方法
【００２４】
（１）データソース
【００２５】
アルゼンチンマツイカは南西太平洋の南部海域に分布し、南極海域に近接し、その産卵場及び索餌場が南極海域の西風海流によるフォークランド海流の影響を受ける。該海域では、南極振動（Ａｎｔａｒｃｔｉｃ  Ｏｓｃｉｌｌａｔｉｏｎ、ＡＡＯ）は重要な気候因子の一つである。南極振動は南半球の大気循環の主なモードであり、複数の尺度で南半球及び北半球の部分地域の気候系に重要な影響を与えている。それは、南半球の中緯度及び高緯度の２つの大気環状活動帯間で大気の質量が変化する全地球規模の「シーソー」構造であり、４０°Ｓ及び７０°Ｓ上の標準化された東西平均海面気圧差を、南極振動変化を測定するための指数とすることができる。南極振動が強いと、南半球の周極低気圧が強くなり、中高緯度の西風が強くなり、その逆も同じである。
【００２６】
南極振動指数はアメリカ海洋大気庁のウェブサイト（ｈｔｔｐ：／／ｗｗｗ．ｃｐｃ．ｎｃｅｐ．ｎｏａａ．ｇｏｖ／ｐｒｏｄｕｃｔｓ／ｐｒｅｃｉｐ／ＣＷｌｉｎｋ／ｄａｉｌｙ＿ａｏ＿ｉｎｄｅｘ／ａａｏ／ｍｏｎｔｈｌｙ．ａａｏ．ｉｎｄｅｘ．ｂ７９．ｃｕｒｒｅｎｔ．ａｓｃｉｉ．ｔａｂｌｅ）からであり、期間が２０００年１月−２０１７年１２月（表１）である。
アルゼンチンマツイカ資源豊度指数ＣＰＵＥ（単位はトン／船）は中国遠洋イカ釣り漁船の年間漁獲量からであり、期間が２００２年−２０１７年（表２）である。
【００２７】
【表１】
2021520575.tif 000004
【００２８】
【表２】
2021520575.tif 000005
【００２９】
（２）研究方法及びステップ
【００３０】
アルゼンチンマツイカ資源豊度はイカ釣り漁船の年間平均漁獲量を指標とし、該指標は生産統計等の要因によって誤差が生じるため、資源豊度指数に対して自然対数ｌｎ（ＣＰＵＥ）を取ることで標準化を行う。
【００３１】
時系列分析方法を利用して、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）と２０００−２０１７年の１−１２月のＡＡＯ値に対して相関性分析を行い、そのうち統計的に有意な相関（Ｐ値＜０．０５）を示す月のＡＡＯ値を選択し、選択したこれらの月のＡＡＯ値を、アルゼンチンマツイカ資源豊度に影響を与える気候因子とし、選択したこれらの月のＡＡＯ値を１、２、３…ｚ…ｍの順で番号付けし、これらの月のＡＡＯ値を順にｘ１、ｘ２、ｘ３…ｘｚ…ｘｍ（ｍは選択した月のＡＡＯ値の数である）とし、
ｘ１、ｘ２、ｘ３…ｘｍのうちの任意１〜ｍ個の気候因子に対して、多変数線形方程式を利用して計２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルを作成し、各予測モデルの統計的Ｐ値を計算し、アルゼンチンマツイカ資源豊度予測モデルの式は、
ｌｎ（ＣＰＵＥ）＝ａ＋ｂ１＊ｘ１＋ｂ２＊ｘ２＋ｂ３＊ｘ３＋…＋ｂｚ＊ｘｚ＋…＋ｂｍ＊ｘｍであり、
式中、ＣＰＵＥは１隻あたりの年間漁獲量、ａは定数、ｂ１、ｂ２、ｂ３、…ｂｚ…、ｂｎはそれぞれｘ１、ｘ２、ｘ３…ｘｚ…ｘｍに対応する係数であり、
上記２ｍ−１個のアルゼンチンマツイカ資源豊度予測モデルのうち、統計的Ｐ値が最も小さいモデルを最適モデルとして選択し、
該最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得して最適モデルに入力し、最適モデルがアルゼンチンマツイカ資源豊度を出力すると、アルゼンチンマツイカ資源量の予測が完了する。
【００３２】
２、研究結果
【００３３】
（１）年間資源豊度ｌｎ（ＣＰＵＥ）の変化
【００３４】
図２からわかるように、アルゼンチンマツイカ資源豊度ｌｎ（ＣＰＵＥ）は大きな年間変化を示し、２００２−２００３年、２００７−２００８年、２０１４−２０１５年では高い資源量レベルであったが、２００５−２００６年、２０１０−２０１２年及び２０１６−２０１７年では低い資源量レベルであった。
【００３５】
（２）資源豊度ｌｎ（ＣＰＵＥ）に影響を与えるＡＡＯ値
【００３６】
資源豊度ｌｎ（ＣＰＵＥ）と過去２年目の月別のＡＡＯ値との相関性分析を行った結果、資源豊度ｌｎ（ＣＰＵＥ）と過去２年目の５月のＡＡＯ値との相関性は有意であり、正の相関を示し、相関係数が０．５４９（Ｐ＜０．０５）であり、
資源豊度ｌｎ（ＣＰＵＥ）と過去１年目の月別のＡＡＯ値との相関性分析を行った結果、資源豊度ｌｎ（ＣＰＵＥ）と過去１年目の８月のＡＡＯ値との相関性は有意であり、負の相関を示し、相関係数が−０．４９７（Ｐ＜０．０５）である。
【００３７】
（３）資源豊度予測モデルの作成
【００３８】
１）予測モデルＩ
気候因子の過去２年目の５月のＡＡＯ値ｘ１に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７３０５＋０．４５５４＊ｘ１であり、
Ｆ値は６．０３４７であり、Ｐ＝０．０２７７＜０．０５であり、
その実際値と予測値の統計表は表３に示され、モデルの分布図は図３（図３中のＡＯＯｔ−２はｘ１である）に示される。
【００３９】
【表３】
2021520575.tif 000006
【００４０】
２）予測モデルＩＩ
気候因子の過去１年目の８月のＡＡＯ値ｘ２に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７２８１−０．４６９４＊ｘ２であり、
Ｆ値は４．５４８０であり、Ｐ＝０．０４９＜０．０５であり、
その実際値と予測値の統計表は表４に示され、モデルの分布図は図４（図４中のＡＯＯｔ−１はｘ２である）に示される。
【００４１】
【表４】
2021520575.tif 000007
【００４２】
３）予測モデルＩＩＩ
気候因子の過去２年目の５月のＡＡＯ値ｘ１及び過去１年目の８月のＡＡＯ値ｘ２に対してアルゼンチンマツイカ資源豊度予測モデルを作成し、具体的には、
ｌｎ  （ＣＰＵＥ）＝６．７０４６＋０．３９０８＊ｘ１−０．３８３６＊ｘ２であり、
Ｆ値は５．５１３５であり、Ｐ＝０．０１８＜０．０５であり、
その実際値と予測値の統計表は表５に示され、モデルの分布図は図５及び図６（図５中のＡＯＯｔ−２はｘ１、図６中のＡＯＯｔ−１はｘ２である）に示される。
【００４３】
【表５】
2021520575.tif 000008
【００４４】
上記３個のモデルを比較分析した結果、予測モデルＩＩＩを最適モデルとして選択し、最適モデルはｌｎ  （ＣＰＵＥ）＝６．７０４６＋０．３９０８＊ｘ１−０．３８３６＊ｘ２である。２００２〜２０１７年の結果に対応するｘ１及びｘ２を最適モデルに入力して予測値を得て（２００４年の結果を例にすると、ｘ１は２００２年５月のＡＡＯ値、ｘ２は２００３年８月のＡＡＯ値である）、その実際値と予測値の資源豊度の変化傾向は図７に示され（図中の正方形は予測値、斜方形は実際値である）、図７からわかるように、予測値と実際値の変化傾向はほぼ同じであり、つまり、本発明の方法によって、アルゼンチンマツイカ資源豊度を効果的に予測できる。
【００４５】
実施例２
【００４６】
電子機器であって、図８に示すように、１つ又は複数のプロセッサ、１つ又は複数のメモリ、１つ又は複数のプログラム、及びデータ収集装置を備え、
データ収集装置は最適モデルに対応するｘ１、ｘ２、ｘ３…ｘｚ…ｘｍを取得するために使用され、１つ又は複数のプログラムはメモリに記憶され、１つ又は複数のプログラムがプロセッサにより実行されると、電子機器に実施例１に記載の南極振動指数に基づくアルゼンチンマツイカ資源量の予測方法を実行させる。
【００４７】
以上、本発明の具体的な実施形態を説明したが、当業者であれば、これらは単に例示的な説明であり、本発明の保護範囲は添付の特許請求の範囲に定められると理解できる。当業者は本発明の原理及び趣旨を逸脱せずに、これらの実施形態に対して種々の変更や修正を行うことができ、これらの変更や修正はすべて本発明の保護範囲に属する。
【図１】
2021520575.tif 000009
【図２】
2021520575.tif 000010
【図３】
2021520575.tif 000011
【図４】
2021520575.tif 000012
【図５】
2021520575.tif 000013
【図６】
2021520575.tif 000014
【図７】
2021520575.tif 000015
【図８】
2021520575.tif 000016
【国際調査報告】
2021520575.tif 000017
2021520575.tif 000018
2021520575.tif 000019
2021520575.tif 000020