(57)【要約】【課題】任意の地点の紫外線強度を高精度かつ細密に予測する。【解決手段】紫外線強度予測部１６は、気象予測モデル演算部１４から雲量の影響を考慮した全天日射量の予測値を取得する。変換式作成部１６２は、予め取得された紫外線強度観測値と全天日射量観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する。波長別紫外線強度算出部１６３は、上記作成された波長毎の変換式に基づいて、任意の地点の緯度とオゾン全量取得部１６１により得られるオゾン全量と、気象予測モデル演算部１４から得られた全天日射量の予測値とをもとに波長別の紫外線強度の予測値を算出する。【選択図】      図２
2008089425.tif 000002
【特許請求の範囲】
【請求項１】
任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得する取得手段と、
予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する作成手段と、
前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出する算出手段と
を具備することを特徴とする紫外線強度予測装置。
【請求項２】
前記全天日射量の予測値は、前記地点の雲量の予測値から求められることをさらに特徴とする請求項１記載の紫外線強度予測装置。
【請求項３】
任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得する取得手段と、
予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する作成手段と、
前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出する算出手段と、
前記算出手段により算出された波長別の紫外線強度をＵＶインデックスに変換する変換手段と
を具備することを特徴とするＵＶインデックス予測装置。
【請求項４】
前記全天日射量の予測値は、前記地点の雲量の予測値から求められることをさらに特徴とする請求項３記載のＵＶインデックス予測装置。
【請求項５】
任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得し、
予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成し、
前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出することを特徴とする紫外線強度予測方法。
【請求項６】
前記全天日射量の予測値は、前記地点の雲量の予測値から求められることをさらに特徴とする請求項５記載の紫外線強度予測方法。
【請求項７】
任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得し、
予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成し、
前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出し、
前記算出された波長別の紫外線強度をＵＶインデックスに変換することを特徴とするＵＶインデックス予測方法。
【請求項８】
前記全天日射量の予測値は、前記地点の雲量の予測値から求められることをさらに特徴とする請求項７記載のＵＶインデックス予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、紫外線強度を予測する紫外線強度予測装置、ＵＶインデックス予測装置、紫外線強度予測方法、及びＵＶインデックス予測方法に関する。
【背景技術】
【０００２】
近年、オゾン層の破壊による地上の紫外線量の増加が人体や一般生物に有害な影響を及ぼすとし、深刻な地球環境問題として取り上げられている。人間や生物にとって太陽光は不可欠である。さらに紫外線は、人間の皮膚内でのビタミンＤの生成や強い殺菌作用など有用な面も少なくない。しかしながら、最近の医学、光生物学における研究により、太陽光とりわけ紫外線の生物に与える有害性の大きさが明らかになっている。紫外線による有害性は、光毒性反応、光感作と光アレルギー、細胞障害、また、皮膚の悪性腫瘍などが挙げられる。オゾン全量の１０％が減少すると、有害なＵＶ−Ｂは世界的にみて年平均約２０％増加すると言われており、紫外線量の増加が懸念される現在において、その高精度な予測技術開発は火急の課題と言える。
【０００３】
既存の紫外線予報として代表とされるのが、気象庁が２００５年５月より開始したＵＶインデックスによる紫外線予報である。ＵＶインデックス（ＵＶ指数）とは、紫外線が人体に及ぼす影響の度合いをわかりやすく示すために、紫外線の強さを指標化したもので、紫外線強度を１３段階にインデックス化した全国の紫外線分布情報が提供されている。気象庁の予測手法では、地上の紫外線量に大きく影響するオゾン量を求めるために、化学輸送モデルを用いて、オゾンに関わる光化学反応を計算している。そこから得られるオゾン量の分布や累年平均エアロゾル観測値を入力値として、放射伝達モデルから２０ｋｍ格子間隔の晴天時紫外線量を予測している。そして最終的には、雲による紫外線減衰効果を考慮するために、過去の観測値と天気（降水、雲天等）の関係から補正係数を算出して、晴天時紫外線量に乗じることで日々の紫外線量を予測している（例えば、非特許文献１又は２を参照。）。
【非特許文献１】「気象庁の紫外線情報について」、安達他、２００６年度日本気象学会春季大会予稿集、ｐ．２６１
【非特許文献２】「紫外線情報の解説」、気象庁、インターネット＜ＵＲＬ：http://www.data.kishou.go.jp/obs-env/uvhp/3-47uvindex_info.html＞
【非特許文献３】「ＵＶインデックスを求めるには」、気象庁、インターネット＜ＵＲＬ：http://www.data.kishou.go.jp/obs-env/uvhp/3-41uvindex_define.html＞
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところが、上述した気象庁の予測手法は極めて理論的な手法といえるが、その反面、放射伝達モデルの計算時間が膨大になるのを避けるため、紫外線量の減衰にもっとも起因する雲量の効果を、天気の区分という大まかな指標から算出している点に問題がある。
【０００５】
この発明は上記事情に着目してなされたもので、その目的とするところは、任意の地点の紫外線強度を高精度かつ細密に予測可能な紫外線強度予測装置、ＵＶインデックス予測装置、紫外線強度予測方法、及びＵＶインデックス予測方法を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的を達成するためにこの発明に係わる紫外線強度予測装置は、任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得する取得手段と、予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する作成手段と、前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出する算出手段とを具備することを特徴とする。
【０００７】
この発明に係わるＵＶインデックス予測装置は、任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得する取得手段と、予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する作成手段と、前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出する算出手段と、前記算出手段により算出された波長別の紫外線強度をＵＶインデックスに変換する変換手段と
を具備することを特徴とする。
【０００８】
また、この発明に係わる紫外線強度予測方法は、任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得し、予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成し、前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出することを特徴とする。
【０００９】
この発明に係わるＵＶインデックス予測方法は、任意の地点の全天日射量の予測値、緯度及びオゾン全量を取得し、予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成し、前記作成された波長毎の変換式に基づいて、前記取得された緯度及びオゾン全量と前記全天日射量の予測値とをもとに波長別の紫外線強度を算出し、前記算出された波長別の紫外線強度をＵＶインデックスに変換することを特徴とする。
【００１０】
上記構成による紫外線強度予測装置、ＵＶインデックス予測装置、紫外線強度予測方法、及びＵＶインデックス予測方法では、気象予測モデルによる雲量の影響を考慮した全天日射量の予測値を用いて紫外線強度を求めることで、従来に比べて高精度かつ細密な紫外線強度を求めることが可能となる。
【発明の効果】
【００１１】
したがってこの発明によれば、任意の地点の紫外線強度を高精度かつ細密に予測可能な気象予測システム及び紫外線強度予測方法を提供することができる。
【発明を実施するための最良の形態】
【００１２】
以下、図面を参照しながら本発明の実施の形態を詳細に説明する。  
図１は、この発明に係わる気象予測システムの一実施形態の構成を示すブロック図である。この気象予測システム１は、ネットワークＮＴを介して気象庁データサーバＤＳ０、及びレーダサイトサーバＤＳ１，ＤＳ２に接続されている。気象予測システム１は、ネットワークＮＴと接続される通信インターフェース１２と、通信処理部１１と、観測データ格納部１３と、気象予測モデル演算部１４と、予測データ格納部１５と、紫外線強度予測部１６と、ＵＶインデックス変換部１７とを備える。
【００１３】
通信処理部１１は、気象庁データサーバＤＳ０やレーダサイトサーバＤＳ１，ＤＳ２から気象予測のもとになる観測データ・予測データ（ＧＰＶデータ）をネットワークＮＴを介して入手する。通信処理部１１で入手された気象観測データは、観測データ格納部１３に格納され、気象予測モデル演算部１４からの要求に応じて選択的に気象予測モデル演算部１４に送られる。
【００１４】
気象予測モデル演算部１４は、気象予測のもととなるデータが観測データ格納部１３に格納されると起動し、気象予測演算を行う。求められた気象予測データは予測データ格納部１５に記憶される。また、観測データ格納部１３に新たな観測データが入力されると、気象予測モデル演算部１４は再び起動し、観測と予測のズレを補正するために気象予測演算を再実行する。本実施形態では、気象予測モデル演算部１４は、雲量を予測し、予測された雲量をもとに予測計算される全天日射量を出力する。
【００１５】
紫外線強度予測部１６は、気象予測モデル演算部１４から出力される全天日射量の予測データを用いて紫外線強度の予測演算を行い、紫外線強度の予測値を出力する。ＵＶインデックス変換部１７は、紫外線強度予測部１６で予測演算された紫外線強度の予測値をＵＶインデックスに変換して出力する。
【００１６】
図２は、紫外線強度予測部１６の機能ブロック図である。紫外線強度予測部１６は、オゾン全量取得部１６１と、変換式作成部１６２と、波長別紫外線強度算出部１６３とを備える。オゾン全量取得部１６１は、各観測ポイントにおけるオゾン情報（オゾンポイント情報）をもとに内挿処理を行うことで、任意の地点のオゾン全量を求める。変換式作成部１６２は、予め取得された全天日射量観測値と紫外線強度観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する。波長別紫外線強度算出部１６３は、上記変換式作成部１６２において作成された波長毎の変換式に基づいて、任意の地点の緯度及びオゾン全量と気象予測モデル演算部１４で予測される全天日射量とをもとに波長毎の紫外線強度を算出する。
【００１７】
次に、このように構成される紫外線強度予測部１６の動作について説明する。図３は、紫外線強度予測部１６の処理手順とその処理内容を示すフローチャートである。  
図３において、紫外線強度予測部１６は、気象予測モデル演算部１４で予測される雲量の影響を考慮した全天日射量を取得する（ステップＳ３ａ）。変換式作成部１６２は、全天日射量を紫外線強度に変換する波長毎の変換式を作成する（ステップＳ３ｂ）。
【００１８】
ここでは、紫外線観測データとして、例えば、気象業務支援センターから配布される「大気・海洋環境観測報告(ＣＤ−ＲＯＭ)」を用いる。２００６年３月現在における最新のもの（第５号平成１５年観測成果）には、５地点における１時間毎の紫外線域日射量が収録されている。収録期間を図４に示す。なお、同ＣＤ−ＲＯＭには、地表に到達するＵＶ−Ｂに影響を及ぼすオゾン全量も収録されている。オゾン全量とは、上空に分布するオゾンを１気圧に圧縮したときの厚さを表す。
【００１９】
図５は、大気・海洋環境観測報告に含まれる２００３年８月のつくばにおける紫外線強度(ＵＶ−Ａ及びＢ)と全天日射量の相関を示す。ＵＶ−Ａ、Ｂともに相関係数が０．９以上であり、全天日射量に対して非常に高い相関を示している。紫外線強度は、基本的には全天日射量と連動しており、日本においては全天日射量の４〜７％となっている。つまり、全天日射量の予測が可能であれば、回帰的に紫外線強度の予測は可能である。このことから、紫外線強度と全天日射量との間には、以下のような回帰関係が成立する。  
［紫外線強度］＝Ｋ×[全天日射量]
よって、この回帰係数Ｋを如何に求めるかが、紫外線強度の予測を行う上で鍵となる。係数Ｋの決定の際には、ＵＶ−Ａ及びＢそれぞれの特徴を踏まえ、以下の点を考慮する。  
ＵＶ−Ａ：緯度を説明変数とした係数Ｋを求める  
ＵＶ−Ｂ：緯度とオゾン全量を説明変数とした係数Ｋを求める
変換式作成部１６２では、図５と同様に各観測点での波長別の回帰係数Ｋを算出する。そして、その結果をもとに任意の地点の係数Ｋを求めるための変換式を作成する。変換式として、以下の関係が成り立つものと仮定した。  
Ｋ＝ａＸ１＋ｂＸ２＋ｃ  
ここで、係数Ｋは、先に求めた各観測点における波長別の回帰係数である。Ｘ１及びＸ２はそれぞれ、緯度およびオゾン全量（参照値）である。参照値とは、例えば１９７１〜２０００年の月平均値で、平均的なオゾン量である。これにより波長毎のａ、ｂ、およびｃの値が決定し、緯度Ｘ１及びオゾン全量Ｘ２を変数とする波長毎の変換式が作成される。
【００２０】
この作成された波長毎の変換式を用いて、波長別紫外線強度算出部１６３は、波長別の紫外線強度を算出する（ステップＳ３ｃ）。波長別紫外線強度算出部１６３は、波長毎の変換式に基づいて、任意の地点の緯度及びオゾン全量取得部１６１により取得されるオゾン全量をもとに波長毎の回帰係数Ｋの値を求め、求められた波長毎の回帰係数Ｋと気象予測モデル演算部１４から得られる全天日射量の予測値とにより波長別の紫外線強度を算出する。
【００２１】
さらに、ＵＶインデックス変換部１７は、上記算出された波長別紫外線強度にＣＩＥ作用スペクトルを乗算し、波長別紅斑紫外線強度を求める（ステップＳ３ｄ）。図６に、波長別赤斑紫外線強度の計算方法を示す（非特許文献３を参考）。紫外線の生物への影響は、国際照明委員会（ＣＩＥ：Commission International de l'Eclairage）が定めるＣＩＥ作用スペクトルを波長別紫外線強度に乗算し、重み付けを施すことで試算される。このように求められた波長別強度は、波長別紅斑紫外線強度と呼ばれる。そして、この波長別紅斑紫外線強度を２８０〜４００ｎｍにわたって波長積分することにより、生物への総合的な影響度が得られる（地上での２９０ｎｍ以下の紫外線は実質的に０と見なせるので、実際には２９０〜４００ｎｍにわたる波長積分で十分に精度よく求められる）。その積分された値を２５ｍＷ／ｍ２で割ったものがＵＶインデックスと呼ばれるものである。
【００２２】
ＣＩＥ作用スペクトルβ（λ）は、以下の式で表される。
【００２３】
β（λ）＝１．０    ｆｏｒ  λ＜２９８ｎｍ，
１００．０９４（２９８−λ）  ｆｏｒ  ２９８＜λ＜３２８ｎｍ，
１００．０１５（１３９−λ）  ｆｏｒ  ３２８＜λ＜４００ｎｍ，
０                  ｆｏｒ  λ＞４００ｎｍ
なお、大気・海洋環境観測報告に３２５〜４００ｎｍの紫外線観測値が含まれていないため、当領域の紫外線強度を一定として扱い、３２５ｎｍの直前の値を代用する。生物への影響度を考慮する際に、波長別紫外線強度にＣＩＥ作用スペクトルをかけると、３２５〜４００ｎｍの波長別紅斑紫外線強度は、非常に小さな値になるため、一定として扱っても全体へ与える影響は非常に小さい。
【００２４】
ＵＶインデックス変換部１７は、上記求められた波長別紅斑紫外線強度を積分して、２５ｍＷ／ｍ２で除算し（ステップＳ３ｅ）、この値をＵＶインデックスとして出力する（ステップＳ３ｆ）。
【００２５】
以上述べたように、上記実施形態では、紫外線強度予測部１６は、気象予測モデル演算部１４から雲量の影響を考慮した全天日射量の予測値を取得する。変換式作成部１６２は、予め取得された紫外線強度観測値と全天日射量観測値との相関関係から全天日射量を紫外線強度に変換する波長毎の変換式を作成する。波長別紫外線強度算出部１６３は、上記作成された波長毎の変換式に基づいて、任意の地点の緯度とオゾン全量取得部１６１により得られるオゾン全量と、気象予測モデル演算部１４から得られた全天日射量の予測値とをもとに波長別の紫外線強度の予測値を算出する。
【００２６】
したがって上記実施形態によれば、気象予測モデルによる雲量の影響を考慮した全天日射量の予測値を用いて紫外線強度を求めることで、従来に比べて高精度かつ細密な紫外線強度を求めることができる。これにより、野外での社会活動を行うにあたって的確なＵＶ指標を提示することが可能となり、有害紫外線に起因する皮膚がんや白内障などの健康被害の予防につながる。
【００２７】
要するにこの発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。
【図面の簡単な説明】
【００２８】
【図１】この発明に係わる気象予測システムの一実施形態を示すブロック構成図。
【図２】図1に示す紫外線強度予測部の構成を示すブロック図。
【図３】図２に示す紫外線強度予測部の処理手順とその内容を示すフローチャート。
【図４】大気・海洋環境観測報告に収録される紫外線域日射量及びオゾン全量の収録期間を示す図。
【図５】紫外線強度（ＵＶ−Ａ及びＢ）と全天日射量の相関を示す図。
【図６】波長別赤斑紫外線強度の計算方法を示す図。
【符号の説明】
【００２９】
１…気象予測システム、１１…通信処理部、１２…通信インターフェース、１３…観測データ格納部、１４…気象予測モデル演算部、１５…予測データ格納部、１６…紫外線強度予測部、１７…ＵＶインデックス変換部、ＮＴ…ネットワーク、ＤＳ０…気象庁データサーバ、ＤＳ１，ＤＳ２…レーダサイトデータサーバ、１６１…オゾン全量取得部、１６２…変換式作成部、１６３…波長別紫外線強度算出部。
【図１】
2008089425.tif 000003
【図２】
2008089425.tif 000004
【図３】
2008089425.tif 000005
【図４】
2008089425.tif 000006
【図５】
2008089425.tif 000007
【図６】
2008089425.tif 000008