三宅島2000年噴火時にGPS観測網が捉えた地殻変動を説明するため、名古屋大学地震火山防災研究センターの方々と、運動の原因となったマグマの移動や断層運動の特定を試みました。私の仕事は、名大の先生方の想像を計算機の中に具現化し、観測された現象を最もよく説明する断層あるいは膨張・収縮源の位置、運動パラメータを得ることでした。一番がんばったのは、一億回の試行を火を吹かずやってのけたワークステーション君でした・・・。


詳しくは、
YAMAOKA, K., M. KAWAMURA, F. KIMATA N. FUJII, and T. KUDO, Dike intrusion associated with the 2000 eruption of Miyakejima volcano, Japan, Bulletin of Volcanology, Vol. 67, No. 3, 231-242, 2005
をご参照ください。



以下解析途中の試行錯誤の軌跡・・・(もったいないので思い出に一部残しときます・・・?)


Clip Board for the MEIDAI MODEL...


Fixed Parameters:
Vp/Vs                    ====>   1.73
###### MAIN Fault(Finite rect.)######
Longitude                ====>  139.31
Latitude                 ====>   34.17
Strike-Angle(from E,c-cw)====>  -40.0
Dip-Angle       (degree) ====>   90.0
Fault Length        (km) ====>   20.0
Strike  -Dislocation (m) ====>    0.0
Dip     -Dislocation (m) ====>    0.0
###### SUB Fault (Point source)#######
Longitude                ====>  139.23
Latitude                 ====>   34.26
Source Depth        (km) ====>   10.0

Model 1: Dike only.


Best Fit Parameters:

### Main F.###
  DEPTH(a):[0.0km]// WIDTH(a):[7.0km]// Tens.DISLOC.(a):[4.0m]//
### Sub F.####
  Not considered.

Model 2: Dike and Aseismic point source



(a) Without constraints for the depth extent of the dike.

Best Fit Parameters:

### Main F.###
  DEPTH(a):[1.0km]// WIDTH(a):[20.0km]// Tens.DISLOC.(a):[2.0m]//
### Sub F.####
  Str.Angle(from E,c-cw)(b):[335.0]// DIP(b):[80.0]
  POT1:[-5.50]// POT2:[-1.50]// POT3:[2.50]// POT4:[-4.17]  (X10^8)


(b) Fixed the depth extent to be 15 km.

Best Fit Parameters:

### Main F.###
  DEPTH(a):[1.5km]// WIDTH(a):[15.0km](fixed)// Tens.DISLOC.(a):[2.0m]//
### Sub F.####
  Str.Angle(from E,c-cw)(b):[340.0]// DIP(b):[85.0]
  POT1:[-6.00]// POT2:[-1.00]// POT3:[2.50]// POT4:[-4.17]  (X10^8)


(c) Fixed the depth extent to be 10 km.

Best Fit Parameters:

### Main F.###
  DEPTH(a):[0.0km]// WIDTH(a):[10.0km](fixed)// Tens.DISLOC.(a):[1.5m]//
### Sub F.####
  Str.Angle(from E,c-cw)(b):[340.0]// DIP(b):[75.0]
  POT1:[-6.00]// POT2:[-2.50]// POT3:[3.00]// POT4:[-5.00]  (X10^8)




Bootstrap method