Graduate School of Science
Kyoto University

Institute for Geothermal Sciences

Annual Report FY2021

Annual Report FY 2021

令和三年度活動報告
Institute for Geothermal Sciences
Graduate School of Science
Kyoto University

京都⼤学
⼤学院理学研究科

附属地球熱学研究施設

序
地球熱学研究施設は、地球物理学研究施設（別府、⼤正１３年設⽴）と⽕⼭研究施設（阿蘇、昭和３年設⽴）とが、平成９年に統合改組された理学研究科附属施設である。地
球上で最⼤規模の⽕⼭・地熱温泉活動域の,と"である中部九州地域を巨⼤な実験装置とみなして、野外観測や室内実験な'を中⼼に、造構運動・⽕⼭活動・地熱温泉活動な'地
球の熱的活動に関する地球熱学の学問体系の構築をめしている。
この基本理念に⽴脚して、専⾨分野の異なる研究者が弾⼒的に協⼒でるように⼤部⾨制を採り、地熱流体論研究分野、地熱テクトニクス研究分野、⽕⼭構造論研究分野、⽕⼭活動論研究分野、地球熱学情報研究分野（外国⼈客員）の５"の研究分野が置かれている。
平成１年の国⽴⼤学法⼈化以降令和３年までの３期１８年にわたる中期計画の間、１２名の教員定員が８名に減少し、運費交付⾦も削減されるな'、施設運にとって厳しい状態が続いている。幸いにして、財政⾯の懸案に対しては、科学研究費な'の競争的資⾦の獲得が顕著であり、その内容は本年報の研究費(Funding)の項にまとめられている。
法⼈化とともに設置された施設運協議会が令和３年も５回開催された。コロナ禍２年⽬となる今年もす-てZoomを⽤いた=RM<R形式での開催となり、この開催形式
は今後定着するともわれる。学内での施設の研究教育⾯での位置づけを明確にすることを⽬的として設置された京都分室に関して、京都勤務者の任期に関する申し合わせが改定されたことにより、分室をこれまで以上に弾⼒的に運⽤することが可能となった。今後も京都分室を中⼼として、学部⽣の教育や研究科内での役割を担っていく所存である。
⻑引くコロナ禍の影響により、これまでの対⾯授業やTV会議を⽤いたセミナーに加え、Zoomを⽤いた=RM<R授業や=RM<Rと対⾯のI<KNCド形式の授業な'、講義・セミナー形態の多様化がさらに進んだ年であった。⼀⽅、遠隔地を活⽤した学部⽣向けの現地 在型実習は、最近２年間は実施出来ていない。緩やかな専⾨化を理念とする
理学部教育にいて、フィールドワークの魅⼒を伝える-実習が実施でないことは懸念す-事態であり、今後、地球惑星科学専攻との協⼒のもと、現地実習再開のための⽅策を検討していたい。
⼈事⾯では、令和３年１２⽉１⽇に澤⼭和貴⽒が助教として施設・本部（別府）に着任された。⼀⽅、⽳井千⾥研究員が令和３年１２⽉１５⽇に⾼知⼤学の特任助教として、梁熙俊研究員（研究 関）が令和４年２⽉１⽇に富⼭⼤学の特命准教授としてそれれ転出した。今後の両⽒の更なる活躍を祈念する。また、地球熱学情報研究分野の外国⼈教員招聘が取りやめになるな'、コロナ禍の影響は施設の⼈事計画にもよんでり、⼀刻も早いコロナ禍の収束を願うのみである。
最後に、２１年熊本地震の本震により被災した⽕⼭研究セRBーの復旧記念式典
が、令和３年５⽉に=RM<Rと対⾯のI<KNCド形式により実施されたことを報告し、これまでの各所からの温かい⽀援に、あらためて感謝申し上げる。
令和４年８⽉
令和３年地球熱学研究施設⻑
⼤倉敬宏

Preface
Institute for Geothermal Sciences (IGS) is an affiliated institute of the
Graduate School of Sciences, Kyoto University, integrated and reorganized fromthe Geophysical Research Laboratory (Beppu, established in 1924) and theVolcanological Research Laboratory (Aso, established in 1928) in 1997.
Central Kyushu, one of the largest volcanic and geothermal areas in theworld, is regarded as a gigantic laboratory, where fieldwork and laboratory
experiments are conducted to establish academic frameworks of geothermalsciences. Based on this fundamental scope of our research, we have the following
five research units in our institute, geothermal fluids, geothermal tectonics, volcanic
structure, volcano-dynamics, and geothermal intelligence section (visiting researchscholars from abroad), so that a variety of research works can flexibly cooperate
within this interdisciplinary geothermal science research system.
During the three periods of the 6-year mid-term plan beginning the 2004
fiscal year, the number of faculty members has been reduced from 12 to 8, andManagement Expenses Grants have also been reduced, resulting in a continuing
difficult situation for the management of the institution. Fortunately, the financial
concerns have been addressed through the acquisition of competitive funds suchas the Grants-in-Aid for Scientific Research, which are summarized in the section on
Research Funding in this annual report.
Meetings of the steering committee of IGS, which was established in 2004,
were held five times in the 2021 fiscal year. The meeting has been held online using
Zoom since the COVID-19 pandemic and we plan to continue this online style in
the future. The Kyoto branch of IGS, which was established in 2007, has been able
to operate flexibly due to the revision of the agreement on the tenure of Kyoto
staff. The Kyoto branch will continue to play an important role in the education ofundergraduates and negotiations with the Graduate School of Sciences.
Due to the prolonged COVID-19 pandemic, various lecture/seminar styles
were established, such as online classes using Zoom and a hybrid of online andface-to-face classes, in addition to conventional face-to-face classes and seminars
using video conferencing. On the other hand, on-site field training courses withaccommodation in Beppu and Aso for undergraduates have not been available
for the last two years. This is a worrisome situation in undergraduate education,
which is based on the principle of gradual specialization in the Faculty of Science.
Therefore, we and faculty members of the Geophysics department are working
together to offer the training course with accommodation safely in the coming
fiscal year.
In terms of personnel, Dr. Kazuki SAWAYAMA was appointed as an assistantprofessor at the Institution Headquarters (Beppu) on December 1, 2021. Meanwhile,
Dr. Chisato ANAI moved to Kochi University as a specially-appointed assistantprofessor on November 15, 2021, and Dr. Heejun YANG moved to Toyama University

as a specially-appointed associate professor on February 1, 2022. The impact of the
COVID-19 pandemic has also affected the personnel plan of the institution, as theinvitation of foreign faculty members in the field of geothermal information
research has been canceled.
Finally, we would like to report that a ceremony to commemorate therestoration of Aso Volcanological Laboratory, which was severely damaged by themain shock of the 2016 Kumamoto earthquake, was held in May 2021 in a hybrid of
online and face-to-face style and would like to thank all of you again for yourheartfelt support.
August 2021
Takahiro OHKURA,
Director of the 2021 fiscal year

⽬次 Contents
序Preface
1． 構成員Member2． 研究
Research
2.1． 研究報告
Scientific Report
2.2． 論⽂
Papers2.3． 講演・⼝頭発表
Presentations
2.4． 書籍等出版物
Books & Other Publications
2.5． 外部資⾦：科学研究費補助⾦
External Funds: KAKENHI
2.6． 外部資⾦：その他
External Funds: Others
3． 教育
Education
3.1． 担当科⽬
Teaching Subjects
3.2． 学位審査
Participation in PhD/Master’s Defense
4． ⼤学運営
Administration
4.1． 全学運営（役職等）
School management
4.2． 部局運営（役職等）
Faculty management
5． 学術・社会貢献
Academic, Social Contribution
5.1． 委員歴
Committee Members
5.2. 学術貢献活動
Academic Contribution
5.3． 社会貢献活動
Social Contribution
6． セミナー
Seminars
7． ⼀般公開報告
Openhouse
8． 来訪者
Visitors
9． 定常観測
Routine Observations
10． 装備・施設
Facilities and Instruments

1
2
2
9
12
13
13
15
16
16
19
21
21
21
22
22
23
23
24
26
30
32
33

．構成員 Members

1

教授 ＊
⼤倉敬宏
⼤信⼆
楠本成寿
＊
施設⻑

⼩澤佑⼈ Yuto Kozawa
菅原崇史 Takashi Sugawara
樋⼝和 Kazuya Higuchi
⽯橋 桜 Sakura Ishibashi
⼤久保美鈴 Misuzu Okubo
梨元 昴 Subaru Nashimoto
服部匠哉 Takuya Hattori
政本⾵⼈ Futo Masamoto
事務補佐員 Assistant Administrative Staff
⼀万⽥ 歩 Ayumi Ichimanda
宮⽥美保 Miho Miyata
⼟井有紀 Yuki Doi
研究⽀援推進員 Research Support Staff
⼩⼭晴⼦ Seiko Koyama
臨時⽤務員 Supply Janitor
⼭代 Sakiyo Yamasaki

Professors
Takahiro Ohkura
Shinji Ohsawa＊
Shigekazu Kusumoto＊Director

准教授
Associate Professors
古川善紹 Yoshitsugu Furukawa
柴⽥智郎 Tomoo Shibata
横尾亮彦 Akihiko Yokoo
助教
Assistant Professors
宇津⽊ 充 Mitsuru Utsugi
澤⼭和貴 Kazuki Sawayama
2021年12⽉1⽇着任
技術専⾨員 Technical Professionals（理学研究科技術部）
⾺渡秀夫 Hideo Mawatari
吉川 慎
Shin Yoshikawa
技術専⾨職員 Specialist Technical Staff（理学研究科技術部）
三島壮智 Taketoshi Mishima
技術職員 Technical Staff（理学研究科技術部）
井上寛之 Hiroyuki Inoue
研究員（研究） Researchers(Research Organization)
梁 熙俊
Heejun Yang
2022年1⽉31⽇離任
⽳井千⾥ Chisato Anai
2021年11⽉15⽇離任
成⽥翔平 Shohei Narita
齋藤 圭
Kei Saitoh
⼤学院⽣ Graduate Student
⽯井佳 Kyoka Ishii
⾕ 協⾄
Kyoji Tani
1．研究 Research

2

2.1.

研究報告 Scientific Report

Temporal change with rock magnetic properties of volcanic ashes ejected during a
one-year eruption event: A case study on the Aso Nakadake 2019–2020 eruption
C. Anai (Kochi Univ.), N. Mochizuki (Kumamoto Univ.),
S. Yoshikawa, T. Ohkura
We investigated temporal changes in the magnetic properties of volcanic ash ejectedfrom the Aso Nakadake volcano during a sequence of ash eruptions from 2019 to 2020.
Titanium-rich titanomagnetite and titanium-poor titanomagnetite were the dominantmagnetic minerals in the samples where titanium-rich titanomagnetite was more
dominant. From the rock magnetic measurements, parameters such as the saturation
remanent magnetization (Mrs), saturation magnetization (Ms), coercivity (Hc), andtitanium content estimated from the Curie temperature (Tc) were extracted andchecked for their temporal changes. The magnetic behavior of the magnetic minerals
was confirmed by the increasing values of Mrs/Ms and Hc at several periods. The samples
with higher values of Mrs/Ms and Hc included titanomagnetite with a low Tc (hightitanium content). The clear increase in Mrs/Ms and Hc suggests that the ratio of thesingle-domain volume fraction increased, indicating that the titanomagnetite particles
became finer in size. Interestingly, the periods of high Mrs/Ms and Hc were synchronizedwith observations of the volcanic glow. These results suggest that changes in themagnetic properties of volcanic ash reflect changes in physical and/or thermal
conditions from the vent to the conduit. (This report has been submitted to EPS.)

Fig. 1 Hysteresis parameters (Mrs/Ms and Hc) in chronological order. Gray squares indicate Hc and blue triangles
denote Mrs/Ms. Red squares indicate the day when the volcanic glow was observed. High Mrs/Ms and Hc
mean the predominance of single domain (SD) particles.

2

Spatio-temporal change in ground deformation of Kuju volcano, as revealed by L-band
InSAR time series (1992-2021)
S. Narita, T. Ohkura
Kuju volcano is one of the active volcanoes of Beppu-Shimabara graben zone.
Iwoyama and Mt. Hossho, those which are at the central part of Kuju volcano, haveexperienced numerous phreatic eruptions in recent 5,000 years (Itoh+ 2014). The latest
eruption occurred in Oct. 1995, forming eruptive vents to 500 m south of Iwoyama.
Since then, significant ground deflation caused by active fumarolic plume discharge
had been observed at Iwoyama and northeast of Mt. Hossho until, at least, 2003
(Nakaboh+ 2003). From 2003, there have been no detailed investigation about the
deformation. To investigate spatio-temporal change in the deformation field aroundIwoyama, we here analyzed L-band Synthetic Aperture Radar (SAR) data of JERS-1(1992-1998), ALOS-1/PALSAR-1 (2006-2011), and ALOS-2/PALSAR-2 (2014-2021).
For SAR interferometry, we used Radar Interferometry Calculation Toolbox (Ozawa+
2016). Then, we applied NSBAS (New Small Baseline subset) mode of GIAnT (Generic
InSAR Analysis Toolbox) to all the interferograms to estimate displacement time series. To
disclose slower deformation which cannot be clearly detected by NSBAS, we applied
StaMPS PS-InSAR to ALOS-2 data after stacking ALOS-2 SLCs (single-look complex
images) by using ISCE2. Finally, we fitted linear function to the time series to extract only
volcanic deformation velocity.
For overall period (1992-2021), line-of-site (LOS) distance increases have been
observed around Iwoyama, indicating long-term ground deflation. However, the spatial
scale clearly differs; 1.5-2 km in JERS-1, 1 km in ALOS and 0.8 km in ALOS-2 results. Thissuggests that the the source depth has been gradually becoming shallower. Although
the spatial extent in the ALOS-1 and ALOS-2 results are similar, quasi east-west (qEW)
and up-down (qUD) components of the velocity fields indicate perfectly different
characteristics. While ALOS-1 results show a simple deflation pattern which can be
explained by simple geometry source, such as a sphere, ALOS-2 results show a more
complex pattern. qEW deformation field shows consistent eastward displacement (2.5cm/yr) over region in which LOS increases are observed (Iwoyama lava dome andnorth-east flank of Mt. Hossho). qUD deformation field shows consistent subsidence (3cm/yr) over the same region. These patterns are not due to simple inflation/deflation of
single pressure source. Focusing close-up of these deformation, we can see two
deformation features; one is landslide-like slope instability, which is seen at the west halfof the deformed area (northeast flank of Mt. Hossho), the other one is local scale (~ 150m) subsidence, which is seen at the east half of the deformed area (the lava dome of
Iwoyama). The subsidence area at Iwoyama overlaps the most active fumarolic area,suggesting that the local subsidence is caused by pore pressure decrease within
near-surface fumarolic conduits.
Other important feature of ALOS-2 results is a uplift-like LOS change over wider areaaround Iwoyama (~2 km diameter). In the results of NSBAS, its signal-to-noise ratio is sosmall that we cannot confirm the uplift-like phase pattern as real deformation signal. To
confirm this, we performed PS-InSAR analysis and succeeded in extracting a very slow
uplift-like deformation (0.4 cm/yr). Spatial pattern and extent of the uplift are similar to
those of the deflation observed during post-eruptive period (1995~1998), indicating thata common hydrothermal reservoir, which is modeled as a deflating spherical source3550-730 m below fumarole area of Iwoyama, had depressurized just after the 1995
eruption due to violent fluid discharge (Nakaboh+ 2003) until at least 2003, and thenstarted to be pressurized in 2014. This is consistent with GNSS baseline length changes,
which show onset of ground inflation around Iwoyama in 2012-2013.
All the observed phenomena from 2012 (pressurization at ~600 m depth anddepressurization at near-surface, and very week thermal activity at surface of
Iwoyama) probably suggest ongoing partial sealing of the fumarolic conduits justbeneath Iwoyama. This is likely due to precipitation of hydrothermal minerals or porecollapse / closing due to long-term depressurization (Hamling+ 2016). To validate thisconceptual model, it is necessary to evaluate mass and heat balance between thepressure sources and surface heat flux.
Prototype geochemical towing observation system for detection of seafloor dischargeof volcanic fluids in coastal waters and its evaluation of performance at volcanic CO2
upwelling area
S. OHSAWA, K. AMITA (Akita University)
T. MISHIMA, K. SAITOH, F. MASAMOTO,
H. TAKAHASHI(AIST), N. MORIKAWA(AIST)
A prototype "Geochemical Towing Observation System" was developed to detect
the presence of volcanic fluids discharging from the coastal seafloor in shallow watersfrom the sea (Fig.1). We confirmed the stable operation of the observation system andevaluated its ability to detect volcanic fluid discharge in one of the coastal waters ofHimeshima Island (Oita Prefecture, Japan), where volcanic CO2 gas has beendischarged from the seafloor. In the observed sea area, as the concentration ofdissolved inorganic carbon (DIC) in seawater increases, its
the estimatedδ C value changes toward13δ C of volcanic CO so we believe that the observation system has
13

2,

successfully captured the seafloor discharge of volcanic fluids. In addition to DIC, asignificant change in redox potential (ORP /Eh) also appeared (Fig.2), which can be
explained by the detection of H2S dissolved in seawater contained in seafloor gases.
Combined measurements of 222Rn, salinity corresponding to electric conductivity, andwater temperature, which have been effective in detecting submarine groundwaterand hot spring water discharges from the seafloor in the Beppu Hot Springs coastal
area in the back of Beppu Bay, showed no expected changes. This may indicate thatsubmarine groundwater discharge from the seafloor along the coast of HimeshimaIsland is not significant and that gas discharge is the main cause.
(Submitted to Journal of Japanese Association of Hydrological Sciences)4

Fig. 1 (a) A view of towing survey by a small boat (a test towing off the east coast of Himeshima Island on
Oct.29, 2019). (b) Layout of prototype geochemical towing observation system composed of two
independent measuring equipment units on the boat in this study

Fig.2 Sequential changes in (a)water temperature, EC (electric
conductivity),
pH, Eh andconcentration of dissolved H2 by
the measuring equipment unit 2
withflow-celland(b)concentrations of dissolved
222
Rn and DIC by the meas equip
unit 1 with spray chamber thatbrings the air and water into
equilibrium
duringthecontinuous towing observationof surface seawater off thecoast Nishiumra. (C) Sequentialchange in DIC concentration inseawatersamplecollected
every 2 minutes from the outletof the flow-line in the meas.
equip unit 2. The highlighted
times P and Q correspondroughly to the times when theobservation boat approachedand left the gas upwelling area,
respectively.

Laboratory experiments of pressure-interference and tracer tests of simulated
geothermal reservoirs involving transport and dead-end pores
K. Sawayama and M. Matsumoto (Kyushu Univ.)
The permeability and porosity mainly control the thermal transport in geothermalareas. In the field experiments, we use pressure-interference and tracer tests, which
individually determine the in-situ value of the permeability and porosity. However, somefield data show a large discrepancy between the estimated value of each method.
One of the possible causes of this discrepancy is the heterogeneous porosity distribution
in nature, such that the pressure-interference (diffusion) and tracer transport

5

(advection) depend on different void areas. Based on the rock physical point of view,a “transport pore” affects both diffusion and advection, whereas a “dead-end pore”
only affects diffusion, not advection. In this study, we introduce these two kinds of pores
to simulate heterogeneous pore structures of actual geothermal reservoirs. In order todevelop a numerical model of comprehensively characterizing pressure changes andadvection in a reservoir, we demonstrate the different void areas ofpressure-interference and tracer tests in the laboratory experiment.
The transport and the dead-end pores were simulated by a tube (an innerdiameter of 4 mm and a total channel length of 2 m) and a tank (a diameter of 100
mm and a height of 100 mm), respectively (Figure 1). We performed a laboratory
fluid-flow test by an apparatus involving these voids to reproduce the pressure andtracer responses. The volume change of the dead-end pore was also simulated bychanging the bulk compressibility (changing the gas-liquid ratio in the tank). The
transient pressure change was driven by the opening of a valve on the higher pressure(higher water level) side, and was measured by three pressure gauges connected to adata logger (GL7000, Graphtec Co., Ltd.). In this experiment, we recorded voltage
values at a sampling rate of 10 kHz, and resampled them to 1 kHz via the Chebyshev
filter for noise reduction. After confirming that the pressure reached a steady state, wethen conducted a tracer test. The tracer was an aqueous KCl solution (10 wt.%), andelectrodes were placed downstream to measure the electrical conductivity between
the electrodes with a multimeter (VOAC7602, Iwasaki Tsushinki Co., Ltd.). The time
variation of the tracer concentration was estimated from the obtained electricalconductivity. This measurement was performed at a sampling rate of 100 Hz, and theraw data was interpolated by a spline method to reduce noise.

Figure1. Outline of the experimental apparatus

Figure 2(a) shows the results of pressure change at Pdown. The horizontal axis shows theelapsed time after the valve opening, and the vertical axis shows the pressure value
converted from the voltage value of the pressure gauge. The pressure change exhibits
different behavior depending on the tank water level (i.e., bulk compressibility). Figure 2(b) shows the results of tracer concentration changes. The horizontal axis shows theelapsed time after tracer injection, and the vertical axis shows the relative tracer
concentration normalized by the injected tracer concentration. In contrast to the6pressure change, tracer concentration shows no significant difference depending on
the tank water level (i.e., bulk compressibility). These results suggest that the pressurechange is affected by the compressibility (=volume) of the dead-end pore anddepends on the total porosity of the transport and dead-end pores. This experimental
result that the volume of dead-end pore contributes to the pressure change but not thetracer concentration change suggests the possibility of explaining the different pressure
and tracer responses observed in geothermal reservoirs.

Figure 2. Plots showing transient changes of (a) pressure change and (b) tracer concentration. Sw
represents the liquid-to-gas ratio in the tank.

Numerical calculation of groundwater flow for local distribution and short-time variation
N. Takita, T. Shibata
Numerical methods are often used to determine hydraulic head distribution of thetarget aquifer. The spatial distributions and temporal variations of the hydraulic headare generally calculated based on the diffusion equation in groundwater flow(1), which
is obtained by the Darcy's law (Pinder and Bredehoeft, 1986),
∂ℎ

∂

(

∂ℎ

)

∂

(

∂ℎ

)

𝑆𝑆𝑆 ∂𝑡𝑡𝑡 = ∂𝑥𝑥𝑥 𝑇𝑇𝑇𝑥𝑥𝑥 ∂𝑥𝑥𝑥 + ∂𝑦𝑦𝑦 𝑇𝑇𝑇𝑦𝑦𝑦 ∂𝑦𝑦𝑦 − 𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊

2where h hydraulic head [L], S storage coefficient, T transmissibility [𝐿𝐿𝐿 /𝑇𝑇𝑇𝑇, W function
describing the net water amount into or out of the system per unit area and unit
time[L/T], respectively. Numerically solving the equation (1) requires the development ofa finite-difference approximation. Actually, we can assume an isotropic andhomogeneous aquifer, the derivative functions are replaced by the second-order
central difference functions at separate discrete points of the equation (2),
∂ℎ
=
∂𝑡𝑡𝑡

2ℎ𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖+ℎ𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖ℎ−2ℎ + ℎ𝑊𝑊𝑊
𝑇𝑇𝑇 ⎡ ℎ𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖−
𝑖
𝑖+ 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖2𝑖 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 ⎤⎥ − 𝑆𝑆𝑆𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 #(2)2
𝑆𝑆𝑆 ⎢𝑥∆𝑥𝑥𝑥𝑥∆𝑦𝑦𝑦𝑦𝑦⎣⎦where i, j, and k are the index of x, y, and t direction respectively, ∆𝑥𝑥𝑥 and ∆𝑦𝑦𝑦 are theincrements in the x and y dimensions, respectively. The general-purpose software,
MODFLOW (Hoffmann et al., 2003), solves the diffusion equation (1) using asecond-order backward finite difference approximation for the time derivative(Hoffmann et al, 2003; Harbaugh, 2005). The method is suitable for obtaining spatial
distributions of hydraulic head in a large scale, but not for local distributions or

7

short-time differentiation. Therefore, to improve time accuracy, we developed acalculation program which solves the equation (2) using forth-order Runge-Kutta
method.
In order to verify the accuracy of the program, we tried to perform two methods. First,
we solved the problems of constant pumping in a well located in the center of aconfined, homogeneous, and isotropic aquifer. The problem can be analytically solved
by Theis (1935),
𝑄𝑄𝑄 ∞ 𝑒𝑒𝑒−𝑧𝑧𝑧

and𝑑𝑑𝑑𝑑𝑑𝑑𝑑ℎ = ℎ0 − 4π𝑤𝑤𝑤𝑇𝑇𝑇 ∫ 𝑧𝑧𝑧 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑
𝑢𝑢𝑢

𝑢𝑢𝑢𝑢𝑢𝑢𝑢 (4𝑇𝑇𝑇𝑇𝑇)2

𝑟𝑟𝑟 𝑆𝑆𝑆

where ℎ0 the initial hydraulic head[L] and 𝑄𝑄𝑄𝑤𝑤𝑤 the pumping rate [L3/T]. In contrast, we

can calculate the problem using the equation (2) with the calculation area of
500[m]×500[m] with hydraulic heads at the edges were constant, and the parametersof S, T, ℎ0, 𝑄𝑄𝑄𝑤𝑤𝑤 ∆𝑡𝑡𝑡 and ∆𝑥𝑥𝑥 of 3.0×10-3, 9.0×10-3[m/s], 0[m], 1.76×102[m3/s], 4.39×10-2[s] and

5[m], respectively (Pinder and Bredehoeft, 1986). The results obtained by the analyticalsolution and the numerical calculation are shown in Fig.1, which are almost the same
except for the time over 103[s]. Although the difference in long time could be due tothe boundary conditions, we are able to confirm that our program is appropriate.

The other method is to confirm whether the pumping test data (Meidai Kougyou Ltd.,
1988) can be reproduced by numerical calculation using our program. The pumping
test was conducted on December 14th, 1987 at the well with a depth of 300m, which islocated in the Institute for Geothermal Sciences, Kyoto University. The decreasing ofwater head by the pumping is fitted by a regression curve of Jacob’s method (Jacob,1950), and the field parameters of storage coefficient, S, and transmissibility, T are
estimated as 2.87×10-3, and 3.66×10-2[m2/s], respectively. Using the parameters, theamount of decrease in water head by pumping can be numerically calculated at any
time. Some decreases in water head estimated from numerical calculation are
corresponded to the observed pumping data (Fig.2).
We created a new program which solves the diffusion equation in groundwater flowusing fourth-order Runge-Kutta method, and suggested that it can calculate for thefields under local distribution and short-time variation.

8

References:
Hoffmann, J., Leake, S.A., Galloway, D.L., and Wilson, A.M. (2003) MODFLOW–2000
ground-water model–User guide to the subsidence and aquifer-systemcompaction (SUB) package. U.S. Geological Survey Open-File Report 03–233.
Harbaugh, A. W. (2005) MODFLOW–2005, the US Geological Survey modular
ground-water model: the ground-water flow process. U.S. Geological Survey
Techniques and Methods 6–A16, 9 chapters.
Jacob, C.E. (1950) Flow of groundwater. In Engineering Hydraulics (ed. Rouse, H.) John
Wiley & Sons, New York, 321–386.
Meidai Kougyou Ltd. (1988) Report of experiment well at Institute for Geothermal
Sciences, Faculty of Science, Kyoto University.
Pinder, G.F. and Bredehoeft, J.D. (1968) Application of the digital computer for aquifer
evaluation. Water Resources Research, 4, 1069–1093.
Theis, C.V. (1935). The relation between the lowering of the piezometric surface and therate and duration of discharge of a well using ground-water storage. Eos, 16,519-524.

2.2.

論⽂ Papers

⼤倉 敬宏

阿蘇⽕⼭中岳第1⽕⼝における2019〜2020年マグマ噴⽕の先駆活動－2019年5⽉3⽇〜5⽇噴⽕を含む⽕⼭活動の特徴－, 宮縁育夫; 飯塚義之; 遠⼊楓⼤; ⼤倉敬宏, ⽕⼭, 2021年09⽉, 査読有り
The 2018 phreatic eruption at Mt. Motoshirane of Kusatsu–Shirane volcano, Japan:

eruption and intrusion of hydrothermal fluid observed by a borehole tiltmeter network,
Akihiko Terada; Wataru Kanda; Yasuo Ogawa; Taishi Yamada; Mare Yamamoto;
Takahiro Ohkura; Hiroshi Aoyama; Tomoki Tsutsui; Shin’ya Onizawa, Earth, Planets and

年12⽉Space, 2021

Locating hydrothermal fluid injection of the 2018 phreatic eruption at Kusatsu-Shirane
volcano with volcanic tremor amplitude, Taishi Yamada; Aika K. Kurokawa; Akihiko
Terada; Wataru Kanda; Hideki Ueda; Hiroshi Aoyama; Takahiro Ohkura; Yasuo Ogawa;

年 ⽉
古地磁気・岩⽯磁気学的⼿法を⽤いた阿蘇⽕⼭中岳⽕⼝周辺域の完新世噴出物の分類と噴⽕年

Toshikazu Tanada, Earth, Planets and Space, 2021 12

9

代の再検討, ⽳井千⾥; 宮縁育夫; 宇津⽊充; 吉川慎; 望⽉伸⻯; 渋⾕秀敏; ⼤倉敬宏, ⽕⼭, 2021年

⼤沢 信⼆

九重硫⻩⼭における⾃然電位観測, ⼤分県温泉調査研究会報告, 網⽥和宏; ⼤信⼆, ⼤分県温泉調査研究会報告, 2021年
分析機器を⽤いない噴気の簡易ガス分析, ⼤信⼆; 三島壮智; 齋藤武⼠; 網⽥和宏, ⼤分県温泉調査研究会報告, 2021年

楠本 成寿

Deformation source revealed from leveling survey in Jigokudani valley, Tateyama
volcano, Japan, Kohei Hotta; Shigekazu Kusumoto; Hidenori Takahashi; Yuichi S.

年1⽉, 査読有り
ブーゲー異常を⽤いた浅部間隙率分布の推定, 楠本成寿; ⾼橋秀徳; 東中基倫; 早川裕弌, 物理探査, 2021年, 査読有り

Hayakawa, Earth, Planets and Space, 2022

柴⽥ 智郎

別府扇状地南部域における不圧地下⽔位の⻑期変動, 梁熙俊; 柴⽥智郎, 地下⽔学会, 2021年8⽉, 査読有り横尾 亮彦

Sampling Volcanic Plume Using a Drone-Borne SelPS for Remotely Determined Stable
Isotopic Compositions of Fumarolic Carbon Dioxide, Urumu Tsunogai; Ryo Shingubara;
Yuhei Morishita; Masanori Ito; Fumiko Nakagawa; Shin Yoshikawa; Mitsuru Utsugi; Akihiko

年 ⽉, 査読有り
熱⾚外線カメラで撮影されたJPEGファイルの解析⽅法, 横尾亮彦; ⽯井佳, ⽕⼭, 2021年9⽉,
査読有り, 筆頭著者, 責任著者

Yokoo, Frontiers in Earth Science, 2022 3

Combined approach to estimate the depth of the magma surface in a shallow conduit

年

at Aso volcano, Japan, Kyoka Ishii; Akihiko Yokoo, Earth, Planets and Space, 2021 12

⽉, 査読有り

Development of a drone-borne volcanic plume sampler, Ryo Shingubara; Urumu
Tsunogai; Masanori Ito; Fumiko Nakagawa; Shin Yoshikawa; Mitsuru Utsugi; Akihiko

年 ⽉, 査読有り

Yokoo, Journal of Volcanology and Geothermal Research, 2021 4

宇津⽊ 充

Aeromagnetic survey in Kusatsu-Shirane volcano, central Japan, by using an
unmanned helicopter, Takao Koyama; Wataru Kanda; Mitsuru Utsugi; Takayuki Kaneko;
Takao Ohminato; Atsushi Watanabe; Hiroshi Tsuji; Taro Nishimoto; Alexey Kuvshinov;

年 ⽉, 査読有り

Yoshiaki Honda, EARTH PLANETS AND SPACE, 2021 7

Magmatic fluid pathways in the upper crust: Insights from dense magnetotelluric

10

observations around the Kuju Volcanoes, Japan, Koki Aizawa; Mitsuru Utsugi; Keigo
Kitamura; Takao Koyama; Makoto Uyeshima; Nobuo Matsushima; Shinichi Takakura;
Haruhiro Inagaki; Hiroki Saito; Yasuhiro Fujimitsu, Geophysical Journal International, 2022

年2⽉
古地磁気・岩⽯磁気学的⼿法を⽤いた阿蘇⽕⼭中岳⽕⼝周辺域の完新世噴出物の分類と噴⽕年代の再検討, ⽳井 千⾥; 宮縁 育夫; 宇津⽊ 充; 吉川 慎; 望⽉ 伸⻯; 渋⾕ 秀敏; ⼤倉 敬宏, ⽕⼭,
2021年Electrical conductive fluid-rich zones and their influence on the earthquake initiation,
growth, and arrest processes: observations from the 2016 Kumamoto earthquake
sequence, Kyushu Island, Japan, Koki Aizawa; Shinichi Takakura; Hisafumi Asaue;
Katsuaki Koike; Ryokei Yoshimura; Ken’ichi Yamazaki; Shintaro Komatsu; Mitsuru Utsugi;
Hiroyuki Inoue; Kaori Tsukamoto; Makoto Uyeshima; Takao Koyama; Wataru Kanda;
Tohru Yoshinaga; Nobuo Matsushima; Kazunari Uchida; Yuko Tsukashima; Takeshi
Matsushima; Hiroshi Ichihara; Dan Muramatsu; Yoshiko Teguri; Azusa Shito; Satoshi

年 ⽉

Matsumoto; Hiroshi Shimizu, Earth, Planets and Space, 2021 12

澤⼭ 和貴

Elastic wave velocity changes due to the fracture aperture and density, and direct
correlation with permeability: an energetic approach to mated rock fractures, K.
Sawayama; T. Ikeda; T. Tsuji; F. Jiang; O. Nishizawa; Y. Fujimitsu, Journal of Geophysical

年 ⽉, 査読有り, 筆頭著者, 責任著者
機械学習による弾性波速度測定実験の⾃動解析技術, 諏訪 雄⼆; 顧 志明; ⼩幡 拓也; 中川 慎⼆;
森下 諒⼀; 澤⼭ 和貴; 辻 健, 電⼦情報通信学論⽂誌 D, 2022年1⽉, 査読有り

Research: Solid Earth, 2022 1

Digital rock physics revealing the relationships between permeability, resistivity andelastic wave velocity of rock fractures, Kazuki Sawayama; Takuya Ishibashi; Fei Jiang;
Tatsunori Ikeda; Takeshi Tsuji; Osamu Nishizawa; Yasuhiro Fujimitsu, Proceedings of the

年11⽉, 査読有り,

14th SEGJ International Symposium, Online, 18–21 October 2021, 2021

筆頭著者, 責任著者

Scale-independent relationship between permeability and resistivity in mated fractureswith natural rough surfaces, K. Sawayama; T. Ishibashi; F. Jiang; T. Tsuji; O. Nishizawa; Y.

年7⽉, 査読有り, 筆頭著者, 責任著者
圧⼒・熱拡散の%>R3ーによる模擬坑井試験 : 熱率計を⽤いた⼀事例, 松本 光央; 澤⼭ 和
貴; 古賀 翼, ⽇本地熱学誌 = Journal of the Geothermal Research Society of Japan,
2021年4⽉, 査読有り

FujimitsuGeothermics, 2021

梁 熙俊

⽔⽂化学的⼿法と多変量解析による 別府温泉の泉質の分類および鉛直分布特性, 梁 熙俊, 地下
⽔学誌, 63, 137-149, 2021年8⽉, 査読有り, 筆頭著者, 責任著者
別府⽕⼭性"状地南部域におる不圧地下⽔位の⻑期変動, 梁 熙俊, 柴⽥智郎, 地下⽔学誌,
11

年8⽉, 査読有り, 筆頭著者, 責任著者

63, 151-157, 2021

Large-scale hydraulic conductivity distribution in an unconfined carbonate aquifer
using the ocean tidal propagation, Heejun Yang, Yasuhiro Tawara, Jun Shimada,

年 ⽉, 査読有り, 筆頭著者,

Makoto Kagabu, Azusa Okumura, Hydrogeology Journal, 2021 6

責任著者

⽳井 千⾥

古地磁気・岩⽯磁気学的⼿法を⽤いた阿蘇⽕⼭中岳⽕⼝周辺域の完新世噴出物の分類と噴⽕年代の再検討, ⽳井 千⾥, 宮縁 育夫, 宇津⽊ 充, 吉川 慎, 望⽉ 伸⻯, 渋⾕ 秀敏, ⼤倉 敬宏, ⽕⼭, 66,
171-186, 2021年, 査読有り, 筆頭著者

齋藤 圭

Mass balance of saline lakes considering inflow loads of rivers and groundwater: thecase of Lake Issyk-Kul, Central Asia, Kei Saitoh; Rysbek Satylkanov; Kenji Okubo, Journal

年 ⽉, 査読有り, 筆頭著者, 責任著者

of Arid Land, 2022 1

Understanding the water environment characteristics in various rivers, Yoichi Morimoto;
Yoshihiro Igari; Kei Saito; Erina Yamagata; Muneo Takemoto; Akihiko Naemura; Koji

年 ⽉, 査読有り

Kodera, International Journal of Human Culture Studies, 2022 1

2.3.

講演・⼝頭発表 Presentations

⼤沢 信⼆

別府湾の沿岸海底温泉湧出の探索, ⼤沢信⼆; 九州⾼等学校理科教育研究会⼤分⼤会, 2021年7⽉, 招待有り
⾼塩分温泉の成因に関する分野横断的研究, ⼤沢信⼆; ⼤分地質学会, 2021年4⽉, 招待有り
⼤分野にる温泉⽔中の%SK?%態窒素安定同位体⽐組成（AHW23-09）, 中⾼志; ⼤
沢信⼆; 網⽥和宏, ⽇本地球惑星科学連合2021年⼤会, 2021年6⽉

柴⽥ 智郎

地震波にり誘発される地下⽔中⼆酸化炭素の気泡成⻑シミュレーシMS, 服部匠哉; 柴⽥智郎,
⽇本地球化学会68回年会, 2021年9⽉
別府南部地域にる不圧地下⽔位の周期変動にい, 梁 熙俊; 柴⽥智郎, ⽇本地下⽔学会
2021年秋季講演会, 2021年11⽉

澤⼭ 和貴

Relationship between permeability and resistivity explored via the numerical simulationof sheared fractures Sawayama, K; Ishibashi T; Jiang; Tsuji, TAmerican Geophysical Union

年12⽉

Fall meeting, 202112

デジタル岩⽯物理とエネルギー論的アプローチによる岩⽯⻲裂の弾性波速度推定, 澤⼭和貴; ⻄澤修; 池⽥達紀; 蒋 ⾶; 辻 健; 藤光康宏, ⽇本地球惑星科学連合2021年⼤会, 2021年06⽉

齋藤 圭

豊後国⽥染荘資料のデータ化による研究基盤の整備と新たな研究に向けた試み, ⾚松秀亮; 齋藤
圭, オンライン国 シンポジウム「宇佐とローマをつなぐ 古代の街道の⽐較研究」, 2022年1⽉
2.4.

書籍等出版物 Books & Other Publications

⼤沢 信⼆

図説⽇本の温泉-170温泉の1イエンス-（⼀般財団法⼈⽇本温泉科学会監修）」分担執筆（Ⅵ九州地⽅：6-11⼭⾹温泉, 6-12雲仙温泉, 6-13阿蘇⽕⼭の温泉，6-15宮崎平野の⼤深度温泉）, 朝
倉書店

楠本 成寿

An introduction to probabilistic fault displacement hazard analysis in site evaluation forexisting nuclear installations,

年11⽉

Agency, 2021

2.5.

共編者(共編著者), 全般, International Atomic Energy

外部資⾦：科学研究費補助⾦ External Funds: KAKENHI

⼤倉 敬宏

年熊本地震は阿蘇カルデラ噴⽕をトリガーするか？, 基盤研究(C), ⼩区分17040:固体地球科学関連, 京都⼤学, ⼤倉敬宏, ⾃ 2018年4⽉1⽇, ⾄ 2022年3⽉31⽇
フィリピン・タール⽕⼭におけるその場観測に基づくマグマシステムの発達過程の研究, 国共同研究加速基,(国 共同研究&化(B)), 中区分17:地球惑星科学およびその関連分野, 名古屋⼤学, 熊⾕ 博之, ⾃ 2019年4⽉1⽇, ⾄ 2023年3⽉31⽇
2016

⼤沢 信⼆

噴気に由来する全物質分析による⽕⼭活動予測：⽔蒸気噴⽕の準備過程の解明, 基盤研究(C),
⼩区分17030:地球⼈間圏科学関連, 信州⼤学, 齋藤 武⼠, ⾃ 2020年4⽉1⽇, ⾄ 2023年3⽉31⽇
カルデラ湖の⽔質を⽤いた⼗和⽥⽕⼭活動モニタリング⼿法の開発, 基盤研究(C), ⼩区分
25030:防災⼯学関連, 秋⽥⼤学, 網⽥ 和宏, ⾃ 2019年4⽉1⽇, ⾄ 2022年03⽉31⽇
酸化還元状態を考慮した冷たい沈み込み帯での吸⽔・脱⽔過程と流体移動経路の解明, 基盤研
究(B), ⼩区分17040: 固体地球科学関連, 京都⼤学, 平島 崇男, ⾃ 2019年4⽉1⽇, ⾄ 2023年3⽉
31⽇
13楠本 成寿重⼒偏差テンソルを⽤いた実⽤的な断層傾斜⾓推定⼿法の研究開発, 基盤研究(C),⾃然災害科
学・防災学, 京都⼤学, 楠本成寿, ⾃ 2017年4⽉1⽇, ⾄ 2022年3⽉31⽇
重⼒異常や重⼒偏差データを⽤いた間隙率分布の推定⼿法についての研究, 基盤研究(C), ⼩区
分31020:地球資源⼯学およびエネルギー学関連, 京都⼤学, 楠本成寿, ⾃ 2020年4⽉1⽇, ⾄2023年3⽉31⽇
フィリピン・タール⽕⼭におる の場観測に基づくマグマシステムの発達過程の研究, 国際共同研究加速基⾦(国際共同研究&化(B)), 中区分17:地球惑星科学および の関連分野, 名古屋⼤学, 熊⾕ 博之, ⾃ 2019年4⽉1⽇, ⾄ 2023年3⽉31⽇

柴⽥ 智郎

⾼性能⽔素吸蔵合⾦を⽤いた質量分析装置内の⽔素除去に関する研究, 基盤研究(C), ⼩区分
17040:固体地球科学関連, 京都⼤学, 柴⽥ 智郎, ⾃ 2021年4⽉1⽇, ⾄ 2024年3⽉31⽇
パッシブサンプラーを⽤いた温泉⽔中溶存ガスの時空間分布の観測, 基盤研究(C), ⼩区分
63010:環境動態解析関連, 東京⼤学, ⾼畑直⼈, ⾃ 2018年4⽉1⽇, ⾄ 2022年3⽉31⽇

横尾 亮彦

⽕⼭噴気ジェット⾳響特性にとづいた新しい⽕⼭ガスSO2放出率推定⼿法の開発, 基盤研究
(C), ⼩区分17030:地球⼈間圏科学関連, 京都⼤学, 横尾 亮彦, ⾃ 2019年4⽉1⽇, ⾄ 2022年3⽉
31⽇

宇津⽊ 充

ドローン搭載型電磁探査による噴⽕発⽣場モニタリングと噴⽕切迫性評価, 基盤研究(A)中区分
25:社会システム⼯学、安全⼯学、防災⼯学および の関連分野, 神奈川県温泉地学研究所, 萬年
⼀剛, ⾃ 2021年4⽉5⽇, ⾄ 2026年3⽉31⽇
⽔蒸気噴⽕に関わる阿蘇⼭地下熱⽔系の観測的研究, 基盤研究(B), ⼩区分17040:固体地球科学関連, 神⼾⼤学, 南 拓⼈, ⾃ 2020年4⽉1⽇, ⾄ 2023年3⽉31⽇
⼩型無⼈機を⽤いた繰り返し空中磁気観測による⽕⼭活動モニタリング, 基盤研究(C), ⼩区分
25030:防災⼯学関連, 京都⼤学, 宇津⽊ 充, ⾃ 2019年4⽉1⽇, ⾄ 2023年3⽉31⽇

⽳井 千⾥

還元化学消磁を⽤いた⽇本の⽩亜系堆積層の古地磁気層序の研究磁気顕微鏡による地球内核形
成'後の地球磁場復元と地球⽣命史への影響の解明 ⼩区分17040:固体地球科学関連, 熊本⼤学,
渋⾕ 秀敏, ⾃ 2020年4⽉1⽇, ⾄ 2023年3⽉31⽇
琉球層群礁性⽯灰岩の古地磁気・岩⽯磁気分析による⾼分解能地球磁場・気候変動の復元, 国際共同研究加速基⾦(国際共同研究&化(B)), 中区分17:地球惑星科学および の関連分野, 国⽴
研究開発法⼈産業技術総合研究所, ⼩⽥ 啓邦, ⾃ 2020年4⽉1⽇, ⾄ 2024年3⽉31⽇
磁気顕微鏡による地球内核形成'後の地球磁場復元と地球⽣命史への影響の解明, 基盤研究(A),中区分17:地球惑星科学および の関連分野, 国⽴研究開発法⼈産業技術総合研究所, ⼩⽥ 啓邦,
14

⾃ 2021年4⽉1⽇, ⾄ 2025年3⽉31⽇

齋藤 圭

流⼊負荷を考慮した塩湖の化学的動態および⽔質形成年代の復元, 若⼿研究, ⼩区分04010:地理学関連, 京都⼤学, 齋藤 圭, ⾃ 2021年4⽉1⽇, ⾄ 2025年3⽉31⽇
2.6.

外部資⾦：その他 External Funds: Others

⼤沢 信⼆

海底カルデラ⽕⼭におけるマグマ起源流体湧出の地球化学的観測, 原⼦⼒規制委員会原⼦⼒規
制庁, 令和３年 原⼦⼒施設等防災対策等委託費（巨⼤噴⽕プロセス 等の知⾒整備に係る研究）⼆次受託研究 受託研究
別府市の地下⽔の現状に関する調査研究, 別府市, 受託研究楠本 成寿

空中重⼒偏差法データを⽤いた断層または地質構造境界の抽出及び傾斜⾓推定⼿法の研究, 受託研究, 住鉱資源株式会社, 受託研究, ⾃ 2021年4⽉15⽇, ⾄ 2021年9⽉15⽇
奨学寄附⾦, 株式会社 地球科学総合研究所, 寄附⾦

宇津⽊ 充を利⽤した在来型地熱貯留層の構造・状態推定, NEDO 国⽴研究開発法⼈ 新エネルギー・産業技術総合開発機構, NEDO 地熱発電導⼊拡⼤研究開発公募, ⾃ 2021年4⽉1⽇, ⾄ 2024年3⽉
31⽇

AI

15

．教育 Education

3

3.1.

担当科⽬ Teaching Subjects

⼤倉 敬宏

セミナー：阿蘇で観る⼤地の営み前期, 全学共通科⽬
固体地球物理学A後期, 理学部
地球惑星科学題演習ＤＡ前期, 理学部
地球惑星科学題演習Ｄ後期, 理学部
地球惑星科学題研究Ｔ３通年, 理学部
地球の物理前期, 全学共通科⽬
フィールド地球科学前期, 全学共通科⽬
地震学ゼミナール Ａ前期, 理学研究科
地震学ゼミナール Ｂ後期, 理学研究科
地震学ゼミナール 前期, 理学研究科
地震学ゼミナール Ｄ後期, 理学研究科
⽕⼭物理学前期, 理学部
⽕⼭物理学・⽕⼭流体学ゼミナールＡ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＢ後期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナール前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＤ後期, 理学研究科
⽕⼭物理学・⽕⼭流体学Ａ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学Ｂ後期集中, 理学研究科
観測地球物理学演習Ａ前期集中, 理学部

ILAS

⼤ 信⼆

地球テクトニクスセミナーIIＡ前期, 理学研究科
地球テクトニクスセミナーIIＢ後期, 理学研究科
地球テクトニクスセミナーIＡ前期, 理学研究科
地球テクトニクスセミナーIＢ後期, 理学研究科
地球惑星科学題演習 ＤＤ後期, 理学部
地球惑星科学題研究Ｔ通年, 理学部
地球熱学前期, 理学部
地球熱学・地熱流体学Ａ前期, 理学研究科
地球熱学・地熱流体学Ｂ後期, 理学研究科
地球熱学・地熱流体学ゼミナールＡ前期, 理学研究科
地球熱学・地熱流体学ゼミナールＢ後期, 理学研究科
地球熱学・地熱流体学ゼミナール前期, 理学研究科
地球熱学・地熱流体学ゼミナールＤ後期, 理学研究科
16

多階層地球変動科学特論：宇宙地球化学後期, 理学研究科
⽔圏地球物理学ゼミナールＩＩＩＡ前期, 理学研究科
⽔圏地球物理学ゼミナールＩＩＩＢ後期, 理学研究科
⽔圏地球物理学ゼミナールＩＩＩＣ前期, 理学研究科
⽔圏地球物理学ゼミナールＩＩＩＤ後期, 理学研究科
観測地球物理学演習Ｂ前期集中, 理学部
陸⽔学前期集中, 理学部

楠本 成寿

固体地球物理学ゼミナールＡ前期, 理学研究科
固体地球物理学ゼミナールＢ後期, 理学研究科
固体地球物理学ゼミナールＣ前期, 理学研究科
固体地球物理学ゼミナールＤ後期, 理学研究科
地球惑星科学題演習 ＤＣ後期, 理学部
地球熱学前期, 理学部
地球熱学・地熱流体学Ａ前期, 理学研究科
地球熱学・地熱流体学Ｂ後期, 理学研究科
地球熱学・地熱流体学ゼミナールＡ前期, 理学研究科
地球熱学・地熱流体学ゼミナールＢ後期, 理学研究科
地球熱学・地熱流体学ゼミナールＣ前期, 理学研究科
地球熱学・地熱流体学ゼミナールＤ後期, 理学研究科
測地学前期, 理学部
測地学ゼミナールＡ前期, 理学研究科
測地学ゼミナールＢ後期, 理学研究科
測地学ゼミナールＣ前期, 理学研究科
測地学ゼミナールＤ後期, 理学研究科
測地学・地殻変動論Ａ前期, 理学研究科
観測地球物理学演習Ｂ前期集中, 理学部

古川 善紹

グロ－バルテクトニクス前期, 理学部
⽕⼭物理学・⽕⼭流体学ゼミナールＡ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＢ後期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＣ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＤ後期, 理学研究科

柴⽥ 智郎

地球惑星科学題演習 ＤＣ後期, 理学部
地球惑星科学題研究Ｔ２通年, 理学部
17地球熱学前期, 理学部
地球熱学・地熱流体学Ａ前期, 理学研究科
地球熱学・地熱流体学Ｂ後期, 理学研究科
地球熱学・地熱流体学ゼミナールＡ前期, 理学研究科
地球熱学・地熱流体学ゼミナールＢ後期, 理学研究科
地球熱学・地熱流体学ゼミナールＣ前期, 理学研究科
地球熱学・地熱流体学ゼミナールＤ後期, 理学研究科
地球の物理前期, 全学共通科⽬
⽔圏地球物理学ゼミナール Ａ前期, 理学研究科
⽔圏地球物理学ゼミナール Ｂ後期, 理学研究科
⽔圏地球物理学ゼミナール Ｃ前期, 理学研究科
⽔圏地球物理学ゼミナール Ｄ後期, 理学研究科
統合科学：閉じた地球で⽣きる（地球 とエネルギー後期, 全学共通科⽬
観測地球物理学演習Ｂ前期集中, 理学部
陸⽔学前期集中, 理学部

横尾 亮彦

セミナー：阿蘇で観る⼤地の営み前期, 全学共通科⽬
固体地球物理学ゼミナールＡ前期, 理学研究科
固体地球物理学ゼミナールＢ後期, 理学研究科
固体地球物理学ゼミナールＣ前期, 理学研究科
固体地球物理学ゼミナールＤ後期, 理学研究科
地球惑星科学題演習 ＤＣ後期, 理学部
地球惑星科学題研究Ｔ３通年, 理学部
⽕⼭物理学前期, 理学部
⽕⼭物理学・⽕⼭流体学ゼミナールＡ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＢ後期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＣ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＤ後期, 理学研究科
⽕⼭物理学・⽕⼭流体学Ａ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学Ｂ後期集中, 理学研究科
観測地球物理学前期, 理学部
観測地球物理学演習Ａ前期集中, 理学部
ILAS

宇津⽊ 充

セミナー：阿蘇で観る⼤地の営み前期, 全学共通科⽬
地球惑星科学題演習 ＤＣ後期, 理学部
地球惑星科学題研究Ｔ通年, 理学部
応⽤地球電磁気学ゼミナールＡ前期, 理学研究科

ILAS

18

応⽤地球電磁気学ゼミナールＢ後期, 理学研究科
応⽤地球電磁気学ゼミナール前期, 理学研究科
応⽤地球電磁気学ゼミナール後期, 理学研究科
⽕⼭物理学前期, 理学部
⽕⼭物理学・⽕⼭流体学ゼミナールＡ前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナールＢ後期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナール前期, 理学研究科
⽕⼭物理学・⽕⼭流体学ゼミナール後期, 理学研究科
観測地球物理学習Ａ前期集中, 理学部

澤⼭ 和貴

地球熱学・地熱流体学ゼミナールＢ後期, 理学研究科
地球熱学・地熱流体学ゼミナール後期, 理学研究科
3.2.

学位審査 Participations in PhD/Master’s Defense

⼤倉 敬宏

阿蘇⽕⼭2014-2015年活動期におけるマグマヘッドの深さの時間変化, ⽯井杏佳, 理学研究科（博⼠）, 副査 2022年3⽉
2016年熊本地震にる阿蘇⽕⼭周辺の地殻変動とマグマだまりへの影響, 菅原嵩史, 理学研究科（修⼠）, 主査 2022年2⽉
多項⽬観測データから⾒た2016年阿蘇⽕⼭爆発的噴⽕の発⽣過程, 樋⼝和也, 理学研究科（修⼠）, 主査 2022年2⽉
空振波形解析にる阿蘇⼭中岳第⼀⽕⼝底物質の粘性推定ー2015年10⽉23⽇マグマ⽔蒸気爆発
発⽣前を例にしー, ⼩澤佑⼈, 理学研究科（修⼠）, 副査 2022年2⽉

⼤ 信⼆

衛星搭載降⽔レーダの⼆周波観測を活⽤した⼭岳域降⽔検出⼿法の開発, 清⽔陸, 理学研究科（修⼠）, 副査 2022年2⽉

楠本 成寿

年熊本地震にる阿蘇⽕⼭周辺の地殻変動とマグマだまりへの影響, 菅原嵩史, 理学研究科（修⼠）, 副査 2022年2⽉
アジョイント法に基づくスロースリップ発⽣域の摩擦特性の推定, ⽥村旭輝, 理学研究科（修⼠）, 副査 2022年2⽉

2016

横尾 亮彦

阿蘇⽕⼭2014-2015年活動期におけるマグマヘッドの深さの時間変化, ⽯井杏佳, 理学研究科
19

（博⼠）, 主査 2022年3⽉
空振波形解析による阿蘇⼭中岳第⼀⽕⼝底物質の粘性推定ー2015年10⽉23⽇マグマ⽔蒸気爆発
発⽣前を例にー, ⼩澤佑⼈, 理学研究科（修⼠）, 主査 2022年2⽉
多項⽬観測データから⾒た2016年阿蘇⽕⼭爆発的噴⽕の発⽣過程, 樋⼝和也, 理学研究科（修⼠）, 副査 2022年2⽉
機械学習を⽤いた桜島に る⽕⼭性地震の分類, 永江航也, 理学研究科（修⼠）, 副査 2022年
2⽉

20

．⼤学運営 Administration

4

4.1.

全学運営（役職等） School management

⼤倉 敬宏

地球熱学研究施設 施設⻑
4.2.

部局運営（役職等） Faculty management

⼤倉 敬宏

地球熱学研究施設 施設⻑
理学研究科地球熱学研究施設運営協議会議員
理学研究科情報セキュリティ委員会委員
理学研究科技術部協議会委員
理学研究科専攻⻑会議構成員
危機管理委員会委員

⼤6 信⼆

理学研究科地球熱学研究施設運営協議会議員

楠本 成寿

理学研究科地球熱学研究施設運営協議会議員
理学研究科将来計画委員会委員

古川 善紹

理学研究科地球熱学研究施設運営協議会議員

柴⽥ 智郎

理学研究科地球熱学研究施設運営協議会議員
理学研究科⾃⼰点検・評価委員会委員
理学部・教育委員会委員
理学部・教育委員会 教務委員会委員
理学研究科環境・安全委員会委員

横尾 亮彦

理学研究科地球熱学研究施設運営協議会議員
理学研究科SACRA広報・社会連携部⾨委員会 広報⼩委員会委員
21

．学術・社会貢献 Academic, Social Contribution

5

5.1.

委員歴 Committee Members

⼤倉 敬宏

⾄ 現在, ⽂部科学省 科学技術・学術審議会測地学分科会 地震⽕⼭部会 臨時委員
⾄ 現在, 東京⼤学地震・⽕⼭噴⽕予知研究協議会委員
⾄ 現在, 熊本⼤学病院 外部評価委員
⾄ 現在, 東京⼯業⼤学地球惑星科学系 外部評価委員
⾄ 現在, 次世代⽕⼭研究・⼈材育成総合プロジェクト⽕⼭噴緊急観測検討作業部会 委員
⾄ 現在, ⽕⼭噴⽕予知連絡会委員
⾄ 現在, ⽕⼭噴⽕予知連絡会 ⽕⼭活動評価検討会委員
⾄ 現在, ⽕⼭噴⽕予知連絡会 あり⽅検討作業部会委員
⾄ 現在, 阿蘇⽕⼭防災会議協議会委員
⾄ 現在, 阿蘇⽕⼭ガス安全対策専⾨委員会委員
⾄ 現在, 阿蘇⽕⼭博物館九⽊⽂化財団学術専⾨委員会委員
⾄ 現在, 阿蘇ジオパーク推進協議会学術顧問

⼤6 信⼆

⾃ 2021年, ⾄ 現在, 委員, 神奈川県温泉地学研究所外部評価委員会
⾃ 2021年, ⾄ 現在, 評議員, ⽇本温泉科学会
⾄ 現在, 委員, 別府市環境保全審議会
⾄ 現在, 委員, ⼤分県温泉監視調査委員会
⾄ 現在, 理事, ⼤分県温泉調査研究会
⾄ 現在, 委員⻑, ⼤分県環境審議会温泉部会
⾄ 現在, 委員, 独⽴⾏政法⼈油天然ガス・⾦属鉱物資源機構地熱資源開発アドバイザリー委員会
⾄ 現在, 評議員, ⽇本温泉科学会
⾄ 現在, 客員教授, ⼤分⼤学減災・復興デザイン教育研究センター

楠本 成寿

⾃ 2013年, ⾄ 現在, 代議員, 物理探査学会
⾃ 2009年8⽉, ⾄ 現在, 「物理探査」編集委員, 物理探査学会
⾃2019年6⽉ ⾄ 2021年5⽉, 評議員, ⽇本測地学会
⾃2016年6⽉ ⾄ 2021年5⽉, 「測地学会誌」編集委員, ⽇本測地学会

柴⽥ 智郎

⾃ 2019年3⽉, ⾄ 現在, 「陸⽔物理学会誌」編集委員, 陸⽔物理学会
22澤⼭ 和貴⾃ 2017年, ⾄ 現在, ⽇本地熱若⼿ネットワーク専⾨部会 広報幹事, ⽇本地熱学会
⾃ 2020年, ⾄ 2021年, Secretary, Student Chapter (Kyushu University), Society of
Exploration Geophysics

齋藤 圭

⾃ 2021年4⽉, ⾄ 現在, 広報委員⻑, 法政⼤学地理学会
5.2.

学術貢献活動 Academic Contribution

齋藤 圭

法政⼤学地理学会 2021年第2回例会（シンポジウム）, パネル司会・セッションチェア等, ⼩原丈明, ⾃ 2021年12⽉．社会貢献活動 Social Contribution

5.3

齋藤 圭

⾝近な⽔環境の全国⼀⻫調査, 助⾔・指導, 調査担当, 国⼟交通省および全国⽔環境マップ実⾏
委員会, ⾃ 2010年6⽉, ⾄ 現在

23

．セミナー Seminars

6

⽕⼭物理学・⽕⼭流体学ゼミナール（地球熱学セミナー
年4⽉20⽇2021年4⽉27⽇
2021年5⽉11⽇
2021年5⽉24⽇
2021年6⽉8⽇
2021年6⽉22⽇
2021年6⽉29⽇
2021年7⽉6⽇
2021年10⽉12⽇
2021年10⽉26⽇
2021年11⽉30⽇
2021年12⽉7⽇
2021年12⽉21⽇
2022年1⽉11⽇
2022年1⽉25⽇
2021

成⽥翔平「熱⾚外線カメラ観測による阿蘇中岳第⼀⽕⼝の放熱率推定(2020-2021年)」
齋藤 圭「最新の温泉泉質データに基づく地下熱⽔流動経路の推定」
⽳井千⾥「熱変質から推定する⽕⼭灰中磁性鉱物の温度履歴」
令和 3 年度 KUS02+HKD04 合同研究集会
梁 熙俊「別府南部における不圧地下⽔位の⻑期変動」
楠本成寿「重⼒異常や重⼒偏差のパCースペクトルとフィルタリングについて」
柴⽥智郎「地下⽔位データから降⽔の影響を取り除く試み」
宇津⽊充「九重⽕⼭における繰り返し空中磁気観測について」
成⽥翔平「九重⼭の最近30年間の地殻変動」
南 拓⼈（神⼾⼤）「⽐抵抗構造から⽰唆される阿蘇⼭の熱⽔系モデルについて」
齋藤 圭「珪藻化⽯分類分析による塩湖の古環境復元」
横尾亮彦「最近の阿蘇⼭中岳第⼀⽕⼝内の地形変化」
⼤信⼆「カルデラ壁からの深部起源CO2の放出」
⼤倉敬宏「阿蘇⽕⼭へのVUI(Volcanic Unrest Index)導⼊の試み」
澤⼭和貴「⽐抵抗と弾性波速度を利⽤し 岩⽯⻲裂内の流体流動挙動の解明」

地球熱学・地球熱流体ゼミナール（⾦曜セミナー
年4⽉23⽇2021年4⽉30⽇
2021年5⽉7⽇
2021年5⽉14⽇
2021年5⽉21⽇
2021年5⽉28⽇
2021年6⽉11⽇
2021年6⽉25⽇
2021年7⽉2⽇
2021年7⽉9⽇
2021年7⽉16⽇
2021年10⽉15⽇
2021年10⽉29⽇
2021

⽯井佳「地震・空振シグナルを⽤い 8グ86-0深さ・⽕道内⾳速の同時推定」
⽯橋 桜「阿蘇2014年噴⽕前後における地下⽐抵抗分布の推移」
服部匠哉「地震波に起因する地下⽔中の⼆酸化炭素気泡の成⻑」
樋⼝和「阿蘇⽕⼭における8グ8⽔蒸気噴⽕および爆発的噴⽕のメカ2
ズム推定に向けて」
⼤久保美鈴「阿蘇カルデラ周辺の地殻変動について」
政本⾵⼈「酸化還元電位測定による海底温泉湧出検出の試み」
⼩澤佑⼈「2015年10⽉23⽇の8グ8⽔蒸気爆発発⽣前にみられ 空振パルスの波形解析」
伊藤良介「九重⽕⼭群における3D磁気インバージョン解析」
宇野幸希「噴⽯軌道の計算」
菅原嵩史「断層運動が阿蘇⽕⼭の8グ8 りに与える影響について」
⼩隼⼈「阿蘇⼭の⽕⼭性地震について」
瀧⽥菜友「地下⽔変化の確率的分析」
梨元 昴「熱⾚外カメラによて得られる温度分布の誤差較正」
服部匠哉「圧⼒変動に起因し 地下⽔中の⼆酸化炭素気泡の成⻑」
菅原嵩史「阿蘇⽕⼭の8グ8 りに対する2016年熊本地震の断層変位に
よる影響」
24

年11⽉5⽇
2021年11⽉12⽇
2021年11⽉26⽇
2021年12⽉3⽇

樋⼝和也「2016年10⽉8⽇阿蘇⽕⼭爆発的噴⽕の直前過程~初期過程について」
⼩澤佑⼈「2015年10⽉23⽇に阿蘇⼭で発⽣したマグマ⽔蒸気爆発の前にみられた空振パルスの波形解析」
⼩隼⼈「阿蘇⽕⼭にる噴⽕準備過程でのA型地震の時間変動」
瀧⽥菜友「地下⽔流動シミュレーションに基づいた⽔位の相互変化」
伊藤良介「スパース磁気インバージョンコードの作成 ーADMMの実装と実データへの適⽤ー」
宇野幸希「2015年9⽉14⽇の阿蘇⼭噴⽕にる噴⽯軌道の計算」
2021年12⽉10⽇
⽯橋 桜「阿蘇噴⽕前後にる地下⽐抵抗分布の推移」
2021年12⽉24⽇
梨元 昴「monoplottingを⽤いた熱⾚外画像の補正⽅法」
2022年1⽉7⽇
政本⾵⼈「強酸性（pH<3）熱⽔に適⽤可能な地化度計年活動期の活動推移」
2022年1⽉21⽇
瀧⽥菜友「地下⽔流動シミュレーションに基づいた⽔位の相互変化」
伊藤良介「スパース磁気インバージョンコードの作成 ーADMMの実装と実データへの適⽤ー」
2022年1⽉28⽇
菅原嵩史「2016年熊本地震にる阿蘇⽕⼭周辺の地殻変動とマグマ 
への影響」
樋⼝和也「多項⽬観測データから⾒た2016年阿蘇⽕⼭爆発的噴⽕の発⽣過程」
⼤久保美鈴「阿蘇カルデラ周辺の⻑期的沈降現象の変動源について」
2022年2⽉4⽇
⼩澤佑⼈「空振波形解析からみた 阿蘇⼭中岳第⼀⽕⼝底物質の粘性推定­2015年10⽉23⽇マグマ⽔蒸気爆発発⽣前を例にして­」
宇野幸希「2014年9⽉14⽇の阿蘇⽕⼭噴⽕にる噴⽯初速度の計算」
⼩隼⼈「阿蘇⽕⼭にいて噴⽕準備過程で発⽣するA型地震」
2021

25

7.

⼀般公開報告 Open House

京⼤ウィークス ⼀般公開報告（別府）京⼤ウィークス2021研究施設⼀般公開事業は、研究施設本館建物のライトアップ展⽰は例年とおり実施したが、施設公開と公開講演会は新型コロナウィルス感染症（COVID-19）の感染防⽌対策のた(、オンデマンドによるデジタルコンテンツとオンラインによる配信で実施した。
施設公開は2021年11⽉5⽇（⾦曜⽇）10時〜20時にオンデマンドによる登録有形⽂化財建物でる研究施設本館建物と研究内容を紹介した。3本の実験動画（計約13分）（ドーム形成とカルデラ形成、七輪マグマ実験，椀がけ実験）、7本の調査動画（計約26分）（蒸気の流量測定，温泉の流量測定、伽藍岳の温泉調査、温泉遊離ガスのサンプリング、海底湧出ガスの調査、海底湧出ガスの調査、噴気地での調査）、4本の施設紹介動画（計約8分）（研修施設紹介、屋上から⾒た景⾊と解説、研究施設の実験室の解説、研究施設の湯）を公開し、閲覧者は72名でった。公開講演会は、11⽉5⽇（⾦曜⽇）15時半〜16時半に実施し、当施設楠本成寿教授による「重⼒異常から'た中部九州の構造とテクトMクス」の講演がり、7名の参加者がった。また、11⽉5⽇（⾦曜⽇）と11⽉6⽇（⼟曜⽇）両⽇の18時30分〜20時に登録有形⽂化財建物でる研究施設本館建物のライトアップ展⽰を⾏った。⼤正時代のレンガ造りの建物内部とライトアップに浮か$⾚レンガの幻想的な⾵景が評判で撮影⽬的で訪れる⼈が多く、両⽇合わせ26名の⽅が参加した。
新型コロナウイルスにより実施の告知活動は、別府市報への掲載、研究施設のTームSージへの掲載の'でったが、多くの⽅の参加がった。

研究施設本館建物のライトアップ

26

オンデマンド施設公開

オンZインによる公開講演会の様⼦

京⼤ウィークス ⼀般公開報告（阿蘇）１．はじめに京⼤ウィークスの⼀環として、2021年7⽉30⽇〜8⽉1⽇の⽇間、「本館のZイトアップ+4m球を使⽤したダジック・アース」「⼀般⾒学会」「建物観覧会」をそれれ実施した。
初⽇のZイトアップでは、本館の照明設備に加えLED電灯を13基使⽤し、それと並⾏してダジック・アースを使ってプレート運動や⽕⼭にいて解説をおこなった。2⽇⽬の⼀般⾒学会では、本館および屋外に展⽰ブースを設け、研究活動の紹介や実験などを⾏った。3⽇⽬の建物観
覧会は、登録有形⽂化財に)指定されている本館の内覧会や講演会を実施した。
また、昨年に引き続き新型コロナウイルス感染拡⼤防⽌対策として、完全予約制を導⼊し、１時間当たりの⼊場者数制限（最⼤10名）を設け、館内でれ違いが起こらないよう動線を設定した。さらに、対応にたる職員および学⽣は、マスクおよびフェイスシールドを着⽤し、来場者に)⼊場時の体温測定やマスクの着⽤をお願いした。
２．内容● ポスター展⽰による研究内容およびBンターの歴史紹介・⽕⼭学の⼀般向け解説● 公開実験・体験
・「七輪マ<マ実験」
・ 「サーモトレーサー体験」
・「地震計を使った振動体験」
・ 「VRウィーヘルト地震計」など● 視覚的展⽰物
・「ガリチン地震計の展⽰」
・ 「伸縮計・⽔管傾斜計の模型展⽰」
・「低周波マイクロフ7ンの展⽰」 ・ 「阿蘇⽕⼭の噴出物薄⽚の展⽰」
・「ダジック・アース展⽰」
・ 「ウィーヘルト地震計の復元作業映像」
27

３．社会告知の⽅法● ポスター掲⽰● 個⼈向けダイレクトメー● 掲載依頼報道関係
朝⽇新聞・熊本⽇⽇新聞・⻄⽇本新聞・読売新聞・毎⽇新聞● Webページ関連
京都⼤学（京⼤ウィークス）・理学研究科・⽕⼭研究センター・地球熱学研究施設● 当⽇の様⼦掲載・放映
朝⽇新聞・熊本⽇⽇新聞・⻄⽇本新聞・テレビ熊本
４．当⽇の様⼦

本館ライトアップの様⼦

岩⽯薄⽚の解説および撮影

低周波マイクロフォンの展⽰

ウィーヘルト地震計

VR

ダジック・アースを使った解説

電気抵抗測定実験
28

ドローンとサーモカメラの展⽰

建物観覧会講義の様⼦

29

．来訪者 Visitors

8

別府
年 4⽉19⽇
5⽉13⽇
5⽉26⽇
7⽉14⽇
11⽉9⽇
11⽉11⽇
12⽉2⽇
12⽉21⽇2022年 1⽉13⽇
3⽉9⽇
3⽉11⽇
3⽉16⽇
3⽉25⽇
2021

別府市職員5名
⼤分県職員（⽣活環境部⾃然保護推進室）
森藤⽒（電⼒中央研究所）
⼤分県職員（⽣活環境部⾃然保護推進室）
⻄島⽒（九州⼤学）
⼤分県職員（⽣活環境部⾃然保護推進室）
照⼭⽒
NHK⼤分取材班
⼤分県職員（⽣活環境部⾃然保護推進室）
⼤分県職員（⽣活環境部⾃然保護推進室）
⻄島⽒（九州⼤学）
森藤⽒（電⼒中央研究所）
NHK⼤分取材班
宗⽒（広島⼤学）

阿蘇
年 4⽉17⽇
⽕⼭研本館⾒学、11名参加
4⽉27⽇
理学部施設安全掛3名
4⽉28⽇
理学部施設安全掛電気担当
5⽉10⽇
⽕⼭研究センター復旧記念式典
5⽉15⽇
⽕⼭研本館⾒学会、8名参加
5⽉16⽇-20⽇
南⽒（神⼾⼤学）
5⽉24⽇-25⽇
⼤⾒⽒（京⼤防災研） 橋本⽒（北海道⼤学）
7⽉7⽇-8⽇
畑⽒（東⼤ 震研）
7⽉31⽇-8⽉1⽇ ⼩⽥⽊⽒（京⼤ 磁気センター） ⾼⾕⽒（京⼤理学技術部）
8⽉4⽇
南阿蘇村教育委員会
8⽉19⽇
池辺⽒（阿蘇⽕⼭博物館）
8⽉21⽇
ジオガイド⾒学会
8⽉31⽇-3⽇
⾵間⽒（京都⼤学）
9⽉2⽇
宮縁⽒（熊本⼤学）
9⽉13⽇-17⽇
南⽒（神⼾⼤学）
9⽉15⽇
⻄村⽒（東北⼤学）
10⽉8⽇
⾦⼦⽒ほか学⽣8名（神⼾⼤学）
10⽉19⽇中⽥⽒（防災科研）ほか10名
10⽉16⽇
南阿蘇⼩学⽣他20⼈
10⽉22⽇中⽥⽒ほか2名（防災科研）
10⽉25⽇森⽥⽒（防災科研）
10⽉29⽇⼤島⽒ほか1名（北海道⼤学）
11⽉25⽇-26⽇
多⽥⽒（JAMSTE） 市原⽒（名古屋⼤学）
11⽉1⽇-4⽇
⼤島⽒（北海道⼤学）

2021

30

⽉⽇
⽉ ⽇
⽉ ⽇
⽉⽇ ⽇
⽉ ⽇ ⽇
⽉ ⽇ ⽇
⽉⽇
⽉⽇
⽉⽇ ⽇
⽉⽇
⽉ ⽇ ⽇
⽉ ⽇
⽉ ⽇
2022年 ⽉ ⽇
⽉ ⽇ ⽇
⽉ ⽇
⽉⽇
⽉⽇
⽉ ⽇
⽉ ⽇
⽉⽇
⽉ ⽇
⽉ ⽇
⽉ ⽇
⽉ ⽇
11 4
11 14
11 15
11 7 -11
11 22 -23
11 23 -27
12 6
12 7
12 7 -8
12 8
12 13 -14
12 15
12 18
1 5
1 11 -13
1 26
2 1
2 2
2 16
2 22
3 7
3 13
3 14
3 23
3 24

阿部⽒（京⼤理学技術部）
ジュニアジオガイド研修、⼩学⽣26名
⾓皆⽒ 伊藤⽒（名古屋⼤学）
STEP10・阿蘇集中講義
DC「マグマから噴⽕まで」受講学⽣7名（京都⼤学）
⽥中⽒ 他5名（北海道⼤学）
⼟井⽒（京⼤防災研）
平松⽒（気象庁）
王⽒（京⼤防災研）
⽵村名誉教授研究集会 松島⽒（防災研）ほか11名
宮縁⽒（熊本⼤学） 原⼝⽒（⼤阪市⽴⼤学）
王⽒ ⼟井⽒（京⼤防災研）
南阿蘇村2名、池辺⽒（阿蘇⽕⼭博物館）
森⽥⽒（産業技術総合研究所） 森⽒（東京⼤学）
環境省
理学部企画課
南島原市市議会議員
松島⽒（産業技術総合研究所）
⻄村⽒（東北⼤学）
⻑⾕⽒（地熱技術開発）
林⽥⽒ほか2名（同志社⼤学）
中尾⽒（⾃然公園財団）
⻄村⽒（東北⼤学）
⽳井⽒（⾼知⼤学）

31

．定常観測 Routine Observations

9

Aso Volcanological Laboratory Permanent Stations
Nakadake monitoring network
Seismic Stations (Short period): HNT, KAF, KAW, SUN, KAK, UMA, TAK2, KAE2
Seismic Stations (Broadband): HND, HNT, SUN, KSM, UMA, KAE2
Tiltmeters: HNT (water tilt 2-comp.), SUN, KSM
Extensometers: HNT (invar 3-comp.)
Microphone: HND, SUN, UMA, KAE2, KAF
Geomagnetic Stations(proton): C1, C3, W1
GPS: ASOM, HOND, MONT, OHJO, SUNS

Seismic, geodetic, geomagnetic and infrasound stations in the central part of Aso. Createdby editing GSI Tiles (Topographic maps).
Central Kyushu regional network
Seismic Stations: AVL, UCM, MAK, NBR, MKN, HDK, TAT, MGR, AKG, CJB, ASJ, HNY,

SKM, BEP, KRK, TSK, AMM, TRM (Online telemetry)

Seismic network in the central Kyushu.

32

．装置・設備 Instruments and Facilities

10

装置 Instruments【別府】【Beppu】
重⾦属アナライザー
Metal monitor
ラドン分析装置
Radon detector
溶存⽔素計
Dissolved hydrogen meter
CO probe
CO 計測器CO , CH fluxmeter
CO -CH フラックスメーター
イオンクロマトグラフ
Ion chromatography
ガスクロマトグラフ
Gas chromatography
⾃動滴定装置Automatic titration system
ピストンシリンダー型⾼圧発⽣装置
Piston cylinder type high pressure apparatus
VLF-EM法電磁探査装置
Very low frequency electromagnetic field
detector
プロトン磁⼒計
Proton magnetometer
接地抵抗計・⼤地抵抗率計
Sounding resistivity meter
⽔準測量システム（⾃動み）Leveling survey system (automatic reading)
【阿蘇】
【Aso】
阿蘇，九重⽕⼭連続地震観測システム
Continuous seismic monitoring system
for Aso and Kuju Volcanoes
地殻変動観測坑道
Observation tunnel for ground deformation
プロトン磁⼒計
Proton magnetometers
フラックス/ート磁⼒計
Fluxgate magnetometers
地磁気絶対測定システム
Geomagnetic absolute measurement system
傾斜計
Tiltmeters
可搬型地震計（帯域, 短周期）
Portable seismometers (broadband short period)
重⼒計
Gravimeters
地磁気地電流測定装置（広帯域型ULF，ELF，VLF型） Magneto-Telluric measurement system(broad-band type, ULF, ELF, and VLF-band)
光波測距儀
Electronic distance measurement system
⽔準測量システム（⾃動み）Leveling survey system (automatic reading)
2

2

2

4

2

4

設備 Facilities【別府】中央実験室：
噴気試料や地熱関連⽔試料，地熱関連鉱物試料の各種化学分析を⾏うための前処理(溶解，濾過，酸化・還元処理)を⾏う．
流体分析室：
地熱関連⽔試料のイオン分析(⾦属元素も含む)を⾏う．分析室にはイオンクロマトグラフィや分光光度計などが設置されている．
機器分析室：
前処理済のガスや溶液の分析を⾏う．分析室にはガスクロマトグラフィやラドン分析装置などが設置されている．
33

【阿蘇】
地下観測坑道（阿蘇⽕⼭地殻変動観測坑道）：
阿蘇中岳第⼀⽕⼝から南⻄1 kmの，地下30 mに設けられた，直⾓三⾓形の⽔平坑道で，1987
年度に竣⼯した．現在は，⽔管傾斜計（25m），伸縮計（20，25 m），短周期地震計，⻑周
期地震計，広帯域地震計，および強震計が設置されている．
⽕⼭研究センター構内地震観測システム：
⽕⼭研究センター構内では，従来から:リ>タイ:による地震観測を⾏ってきたが，平成13
年度に，ノイズ低減の為，約200 mのボーリング孔を4本掘削し，孔底に地震計を導⼊した．
これにより，S/N⽐は⼤幅に改善され，従来識別できなかった中岳の⻑周期微動が検出されるようになった．また，ボーリング0(を採取したこにより，研究センターの丘，⾼野尾⽻根（たかのおばね）⽕⼭について地質学的に新たな知⾒が得られつつある．これは，阿蘇中央⽕⼝丘の噴⽕史を研究する上で 貴重な資料である．

34

Institute for Geothermal Sciences
Graduate School of Science, Kyoto University

京都⼤学⼤学院理学研究科 附属地球熱学研究施設
Beppu Geothermal Research Laboratory
Noguchibaru, Beppu, Oita 874-0903
Japan
Telephone: +81-977-22-0713
Facsimile: +81-977-22-0965【 別府 】〒874-0903 ⼤分県別府市野⼝原
電話: 0977-22-0713
ファックス: 0977-22-0965
Homepage: http://www.vgs.kyoto-u.ac.jp

Aso Volcanological Laboratory
Kawayo, MinamiAso, Kumamoto869-1404, Japan
Telephone: +81-967-22-0022
Facsimile: +81-967-22-2153【 阿蘇（⽕⼭研究センター）】〒869-1404熊本県阿蘇郡南阿蘇村河陽
電話: 0967-22-0022
ファックス: 0967-22-2153
Homepage:http://www.aso.vgs.kyoto-u.ac.jp/

Front Cover Image: A strombolian eruption at the Nakadake 1st crater of Aso volcano in April 2015.
(Photo by A. Yokoo): 2015 4

表紙の写真

年 ⽉の阿蘇中岳第⼀⽕⼝の式噴⽕の様⼦（横尾亮彦撮影）

Editorial compilation by A. Yokoo and Y. Doi

令和４年９月発行
編集兼発行者

京都大学大学院理学研究科

附属地球熱学研究施設

別府〒874－0903 大分県別府市野口原
電話：0977－22－0713
ファックス：0977－22－0965

阿蘇（火山研究センター）〒869－1404 熊本県阿蘇郡南阿蘇村河陽
電話：0967－67－0022
ファックス：0967－67－2153
印刷所
ホープ印刷株式会社〒861－8007 熊本県熊本市北区龍田弓削１丁目４－12
TEL 096－338－0500 FAX 096－386－3001
E-mail：mb@hope-printing.co.jp http://www.hope-printing.co.jp/