2015 URSI-Japan Radio Science Meeting(URSI-JRSM 2015)
Tokyo Institute of Technology, Tokyo, Japan, September 3-4, 2015

Program and Abstracts
Book

Sponsored by:- The Institute of Electronics, Information and Communication Engineers (IEICE)Technically supported by:- International Union of Radio Science (URSI)
In cooperation with:- Japan Geoscience Union- Science Council of Japan- The Astronomical Society of Japan- The Institute of Electrical Engineers of Japan- The Laser Society of Japan- The Remote Sensing Society of Japan

Table of Contents
Message from URSI-JRSM 2015 General Chair ………..…………………………………………….. 1
Sponsorships …………………………………………………………………………………………... 3
Conference Organization ………………………………………………………...……………………. 3
Program and Time Table ………………………………………………………...…………………….. 5
List of Poster Papers …..………………………………………………………...…………………….. 8
Abstracts (Keynote Lectures, Special Lectures and Invited Papers) ……………………………….. 14
Biographical Sketches of Speakers …………………………………………………………………... 29
Abstracts (Contributed Poster Papers) …………………………..………………………………….. 37
Conference Events ……………………………………………………………...…………………. 100
Opening Ceremony ……………………………………………………….………………….. 100
Poster Session …………………………………………………………….………………….. 100
Exhibitions ……….……………………………………………………….………………….. 100
Technical Tour …………………………………………………………….………………….. 100
URSI-Japan Commission Business Meetings ……………………………….……………………. 101
Banquet……………………………………………………………………….…………………... 101

Closing Ceremony ….…………………………………………………….………………….. 101
URSI-JRSM 2015 Steering Committee Meeting (members only) ……………………………101
SPC Committee Meetings (members only) ……………………………….………………….. 101
Floor Guide ……………………………………………………………………...…………………. 102
List of Exhibitors……………..……………………………………………………...…………… 105

ISBN978-4-88552-299-4 C3855
Copyright © 2015 by The Institute of Electronics, Information and Communication Engineers (IEICE)All rights reservedMessage from URSI-JRSM 2015 General Chair
It is our great pleasure to welcome all the participants to the “2015
URSI-Japan Radio Science Meeting” (URSI-JRSM 2015) held at
O-okayama Campus of Tokyo Institute of Technology, Tokyo, Japan on
September 3-4, 2015.
The “URSI-Japan Radio Science Meeting” (URSI-JRSM) is aregional URSI conference organized by the Japan National Committee of
URSI, and is held in Japan. The URSI-JRSM provides a scientific forum
for radio scientists and engineers in Japan and in the countries of Asia. The
objective of the Conference is to review current research trends, present
new discoveries, and make plans for future research and special projects in all areas of radio science, and aparticular emphasis is placed on promoting URSI activities in the Asian countries. At present, there are only
6 URSI member countries in Asia (China, India, Japan, Korea, Singapore, and Taiwan) and hence, it isdesirable to create a cooperative relationship among the URSI national committees in these countries andenhance the visibility of URSI in the Asian region. Therefore the URSI-JRSM also aims at establishingclose ties among the URSI Member Committees in Asia.
The URSI-JRSM was held for the first time in Tokyo in September 2014, and it was a success with aseries of presentations including 3 Keynote Lectures and 10 Invited Talks by outstanding scientists fromJapan and the Asian countries. The URSI-JRSM 2015 is the second URSI-JRSM, and the Conference has
been expanded extensively from the first edition from various aspects. The URSI-JRSM 2015 is sponsoredby The Institute of Electronics, Information and Communication Engineers (IEICE) and technicallysupported by the International Union of Radio Science (URSI), in cooperation with various academicorganizations including, Japan Geoscience Union, Science Council of Japan, The Astronomical Society ofJapan, The Institute of Electrical Engineers of Japan, The Laser Society of Japan, and The Remote Sensing
Society of Japan. Subject areas covered by the URSI-JRSM 2015 are broad including topics related to the
URSI Commissions A-K. The Conference consists of an oral session with 2 Keynote Lectures, 2 Special
Lectures, and 10 Invited Papers by prominent scientists in the radio science community as well as a poster
session with a total of 62 contributed papers. In addition to the scientific program, we are also organizing
the URSI-Japan Commission Business Meetings during lunch breaks on September 3 and 4, which areopen to all the Conference participants. For those of you who are interested in the Commission activities inJapan, please plan on attending these business meetings.
There are several unique features at this year’s URSI-JRSM. We have organized the “Student Paper
Competition” (SPC) as a special program for young scientists, which is similar to the program arranged atthe “URSI General Assembly and Scientific Symposium” (URSI GASS) and the “URSI Atlantic Radio
Science Conference” (URSI AT-RASC). Our SPC program is aimed at encouraging and supporting
full-time university students in degree programs so that they will be able to contribute to the URSI
community actively. We have received a total of 13 applications for the SPC program. The SPC review-1-

process is based on comprehensive evaluations of full-length papers submitted by the applicants and theirposter/oral presentations. On Thursday afternoon, September 3, we will select the three finalists after poster
presentations by the applicants, and these finalists will be announced at the Banquet held on the same day.
Subsequently on Friday morning, September 4, we will judge oral presentations by the three finalists, andselect the First, Second, and Third Prize winners. The selection results will be announced at the Award
Ceremony held on Friday afternoon, September 4. This SPC program is financially supported by URSI
central funds, and the support coming from URSI is greatly acknowledged.
I would also like to mention that we are organizing two special issues based on the papers presented at
the URSI-JRSM 2015. One of them is the “Special Issue of the 2015 URSI-Japan Radio Science Meeting”,
which is planned to be published in the journal “Radio Science” in 2016. This special issue will provide acollection of papers presented at the URSI-JRSM 2015, and will include invited papers (original andreview articles) and contributed papers (original articles). It is aimed at reporting progress and recent
advances in the radio science research ranging from fundamental theories to various applications, andcovers the areas related to all the URSI Commissions A to K. For those of you who are interested in thespecial issue in Radio Science, please submit a paper based on your presentation at URSI-JRSM 2015. The
other special issue is the “Special Issue on the URSI-JRSM 2015 Student Paper Competition”, which isscheduled for publication in the March 2016 issue of the URSI journal “Radio Science Bulletin”. This
special issue will provide a collection of papers by the three winners of the SPC program.
We hope that the “URSI-Japan Radio Science Meeting” (URSI-JRSM) will serve as one of the major,
regional URSI conferences in the international radio science community. We would be happy to receive anysuggestions on future directions of the URSI-JRSM conferences from all of you attending the URSI-JRSM
2015.
On behalf of the URSI-JRSM 2015 Committees, I would like to express a hearty welcome to all of
you working in the area of radio science for your participation in the URSI-JRSM 2015. After the scientificprogram of the first day is over, the Banquet will be held on campus of Tokyo Institute of Technology in theevening of September 3. We wish that you enjoy the Conference itself as well as the Banquet, and have apleasant stay in Tokyo. We are hoping that the URSI-JRSM 2015 will lead to a great success and the
URSI-JRSM will be held regularly in future.
Welcome to Tokyo, and welcome to the URSI-JRSM 2015!
September 3, 2015

Kazuya Kobayashi
General Chair, 2015 URSI-Japan Radio Science Meeting (URSI-JRSM 2015)
President, Japan National Committee of URSIURSI Assistant Secretary-General (AP-RASC)
Vice-Chair, URSI Commission B-2-

Sponsorships
This conference is sponsored by:- The Institute of Electronics, Information andCommunication Engineers (IEICE)Technically supported by:- International Union of Radio Science (URSI)
In cooperation with:- Japan Geoscience Union- Science Council of Japan- The Astronomical Society of Japan- The Institute of Electrical Engineers of Japan- The Laser Society of Japan- The Remote Sensing Society of Japan(in alphabetical order)

Conference Organization
General Chair
Kazuya Kobayashi, Chuo University (President,
Japan National Committee of URSI; URSI
Assistant Secretary-General (AP-RASC);
Vice-Chair, URSI Commission B)
Steering Committee
Chair:
Kazuya Kobayashi, Chuo University (President,
Japan National Committee of URSI; URSI
Assistant Secretary-General (AP-RASC);
Vice-Chair, URSI Commission B)
Secretary (General Affairs and Technical
Program):Satoshi Yagitani, Kanazawa University(Secretary, Japan National Committee of
URSI)
Secretary (Local Arrangements and Publicity):
Jun-ichi Takada, Tokyo Institute of Technology(Assistant Secretary, Japan National
Committee of URSI)
Secretary (Finance):
Tsuneki Yamasaki, Nihon University (Assistant
Secretary, Japan National Committee of
URSI)

-3-

Members:
Mitsuru Musha, The University of ElectroCommunications (Commission A Official
Member, Japan National Committee of
URSI)
Toru Uno, Tokyo University of Agriculture andTechnology (Commission B Official
Member, Japan National Committee of
URSI)
Nobuyoshi Kikuma, Nagoya Institute ofTechnology (Commission C Official
Member, Japan National Committee of
URSI)
Katsutoshi Tsukamoto, Osaka Institute ofTechnology (Commission D Official
Member, Japan National Committee of
URSI)
Osami Wada, Kyoto University (Commission E
Official Member, Japan National Committeeof URSI)
Yasuyuki Maekawa, Osaka ElectroCommunication University (Commission F
Official Member, Japan National Committeeof URSI)
Mamoru Yamamoto, Kyoto University(Commission G Official Member, Japan
National Committee of URSI)
Yoshiya Kasahara, Kanazawa University(Commission H Official Member, Japan
National Committee of URSI)
Kenta Fujisawa, Yamaguchi University(Commission J Official Member, Japan
National Committee of URSI)
Koichi Ito, Chiba University (Commission K
Official Member, Japan National Committeeof URSI)
Makoto Ando, Tokyo Institute of Technology(URSI Vice-President)
Yasuhiro Koyama, National Institute ofInformation and Communications
Technology (Chair, URSI Commission A)
Masao Taki, Tokyo Metropolitan University(Past Chair, URSI Commission K)
Yoshiharu Omura, Kyoto University (Past
Chair, URSI Commission H)
Members (Overseas):
Archana Bhattacharyya, Indian Institute ofGeomagnetism, India (President, India
National Committee of URSI)
Lou-Chuang Lee, Academia Sinica, Taiwan

(President, China (SRS) National Committeeof URSI)
Sangwook Nam, Seoul National University,
Korea (President, South Korea National
Committee of URSI)
Jian Wu, China Electronics Technology Group
Corporation, China (President, China (CIE)
National Committee of URSI)

Commission E:
Osami Wada, Kyoto University (Commission E
Official Member, Japan National Committeeof URSI)
Commission F:
Yasuyuki Maekawa, Osaka ElectroCommunication University (Commission F
Official Member, Japan National Committeeof URSI)
Commission G:
Mamoru Yamamoto, Kyoto University(Commission G Official Member, Japan
National Committee of URSI)
Commission H:
Yoshiya Kasahara, Kanazawa University(Commission H Official Member, Japan
National Committee of URSI)
Commission J:
Kenta Fujisawa, Yamaguchi University(Commission J Official Member, Japan
National Committee of URSI)
Commission K:
Koichi Ito, Chiba University (Commission K
Official Member, Japan National Committeeof URSI)

Technical Program Committee
Chair:
Satoshi Yagitani, Kanazawa University(Secretary, Japan National Committee of
URSI)
Members (in charge of specific Commissions):
Commission A:
Mitsuru Musha, The University of ElectroCommunications (Commission A Official
Member, Japan National Committee of
URSI)
Commission B:
Toru Uno, Tokyo University of Agriculture andTechnology (Commission B Official
Member, Japan National Committee of
URSI)
Commission C:
Nobuyoshi Kikuma, Nagoya Institute ofTechnology (Commission C Official
Member, Japan National Committee of
URSI)
Commission D:
Katsutoshi Tsukamoto, Osaka Institute ofTechnology (Commission D Official
Member, Japan National Committee of
URSI)

Secretariat
Prof. Tsuneki Yamasaki
College of Science and Technology,
Nihon University
1-8-14, Kanda-Surugadai, Chiyoda-ku,
Tokyo 101-8308
Tel/Fax：+81-3-3259-0771
E-mail：yamasaki@ele.cst.nihon-u.ac.jp-4-

Program and Time Table
Thursday, September 3, 2015
9:00

Opening Ceremony

9:15

Special Lecture 1
Chair: Kazuya Kobayashi (President, Japan National Committee of URSI;URSI Assistant Secretary-General (AP-RASC); Vice-Chair, URSI Commission B)
Makoto Ando (URSI Vice-President; Tokyo Institute of Technology, Japan)“International Union of Radio Science (URSI): Its Mission, Structure and Activities”

9:45

Invited Paper 1 (Commission A)
Chair: Mitsuru Musha (Commission A Official Member, Japan National Committee of URSI)
Nozomi Ohtsubo (National Institute of Information and Communications Technology, Japan)“Present Ion Optical Frequency Standards and New Challenges”

10:15

Coffee Break

10:45

Invited Paper 2 (Commission B)
Chair: Toru Uno (Commission B Official Member, Japan National Committee of URSI)
Toru Sato (Kyoto University, Japan)“Program of the Antarctic Syowa (PANSY) MST/IS Radar”

11:15

Invited Paper 3 (Commission F)
Chair: Yasuyuki Maekawa (Commission F Official Member, Japan National Committee of URSI)
Yoshio Karasawa (The University of Electro-Communications, Japan)“On Physical Limit of Wireless Data Transmission from Radiowave Propagation Viewpoint”

11:45

Special Lecture 2
Chair: Makoto Ando (URSI Vice-President)
Kenichi Iga (Professor Emeritus / Former President, Tokyo Institute of Technology, Japan)“Temperature-Insensitive Quartz Oscillator Enabled Highly Stable Communications
and Clocks – Review of Issac Koga’s Works –”

12:15
13:45

Lunch / Technical Tour / Commission Business Meetings
Keynote Lecture 1
Chair: Kazuya Kobayashi (President, Japan National Committee of URSI;URSI Assistant Secretary-General (AP-RASC); Vice-Chair, URSI Commission B)
Alexander I. Nosich (National Academy of Sciences of Ukraine, Ukraine)“Method of Analytical Regularization in Computational Photonics”-5-

14:30Invited Paper 4 (Commission K)
Chair: Masao Taki (Past Chair, URSI Commission K)
Koichi Ito (Chiba University, Japan)“Recent Medical Applications of Electromagnetic Waves”

15:00

Coffee Break / Exhibition

15:30

Poster Session / Exhibition

18:30

BanquetFriday, September 4, 2015
9:00

Invited Paper 5 (Commission D)
Chair: Katsutoshi Tsukamoto (Commission D Official Member, Japan National Committee of URSI)
Yuichi Kado (Kyoto Institute of Technology, Japan)“Multi-port Power Router and Its Impact on Future Smart Grid”

9:30

Invited Paper 6 (Commission G)
Chair: Mamoru Yamamoto (Commission G Official Member, Japan National Committee of URSI)
Toshitaka Tsuda (Kyoto University, Japan)“Study of Coupling Processes in the Solar-Terrestrial System”

10:00

Invited Paper 7 (Commission H)
Chair: Yoshiya Kasahara (Commission H Official Member, Japan National Committee of URSI)
Yohei Miyake (Kobe University, Japan)“Spacecraft-Plasma Interaction Effects on In-Space Electric Field Measurements”

10:30
11:00
12:15
13:45

Coffee BreakSPC (Student Paper Competition) Special Session
Lunch / Technical Tour / Commission Business Meetings
Keynote Lecture 2
Chair: Jun-ichi Takada (Assistant Secretary, Japan National Committee of URSI)
Ryuji Kohno (Yokohama National University, Japan / University of Oulu, Finland)“Harmonization Based on Regulatory Science for Scientific and Commercial Radio Uses”

14:30

Invited Paper 8 (Commission C)

Chair: Nobuyoshi Kikuma (Commission C Official Member, Japan National Committee of URSI)
Hiroyuki Tsuji (National Institute of Information and Communications Technology, Japan)“R&D of Mobile Communication Systems Using Millimeter Wave”-6-

15:00Coffee Break

15:30

Invited Paper 9 (Commission E)
Chair: Yasuhide Hobara (Commission E Member, Japan National Committee of URSI)
Yuichi Hayashi (Tohoku Gakuin University, Japan)“EM Information Leakage from Smart Devices”

16:00

Invited Paper 10 (Commission J)
Chair: Kenta Fujisawa (Commission J Official Member, Japan National Committee of URSI)
Masatoshi Ohishi (National Astronomical Observatory of Japan, Japan)“The Universe, Life and the Radio Sciences”

16:30

SPC Award Ceremony / Closing Ceremony

Conference Room Information:
Opening Ceremony …………………………………………. Multi-Purpose Digital Hall, West Building 9
Keynote Lectures, Special Lectures and Invited Papers …….. Multi-Purpose Digital Hall, West Building 9
Closing Ceremony …………………………………………. Multi-Purpose Digital Hall, West Building 9
Poster Session……………………………………….Media Hall / Collaboration Room, West Building 9

Exhibition ……………………………………………………………………Media Hall, West Building 9Technical Tour
Banquet………………………………………………..Tokyo Tech Museum, Centennial Hall………………………………………………………………………………………Student Hall

Commission Business Meetings …………………West Building 3 / West Building 9 (Details on page 101)-7-

List of Poster Papers
Commission A: Electromagnetic Metrology
A1.

Observation of the Magnetic-Insensitive Clock Transition in the 2S1/2(F=0)−2D3/2(F=2) Transitions inSingle 171Yb+
Yasutaka Imai, Kazuhiko Sugiyama and Masao Kitano (Kyoto University, Japan)

Commission B: Fields and Waves
B1.

The Physics of Angular Momentum RadioBo Thidé (Swedish Institute of Space Physics & Acreo Swedish ICT AB, Sweden); Fabrizio
Tamburini (Twist Off SRL and University of Padua, Italy)

B2.

Asymptotic-Numerical Analysis for Transient Scattered Field by a Coated Cylinder with a Thin
Lossy Medium Excited by a UWB Pulse Wave
Keiji Goto, Oki Okawa and Naoki Kishimoto (National Defense Academy, Japan)

B3.

Preliminary Study of Physical Optics for the Prediction of Scattering from Rough Surface inmmWaveRieko Tsuji, Kentaro Saito and Jun-ichi Takada (Tokyo Institute of Technology, Japan)

B4.

Energy Distribution of Air-hole Waveguides with Dielectric Cylinder Outside of the Defect Layer
Ryosuke Ozaki and Tsuneki Yamasaki (Nihon University, Japan)

B5.

Wiener-Hopf Analysis of the Plane Wave Diffraction by a Thin Material Strip
Takashi Nagasaka (Chuo University, Japan)

B6.

Maxwell-Schrödinger Hybrid Simulation for Optically Controlling Quantum States: A Two-Level
System Manipulated by a Light Pulse Pair
Takashi Takeuchi, Shinichiro Ohnuki and Tokuei Sako (Nihon University, Japan)

Commission C: Radiocommunication Systems and Signal Processing
C1.

Twin is Better than Single: A Hybrid Localization Method Using RSS and Time of Flight
Nopphon Keerativoranan (National Electronics and Computer Technology Center, Thailand);Chong-kwon Kim (Seoul National University, Korea)

C2.

Localization of Coherent Targets in a Multiple-Input Multiple-Output Radar Employing Spatial
Smoothing Between Virtual Antennas Extrapolated From Real Antennas
Ryo Nishikawa and Takehiko Kobayashi (Tokyo Denki University, Japan)

C3.

Human Motion Classification Utilizing Radio Signal Strength in Wireless Body Area Network
Sukhumarn Archasantisuk, Takahiro Aoyagi (Tokyo Institute of Technology, Japan); Tero
Uusitupa (Aalto University School of Electrical Engineering, Finland); Minseok Kim(Niigata University, Japan); Jun-ichi Takada (Tokyo Institute of Technology, Japan)

C4.

Localization of Illegal Radios Utilizing Phase-Difference as Location Fingerprints-8-

Azril Haniz, Gia Khanh Tran, Kei Sakaguchi and Jun-ichi Takada (Tokyo Institute ofTechnology, Japan); Daisuke Hayashi, Toshihiro Yamaguchi and Shintaro Arata (Koden
Electronics Co., Ltd., Japan)

Commission D: Electronics and Photonics
D1.

Resonant Tunneling Diode Receivers for 300-GHz-band Wireless Communications
Kousuke Nishio (Osaka University, Japan); Sebastian Diebold (Karlsruhe Institute ofTechnology, Germany); Shunsuke Nakai (Osaka University, Japan); Kazuisao Tsuruda,
Toshikazu Mukai and Jaeyoung Kim (Rohm Co, Ltd., Japan); Masayuki Fujita and Tadao
Nagatsuma (Osaka University, Japan)

D2.

Study on Hybrid RAN Using RoF and RoR for Distributed Small Cell Configurations
Katsutoshi Tsukamoto and Kazuo Kumamoto (Osaka Institute of Technology, Japan)

D3.

Terahertz Wireless Transmission Enabled by Photonics Using Binary Phase-shift Keying at 300 GHz
Yu Yasuda, Yusuke Fujita and Shintaro Hisatake (Osaka University, Japan); Shigeru Kuwano,
Jun Terada and Akihiro Otaka (NTT Corporation, Japan); Tadao Nagatsuma (Osaka
University, Japan)

Commission E: Electromagnetic Environment and Interference
E1.

Electromagnetic Imaging and Monitoring of Active Volcanoes: a Case Study and Future EM
Monitoring Plan for Mt. Kusatsu-Shirane Volcano, Japan
Yasuo Ogawa and Mineo Kumazawa (Tokyo Institute of Technology, Japan)

E2.

Development of an Upgraded BOLT System to Improve Location Accuracy
Yasuhiro Akiyama, Hiroshi Kikuchi, Ting Wu and Michael Stock (Osaka University, Japan);
Yoshitaka Nakamura (Kobe City College of Technology, Japan); Satoru Yoshida(Meteorological Research Institute, Japan); Tomoo Ushio (Osaka University, Japan)

E3.

Atmosphere and Ionosphere Connection as Revealed by the Network Observation of Very Low
Frequency Radio Signals
Sujay Pal and Yasuhide Hobara (The University of Electro-Communications, Japan)

E4.

A Statistical Study of Sub-ionospheric VLF Signal Anomaly Due to Geomagnetic Storms
Kenshin Tatsuta and Yasuhide Hobara (The University of Electro-Communications, Japan)

E5.

Statistical Spatial Distributions of Lightning with Charge Moment Changes Over Northern Area ofJapan by ELF and LLS Observations
Junpei Yamashita and Yasuhide Hobara (The University of Electro-Communications, Japan)

E6.

Modeling and Possible Determination of D-region Ionospheric Perturbation during Annular Solar
Eclipse on May 21, 2012, From VLF Signal Simulation and Multi-Propagation Path Observation
From UEC-VLF Network: Preliminary Results
Tamal Basak and Yasuhide Hobara (The University of Electro-Communications, Japan)

E7.

Bulk Current Injection Simulation Using Integrated Circuit Immunity Macro Model-9-Yosuke Kondo, Shinichiro Ueyama and Masato Izumichi (DENSO Corporation, Japan);
Osami Wada (Kyoto University, Japan)
E8.

Estimation of Packet Error Rate Considering Pulse Duration of Burst DisturbanceKazuhiro Takaya (NTT Corporation, Japan); Daisuke Tomita, Kouki Umeda, Tohlu
Matsushima, Takashi Hisakado and Osami Wada (Kyoto University, Japan)

E9.

Wideband Characterization of Multilayer Wave Absorber Using FDTD Method
Kosorl Thourn, Takahiro Aoyagi and Jun-ichi Takada (Tokyo Institute of Technology, Japan)

E10. A Novel Prediction for Very Low Frequency Transmitter Signal Amplitude Using NARX Neural
Network
Hendy Santosa (The University of Electro-Communications, Japan & University of
Bengkulu, Indonesia); Yasuhide Hobara (The University of Electro-Communications, Japan)

Commission F: Wave Propagation and Remote Sensing
F1.

Detections of Electromagnetic Pulses During and Prior to Earthquakes
Minoru Tsutsui (Kyoto Sangyo University, Japan)

F2.

A Study on Rain Area Motion along the Propagation Path and Ground Wind Velocity Around theEarth Station in Ku-Band Satellite Communications Links
Yasuyuki Maekawa and Yoshiaki Shibagaki (Osaka Electro-Communication University,
Japan);

F3.

Effects of Rain Area Motion on the Ku-Band Satellite Communications Links in Tropical Wet
SeasonAkihiro Tama, Yoshiaki Shibagaki and Yasuyuki Maekawa (Osaka Electro-Communication
University, Japan)

F4.

Effects of Rain Area and Rain Front Velocities on Rain Attenuation Characteristics in Ku-Band
Satellite Communications Links
Naoki Kubota, Yoshiaki Shibagaki and Yasuyuki Maekawa (Osaka Electro-Communication
University, Japan)

F5.

Effects of Rain Area Motion on the Ku-Band Satellite Communications Links in Tropical Dry
SeasonKeigo Takemoto, Yoshiaki Shibagaki and Yasuyuki Maekawa (Osaka
Electro-Communication University, Japan)

F6.

Channel Sounding in 2GHz and 60GHz Band for Multi-band Wireless LAN System
Tianyang Min, Kentaro Saito and Jun-ichi Takada (Tokyo Institute of Technology, Japan)

F7.

Application of Adaptive Digital Beam Forming for Polarimetric Phased Array Weather Radar
Hiroshi Kikuchi, Takuro Tashima, Ting Wu, Gwan Kim and Tomoo Ushio (Osaka University,
Japan); Hideto Goto and Fumihiko Mizutani (Toshiba Corporation, Japan)

F8.

Large Scale Subsurface Velocity Estimation with Array GPR System YAKUMO
Li Yi, Kazunori Takahashi and Motoyuki Sato (Tohoku University, Japan)- 10 -

F9.

Enhanced Clustering of Multipath Components Utilizing Geometrical Parameters for Dynamic
Indoor Double-Directional Propagation Channel
Panawit Hanpinitsak, Kentaro Saito and Jun-ichi Takada (Tokyo Institute of Technology,
Japan); Minseok Kim (Niigata University, Japan)

F10. Fundamental Study of GB-SAR Imaging by Compressive Sensing
Lilong Zou and Motoyuki Sato (Tohoku University, Japan)
F11. Leaf Shape Effects on the Scattering From Simplified Deciduous Leaf Structures Using Spherical
Wave Harmonics
Paul Jason Co and Jun-ichi Takada (Tokyo Institute of Technology, Japan)
F12. Experimental Evaluation of UWB Propagation in a Closed Box for Replacing Wired Interface Busesin Spacecrafts
Miyuki Hirose and Takehiko Kobayashi (Tokyo Denki University, Japan)

Commission G: Ionospheric Radio and Propagation
G1.

Measurement of Ionosphere over the Western Pacific Ocean
Mamoru Ishii, Hidekatsu Jin, Tatsuhiro Yokoyama, Takuya Tsugawa, Michi Nishioka andTakashi Maruyama (National Institute of Information and Communications Technology,
Japan)

G2.

Ionospheric Scintillation Observations by a Digital Beacon Receiver in Tromsø
Yasunobu Ogawa (National Institute of Polar Research, Japan); Yuichi Otsuka (Nagoya
University, Japan); Mamoru Yamamoto (Kyoto University, Japan)

G3.

Development of Digital Receiver for LF/MF Band Radio Wave Onboard Sounding Rocket
Takafumi Mizuno (Toyama Prefectural University, Japan); Yuki Ashihara (Nara National
College of Technology, Japan); Keigo Ishisaka (Toyama Prefectural University, Japan)

G4.

Electron Density Estimation by Using Simultaneous Observations of Lightning Optical Emissionsand Whistlers from ISS GLIMS Mission
Katsunori Suzuki, Yasuhide Hobara and Kanata Kakinuma (The University ofElectro-Communications, Japan); Ivan Linscott (Stanford University, USA); Umran Inan
(Koc University, Turkey); Mitsuteru Sato and Yukihiro Takahashi (Hokkaido University,
Japan); Tomoo Ushio and Zen Kawasaki (Osaka University, Osaka, Japan); Takeshi
Morimoto (Kinki University, Japan); Atsushi Yamazaki and Makoto Suzuki (Japan
Aerospace Exploration Agency, Japan)

G5.

VHF Lightning Observations and DOA Estimation From ISS / JEM-GLIMS
Takeshi Morimoto (Kinki University, Japan); Hiroshi Kikuchi and Tomoo Ushio (Osaka
University, Japan); Mitsuteru Sato (Hokkaido University, Japan); Atsushi Yamazaki and
Makoto Suzuki (Japan Aerospace Exploration Agency, Japan); Yasuhide Hobara (The
University of Electro-Communications, Japan)

G6.

Beacon Experiment of Low-latitude Ionosphere from Southeast Asia
Mamoru Yamamoto (Kyoto University, Japan)- 11 -

G7.

GNSS Scintillations and Ambient Ionization near the Northern EIA Crest Under Superstorm Eventof March, 2015
Debasis Jana and Shyamal Chakraborty (Raja Peary Mohan College, India)

G8.

Estimation of Electron Density Using the Propagation Characteristics of Radio Waves by S-520-29
Sounding Rocket
Keita Itaya and Keigo Ishisaka (Toyama Prefectural University, Japan); Yuki Ashihara (NaraNational College of Technology, Japan); Junichi Kurihara (Hokkaido University, Japan);
Takumi Abe (Japan Aerospace Exploration Agency, Japan)

Commission H: Waves in Plasmas
H1.

A Simple Kinetic-Fluid Model for Ion-Scale Waves in Collisionless Plasmas
Yasuhiro Nariyuki (University of Toyama, Japan); Shinji Saito and Takayuki Umeda(Nagoya University, Japan); Takuma Nakamura (Space Research Institute, Austrian Academyof Sciences, Austria)

H2.

Collective Thomson Scattering Diagnostics in Non-Equilibrium Plasma: Application to High Power
Laser Experiment
Shuichi Matsukiyo (Kyushu University, Japan); Youichi Sakawa (Osaka University, Japan);
Kentaro Tomita (Kyushu University, Japan); Yasuhiro Kuramitsu (National Central
University, Taiwan); Hideaki Takabe (Osaka University, Japan)

H3.

Dependencies of the Generation Process of Whistler-Mode Emissions on Temperature Anisotropy ofEnergetic Electrons in the Earth's Inner Magnetosphere
Yuto Katoh (Tohoku University, Japan); Yoshiharu Omura (Kyoto University, Japan); Yohei
Miyake and Hideyuki Usui (Kobe University, Japan); Hiroshi Nakashima (Kyoto University,
Japan)

H4.

Study on Improvement of Direction Finding Method for VLF Waves
Mamoru Ota, Yoshiya Kasahara and Yoshitaka Goto (Kanazawa University, Japan)

H5.

Effect of Electron Core-Halo Drift Speed on the Strahl Formation in the Solar Wind: Particle-in-cell
Simulations
Jungjoon Seough and Yasuhiro Nariyuki (University of Toyama, Japan); Peter Yoon
(University of Maryland, USA); Shinji Saito (Nagoya University, Japan)

H6.

Low Frequency Characteristics of a Wire Antenna aboard a Scientific Spacecraft
Tomohiko Imachi (Kanazawa University, Japan); Ryoichi Higashi (Ishiakwa National
College of Technology, Japan); Mitsunori Ozaki and Satoshi Yagitani (Kanazawa University,
Japan)

H7.

Effect of Lunar Surface Reflection on Polarizations of Auroral Kilometric Radiation Observed by
KAGUYA
Yoshitaka Goto, Ryota Kimura and Yoshiya Kasahara (Kanazawa University, Japan); Atsushi
Kumamoto (Tohoku University, Japan)

H8.

Development of Automatic Detection of Omega Signals Captured by PFX Onboard Akebono- 12 -I Made Agus Dwi Suarjaya, Yoshiya Kasahara and Yoshitaka Goto (Kanazawa University,
Japan)

Commission J: Radio Astronomy
J1.

Feasible Spectral Indices Measurement by Broad-band VLBI in VGOS Era
Kazuhiro Takefuji (National Institute of Information and Communications Technology,
Japan)J2.

Development of Wide-Band Bandwidth Synthesis Software for Geodetic VLBI
Tetsuro Kondo and Kazuhiro Takefuji (National Institute of Information andCommunications Technology, Japan); Yoshihiro Fukuzaki (Geospatial Information Authorityof Japan, Japan)

J3.

Development of Wideband Feed
Hideki Ujihara (National Institute of Information and Communications Technology, Japan)

J4.

New X-band Receiver System of Usuda 64m Antenna
Yasuhiro Murata, Masato Tsuboi, Kentaro Yamaguchi, Kenta Uehara, Hiroshi Takeuchi,
Kiyoshi Nakajima, Zen'ichi Yamamoto and Nanako Mochizuki (Japan Aerospace
Exploration Agency, Japan); Yusuke Kono, Masahiro Kanaguchi, Syunsaku Suzuki andTomoaki Oyama (National Astronomical Observatory of Japan, Japan)

J5.

Source Frequency Phase Referencing or 22 GHz H2O Masers and 43 and 86 GHz SiO Masers of the
Red Supergiant, S Per
Yoshiharu Asaki (National Astronomical Observatory of Japan, Japan)

J6.

Hours Timescale Radio Transients Detected by the Nasu Radio Telescope Array
Takahiro Aoki (Waseda University, Japan)

Commission K: Electromagnetics in Biology and Medicine
K1.

Microwave Induced Oxidative Stress Mediated Toxicity on Male Fertility Pattern and Possible
Protective Measures
Kavindra Kesari (University of Eastern Finland, Finland)

K2.

Temperature and Frequency Dependent on Electrical Properties of Swine Liver
Shohei Kon, Kazuyuki Saito and Koichi Ito (Chiba University, Japan)

K3.

Development of Biological Tissue Coagulating and Cutting Device Combining Microwave andRadio Frequency Current
Sho Suzuki, Yuta Endo, Kazuyuki Saito and Koichi Ito (Chiba University, Japan)

K4.

Study on Cytotoxicity of 0.07-0.30 THz Wave Exposure at Cultured Cells
Noriko Yaekashiwa (RIKEN, Japan); Shin'ichiro Hayashi (RIKEN, Japan); Kodo Kawase(Nagoya University & RIKEN, Japan)- 13 -

Abstracts(Keynote Lectures, Special Lecturesand Invited Papers)- 14 -

Keynote Lecture 1

Method of Analytical Regularization in Computational Photonics
Alexander I. Nosich
Laboratory of Micro and Nano Optics, Institute of Radio-Physics and Electronics NASU
vul. Proskury 12, Kharkiv 61085, Ukraine
Email: anosich@yahoo.com
Thanks to micro and nano technologies, the size of photonic-circuit components is now comparable andeven smaller than the length of the electromagnetic wave in the optical range. This makes traditional
ray-tracing methods of their analysis and design heavily inaccurate and calls for the full-wave study of theboundary-value problems (BVP) for the time-harmonic-Maxwell equations. Uniqueness of their solution isnormally guaranteed by a sufficient set of boundary, edge, and radiation conditions. Although one may try
to find it by using finite-difference discretizations of the partial differential equations, associated problemsof domain truncation, "good" exterior meshing, and solving enormous matrices are hardly compatible withhigh accuracy. To avoid these, BVP can be cast to an equivalent integral equation (IE). Then the radiation
condition is automatically satisfied and only the finite domains or boundaries need to be discretized.
Many textbooks and journal papers on computational electromagnetics deal with various IEs fordetermining the surface or polarization currents on metallic or dielectric scatterers, given the incident field.
Such IEs are obtained from the Green’s formulas and always have logarithmic or higher singular kernels;
many types of them are first-kind IEs. Additionally, some types of IEs loose the equivalence to BVP on
discrete sets of frequencies. IEs are frequently discretized for numerical solution with Method-of-Moments(MoM) using subdomain or entire-domain basis functions. Although this commonly brings meaningful anduseful results, unfortunately, there are not any theorems proving general MoM convergence, or even theexistence of exact solutions, for such IEs [1]. A rule-of-thumb of taking at least 10 mesh points per
wavelength is only a rule-of-thumb, and by no means does it guarantee any number of correct digits.
Meanwhile, there exists general approach to obtaining the second kind IEs of the Fredholm type, with asmooth enough kernel, from first-kind equations. Discretization of these new equations generates matrix
equations with condition numbers, which remain small if the truncation order is progressively increased.
The approach mentioned is collectively called the Method of Analytical Regularization (MAR) [2];sometimes semi-inversion is used as a synonym. It is based on the identification and analytical inversion of
the whole singular part of the original IE or its most singular component.
Theoretical merits of the MAR are numerous: exact solution exists, convergence to it is guaranteed, andeven asymptotic analytical solution is possible. Computationally, the MAR results in small matrix size for apractical accuracy, and sometimes no numerical integrations are needed for filling in the matrix. Thus, thecost of MAR algorithms is low in terms of both CPU time and memory. A frequent feature is that both
power conservation and reciprocity are satisfied at the machine-precision level, whatever is the number of
equations. The fact that the condition number is small and stable means that conjugate-gradient numericalalgorithms are very promising as emphasized in [3]. It is also important that the accuracy of the MAR is
uniform including sharp resonances known as whispering-gallery and lattice ones, which are inaccessible
by rough numerical approximations. By now, quite impressive number of optical-wave problems associatedwith dielectric, metal, semiconductor, and graphene objects has been analyzed with MAR, and interesting
effects have been studied in detail and with guaranteed accuracy [4].
References
1. D.G. Dudley, Electromagnetics, vol. 5, no 2-3, pp. 89-97, 1985.
2. A.I. Nosich, IEEE Antennas Propagat. Magazine, vol. 41, no 3, pp. 25-49, 1999.
3. V. Rokhlin, J. Computat. Physics, vol. 86, no 2, pp. 414-439, 1990.
4. LMNO IRE NASU Website, http://www.ire.kharkov.ua/LMNO/.

- 15 -

Keynote Lecture 2

Harmonization Based on Regulatory Science for Scientific and Commercial Radio Uses
Ryuji KOHNO
Ph.D., IEEE & IEICE Fellows
Professor, Graduate School of Engineering, Yokohama National University, Japan*
Finnish Distinguished Professor (FiDiPro), University of Oulu, Finland*** 79-5 Tokiwadai, Hodogaya-ku, Yokohama, 240-8501 Japan**Linnanmaa, P.O. Box 4500, FIN-90570 Oulu, FI-90014 Finland
Email: kohno@ynu.ac.jp, ryuji.kohno@oulu.fi
Harmonization for non-commercial scientific radio use such as radio astronomy, medical radio imaging
and commercial business radio use such as wireless communications and radio broadcasting must be suchan important inherent subject that URSI can discuss and contribute on global rule making from various
points of view by a fair manner. For these decades while demand for commercial radio uses has been
increasing, commercial radio uses have been always argued to harmonize with non-commercial scientificuses in radio regulation, because it seems to be a difficult problem how the same frequency of radio
resources can be shared by both scientific and commercial applications such as radio astronomy
observatory and wireless mobile communications. I have also committed in regional and international
regulatory committees. For instance, lately ultra wideband (UWB) radio in micro wave and millimeterwave bands have been quite attractive in commercial while radio astronomy observatory using the sameband may be interfered. When chairing the committees, it must be quite tough to prove a reasonable
judgment while keeping fair balance of long term scientific importance and urgent commercial benefits.
For this discussion on harmonization for non-commercial scientific and commercial business radio uses,
we have been taken care in some events of URSI, i.e. special session on URSI-ITU relations in
URSI-GASS2011, special workshop on harmonization of scientific and commercial radio uses in
AP-RASC2013. However, such argument has not been converged yet while increasing requirement to solve
such a problem. The reason why URSI is a best venue to discuss on the harmonization is obvious, because
URSI is community of mature representatives of both scientific and commercial radio uses. Scientific use
side could be represented by radio astronomy observatory, space agency (ETA, DASA, JAXA, NASA), andradio armature association in commissions F, G, H and J. In commercial use side, there are many
representatives from mobile cell phone operators, WiFi or Zigbee industry alliance, and radio broadcasting
association in commissions A, B, C, D, E, F, and H.
In this keynote speech in URSI-Japan Radio
Science Meeting (URSI-JRSM), I would like toapproach to resolve remained issues in theharmonization with “regulatory science,” that is amultidisciplinary concept and subject. Regulatory
science has been applied for medicine and medical
devices for instance because risk versus benefit to
apply newly invented medical devices to medical
healthcare services can be scientifically analyzed
with numerical evaluation. Then a rule or regulation for the devices can be made with scientific manner.
Regulatory body will examine the devices and systems compliant to the regulation. A plan of establishing a
regulatory science center by Science Council of Japan in Cabinet Office may be also introduced.
References
1. Ryuji Kohno, “Harmonization in Ultra Wide Band Regulation for both Scientific and Commercial
Radio Usages,” Special session on URSI-ITU relations in URSI-GASS2011, Istanbul, Aug.18, 2011
2. Ryuji Kohno, Special workshop in AP-RASC2013, Taipei, Sept. 4, 2013- 16 -

Special Lecture 1

International Union of Radio Science (URSI)
Its Mission, Structure and Activities
Makoto Ando
Department of Electrical and Electronic Engineering
Tokyo Institute of Technology
S3-19, 2-12-1, Oookayama, Meguro, Tokyo 152-8552, Japan
Email: mando@antenna.ee.titech.ac.jp
URSI: Mission and Structure
URSI is a non-governmental and non-profit union under the International Council for Science(ICSU).
URSI also consists of 44 member committees and 10 scientific commissions as his warp and woof. Its
mission, as appears on the Web, is to stimulate and co-ordinate, on an international basis, studies, research,
applications, scientific exchange, and communication in the fields of “radio science” or all aspects ofelectromagnetic fields and waves. Two important objectives are served by the following structures.
1. to represent radio science to the general public and to disseminate the studies and activities in radio
science and its applications regionally and globally. Currently, full and associate membershipcommittees are 44 countries/regions.
2. to stimulate and co-ordinate studies of the scientific aspects of electromagnetic waves, that is,
telecommunications using guided and unguided, the generation, emission, radiation, propagation,
reception, and detection of fields and waves, and the processing of the signals embedded in them.
Currently 10 scientific commissions cover the topics in radio science.
URSI: Conference Activities
Historically, URSI has general assemblies (GA) since 1922 and currently every triennial after 1954. GA
circulates over the world and the latest one was in 2014 in Beijing China. In GA, the council meeting withmembership committees as well as the business meetings with 10 commissions are held and the research
trends, new discoveries during the triennium are exchanged and the future planning in the next triennium
are discussed. As one of the URSI mission (1), dissemination into developing regions in Asia, Africa andsouth America should be strengthened. AP-RASC is the regional meeting established in 2001 enhance thevisibility of URSI in the Asian countries.
As another mission of URSI (2), since each 10 commission has their own commission specificmeetings and academic communities, the unique and distinguishing feature of URSI is to focus upon theinterdisciplinary studies between the multi commissions. In addition, to encourage young scientists tocontribute to various URSI activities is also an important task of URSI.
URSI: Latest Initiatives
Flagship meetings The Board decided that URSI should focus on yearly URSI “Flagship meetings”.
These URSI Flagship meetings should allow to maintain the momentum created at the GASS during thetriennium. The coming triennium (20115-2017) will be a transition period during which the focus will beon the 2015 AT-RASC meeting and the 2017 GASS. AP-RASC will also be organized in its current format
while the transition to an URSI Flagship meeting will be discussed. The commission specific meetings,
EMTS of Commission B and ISSSE of Commission C and D for example, can still be held if theCommission wishes to do so. Early carrier representatives (ECR) The Board has decided to have two
ECR for each commission, which invite the young generation into URSI activities. Now each commission
has three officers, the Chair, Vice-Chair and ECR; another ECR will be appointed as the fourth officer.
References
1. URSI Website, http://www.ursi.org/en/home.asp, ICSU Website, http://www.icsu.org/
2. URSI AT-RASC2015 Website, http://www.at-rasc.com/
3. URSI AP-RASC2016 Website, http://aprasc2016.org/- 17 -

Special Lecture 2

Temperature-Insensitive Quartz Oscillator
Enabled Highly Stable Communications and Clocks- Review of Issac Koga's WorksKenichi IGATokyo Institute of Technology, Professor Emeritus /Former President
2-33-10, Tsukushino, Machida, Tokyo, Japan 194-0001
E-mail Address: kiga@k00.itscom.net
On April 3, 1933, Issac Koga of Tokyo Institute of Technology, reported quartz crystal cut having
the zero temperature coefficient of frequency; R-cut, as shown in Fig. 1. The crystal quartz plate based on
this invention was used for transmitters at first and then for clocks. Currently, this type oftemperature-insensitive quartz crystal oscillator has become indispensable for radio communication systems

and information systems both for professional and consumer electronics.

Fig. 1 Temperature coefficients vs. crystal cut angles
References[1] I. Koga, "Longitudinal Vibration of Short Circular Cylinders (in Japanese)", JIEE Japan, Vol.50, No.508,
pp.1209-1224, Nov. 1930.
[2] I. Koga, "Thickness Vibration of Piezoelectric Oscillating Crystal", Physics, Vol. 3, No.2, pp. 70-80,
Aug. 1932.
[3] I. Koga, "R-cut Quartz Oscillating Plates and Harmonic Oscillation (in Japanese)", Proc. of 2nd
Conference on Engineering, Electrical Engineering Section, No. 102, p. 170, April 1932.
[4] I. Koga and K. Ichinose, "Quartz Oscillating Plates with Small Temperature Coefficients for Short-Wave(in Japanese)", Proc. of 8th Joint Conference on Electrical Engineering, No.135, pp. 205-206, April 3,
1933.
[5] I. Koga and N. Takagi, "Piezoelectric Oscillating Quartz Plates with Temperature Coefficients less than
10-7 /℃ (in Japanese)", JIEE Japan, Vol. 54, No. 543, p.940, Oct. 1, 1933.

- 18 -

Invited Paper 1

Present ion optical frequency standards and new challenges
Nozomi Ohtsubo
National Institute of Information and Communications Technology
4-2-1 Nukui-Kitamachi, Koganei, Tokyo 184-8795, Japan
Email: ohtsubo@nict.go.jp
Optical frequency standards realize frequency uncertainty much smaller than those of microwave
frequency standards. Relative frequency uncertainty smaller than 10-17 has been reported with single-ion
optical clocks [1-3] as well as with optical lattice clocks [4]. The optical clocks have their origin in the
proposal by H. Dehmelt [5], in which a trapped ion in the alkaline-earth electron configuration (B+, Al+,
Ga+, In+, Tl+ ;ions of group IIIA) supplies a narrow optical transition free from the quadrupole shift andwith negligibly small black body radiation (BBR) shift. This feature enables an inaccuracy in the 10-18
level.
Attempts have been made to implement the single-ion optical clock using ions with the alkaline-metalelectron structure (Ca+, Sr+, Yb+, Hg+) instead of those ions in the IIIA group due to difficulty in generating
a continuous-wave vacuum ultraviolet (VUV) radiation to laser-cool and to observe the ions. Although theclock transitions of alkaline-metal-like ions are not free from the quadrupole shift and the BBR shift, suchprogresses as discovery of exceptionally small BBR shift [2], engineering of the trap drive frequency tonull the shifts [3] have brought the inaccuracies of the alkaline-metal-like-ion optical clocks down to the
10-18 level. Matureness of the quantum information technologies with trapped ions has enabled cooling anddetection of IIIA ions without using the VUV radiation. A small-scale quantum computation method called
quantum logic spectroscopy (QLS) has been employed to implement a 27Al+ optical lock with a relative
uncertainty of the 10-18 level by NIST group.
Stability of the optical clocks grows as N1/2, where N is the number of atoms. The stability of ion clocks
can be improved by simultaneous trapping of multiple ions, but this strategy excludes thealkaline-metal-like ions due to the quadrupole shift. The Al+ optical clock cannot be extended to multiple
ions by the limitation of present QLS. We found that the most promising candidate reachable with present
technologies is 115In+. In our scheme the ions is sympathetically cooled withlaser-cooled 40Ca+ ions in a linear trap. The clock transition (1S0-3P0) is
detected by electron shelving using the 1S0-3P1 transition. This optical clock
scheme is directly extended to the multi-ion clock. We have alreadydeveloped fundamental technologies for implementing the single-ion In+
clock including resonant photo-ionization, sympathetic cooling, ion sorting[6], deep UV generation with high efficiency, stabilization of the clock laser
and generation of the VUV radiation at 159 nm [7]. We will present thedetails of the experiment and the future prospects of extending the clock to
the multi-ion version.
Fig1. Energy diagram of In+ ion.
References
1. C. W. Chou, et. al., Phys. Rev. Lett. 104, 070802(2010)
2. N. Huntemann, M. Okhapkin, B. Lipphardt, S. Weyers, Chr. Tamm, E. Peik, Phys. Rev. Lett. 108,090801 (2012)
3. Alan A. Madej, Pierre Dubé, Zichao Zhou, John E. Bernard, and Marina Gertsvolf, Phys. Rev. Lett.
109, 203002 (2012)
4. I. Ushijima, M. Takamoto, M. Das, T. Ohkubo, H. Katori, Nature Photon. 9, 185 (2015)
5. H. Dehmelt, IEEE Trans. Instrum. Meas., IM-31, 83(1982)
6. K. Hayasaka, Appl. Phys B 107, 965(2012)
7. K. Wakui, K. Hayasaka, T. Ido, Appl. Phys. B 117, 957(2014)- 19 -

Invited Paper 2

Program of the Antarctic Syowa (PANSY) MST/IS Radar
Toru Sato
Department of Communications and Computer Engineering
Graduate School of Informatics
Kyoto University
Yoshida-Honmachi, Sakyo-ku, Kyoto 606-8501, Japan
Email: sato.toru.6e@kyoto-u.ac.jp

The PANSY radar is the first Mesosphere-Stratosphere-Troposphere/Incoherent Scatter(MST/IS) radar in the Antarctic region. It is a large VHF monostatic pulse Doppler radar
operating at 47 MHz, consisting of an active phased array of 1,045 Yagi antennas and an
equivalent number of transmit-receive (TR) modules with a total peak output power of 500
kW [1]. The construction of the radar system started at Syowa Station (69°00’S, 40°35’E) in2010, and after operations with partial systems, the full operation was started in March 2015.
Here we briefly review the project’s scientific objectives, technical descriptions, and thepreliminary results of observations made to date. The radar is designed to clarify the role ofatmospheric gravity waves at high latitudes, and to explore the dynamical aspects of uniquepolar phenomena such as polar mesospheric clouds (PMC) and polar stratospheric clouds(PSC). The katabatic winds as a branch of Antarctic tropospheric circulation and as animportant source of gravity waves are also of special interest. Moreover, strong and sporadicenergy inputs from the magnetosphere by energetic particles and field-aligned currents can bequantitatively assessed by the broad height coverage of the radar which extends from thelower troposphere to the upper ionosphere.
From engineering points of view, the radar had to overcome restrictions related to thesevere environments of Antarctic research, such as very strong winds, limited power
availability, short construction periods, and limited manpower availability. We resolved these
problems through the adoption of specially designed class-E amplifiers, lightweight and toughantenna elements, and versatile antenna arrangements with flexible electronic control.
Although the radar has so far been operating mainly with partial systems, we have alreadyobtained interesting results on the Antarctic troposphere, stratosphere and mesosphere, suchas gravity waves, multiple tropopauses associated with a severe snow storm in the troposphereand stratosphere [1], and polar mesosphere summer echoes (PMSE) and polar mesosphere
winter echoes (PMWE) [2].
References
1. K. Sato, M. Tsutsumi, T. Sato, T. Nakamura, A. Saito, Y. Tomikawa, K. Nishimura, M.
Kohma, H. Yamagishi, and T. Yamanouchi, Program of the Antarctic Syowa MST/IS Radar(PANSY), J. Atmos. Solar-Terr. Phys., Vol.105, No.DOI:10.1016/j.jastp.2013.08.022,
pp.2-15, 2014..
2. T. Nishiyama, K. Sato, T. Nakamura, M. Tsutsumi, T. Sato, M. Kohma, K. Nishimura, Y.
Tomikawa, M.K. Ejiri, and T.T. Tsuda, Height and time characteristics of seasonal anddiurnal variations in PMWE based on 1 year observations by the PANSY radar (69.0S,
39.6E), Geophys. Res. Let., Vol.42, No.7, pp.2100-2108, 2015.

- 20 -

Invited Paper 3

On Physical Limit of Wireless Data Transmission
from Radiowave Propagation Viewpoint
Yoshio Karasawa
Advanced Wireless and Communication Research Center (AWCC)The University of Electro-Communications (UEC Tokyo)
1-5-1 Chofugaoka, Chofu, Tokyo 182-8585, Japan
Email: karasawa@ee.uec.ac.jp

Under the steady state condition with thermal noise, the physical limit of information transmission isgoverned by Shannon’s channel capacity theorem. However, in this formula, it does not contain factors ofradiowave propagation. In other words, for the ultimate information transmission, a sufficiently long time
for the coding and signal processing is expected
However, since wave propagation prevents its premise, there is another physical limit for informationtransmission in a different perspective with Shannon's channel capacity theorem. Even if the SN ratio issufficiently high, there is the limit for information transmission. This presentation deals with this matter
concerning physical limit of wireless transmission from a radiowave propagation viewpoint.
In wireless data transmission under multipath environment, a condition for not occurring bit errors due
to inter-symbol interference (ISI) is /Ts <<1 (: delay spread of multipath channel, Ts: symbol period ofdigitally modulated signal). Is this condition truly the physical limit? The answer is NO because it can besolved by employing a waveform equalizer or OFDM transmission scheme based on signal processing
using sufficiently long time period than Ts. On the other hand, a condition for not occurring bit errors due to
unexpected phase variations by Doppler frequency spread (f) is fTs <<1. Is this condition truly thephysical limit? The answer is NO because it can be solved by widening the signal bandwidth adopting suchas spread spectrum scheme.
However, for the avoidance, there are the following conditions. To release the condition of /Ts << 1,
it is necessary for keeping fTs <<1 so that there is no temporal variation during the signal processing
period. To release the condition of fTs <<1, it requires frequency-flat fading condition, namely, thecondition of /Ts << 1. Therefore, the physical limit is truly in the condition of f<<1 which iscomposed of two major propagation factors: delay spread () and Doppler spread (f). The image is shownin Fig. 1.
In the presentation, we will review this condition to taking future wideband wireless communication
systems into consideration. Finnaly, importance of collaboration among antennas (A), propagation (P) andsystems (S), namely, APS collaboration, will be stressed.

Fig. 1 Data transmission available/unavailable area

- 21 -

Invited Paper 4

Recent Medical Applications of Electromagnetic Waves
Koichi Ito
Center for Frontier Medical Engineering
Chiba University
1-33 Yayoi-cho, Inage-ku, Chiba-shi, 263-8522, Japan
Email: ito.koichi@faculty.chiba-u.jp

In recent years, various types of medical applications of electromagnetic waves have widely been
investigated and reported. Typical recent applications include:(1) Information transmission:- RFID (Radio Frequency Identification) tag system / Wearable or Implantable sensor- Wireless telemedicine / Mobile health system- Wireless capsule endoscopy(2) Diagnosis:- High intensity MRI (Magnetic Resonance Imaging)- Microwave CT (Computed Tomography) / Radiometry- Sleep monitor(3) Treatment:- Thermal therapy (Hyperthermia, ablation, etc)
- Transcranial magnetic stimulation (TMS)- Surgical device (Coagulation device, microwave knife, etc)
In this presentation, some practical medical applications of electromagnetic waves which have beenstudied in our laboratory are introduced. Firstly, a wearable dual-mode antenna [1] for vital data monitoring
systems is presented. A key technology for the antenna is body-centric wireless communications.
Secondly, an X-band antenna for a microwave sleep monitor is demonstrated with human phantom
experiments [2]. A “dynamic” phantom played an important role for the study. Thirdly, after a brief
description of thermal therapy and microwave heating, a coaxial-slot antenna and an array applicator
composed of several coaxial-slot antennas for minimally invasive microwave thermal therapy are
overviewed. A few results of actual clinical trials by use of coaxial-slot antennas are demonstrated from atechnical point of view. Then, as a new therapeutic application of coaxial-slot antennas, intracavitary
hyperthermia for bile duct carcinoma is briefly introduced [3]. Finally, a few different types of surgicaldevices using high power microwave energy, including a new coagulation device which can detect thecomplete coagulation, are introduced [4]. Heating characteristics of such microwave surgical devices areevaluated by numerical calculation as well as some experiments using phantoms, meat and animals.
References[1] C-H. Lin, and K. Ito, “Compact Wearable Dual-mode Antennas for Body-centric WirelessCommunications,” IEEE Int. Symp. Antennas and Propagat., pp.307-308, Memphis, USA, July 2014.
[2] J. Yonebayashi, et. al., “Evaluation on Performance of Doppler Radar for Breath Detection by Dynamic
Phantom,” AP-RASC’10, Toyama, Japan, Sep. 2010.
[3] K. Saito, et. al., “Estimation of Heating Performances of a Coaxial-Slot Antenna With Endoscope forTreatment of Bile Duct Carcinoma,” IEEE Trans. MTT, Vo. 54, No. 8, Aug. 2006.
[4] N. Ogasawara, et. al., “Development of Tissue Coagulation Forceps by Microwave Energy,” IEEE Int.
Symp. Antennas and Propagat., pp.723-724, Vancouver, Canada, July 2015.

- 22 -

Invited Paper 5

Multi-port Power Router and Its Impact on Future Smart Grid
Yuichi Kado
Graduate School of Science and Technology
Kyoto Institute of Technology
Matsugasaki, Sakyo-ku, kyoto 606-8585, Japan
Email: kado@kit.ac.jp
Through the efforts in achieving security of self-sufficient energy supplies and the prevention of global
warming, the effective utilization of electric energy is a critical issue. Therefore, at every stage from power
generation to power consumption, it is essential to develop the power electronics technologies that enableefficient power conversion and control. Especially, it’s necessary to establish the power interchanging
system that can flexibly interchange power among many energy resources, loads, and energy storage
devices. Several studies about grid systems have been reported [1], [2], but there is not enoughexperimental verification about power routers that perform the power interchanging function regardless of
DC and AC power and easily construct power networks.
Under such situations, we propose the multi-port power router that can be used as a unit cell to easily
build a power delivery system in order to meet many kinds of user’s requirements and reduce the systemdevelopment costs as shown in Fig. 1. The router is composed of a three-way isolated DC/DC converter asa core unit, AC/DC converter units, and a software-defined control unit. The router can distribute power inthree directions and perform the power interchanging function regardless of DC and AC. Furthermore,
many kinds of loads, energy resources, and energy storage devices are connectable to the routers by only
rewriting the control software.
For multiple energy routers installed at various locations to work cooperatively to implement flexible
power interchange in real time, as shown in Fig. 2, it is necessary to measure the current and voltage
amplitude at each node and reflect the results in the phase-shift modulation (PSM) control of each energy
router. We consider the communication requirements of a control framework for enabling dynamic
adjustments of power flow through the coordinated operation of multiple energy routers placed in different
locations. That requires highly accurate time synchronization of the control modules of the energy routers.
A strong candidate for the control network is an optical fiber network access system. When optical fiber isused for high-voltage power converter control, electrical isolation, robustness against lightning strikes, andexcellent robustness against electrical noise can be expected.

Fig.1. Application of the multi-port power router to future grid systems.

Fig.2. Control layers between energy routers

References
1. A. Huang, et al., “The Future Renewable Electric Energy Delivery and Management (FREEDM)
System: The Energy Internet,” Proceeding of the IEEE, Vol.99, pp. 134-148, 2011.
2. R. Abe, et al., “Digital Grid: Communicative Electrical Grids of the Future,” IEEE TRANSACTIONSON SMART GRID, Vol. 2, No. 2, pp. 399-410, 2011.

- 23 -

Invited Paper 6

Study of Coupling Processes in the Solar-Terrestrial System
Toshitaka Tsuda
Research Institute for Sustainable Humanosphere (RISH)
Kyoto University
Gokasyo, Uji, Kyoto, 611-0011, Japan
Email: tsuda@rish.kyoto-u.ac.jp"Coupling process in the solar-terrestrial system" aims to study the solar energy inputs into the Earth, andthe response of Geospace (magnetosphere, ionosphere and atmosphere) to the energy input. The solar energy
can mainly be divided into two parts - the solar radiation involving infra-red, visible, ultra-violet and X-ray,and solar wind which is a high-speed flow of plasma particles. The solar radiation becomes maximum atthe equator; and atmospheric disturbances are actively generated near the Earth's surface. They further excite
various types of atmospheric waves which propagate upward carrying energy and momentum. On the other
hand, the energy associated with the solar wind converges into the polar regions where disturbances are
generated. A part of the energy is transported toward lower latitudes and lower atmospheric regions. We
propose to establish large atmospheric radars with active phased array antenna at the equator and the Arctic
region. Among the equatorial regions, we focus on the Indonesian region where atmospheric disturbances
are most intense. We will establish a comprehensive observatory in Indonesia with the Equatorial MU (EMU)
radar as its main facility. Alongside, we will take part in the construction of the state-of-the-art radar, called
EISCAT_3D, in northern Scandinavia under international collaborations. We will also develop the global
observation network of portable equipment from the equator to both polar regions, and study the flow of theenergy and materials in the whole atmosphere.

- 24 -

Invited Paper 7

Spacecraft-Plasma Interaction Effects on In-Space Electric Field Measurements
Yohei Miyake
Graduate School of System Informatics
Kobe University
1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
Email: y-miyake@eagle.kobe-u.ac.jp
Space exploration and exploitation have been rapidly increasing, and a strong demand arises regarding comprehensive understanding of spacecraft-plasma (SP) interactions. This is clearly required
to ensure survivability and proper operations of space-based systems, and also for correct interpretationof observation data collected in situ by scientific satellites. In space environments, not only the spacecraft charging but also surrounding plasma disturbances such as caused by the photo-/secondary electron
emissions and sheath/wake formations may interfere directly with in-situ electric field measurements
using the double probes/electric antennas.
In DC to a low frequency range, the double probe instrument basically measures an electric field as
the potential difference between two probe positions. An asymmetric electrostatic environment creatednear the spacecraft due to sheath/wake formation and photo-/secondary electron emission is thus a potential source of errors in the field measurements. The near-spacecraft potential structure is relevant to
the process of spacecraft charging, and the significance of such spurious effects depends strongly on a
spacecraft potential.
In a higher frequency range, the response of an antenna to alternating electric field should be characterized by its effective length and RF impedance. Particularly, the antenna impedance in plasma is known
to differ from that in free space. Further difficulty is the plasma non-uniformity due to the formation of
ion or photoelectron sheaths around the antenna. Since a full treatment of the impedance in magnetizedplasma is quite complicated, a simplified model represented by the sheath capacitance and resistance isused commonly. Even using the simplified model, an actual sheath structure should be determined by thecomplex SP interactions. Hence, the self-consistent solution of such phenomena is necessary to assess
these effects.
We apply the plasma particle simulations based on a particle-in-cell (PIC) method to the problems.
Our original simulation code EMSES[1] is based on the standard electromagnetic PIC method, and also
has the capability to include the conducting bodies of a spacecraft, based on the capacitance matrix
method. In addition, a number of physics such as the photoelectron emission and the secondary emission
are modeled numerically in the latest version of EMSES. The code has been applied to some specificspacecraft, e.g., Geotail, Cluster[2, 3], BepiColombo/MMO[4], and Solar Probe Plus. In this paper, we
present selected examples of the EMSES applications and overview the SP interaction effects on in-spaceelectric field measurements.
References
1. Y. Miyake and H. Usui, “New electromagnetic particle simulation code for the analysis of spacecraftplasma interactions,” Phys. Plasmas, Vol. 16, p. 062904, 2009.
2. Y. Miyake, C. M. Cully, H. Usui and H. Nakashima, “Plasma particle simulations of wake formation
behind a spacecraft with thin wire booms,” J. Geophys. Res., Vol. 118, pp. 5681-5691, 2013.
3. Y. Miyake, Y. Nishimura and Y. Kasaba, “Asymmetric electrostatic environment around spacecraft in
weakly streaming plasmas,” J. Geophys. Res., in press.
4. Y. Miyake, H. Usui, H. Kojima and H. Nakashima, “Plasma particle simulations on stray photoelectron current flows around a spacecraft,” J. Geophys. Res., Vol. 117, p. A09210, 2012.

- 25 -

Invited Paper 8

R&D of Mobile Communication Systems Using Millimeter Wave
Hiroyuki Tsuji
Wireless Network Research InstituteNational Institute of Information and Communications Technology
4-2-1, Nukui-Kitamachi, Koganei, Tokyo 184-8795, Japan
Email: tsuji@nict.go.jp
Demand has been increasing for better mobile phone and wireless local area network (LAN) access forpeople on board aircraft or trains. Now, many airline companies have introduced in-cabin use of cellular
phones and Internet accesses with a system involving satellites.
The National Institute of Information and Communications Technology (NICT) and Mitsubishi Electrichave conducted the research and development of mobile communication systems using high-frequency bands
such as the millimeter wave band (40 GHz band), which can provide broadband communication, in order tomeet the recent demand for Internet connectivity and high-speed data communication on board fast-moving
vehicles such as aircraft and bullet trains [1]. The proposed communication system requires several key
technologies such as tracking the target terminal antennas. The ground-based tracking antenna has tocontinuously track the moving vehicle with a high degree of accuracy. Meanwhile, the onboard antenna hasto track the ground-based antenna by considering the moving vehicle attitude and its environment as well.
We, therefore, developed a high-accuracy mono-pulse tracking antenna as the ground-based antenna and adownsized active phased array antenna as the on-board antenna.
The ground tracking antenna with a mechanically controlled reflector (Fig. 1) has a feature that enables
it to direct its antenna beam in a specific direction by tilting a reflection disk mechanically and to provide a
cost-effective and power-efficient tracking antenna. As for the onboard antenna, any onboard antenna
equipment must be compact and lightweight. A tracking antenna that is electrically steered is considered
suitable for use on an onboard antenna because the phased-array antenna has no moving parts and it can also
be resized. We have succeeded in reducing the array antenna's weight by using plastic resin array elements(Fig. 2). We also carried out some experiments using actual airplanes and confirmed that over 100 Mbps
transmission speeds could be achieved between the ground and the air.
The results of this research were reflected in the new report ITU-R M.2282: "Systems for public mobile
communications with aircraft" in ITU-R SG5 WP5A [2].
We are also trying to deploy this technology on railways, and are promoting research and development
to provide broadband access to high-speed trains. To reflect this achievement, we have proposed a new ITUR report that summarizes the relevant technology in ITU-R SG5 WP5A, and are currently performing studies
for the completion of the report.

Figure 1. Ground tracking antenna with lens antenna
Figure 2. Onboard active phased array antenna.
References
1. H. Tsuji, M. Suzuki, and T. Morisaki, “Millimeter Broadband Wireless System Aims To Accommodate
Inflight Wireless Communications,” The 13th International Symposium on Wireless Personal
Multimedia Communications (WPMC 2010), 2010.
2. ITU-R Website, http://www.itu.int/md/R12-SG05-C-0072- 26 -

Invited Paper 9

EM Information Leakage from Smart Devices
Yuichi Hayashi
Department of Electrical Engineering and Information Technology, Faculty of Engineering
Tohoku Gakuin University
1-13-1 Chuo, Tagajo, Miyagi 985-8537, Japanyu-ichi@mail.tohoku-gakuin.ac.jp

At present, ICT devices at risk for EM information leakage include many commercial products, such as
CRT and LCD monitors [1, 2], printers [3], keyboards [4], cryptographic modules [5], and so on. In thispresentation, we focus on EM information leakage from smart devices (e.g. Tablet PCs, Smartphones). The
use of smart devices is spreading rapidly, and accordingly Users browsing and inputting personal
information in public spaces can often be seen by third parties. Unlike conventional mobile phones andnotebook PCs equipped with distinct input devices (e.g., keyboards), tablet PCs and smartphones havetouchscreen keyboards for data input. Such integration of display and input device increases the potential
for harm when malicious attackers capture the screen. This presentation introduces a new threat: a captureof tablet PC displays via measurement of electromagnetic (EM) emanation. In conventional studies, such
EM display capture has been achieved by using non-portable setups. Those studies also assumed that asignificant amount of time was available in advance of arrest to obtain the electrical parameters of thetarget display. In contrast, this presentation demonstrates that such EM display capture is feasible in realtime by a setup that fits in an attaché case. The EM display capture does not require the conventionalparameter tuning, relying instead on prior profiling and real-time signal processing. The attack distancemakes this method a practical threat to general tablet PCs in public places. We discuss possible attack
scenarios based on the setup described above. Also, we describe a mechanism of EM emanation from tablet
PCs and a countermeasure against such EM display capture.

References
1. M.G. Kuhn, “Compromising emanations of LCD TV sets,” IEEE International Symposium on
Electromagnetic Compatibility, pp. 931-936, Long Beach, California, USA, August, 2011.
2. H. Sekiguchi and S. Seto, “Estimation of receivable distance for radiated disturbance containing information
signal from information technology equipment,” IEEE International Symposium on Electromagnetic
Compatibility, pp.942-945, Long Beach, California, USA, August, 2011.
3. T. Tosaka, K. Taira, Y. Yamanaka, A. Nishikata, M. Hattori, “Feasibility study for reconstruction of
information from near field observations of the magnetic field of laser printer,” 17th International Zurich
Symposium on Electromagnetic Compatibility, 2006. EMC-Zurich 2006, pp.630-633, Singapore, February,
2006.
4. M. Vuagnoux and S. Pasini, “An improved technique to discover compromising electromagnetic
emanations,” IEEE International Symposium on Electromagnetic Compatibility, pp.121-126, Fort
Lauderdale, FL, USA, July 2010.
5. D. Agrawal, B. Archambeault, R. Rao, and P. Rohatgi, “The EM Sidechannel(s),” CHES 2002, Lecture
Notes in Computer Science, pp. 29-45, San Francisco, CA, August, 2002.

- 27 -

Invited Paper 10

The Universe, Life and the Radio Sciences
Masatoshi Ohishi
Astronomy Data Center, National Astronomical Observatory of Japan
2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
Email: masatoshi.ohishi@nao.ac.jp
It is widely accepted that prebiotic chemical evolution from small to large and complex molecules would
have resulted in the Origin of Life. On the other hand there are two conflicting views where inorganic
formation of complex organic molecules (hereafter COMs) occurred in the early Earth, on the Earth or outof the Earth. Ehrenfreund et al.[1] indicated that exogenous delivery of COMs by comets and/or asteroids tothe early Earth could be larger than their terrestrial formation by three orders of magnitude. If amino acids
are formed in interstellar clouds, significant amount of them may be delivered to planets. Detection of aminoacids would accelerate the discussion concerning the universality of “life”.
So far, many trials by means of radio telescopes to detect the simplest amino acid, glycine (NH2CH2COOH),
were made towards Sgr B2 and other high-mass forming regions, but none of them were successful. One idea
to overcome this situation would be to search for precursors to glycine. Although the chemical evolution ofinterstellar N-bearing COMs is not well known, methylamine (CH3NH2) is proposed as one precursor toglycine. CH3NH2 may be formed through hydrogenation to HCN on dust surface [2][3]: HCN → CH2NH→ CH3NH2. Since it is well known that CO2 exists in most of molecular clouds, glycine would be expected
where CH3NH2 is rich.
In April 2013, we made a survey to find CH2NH-rich sources byusing the Nobeyama 45 m radio telescope. We succeeded to detectfour new CH2NH sources. The derived fractional abundances of
CH2NH relative to H2 are as high as 6x10-8, implying that CH2NH
may exist widely in the Universe. Followed by this success, wemade a survey of CH3NH2 in 2014, since further hydrogenation to
CH2NH would efficiently produce CH3NH2. We succeeded todetect CH3NH2 towards two sources as shown in the Figure on theleft. The estimated fractional abundance of CH3NH2 to H2 was ~10-8, about 10 times higher than that reported towards SgrB2(N) [4].
This result would mean that interstellar glycine is highly expected
to exist in various sources. We would be able to find more
CH3NH2-rich sources by observing the 79 GHz region.
However, a serious problem may occur in the 76 to 81 GHz range.
Car manufacturers have strongly pushed to use this frequency range for automobile radars. Improving safety
is very important, however, unless an exclusion zone is established around each radio observatory operatingin the 76 to 81 GHz range, highly sensitive radio astronomical observations would become impossible !
URSI is one of the right fora to find measures for various radio uses to coexist.
References[1] Ehrenfreund et al. Rep. Prog. Phys., 65, 1427, 2002;[2] Dickens et al., ApJ, 479, 307, 1997;[3] Theule et al., 2011, A& A, 534, A64;
[4] Halfen et al., ApJ, 767, id.66- 28 -

Biographical Sketches of Speakers- 29 -

Keynote Lecture 1: “Method of Analytical Regularization in Computational Photonics”
Prof. Alexander I. Nosich, National Academy of Sciences of Ukraine, Ukraine
Alexander I. Nosich is Professor and Principal Scientist of the Institute ofRadio-Physics and Electronics of the National Academy of Sciences of Ukraine,
Ukraine. He was born in Kharkiv, Ukraine, in 1953. He earned his M.S., Ph.D. and
D.Sc. degrees, in Radio Physics, from the Kharkiv National University in 1975,
1979 and 1990, respectively. Since 1979, he has been with the Institute ofRadio-Physics and Electronics of the National Academy of Sciences of Ukraine in
Kharkiv, currently as Principal Scientist, Professor, and Head of the Laboratory ofMicro and Nano Optics. Since 1992, he has held many guest fellowships andprofessorships in the EU, Japan, Singapore, and Turkey. In 1990 he was one of theinitiators and then technical program committee chairman of the international
conference series on Mathematical Methods in Electromagnetic Theory (MMET)
held in Ukraine. In 1995, he organized the IEEE Antennas and Propagation Society East Ukraine Chapter, thefirst one in the former USSR. In 2001-2003, he represented Ukraine, Poland and the Baltic States in theEuropean Microwave Association. Since 2008, he represents Ukraine and Georgia in the European
Association on Antennas and Propagation. He was elected Fellow of IEEE in 2004, via the Antennas andPropagation Society, and Senior Member of the Optical Society of America in 2012. In 2015, he was awarded
the title of Doctor Honoris Causa of the University of Rennes 1, France. His research interests include themethod of analytical regularization, propagation and scattering of waves, open resonators, open waveguides,
optical antennas, linear modeling of lasers, and the history of microwaves.

Keynote Lecture 2: “Harmonization Based on Regulatory Science for Scientific andCommercial Radio Uses”
Prof. Ryuji Kohno, Yokohama National University, Japan / University of Oulu, Finland
Ryuji Kohno is IEEE and IEICE fellows. He received the Ph.D. degree inDept.EE from the Univ. of Tokyo in 1984. Since 1998 he has been a Professor in
YNU. During 1984-1985 he was a Visiting Scientist in Dept.EE, Univ. of Toronto.
Since 2007, he is also a Finnish Distinguished Professor (FiDiPro) in Univ. ofOulu, Finland. The meanwhile, he was also a director of SONY CSL/ATL during
1998-2002 and was a director of the UWB Tech. Inst. and a program coordinatorof Medical ICT Inst. of National Institute of Information and Communications
Technology (NICT) during 2002-2011．Currently he is the CEO of the Universityof Oulu Research Institute Japan – CWC-Nippon Inc. Ltd. since March 2012. He
was a principal leader of MEXT 21st century and Global COE programs during
2002-2007 and 2008-2013, respectively. Since 2003, he is currently a director ofboth Medical ICT Center and Kanagawa Medical Device Regulatory Science Center in YNU. He is anassociate member of the Science Council of Japan since 2006. Prof. Kohno was elected to be a BoG
member of the IEEE Information Theory Society three times on 2000, 2002, and 2006. He was awarded
IEICE Greatest Contribution Award and NTT DoCoMo Mobile Science Award in 1999 and 2002,
respectively.

- 30 -

Special Lecture 1: “International Union of Radio Science (URSI): Its Mission, Structure
and Activities”
Prof. Makoto Ando, URSI Vice-President; Tokyo Institute of Technology, Japan
Makoto Ando received the D.E. degrees in electrical engineering from Tokyo
Institute of Technology, Tokyo, Japan in 1979. From 1979 to 1983, he worked at
Yokosuka Electrical Communication Laboratory, NTT, and was engaged indevelopment of antennas for satellite communication. He moved to Tokyo
Institute of Technology in 1983 and is currently a Professor and the Executive
Vice President for research of the Tokyo Institute of Technology. His maininterests have been high frequency diffraction theory such as Physical Optics andGeometrical Theory of Diffraction. His research also covers the design ofwaveguide planar arrays and millimeter-wave antennas. He received theAchievement Award and the Paper Awards from IEICE Japan in 1993 and 2009.
He also received the 8th Inoue Prize for Science in 1992, the Meritorious Awardof the Minister of Internal Affairs and Communications and the Chairman of the Broad of ARIB in 2004
and the Award in Information Promotion Month 2006, the Minister of Internal Affairs and Communications.
He received IEICE Distinguished Achievement and Contributions Award in 2014.
He served as the guest editor-in-chief of more than six special issues in IEICE, Radio Science andIEEE AP. He was the general chair of the 2004 URSI EMT symposium in Pisa and of the ISAP 2007 in
Niigata. He was the Program Officer for engineering science group in Research Center for Science Systems,
JSPS in 2006-2009. He served as the member of Scientific Council for Antenna Centre of Excellence ACE in EU’s 6’th framework programme since 2004. He served as the Chair of Commission B of URSI
2002-2005. He was the 2007 President of Electronics Society IEICE, the 2009 President of IEEE Antennas
and Propagation Society and is currently the Co-chair of ISAP International Steering Committee andvice-president of URSI and IEICE. He is the Fellow IEEE and IEICE.

Special Lecture 2: “Temperature-Insensitive Quartz Oscillator Enabled Highly Stable
Communications and Clocks – Review of Issac Koga’s Works –”Prof. Emeritus Kenichi Iga, Professor Emeritus / Former President, Tokyo Institute ofTechnology, Japan
Kenichi Iga received the B.E., M.E., and Dr. Eng. degrees from the Tokyo
Institute of Technology, in 1963, 1965 and 1968, respectively. From 1968, hejoined the P&I Lab., Tokyo Institute of Technology, became Associate Professor
in 1974 and Professor in 1984. He retired Tokyo Institute of Technology inMarch 2001 and was awarded by Professor Emeritus.
Dr. Iga joined Japan Society for the Promotion of Science (JSPS) as
Executive Director from April 2001 to September 2007. He has been serving as
the President of Tokyo Institute of Technology from October 2007 till
September 2012. From 1979 to 1980, he stayed at Bell Laboratories as Visiting
Technical Staff Member.
Prof. Iga first proposed and pioneered the research of surface emitting laser(VCSEL) and microoptics. He received 1992 IEEE/LEOS William Streifer Award, 1998 IEEE/OSA John
Tyndall Award, 2003 IEEE Daniel E. Noble Award, 2002 Rank Prize, 2001 Purple Ribbon Prize, 2003
Fujiwara Award, 2007 C&C Prize, and 2009 NHK Broadcast Cultural Award. In 2013, he was awarded
2013 Franklin Medal with the Bower Award and Prize in Science.

- 31 -

Invited Paper 1 (Commission A):“Present Ion Optical Frequency Standards and New Challenges”
Dr. Nozomi Ohtsubo, National Institute of Information and Communications Technology,
Japan
Nozomi Ohtsubo received the B. S., M. S. degrees from Keio University, Japan,
in 2008 and 2010, the Ph.D. degree from University of Tokyo, Japan, in 2013.
He joined National Institute of Information and Communications Technology(NICT) as a post-doctoral researcher in 2013. He is engaged in the ion opticalfrequency standard.

Invited Paper 2 (Commission B):“Program of the Antarctic Syowa (PANSY) MST/IS Radar”
Prof. Toru Sato, Kyoto University, Japan
Toru Sato received his B.E., M.E., and Ph.D. degrees in electrical engineeringfrom Kyoto University, Kyoto, Japan in 1976, 1978, and 1982, respectively. He
has been with Kyoto University since 1983 and is currently a Professor in the
Department of Communications and Computer Engineering, Graduate School ofInformatics. His major research interests include system design and signal
processing aspects of UWB radars, atmospheric radars, radar remote sensing of
the atmosphere, and biomedical imaging. He is a fellow of the Institute ofElectronics, Information, and Communication Engineers of Japan, and a memberof the Society of Geomagnetism and Earth, Planetary and Space Sciences, theJapan Society for Aeronautical and Space Sciences, the Institute of Electrical andElectronics Engineers, and American Meteorological Society.

- 32 -

Invited Paper 3 (Commission F):“On Physical Limit of Wireless Data Transmission from Radiowave
Propagation Viewpoint”
Prof. Yoshio Karasawa, The University of Electro-Communications, Japan
Yoshio Karasawa received B.E. degree from Yamanashi University in 1973 andM.E. and Dr. Eng. Degrees from Kyoto University, Japan, in 1977 and 1992,
respectively. He joined KDD R&D Labs. in 1977. From July 1993 to July 1997,
he was a Department Head of ATR Optical and Radio Communications Res.
Labs. and ATR Adaptive Communications Res. Labs. both in Kyoto. Since 1999,
he has been a professor in the University of Electro-Communications (UEC),
Tokyo, and acted as the first director of Advanced Wireless Communication
research Center (AWCC) in UEC from 2005 to 2007. He has engaged in studies
on wave propagation and antennas, particularly on theoretical analysis andmeasurements for wave-propagation phenomena, such as multipath fading in
mobile radio systems, tropospheric and ionospheric scintillation, and rain
attenuation. His recent interests are in frontier regions bridging “wave propagation” and “digitaltransmission characteristics” in wideband mobile radio systems such as MIMO. In the field of wirelesscommunications, particularly of radiowave propagation, he authored/co-authored 10 text books, more than
200 journal papers, and 200 international conference papers.
Dr. Karasawa received the Young Engineer Award from IECE of Japan in 1983, the Meritorious Award
on Radio from the Association of Radio Industries and Businesses (ARIB, Japan) in 1998, Research Award
from ICF in 2006, two Paper Awards from IEICE in 2006, and Best Tutorial Paper Awards in 2007 and2008 from ComSoc of IEICE. He is a fellow of the IEEE and IEICE.

Invited Paper 4 (Commission K):“Recent Medical Applications of Electromagnetic Waves”
Prof. Koichi Ito, Recent Medical Applications of Electromagnetic Waves
Koichi Ito received the B.S. and M.S. degrees from Chiba University, Japan,
and the D.E. degree from Tokyo Institute of Technology, Japan. He is currently
a Professor at the Center for Frontier Medical Engineering, Chiba University.
From 2005 to 2009, he was Deputy Vice-President for Research, Chiba
University. From 2009 to 2015, he served as Director of the Center for Frontier
Medical Engineering, Chiba University.
His main research interests include small antennas for mobile
communications, microwave antennas for medical applications such as cancer
treatment, research on evaluation of the interaction between electromagnetic
fields and the human body by use of phantoms, and antenna systems forbody-centric wireless communications.
Professor Ito is a Fellow of the IEEE, a Fellow of the Institute of Electronics, Information andCommunication Engineers, Japan (IEICE). He served as Chair of the Technical Committee on Human
Phantoms for Electromagnetics, IEICE, from 1998 to 2006, Chair of the Technical Committee on Antennas
and Propagation, IEICE, from 2009 to 2011, Chair of the IEEE AP-S Japan Chapter from 2001 to 2002,
General Chair of the 2008 IEEE International Workshop on Antenna Technology (iWAT2008), an AdCom
member for the IEEE AP-S from 2007 to 2009, an Associate Editor for the IEEE Transactions on Antennas
and Propagation from 2004 to 2010, a Distinguished Lecturer for the IEEE AP-S from 2007 to 2011,
General Chair of the 2012 International Symposium on Antennas and Propagation (ISAP2012), and amember of the Board of Directors, BEMS, from 2010 to 2013. He currently serves as Chair of the IEEE
AP-S Committee on Man and Radiation (COMAR) and a Councilor to the Asian Society of Hyperthermic
Oncology (ASHO). He has been elected as a delegate to the European Association on Antennas andPropagation (EurAAP) since 2012 and Chair of Commission K, Japan National Committee of URSI
(International Union of Radio Science) since 2015.

- 33 -

Invited Paper 5 (Commission D):“Multi-port Power Router and Its Impact on Future Smart Grid”
Prof. Yuichi Kado, Kyoto Institute of Technology, Japan
Yuichi Kado received M.S. and Ph.D. degrees in electronics from Tohoku
University, Miyagi, Japan, in 1983 and 1998, respectively. In 1983 he joined theElectrical Communication Laboratories of Nippon Telegraph and Telephone
Public Corporation (now NTT), Kanagawa, Japan, where he was engaged inresearch on SOI structure formation by hetero-epitaxial growth. From 1989 to
1998 he worked on the development of fully depleted CMOS/SIMOX LSIs andultra-low-power CMOS circuits. From 1999 he was engaged in R&D on
compact network appliances using ultralow-power CMOS circuit technologies
for ubiquitous communications. He led research and development projects on
ultra-low-power
network
appliances,sub-terahertz-wave
wireless
communication, and intra-body communication as a director of Smart Devices
Laboratory at NTT Microsystem Integration Laboratories (2003-2010). In July 2010, he joined theDepartment of Electronics, Kyoto Institute of Technology, Kyoto, Japan. His current research interests
include multi-port power router to build energy interchanging systems and human body communication. He
has been the recipient of awards including the 2009 Nikkei BP Technology Award, the 2009 Radio wave
Achievement Award presented by the ARIB, and the 2009 Telecom System Technology Award presented by
the Telecommunications Advancement Foundation. He is a board member of NPERC-J(http://www.nperc-j.or.jp) and a member of IEEE.

Invited Paper 6 (Commission G):“Study of Coupling Processes in the Solar-Terrestrial System”
Prof. Toshitaka Tsuda, Kyoto University, Japan
Toshitaka Tsuda received B. Eng., M. Eng., and Dr. Eng. degrees in Electronics
from Kyoto University in 1975, 1977, and 1982, respectively. He was AssistantProfessor in Electronics, Kyoto University since 1977 and Associate Professor inRadio Atmospheric Science Center (RASC), Kyoto University since 1987. From
1995, he became Professor at RASC, Kyoto University. He is currently Director
and Professor (Laboratory of Atmospheric Sensing and Diagnosis) in ResearchInstitute for Sustainable Humanosphere (RISH), Kyoto University, and Vice
Executive Director, Kyoto University. His research area is observational studiesof the middle and upper atmosphere dynamics and meteorology withatmospheric radars and GPS meteorology techniques, and development of new
radio and optical remote-sensing techniques. He is the President of Japan
Geoscience Union (JpGU) since May, 2012, and a liaison-member of Science
Council of Japan since September, 2006. He served Co-Chair of Climate and Weather of the Sun-Earth
System (CAWSES) Ⅱ which is an international program sponsored by Scientific Committee on
Solar-Terrestrial Physics (SCOSTEP) of ICSU in 2009-2013. He was awarded Tanakadate Award(SGEPSS) in 1985, Horiuchi Award (Meteorological Society of Japan) in 1994, and Society Award(Meteorological Society of Japan) in 2003. He is also an ISI highly cited scientist in Geoscience.

- 34 -

Invited Paper 7 (Commission H):“Spacecraft-Plasma Interaction Effects on In-Space Electric Field
Measurements”
Assist. Prof. Yohei Miyake, Kobe University, Japan
Yohei Miyake received his B.E., M.E., and Dr. Eng. degrees in electricengineering from Kyoto University in 2004, 2006, and 2009, respectively. He
was a post-doctoral researcher of the Academic Center for Computing andMedia Studies, Kyoto University from 2009 to 2011. Since 2011, he has beenwith the Graduate School of System Informatics and the Education Center on
Computational Science and Engineering, Kobe University. In 2014, he became
an assistant professor in the Graduate School of System Informatics, Kobe
University. He has been engaged in the study of spacecraft-plasma interactions
and high-performance computing in plasma physics. He was awarded URSI
Young Scientist Award in 2005. He is a member of Society of Geomagnetismand Earth, Planetary and Space Sciences (SGEPSS), Japan Geoscience Union(JpGU), American Geophysical Union (AGU), and Japan Society for Computational Engineering andScience (JSCES).

Invited Paper 8 (Commission C):“R&D of Mobile Communication Systems Using Millimeter Wave”
Dr. Hiroyuki Tsuji, National Institute of Information and Communications Technology,
Japan
Hiroyuki Tsuji received B.E., M.E., and D.E. degrees in Electrical Engineeringfrom Keio University, Tokyo, Japan, in 1987, 1989, and 1992, respectively. In
1992, he joined the Communications Research Laboratory (CRL, now part of the
National Institute of Information and Communications Technology, or NICT).
From 1999 to 2000, he was a visiting research fellow of the University of
Minnesota. Since 2002, he has been working at Yokohama National University,
Yokohama, Japan as a visiting professor, which is an additional post. His
research interests are array antennas for wireless communications and signal
processing for communications and other systems. Dr. Tsuji received the Young
Engineer Award from IEICE (The Institute of Electronics, Information andCommunication Engineers) of Japan in 1995, the Electrical Science andEngineering Promotion Award from the Promotion Foundation for Electrical Science and Engineering in
2010, and ITU-AJ Encouragement Award from the ITU Association of Japan (ITU-AJ) in 2014,
respectively. He is a member of IEICE of Japan.

- 35 -

Invited Paper 9 (Commission E):“EM Information Leakage from Smart Devices”
Assoc. Prof. Yuichi Hayashi, Tohoku Gakuin University, Japan
Yuichi Hayashi received the B.E. degree in computer science and engineeringfrom Aizu University, Aizuwakamatsu, Japan, in 2003, and the M.S. and Ph.D.
degrees in information sciences from Tohoku University, Sendai, Japan, in 2005
and 2009, respectively. He is currently an Associate Professor at the Faculty ofEngineering, Tohoku Gakuin University.
His research interests include electromagnetic compatibility andinformation security. He is the Chair of Electromagnetic Information Leakage
Subcommittee in IEEE Electromagnetic Compatibility Technical Committee 5.
He received IEEE EMC International Symposium on Electromagnetic
Compatibility Best Symposium Paper Award in 2013, Symposium on
Cryptography and Information Security Innovation Paper Award in 2014 and2015, Research Institute of Electrical Communication (RIEC) Award at Tohoku University in 2014, andWorkshop on Cryptographic Hardware and Embedded Systems (CHES) Best Paper Award in 2014.

Invited Paper 10 (Commission J):“The Universe, Life and the Radio Sciences”
Assoc. Prof. Masatoshi Ohishi, National Astronomical Observatory of Japan, Japan
Masatoshi Ohishi, the Director of the Astronomy Data Center, the National
Astronomical Observatory of Japan (NAOJ), is an active radio astronomer with abroad scientific background.
He has deep knowledge and experience on the molecules in space(interstellar molecules), and he and his colleague discovered many new
interstellar molecules by using the 45m radio telescope of the Nobeyama Radio
Observatory, NAOJ. His current scientific interest is on exogenous delivery ofprebiotic material from interstellar space to planets, as an origin of life.
He has variety of international activities: he was the president ofCommission 5 (Astronomical Data) of the International Astronomical Union(IAU) between 2009 and 2012. He serves the vice president of an IAU
Commission, Astrobiology, between 2015 and 2018. Since 1993 he joined the International
Telecommunication Union (ITU) in order to protect the radio astronomy observations against man-madeinterference. He was the chairman of Working Party 7D (radio astronomy) of the Radiocommunications
Sector of the ITU, between 2000 and 2010.
He is a member of the IAU, the URSI and the Committee on Space Research (COSPAR).

- 36 -

Abstracts(Contributed Poster Papers)- 37 -

A1

Observation of the Magnetic-Insensitive Clock Transition in the2
S1/2 (F = 0) −2 D3/2 (F = 2) Transitions in Single 171 Yb+
Yasutaka Imai, Kazuhiko Sugiyama, and Masao Kitano
Graduate School of Electronic Science and Engineering
Kyoto University
Kyotodaigaku-katsura, Nishikyo-ku, Kyoto-shi, Kyoto 615-8510, Japan
Email: imai@giga.kuee.kyoto-u.ac.jp

Optical clocks referenced to long-lifetime transitions in isolated atoms or ions in optical region
have a potential for better uncertainty than that of the current frequency standard defined in
microwave region. Extremely low uncertainties in optical clocks enable us to conduct precise
tests in the field of fundamental physics such as search for a time variation of the fine structure
constant α [1].
We aim at realization of optical clocks based on single 171 Yb+ ions. 171 Yb+ has the mF =0 − mF ′ = 0 clock transitions, which are free from the 1st-order Zeeman shift. Also, 171 Yb+ has
relatively simple hyperfine structures. A large diﬀerence in the sensitivity to the variation of α
between the quadrupole and octupole clock transitions enables us to search for a time variationof α by using only 171 Yb+ ions [2, 3]. In this work, we conduct single-ion spectroscopy of the
2S2
171 Yb+ , and identify the m = 0 − m ′ = 0
F
1/2 (F = 0)− D3/2 (F = 2) clock transitions inFcomponent.
The 2 S1/2 (F = 0)−2 D3/2 (F = 2) transition has five Zeeman components corresponding to∆mF = 0, ±1, ±2. In order to identify the mF = 0 − mF ′ = 0 component, we probed all
Zeeman components by frequency sweep of the clock laser by 20-kHz intervals. The result isshown in Fig. 1. We identified the mF = 0 − mF ′ = 0 component. We conducted spectroscopy
with higher resolution and resolved motional sidebands. Then, we probed the carrier spectrumwith even higher resolution. We detected the spectra two times by changing the sweep directionof the clock-laser frequency, in order to measure and remove a linear frequency drift in the clock
laser. We obtained a spectrum of the full width at half maximum of 380 Hz as shown Fig. 2.
0.5

0 - -20 - -1

0-0

0-1

0-2

Excitation probability

Excitation Probability

0.5
0.4
0.3
0.2
0.1
0

-20

-10

0

Detuning/ MHz

10

0.4
0.3
0.2
0.1
0

20

γ= 380 Hz

-2

-1

0
Detuning/ kHz

1

2

Figure 2. Carrier spectrum of the mF =0 − mF ′ = 0 component. Frequency is
swept in 20-Hz intervals.

Figure 1. All zeeman components of the2S2
1/2 (F = 0)− D3/2 (F = 2) transition detected by using a single 171 Yb+ ion.

References
1. T. Rosenband et al., Science 319, 1808, (2008).
2. R. M. Godun et al., Phys. Rev. Lett. 113, 210801 (2014).
3. N. Huntemann et al., Phys. Rev. Lett. 113, 210802 (2014).

- 38 -B1

The Physics of Angular Momentum RadioBo Thidé(1) and Fabrizio Tamburini(2)(1) Swedish Institute of Space Physics, Uppsala, Sweden
Acreo Swedish ICT AB, Kista, Stockholm, Sweden(2) Twist Off S.R.L., Padua, Italy
Radio science and technology, be it fundamental physical or information-theoretical research, implementations such as radio astronomy and radar probing, or wireless communication and navigation applications, currently use predominantly techniques that involve only the Poynting (linear
momentum density) vector. I.e., the methods used are almost entirely based on the electromagnetic
linear (translational) momentum physical layer.
As a complement and/or supplement to this conventional approach, innovative techniques rooted
in the electromagnetic angular (rotational) momentum physical layer have been advocated, andpromising results from proof-of-concept laboratory and real-world angular momentum wirelessand fiber communication experiments were recently reported. A physical observable in its own
right, albeit much more sparingly used than other electromagnetic observables such as the energy
and linear momentum, the angular momentum exploits the rotational symmetry of the electromagnetic field and the rotational (spinning and orbiting) dynamics of the pertinent charge and current
densities.
We have shown that the volume integrated angular momentum density, i.e., the total angular momentum emitted by any source, always tends asymptotically to a constant when the distance fromthe source volume tends to infinity. This proves that a radiation arrow of time exists also for angular momentum, as it does for linear momentum (the volume integrated Poynting vector) andthat therefore angular momentum physics (torque action) offers an alternative or a supplement toconventional linear momentum physics (force action) as a means of transferring information wirelessly over very long distances. This proof opens possibilities for, among other things, a more
flexible utilization of the radio frequency spectrum and paves the way for new information transfer methodologies, technologies and techniques. We discuss implementation aspects and illustratethem by examples based on analytic and numerical analyses.
We report results from experiments demonstrating the feasibility of using electromagnetic angularmomentum (vortices) in real-world applications, and how vortical (twisted) beams can be shaped
and their divergence controlled. A scenario with angular momentum transducers of new types,
including optomechancial ones and those based on the interplay between charge, spin and orbital
angular degrees of freedom, heralded by recent advances in spintronics/orbitronics, condensedmatter skyrmion, and quantum ring physics and technology, is briefly delineated.

- 39 -

B2

Asymptotic-Numerical Analysis for Transient Scattered Field by a Coated Cylinder witha Thin Lossy Medium Excited by a UWB Pulse Wave
Keiji Goto, Oki Okawa, and Naoki Kishimoto
Department of Communications Engineering
National Defense Academy
1-10-20 Hashirimizu, Yokosuka-shi, Kanagawa 239-8686, Japan
Email: keigoto@nda.ac.jp

The studies on high-frequency scattering analysis by a smooth convex body have been an importantresearch subject in the area of the analysis of radiation patterns of antennas mounted near curved surfaces
such as an aircraft fuselage and the radar cross section of an airborne vehicle [1].
We have derived a frequency-domain (FD) uniform asymptotic solution (FD-UAS), which is useful forengineering applications, for a two-dimensional scattering problem by a coated conducting cylinder
covered with a thin lossy medium [2]. The FD-UAS is uniform in the sense that it remains valid within thetransition region adjacent to the shadow boundary, and it smoothly connects a geometric optical ray (GO)solution and a geometrical theory of diffraction (GTD) solution exterior to the transition region,
respectively. The FD-UAS is valid for a source point and an observation point either near the coating
surface or in the far-zone. The FD-UAS is represented by a combination of scattered field component
solutions, namely, the geometric optical ray (GO) solution composed of a direct GO (DGO) and a reflected
GO (RGO), the extended uniform GTD (UTD) solution made up of a DGO and a pseudo surface diffracted
ray (pseudo SD), the modified UTD solution representing SD series, and the lowest order SD solution.
In this paper, by extending the FD-UAS in [2], we derive a TD asymptotic-numerical solution(TD-ANS) for the transient scattered field when a modulated ultra-wideband (UWB) pulse wave [3] isincident on a coated cylinder with a thin lossy medium [4]. We assume that the thickness of the medium isthin as compared with the wavelength of a central angular frequency of a UWB pulse wave. The TD-ANS
is represented by a combination of transient scattered field component solutions, namely, the GO solution,
the extended UTD solution, the modified UTD solution, and the lowest order SD solution.
The validity and computation rate of the TD-ANS is confirmed by comparing with a reference solution
which is computed numerically by applying the FD exact solution and the FFT numerical code to atransient scattered field integral. We show that the TD-ANS is effective in understanding the transient
scattering phenomena of a coated cylinder, because it can extract and observe each pulse wave elementfrom the response waveform even when the pulse wave elements overlap mutually.
Acknowledgments
This work was supported in part by Japan Society for the Promotion of Science (JSPS) KAKENHI Grant
Number 15K06094.
References[1] H. T. Kim and N. Wang, IEEE Trans. Antennas Propag., vol. 37, no. 11, pp. 1463-1472, Nov. 1989.
[2] K. Goto et al., IEICE Electron. Express, vol. 12, no. 3, 20141110, pp. 1-10, Feb. 2015.
[3] Federal Communications Commission (FCC), First report and order: Revision of part 15 of thecommission’s rules regarding ultra-wideband transmission systems, FCC 02-08, Washington, D.C., 22
April 2002.
[4] K. Goto et al., Proc. of 2007 International Conference on Electromagnetics in Advanced Applications(ICEAA2015), Torino, Italy (to be appeared).

- 40 -

B3

Preliminary Study of Physical Optics for the Prediction of Scattering from RoughSurface in mmWave
Rieko Tsuji, Kentaro Saito and Jun-ichi Takada
Department of International Development Engineering
Tokyo Institute of Technology
2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
Email: {tsujirieko, saitoken, takada}@ap.ide.titech.ac.jp

As the demands for high speed wireless communications increase, the utilization of mmWave in themobile wireless networks have become a topic of great interest. In those higher frequency band, it is thoughtthat the wave scattering becomes dominant propagation mechanism since the wavelength is relatively same
size to the roughness of objects’ surface [1]. However, scattering phenomenon is not fully clarified. Physical
Optics (PO) is an effective method to obtain electromagnetic fields of scattering with approximating theinduced current focusing on a locality of the scattering area. In this paper, we present the received power
characteristics scattered from a random rough surface which is characterized by rms height h with Gaussian
distribution and correlation length l with Gaussian correlation [2], calculated by PO. We also discuss theappropriate integral area introducing Fresnel Zone, which are transformed as an ellipse depends on theincident angle.
Figure 1. shows simulation setting. As we consider mmWave and indoor environment, such as a wallwith random objects like pipes, h is set as 1λ (10 to 1 mm) and l is set as 5λ (50 to 1 mm).

Figure 1. Simulation Setting.

Figure 2. Received Power Comparison (θcomponent).

Figure 2. shows the result of the received power comparison in terms of the integral area. In this figure,the integral areas are represented by the number of fresnel zone and the index of S refers the number in thelegend. Black line represent the received power from a flat surface by integral area of S1. You can observethat the peak power of S1 is almost 20 dB less than S1-flat, and the power spread much wider than S1. Morespecific, rms angular spread of S1 is 12.09 degree while it is 3.85 degree in case of S1-flat. Switching ourattention to the integral area, the peak received power difference between S3 and S4 is 2.04 dB, while it is0.25 dB in case of S4 and S5. It seems that the difference is suddenly decreased after S4. For the future plan,we will simulate in the case of variable range of roughness and compare simulation and theory.
References
1. S. Salous, M. Nekovee, V. Degliesposti, and S. Hur, “Millimeter-wave Propagation Characterization andModelling Towards 5G Systems,” 9th Meeting COST IC 1004, Aalborg, Denmark, 2014.
2. C. Corbel, C. Bourlier, N. Pinel, and J. Chauveau, “Rough Surface RCS Measurements and Simulations
Using the Physical Optics Approximation,” vol. 61, no. 10, pp. 5155–5165, 2013.

- 41 -

B4

Energy Distribution of Air-hole Waveguides with Dielectric Cylinder Outside of the
Defect Layer
Ryosuke Ozaki and Tsuneki Yamasaki
Department of Electrical Engineering, College of Science and Technology,
Nihon University
1-8-14 Surugadai, Kanda, Chiyoda-ku, Tokyo, 101-8308, Japan
Email: ozaki@ele.cst.nihon-u.ac.jp, yamasaki@ele.cst.nihon-u.ac.jp
Light or electromagnetic waves within a certain wavelength band are confined a relatively small region,
cannot be propagated in a particular direction. A photonic crystal such as optical nano periodic structure is wellknown as optical device for controlling the propagation of light by interaction of both the wave nature of lightand periodicity. In particular, guiding problems of electromagnetic waves are very interesting in many areas ofphysics and engineering. In order to understand the properties of the periodic structure waveguide, it is necessary
to examine the photonic band gaps or stop-pass band regions. But, numerical results of propagation properties
are not indicated in detailed. And also, in the numerical analysis of photonic crystals waveguide, it is employed
the FDTD method[1], FEM method based on the variational method[2], and another numerical techniques[3].
In previous paper[4], we have analyzed the periodically dielectric waveguides with different radius composedof dielectric circular cylinders array with air-hole cylinder, and investigated the influence of complex
propagation constants at the first stop band region and distribution of energy density in the defect area and byusing the combination of improved Fourier series expansion method and multilayer method for TE mode.
In this paper, we analyze the energy distribution of air-hole waveguide with dielectric cylinder outside of thedefects layer as shown in Fig.1. Consequently, the aim of this paper is to obtain optimum dielectric waveguidestructure to concentrate the energy into the defect area for TE and TM mode. Numerical results indicate for thecomplex propagation constants at the first stop band region and the energy density for TE mode compared withprevious dielectric waveguide structures[4].
x0 3( m ) b

rn

S1 ( 0 , 0 )

c

0

z
d1 d

1

d2 d

S2

D
d

E y or H y
1

d2 d

0

d1 d

p

p

S3 ( 0 , 0 )tn

Figure 1. Structure and coordinate system of air-hole waveguides with dielectric cylinder outside ofthe defects layer
References
1. A. Adibi, Y. Xu, R. K. Lee, and A. Yariv, “Properties of the Slab Modes in Photonic Crystal Optical
Waveguide,” J. Light Tech., vol.18, no.11, pp. 1554-1564, 2000.
2. M. Koshiba, Y.Tsuji, and S. Saeki, “High-Performance Absorbing Boundary Conditions for Photonic Crystal
Waveguide Simulations,” IEEE Micro. Wireless Lett., vol.11, no.4, pp.152-154, 2001.
3. K. Watanabe and Y. Nakatake, “Floquet-Mode Analysis of Two-Dimensional Photonic Crystal Waveguides
Formed by Circular Cylinders Using Periodic Boundary Conditions,” IEICE Trans. Electron., vol.E93-C,
no.1, pp.24-31, 2010.
4. R. Ozaki and T. Yamasaki, “Distribution of Energy Flow by Dielectric Waveguides with Different
Radius Composed of Dielectric Circular Cylinder and Air Cylinder”, Proc. ICEAA, pp.97-100, 2014.

- 42 -

B5

Wiener-Hopf Analysis of the Plane Wave Diffraction by a Thin Material Strip
Takashi Nagasaka
Department of Electrical, Electronic, and Communication Engineering
Chuo University
1-13-27 Kasuga, Bunkyo-ku, Tokyo 112-8551, Japan
E-mail: a05.jne3@g.chuo-u.ac.jp
The analysis of the scattering by material strips is an important subject in electromagnetic theory andradar cross section (RCS) studies. Volakis (Volakis, J. L., Radio Sci., Vol. 23, pp. 450-462, 1988) analyzed
the diffraction problem involving a thin material strip using the dual integral equation approach and theextended spectral ray method together with approximate boundary conditions. Volakis has obtained a highfrequency asymptotic solution to the original strip problem by superposing the diffracted fields from the two
independent material half-planes and the doubly/triply diffracted fields between the edges of the two halfplanes. Therefore, his analysis may not be applicable unless the strip width is large compared with thewavelength. In this paper, we shall consider the same strip geometry as in the paper by Volakis, and analyze
the H-polarized plane wave diffraction using the Wiener-Hopf technique together with approximate boundary
conditions.
The geometry of the material strip is shown in Figure 1, where r and r denote the relative
permittivity and the relative permeability of the material, respectively, and i is the incident field of H
polarization. Introducing the Fourier transform of the scattered field and applying approximate boundary
conditions in the transform domain, the problem is formulated in terms of the simultaneous Wiener-Hopf
equations, which are solved exactly via the factorization and decomposition procedure. However, the solution
is formal in the sense that infinite branch-cut integrals with unknown integrands are contained. By using arigorous asymptotic method together with a new special function defined in this paper, we shall derive ahigh-frequency solution of the Wiener-Hopf equations, which is expressed in terms of an infinite asymptotic
series. Our solution is valid for the strip width greater than about the incident wavelength and accounts for
all the higher order multiple diffraction effects rigorously. The scattered field is evaluated asymptotically bytaking the inverse Fourier transform and applying the saddle point method. It is to be noted that our final
solution is uniformly valid for incidence and observation angles. Numerical examples of the RCS arepresented for various physical parameters and far field scattering characteristics of the strip are discussed indetail. Some comparisons with Volakis are also provided.
Acknowledgment
The author would like to thank Professor K. Kobayashi for helpful discussion in the preparation of thismanuscript.
xi ( H yi )

b/2

a

r , ry0ab / 2

Figure 1. Geometry of the problem.

- 43 -z

B6

Maxwell-Schrödinger Hybrid Simulation for Optically Controlling Quantum States:
A Two-Level System Manipulated by a Light Pulse Pair
T. Takeuchi*, S. Ohnuki*, and T. Sako** College of Science and Technology, Nihon University
1-8-14 Surugadai, Kanda, Chiyoda-ku, Tokyo 101-8308, Japan
E-mail: csts13001@g.nihon-u.ac.jp

Recent remarkable advancements in laser technology have enabled us to control quantum states ofatoms, molecules, and quantum dots by laser pulses. This innovative technology of controlling quantum
systems by external electromagnetic fields would be a promising seed for the next generation key
technology in various applications of electrical and nano engineering. Pioneering efforts in this direction
have been already made, aiming at, such as highly-efficient state-selective photochemical reactions,arbitrary q-bits operation in quantum computation, and so on.
From theoretical points of view dealing with such control problems boils down to the issue ofdesigning an optimal spatiotemporal profile of the applied laser pulse (called light control pulse hereafter).
Recent experimental studies have demonstrated that adroit modulation to the laser pulses through aliquid-crystal phase modulator driven by a genetic algorithm enables them to control product yields incertain photochemical reactions. At the same time a significant theoretical progress in the last decade has
added a new dimension to the design of light control pulses, namely, designing light control pulses basedon optimal control theory.
In the present study we have focused on a quantum two-level system of an electron confined in aquasi-one-dimensional nanostructure modelling quantum dots or nanowires, which manifests an ideal
system for exploring basic q-bit operations in quantum computation. We have studied switching propertiesof this system manipulated by laser pulses, that is, an arbitrary and complete transfer of probability
densities over the intended two quantum levels by using theoretically-designed laser pulses. Light-control
pulses have been generated by our recently proposed designing method that relies on highly-accurate
solutions of the coupled Maxwell-Schrödinger equations [T. Takeuchi, et al., IEEE J. Quantum Electron.,50, 2014; T. Takeuchi, et al., PIER (Invited paper), 148, 2014]. The Maxwell-Schrödinger hybrid
simulation has been employed here since, as we have demonstrated recently, the incident laser pulse in thevicinity of the target system can be disturbed significantly by the induced radiation from the excitedelectrons of the target system itself [T. Takeuchi, et al., Phys. Rev. A, 91, 2015]. This effect of the local
modification of the incident laser fields has been thought to be negligibly small and thus has been
neglected, to our best knowledge, in all previous theoretical studies.
As a demonstration we have examined control ability of the light control pulses with respect to theswitching operations designed by our proposed method solving the coupled Maxwell-Schrödinger
equations and by the conventional one [Y. Ohtsuki, et al., J. Chem. Phys., 109, 1998], that solves only thetime-dependent Schrödinger equation, thus neglecting the local modification of the incident laser pulses by
the electron excitation. A comparison between the results employing these two distinct pulses has clearlyshown that the light control pulses obtained by the conventional scheme could hardly control the systemstably while does our proposed pulses thanks to incorporating the effect of this local modification. The
present study suggests that our proposed scheme of designing light control pulses can be an indispensable
technology for light-driven ultrafast computation.

- 44 -

C1

Twin is Better Than Single: A Hybrid Localization Method Using RSS and Time of Flight
Nopphon Keerativoranan
National Electronics and Computer Technology Center*
112, Khlong Neung, Khlong Luang, Pathum Thani
12120, Thailand
Email: ke.nopphon@gmail.com

Chong Kwon Kim
Department of Computer Science and Engineering
Seoul National University
1 Gwanak-ro, Gwanak-gu, Seoul 151-742, South Korea
Email: ckim@snu.ac.kr

A positioning system allows user to track their position in real-time by using the knowledge of radio frequency signal(RF). Global Positioning System (GPS) is the most widely used outdoor localization system; however, it is infeasible inindoor environment due to the rich of obstacles and multipath effect. General ranging-based localization technique that ismainly used in GPS, turns out to be difficult for applying in indoor environment because it requires a very precise timing
and synchronization which is unrealistic in distributed systems. To circumvent the difficulty of time synchronization,
ToF(Time of Flight) based ranging is adopted; however its accuracy depends on sampling frequency of ADCs. Anotherfeasible Indoor Positioning System (IPS) method is RSS-based fingerprinting where location is estimated by matching
Received Signal Strength (RSS) profile with the database that stored RSS fingerprints at reference locations. Again, RSS
measurements are not robust and highly time-varying, causing a large error in some applications.
This work presented a new way of indoor localization by using Time-of-Flight (ToF) as well as RSS. In theory, ToF is
better than RSS in term of reliability and accuracy because their temporal and spatial dynamicity are less than those of RSS.
However, ToF based methods are applicable only when many LoS (Line of Sight) paths are available. Instead of using ToF
as ranging information we use ToF as fingerprinting information. Also we propose to combine ToF and RSS such that the
proposed methods can be applicable to indoor environments with many obstacles. We have experimented various
estimation method in order to find a better solution than basic fingerprinting algorithm where accuracy is limited by a gap
between fingerprinting points. We found that by using the knowledge of neighbor fingerprinting locations, it cansignificantly improve the accuracy. We evaluated the performance of the proposed method via simulation and experiments,
and the result has shown that our work outperformed previous methods.
References
1.

Gunther, and C. Hoene, “Measuring Round Trip Times to Determine the Distance between WLAN Nodes,” in Proc. Of
Networking 2005.
2. Z. Li, D. C. Dimitrova, D. H. Raluy, and T. Braun, “TDOA for Narrow-band Signal with Low Sampling Rate and Imperfect
Synchronication,” Wireless and Mobile Networking Conference, WMNC 2014.
3. D. j Geiger, “High Resolution Time Difference of Arrival Using Timestamps for Localization in 802.11b/g Wireless
Networks,”Wireless Communication and Networking Conference, WCNC 2010.
4. T. Ohlemueller, F. Winkler, and E. Grass, “Radio Localization in OFDM networks using the Round Trip Phase,” 7thworkshop on Positioning Navigation and Communication, WPNC 2010
5. M. Ciurana, F. Barcelo-Arroyo, and F.Izquierdo, “A ranging method with IEEE 802.11 data freames for indoor localization,”
WCNC 2007.
6. K. Liu, X. Liu, and X. Li, “Guoguo: Enabling Fine-grained Indoor Localization via Smartphone,” MobiSys 2013.
7. N. B. Priyantha, A, Chakraborty, and H. Balakrishanan, “The Cricket Location-Support System,” MOBICOM 2000.
8. P. Bahl, and V. N. Padmanabhan, “RADAR: And In-Building RF-Based User Location and Tracking System,” IEEE
INFOCOMM 2000.
9. J. Wang, and D. Katabi, “Dude, Where’s My Card? RFID Positioning That Works with Multipath and Non-Line of Sight,”
SIGCOMM 2013.
10. J. Xiong, and K. Jamieson, “ArrayTrack: A Fine-Grained Indoor Location System,” NSDI 2013.
11. K. Joshi, S. Hong, and S. Katti, “PinPoint: Localizing Interfering Radios,” NSDI 2013.
12. S. Sen, J. K. Lee, K. H. Kim, and P. Congdon, “Avoiding Multipath to Revive Inbuilding WiFi Localization,” Mobisys 2013.
13. Q. Pu, S. Gupta, S. Gollakota, and S. Patel, “Whole-Home Gesture Recognition Using Wireless Signals,” MobiCom 2013.
14. F. Adib, Z. Kabelac, D. Katabi, and R. C. Miller, “3D Tracking via Body Radio Reflections,” NSDI 2014.
15. G. Wang, Y. Zou, Z. Zhou, K. Wu, and L. M. Ni, “ We Can Hear You with Wi-Fi,” MobiCom 2014.
16. S. Goswami, “Indoor Location Technologies,” Springer, 2013.
17. B. R. Mahafza, “Radar Systems Analysis and Design Using MATLAB,” CRC Press, Third Edition, 2013*This work was done while the author was a graduate student at Department of Computer Science and Engineering, Seoul
National University, 1 Gwanak-ro, Gwanak-gu, Seoul 151-742, South Korea- 45 -

C2

Localization of Coherent Targets in a Multiple-Input Multiple-Output Radar
Employing Spatial Smoothing between Virtual Antennas
Extrapolated from Real Antennas
Ryo Nishikawa and Takehiko Kobayashi
Wireless Systems Laboratory
Tokyo Denki University
5 Senju-Asahi-cho, Adachi-ku, Tokyo 120-8551, Japan
Email: nishikawa@grace.c.dendai.ac.jp

A multiple-input multiple-output (MIMO) radar system utilizes multiple antennas for both transmission
and reception. The system transmits signals that are orthogonal to each other and combines the received
signals at a signal processing unit. It can localize multiple targets through joint estimation of the directionsof departure and directions of arrival of the targets. The estimation is commonly based on an
eigenspace-based algorithm, such as multiple signal classification (MUSIC), Capon, or estimation of signal
parameters via rotational invariance techniques (ESPRIT). The ‫ܯ‬௧ ×‫ܯ‬௥ MIMO system, consisting of ‫ܯ‬௧transmitting and ‫ܯ‬௥ receiving antennas, can localize ‫ܯ‬௧ ‫ܯ‬௥ – 1 targets, when the received signals fromthe targets are mutually incoherent, and thus the rank of the correlation matrix of the signals is full. Whenthe signals are coherent, however, the rank deteriorates and the localization cannot be achieved.
This problem has been mitigated by employing a spatial smoothing technique [1], which defines
subarrays from the transmitting and receiving arrays, then averages the correlation matrices among thesubarrays, and recovers the rank of the matrix. The number of the simultaneous localizable targets, however,decreases to ܰ௧ ܰ௥ – 1, where ܰ௧ and ܰ௥ are respectively the number of subarray elements fortransmission and reception.
This paper proposes an amended spatial smoothing technique. It extrapolates virtual elements outside thereal transmitting and receiving arrays by estimating the arriving phase at those elements, and then averages
the correlation matrices of the subarrays consisting of the real and virtual elements. Assume transmitting
signals possessing (ultra) wide bandwidth yielding a high time resolution. Let ߬௞,௠,௡ be the echo arrival
time from the m-th (1 ≦ ݉ ≦ ‫ܯ‬௧ ) transmitting element to the n-th (1 ≦ ݊ ≦ ‫ܯ‬௥ ) receiving element via
k-th (1≦k≦K) target, where K is the total number of the targets. The ߬௞,௠,௡ can be approximated by eithera plane on a sphere against (m, n) for each of k, when the targets are respectively located in the far or thenear distance with respect to the sizes of the transmitting and receiving arrays. The receiving phases at thevirtual elements are then extrapolated from the approximated arrival times. The number of localizable
targets can be maintained, even when the received signals are coherent, by spatially smoothing thesubarrays clipped out from the extended arrays. The performance was evaluated through simulation: a 3×3
MIMO radar system, unable to localize five coherent targets, was extended to 5×5 and 7 × 7 systems,
which can localize the coherent targets.
Reference[1] W. Zhang, W. Liu, J. Wang, and S. Wu, “DOA estimation of coherent target in MIMO radar,” ICASSP
2013, pp. 3929-3933, Vancouver, Canada, May 2013.

- 46 -

C3

Human Motion Classification Utilizing Radio Signal Strength in Wireless Body Area
NetworkSukhumarn Archasantisuk*, Takahiro Aoyagi*, Tero Uusitupa, Minseok Kim**, and Jun-ichi Takada***Graduate School of Decision Science and Technology*
Graduate School of Engineering***
Tokyo Institute of Technology
2-12-1 O-okayama, Meguro-ku, Tokyo 152-8552, Japan
Faculty of Engineering, Niigata University**
8050 Ikarashi 2-no-cho, Nishi-ku, Niigata-shi, 950-2181, Japan
Human motion classification has been studied intensively since information of human activities is very
useful for many kinds of applications such as activity monitoring, physical rehabilitation, abnormal activity
detection and so forth. In wireless body area network (WBAN), human activity-related information is not
only beneficial for application level, but it can also be used to design reliable communication protocol, forinstance, in MAC layer. Generally, the human motion classification system has been researched using
camera-embedded equipment or body-mounted sensors such as accelerometer or gyroscope. However,
although we notice that human motions highly affect the characteristic of wireless channel, only few
researches considered using the radio signal strength to classify human motions. Therefore, we aim toinvestigate whether human motions could be recognized by the pattern of received radio signal strength.
The WBAN radio signal strength in several human motions was obtained from the numerical simulation
using FDTD generated at 403.5 MHz [1] and 2.45 GHz [2]. The WBAN channel propagation from thesimulation results was used in this research to construct the human motion classifier. The human motion
classification consisted of two main stages: one was pre-processing stage; the other one was classification
stage. Pre-processing stage was responsible for preparing input vector for the classification stage.
Time-series of received signal levels were extracted from the simulation result. In current research,time-domain features were computed from the time-series of received signal levels. The time-domainfeatures were then normalized and inputted to the classifier. The classifier was trained by four different
algorithms, which were back propagation (BP), k-nearest neighbor (kNN), support vector machine (SVM),
and decision tree (DT). 3-class classifier was constructed to classify 3 human motions including running,
walking and weakly walking. 5-class classifier was constructed to classify 5 human motions including 3
motions mentioned earlier, sitting/standing and sleeping motion. The classifier was constructedindividually for each particular sensor link. Performance of the classification was measured by theaccuracy rate, which refers to the percentage of correct classifications and total classifications. BP and
SVM performed well in this classification task with the accuracy rate of 73.9–95.7 and 75.4–93.7 percent
respectively, while kNN and DT obtained slightly lower accuracy rates, which were 68.1–91.3 and 63.8–
89.4 percent respectively. Among six data links, Navel-Chest (Tx-Rx) link provided the most recognizable
information, which led to 90.1–93.7 percent of accuracy rate when classifier was SVM. These results
suggest that the human motion classification can be employed in WBAN utilizing only the continuous radio
signal strength.
References
1. T. Aoyagi, M. Kim, J. Takada, K. Hamaguchi, and R. Kohno, “Numerical simulations for wearable
BAN propagation channel during various human movements,” IEICE transactions on communications, Vol. 94, No. 9, pp. 2496–2500, 2011.
2. T. Uusitupa and T. Aoyagi, “Analysis of dynamic on-body communication channels for various
movements and polarization schemes at 2.45 GHz,” IEEE Transactions on Antennas andPropagation, Vol.61, No.12, pp.6168–6179, 2013.

- 47 -

C4

Localization of Illegal Radios Utilizing Phase-Difference as Location Fingerprints
A. Haniz*, G.K. Tran**, K. Sakaguchi**, J. Takada*,
D. Hayashi***, T. Yamaguchi*** and S. Arata***✄ Dept. of International Development Engineering, Tokyo Institute of Technology
S6-4, 2-12-1 O-okayama, Meguro-ku, Tokyo, 152-8552 Japan
E-mail: ❢azril,takada❣@ap.ide.titech.ac.jp✄✄ Dept. of Electrical & Electronic Engineering, Tokyo Institute of Technology
2-12-1 O-okayama, Meguro-ku, Tokyo, 152-8552 Japan
E-mail: ❢khanhtg,sakaguchi❣@mobile.ee.titech.ac.jp✄✄✄ Koden Electronics Co., Ltd., 2-13-24 Tamagawa, Oota-ku, Tokyo, 146-0095 Japan
E-mail: ❢d-hay,ymgct-t,arata-s❣@koden-electronics.co.jp

Illegal radios have the potential to cause harmful interference to licensed radio systems, and may be
a large threat to those used by emergency services. Statistics of illegal radios by the Ministry of Internal
Affairs and Communications (MIC) in Japan indicate that there is still a large number of illegal radio
appearances [1], and it is crucial to localize and tackle these illegal radios to ensure security especiallyin large events such as the upcoming 2020 Tokyo Olympics.
Conventional localization techniques utilizing angle of arrival (AOA) and received signal strength(RSS) suffer from low accuracy in dense urban environments due to a large number of multipaths andthe lack of line-of-sight (LOS) between the target and Rx sensors. Fingerprint-based techniques havebeen proposed to utilize the richness of information contained in the multipaths [2], but generally theyrequire knowledge of the target’s signal parameters when collecting fingerprints in the training phase.
To solve the abovementioned problem, [3] proposed a fingerprint-based localization technique using thecross-correlation of channel impulse responses (CIR) as location fingerprints, and these fingerprints areinterpolated in bandwidth, frequency and spatial domains. However, they only utilized the magnitude offingerprints while the phase information was discarded.
In this paper, the original framework proposed in [3] is employed, and a novel technique which utilizes phase difference between elements of an antenna array as location fingerprints is proposed. In orderto perform frequency domain interpolation of location fingerprints, an AOA estimation technique utilizing data from several frequency bands simultaneously is proposed. The location of the illegal radio is
estimated by performing pattern matching between the training fingerprints and the actual illegal radio’s
fingerprints based on the minimum squared error criterion. The accuracy of the proposed algorithm wasevaluated through ray-tracing simulations using a 3-D model of the urban area surrounding Shinjuku
station, Tokyo. Results show that the proposed algorithm is comparable with results in the previous
literature, and a hybrid algorithm is expected to achieve a greater improvement in localization accuracy.
References
1. Ministry of Internal Affairs and Communications Japan, Number of Unlicensed Radio Stations
Detected and Dealt With [Online] Available: http://www.tele.soumu.go.jp/e/adm/
monitoring/summary/ad_pro/index.htm
2. Y. Jin, W. Soh, and W. Wong, “Indoor localization with channel impulse response based fingerprint
and nonparametric regression,” IEEE Trans. Wirel. Commun., Vol.9, No.3, pp. 1120-1127, 2010.
3. A. Haniz, G.K. Tran, K. Sakaguchi, J. Takada, D. Hayashi, T. Yamaguchi, and S. Arata, “Localization of Illegal Radios Utilizing Cross-correlation of Channel Impulse Response with Interpolation in
Urban Scenarios,” to appear in Proc. of IEEE MILCOM, 2015.

- 48 -

D1

Resonant Tunneling Diode Receivers for 300-GHz-band Wireless Communications
Kousuke Nishio1, Sebastian Diebold1, Shunsuke Nakai1, Kazuisao Tsuruda1,2, Toshikazu Mukai2,
Jaeyoung Kim2, Masayuki Fujita1, Tadao Nagatsuma1
1Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
2 ROHM Co., Ltd. Sensor Business Strategy, Optical Device R&D Div, 21 Saiin Mizosaki, Ukyo, Kyoto 615-8585, Japan
Email: kousukenishio118@s.ee.es.osaka-u.ac.jp, fujita@ee.es.osaka-u.ac.jp

Resonant Tunneling Diodes (RTDs) are promising devices for terahertz (THz) wireless communication systems
not only for a transmitter as an oscillator but also for a receiver as a detector [1]. There have so far been several
reports on the increase of the data rate with respect to the modulation bandwidth [2] and antenna bandwidth [3]. As
for the receiver using RTDs, the highest data rate is 13 Gbps [4]. One of the bottlenecks which determine thereceiver bandwidth is a baseband (BB) circuit after the demodulation. In this paper, we propose a technique toenhance the bandwidth of the BB circuit and demonstrate an error-free (bit error rate < 10-11) transmission at a data
rate up to 17 Gbps using the RTD receiver at 300-GHz band.
We examined two types of BB circuits as shown in Figure 1, which are connected with the RTD chip [3, 4] bybonding wires. One BB circuit (type A) consists of a coplanar strip line (CPS), a grounded coplanar waveguide(GCPW) with via holes and a SMA connector as shown in Figure 1-(a). The CPS-GCPW transition causes animpedance mismatch, which limits the available 3dB-bandwidth to 8 GHz, as shown in Fig. 2-(a). The other one(type B) shown in Figure 1-(b) is only constructed from a tapered coplanar waveguide (CPW) and a SMA connector.
The tapered structure ensures the impedance matching. In addition, two ground conductors of the CPW areconnected at the end. The 3dB-bandwidths are 8 GHz and 11 GHz for type A and type B, respectively.
We performed a wireless transmission experiment using the RTD receiver with a type B circuit, together with aphotonics-based transmitter [5]. The carrier frequency and the transmission distance were set to be 297 GHz and 30
mm, respectively. Figure 2-(b) shows an eye diagram at a data rate of 17 Gbps, and an error-free transmission was
also confirmed by the bit error rate tester.

Fig. 1 Pictures and schematics of BB circuits. (a) Type A, (b) Type B

Fig. 2 (a) Frequency dependence of detected power.
(b) Eye diagram

This work was supported in part by the Strategic Information and Communications R&D Promotion Program (SCOPE), from the
Ministry of Internal Affairs and Communications, Japan and JSPS Postdoctoral Fellowship for Overseas Researchers.

References[1] T. Nagatsuma et al., Int. Conf. Telecommunications Remote Sensing, 41 (2014).
[2] Y. Ikeda et al., IEICE Electron. Express, 12, 20141161 (2015).
[3] K. Tsuruda et al., Global Symposium on Millimeter Wave, T1B-3 (2014).
[4] S. Diebold et al., IEICE Tech. Rep, 115, MWP2015-14, 17 (2015).
[5] T. Nagatsuma, XXXIth URSI General Assembly and Scientific Symposium, DFC01.3 (2014)- 49 -

D2

Study on Hybrid RAN Using RoF and RoR for Distributed Small Cell Configurations
Katsutoshi Tsukamoto1 and Kazuo Kumamoto2
1
Faculty of Information Science and Technology and 2 Faculty of Engineering,Osaka Institute of Technology, Japankatsutoshi.tsukamoto@oit.ac.jp
Wireless access traffic is rapidly increasing to access variety cloud services, and the access methods
diversify in various types of radio air interfaces, 3.5G, 3.9G, and 4G and future 5G, or commercial / private
WLAN. Moreover, as IoT will be spread, various wireless sensor networks will be involved in increasing
air interface diversification and the traffic with various types of QoS. This trend requires the more andmore efficient use of radio frequency, and the size reduction of radio cell offering heterogeneous wirelessaccess services. An effective solution will be small cell architecture, however, front-haul networks for alarge number of small cells become important in their cost, construction easiness, universality for various
types of air interfaces. Then, Hybrid RAN (radio access networks) using RoF (Radio over Fiber) and RoR(Radio on Radio) are effective to realize such front-haul networks.
RoF applied for WDM-PON optical access network has been widely studied [1] as a next generation
mobile backhaul because of its high capacity, low transmission loss and easy construction of a large
number of RBSs. As a broadband and universal front-haul network, we have proposed Radio over Fiber
Distributed Antenna System in WDM-PON [2]. When implementing DAS on the WDM-PON architecture
as shown in figure 1, RoF technologies enable the efficient and economical antenna distribution for several
different types of MIMO RF signals. Some evaluation results in the capacity achieved by the system will bepresented.
The transparency of the front-haul networks can be also realized by RoR (Radio on Radio) networks [3,4], which can provide a free space for heterogeneous wireless services in millimeter-wave radio. Figure 2
illustrates the configuration of the proposed cascaded RoF and Digital MMW RoR front-haul networks forsmall cells. RoR using millimeter-wave has a potential broadband to transport microwave signals keeping
their signal formats and a distributing ability to some small cells. To combat against lower transmission
quality of MMW RoR compared with RoF, we employ digital MMW RoR link cascaded with RoF to
expect FEC effect. In Fig. 2, point-to-point MMW RoR is shown, but RoR can be easily extended to P-MP
system. Some basic experimental results in link design and signal quality will be presented.

Fig 1: RoF distributed antenna system.

Fig 2: Cascaded RoF and digital MMW RoR front-haul.

Acknowledgments: This research was partially supported by JSPS KAKENHI Grant Number 25330122
and also SCOPE, Project No. 135007101(finished in 2013), provided by MIC. The authors are also grateful
thank to Prof. K. Iwatsuki, Tohoku Univ., Dr. S. Kuwano and the researchers in ASNL of NTT for theirfruitful discussions and co-operations.
References:
1. C. Choi, et.al., “Mobile WDM Backhaul Network Design for LTE-Advanced and Beyond Systems”, APMP
2012, paper FC-2, April 2012.
2. T. Iwakuni, et. al., ''Analysis of Wireless Channel Capacity in RoF-DAS over WDM-PON System'', IEICE
Transactions on Electronics, Vol.E96-C, No.2, pp.171-179, Feb. 2013.
3. Kanno A., et. al., “90-GHz radio-on-radio-over-fiber system for linearly located distributed antenna systems”,
Photonics Global Conference (PGC), Vol.1, No.1, pp.1-5, 2012.
4. K. Kumamoto, er. al., '' Development of Entrance Link using Millimeter-wave Digital Radio-on-Radio for
Small Cell Mobile Access'', Proc. of XXXI URSI GAAS, Beijing, China, DC01.1, Aug. 2014.

- 50 -

D3

Terahertz Wireless Transmission Enabled by Photonics
Using Binary Phase-shift Keying at 300 GHz
Y. Yasuda*, Y. Fujita*, S. Hisatake*, S. Kuwano**, J. Terada**, A. Otaka**, and T. Nagatsuma** Graduate school of Engineering Science, Osaka University
1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
E-mail: yuuyasuda101@s.ee.es.osaka-u.ac.jp** NTT Access Network Service Systems Laboratories, NTT Corporation,
1-1 Hikari-no-oka, Yokosuka, Kanagawa 239-0847, Japan

Recently, there has been a growing interest in the application of terahertz (THz) waves (0.1 THz~10
THz) to ultra-fast wireless communications [1]. This paper presents the first real-time wireless transmission
at 300 GHz using a binary phase-shift keying (BPSK) modulation. A coherent transmitter based on
photonics has enabled an error-free (bit error rate: BER < 10-11) transmission at a record data rate of 40
Gbit/s.
Figure 1 shows the schematic diagram of wireless coherent transmission based on BPSK modulation.
In the transmitter, we proposed a new approach to actively stabilizing the phase of carrier signals by using
the optical frequency comb (OFC) [2]. Two optical sub-carrier signals at different frequencies are extracted
from the OFC using an optical filter, and combined with an optical coupler (OC). The phase fluctuations of
optical sub-carriers in the optical fibers between the optical filter and the OC is stabilized by locking thephases of the two optical sub-carriers to those of the OFC signals by negative feedback for phase shifters(PS). The phase stabilized optical sub-carriers are input to the photodiode to generate the THz waves at 330
GHz which corresponds to the difference between two optical sub-carriers. In the receiver, thesub-harmonic mixer demodulated the THz waves by mixing it with the local oscillator (LO) signal, andoutput the recovered data signal.
Figure 2 shows BER characteristics at 40 and 45 Gbit/s and an eye diagram at 40 Gbit/s. The
error-free condition was confirmed at 40 Gbit/s, when the transmitter power was -21 dBm. When we
increased the data rate to 45 Gbit/s, the BER was 1.6×10-6. The main reason of the BER deterioration at 45
Gbit/s is the band-limitation of the electro-optic modulator (EOM, electrooptic bandwidth of 25 GHz).

Figure 1. Schematic diagram of experimental setup for a  Figure 2. BER characteristic at 40 and 45
Gbit/s and an Eye diagram at 40 Gbit/s.
wireless coherent transmission based on BPSK modulation.
References
1. J. Federici and L. Moeller, “Review of terahertz and subterahertz wireless communications,” J. Appl.
Phys., 107(11), 2010, p.111101.
2. Y. Yoshimizu, S. Hisatake, S. Kuwano, J. Terada, N. Yoshimoto, and T. Nagatsuma, “Wirelesstransmission using coherent terahertz wave with phase stabilization,” IEICE Electronics Express,
10(18), 2013, pp. 1–8.

- 51 -

E1

Electromagnetic Imaging and Monitoring of Active Volcanoes: a Case Study and Future
EM Monitoring Plan for Mt. Kusatsu-Shirane Volcano, Japan
Yasuo Ogawa
Volcanic Fluid Research center, Tokyo Institute of Technology
2-12-1 Ookayama, Meguro, Tokyo 152-8551, Japan
Email: oga@ksvo.titech.ac.jp
Mineo Kumazawa
Earth and Life Science InstituteTokyo Institute of Technology
2-12-1 Ookayama, Meguro, Tokyo 152-8551, Japan
Email: kumazawa@elsi.jp

Diagnostic understanding of active volcanoes is one of the targets of strong social demands, as indicated by
the recent phreatic eruption of Ontake volcano on 27 September 2014 and also by many other volcanic
activities in Japanese Islands. Kusatsu-Shirane Volcano, in particular, is the most important target for the
volcanologists at Tokyo Institute of Technology, and we have been maintaining the volcano observatory
since 1986.
In the last decade, we have had successfully imaged the edifice of the volcano in three dimensions with
200m horizontal resolutions to a depth of 2km by measuring impedances using natural electromagnetic
fields in the frequency range between 0.1Hz and 10kHz (Nurhasan 2006; Nurhasan et al 2006). We havefound that the peak of the volcano has a bell-shaped electrical conductor, consisting of clay minerals(smectite), which works as an impermeable cap to trap underlying vapor and fluids. The micro-seismic
hypocenters of the volcano are consistently located under the bell-shape conductor. The pressure source
location for the recent expansion of the volcanic edifice since March 2014, is also located within thebell-shaped conductor. We think that the increase in the resistivity under the clay cap and the breakdown of
the bell-shaped cap itself will lead to a phreatic eruption.
The continuous igaming of the cap structure by the electromagnetic method will thus be a key tosuccessfully monitor the potential phreatic eruption. We now propose a controlled source electromagnetic
imaging and monitoring system, consisting of multiple current loops sources and magnetic impedance(MI) sensor receiver arrays to be installed around the peak of Kusatsu-Shirane volcano.
Following the concept of ACROSS (Accurately Controlled Routinely Operated Signal System) (Kumazawa
and Takei 1994), transmitting signals consist of frequency combs and receiving the signal by a large
number of stacking will enable detections of small magnetic signals and their temporal changes originatingfrom the deep geothermal system prior to future phreatic eruptions.
References
1. Nurhasan,.PhD Thesis, Tokyo Institute of Technology
2. Nurhasan , Y. Ogawa , N. Ujihara, S.B. Tank , Y. Honkura, S. Onizawa,, T. Mori , and M. Makino,
Two electrical conductors beneath Kusatsu-Shirane volcano, Japan, imaged by audiomagnetotellurics
and their implications for hydrothermal system, Earth Planets Space, 58, 1053-1059, 2006.
3. Kumazawa M, Takei Y, Active method of monitoring underground structures by means of accurately
controlled rotary seismic source (ACROSS) 1. Purpose and principle. In: Abstracts of fall meeting of
the Seismological Society of Japan, p 158, 1994.

- 52 -

E2

Development of Upgrade BOLT System for Improvement of Location Accuracy
Yasuhiro Akiyama, Hiroshi Kikuchi, Ting Wu, Michael Stock, Yoshitaka Nakamura,Satoru Yoshida, Tomoo
Ushio
Graduate School of Engineering, Division of Electrical, Electronic and Information Engineering
Osaka University
2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan
Email: kikuchi.hiroshi@comf5.comm.eng.osaka-u.ac.jp
As is well known, the lightning discharge process produces strong impulsive radiation, witha strong broadband spectrum. Observation sensors for lightning discharges senseelectromagnetic waves, mainly in the ELF to UHF range, and especially in the LF and VHF
bands. VHF band sensor sensors can observe lightning discharge process in detail but its
observation coverage is limited in about 20 km. On the other hand, LF band sensor can observe
lightning at much great distances in about several hundred kilometers. The scale of athunderstorm is about several tens of kilometers maximally. Therefore, LF sensors are useful
to observe lightning throughout a thunderstorm’s life cycle.
Our research group has been developing the Broadband Observation network for Lightning
and Thunderstorm (BOLT), which locates radiation sources associated with lightning
discharges in three spatial dimensions[1]. BOLT consists of 11 LF band sensors which detect
lightning pulses wide frequency range from 5 kHz to 500 kHz. We have been operating BOLT
in Kansai area of Japan, locating both cloud-to-ground and intracloud discharges. Currently,
the BOLT system observes about 100 to 1000 lightning pulses per a flash, but we are striving
to improve both the detection efficiency and the location accuracy. Figure 1 shows an exampleof the lightning location results using BOLT system on 28 Jul., 2015. The white circle indicatesthe position of the sensor site. The plus sign is the position of the detected lightning discharges.
In order to improve the location accuracy of the BOLT system, preliminary investigation showthat the number of sources located, increases dramatically when only the highest portion of
the BLOT frequency band is used far location. So, our research group has proposed improving
a new “DDT” antenna sensor design to improve the high frequency sensitivity of the antenna.
The DDT antenna consists of a modified charge amplifier circuit. In this research, we will show
the observational results with the BOLT system and discuss the advantages of the proposed
DDT antenna.

Figure 1. Lightning location mapping by BOLT system.
References
1. Y. Takayanagi, M. Akita, Y. Nakamura, S. Yoshida, T. Morimoto, T. Ushio, Z.-I. Kawasaki, and K.
Yamamoto, “Development of VLF/LF bands interferometer and its initial observations,” IEEJ.Trans.FM.,
No.131, pp. 716-722, 2011

- 53 -

E3

Atmosphere and ionosphere connection as revealed by the network observationof Very Low Frequency radio signals
Sujay Pal1 and Yasuhide Hobara1,2
Department of Communication Engineering and Informatics, The University of Electro-Communications,
Chofu, Tokyo, Japan 182-8585.
2
Earth Environment Research Station, The University of Electro-Communications, Tokyo, Japan,
Email: myselfsujay@gmail.com, hobara@ee.uec.ac.jp

1

In this paper, we reported short term to long term connection between the atmosphere and lower
ionosphere using the Very Low Frequency (VLF) radio signals observed at several receiving stations within
Japan. These VLF receivers are part of the UEC’s VLF/LF observation network and continuously measure
the amplitude and phase of narrow band VLF signals from the distant transmitters around the globe.
Measurements of VLF signals allow us to track the changes occurred in the Upper-Mesosphere and Lower
Ionosphere (UMLI) region as the VLF waves reflected back from this region while propagating a long
distance within the earth-ionosphere waveguide. To investigate the long term relationship between thelower atmosphere and UMLI region, we compared the amplitude of VLF signals with the atmosphericparameters such as Total Column Ozone (TCO) density and stratospheric temperature at various heights forthe first time for three different latitudinal regions. We show that the VLF amplitude is strongly correlatedwith the TCO density, stratospheric temperatures for mid-latitude propagation paths throughout the years.
For high and low latitude regions, this correlation between the VLF amplitude and atmospheric parameters
is small and not significant. This study indicates the experimental observation of latitudinal dependence of
atmospheric influence on the upper mesosphere.
We also analyzed the VLF signals corresponding to a major Sudden Stratospheric Warming (SSW)
event occurred in 2009 to investigate the short term atmosphere and ionosphere connection. As a result, wefound the effect of the SSW event on the VLF signals in mid-latitude and high-latitude propagation paths.
We found an enhancement of the nighttime VLF amplitude during the SSW event for all propagation paths.
The strength of quasi 16-day planetary waves was also increased during the SSW event consistent withstratospheric temperature enhancement.

References:
1. Correia, E., Kaufmann, P., Raulin, J.P., Bertoni, F., and Gavilan, H.R., 2011, “Analysis of daytime
ionosphere behavior between 2004 and 2008 in Antarctica”, J. Atmos. Sol. Terr. Phys., 73, 2272–2278.
2. Lastovicka, J., 2006, “Forcing of the ionosphere by waves from below”, J. Atmos. Sol. Terr. Phys., 68,
479–497.
3. Pal, S., Chakrabarty, S., and Chakrabarti, S.K., 2015, “On the use of Very Low Frequency transmitter
data for remote sensing of atmospheric gravity and planetary waves”, Adv. Space Res., 55, 1190-1198.
4. Silber, I., Price, C., Rodger, C.J., and Haldoupis, C., 2013, “Links between mesopause temperatures andground-based VLF narrowband radio signals”, J. Geophys. Res., Vol. 118, 1–12.

- 54 -

E4

A statistical study of sub-ionospheric VLF signal anomaly due to geomagnetic storms
K. Tatsuta1, Y. Hobara1,2
1
Department of Communication Engineering and Informatics,
The University of Electro-Communications, Tokyo, Japan.
2
Earth Environment Research Station,
The University of Electro-Communications, Tokyo Japan.
1-5-1 Chofugaoka, Chofu city, Tokyo 182-8585, Japan
Email: t1431069@edu.cc.uec.ac.jp, hobara@ee.uec.ac.jp

We investigate quantitatively the effect of geomagnetic storms on the sub-ionospheric VLF/LF propagations
for different latitudes based on 2-years nighttime data from VLF/LF observation network operated by thegroup of the University of Electro-Communications. Three statistical parameters such as average signal
amplitude, variability of the signal amplitude, and nighttime fluctuation were calculated daily during two
years for 16~21 independent VLF/LF transmitter – receiver propagation paths consisting of 3 transmittersand 7 receiving stations. These propagation paths are suitable to simultaneously study high latitude, low-mid
latitude and mid-latitude D/E-region ionospheric properties. We found that these three statistical parametersindicate significant anomalies exceeding at least 2 times of their standard deviation from mean value during
the geomagnetic storm time period in the high-latitude paths with an occurrence rate of anomaly between 40
and 50 % presumably due to the auroral energetic electron precipitation. The mid-latitude and low-mid
latitude paths have smaller influence from the geomagnetic activity because of smaller occurrence rate ofanomalies even during the geomagnetically active time period (from 20 to 30 %). The anomalies except
geomagnetic storm periods may be caused by atmospheric and/or lithospheric origins. The obtained statisticaloccurrence rates of ionospheric anomalies for different latitudinal paths during geomagnetic storm and nonstorm time periods are basic and important information not only to identify the space weather effects toward
the lower ionosphere depending on the latitudes but also to separate various external physical causes of lower
ionospheric disturbances.

- 55 -

E5

Statistical Spatial Distributions of Lightning with Charge Moment Changes
Over Northern Area of Japan by ELF and LLS Observations
Junpei Yamashita1, Y. Hobara1, 2,
1
Department of Communication Engineering and Informatics,
The University of Electro-Communications, Tokyo, Japan.
2
Earth Environment Research Station,
The University of Electro-Communications, Tokyo, Japan.
1-5-1 Chofugaoka, Chofu city, Tokyo 182-8585, Japan
Email: y1431106@edu.cc.uec.ac.jp, hobara@ee.uec.ac.jp
Electrical characteristics of Cloud to Ground Flashes (CGFs) over northern area of mainland Japan were
investigated by focusing on the information of Charge Moment Changes (CMCs). Peak current (Ip), polarity
information and lightning geolocation were detected from the conventional lightning location system(LLS), while corresponding lightning CMCs were calculated by using transient electromagnetic radiationsin the Extremely Low Frequency (ELF) range so-called ELF transients observed by the group of the
University of Elecro-Communications (UEC) in Moshiri, Hokkaido, Japan. Statistical spatial distributions of
CGFs with their CMCs and peak currents were obtained. In this work, we analyzed the lightning data fromthree successive years from the summer of 2011 to the end of 2014 including several data gaps in the ELF
data due to the system failure and/or electrical power cut. Therefore practical number of days for the data
analysis is about 800. Remarkable difference between CMCs and Ip was found in the spatial distributions of
CGFs. Negative CGs with a large CMCs were inferior to those for +CGs and were mainly distributed overinland regions. This tendency represents characteristics of lightning activity during the summer thunderstorm
activity in Japan. On the contrary, the number of positives with a large CMC was small and seen over the Sea
of Japan and along the Hokuriku coast, where winter lightning with a large positive CMC are expected [1].
Negative CGFs with a large Ip were distributed more over the coastal areas and the Pacific and Sea of Japan,
but positives were observed mostly over the Sea of Japan. These statistical distributions provide the basic
information of not only about the lightning physics but also to prevent/mitigate potential damages to thepower grid systems and reusable power generating systems because CGFs with a large amount of charge
amount play a critical role on incidents to these systems.

References
1. Y.Hobara, T.Inoue, M.Hayakawa, K.Shiokawa, “Deducing Locations and Charge Moment Changes ofLightning Discharges by ELF Network Observations in Japan”, IEEJ Transactions on Power and Energy,
Vol.133, No.12, pp.994–1000 DOI: 10.1541/ieejpes.133.994,2014.

- 56 -

E6

Modeling and possible determination of D-region ionospheric perturbation during
annular solar eclipse on May 21, 2012, from VLF signal simulation andmulti-propagation path observation from UEC-VLF network: preliminary results
Tamal Basak, Yasuhide Hobara
Department of Communication Engineering and Informatics, Graduate School of Informatics andEngineering
The University of Electro-Commuications
1-5-1 Chofugaoka, Chofu, Tokyo 182-8585, Japan
Email: tamalbasak@gmail.com

A major part of the path of the annular solar eclipse of May 21, 2012 (magnitude 0.9439) was oversouthern Japan. The D-region ionospheric changes associated with that eclipse, led to several degree ofobservable perturbations of sub-ionospheric very low frequency (VLF) radio signal. The University ofElectro-Communications (UEC) operated VLF network is distributed over Japan and the eclipse associated perturbations were recorded in several receiving stations (Rx ) simultaneously for the VLF signals
coming from NWC/19.8kHz, JJI/22.2kHz, JJY/40.0kHz, NLK/24.8kHz and other VLF transmitters (Tx ).
The objective of this work is to locate and quantify the amount of lower ionospheric perturbation fromobserved VLF signal changes during eclipse. Inui and Hobara (2014) has already analyzed the same
event using 2D-Finite Difference Time Domain (FDTD) method and reproduced most of the observed
VLF signal deviations. For simulation, they used same set of signal propagation path parameters for
all paths. We include respective parameters like ground conductivity map, anisotropic ionosphere fordifferent propagation paths from well known Long Wave Propagation Capability (LWPC) code. In thiswork, we considered the 2-parameter exponential D-region ionospheric model and assumed a generalized
space-time dependent 2-dimensional elliptical Gaussian distribution type model perturbation for ionospheric parameters, such as, effective reflection height (h′ ) and sharpness factor (β). The depth (△hmax ,△ βmax ) and half widths (σθ,φ ) of that Gaussian distribution are being used as optimizing simulation
parameters, keeping consistency with real eclipse shadow pattern. In the vicinity of the eclipse zone, for
VLF amplitude and phase simulation over several signal propagation paths of UEC-VLF network, weused the LWPC code as a tool (Pal et al. 2012). Results have both similarities and dissimilarities withactual observed signal perturbations according to the limitations of the model. Based on it, a preliminarymodel for determination of spatial scale and degree of lower ionospheric perturbation has been proposed.
Hence, the initial simulation results are presented here (Basak and Hobara, 2015).
References
1. T. Basak, Y. Hobara, “Modeling and possible determination of D-region ionospheric perturbation
during annular solar eclipse on May 21, 2012, from VLF signal simulation and multi-propagation
path observation from UEC-VLF network”, (paper in preparation), 2015.
2. S. Pal, S. K. Chakrabarti, S. K. Mondal, “Modeling of sub-ionospheric VLF signal perturbations
associated with total solar eclipse, 2009 in Indian subcontinent,” Advances in Space Research, Vol.
50, pp. 196-204, July 2012.
3. D. Inui, Y. Hobara, “Spatio-temporal characteristics of sub-ionospheric perturbations associated withannular solar eclipse over Japan: Network observations and modeling”, URSI-GASS, 2014 XXXIthURSI, 16-24 Aug. 2014, DOI: 10.1109/URSIGASS.2014.6929555.

- 57 -

E7

Bulk Current Injection Simulation Using Integrated Circuit Immunity Macro Model
Yosuke Kondo#*, Shinichiro Ueyama#, Masato Izumichi#, and Osami Wada*# Semiconductor Circuit R&D Division, DENSO CORPORATION
1-1 Showa-cho, Kariya, Aichi 448-8661, Japan
Email: yousuke_kondou@denso.co.jp* Department of Electrical Engineering, Graduate School of Engineering, Kyoto University
Nishikyo-ku Kyoto-daigaku-katsura, Kyoto, 615-8510, Japan

RF voltage ratio

As automotive semiconductors and electronic products have increased, electromagnetic susceptibility
has become more important for integrated circuit (IC) reliability in terms of safety. For conducted immunity,
the bulk current injection (BCI) method is widely used for testing automotive components. The difficulty ofpredicting the results of BCI tests causes a long term development of automotive components. Hence, aprediction method for BCI test such as numerical simulation is highly desired. This paper provides a new
simulation method for BCI tests. Injected RF disturbance that reaches to an integrated circuit (IC) is
calculated by electromagnetic analysis with a high accuracy injection probe model [1]. In order to define
the failure criterion, the direct power injection (DPI) based IC immunity macro model for conducted
immunity (ICIM-CI [2]) is adopted.
The DUT is an automotive IC chip which has a pair of power supply terminals (Vcc, Gnd) and one
analog output terminal (Vout). DPI tests are performed in order to determine the voltage criteria of failure(Vout voltage offset larger than 50 mV). Figure 1 shows the BCI setup which is simulated byelectromagnetic analysis. The IC is mounted on the EUT. The RF disturbance voltage injected by BCI
method is calculated and compared with the voltage criteria of failure. Figure 2 shows the calculated andsimulated results. If the RF disturbance voltage at the terminal exceeds the voltage criteria of failure, the10
BCI test result is judged as failure in the simulation. The
Vcc
experimental failures at 120 MHz, 200 MHz, 270 MHz
Vout
1
and 300-350 MHz are well reproduced in the simulation.
The results indicated that the proposed method has the
0.1
possibility to predict BCI test results.
Electric load(a) Voltage ratio of RF disturbance and IC voltage criteria

Injection probe

Output deviation mV

LISN

0.01

EUT

Wires

QTEM port

200
100
0

-100
-200
0

LISN & Electric load

Injection probe

EUT

100

200
300
Frequency MHz

400(b) Measured DC Vout voltage deviation

Figure 2. Results of BCI Simulation and experiment.

Figure.1 BCI simulation model

References
1. Y. Kondo, M. Izumichi, K. Shimakura, O. Wada, “Modeling of Bulk Current Injection Setup for
Automotive Immunity Test Using Electromagnetic Analysis,” IEICE Trans. Commun., Vol. E98-B,
no.07, pp.1212-1219, July, 2015.
2. F. Lafon, F. D. Daran, M. Ramdani, R. Perdriau, M. Drissi, “Immunity modeling of integrated circuit:an industrial case,” IEICE Trans. Commun., Vol. E93-B, no.7, pp.1723-1730, July, 2010.

- 58 -

E8

Estimation of Packet Error Rate Considering Pulse Duration of Burst DisturbanceKazuhiro Takaya#*, Daisuke Tomita*, Kouki Umeda*, Tohlu Matsushima*, Takashi Hisakado*, and Osami Wada*#: NTT Network Technology Laboratories*: Department of Electrical Engineering, Kyoto University#: 3-9-11 Midori-cho, Musashino-shi, Tokyo 180-8585, Japan
Email: takaya.kazuhiro@lab.ntt.co.jp

As one of the evaluation indexes of electromagnetic disturbances, the amplitude probability distribution(APD) has been specified in CISPR 16-1-1. In some previous researches, the estimation methods of the
impact on communication quality from APD of disturbances were proposed on the assumption that each bit
error occurs independently 1, 2. However, these estimation methods cannot come into effect for continuous
burst disturbances because the assumption of independent bit error does not hold good. Moreover, the datais transferred not by bits but by a packet in IP-based communication networks.
Kth burst disturbance
We have studied the impacts on packet error rate (PER)
Ts
of a direct-sequence spread-spectrum (DSSS) system by X パケット
Td,k
パケット
M
パケット
M
M
M
・・・
burst disturbances, and we propose the estimation method
nTp
Tp
of PER considering the collision probability between a
Figure 1. Communication model.
packet and a burst disturbance with various pulse durations
To simplify the relationship between PER and pulse
duration, the model with constant transmission period Tp,
constant packet length M, and constant transmission time
Pulse
duration
Disturbanceof a packet Ts is considered as shown in Fig. 1. Assuming
generatorthat the duration in which the amplitude of K-th burst
Burst disturbance
disturbance is greater than the level X is Td,k, the PER PPER
DSSS
Directional
DSSS
Attenuatortransmitter
coupler
receivercan be estimated by Ts, Td,k, and measurement time Tall.
Disturbance level [dBm]

-70

-80

-90

0.44

0.442

0.444 0.446
Time [s]0.448DSSS signal
Transmitted packet

ACK

(1)

Figure 2. PER measurement system.
Meanwhile, the APD PAPD is the fraction of time that the
disturbance amplitude is greater than the level X.

Figure 2 indicates the measurement system to confirm theavailability of our estimation, and Figure 3 compares theestimated results by proposed method and from APD. The
estimated result by proposed method agrees with measured
one, and the difference in PER between our proposal and
APD is Ts just as shown in the difference between Eqs. (1)and (2).

Packet error rate

(2)

0.5
0.4
0.3

Measured
Estimated
APD
Ts

0.2
0.1
0
0

10

20

Pulse duration Td [ms]Figure 3.

Measured and Estimated PER
for various pulse durations.

References
1. K. Wiklundh, “Relation between the Amplitude Probability Distribution of an Interfering Signal and its
Impact on Digital Radio Receivers,” IEEE Trans. EMC, vol.48, no.3, pp.537-544, Aug. 2006.
2. K. Gotoh, T. Shinozuka, S. Ishigami, and Y. Matsumoto, "Electromagnetic Disturbance Evaluation
using Pulse Duration Distribution," EMC Europe 2012, Roma, P3-3-8, Italy, Sep. 2012.

- 59 -

0.45

E9

Wideband Characterization of Multilayer Wave Absorber Using FDTD Method
Kosorl Thourn† , Takahiro Aoyagi†† , and Jun-ichi Takada†† Department of International Development Engineering
Graduate School of Science and Engineering, Tokyo Institute of Technology
2-12-1 O-okayama, Megoru-ku, Tokyo, 152-8552, Japan†† Department of Human System Science
Graduate School of Decision Science and Technology, Tokyo Institute of Technology
2-12-1-W9-110 O-okayama, Megoru-ku, Tokyo, 152-8552, Japan†
E-mail: kosorl.t.aa@m.titech.ac.jp, †† aoyagi.t.aa@m.titech.ac.jp, and † takada@ide.titech.ac.jp
Most of absorbing materials are the frequency dispersive materials which can be effectively modeled by the Finite-Difference Time-Domain (FDTD) method. To facilitate the FDTD simulation for
the wideband response, the appropriate parametric model of material as the function of frequency isutilized. In previous work [1], the frequency characteristic of the absorbing material made of foamed
polystyrene was approximated by the multiple poles Debye model. To widen the absorption bandwidth,
the multilayer absorbers are frequently used. However, their performance depends on the combination of
thickness and electrical properties (permittivity and/or permeability) of the material in each layer. In thispaper, based on the optimum design parameters from [2], the frequency characteristic of multilayer absorbers (from one layer to eight layers) are investigated and modeled by 1D-FDTD method with various
orders of Debye model for material property approximation. Figure 1 depicts the calculation model for
multilayer wave absorber using 1D-FDTD technique. Each layer is characterized by thickness dn and
Debye parameters computed from carbon content Gn . A modulated Gaussian pulse is used as the pulse
excitation source with a center frequency 5.2 GHz and a bandwidth 4 GHz. This pulse is injected into
the simulation domain at distance λ0 = 0.75 m from the front layer of absorber.

Pulse

0.25
Tran. Line Theory [2]1D-FDTD/Debye 1st Order
1D-FDTD/Debye 2nd Order
1D-FDTD/Debye 3rd Order

0.2

G 1 G2P
GN E
C

Reflectivity

C
P
M
L

Observation Plan

2λ0 (λ0 : max. wavelength in simulation) x

0.15

0.1

0.05

d1

d2

dN

0
0

1

2

3

4

5

6

7

8

9

10

Frequency [GHz]

Figure 1. Calculation model for multilayer absorber using 1D-FDTD technique.

Figure 2. Reflectivity of eight-layer absorber.

The performance of absorbers is evaluated by reflection coefficient. To validate this coefficient, it
is compared with one computed by transmission line theory. Figure 2 shows the reflection coefficientcomputed for the eight-layer wave absorber using tran. line theory and 1D-FDTD method. As a result,
the reflection coefficients calculated by both methods showed good agreement, and the effectiveness of
the Debye approximation by this research, which can be used in wideband FDTD simulation, is shown.
For future works, 1D, 2D and 3D design of wave absorber in time domain will be taken into account.
References
1. K. Thourn, T. Aoyagi, and J. Takada, “Characterization of Absorbing Material Using Multiple Poles Debye
Model and Its Application for Wideband Simulation in FDTD Method”, IEICE Technical Report, Vol. 115, No.
131, EMCJ2015-41, pp. 61-66, July 2015.
2. Y. Shimizu, “A Design of Multi-Layer Wave Absorber by Point Frequency Matching Method”, IEICE Trans.
Vol. J62-B, No. 4, pp. 428-434, April 1979 (in Japanese).

- 60 -E10

A Novel Prediction for Very Low Frequency Transmitter Signal Amplitude
using NARX Neural Network
Hendy Santosa* and Yasuhide Hobara**Department of Communication Engineering and Informatics,
The University of Electro-Communications, 1-5-1 Chofu, Tokyo 182-8585, Japan
E-mail: hobara@ee.uec.ac.jp

In the recent three decade Very Low Frequency (VLF) waves (3-30 kHz) from powerful transmitters havebeen used to study and monitor the lower ionospheric conditions. The ionospheric perturbations are identified
in relation with e.g. thunderstorm activity [1], in association seismic activity [2] and others. However time
temporal dependence of VLF amplitude has a complicated and large daily variabilities in general due tocombinations of both effects from above and below (atmospheric and crustal activities). And quantitative
contributions from different external sources have not known well yet. Thus the modelling and prediction of
VLF wave amplitude is an important issues to study the lower ionospheric responses from various externalparameters and to also detect the anomalies from unknown parameters.
A NARX (Nonlinear Autoregressive with Exogenous Input) neural network is used as a novel method
for predicting daily nighttime averaged amplitude of Very Low Frequency (VLF) transmitter signals
indicating the ionospheric perturbation around the transmitter-receiver path. The NARX neural networks has
a good accuracy in predicting time series data [3] and thus are more suitable for dynamic modelling [4]. The
NARX model, which was built based on daily input variables of various physical parameters such as
stratosphere temperature, cosmic rays and total column ozone, possessed good accuracies during the modelbuilding. The aim of the study is to model and predict daily nighttime average amplitude of VLF wave
propagation from Hawaii (NPM) to Chofu (CHO) Tokyo, Japan path using the NARX neural network.
As the result, the constructed models are capable of performing accurate one step (1 day) ahead
predictions with the Pearson r is 0.92 and RMSE is 2.0 dB as shown in Figure 1. Furthermore, our results
indicate the feasibility and reliability of predicting lower ionospheric properties by NARX neural network
approach, and provide physical inside on the responses of lower ionosphere due to various external forcings.

Figure 1. One Step (1 day) Ahead predictions of NARX N.N. model daily nighttime average amplitude of
VLF waves NPM-CHO path over the time interval from January 2011 to December 2013 (Observed VLFblue; OSA-red).
References
1. U. S. Inan, S. A. Cummer, and R. A. Marshall, “A survey of ELF and VLF research on lightningionosphere interactions and causative discharges,” J. Geophys. Res., vol. 115, p. A00E36, Jun. 2010.
2. M. Hayakawa and Y. Hobara, “Current status of seismo-electromagnetics for short-term earthquake
prediction,” Geomatics, Nat. Hazards Risk, vol. 1, no. 2, pp. 115–155, 2010.
3. J. Connor, R. Martin, and L. Atlas, “Recurrent neural networks and robust time series prediction,” Neural
Networks, IEEE …, vol. 5, no. 2, pp. 240–254, 1994.
4. N. Bhat and T. J. McAvoy, “Use of neural nets for dynamic modeling and control of chemical processsystems,” Comput. Chem. Eng., vol. 14, no. 4–5, pp. 573–583, 1990.

- 61 -

F1

Detections of Electromagnetic Pulses During and Prior to Earthquakes
Minoru Tsutsui
Professor Emeritus, Kyoto Sangyo University
3-1-47, Nan-Ryo-Cho, Uji, Kyoto 611-0028, Japan
Email: tsutsui@cc.kyoto-su.ac.jp

I have been observing electromagnetic (EM) waves in the earth and above the ground using sensor
systems composed of tri-axial magnetic search coils and a vertical electric field antenna. One sensor system
was inserted into a bore hole of 100 m in depth, and another was installed on the ground. When a big
earthquake (M 6.3) occurred at 115 km distance from the EM observation site on April 13, 2013, the sensor
system on the ground first detected a magnetic pulse, and 13 sec later, another magnetic pulses was
detected in the borehole together with a seismic wave detected on the ground surface. This result suggeststhat the EM pulses were excited due to piezo-electric effect in the earth crust caused by vibration ofaseismic wave, and they were radiated out of the ground surface [1].
In the further analysis of waveforms of EM pulses simultaneously detected above and under the groundsurface when a seismic wave arrived at the EM observation site, the EM pulse detected above the groundshowed an elliptic polarization whereas that in the earth showed a linear polarization. This result has
proved that the EM pulse wave radiated from the deep earth, passing through a boundary (ground surface)of two media with different refractive indices. Based on this fact, we can regarded that the ellipticallypolarization of EM wave detected above the ground is the result of radiation of EM wave in the earth. Next,
I proposed a hypothesis on an excitation mechanism of EM pulses in the earth. Then I proved the
hypothesis by a laboratory simulation experiment on EM radiations form earth’s crust.
I found an outstanding EM pulse detected above the ground at about 7 hour prior to an earthquake(M3.9) occurred at a distance of 24 km from the EM observation site(see Fig. 1), which showed ellipticalpolarization andcontained electric fieldcomponent. The
waveform of its electriccomponent was similar
to that obtained by thelaboratory experiment.
This kind of EM
pulse could be regarded
as a precursory signal of
the earthquake.
Figure 1 Spectrograms of EM pulses detected during earthquake(at 21:02) and prior to the earthquake (at 14:08) on
Jan. 14, 2015.
References
1. Minoru Tsutsui, Behaviors of Electromagnetic Waves Directly Excited by Earthquakes, IEEE
Geoscience and Remote Sensing Letters, Vol. 11, No. 11, pp-1961-1965, 2014.

- 62 -F2

A Study on Rain Area Motion Along the Propagation Path and Ground Wind Velocity
Around the Earth Station in Ku-Band Satellite Communications Links
Yasuyuki Maekawa and Yoshiaki Shibagaki
Department of Telecommunications and Computer Networks
Osaka Electro-Communication University
18-8 Hatsucho, Neyagawa, Osaka 572-8530, Japan
Email: maekawa@maelab.osakac.ac.jp
Recently, satellite communications links use frequency bands of higher than 10 GHz such as Ku band
(14/12 GHZ), in addition to C band (6/4 GHz). In these higher frequency bands, however, the rain
attenuation of radio waves becomes one of the severe problems of satellite communications. Thus, detailed
studies on rain attenuation characteristics are necessary for reliable operations of satellite communication
links, using site diversity techniques and so on. In this study, the rain attenuation characteristics of
Ku-band satellite communications links are investigated from the broadcast satellite (BS) signal levels
simultaneously obtained at Neyagawa campus of Osaka Electro-Communication University (Neyagawa,
Osaka), as well as its Shijonawate campus (Shijonawate, Osaka) and high school (Moriguchi, Osaka)1.
The direction and speed of rain area motions are estimated from the time differences in the attenuation
detected at these three locations in light of the operations of site diversity techniques. These results arecompared with the ground wind velocities observed by nearby AMeDAS stations2 at three locations in
Hirakata, Osaka, and Toyonaka cities from 2008 to 2011.
Figure 1 shows (a) the directions and speeds of rain area motion obtained from BS signal levelmeasurements at the three locations at Neyagawa, Shijonawate, and Moriguchi, together with (b) the
ground wind velocities measured at Hirakata AMeDAS station that is about 6.7 km northeast from the
Neyagawa campus, respectively, in the form of circular diagrams. It is found from Fig.1 that the speeds ofthe rain area motions in the sky are, as a whole, several times larger than the wind speeds measured by
AMeDAS on the ground. However, a considerable relationship is found between the directions of the rain
area motions and ground winds, although the directions of the ground winds tend to be aligned in thenortheast-southwest direction along the area of Osaka plains. Table 1 summarizes cross-correlation
coefficients between the directions of the rain area motions and ground winds obtained at Hirakata, Osaka,
and Toyonaka AMeDAS stations, respectively. It is found that there are still correlations that are notnegligible between the directions of rain areal motions and ground winds. Thus, the AMeDAS data that is
easily available in Web site seems to be hopeful to determine the direction of satellite stations when the site
diversity techniques are conducted.
(a)

(b)

Table1. Cross-correlation coefficients between
directions of rain area motion and ground winds
Hirakata

Osaka

Toyonaka

0.52

0.51

0.52

Figure 1. (a) Rain area motion and (b) ground wind velocities.
References
1. Y. Maekawa, K. Sawai. Y. Shibagaki, T. Takami and T Hatsuda, “Site diversity effects on Ku-bandsatellite signal attenuation related to rain area motion,” Proc. ISAP2008, 3A07-3, Taipei, Taiwan,
Oct.2008.
2. http://www.jma.go.jp/jp/amedas/- 63 -

F3

Effects of Rain Area Motion on the Ku-Band Satellite
Communications Links in Tropical Wet Season
Akihiro Tama, Yoshiaki Shibagaki and Yasuyuki Maekawa
Department of Telecommunications and Computer Networks
Osaka Electro-Communication University
18-8 Hatsucho, Neyagawa, Osaka 572-8530, Japan
Email: mf14a007@oecu.jp
The effects of precipitating clouds on satellite communications links are significant as the frequency of
the used radio waves becomes higher than 10GHz, especially in heavy rain areas such as Southeast Asia.
Osaka Electro-Communication University conducted the rain attenuation measurements of Ku-band upand down-link radio waves (14/12/GHz) in the tropical region, using Superbird C that connects ResearchInstitute of Sustainable Humanosphere (RISH, Uji) and Equatorial Atmosphere Radar (EAR, West Sumatra,Indonesia) in collaboration with Kyoto University1,2. In this study, relationship between rain attenuation
characteristics of Ku-band satellite signals and rain cloud motion is investigated in wet season of Indonesia,
with aids of the X-band meteorological radar located at the EAR by the Shimane University group.
Figure 1 shows the rain attenuation of the satellite up-link (14GHz) signals measured at EAR from12:00 to 18:00 LT on January 7, 2006. It is found that the rain attenuation as large as 4 dB is detectedduring about 20 min from 14:30. In this period, the horizontal display of the X-band radar echo
simultaneously indicates that the core of the precipitating clouds passes over EAR at the speed of 5.73 m/s.
Considering the core size of the precipitating clouds, the duration time of the rainfall at EAR is estimated to
be about 15 min, being in good agreement with the duration of the rain attenuation as was shown in Fig.1.
Table 1 summarizes yearly averages of the speed of cloud motions and the duration time of theattenuation for the satellite communication links obtained during the wet season in the years of 2005 and2006. Thus, the average duration time is found to be longer due to much slower moving speeds than thosein Japan, which are typically 10-30 m/s. These characteristics seem to be carefully taken into account indesign of satellite communications links in the tropical region.
Table 1. Yearly averages of speed of cloud motion
and duration time of the attenuation.
Year
Speed [m/s]Duration [min]
2005
4.485
39.3
2006
5.578
62.3

Figure 1. Rain attenuation of up-link signals.
Acknowledgments

The authors thank the Shimane University group that supplied the X-band radar data.

References
1. Y. Maekawa, T. Fujiwara, Y. Shibagaki, T. Sato, M. Yamamoto, H.Hashiguchi, S. Fukao, “Effects oftropical rainfall to the Ku-band satellite communications links at the Equatorial Atmosphere Radar
Observatory,” J. Meteor. Soc. Japan, vol. 84A, pp.211-225, 2006.
2. Y. Maekawa, S. Miyamoto, K. Sawai. Y. Shibagaki, T. Sato, M. Yamamoto, H. Hashiguchi, and S.
Fukao, “Estimation of Rain Attenuation Characteristics of Satellite Communication Links using
X-band Meteorological Radars,” Proc. of ICCAS-SICE09. 2A16-4, Fukuoka, Aug. 2009.

- 64 -

F4

Effects of Rain Area and Rain Front Velocities on Rain Attenuation
Characteristics in Ku-Band Satellite Communications Links
Naoki Kubota, Yoshiaki Shibagaki and Yasuyuki Maekawa
Department of Telecommunications and Computer Networks
Osaka Electro-Communication University
18-8 Hatsucho, Neyagawa, Osaka 572-8530, Japan
Email: mf15a002@oecu.jp
The recent satellite communications tend to use frequency bands of higher than 10 GHz, in addition to
the traditional C band of 6/4 GHz. In these higher frequency bands, however, the rain attenuation of radio
waves becomes one of the severe problems of satellite communications. Thus, detailed observations andstatistical analyses of rain attenuation characteristics are necessary for reliable operations of satellite
communication links, using site diversity techniques and so on. In this study, Ku-band satellite signal
levels have been simultaneously observed at Osaka Electro-Communication University (OECU, Neyagawa,
Osaka), the Research Institute for Sustainable Humanosphere (RISH, Uji, Kyoto) and the Shigaraki MU
radar observatory (Koga, Shiga) of Kyoto University since September 20021. The direction and speed ofrain area motions are estimated from the time differences in the attenuation detected at these three locations
in light of the effects of site diversity techniques. These results are compared with those of rain front
motions inferred from the weather charts published by Japan Meteorological Agency in Web sites2.
Figure 1 shows the motion of rain fronts extracted from the weather charts issued at 3 and 6 LT on
March 30, 2007. During this period, the cold front is found to move southeastward at a speed of 21.5 m/s,
with the direction of 130.8°clockwise from the north over Osaka area, as indicated by an arrow. The speedof the cold front is estimated from reduction rate of the weather chart. On the other hand, Fig.2 shows the
motion of the rain area estimated from the time differences in rain attenuation occurrence among the threelocations of OECU, RISH, and MU during the same period. As shown by thin and thick dashed lines thatindicate the positions of the rain front inferred at RISH and MU, respectively, the rain area associated withthe front is found to move southeastward at a speed of 22.9 m/s, and the direction of the motion is 123.1°
clockwise from the north, being in good agreements with those estimated from the actual rain front motionin the weather chart in Fig.1. Thus, the speeds and directions of the motion of rain area that are importantfor site diversity operations are shown to be easily inferred from the weather charts published in Web sites,
using simple image processing methods.

Figure 1. Rain front motion in weather charts

Figure 2. Rain area motion observed by three locations.

References
1. Y. Maekawa, T. Nakatani, Y. Shibagaki, and T. Hatsuda, “A study on site diversity techniques related
to rain area motion using Ku-band satellite signals,” IEICE Trans. on Commun., Vol. B91-B, No. 6,
pp.1812-1818, Jun. 2008．
2. http://www.jma.go.jp/jp/g3/- 65 -

F5

Effects of Rain Area Motion on the Ku-Band Satellite
Communications Links in Tropical Dry Season
Keigo Takemoto, Yoshiaki Shibagaki and Yasuyuki Maekawa
Department of Telecommunications and Computer Networks
Osaka Electro-Communication University
18-8 Hatsucho, Neyagawa, Osaka 572-8530, Japan
Email: mf14a006@oecu.jp
The effects of precipitating clouds on satellite communications links are significant as the frequency of
the used radio waves becomes higher than 10GHz, especially in heavy rain areas such as Southeast Asia.
Osaka Electro-Communication University conducted the rain attenuation measurements of Ku-band upand down-link radio waves (14/12 GHz) in the tropical region, using Superbird C that connects ResearchInstitute of Sustainable Humanosphere (RISH, Uji) and Equatorial Atmosphere Radar (EAR, West Sumatra,Indonesia) in collaboration with Kyoto University1,2. In this study, relationship between rain attenuation
characteristics of Ku-band satellite signals and rain cloud motion is investigated in dry season of Indonesia,
using of the X-band meteorological radar located at the EAR by the Shimane University group.
Figure 1 shows the rain attenuation of the up-link (14GHz) satellite signals measured at EAR from14:00 to 21:00 LT on September 16, 2006. It is found that the rain attenuation as large as nearly 12 dB is
detected during about 15 min around 17:10. In this period, the horizontal display of the X-band radar echo
indicates that the core of the precipitating clouds passes over EAR at the speed of about 4 m/s. Because
the core size of the precipitating clouds is about 4 km, the duration time of the rainfall at EAR is estimatedto be about 17 min, being in good agreement with the duration of the rain attenuation in Fig.1.
Table 1 summarizes the average speed of cloud motions and the average duration time of theattenuation for the satellite communication links obtained in the case of the peak rainfall rates above andbelow 50 mm/h, respectively, during the years of 2005 and 2006. It should be noted that the duration timeis still much longer in the case of peak rain fall rates larger than 50 mm/h. These characteristics seem tobe carefully taken into account in design of satellite communications links in the tropical region.
Table 1. Yearly averages of speed of cloud motion
and duration time of the attenuation.
Peak rain fall
Speed [m/s]Duration [min]rate [mm/h]< 50 mm/h
5.29
61.9
4.26
70.7
≧50 mm/h

Figure 1. Rain attenuation of up-link signals.
Acknowledgments

The authors thank the Shimane University group that supplied the X-band radar data.

References
1. Y. Maekawa, T. Fujiwara, Y. Shibagaki, T. Sato, M. Yamamoto, H.Hashiguchi, S. Fukao, “Effects oftropical rainfall to the Ku-band satellite communications links at the Equatorial Atmosphere Radar
Observatory,” J. Meteor. Soc. Japan, vol. 84A, pp.211-225, 2006.
2. Y. Maekawa, S. Miyamoto, K. Sawai. Y. Shibagaki, T. Sato, M. Yamamoto, H. Hashiguchi, and S.
Fukao, “Estimation of Rain Attenuation Characteristics of Satellite Communication Links using
X-band Meteorological Radars,” Proc. of ICCAS-SICE09. 2A16-4, Fukuoka, Aug. 2009.

- 66 -

F6

Channel Sounding in 2GHz and 60GHz Band for Multi-band Wireless LAN System
Tianyang MIN , Kentaro SAITO , Jun-ichi TAKADA
Department of International Development of Engineering
Tokyo Institute of Technology
2-12-1 O-okayama, Meguro-ku, Tokyo 152-8550, Japan
Email: mty@ap.ide.titech.ac.jp

Recent years, with the increasing demand of higher data rate, the promising millimeter wave (60GHz)
band which provides a huge bandwidth (57-66GHz) is widely being investigated. Compared to low
frequency signals in 2.4GHz, the millimeter wave signals are more fragile, and the typical coverage
is very limited due to significantly higher propagation loss. Also, because of the smaller wavelength,60GHz directional links are far less capable to diffract around obstacles than 2.4GHz. Hence, 60GHz
signals can be easily blocked by an obstacle, such as human, and observe a high penetration loss. Becauseof the reasons discussed above, to combine the merits of lower frequency band such as the large coverageand the merits of higher frequency band such as the high transmission rate, the viability of multi-band
wireless LAN that are capable of operating at both 2.4GHz and 60GHz bands becomes very attractive.
The multi-band channel sounding system consists of two channel sounders. For 2.4GHz band measurement, a low-cost and self-developed channel sounder which was developed based on USRP (Universal Software Radio Platform)[1] was used. We implemented the 60GHz channel sounding measurement
by using a channel sounder whose details was presented in[2].
In context of his multiband wireless LAN, we implemented channel sounding measurements at both
2.4 and 60GHz with the same transceiver location in indoor environment. The comparison results show
propagation path characteristics of 2.4 and 60GHz indoor radio channels to discuss their differences interms of the diffuse scattering power and the delay spread to provide insights on a possible channel
modeling strategy at the two frequencies.
By assuming that the channel impulse response is obtained under the wide-sense stationary condition
where the path loss and shadowing remains constant, we rotate both the Tx and Rx antennas over 360 ◦azimuth angle and acquire the channel transfer function at every azimuth angle. We derive power angularprofile by the squared average of channel impulse responses obtained with the Tx and Rx antennas
pointing at different azimuth angles. Figure shows the power angular profile of 2GHz band when Tx ispointing at 0 ◦ and Rx is pointing over 360 ◦ azimuth angle.

References
1. T. Min, M. Kim, J. Takada ,“Development of Channel Sounder using GNU Radio/USRP: Timing
Synchronization and Phase Alignment ,” IEICE technical report, vol. 114, no. 491, SR2014-131, pp.
121-126, Mar. 2015.
2. K. Wangchuk, P.K. Hung, M. Kim, J. Takada ,“mmWave Channel Sounding for Future Heterogeneous
Cellular Network ,” IEICE technical report, vol. 114, no. 492, SRW2014-52, pp. 31-36, Mar. 2015.

- 67 -

F7

Application of Adaptive Digital Beam Forming for Polarimetric Phased Array Weather
Radar
Hiroshi Kikuchi, Takuro Tashima, Ting Wu, Gwan Kim, Tomoo Ushio, Hideto Goto, Fumihiko Mizutani
Graduate School of Engineering, Division of Electrical, Electronic and Information Engineering
Osaka University
2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan
Email: kikuchi.hiroshi@comf5.comm.eng.osaka-u.ac.jp

In 2012, we have developed a rapid scanning phased array radar and installed at the campus in the Osaka
University. The phased array radar consists of the 128 antenna elements in tandem arrangements and achieves
the 3-D precipitation measurements in less than 10 or 30 seconds in a range of 20 or 60 km, respectively. The
phased array radar is also achieved a high density observation in elevation angles, which consists of about
100 sharped beams. Within a few years, we plan to develop a polarimetric phased array weather radar as anext generation rapid scanning weather radar, which has a dual-polarized antenna with two-dimensional
circular planar phase-array elements. The polarimetric measurements has the following advantages. First, it
is possible to discriminate the types of hydrometeors, which is an important step prior to rainfall estimates.
Seconds, an estimation of rainfall rates not only on the ground is expected to be available, which means thata vertical structure gives insight in precipitation processes. The under development radar is capable ofpolarimetric measurements for the 3-D rainfall distribution in less than 10 or 30 seconds in a range of 20 or
60 km, respectively. For rapid scanning, a digital beam forming (DBF) method is used as an electronic
scanning system. It is one of the main components to determine an observation accuracy of the phased arrayradar. This paper is focused on a DBF method for the under development radar. DBF is a signal processing
method that makes a directional beam pattern. We proposed to apply minimum mean square error (MMSE)
method. In [1], MMSE method is proposed as one of the adaptive digital beamforming method. The method
is useful for the estimation of the distribution targets as a precipitation. The under development planar phasedarray antenna is based on 6992 array elements that are placed in a circular shape. We conducted a numericalsimulation for precipitation measurements to evaluate an impact of a beamforming method. We comparedestimation results with MMSE to those with the conventional DBF method such as Fourier beam forming(FR), to evaluate the effect on DBF. From the results, MMSE’s performance was superior to DBF whencomparing all polarimetric variables.
In order to evaluate the effect on the MMSE’s performance, we applied the observational data using
the phased array radar installed at the Osaka University. From comparison with the conventional method,
the mean error is improved in the estimation results.
References
1. E. Yoshikawa, T. Ushio, Z. Kawasaki, S. Yoshida, T. Morimoto, F. Mizutani, and M. Wada, “MMSE
Beam Forming on Fast-Scanning Phased Array Weather Radar,” IEEE Transactions on Geoscience andRemote Sensing , vol. 51, no. 5, pp. 3077-3088, May 2013.

- 68 -

F8

Large Scale Subsurface Velocity Estimation with Array GPR System YAKUMO
Li Yi*, Kazunori Takahashi** and Motoyuki Sato*** Graduate School of Environmental Studies, Tohoku University,
41 Kawauchi, Aoba-Ku, Sendai 980-8576, Japan
E-mail: yi_li@cneas.tohoku.ac.jp** Center of Northeast Asian Studies, Tohoku University
41 Kawauchi, Aoba-Ku, Sendai 980-8576, Japan
E-mail: kazunori.takahashi@cneas.tohoku.ac.jp
sato@cneas.tohoku.ac.jp

We demonstrate an approach to estimate velocity of electromagnetic wave propagation in subsurface
with an array ground penetrating radar (GPR) system “YAKUMO” [1]. The common midpoint (CMP) dataacquired by YAKUMO at each position can be used to estimate the velocity at different depth at that
position [2]. As it shown in Figure 1(a), the CMP data acquired by the YAKUMO system only include few
traces and it may introduce some artifacts to the velocity spectrum which is used to estimate the velocity, inorder to estimate the velocity automatically, further processing is necessary.
In this paper, we introduced a transform based method to remove the aliasing in the CMP data, so thatthe irregular CMP data can be interpolated properly even there are only 8 traces are available in a CMP data.
The results are shown in Figure 1(b). After the proposed processing the quality of the velocity spectrum canbe enhanced. Hence the velocity profile can be generated by picking the velocity at each position
automatically.

(a)(b)
Figure 1. The simulated CMP data and its velocity spectrum; (a) Raw data; (b) Reconstructed data with theproposed processing.
We also applied this approach to the data obtained at a simulated airport runway. The results are shownin Figure 2. Beside man-made voids can easily be detected, we also found that the area with velocity about
0.03 m/ns lower than surrounding areas indicates the place where man-made grooves exist. The velocitychanges may be used to detect the damaged pavements.
References
1. M. Sato, "Array GPR “Yakumo” and its application to archaeological survey and environmental
studies," in Proc. General Assembly and Scientific Symposium (URSI GASS) 2014, Beijing , 2014, pp.
1-2.
2. Yi, K. Takahashi, and M. Sato, “A fast iterative interpolation method in f-k domain for 3D irregularly
sampled GPR data,” submitted to IEEE Journal of Selected Topics in Applied Earth Observations andRemote Sensing.
3. O. Yilmaz, “Velocity Analysis and Statics Corrections,” in Seismic Data Processing, Tulsa, Society of
Exploration Geophysicists, 1987.

- 69 -

F9

Enhanced Clustering of Multipath Components Utilizing Geometrical Parameters for
Dynamic Indoor Double-Directional Propagation Channel
P. Hanpinitsak*, K. Saito*, J. Takada* and M. Kim***Department of International Development Engineering, Tokyo Institute of Technology
2-12-1, O-okayama,Meguro-ku, Tokyo 152-8552, Japan
E-mail: hanpinitsak.p@ap.ide.titech.ac.jp, saitouken,takada@ide.titech.ac.jp**Department of Electrical and Electronic Engineering, Niigata University
8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-2181, Japan
E-mail: mskim@eng.niigata-u.ac.jp

-10
0

Tx
Rx

-20
-10
0

Tx
Rx

y-coordinate [m]-20y-coordinate [m]y-coordinate [m]This paper extends previously proposed clustering approach [1] which utilizes scattering points obtained from measurement-based ray tracer (MBRT) from static into dynamic double-directional propagation channel. Similar to the previously proposed approach, multipath components (MPCs) were extracted
from the measurement data obtained from wideband MIMO channel sounder at 11 GHz [2] which is a
future frequency band. After that, KPowerMeans algorithm (KPM) was used to group MPCs where it
was extended to dynamic clustering. In this paper, two major contributions was made: 1. MBRT withangular correction was proposed to overcome the limitation of array size in channel sounding in highfrequency bands by taking the physical environment into consideration and 2. KPM algorithm was extended to capture the dynamic behaviour of cluster.
Figure 1 shows the clustering results of three consecutive snapshots in an indoor environment where
the Tx is moving and Rx is static. It can be observed that the algorithm can capture dynamic behaviourof clusters accurately.

-20
-10
0

Tx
Rx

30 20 10 0

30 20 10 0

30 20 10 0

x-coordinate [m]x-coordinate [m]x-coordinate [m]Figure 1. Clustering results of three consecutive snapshots in an indoor environment.

References
1. P. Hanpinitsak, M. Kim, and J. Takada, “Clustering of multipath components utilizing geometricalparameters for indoor double-directional propagation channel,” in 2014 Asia-Pacific Microwave Conference (APMC), 2014, p. 1057.
2. Y. Konishi, Y. Chang, M. Kim, and J. Takada, “Versatile Radio Channel Sounder for Double Directional and Multi-link MIMO Channel Measurements at 11 GHz, ” IEICE Trans. Electron., E97-C,
10, 2014, pp. 994-1004.

- 70 -

F10

Fundamental Study of GB-SAR Imaging by Compressive Sensing
Lilong Zou, Motoyuki Sato
Graduated School of Environmental Studies
Tohoku University
41 Kawauchi, Aoba-ku, Sendai 980-8576, Japan
E-mail: zou_ll@cneas.tohoku.ac.jp

In recent years, Ground-Based Synthetic Aperture Radar (GB-SAR) has been widely applied formonitoring of landslides and infrastructures thanks to its capability of acquiring short ranges scene in justfew minutes. The transmitted and received antennas are moving along a rail to synthetize a 2D image and aﬁxed sampling of the system bandwidth and of the azimuth dimension is performed to obtain the GB-SAR
images that shows electromagnetic properties of the object. However, depending on the applications, themonitoring requirements and conditions and the system storage capabilities, a less than regular sampling ofreturned echos could be employed. One solution for this data acquisition problem is compressive sensing(CS). It is a recently introduced theory that permits the reconstruction of sparse signals even if a number ofsamples much lower than those required by the Shannon-Nyquist theorem is available.
In recent years, compressive sensing has been applied in the field of Ground Based Synthetic Aperture
Radar (GB-SAR) imaging and shows great potential. The existing models are, however, based on
application of the sensing matrix acquired by the exact observation functions. As a result, the corresponding
reconstruction algorithms are much more time consuming than traditional focusing methods, especially inhigh resolution systems. In this paper, we formulate a new GB-SAR imaging model based on the use of theapproximated GB-SAR observation deducted from the inverse of focusing procedures called Far-Field
Pseudopolar Format Algorithm (FPFA) [1]. The proposed model forms a new compressive sensing
GB-SAR imaging method that can be applied to high-quality and high-resolution imaging under
sub-Nyquist rate sampling, while saving the computational cost substantially in time. In order to do this, itis necessary to use a reconstruction algorithm: spectral projected gradient (SPGL1) algorithm [2] has been
employed because of its speed, transparency and guarantee of sparse recovery if the measurement matrix
satisﬁes the restricted isometry condition. Real GB-SAR data application supports that the proposed
method can perform GB-SAR imaging effectively and efficiently under Nyquist rate, especially for large
scale applications.

GB-SAR System(a)(b)
Fig.2. (a) Scene Map of North-Kawauchi Campus, Tohoku University; (b) Reconstructed image by
proposed method with 25% random samples data.
References
1. J. Fortuny-Guasch, “A fast and accurate far – field pseudo - polar format radar imaging algorithm,”
IEEE Trans. Geosci. Remote Sens., Vol. 47, No. 4, pp. 1187-1196, 2009.
2. E. van den Berg and M. P. Friedlander, “Sparse optimization with least-squares constraints,” SIAM J.
on Optimization, SIAM J. on Optimization, Vol.21, No.4, pp. 1201–1229, 2011.

- 71 -

F11

Leaf Shape Effects on the Scattering from Simplified Deciduous Leaf
Structures Using Spherical Wave Harmonics
P.J. Co and J. Takada
Department of International Development Engineering, Graduate School of Science andEngineering, Tokyo Institute of Technology
S6-4, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
E-mail: pjco@ap.ide.titech.ac.jp
E-mail: takada@ide.titech.ac.jp

Abstract— Leaf scatterers are a type of obstacle that is abundant in rural propagation environments. These leaf scatterers interact with propagating radio waves which may affect the performance of wireless systems operating in such environments. These leaf scatterers have complex
geometries, which provide a significant challenge in modelling their electromagnetic interactions with propagating radiowaves. One common approach is to simplify the leaf structures into
simpler geometries such as thin dielectric disks, cylinders and ellipsoids. In this paper, deciduous leaves are approximated as thin dielectric disks as well as thin dielectrics with realistic leaf
shapes. The realistic leaf shape is approximated using a polygonal model which can estimate
the leaf shape of real world leaves. The electromagnetic scattering is modelled as a T-Matrix
evaluated from the least squares solution of numerical results, which describes the scattering interms of the spherical wave harmonics. The scattered fields from similarly sized leaf shapes areevaluated from the T-matrix model and is compared with the simple disk shape approximation.
It was found that the simple disk shape model underestimates the cross-polarized scattering
when compared to more realistic leaf shapes.
References
1. R. Schiffer and K.O. Thielheim, “Light Scattering by Dielectric Needles and Disks,” Journalof Applied Physics, Vol. 50, No. 4, 1979, pp. 2476–2483.
2. M.A. Karam, A.K. Fung and Y.M.M. Antar, “Electromagnetic Wave Scattering from Some
Vegetation Samples,” IEEE Trans. Geosci. Remote Sensing, Vol. 26, No. 6, 1988, 799–808.
3. M.A. Karam and A.K. Fung, “Leaf Shape Effects in Electromagnetic Scattering from Vegetation,” IEEE Trans. Geosci. Remote Sensing, Vol. 27, No. 6, 1988, pp. 687–697.
4. Y. Oh and J.Y. Hong, “Re-examination of Analytical Models for Microwave Scattering from
Deciduous Leaves,” IET Microw. Antennas Propag., Vol. 1, No. 3, 2007, pp. 617–623.
5. G. Cerutti, L. Tougne, A. Vacavant and D. Coquin, “A Parametric Active Polygon for Leaf
Segmentation and Shape Estimation,” Advances in Visual Computing, Springer Berlin Heidelberg, 2011, pp. 202-213
6. P.C. Waterman, “Matrix Formulation of Electromagnetic Scattering,” Proceedings of theIEEE, Vol. 53, No. 8, 1965, pp. 805-812.
7. Hansen, J. E., Spherical Near-Field Measurements, IEE Electromagnetic Wave Series 26,
London U.K., 1988.

- 72 -

F12

Experimental Evaluation of UWB Propagation in a Closed Box
for Replacing Wired Interface Buses in Spacecrafts
Miyuki Hirose and Takehiko Kobayashi
Wireless Systems Laboratory
Tokyo Denki University
5 Senju-Asahi-cho, Adachi-ku, Tokyo 120-8551, Japan
E-mail: miyuki@wsl.c.dendai.ac.jp
A use of ultra wideband (UWB) technology within spacecrafts has been proposed with a view to(partly) replacing wired interface buses with wireless connections [1]. Adoption of wireless technologies
within the spacecrafts could contribute to: (i) reduction of cable weight and launching cost as a results, (ii)reduction in the cost of design, manufacture, and tests, (iii) ensuring reliable robust communication due toeffect of route diversity, (iv) more flexibility in layout of spacecraft subsystems, and (v) more reliable
connections at rotary, moving, and sliding joints. However, multipath propagation in semi-closed
conductive enclosures, such as spacecrafts, restricts digital transmission performance. Narrowband wirelesscommunication systems suffer flat fading in such an environment, hence need a substantial amount offading margin, and only yield a low data rate. On the contrary, UWB is a promising technology because of
its robustness against multipath propagation and low power consumption.
UWB (3.1 – 10.6 GHz, 4.2 – 4.8 GHz, and 7.4 – 7.9 GHz,) propagation and transmission were
measured and characterized inside a 90-liter shield box made of stainless-steel, simulating a miniaturescientific spacecraft. While the continuous wave (6.85 GHz) resulted in nearly 35-dB fading at several“dead spots” caused by multipaths, UWB practically yielded none. The UWB systems have therefore anadvantage over narrowband from the viewpoint of reducing fading margins. A bandwidth over 500 MHz
was capable of reducing the fading depth by approximately 3 dB in conductive semi-closed spaces likespacecraft. No discernible dependence on the distance was observed for UWB propagation gain, delay
spread, or throughput, and no apparent spatial correlation between them.
The effects of inner volume, a piece of radio absorber (paneled on the inner surface), and apertures(perforated on the outer surface of satellites) on UWB propagation and link throughputs were evaluatedinside the box. Increase in the inner volume of boxes resulted in higher UWB propagation gains, but wider
delay spreads. Propagation gain decreased and the fluctuation range of the gain increased when increasing
the area of radio absorber attached on an inner surface. On the other hand, propagation gains were almost
invariable for the total area of apertures normalized by the total area of the inner surface between 0.01 and
0.1%, and gradually decreased with the area beyond 0.1%. The higher UWB propagation gains in thehigh-band UWB were ascribable to higher free space propagation losses. The delay spreads can besuppressed with the use of a small patch of radio absorber and/or with apertures. The throughput for lowand high-band UWB was up to 96 and 100 Mb/s, when the absorber panel covered 4 and 8% of the total
inner surface area, respectively. With apertures, the throughputs were almost invariable for the normalized
area of apertures between 0.05 and 0.4%, and gradually increased with the area beyond 0.4%.
References
1. A. Matsubara, A. Tomiki, T. Toda, and T. Kobayashi, “Measurements and characterization of ultra
wideband propagation within spacecrafts—proposal of wireless transmission for replacing wired
interface buses,” Advances in Spacecraft Technologies, pp. 61-74, IN-TECH, 978-953-307-551-8,
Vienna, Austria, Feb. 2011.

- 73 -

G1

Measurement of Ionosphere Over the Western Pacific Ocean
Mamoru Ishii, Hidekatsu Jin, Tatsuhiro Yokoyama, Takuya Tsugawa, Michi Nishioka,
and Takashi Maruyama
Space Weather and Environment Informatics Laboratory, National Institute of Information andCommunications Technology
4-2-1 Nukui-kita, Koganei, Tokyo 184-8795, Japan
Email: mishii@nict.go.jp

Recently it becomes more important to measure the ionosphere over the ocean. It is huge unobserving
area for a long time and prevents the scientists from establishing precise global model of ionosphere.
Now many of social requirements, e.g., international aviation glows to know the information over theocean.
NICT has a long history of operational ionospheric observation with ionosondes since IGY 1957. On
the beginning, we had four domestic observatories, Wakkanai, Akita, Kokubunji and Yamagawa. Afterthat Akita was closed and Okinawa joined and we operate these four observatories continuously. In
addition, Syowa station in Antarctica has been observing ionosphere by NICT since IGY, too. In
addition as the World Data Center forionosphere, we have a lot of number ofionospheric data obtained by foreign
institutes. There are several methods toobserve/estimate

theionosphericcondition over the ocean; (1) oblique
sounding using ionosondes, (2) transequatorial propagation, (3) occultation
with low orbital satellites, and (4)GNSS buoy. Now we are planning tohaveinternational

cooperationtoestablish the oblique sounding using
ionosondes

andtrans-equatorialpropagation. In addition, we aredeveloping global simulation code
Figure 1. The image of Trans-equatorial Propagation

named “GAIA”, which calculates

between Australia and Japan

climate conditions from the ground to
the ionosphere. In future, we plan tobuild data assimilation system using these simulation code and the ionospheric observation system.

- 74 -

G2

Ionospheric Scintillation Observations by a Digital Beacon Receiver in Tromsø
Yasunobu Ogawa1,2, Yuichi Otsuka3, Mamoru Yamamoto4, and Yoshiyuki Hamaguchi3
1

National Institute of Polar Research, Japan
Email: kazuya@tamacc.chuo-u.ac.jp
2
Graduate University for Advanced Studies, Japan
3
STEL, Nagoya University, Japan
4
RISH, Kyoto University, Japan
Digital beacon receiver observations in Tromsø Norway have been conducted since October 2011. The
digital beacon receiver consists of the open-source hardware called Universal Software Radio Peripheral(USRP) and the open-source software toolkit for the software deﬁned radio (GNU Radio) [Yamamoto,2008]. The digital beacon receiver receives radio beacon signals transmitted from Low Earth Orbit
Satellites (LEOS) and enables us to measure the ionospheric total electron content (TEC) and ionosphericscintillation. We have investigated a detailed relationship between ionospheric scintillation and auroral arcs
by using the digital beacon receiver and all-sky/narrow field-of-view imagers with high time resolution(1-4 Hz). Our initial results show that the ionospheric scintillation frequently occurs at the boundary of anauroral arc and in the region between auroral arcs, not in an auroral arc. The results imply that enhanced
electric fields near auroral arcs may contribute to the formation of the ionospheric scintillation.

Figure 1: Time variation of ROTI (the standard
deviation of the rate of change of TEC) derived from
TEC data of the digital beacon receiver in Tromsø on
January 23, 2012.

Figure 2: A selected all-sky image together
with the locations of the beacon satellite(marked ‘+’) and GPS satellites (marked ‘×’)obtained at 01:52:52.1 UT on January 23, 2012.
The ROTI increased at 01:52:48 and 01:52:52
UT when the beacon satellite was located nearauroral arcs.

Reference
Yamamoto, M., Digital beacon receiver for ionospheric TEC measurement developed with GNU Radio,
Earth Planets Space, 60, e21–e24, 2008.

- 75 -

G3

Development of Digital Receiver for LF/MF Band Radio Wave Onboard Sounding
RocketTakafumi Mizuno
Faculty of Engineering, Toyama prefectural University
5180 Kurokawa, Imizu-shi, Toyama 939-0398, Japan
Email: mizuno@rdw.pu-toyama.ac.jp
Yuki Ashihara
Department of Electrical Engineering, Nara National College of Technology
22 Yata-cho, Yamatokoriyama, Nara, Japan
Email: ashihara@elec.nara-k.ac.jp
Keigo Ishisaka
Faculty of Engineering, Toyama prefectural University
5180 Kurokawa, Imizu-shi, Toyama 939-0398, Japan
Email: ishisaka@pu-toyama.ac.jp

MF absorption method is the only technique to obtain electron density profiles such as an ionospheric D
region. This method estimates an electron density profile by analyzing the propagation characteristics ofradio waves in the region from the ground to the lower ionosphere. The LF/MF band radio receiver isnecessary for observing the broadcast waves in the ionosphere. The LF/MF band radio receiver has everconsisted of analog circuit using super-heterodyne system. This receiver has been loaded on many sounding
rockets and observed the propagation characteristics of radio waves in the ionosphere. And we could obtainthe frequency components indicated by the propagation characteristics, such as the downward wave and the
upward wave and the left and right handed waves by using the frequency analysis.
In this study we present how to obtain the specific frequencies of radio waves in the ionosphere using thedigital radio receiver with the FPGA so that the receiver can be miniaturization and high-performance.
Especially, we investigate that we obtain the specific frequencies, such as Doppler shift frequency, upwardand downward wave frequency and so on, of many radio waves with the different frequency each other byusing the digital radio receiver onboard the sounding rocket. Then it's investigated whether it's possible toobtain the horizontal structure of the electron density in the lower ionosphere by using a lot of radio data ofdifferent path observed along the rocket trajectory.

- 76 -

G4Electron Density Estimation by Using Simultaneous Observations of Lightning
Optical Emissions and Whistlers from ISS GLIMS Mission
Katsunori Suzuki1, Y.Hobara1, 2, K.Kakinuma1, I.Linscott3, U.Inan3, M.Sato4, Y.Takahashi4,
T.Ushio5, Z.Kawasaki5, T.Morimoto6, A.Yamazaki7, M.Suzuki7
1

Department of Communication Engineering and Informatics,
The University of Electro-Communications, Tokyo, Japan.
2

Earth Environment Research Station,

The University of Electro-Communications, Tokyo, Japan
1-5-1 Chofugaoka, Chofu city, Tokyo 182-8585, Japan.
3Department of Electrical Engineering, Stanford University.
4Department of Cosmosciences, Graduate School of Science, Hokkaido University.
5Information and communication engineering department, Osaka University.
6Faculty of Science and Engineering, Kinki University.
7Institute of Space and Astronautical Science, Japan Aerospace Exploraton Agency.

Email: s1531053@edu.cc.uec.ac.jp, hobara@ee.uec.ac.jp
The atmospherics from the lightning penetrating through the ionosphere and are observed as socalled lighting whistlers in the magnetosphere. Since ionosphere is the dispersive medium, groupdelays of whistler waves with different frequencies contain the information on the electron density
along the propagation path.In this paper we estimate the maximum electron density of ionospheric F2
layer (NmF2) by using the time delay between the lightning optical emissions and associated whistlers,
which were simultaneously measured by Global Lightning and SprIte MeasurementS (GLIMS)
mission onboard ISS. We found that the calculated NmF2 by using group time delay are in good
agreement with those estimated from classical dispersion analysis of lightning whistlers, IRI modeland ionosonde for relatively small dispersion events. This method will be useful to cover the area
where no ground-based measurements are available such as over the ocean and remote areas.

- 77 -G5

VHF Lightning Observations and DOA Estimation from ISS / JEM-GLIMS
Takeshi Morimoto
Department of Electric and Electronic Engineering, Faculty of Science and Engineering, Kinki University
3-4-1 Kowakae, Higashiosaka City, Osaka 577-8502, Japan
Email: morimoto@ele.kindai.ac.jp
Hiroshi Kikuchi, Tomoo Ushio
Division of Electrical, Electric and Information Engineering, Graduate School of Engineering, Osaka
University
Mitsuteru Sato
Department of Cosmosciences, Graduate School of Science, Hokkaido University
Atsushi Yamazaki, Makoto Suzuki
Department of Solar System Science, Institute of Space and Astronautical Science, JAXA
Yasuhide Hobara
Department of Communication Engineering and Informatics, Graduate School of Informatics andEngineering, The University of Electro-Communications

Global Lightning and sprIte MeasurementS (GLIMS) mission has been conducted on Exposed Facilityof Japanese Experiment Module (JEM-EF) of the International Space Station (ISS) for more than 2 and halfyears. It is a space mission to observe global distributions of lightning and lightning-associated TLEs by
combining observations with radio and optical sensors. This paper focuses on an EM payload of
JEM-GLIMS mission, VITF, and its primary observation results.
VITF consists of two sets of antennas, band-pass filters, amplifiers, and 2-channel-AD-converter.
Impulsive EM radiations received by the antennas are digitized by the AD converter synchronizing withanother channel through the filters and the amplifiers. A patch type antenna is developed within the size of
200*200 mm. It is mounted on the antenna base made of aluminum alloy and Teflon block with the totalheight of 100 mm to gain its bandwidth and to reduce the interference from other structural objects. The
same two units of antennas are installed with the separation of 1.6 m. Their bandwidths with the higher
return loss than -3 dB are from 70 to 100 MHz. The signals received by the antenna are transmitted alongcables with the same lengths to the electronics. The AD converter records 130 waveforms as maximum of
one dataset with the duration of 2.5 μs with 200 MS/s. The developments of VITF are based on the heritageof VHF sensor on Maido-1 satellite. It is designed to estimate the direction-of-arrival (DOA) with about 10
km resolution that is equivalent to the scale of a thundercloud. It means that VITF is able to monitor
thunderclouds with lightning activities globally. Comprehensive observations with the optical payloads of
JEM-GLIMS mission for lightning activities and TLEs are expected to give us many scientific impacts tothe field.
JEM-GIMS mission payload was successfully launched at the end of July 2012, and transported andinstalled to the ISS. After the initial checkout and maintenance, its nominal operation has been continuingfrom December 2012. Through the operation period, VITF corrects numerous VHF EM data synchronizedwith optical signals. 2,904 events with active VHF radiations out of 4,173 optical observations, about 70%,
are observed for the 20 months in 2013 and 2014. The DOA estimations of the received VHF pulses areattempted using the broadband digital interferometry. Some results agree with the optical observations,even though DOA estimation has difficulties caused by its very short baseline of the antennas and multiple
pulses in short time, namely burst-type EM waveforms. It is the worlds’ first space-borne lightning location
by EM DOA estimation. The DOA estimation via ionospheric dispersion is also applied the received VHF
EM pulses. The simulated dispersion curves for various DOA with lay-tracing method and for received
waveform enables to DOA estimation.

- 78 -

G6

Beacon Experiment of Low-latitude Ionosphere from Southeast Asia
Mamoru Yamamoto
Research Institute for Sustainable Humanosphere (RISH)
Kyoto University
Gokasho, Uji, Kyoto 611-0011, Japan
Email: yamamoto@rish.kyoto-u.ac.jp

We have been successfully conducted observations of total-electron content (TEC) of the ionosphere by
the satellite-ground beacon experiment. An unique dual-band (150/400MHz) digital receiver GRBR (GNURadio Beacon Receiver) were developed based on the recent digital signal processing technologies. The
GRBR receivers were deployed first in Japan, and then in southeast Asia. Data from the GRBR network has
been used for the investigations of variety of ionospheric phenomena. Especially in low-latitude region,
longitudinal “large-scale wave structures (LSWS)” were studied in detail as a source of equatorial Spread-F(ESF) events. Also we were successful to measure the equatorial ionospheric anomaly (EIA) near 100E
longitude in large latitudinal extent of at most +/-20 degrees around the geomagnetic equator. The
technique is utilized for sounding rocket-ground experiment as well. Now there are several plans to launchseveral beacon satellites in the low-latitude region in 2015/2016, and a new project to develop new series ofdigital receivers has started this year. We review the ionospheric studies with the GRBR network, anddiscuss future direction of the related studies.

- 79 -

G7

GNSS scintillations and ambient ionization near the northern EIA crest
under superstorm event of March, 2015
D. Jana*2 and S. K. Chakraborty*1*

Raja Peary Mohan College, Uttarpara, Hooghly, West Bengal-712258, India
1

Also at: Maharaja Srischandra College, Kolkata - 700003

Email- *1skchak2003@yahoo.com, *2janadebasis69@gmail.com

Abstract - Observations of GNSS scintillations and TEC from RPM College centre(22.655oN, 88.348oE), situated near the northern EIA crest, in the main phase of recent
superstorm event of March 17-19, 2015 reveal opposite features contrary to quiet days, in thepost sunset variability of TEC at two nearby longitudes 88 oE and 90oE and occurrence of
scintillation. On the 2nd day TEC at both the sectors exhibits extremely low diurnal values in
conjunction with the whole day counter-electrojet event. PCA analysis of recovery phase
TEC and EEJ reflect a good correspondence suggesting dominant contribution of disturbance
dynamo field in ionospheric variability.

- 80 -

G8

Estimation of Electron Density Using the Propagation Characteristics ofRadio Waves by S-520-29 Sounding Rocket
Keita Itaya, Keigo Ishisaka
Faculty of Engineering, Toyama prefectural University
5180 Kurokawa, Imizu-shi, Toyama 939-0398, Japan
Email: itaya@rdw.pu-toyama.ac.jp, ishisaka@pu-toyama.ac.jp
Yuki Ashihara
Department of Electrical Engineering, Nara National College of Technology
22 Yata-cho, Yamatokoriyama, Nara, Japan
Email: ashihara@elec.nara-k.ac.jp
Junichi Kurihara
Hokkaido University
5 Hatijonishi, Sapporo-shi, Hokkaido, 060--0808 Japan
Email: kurihara@mail.sci.hokudai.ac.jp
Takumi Abe
Japan Aerospace Exploration Agency (JAXA), Institute of Space and Astronautical
Science (ISAS)
3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa, 252-5210 Japan
Email: abe.takumi@jaxa.jp
S-520-29 sounding rocket experiment was carried out at Uchinoura Space Center (USC) at
19:10 JST on 17 August, 2014. The purpose of this sounding rocket experiments isobservation of sporadic E layer that appears in the lower ionosphere at near 100km. This
rocket is equipped with LF/MF band radio receiver in order to observe of propagation
characteristics of LF/MF band radio waves in rocket flight. LF/MF band radio receiverreceives three radio waves of 873kHz (JOGB), 666kHz (JOBK), 60kHz (JJY) from theground. In the sounding rocket experiment, LF/MF band radio receiver was working
properly. We have completed the observation of radio wave intensity. We analyze theobservation results using a Doppler shift calculations by frequency analysis. From theanalysis, we obtained the propagation characteristics of LF/MF band radio waves andestimated the electron density of the ionosphere. As a result, 873kHz, 666kHz radio waves
are reflected by the ionosphere. 60kHz wave was able to propagate in ionosphere because
wavelength of 60kHz was longer than the thickness of the sporadic E layer.
In this study, we explain the result of LF/MF band radio waves observations and theelectron density of the ionosphere using frequency analysis by S-520-29 sounding rocket.

- 81 -

H1

A simple kinetic-fluid model for ion-scale waves in collisionless plasmas
Yasuhiro Nariyuki
Faculty of Human Development, University of Toyama
3190, Toyama City, Toyama 930-8555, Japan
Email: nariyuki@edu.u-toyama.ac.jp
Sinji Saito, Takayuki Umeda
Solar-Terrestrial Environment Laboratory, Nagoya University
Takuma Nakamura
Space Research Institute (IWF), Austrian Academy of Sciences (OeAW)In this presentation, we discuss a simple kinetic-fluid model for obliquely propagating ion scale
waves, which approximately includes the effects of ion kinetics and electron inertia. In the present
study, the specific heat ratio of the Vlasov-MHD model[1] is applied to a two-fluid model with isotropicpressure[2]. The linear dispersion relation of the system is[(M ω 2 − k 2 )T − CS2 k 2 ky2 ](M ω 2 − k 2 cos2 θbk ) = ω 2 k 4 cos2 θbk T .

(1)where T = ω 2 − CS2 k 2 , CS2 = βf = γi Ti + Te , M = 1 + me k 2 , and k 2 = kx2 + ky2 . γi is the specificheat ratio defined by the past study[1].
Figure 1 shows (a)ωr = Re[ω] and (b)−ωi of the shear mode (kinetic Alfvén waves ([4])) in the
Vlasov-MHD model (1) and the full-Vlasov theory ([3]) in the case of me = 1/1836, θbk = 89◦ andβe = βi = 2.0. The agreement of the phase velocity (ωr /k) and the growth rate among the models
is observed when k <∼ 5. On the other hand, in contrast to the Landau fluid model ([4]), −ωi of the
Vlasov-MHD model decreases with increasing k at k >∼ 5. This may be due to the lack of the finite ion
Larmor radius effects and correct evaluation of the transit time damping.
This work was supported by MEXT/JSPS under Grant-in-Aid for Scientific Research (Nos. 26287119,
26287041, 15K17770, and 15K13572).

Figure 1. (a)ωr = Re[ω] and (b)−ωi in the Vlasov-MHD model (1) (solid lines) and the full-Vlasov
theory ([3]) (dashed lines).

References
1. Y. Nariyuki, S. Saito, and T. Umeda, Phys. Plasmas, 20, 072115, 2013.
2. T. K. M. Nakamura, D. Hayashi, M. Fujimoto, and I. Shinohara, Phys. Rev. Lett., 92, 145001, 2004.
3. Stix, T. H., Waves in plasmas, by American Institute of Physics, New York, 1992.
4. T. Passot,P. L. Sulem, and P. Hunana, Phys. Plasmas, 19, 082113, 2012.

- 82 -

H2

Collective Thomson scattering diagnostics in non-equilibrium plasma:
Application to high power laser experiment
Shuichi Matsukiyo1, Youichi Sakawa2, Kentaro Tomita1, Yasuhiro Kuramitsu3, Hideaki Takabe2
1
Faculty of Engineering Sciences, Kyushu University
6-1 Kasuga-Koen, Kasuga, Fukuoka 816-8580, Japan
Email: matsukiy@esst.kyushu-u.ac.jp
2
Institute of Laser Engineering, Osaka University, Japan
3
National Central University, Taiwan

The collective Thomson scattering is widely used to measure local quantities of a plasma in laboratory
experiment. It is the application of the elastic scattering of low frequency electromagnetic waves by freeelectrons in a plasma [1]. In the measurement one observes scattered waves which are radiated through theinteraction between an incident probe laser light and collective motions of electrons in a plasma. The
characteristics of the scattered waves enable us to infer local electron density, plasma temperature, andvalence of ions. Details of the scattering theory is complex. In particular the theory of the collective
Thomson scattering in a non-equilibrium plasma has not been established well. Here, a high power laser
experiment of a collisionless shock performed at the Institute of Laser Engineering, Osaka university, istaken as an example (Figure 1). The wave equation of the scattered waves interacting with a plasma
containing arbitrary density fluctuations is numericallysolved to discuss how the non-equilibrium plasma in a
shock transition region scatters an incident probe laser
light. In the past collective Thomson scattering
measurement in a laboratory plasma the so-called ion
feature is mainly paid attention, since the spectral
intensity of scattered waves is relatively high. Here, wesuggest that the electron feature, whose spectral intensity
is usually very weak, can be observed in a highly
Figure 1. Collisionless shock experiment at ILE.
non-equilibrium plasma as in the shock transition region.

Reference
1. D. H. Froula, S. H. Glenzer, N. C. Luhmann Jr., and J. Sheffield, "Plasma Scattering ofElectromagnetic Radiation", Elsevier, 2014.

- 83 -

H3

Dependencies of the Generation Process of Whistler-Mode Emissions on Temperature
Anisotropy of Energetic Electrons in the Earth’s Inner Magnetosphere
Y. Katoh1, Y. Omura2, Y. Miyake3, H. Usui3, and H. Nakashima4
Graduate School of Science, Tohoku University, 6-3 Aramaki-Aza-Aoba, Aoba, Sendai, Miyagi, Japan,
yuto@stpp.gp.tohoku.ac.jp
2
Research Institute for Sustainable Humanosphere, Kyoto University, Gokasho, Uji, Kyoto 611-0011, Japan
3
Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
4
Academic Center for Computing and Media Studies, Kyoto University, Yoshida-Honmachi, Sakyo-ku, Kyoto
606-8501, Japan
1

Whistler-mode chorus emissions are electromagnetic plasma waves commonly observed in planetary
magnetospheres. In the Earth's inner magnetosphere, chorus emissions are observed mostly on the dawnside and are enhanced during geomagnetically disturbed periods. Chorus emissions appear in the typicalfrequency range from 0.2 to 0.8 fce0 with a gap at the half fce0, where fce0 represents the electron
gyrofrequency at the magnetic equator. Recent in situ observation in the magnetosphere revealed thepresence of whistler-mode hiss-like emissions, whose wave amplitude is comparable to those of chorus
emissions.
By a series of electron hybrid simulations, we study dependencies of the generation process ofwhistler-mode chorus and hiss-like emissions on temperature anisotropy of energetic electrons. The
generation process of chorus has been reproduced in electron hybrid simulations and has been explained by
the nonlinear wave growth theory [see review by Omura et al., 2012]. The generation mechanism ofhiss-like emissions is also explained by the nonlinear wave growth theory and has been reproduced bysimulations [Katoh and Omura, 2013]. In the present study, by an improved electron hybrid code with
OhHelp library [Nakashima et al., 2009], we conduct a series of electron hybrid simulations for different
temperature anisotropy (AT) of the initial velocity distribution function of energetic electrons. We vary AT
in the range from 3 to 9 with changing the number density of energetic electrons (Nh) so as to study
whether distinct rising-tone chorus emissions are reproduced or not in the assumed initial condition.
Simulation results reveal that Nh required for the chorus generation decreases as AT of energetic electrons
increases and that reproduced spectra becomes hiss-like for large Nh cases. These simulation results clarify
the validity of the nonlinear wave growth theory and suggest that the nonlinear wave-particle interactions
play crucial roles in the generation process of whistler-mode chorus and hiss-like emissions in themagnetosphere.

References
1. Y. Katoh and Y. Omura, “Effect of the background magnetic field inhomogeneity on generation
processes of whistler-mode chorus and hiss-like broadband emissions”, J. Geophys. Res. Space Physics,
118, pp. 4189-4198, 2013.
2. H. Nakashima, Y. Miyake, H. Usui, and Y. Omura, "OhHelp: A Scalable Domain-Decomposing
Dynamic Load Balancing for Particle-in-Cell Simulations", Proceedings of the 23rd International
Conference on Supercomputing, New York, 90-99, 2009.
3. Y. Omura, D. Nunn, and D. Summers, Generation process of whistler mode chorus emissions: Current
status of nonlinear wave growth theory, in Dynamics of the Earth's Radiation Belts and Inner
Magnetosphere, edited by D. Summers et al., 243–254, Geophys. Monogr. Ser., vol. 199, AGU,
Washington, D. C., 2012.

- 84 -

H4

Study on Improvement of Direction Finding Method for VLF Waves
Mamoru Ota ; Yoshiya Kasahara ; Yoshitaka Goto
Graduate School of Natural Science and Technology
Kanazawa University
Kanazawa-shi, Ishikawa 920-1192, Japan
Email: ota@cie.is.t.kanazawa-u.ac.jp

Investigating characteristics of plasma waves observed by scientific satellites in the Earth’s
magnetosphere and plasmasphere is the effective key to understand not only generation mechanisms of the
waves but also a plasma environment which influences its generation and propagation conditions. In
particular, direction finding of plasma waves is most important for understanding these characteristics.
The wave distribution function (WDF) method was proposed for VLF waves in the Earth’s
magnetosphere/plasmasphere [1]. This method assumes that the observed signals are combinations of acontinuum of superposed plane waves of different frequencies, propagating in different directions with nomutual phase coherence. In addition, this method also assumes that the observed signals are stationary andfollow ergodic Gaussian random process with zero mean. Under these assumptions, the WDF method canestimate a WDF as directional distribution of wave energy density by using a spectral matrix which
composes by cross spectra of observed signals. The WDF method is preferred when a wave source does notsatisfy the plane-wave approximation or when it is widely extended. However, the WDF estimation is
ill-posed problem, that is, the solution is not determined uniquely. Models as additional information for
WDF must be needed to determine the solution uniquely. Many models have been proposed until now suchas the Gaussian distribution model, and Markov random field model [2]. The estimation using these models
works well if the sample number of observed signals is large enough to calculate spectral matrices exactly.
Actually the number of sample observed by satellites is very few. We therefore must take into account thatthe spectral matrix which can be used for WDF estimation contains uncertainty.
We studied two issues for finding solution of above-mentioned problems. First, to treat the uncertainty,
we used the Bayesian inference, and we introduced probability density distribution which determines
relationship between observed and theoretical spectral matrices. This introduction makes the WDF
estimation enable to take into account the effect of sample number. Second, we studied Markov chain
Monte Carlo (MCMC) methods as a Bayesian inference method for WDF estimation. Variational Bayesian
method is also well known as a Bayesian inference method for complicated probability model. In thismethod, we need to get a proper class of approximate probability distribution. However, we cannot get it
due to the complexity of probability distribution we set. We therefore must select MCMC methods. MCMC
methods have many types such as Metropolis-Hastings, Gibbs sampler, Hamiltonian Monte Carlo (HMC),
and so on [3]. We investigated which type is better in terms of reliability and efficiency.

References
1. Storey, L. R. O., and F. Lefeuvre, “The analysis of 6-component measurements of a random
electromagnetic wave field in a magnetoplasma I. The direct problem,” Geophysical JournalInternational, 56(2), pp. 255-268, 1979.
2. M. Ota, Y. Kasahara, and Y. Goto, “A new method for direction finding based on Markov random fieldmodel,” Radio Science, July 2015.
3. Bishop, C. M., “Pattern Recognition and Machine Learning,” Springer-Verlag New York, Inc.,
Secaucus, NJ, USA, 2006.

- 85 -

H5

Effect of Electron Core-Halo Drift Speed on the Strahl Formationin the Solar Wind: Particle-in-cell Simulations
Jungjoon Seough
Faculty of Human Development
University of Toyama
3190 Gofuku, Toyama city, Toyama 930-8555, Japan
Email: seough@edu.u-toyama.ac.jp
The halo is present at all pitch angles, whereas the strahl appears as a magnetic-field-aligned beam
streaming away from the Sun. Recently, it was found that the strahl could be generated through nonlinear
pitch angle scattering resulting from the interactions between the whistler waves excited by anisotropic
core and the drifting halo electrons. The new scenario of local formation of the strahl was demonstrated by
the particle-in-cell simulation. One of key parameters used in the simulation model is a relatively large
core-halo drift speed. Such a drift is typically observed in the collisionless solar wind plasmas at 1 AU andits speed extends up to a few tens of the local Alfvé n speed. In the present work, a parametric study isconducted to investigate how the core-halo drift speed affects the formation of the strahl. It is found that thethreshold drift speed is necessary for effectively generating the stahl through pitch angle scattering ofdrifting halo when the whistler instability occurs in interplanetary space.

- 86 -

H6

Low Frequency Characteristics of a Wire Antenna aboard a Scientific Spacecraft
Tomohiko Imachi1, Ryoichi Higashi2, Mitsunori Ozaki1, Satoshi Yagitani1
1
Kanazawa University
Kakuma-machi, Kanazawa, Ishikawa 920-1192, Japan
Email: imachi@imc.kanazawa-u.ac.jp
2

Ishiakwa National College of Technology

Observation of plasma waves in space is an important task for scientific satellites. In order to accurately
determine the magnitude of the electric field of plasma waves observed by a scientific satellite, the calibrationof the sensor must be accurate. For the observation of the electric field, a wire antenna which is a dipole
antenna using a wire is often used. The effective length is one of the needed characteristics, however, it isdifficult to acquire it in a ground test.
In order to study the estimation of the effective length, we have performed a rheometry experiment. In
this experiment, a signal applied to two metal plate in water generates a quasi-static field between them
(Figure 1.). The field is considered as the electric field of an electromagnetic wave and receiving it by anantenna, we can measure its output voltage. The magnitude of the electric field is known, so that we cancalculate the effective length of the antenna from its output voltage. The frequency range of the experiment
we have performed is from 10 Hz to 100 kHz.
And we have made a theoretical calculation method which is called DPS (Distributed Parallel Sources)
method. In this calculation method, the electric field is approximated to a parallel circuit of potential sources.
According of this method, we can make an equivalent circuit of the experiment, and calculate the value ofthe effective length as a function of the frequency. We have applied to this method to the experiment we have
performed, and the calculation result is consistent to the experiment result.

Figure 1. Overview of rheometry experiment

Figure 2. Basic concept of the equivalent circuit

In the present research we have investigate frequency dependence of the effective length of simple wire
antennas, and found that the insulator which covers the side of wire cause frequency dependence at low
frequencies. In this time we made a structure model of the actual satellite “ERG”, and investigate the effectof the deployment structure such as solar battery panels to the wire antenna. We will show the detail at themeeting.
References
1. T. Imachi, et al., “Rheometry Experiment for a Wire Antenna Aboard Spacecraft at Low
Frequencies”, Vol. J89-B, No. 4, pp. 552–559, April 2006.

- 87 -

H7

Effect of Lunar Surface Reflection on Polarizations of Auroral Kilometric Radiation
Observed by KAGUYA
Yoshitaka Goto1, Ryota Kimura1, Yoshiya Kasahara1, Atsushi Kumamoto2
1
School of Electrical and Computer Engineering, Kanazawa University
Kakuma-machi, Kanazawa 920-1192, Japan
Email: ygotou@is.t.kanazawa-u.ac.jp
2
Department of Geophysics, Graduate school of Science, Tohoku University
6-3, Aramaki Aza-Aoba, Aoba-ku, Sendai 980-8578, JapanThe lunar orbiter KAGUYA frequently observed an electromagnetic wave from the Earth that is called
AKR (auroral kilometric radiation). One of the interesting features of AKR observation from KAGUYA is
that the observed AKR consists of not only waves reached directly from the Earth's auroral region but also
waves reflected on the lunar surface [1, 2]. Reflection characteristics are determined by electricalparameters of the lunar surface. These parameters are, however, still unknown especially for low
frequencies (LF band) although they may bring us a new knowledge about materials of the lunar surface.
The reflection influences polarizations of the observed AKR. Since KAGUYA is a three-axis stabilized(non-spinning) satellite and obtains only electric field on an antenna plane, measured polarizations dependon arrival direction of the AKR. When k-vector of the arrival AKR is perpendicular to the antenna plane,
the measured polarization is equivalent to the true AKR polarization. On the other hand, when the k-vector
is parallel to the antenna plane, the measured polarization is always linear regardless of the wave
polarization. KAGUYA always keeps its antenna plane parallel to the lunar surface along the polar orbit of
the Moon at an altitude of about 100 km. Thus, polarization should show a sinusoidal variation along thesatellite movement if the reflection did not occur. In the statistical analysis of the observed polarization,
however, the variation does not show pure sinusoidal curves due to the effect of the reflection.
In the present study, we examined whether the effect of the lunar surface reflection on the AKR
polarization can be explained by a theoretical reflection model that is based on Fresnel equations. From thecomparison between observed polarizations and those calculated for various parameter sets of original
AKR polarization and the lunar surface, we obtained following results; (1) the original AKR polarization isperfectly circular, and (2) dependence of observed polarizations on incident angle is reconstructed under
the condition that the permittivity of the lunar surface takes values from 4 to 8. The former result is
consistent with observations of AKR from Cassini. As for the permittivity, it should vary depending on
regions while we assume that it is constant all over the regions of the same incident angle. In that sense, thepermittivity is an average value for each incident angle. In the future study, we will examine its regional
dependences through detailed data analysis and compare them with permittivity map for high frequencies(HF band) that are derived from measurements by the lunar radar sounder instrument (LRS) on KAGUYA.
References
1. T. Ono, A. Kumamoto, Y. Kasahara, Y. Yamaguchi, A. Yamaji, T. Kobayashi, S. Oshigami, H.
Nakagawa, Y. Goto, K. Hashimoto, Y. Omura, T. Imachi, H. Matsumoto, and H. Oya, “The Lunar
Radar Sounder (LRS) Onboard the KAGUYA (SELENE) Spacecraft,” Space Sci. Rev., 154, 145–192,
2010.
2. Y. Goto, T. Fujimoto, Y. Kasahara, A. Kumamoto, and T. Ono, “Lunar ionosphere exploration method
using auroral kilometric radiation,” Earth Planets Space, 63, 47–56, 2011.

- 88 -

H8

Development of Automatic Detection of Omega Signals
Captured by PFX Onboard Akebono
IMAD. Suarjaya*, Y. Kasahara** and Y. Goto***Graduate School of Natural Sci. and Tech., Kanazawa University,
Kakuma-machi, Kanazawa 920-1192, Japan
E-mail: *suarjaya@cie.is.t.kanazawa-u.ac.jp,**kasahara@is.t.kanazawa-u.ac.jp, ***ygotou@is.t.kanazawa-u.ac.jp
The Akebono satellite was launched in 1989 to observe the Earth’s magnetosphere andplasmasphere. Wave normal and Poynting flux Analyzer (PFX) is a waveform receivermeasuring two components of electric field and three components of magnetic field with bandwidth of 50 Hz in a frequency range from 100 Hz to 12.75 kHz [1]. Omega was a navigation
system transmitted from 8 different ground stations with transmission pattern every 10 s. Eachstation transmitted a different pattern of frequency but has common frequency at 10.2 kHz. The
Omega system was terminated in 1997 giving way to the GPS system. The PFX onboard
Akebono has received signals at 10.2 kHz from these stations from 1989 to 1997, and hugeamounts of PFX data is valuable to study propagation characteristics of VLF waves in theionosphere and plasmasphere around the Earth. In the present paper, we introduce detection
method of Omega signals from the PFX data in a systematic way and show statistical featuresof Omega signals observed in the plasmasphere.
We first detect the raise time of each signal by comparing average intensity of specifictime frame with threshold level expecting sudden increase of intensity. Second, we determine
the transmission station by comparing the raise time and the transmission patterns of 8 Omega
stations. Then we calculate the delay time of the signal subtracting the transmission time of thestation. Next we calculate the intensity of the signal both for electric and magnetic components.
The software we develop can analyze and show result of intensity and delay time map in
geomagnetic or geographic map for both magnetic and electric field. Analyses of local timedependence is also available. In this study, we developed advanced method for automatic
detection to process huge amounts of several years PFX’s data. This method enables us todistinguish noises and also handle the error detection to produce accurate result and finding.
References
1. I. Kimura, K. Hashimoto, I. Nagano, T. Okada, M. Yamamoto, T. Yoshino, H. Matsumoto,
M. Ejiri and K. Hayashi, “VLF Observation by the Akebono (EXOS-D) Satellite”, J.
Geomag. Geoelectr., 42, 1990, pp. 459-478.

- 89 -

J1

Feasible Spectral Indices Measurement by Broad-band VLBI in VGOS Era
Kazuhiro Takefuji
National Institute of Informations and Communication Technology
893-1 Hirai,Kashima, Ibaraki 314-8501, Japan
Email: takefuji@nict.go.jp

To establish next generation geodetic VLBI, called VGOS ( VLBI Global Observing System ) which isbased on the so-called broadband delay which uses four or more frequency bands in the range from 2.5
GHz to 14 GHz(1), we have been developing a broad-band system for Kashima 34 meter radio telescope
which is the third largest radio telescope with cassegrain optics in Japan.
Currently we have installed two type of broadband feed (IGUANA-H and NINJA). Firstly for 6.5 GHz to
15 GHz since December 2013 and secondary 3.2 GHz to 14.4 GHz since July 2015. VGOS specification is
2.5GHz, but we designed the NINJA feed from 3.2 GHz to prevent highly RFI circumastances in Japan.
Our broad-band system with ambient receiver system having the system temperature is about 150 Kelvin.
The signal after the LNA transfers via optical fiber, then high-speed direct sampler called K6/GALAS
converted analog to digital without frequency conversion.
By the way, full VGOS radio telescope is built by The Geospatial Information Authority of Japan in
Ishioka, Ibaraki, Japan. The feasible broad-band VLBI sessions with Kashima 34 m and new VGOS-type
Ishioka 13 m were carried out on January and July 2015. We carried out intensive sessions for 6 to 14GHz(6-bands, 2048Msps, 1bit). Fringes from every six frequencies could be simultaneously obtained after
software correlation with GICO3 (2). Then a coherent phase connection over six frequencies were
performed. An error of delay resolution function by the super bandwidth synthesis is estimated only 0.1 ps
in 1 sec integration. Since we obtain 6 GHz bandwidth correlation result, we try to measure the spectral
indices on broad-band correlation by comparing radio sources in these sessions. We will report more detail
information in the meeting.

References
1. International VLBI Service, NASA Website, http://ivscc.gsfc.nasa.gov/technology/vgos-general.html.
2. M.Kimura,”Development” of the software correlator for the VERA system II”, IVS NICT-TDC News.
No.28, pp.22-25, 2007.

- 90 -

J2

Development of Wide-Band Bandwidth Synthesis Software for Geodetic
VLBI
Tetsuro Kondo1 , Kazuhiro Takefuji1 , and Yoshihiro Fukuzaki2
1 Kashima Space Technology Center
National Institute of Information and Communications Technology
893-1 Hirai, Kashima, Ibaraki 314-8501, Japan
Email: kondo@nict.go.jp
2 Geospatial Information Authority of Japan

Wide-band bandwidth synthesis (WBWS) is a technique to combine multi-band VLBI dataof which bandwidth exceeds 10 GHz or more (band width is about 10 times wider than that ofconventional bandwidth synthesis technique) not only to improve a delay resolution but also toincrease a sensitivity of fringe detection. We have been developing WBWS software and havesuccessfully applied it to actual wide-band VLBI data.
Three corrections considered in the development of WBWS software are as follows: 1) innerband phase correction, 2) inter-band delay correction, and 3) ionospheric delay correction. As
for corrections 1) and 2), a realistic method has been developed which corrects the phase anddelay based on those obtained by a reference scan (which is usually defined as a scan observing astrong radio star). As for an ionospheric correction, it is not implemented yet in the software, but
the total electron content (TEC) will be estimated by using the phase characteristics against thefrequency. WBWS is applied to a wide-band VLBI observation data obtained by an experimentconducted on Kashima-Ishioka baseline in Jan. 16, 2015 and got a successful result (Figure 1).

Figure 1: An example of so-called search function after WBWS for an observation on Jan.16,
2015 on Kashima-Ishioka baseline.

- 91 -

J3

Development of Wideband Feed
Hideki Ujihara
National Institute of Information and Communications Technology
893-1 Hirai, Kashima, Ibaraki 314-8501, Japan
Email: ujihara@nict.go.jp
The Author developed wideband feeds for Kashima 34m antenna and small movable VLBI stations“MARBLE”. These feeds are used for our study on Time and Frequency transfer with VLBI, which
compares atomic clocks separated in long distance. However our wideband antenna can be used for otherapplications such as radio astronomy, geodesy or other observation needs wider bandwidth never done
before.
This project of Kashima VLBI group was named Gala-V and covers 2.2-14.4GHz or more over.
Frequency coverage of our Gala-V system, shown in Table 1, was defined to avoid RFI and minimize thedata set of our VLBI observation. We hope our new technologies will be used in SKA(Square Kilometer
Array), VGOS(VLBI Global Observing System) or other next generation projects. SKA and VGOS uses
smaller wideband feed than us now, because they cannot develop wideband feeds for conventionalcassegrain antennas which needs feeds with narrower beam width and larger aperture size. Thus they have
limited design strategy, however we can arrange the beam size of feed for cassegrain antenna or center
focus of the parabola or other optics in our design capabilities.
Now two types of wideband feeds were installed in 34m antenna. The first is a prototype of IGUANA
feed designed for 2.2-18 or up to 22GHz. IGUANA feed was named after Gala-V, Gala-V was named fromunique evolutions in Galapagos, off course, V is for VLBI. This prototype is designed for higher frequency
band of IGUANA, thus it is called IGUANA-H and now used in 6.5-15.0GHz with 30-50% aperture
efficiency of 34m antenna, which is shown in right of Figure 2. First prototype was set in 34m antenna at
the end of 2013, and the second replaced it in 2014 and is used still now. The Author previously designed
multimode horns for ASTRO-G/VSOP-2 satellite, 6.7GHz feeds for 20m antennas of VERA project in
National Astronomical Observatory and 25m antenna in Shanghai, and now, they were arranged for Gala-V.
The other new Wideband feed is designed for 3.2-14.4GHz, named NINJA Feed, is shown in left of
Figure 2. Aperture efficiency of 34m antenna with NINJA is shown in Figure 3. however, it will be
improved after adjusting the feed positions. They replaced in previous two C-band feeds in 34m, which
both are corrugated horns for 5GHz and 10GHz. These new wideband feeds brought us more wider
simultaneous bandwidth which no other large radio telescopes have. The Author supposes the IGUANA
feed with 2.2-18 or 22GHz coverage will be manufactured in 2015 or early 2016.
Table 1. Frequencies of Gala-V

Figure 1. Kashima 34m antenna(photo by T.Kondo).

60
50
40
30
20
10
0
3

5

7

9

11

13

15

Figure 2. Wideband feeds in

Figure 3. Aperture efficiency of 34m antenna with

Kashima 34m antenna.

NINJA feed(Aug.6,2015, measured by K.Takefuji).

Reference
1. Delopement of Mutimode Horns and Wideband Feed for Radio Telescopes, H.UJIHARA et.al.
http://www.ursi.org/proceedings/procGA11/ursi/JP2-16.pdf- 92 -J4

New X-band Receiver System of Usuda 64m Antenna
Yasuhiro Murata, Masato Tsuboi, Kentaro Yamaguchi, Kenta Uehara, Hiroshi Takeuchi, Kiyoshi Nakajima,
Zen-ichi Yamamoto, Nanako Mochizuki
Institute of Space and Astronautical Science, JAXA
1-1-3 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 252-5210, Japan
Email: murata@vsop.isas.jaxa.jp
Yusuke Kono, Masahiro Kanaguchi, Shunsaku Suzuki, Tomoaki Oyama
National Astronomy Observatory of Japan
2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
We have used the S/X-band common feed port for the spacecraft
operation, for the radio astronomy using Usuda 64m antenna (Figure1). This X-band receiving system connected to this port has higher
system noise temperature because it has some microwave circuits,
such as diplexer, filters, and waveguide switch in room temperature,
which are all required for the spacecraft, and has narrower (~ 300
MHz) band pass for radio astronomy. Because of the improvementof backend of the radio-astronomy system, more than 500 MHz
bandwidth is required for getting better sensitivity to observe
continuum sources.
We have paid attention to the unused X-band receiving-only
feed (#6 port) of Usuda 64m, and developed new receiver (Figure 2)
to attach it. The new receiver is consist of thermal insulating
waveguide, septum type low loss circular polarizer, coupler forcalibration signal insertion, and 2-low noise amplifiers (LNA’s) for
both LHCP and RHCP detection.. All of those are operated in the
dewar cooled less than 10 K. 2-LNA’s are LNF-LNC7_10A made by
Low Noise Factory, and has 36 dB gain, and 4 K receiver noise
temperature in 8.0-8.5 GHz at the physical temperature of 9 K
environment[1].
We installed this receiver at the port #6 of Usuda 64m andmeasure the system temperature and antenna efficiency though port
#6. We have got the typical system noise temperature at zenith25-30 K, and 45 % antenna efficiency, which improved the noiselevel of 2-3 times better. Furthermore, this new receiver allows us toextend the observing bandwidth for the continuum sources. Total
improvement of the sensitivity is factor 4-5, which allow us toobserve new category of sources.
Using this receiver we made VLBI observation with Tsukuba
32m (GSI), Yamaguchi 32m (Yamaguchi Univ. and NAOJ), andIbaragi 32m (Ibaragi Univ. and NAOJ) on May 5th, 2015, andsuccessfully confirmed we can use this new system for X-band
VLBI observation. This receiver allows us to do high sensitivity
VLBI observation and signal detection in X-band.
References
1. Low noise factory Website, http://www.lownoisefactory.com/- 93 -

Figure 1 Picture of Usuda 64m Antenna

Figure 2 Inside of new X band receiver.
(1) Thermal insulated wave guide, (2)Circular

Polarizer,

(3)

Low

Amplifier, LNF-LNC7_10A

Noise

J5

Source/Frequency Phase Referencing
for 22 GHz H2 O Masers and 43 and 86 GHz SiO Masersof the Red Supergiant, S Per
Yoshiharu Asaki
National Astronomical Observatory of Japan (NAOJ) Chile Observatory /Joint ALMA Observatory (JAO),
Alonso de Cordova 3107, Vitacura 763 0355, Santiago Chile
Email: yoshiharu.asaki@nao.ac.jp

We report on astrometric study results of quadruple maser lines of the red supergiant S Per with anovel VLBI astrometric analysis method of source/frequency phase referencing (SFPR) using the joint
VLBI array, KaVA, consisting of a Japanese VLBI array of the VERA and a Korean VLBI array of the
KVN. We observed microwave to millimeter wave maser lines of H2 O at 22 GHz, v = 1 and 2, J = 1 →0 SiO at 43 GHz, and v = 1, J = 2 → 1 SiO at 86 GHz of the circumstellar envelop (CSE) of S Per anda nearby calibrator, ICRF J0244+6228, as a positional calibrator with the separation angle of 4◦ . The
fringe phases of the 43 and 86 GHz SiO masers (v = 1 and 2, J = 1 → 0, and v = 1, J = 2 → 1) were
calibrated with those of the H2 O maser with the radial velocity of −50 km s−1 using multi-frequency
phase referencing (MFPR) in order to improve coherence at the higher frequencies. In addition, tocompare the relative positions between those four maser lines without astrophysical assumptions, MFPR
for the calibrator was also conducted to calibrate frequency band dependent phase errors such as phasedifferences of receiver instruments and ionospheric phase delays. The H2 O maser image of S Per and
ICRF J0244+6228 images were made with the KaVA data with good quality while the SFPR data analysis
has been made with the KVN data. Because of a trouble during the observation, only one baseline of the
KVN was available for the SFPR data analysis.
The SFPR imaging of the J = 1 → 0 SiO masers at 43 GHz was successful, and we could achieve
astrometric analysis between the 22 GHz H2 O and the 43 GHz SiO masers by a simple superposition of
the images. The 43 GHz SiO masers is located at around the center of the H2 O maser distribution along anortheast–southwest direction. At 86 GHz, the SFPR analysis was also successful, and we could confirm
the maser emission in the cross power spectrum after SFPR while it can hardly seen before SFPR due tovery rapid fringe phase fluctuations. On the other hand, a synthesis image at 86 GHz could not be madewith good quality because the cross power spectrum amplitude and phase has a rapid time variation. This
indicates that the 86 GHz SiO maser emission has more complicated structure which could not retrievewith the single KVN baseline.
References
1. Dodson, R. and Rioja, M., “Astrometric calibration of mm-VLBI using Source/Frequency Phase Referenced observations”, VLBA Scientific Memorandum, No. 31, 2009.
2. Dodson, R et al., “Astrometrically Registered Simultaneous Observations of the 22 GHz H2 O and 43
GHz SiO Masers toward R Leonis Minoris Using KVN and Source/Frequency Phase Referencing”,
The Astronomical Journal, 148, p. 97, 2014.
3. Rioja, M et al., “Verification of the Astrometric Performance of the Korean VLBI Network, using
comparative SFPR studies with the VLBA at 14/7 mm”, The Astronomical Journal, 148, p. 15, 2014.

- 94 -

J6

Hours Timescale Radio Transients Detected by the Nasu Radio Telescope Array
Takahiro Aoki
Faculty of Science and Engineering
Waseda University
3-4-1 Ohkubo, Shinjuku-ku, Tokyo 169-8555, Japan
Email: takahiro aoki@aoni.waseda.jp

We report on the discovery rate of transient radio sources with a timescale of hours. The Nasu radio
telescope array surveyed an area of 0.5 million square degree and discovered a transient object, labeled
WJN J1443+3439 [1], from an unknown origin. Its flux variability is thought to have lasted from minutes
to days (convincingly, hours). On the basis of this single detection, the sky surface density, Σ, of live
>3 Jy
−6 deg−2 at a flux density above 3 Jy and a freradio transients was estimated to be Σ1.42 GHz = 2+9−1.9 ×10
quency of 1.42 GHz [2]. Defining the event rate, R, as the sky surface density divided by the flux variabil>3 Jy−3−1 deg−2 year−1 .
ity timescale, τ , we estimated the event rate to be R1.42 GHz = 1+7−1 × 10 (τ /0.5 day)The left panel of Figure 1 shows the consistency between our result and others. The discovered
transient, WJN J1443+3439, locates within the shaded region in the right panel of Figure 1. From thatpoint of view, the transient might originate from a flare star or an active galactic nucleus.

Figure 1. Left: Log-log plot of sky surface density Σ >S vs. flux density S for radio transients. Several
survey results are plotted disregarding observation frequencies. Right: Parameter space of variable andtransient radio sources. The vertical axis indicates luminosity, SD2 , for flux density S and distance D.
The horizontal axis indicates the product of emission frequency ν and variability timescale τ .
>3 JyAssuming a power-law relation between source count and flux density, the surface density of Σ1.42 GHz =>S−3/2 deg−2 . From a statistical perspective, to verify
2 × 10−6 deg−2 yields Σ1.42
GHz = 0.3 × (S/mJy)>S+3/2 deg2
our survey result of Σ1.42
GHz with 95% confidence, a survey area needs at least 10 × (S/mJy)on the sky. If the WJN J1443+3439 shined on a timescale of minutes, it seems too difficult to detect such
events for the existing radio facilities; therefore, the Square Kilometre Array (SKA) with instantaneously
high sensitivity is needed. Our survey result will be completely verified by the SKA1-MID array under
its system baseline design in 2013 (0.83 deg2 field of view, 1.69 Jy system equivalent flux density, and
800 MHz bandwidth).
References
1. Niinuma, K., et al., The Astrophysical Journal Letters, 657, L37, 2007.
2. Aoki, T., et al., The Astrophysical Journal, 781, 10, 2014.

- 95 -

K1

Microwave Induced Oxidative Stress Mediated Toxicity on Male Fertility
Pattern and possible protective measures
Kavindra Kumar Kesari1,2, PhD
1

Department of Environmental Science, University of Eastern Finland, Kuopio, Finland
2
School of Life Science, Jaipur National University, Jaipur, India
kavindra.kesari@uef.fi, kavindra_biotech@yahoo.co.in

ABSTRACTThe present evidence on mobile phone radiation exposure is based on scientific research andpublic policy initiative to give an overview of what is known of biological effects that occur
at radiofrequency (RF)/ electromagnetic fields (EMFs) exposure. The reviews of last few
decades on health endpoints reported to be associated with RF include childhood leukaemia,
brain tumours, genotoxic effects, neurological effects and neurodegenerative diseases,
immune system deregulation, allergic and inflammatory responses, infertility and somecardiovascular effects. Most of the reports conclude a reasonable suspicion of mobile phonerisk based on clear evidence of bio-effects. Present study has been focused on bio-interaction
mechanism between mobile phone exposure and possible health effects which may be occur
due to reactive oxygen species (ROS) formation. Available data suggests that increased ROS
play an important role by enhancing the effect of radio frequency radiations which may cause
infertility as well as neurodegenerative diseases.
Several studies have reported significant (P< 0.05) causative effect of radio frequency
electromagnetic field on different organs with different parameters like; DNA strand break,apoptosis, melatonin, creatine kinase, PKC etc. But in the present scenario, very less oralmost none has found a solution against protection to these radiations. Here we are going todiscuss the possibility of therapeutic or retrieval impact of herbal plants and antioxidant
enzymes in response to RF-EMF radiations. There are several known herbal plants available
with known compounds which can act as antioxidants, anti cancerous leads. Radioprotective
potential of caffeic acid (green tea) and its extracts are rich source of polyphenolic
compounds. The polyphenols present in green tea (Camellia sinensis) have numerous healthbenefits due to their anti-inflammatory and anti-oxidant activities which may also be usefulagents against microwave exposure. Moreover, pulsed electromagnetic field (PEMF) therapy
is a non-invasive, simple technique that can be used as a scavenger agent to combat oxidative
stress. Pulsed electromagnetic field therapy provides significant protection by controlling
ROS production. Also another compound, melatonin has strong antioxidative potential
against microwave radiation induced oxidative stress mediated DNA damage in barin andtesticular cells.
The above discussed studies are an important issue in the present scenario. Because most of
the population (more than 95%) are using cell phones and also microwave appliances. The
whole study will explore the possible mechanism of microwave radiation interaction and also
protection against these radiations.
Key words: Cell phone, Apoptosis, ROS, Antioxidant- 96 -

K2

Temperature and frequency dependent onelectrical properties of swine liver
Shohei Kon, Kazuyuki Saito, and Koichi Ito
Chiba University
1-33 Yayoi-cho, Inage-ku, Chiba 263-8522, Japan
Email: s.kon@chiba-u.jp
In recent years, electromagnetic waves (EM) are indispensable for our life because they were used inmany situations, for example, telecommunications, home electrical appliance and medical applications (e.g.
hyperthermia, body-centric wireless communications). With the spread of the mobile radio terminals, a lotof studies in the ultra-high frequency (UHF) band and above have been made [1]. In particular, it has beenstudied about the effect of EM waves in the UHF for biological tissue by using of mobile phones. Therefore,
there are a lot of data about electrical properties of biological tissue in the UHF [2]. Those data can be used
for making phantoms and numerical calculations. There are some frequency bands like the high frequency(HF), and the medium frequency (MF) that have lower frequency than the UHF band. Although, thesefrequency band (range from hundreds kilohertz to tens megahertz) are used to some medical devices such as
an electro surgical knife and a hyperthermia, few studies have been reported [3]. As a result, there is notenough data about electrical properties of biological tissue in the MF and HF band, especially the data taking
into consideration about changes of electrical properties caused by temperatures and coagulations. Therefore,
it is necessary to study about electrical properties of biological tissue in the MF and HF band. I
In this study, we heated the fresh swine liver gently until coagulation and measured its electrical properties(relative permittivity and conductivity).We used the two terminal method using a LCR meter (6530B, Wayne
Kerr Electronics, London) for measurement. We investigated temperature dependent on electrical propertiesof swine liver from 0.3 MHz to 30 MHz with these data. Figure shows the result of mesurement at 0.3MHz.
Relative permittivity and conductivity greatly changed between 60 oC to 70 oC. It appears these great change ofelectrical properties between 60 oC and 70 oC are caused by the destruction of cell membrane.
5000

0.9
0.8
0.7
0.6

Relative permittivity

3000

0.5

Conductivity

2000

0.4
0.3

Conductivity [S/m]Relative permittivity

4000

0.2

1000

0.1
0

0
20

30

40

50

60

70

80

90

100

Temperature [oC]

Figure Electrical properties at 0.3 MHz
References
1. F. Sterzer,“Microwave medical devices,” IEEE Microwave Mag., vol. 3, no. 1,pp. 65-70, Mar. 2002.
2. C. Gabriel, “Compilation of the Dielectric Properties of Body Tissues at RF and Microwave Frequencies,”
Report N.AL/OE-TR- 1996-0037, Occupational and environmental health directorate, Radiofrequency
Radiation Division, Brooks Air Force Base, Texas (USA), June 1996.
3. K.R. Foster and H.P. Schwan, “Dielectric properties of selected tissues”, In: “Handbook of biologicaleffects of electromagnetic fields” edited by C. Polk and E. Postow, CRC Press, pp.27-96, 1986.

- 97 -

K3

Development of Biological Tissue Coagulating and Cutting Device
Combining Microwave and Radio Frequency Current
Sho Suzuki, Yuta Endo, Kazuyuki Sato, and Koichi Ito
Chiba University
1-33 Yayoi-cho, Inage-ku, Chiba 263-8522, Japan
Email: s.sho@chiba-u.jp
1. Introduction
Various types of surgical devices using microwave energy in order to stop bleeding from the organs andblood vessels are investigated and reported [1]. For smooth progress of the surgery, surgical devices are
required to have multi functions: tissue coagulation and tissue resection. Generally, microwave energy is used
for tissue coagulation and unsuitable for tissue resection because of the mild heating effect. Here, devices
using joule heating by discharging radio frequency (RF) current are employed for tissue resection. In thisstudy, a hybrid surgical device for tissue coagulation and resection with the microwave energy which has agood coagulation ability and radio frequency current which has a good resection ability is proposed. In thisabstract, the heating characteristic of the designed microwave heating antenna for tissue coagulation loadedwith the RF electrode for cutting device is simulated.
2. The proposed device and the calculation model
Figure 1 shows the designed heating antenna. In Fig. 1, antenna elements are connected to an inner conductor and an outer conductor of the coaxial cable. Moreover, an active electrode which emits RF current isloaded with the proposed device. In analysis, a copper plate is assumed as the active electrode. Each element
is covered with polytetrafluoroethylene (PTFE) in order to prevent adhesion of the coagulated tissue. However, the length of 2 mm from the tip of the other heating antenna connected to the outer conductor is notcovered with PTFE in order to emit RF current. Figure 2 shows the calculation model simulating the situation
when the tip of the proposed device is pressed to the surface of liver tissue. The operating frequency, input
power to the device, and heating duration were set to 2.45 GHz, 30 W, and 10 s, respectively.
3. Conclusion
Figure 3 shows a temperature distribution on the surface of liver tissue attached the proposed device (x-zplane). In Fig. 3, the parts of blood vessel near the heating antenna elements are heated more than 60 ºC.
Moreover, the reflection coefficient at the feeding point of the device is approximately -12 dB. As theseresults, it is confirmed that the designed antenna is useful as the heating antenna loaded with the RF electrode.
As further study, the proposed device will be produced and its heating characteristic will be evaluatedexperimentally.
4. Acknowledgement
This work was supported by Grant-in-Aid for Scientific Research (C) Grant Number 15K06010.
References
1. Y. Endo, K. Saito, and K. Ito, “The development of forceps-type microwave tissue coagulator for sugicaloperation” IEEE Trans. Microwave Theory Tech., vol.63, no.6 pp.2041-2049, June 2015.

Fig.1 Heating antenna.

Fig.2 Calculation model.

- 98 -

Fig.3 Temperature distribution.

K4

Study on Cytotoxicity of 0.07-0.30 THz Wave Exposure at Cultured Cells
Noriko Yaekashiwa 1, Shin'ichiro Hayashi 1, and Kodo Kawase 1, 2
1 RIKEN, RIKEN Center for Advanced Photonics, 519-1399 Aramaki-aoba, Aoba, Sendai 980-0845, Japan
Email: yaekashiwa@riken.jp
2
Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8603 Japan
Email: kawase@nuee.nagoya-u.ac.jp

We study on biological effects of exposure to THz radiation and explore the resonant vibration of cell
membrane based on the Frohlich hypothesis. Terahertz (THz) wave has characteristics of both radio waves
and visible light. This frequency will be expected early commercial use. However there are few
experimental data about safety of human body in this frequency range. So we need more investigate
regarding possible health effects.
We have carried out culture experiments to examine the cytotoxicity of 0.07-0.30 THz exposure in two
cultured human cell lines, NB1RGB (fibroblast) and A172 (glioblastoma). The Uni-Traveling Carrier
Photodiode (UTC-PD, NTT Electronics Incorporation) was used as a widely tunable THz/MMW sourcewith several micro watt output power. The distance between UTC-PD head and the cultured plate was 100
mm. The frequency was tuned within the 0.07-0.30 THz range by 1.0 GHz step by changing the difference
frequency between two pump CW fiber lasers. We have measured the cell viability assays using MTS assay
which was colorimetric method to determine the number of viable cells in proliferation or cytotoxicity
assays. This assay was based on metabolic activity to uptake and the reduction reaction of MTS reagent by
mitochondrial dehydrogenase. The absorbance was recorded at 490 nm with Microplate reader. To researcheffect of the THz exposure at cell proliferation, the cells cultured 70 hours with radiated 0.07-0.30 THz
wave. Then colorimetric reaction for MTS assay was three hours without irradiation, and absorbance wasmeasured. Another experiment, after 72 hours cultured cells without irradiation, MTS assay was carried out
three hours with irradiation for researching the cytotoxicity. Experiments were replicated three times.
The results show absorbance of irradiated cells was not decreased compared to sham cells. So we did
not find cytotoxicity of exposure at 0.07-0.30 THz in MTS assay. This suggests 0.07-0.30 THz exposure
will not cause human skin damage by ultra-weak power. However we need more detail research in this THz
range.

References
1. H. Frohlich, “Long-range coherence and energy storage in biological systems,” Int. J. Quant. Chem 2
1968, p. 641 - 649.
2. H. Ito, T. Nagatsuma, and T. Ishibashi. “Uni-Traveling-Carrier Photodiodes for High-Speed Detection
and Broadband Sensing,” Proc. of SPIE, 2007, Vol. 6479, 64790X.
3. Y. Niwano, T. Kanno, A. Iwasawa, M. Ayaki, and K. Tsubota, “Blue light injures corneal epithelial
cells in the mitotic phase in vitro,” Br J Ophthalmol, 98(7), 990-2, 2014.

- 99 -

Conference Events
Opening Ceremony

Exhibition
15:00-18:00 (September 3, 2015)
Media Hall, West Building 99:00-9:15 (September 3, 2015)
Multi-Purpose Digital Hall, West Building 9

List of Exhibitors:- AET, INC.
- MITSUBISHI ELECTRIC CORPORATION

MC: Satoshi Yagitani (Chair, URSI-JRSM 2015
Technical Program Committee)
Opening Address:
Kazuya Kobayashi (URSI-JRSM 2015 General
Chair; President, Japan National Committeeof URSI; URSI Assistant Secretary-General(AP-RASC); Vice-Chair, URSI Commission
B)
Makoto Ando (URSI Vice-President; Executive
Vice President for Research, Tokyo Instituteof Technology, Japan)

Technical Tour
10:30-16:30 (September 3, 2015)
10:30-16:30 (September 4, 2015)
Tokyo Tech Museum, Centennial Hall

Poster Session
15:30-18:00 (September 3, 2015)
Media Hall / Collaboration RoomWest Building 9
Instructions for Poster Presentation:
The board size for the poster presentation is1700 mm (height) and 1150 mm (width). The
poster size should be less than A0 size (1189 mm
height and 841 mm width). Authors are required touse only the boards corresponding to their posters.
Each poster board is marked with a poster
ID-number, which can be found in the poster
program.
Please keep in mind that the poster message
should be clear and understandable without oral
explanation. Please use an appropriate font size
allowing posters to be readable by participantsfrom 2 m away.
Posters can be put up from 12:15 on
September 3, 2015. After the poster session is over,
all the posters should be removed immediately.

Participants are free to explore Tokyo Tech
Museum in Centennial Hall on their own any time
during USRI-JRSM 2015. An exhibition room on
the 2nd floor highlights the evolution of electronics
and the development of photonics at Tokyo Tech
(Tokyo Institute of Technology), including theinvention of the high-stability quartz oscillator by
Prof. Issac Koga. During URSI-JRSM 2015, Issac
Koga Gold Medal, which is presented to an
outstanding young scientist under the age of 35 by
the Japan National Committee of URSI on theoccasion of the “URSI General Assembly and
Scientific Symposium” (URSI GASS), will also be
displayed.

- 100 -URSI-Japan
Commission Business Meetings(open to all participants)

Closing Ceremony
16:30-17:00 (September 4, 2015)
Multi-Purpose Digital Hall, West Building 9

13:00-13:45 (September 3, 2015):
Commission B: Room W935, West Building 9
Commission E: Room W936, West Building 9

MC: Jun-ichi Takada (Assistant Secretary,
Japan National Committee of URSI)

12:15-13:25 (September 4, 2015):
Commission A: Room W321, West Building 3
Commission C: Room W934, West Building 9
Commission D: Room W935, West Building 9
Commission F: Room W936, West Building 9
Commission G: Room W932, West Building 9
Commission H: Room W933, West Building 9
Commission J: Room W322, West Building 3

SPC Award Ceremony:
Winners of SPC will be awarded the 1st, 2nd
and 3rd prizes.
Closing Address:
Kazuya Kobayashi (URSI-JRSM 2015 General
Chair; President, Japan National Committeeof URSI; URSI Assistant Secretary-General(AP-RASC); Vice-Chair, URSI Commission
B)

Banquet18:30-20:30 (September 3, 2015)
Student Hall

URSI-JRSM 2015 Steering Committee
Meeting (members only)
12:15-13:15 (September 3, 2015)
Room W931, West Building 9

MC: Jun-ichi Takada (Assistant Secretary,
Japan National Committee of URSI)
Welcome Address:
Satoshi Yagitani (Chair, URSI-JRSM 2015
Technical Program Committee)
Alexander I. Nosich (Professor / Principal
Scientist, National Academy of Sciences of
Ukraine, Ukraine)

SPC Committee Meetings (members only)
SPC Committee Meeting 1:
17:25-18:10 (September 3, 2015)
Room W931, West Building 9

Toast:
Yasuhiro Koyama (Chair, URSI Commission A)
During the Banquet, three finalists for Student
Paper Competition (SPC) will be announced.

SPC Committee Meeting 2:
13:25-13:45 (September 4, 2015)
Room W931, West Building 9- 101 -

Floor Guide
Campus map, Tokyo Institute of Technology
West Bldg. 8(Smoking Area is locatedoutside this building)Okusawa Sta.

Midorigaoka Sta.
West Bldg. 9(Main Venue, andCommittee Meeting)West Bldg. 3(Committee Meetings)

Student Hall
(Banquet)

Centennial Hall(Museum)Ōokayama Sta.

2nd Floor, West Building 9
Collaboration Room

Poster session
Entrance

Cloak
Reception desk

Media Hall
Rest
Area

Exhibition & Poster
SessionRefreshment Corner

Multi-Purpose
Digital Hall

Registration
Oral Session- 102 -

3rd Floor, West Building 9
URSI-JRSM 2015 Steering Committee Meeting
SPC Committee Meetings
Commission G

W931

W936

W932

Commission H

W933

W935

W934

Commission E
Commission F

Commission C
Commission B
Commission D

3rd Floor, West Building 3

West Bldg. 2

W323

W322

Commission J

W321
Commission A

The front entrance of the building is located on the 3rd floor.
Please go down to the 2nd floor to find these meeting rooms.

- 103 -

Posters and Exhibition
Media Hall, West Building 9
Exhibitionbooth area

Media Hall
A1, B1 – B3
B4 – B6

C1 - C4

D1 – D3

Reception desk

Poster area
Coffee service

Multi-Purpose
Digital Hall

Collaboration Room, West Building 9
Collaboration Room

E1 – E4
E9, E10, F1, F2

E5 – E8
F3 – F6

F7 – F10
G3 – G6

F11, F12, G1, G2
G7, G8, H1, H2

H3 – H6
J3 – J6

H7, H8, J1, J2
K1 – K4

To the Reception desk
To the Media Hall- 104 -

List of Exhibitors- AET, INC.
- MITSUBISHI ELECTRIC CORPORATION- 105 -