SNRの利用法

前にも書いたと思うが、Colorado大学のLarsonが、SNRを見ているとアンテナへの着雪が見える、という論文を出しているので、ちょっと自分でもやってみた。
Larson, K.M., A Methodology to Eliminate Snow and Ice-Contaminated Solutions from GPS Coordinate Time Series, J. Geophys. Res, Vol. 118(8), 4503-4510, doi:10.1002/jgrb.50307.

やったこと。
１）SNRのデータをまず仰角のsinの多項式で展開して、アンテナゲインパターンを調べる。
２）その関数を使って、指定する仰角(55度）のSNRを計算して、指標とする。

時系列を見てみると、Larsonの言う通り、なにもなければきれいな温度相関する変動がみられる。レドームの着雪が問題となっている電子基準点でみると、予想通り冬場に乱れて
それが座標の乱れとよく対応している。これは使えそう。

黒潮による荷重変形？

せっかくRTKLIBのPPPで全電子基準点の座標時系列を求めたので、鉛直座標のばらつきを調べてみた。まだちゃんと見てないけど、やっぱり北海道、東北脊梁部や日本海側でばらつきが大きい。多分積雪による荷重変形が見えているのだな。もちろんつくばの年周変動もきれいに見えている。あと、面白いのが、四国から九州にかけての太平洋岸に、帯状にばらつきが大きい点が分布していること。これ、なんだろ。黒潮による荷重変形かなぁ。だとするとちょっと面白いけど。ECCOでも見てみるかなぁ。

negative pseudorange

RTKLIBで昔のデータを解析していると、途中で解析ができなくなる事例があって、なんでかなぁ、と思っていたのだった。で、ちょっとトレースを見てみると、どうも擬似距離がおかしいらしい。見ると確かに負の値が入っている。これをはじいているわけだな。でもこれ、GIPSYでは解けるんだよね。なんでだろ？と思って調べてみたら、昔のデータでは、擬似距離が負になることもあるらしい。
http://postal.unavco.org/pipermail/teqc/2011/001150.html
本当は不要なんだけど、無理やりclockprepをかけて、time tagを直して解くことにしたら無事解けたのでよしとする。

デジタルDEM

GMTで三宅島の地図を描いてみたのだが、海岸線があまりに荒くてかっこよくない。なんかいいのないかなぁ、と探してみた。

簡単なのはSRTM3。90mメッシュ。
http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/

こと三宅島に関しては、ここから10mメッシュDEMが手に入るらしい。
http://fgd.gsi.go.jp/download/DEMkind.htm

ま、SRTM3でも当座十分なのだが、10mメッシュも気になるので、ちょっと時間ができたら使ってみよう。

あれ？よく調べたら、実は10mメッシュ、全国で手に入るのか。精度は落ちるし、50mメッシュから発生させている地域もあるようだが。これが昔の50mメッシュの代わりらしい。灯台もと暮らし。アホだorz例えばこれをこちらで公開していただいているソフトウェアで変換すればGMTで読めるみたい。
https://staff.aist.go.jp/a.noda/programs/jpgis2grd/jpgis2grd.html

時間があったら比較してみたい。

遅ればせながらの三宅島

ここのところ三宅島のデータをいじっていることもあって、集中的に過去の論文を読んでいるが、これがとても面白い。もう14年近く立っているせいもあって、地球物理、地球化学、岩石学などいろいろな分野で決定的な論文が出そろっており、ダイクの西方貫入からカルデラ形成、脱ガスまで、なにが起こったのか手に取るようにわかる。特に岩石学の論文は秀逸。同期のG氏や元同僚N氏の論文群も大変参考になった。うーん、こんなことやっていたのか（大恥。今回自分が得た結果は、残念ながらそういったビッグピクチャーをひっくり返すものではない（pico STAP程度？）のだが、その物理的な背景を考える材料となると思っているので、個人的には気に入っている。

しかし、これ、リアルタイムで貢献できていたらもっと楽しかっただろうなぁ。当時、ポスドクとして一応ちょっとは関わっていたはずなのだが、一体何やっていたんだろう。はぁ。