Showing papers in "Earth, Planets and Space in 2018"

Geospace exploration project ERG

Abstract: The Exploration of energization and Radiation in Geospace (ERG) project explores the acceleration, transport, and loss of relativistic electrons in the radiation belts and the dynamics for geospace storms. This project consists of three research teams for satellite observation, ground-based network observation, and integrated data analysis/simulation. This synergetic approach is essential for obtaining a comprehensive understanding of the relativistic electron generation/loss processes of the radiation belts as well as geospace storms through cross-energy/cross-regional couplings, in which different plasma/particle populations and regions are strongly coupled with each other. This paper gives an overview of the ERG project and presents the initial results from the ERG (Arase) satellite.

The Plasma Wave Experiment (PWE) on board the Arase (ERG) satellite

Abstract: The Exploration of energization and Radiation in Geospace (ERG) project aims to study acceleration and loss mechanisms of relativistic electrons around the Earth. The Arase (ERG) satellite was launched on December 20, 2016, to explore in the heart of the Earth’s radiation belt. In the present paper, we introduce the specifications of the Plasma Wave Experiment (PWE) on board the Arase satellite. In the inner magnetosphere, plasma waves, such as the whistler-mode chorus, electromagnetic ion cyclotron wave, and magnetosonic wave, are expected to interact with particles over a wide energy range and contribute to high-energy particle loss and/or acceleration processes. Thermal plasma density is another key parameter because it controls the dispersion relation of plasma waves, which affects wave–particle interaction conditions and wave propagation characteristics. The DC electric field also plays an important role in controlling the global dynamics of the inner magnetosphere. The PWE, which consists of an orthogonal electric field sensor (WPT; wire probe antenna), a triaxial magnetic sensor (MSC; magnetic search coil), and receivers named electric field detector (EFD), waveform capture and onboard frequency analyzer (WFC/OFA), and high-frequency analyzer (HFA), was developed to measure the DC electric field and plasma waves in the inner magnetosphere. Using these sensors and receivers, the PWE covers a wide frequency range from DC to 10 MHz for electric fields and from a few Hz to 100 kHz for magnetic fields. We produce continuous ELF/VLF/HF range wave spectra and ELF range waveforms for 24 h each day. We also produce spectral matrices as continuous data for wave direction finding. In addition, we intermittently produce two types of waveform burst data, “chorus burst” and “EMIC burst.” We also input raw waveform data into the software-type wave–particle interaction analyzer (S-WPIA), which derives direct correlation between waves and particles. Finally, we introduce our PWE observation strategy and provide some initial results.

The ARASE (ERG) magnetic field investigation

Abstract: The fluxgate magnetometer for the Arase (ERG) spacecraft mission was built to investigate particle acceleration processes in the inner magnetosphere. Precise measurements of the field intensity and direction are essential in studying the motion of particles, the properties of waves interacting with the particles, and magnetic field variations induced by electric currents. By observing temporal field variations, we will more deeply understand magnetohydrodynamic and electromagnetic ion-cyclotron waves in the ultra-low-frequency range, which can cause production and loss of relativistic electrons and ring-current particles. The hardware and software designs of the Magnetic Field Experiment (MGF) were optimized to meet the requirements for studying these phenomena. The MGF makes measurements at a sampling rate of 256 vectors/s, and the data are averaged onboard to fit the telemetry budget. The magnetometer switches the dynamic range between ± 8000 and ± 60,000 nT, depending on the local magnetic field intensity. The experiment is calibrated by preflight tests and through analysis of in-orbit data. MGF data are edited into files with a common data file format, archived on a data server, and made available to the science community. Magnetic field observation by the MGF will significantly improve our knowledge of the growth and decay of radiation belts and ring currents, as well as the dynamics of geospace storms.

The ERG Science Center

Abstract: The Exploration of energization and Radiation in Geospace (ERG) Science Center serves as a hub of the ERG project, providing data files in a common format and developing the space physics environment data analysis software and plug-ins for data analysis. The Science Center also develops observation plans for the ERG (Arase) satellite according to the science strategy of the project. Conjugate observations with other satellites and ground-based observations are also planned. These tasks contribute to the ERG project by achieving quick analysis and well-organized conjugate ERG satellite and ground-based observations.

High Frequency Analyzer (HFA) of Plasma Wave Experiment (PWE) onboard the Arase spacecraft

Abstract: The High Frequency Analyzer (HFA) is a subsystem of the Plasma Wave Experiment onboard the Arase (ERG) spacecraft. The main purposes of the HFA include (1) determining the electron number density around the spacecraft from observations of upper hybrid resonance (UHR) waves, (2) measuring the electromagnetic field component of whistler-mode chorus in a frequency range above 20 kHz, and (3) observing radio and plasma waves excited in the storm-time magnetosphere. Two components of AC electric fields detected by Wire Probe Antenna and one component of AC magnetic fields detected by Magnetic Search Coils are fed to the HFA. By applying analog and digital signal processing in the HFA, the spectrograms of two electric fields (EE mode) or one electric field and one magnetic field (EB mode) in a frequency range from 10 kHz to 10 MHz are obtained at an interval of 8 s. For the observation of plasmapause, the HFA can also be operated in PP (plasmapause) mode, in which spectrograms of one electric field component below 1 MHz are obtained at an interval of 1 s. In the initial HFA operations from January to July, 2017, the following results are obtained: (1) UHR waves, auroral kilometric radiation (AKR), whistler-mode chorus, electrostatic electron cyclotron harmonic waves, and nonthermal terrestrial continuum radiation were observed by the HFA in geomagnetically quiet and disturbed conditions. (2) In the test operations of the polarization observations on June 10, 2017, the fundamental R-X and L-O mode AKR and the second-harmonic R-X mode AKR from different sources in the northern polar region were observed. (3) The semiautomatic UHR frequency identification by the computer and a human operator was applied to the HFA spectrograms. In the identification by the computer, we used an algorithm for narrowing down the candidates of UHR frequency by checking intensity and bandwidth. Then, the identified UHR frequency by the computer was checked and corrected if needed by the human operator. Electron number density derived from the determined UHR frequency will be useful for the investigation of the storm-time evolution of the plasmasphere and topside ionosphere.

Understanding and forecasting phreatic eruptions driven by magmatic degassing

Abstract: This paper examines phreatic eruptions which are driven by inputs of magma and magmatic gas. We synthesize data from several significant phreatic systems, including two in Costa Rica (Turrialba and Poas) which are currently highly active and hazardous. We define two endmember types of phreatic eruptions, the first (type 1) in which a deeper hydrothermal system fed by magmatic gases is sealed and produces overpressure sufficient to drive explosive eruptions, and the second (type 2) where magmatic gases are supplied via open-vent degassing to a near-surface hydrothermal system, vaporizing liquid water which drives the phreatic eruptions. The surficial source of type 2 eruptions is characteristic, while the source depth of type 1 eruptions is commonly greater. Hence, type 1 eruptions tend to be more energetic than type 2 eruptions. The first type of eruption we term “phreato-vulcanian”, and the second we term “phreato-surtseyan”. Some systems (e.g., Ruapehu, Poas) can produce both type 1 and type 2 eruptions, and all systems can undergo sealing at various timescales. We examine a number of precursory signals which appear to be important in understanding and forecasting phreatic eruptions; these include very long period events, banded tremor, and gas ratios, in particular H2S/SO2 and CO2/SO2. We propose that if these datasets are carefully integrated during a monitoring program, it may be possible to accurately forecast phreatic eruptions.

Onboard software of Plasma Wave Experiment aboard Arase: instrument management and signal processing of Waveform Capture/Onboard Frequency Analyzer

Abstract: We developed the onboard processing software for the Plasma Wave Experiment (PWE) onboard the Exploration of energization and Radiation in Geospace, Arase satellite. The PWE instrument has three receivers: Electric Field Detector, Waveform Capture/Onboard Frequency Analyzer (WFC/OFA), and the High-Frequency Analyzer. We designed a pseudo-parallel processing scheme with a time-sharing system and achieved simultaneous signal processing for each receiver. Since electric and magnetic field signals are processed by the different CPUs, we developed a synchronized observation system by using shared packets on the mission network. The OFA continuously measures the power spectra, spectral matrices, and complex spectra. The OFA obtains not only the entire ELF/VLF plasma waves’ activity but also the detailed properties (e.g., propagation direction and polarization) of the observed plasma waves. We performed simultaneous observation of electric and magnetic field data and successfully obtained clear wave properties of whistler-mode chorus waves using these data. In order to measure raw waveforms, we developed two modes for the WFC, ‘chorus burst mode’ (65,536 samples/s) and ‘EMIC burst mode’ (1024 samples/s), for the purpose of the measurement of the whistler-mode chorus waves (typically in a frequency range from several hundred Hz to several kHz) and the EMIC waves (typically in a frequency range from a few Hz to several hundred Hz), respectively. We successfully obtained the waveforms of electric and magnetic fields of whistler-mode chorus waves and ion cyclotron mode waves along the Arase’s orbit. We also designed the software-type wave–particle interaction analyzer mode. In this mode, we measure electric and magnetic field waveforms continuously and transfer them to the mission data recorder onboard the Arase satellite. We also installed an onboard signal calibration function (onboard SoftWare CALibration; SWCAL). We performed onboard electric circuit diagnostics and antenna impedance measurement of the wire-probe antennas along the orbit. We utilize the results obtained using the SWCAL function when we calibrate the spectra and waveforms obtained by the PWE.

Medium-energy particle experiments—electron analyzer (MEP-e) for the exploration of energization and radiation in geospace (ERG) mission

Abstract: The medium-energy particle experiments—electron analyzer onboard the exploration of energization and radiation in geospace spacecraft measures the energy and direction of each incoming electron in the energy range of 7–87 keV. The sensor covers a 2π-radian disklike field of view with 16 detectors, and the full solid angle coverage is achieved through the spacecraft’s spin motion. The electron energy is independently measured by both an electrostatic analyzer and avalanche photodiodes, enabling significant background reduction. We describe the technical approach, data output, and examples of initial observations.

Detection of water vapor time variations associated with heavy rain in northern Italy by geodetic and low-cost GNSS receivers

Abstract: GNSS atmospheric water vapor monitoring is not yet routinely performed in Italy, particularly at the regional scale. However, in order to support the activities of regional environmental protection agencies, there is a widespread need to improve forecasting of heavy rainfall events. Localized convective rain forecasts are often misplaced in space and/or time, causing inefficiencies in risk mitigation activities. Water vapor information can be used to improve these forecasts. In collaboration with the environmental protection agencies of the Lombardy and Piedmont regions in northern Italy, we have collected and processed GNSS and weather station datasets for two heavy rain events: one which was spatially widespread, and another which was limited to few square kilometers. The time variations in water vapor derived from a regional GNSS network with inter-station distances on the order of 50 km were analyzed, and the relationship between the time variations and the evolution of the rain events was evaluated. Results showed a signature associated with the passage of the widespread rain front over each GNSS station within the area of interest. There was a peak in the precipitable water vapor value when the heavier precipitation area surrounded the station, followed by a steep decrease (5–10 mm in about 1 h) as the rainclouds moved past the station. The smaller-scale event, a convective storm a few kilometers in extent, was not detected by the regional GNSS network, but strong fluctuations in water vapor were detected by a low-cost station located near the area of interest.

Mean winds at the cloud top of Venus obtained from two-wavelength UV imaging by Akatsuki

Abstract: Venus is covered with thick clouds. Ultraviolet (UV) images at 0.3–0.4 microns show detailed cloud features at the cloud-top level at about 70 km, which are created by an unknown UV-absorbing substance. Images acquired in this wavelength range have traditionally been used to measure winds at the cloud top. In this study, we report low-latitude winds obtained from the images taken by the UV imager, UVI, onboard the Akatsuki orbiter from December 2015 to March 2017. UVI provides images with two filters centered at 365 and 283 nm. While the 365-nm images enable continuation of traditional Venus observations, the 283-nm images visualize cloud features at an SO2 absorption band, which is novel. We used a sophisticated automated cloud-tracking method and thorough quality control to estimate winds with high precision. Horizontal winds obtained from the 283-nm images are generally similar to those from the 365-nm images, but in many cases, westward winds from the former are faster than the latter by a few m/s. From previous studies, one can argue that the 283-nm images likely reflect cloud features at higher altitude than the 365-nm images. If this is the case, the superrotation of the Venusian atmosphere generally increases with height at the cloud-top level, where it has been thought to roughly peak. The mean winds obtained from the 365-nm images exhibit local time dependence consistent with known tidal features. Mean zonal winds exhibit asymmetry with respect to the equator in the latter half of the analysis period, significantly at 365 nm and weakly at 283 nm. This contrast indicates that the relative altitude may vary with time and latitude, and so are the observed altitudes. In contrast, mean meridional winds do not exhibit much long-term variability. A previous study suggested that the geographic distribution of temporal mean zonal winds obtained from UV images from the Venus Express orbiter during 2006–2012 can be interpreted as forced by topographically induced stationary gravity waves. However, the geographic distribution of temporal mean zonal winds we obtained is not consistent with that distribution, which suggests that the distribution may not be persistent.

Difference of horizontal-to-vertical spectral ratios of observed earthquakes and microtremors and its application to S-wave velocity inversion based on the diffuse field concept

Abstract: We have been discussing the validity of using the horizontal-to-vertical spectral ratios (HVRs) as a substitute for S-wave amplifications after Nakamura first proposed the idea in 1989. So far a formula for HVRs had not been derived that fully utilized their physical characteristics until a recent proposal based on the diffuse field concept. There is another source of confusion that comes from the mixed use of HVRs from earthquake and microtremors, although their wave fields are hardly the same. In this study, we compared HVRs from observed microtremors (MHVR) and those from observed earthquake motions (EHVR) at one hundred K-NET and KiK-net stations. We found that MHVR and EHVR share similarities, especially until their first peak frequency, but have significant differences in the higher frequency range. This is because microtremors mainly consist of surface waves so that peaks associated with higher modes would not be prominent, while seismic motions mainly consist of upwardly propagating plain body waves so that higher mode resonances can be seen in high frequency. We defined here the spectral amplitude ratio between them as EMR and calculated their average. We categorize all the sites into five bins by their fundamental peak frequencies in MHVR. Once we obtained EMRs for five categories, we back-calculated EHVRs from MHVRs, which we call pseudo-EHVRs (pEHVR). We found that pEHVR is much closer to EHVR than MHVR. Then we use our inversion code to invert the one-dimensional S-wave velocity structures from EHVRs based on the diffuse field concept. We also applied the same code to pEHVRs and MHVRs for comparison. We found that pEHVRs yield velocity structures much closer to those by EHVRs than those by MHVRs. This is natural since what we have done up to here is circular except for the average operation in EMRs. Finally, we showed independent examples of data not used in the EMR calculation, where better ground structures were successfully identified from pEHVRs again. Thus we proposed here a simple empirical method to estimate S-wave velocity structures using single-station microtremor records, which is the most cost-effective method to characterize the site effects.

A comprehensive model of Earth’s magnetic field determined from 4 years of Swarm satellite observations

Abstract: The European Space Agency’s three-satellite constellation Swarm, launched in November 2013, has provided unprecedented monitoring of Earth’s magnetic field via a unique set of gradiometric and multi-satellite measurements from low Earth orbit. In order to exploit these measurements, an advanced “comprehensive inversion” (CI) algorithm has been developed to optimally separate the various major magnetic field sources in the near-Earth regime. The CI algorithm is used to determine Swarm Level-2 (L2) magnetic field data products that include the core, lithospheric, ionospheric, magnetospheric, and associated induced sources. In addition, it has become apparent that the CI is capable of extracting the magnetic signal associated with the oceanic principal lunar semidiurnal tidal constituent $$M_2$$ to such an extent that it has been added to the L2 data product line. This paper presents the parent model of the Swarm L2 CI products derived with measurements from the first 4 years of the Swarm mission and from ground observatories, denoted as “CIY4,” including the new product describing the magnetic signal of the $$M_2$$ oceanic tide.

Ultraviolet imager on Venus orbiter Akatsuki and its initial results

Abstract: The ultraviolet imager (UVI) has been developed for the Akatsuki spacecraft (Venus Climate Orbiter mission). The UVI takes ultraviolet (UV) images of the solar radiation reflected by the Venusian clouds with narrow bandpass filters centered at the 283 and 365 nm wavelengths. There are absorption bands of SO2 and unknown absorbers in these wavelength regions. The UV images provide the spatial distribution of SO2 and the unknown absorber around cloud top altitudes. The images also allow us to understand the cloud top morphologies and haze properties. Nominal sequential images with 2-h intervals are used to understand the dynamics of the Venusian atmosphere by estimating the wind vectors at the cloud top altitude, as well as the mass transportation of UV absorbers. The UVI is equipped with off-axial catadioptric optics, two bandpass filters, a diffuser installed in a filter wheel moving with a step motor, and a high sensitivity charge-coupled device with UV coating. The UVI images have spatial resolutions ranging from 200 m to 86 km at sub-spacecraft points. The UVI has been kept in good condition during the extended interplanetary cruise by carefully designed operations that have maintained its temperature maintenance and avoided solar radiation damage. The images have signal-to-noise ratios of over 100 after onboard desmear processing.

Venus looks different from day to night across wavelengths: morphology from Akatsuki multispectral images

Abstract: Since insertion into orbit on December 7, 2015, the Akatsuki orbiter has returned global images of Venus from its four imaging cameras at eleven discrete wavelengths from ultraviolet (283 and 365 nm) and near infrared (0.9–2.3 µm), to the thermal infrared (8–12 µm) from a near-equatorial orbit. The Venus Express and Pioneer Venus Orbiter missions have also monitored the planet for long periods but from polar or near-polar orbits. The wavelength coverage and views of the planet also differ for all three missions. In reflected light, the images reveal features seen near the cloud tops (~ 70 km altitude), whereas in the near-infrared images of the nightside, features seen are at mid- to lower cloud levels (~ 48–60 km altitude). The dayside cloud cover imaged at the ultraviolet wavelengths shows morphologies similar to what was observed from Mariner 10, Pioneer Venus, Galileo, Venus Express and MESSENGER. The daytime images at 0.9 and 2.02 µm also reveal some interesting features which bear similarity to the ultraviolet images. The nighttime images at 1.74, 2.26 and 2.32 µm and at 8–12 µm reveal features not seen before and show new details of the nightside including narrow wavy ribbons, curved string-like features, long-scale waves, long dark streaks, isolated bright spots, sharp boundaries and even mesoscale vortices. Some features previously seen such as circum-equatorial belts (CEBs) and occasional areal brightenings at ultraviolet (seen in Venus Express observations) of the cloud cover at ultraviolet wavelengths have not been observed thus far. Evidence for the hemispheric vortex organization of the global circulation can be seen at all wavelengths on the day- and nightsides. Akatsuki images reveal new and puzzling morphology of the complex nightside cloud cover. The cloud morphologies provide some clues to the processes occurring in the atmosphere and are thus, a key diagnostic tool when quantitative dynamical analysis is not feasible due to insufficient information.

Prediction of global ionospheric VTEC maps using an adaptive autoregressive model

Abstract: In this contribution, an adaptive autoregressive model is proposed and developed to predict global ionospheric vertical total electron content maps (VTEC). Specifically, the spherical harmonic (SH) coefficients are predicted based on the autoregressive model, and the order of the autoregressive model is determined adaptively using the F-test method. To test our method, final CODE and IGS global ionospheric map (GIM) products, as well as altimeter TEC data during low and mid-to-high solar activity period collected by JASON, are used to evaluate the precision of our forecasting products. Results indicate that the predicted products derived from the model proposed in this paper have good consistency with the final GIMs in low solar activity, where the annual mean of the root-mean-square value is approximately 1.5 TECU. However, the performance of predicted vertical TEC in periods of mid-to-high solar activity has less accuracy than that during low solar activity periods, especially in the equatorial ionization anomaly region and the Southern Hemisphere. Additionally, in comparison with forecasting products, the final IGS GIMs have the best consistency with altimeter TEC data. Future work is needed to investigate the performance of forecasting products using the proposed method in an operational environment, rather than using the SH coefficients from the final CODE products, to understand the real-time applicability of the method. The ionosphere affects national defense, aerospace, economic development, and human life. Short-term predictions of global ionospheric vertical TEC (VTEC) maps are important for the scientific analysis of the ionosphere and practical applications, such as satellite navigation. In this contribution, an adaptive autoregressive model is proposed and developed for the prediction of global ionospheric vertical TEC maps. The forecasting products are compared with final CODE and IGS GIMs, as well as altimeter TEC from JASON data under both low and mid-to-high solar activity. The results indicate that forecasting products exhibit good consistency with final GIMs.

Chronology of the 2015 eruption of Hakone volcano, Japan: geological background, mechanism of volcanic unrest and disaster mitigation measures during the crisis

Abstract: The 2015 eruption of Hakone volcano was a very small phreatic eruption, with total erupted ash estimated to be in the order of only 102 m3 and ballistic blocks reaching less than 30 m from the vent. Precursors, however, had been recognized at least 2 months before the eruption and mitigation measures were taken by the local governments well in advance. In this paper, the course of precursors, the eruption and the post-eruptive volcanic activity are reviewed, and a preliminary model for the magma-hydrothermal process that caused the unrest and eruption is proposed. Also, mitigation measures taken during the unrest and eruption are summarized and discussed. The first precursors observed were an inflation of the deep source and deep low-frequency earthquakes in early April 2015; an earthquake swarm then started in late April. On May 3, steam wells in Owakudani, the largest fumarolic area on the volcano, started to blowout. Seismicity reached its maximum in mid-May and gradually decreased; however, at 7:32 local time on June 29, a shallow open crack was formed just beneath Owakudani as inferred from sudden tilt change and InSAR analysis. The same day mud flows and/or debris flows likely started before 11:00 and ash emission began at about 12:30. The volcanic unrest and the eruption of 2015 can be interpreted as a pressure increase in the hydrothermal system, which was triggered by magma replenishment to a deep magma chamber. Such a pressure increase was also inferred from the 2001 unrest and other minor unrests of Hakone volcano during the twenty-first century. In fact, monitoring of repeated periods of unrest enabled alerting prior to the 2015 eruption. However, since open crack formation seems to occur haphazardly, eruption prediction remains impossible and evacuation in the early phase of volcanic unrest is the only way to mitigate volcanic hazard.

Low-energy particle experiments–ion mass analyzer (LEPi) onboard the ERG (Arase) satellite

Abstract: Low-energy ion experiments–ion mass analyzer (LEPi) is one of the particle instruments onboard the ERG satellite. LEPi is an ion energy-mass spectrometer which covers the range of particle energies from < 0.01 to 25 keV/q. Species of incoming ions are discriminated by a combination of electrostatic energy-per-charge analysis and the time-of-flight technique. The sensor has a planar field-of-view, which provides 4 $$\pi$$ steradian coverage by using the spin motion of the satellite. LEPi started its nominal observation after the initial checkout and commissioning phase in space.

Ground geomagnetic field and GIC response to March 17, 2015, storm

Abstract: The St. Patrick’s Day geomagnetic storm on March 17, 2015, has been chosen by the space community for synergetic analysis to build a more comprehensive picture of the storm’s origin and evolution. This storm had an unusually long (~ 17 h) main phase. During this period, many substorm-like activations occurred. These activations resulted in bursts of geomagnetically induced currents (GICs) in power lines on the Kola peninsula. To examine the substorm activations in more detail, we apply various data processing techniques for the world-wide array of magnetometers: the virtual magnetograms, magnetic latitude–local time (MLT) snapshots, and magnetic keograms. These techniques are simple tools that are supplementary to more advanced facilities developed for the analysis of SuperDARN, IMAGE, and CARISMA arrays. We compare the global spatial localization and time evolution of the geomagnetic X-component disturbance and magnetic field variability measured by the Hilbert transform of time derivative dB/dt. The latitude-MLT mapping of these magnitudes shows that very often a region with highest magnetic variability does not overlap with a substorm “epicenter” but is shifted to its poleward or equatorward boundaries. Highest variability of the geomagnetic field, and consequently intense GICs, are caused by medium-scale fast varying structures. There is no one-to-one correspondence between substorm intensity and GIC magnitude.

Mass-related excitation of polar motion: an assessment of the new RL06 GRACE gravity field models

Abstract: The new Release-06 (RL06) Gravity Recovery and Climate Experiment (GRACE) gravity field solutions are evaluated by converting them into equatorial effective angular momentum functions (so-called excitation functions) for polar motion and comparing these to respective time series based on space-geodetic observations (geodetic excitation). The same is performed for the older RL05 solutions using identical processing. Maps of equivalent water heights derived from both releases show that the signal-to-noise ratio is significantly improved in RL06. The derived polar motion excitation functions from RL05 and RL06 differ by about 15 $$\%$$ . An analysis of the contributions of different Earth subsystems revealed that the release update mainly influenced the hydrological (12 $$\%$$ ) and oceanic excitations (17 $$\%$$ ), but it has a relatively small impact on the cryospheric excitations related to Antarctica (4 $$\%$$ ) and Greenland (1 $$\%$$ ). The RL06 data from different GRACE processing centers are more consistent among each other than the previous RL05 data. Comparisons of the GRACE-based excitation functions with the geodetic and model-based oceanic excitations show that the latest release update improved the agreement by about 2 to 15 percentage points.

Relating gas ascent to eruption triggering for the April 27, 2016, White Island (Whakaari), New Zealand eruption sequence

Abstract: The April 27, 2016 eruption sequence at White Island was comprised of 6 discrete eruptive events that occurred over a 35-min period. Seismicity included three episodes of VLP activity: the first occurring ~ 2 h and a second occurring 10 min prior to the first eruption. A third larger VLP event occurred just prior to the fourth eruption. A VLP source depth of 800–1000 m below the vent is obtained from an analysis of the waveform semblance, and a volumetric source is obtained from waveform inversion of the largest VLP event. Lag times between VLP occurrence and eruption onsets provide an opportunity to examine gas migration and stress transfer models as potential triggers to the eruptive activity. Plausible lag times for a deep gas pulse to the surface are obtained by application of a TOUGH2 computational model which suggests propagation times of 0.25–1.9 m/s and are informed by previously measured White Island rock porosities and permeabilities. Results suggest that pre-eruption VLP may be plausibly linked to advection of gas from the VLP source at a magmatic carapace located ~ 800–1000 m depth. Alternatively, the large VLP that occurred just prior to the fourth eruption may be linked to a quasi-dynamic or quasi-static stress perturbation.

Magnetic Search Coil (MSC) of Plasma Wave Experiment (PWE) aboard the Arase (ERG) satellite

Abstract: This paper presents detailed performance values of the Magnetic Search Coil (MSC) that is part of the Plasma Wave Experiment on board the Arase (ERG) satellite. The MSC consists of a three-axis search coil magnetometer with a 200-mm-long magnetic core. The MSC plays a central role in the magnetic field observations, particularly for whistler mode chorus and hiss waves in a few kHz frequency range, which may cause local acceleration and/or rapid loss of radiation belt electrons. Accordingly, the MSC was carefully designed and developed to operate well in harsh radiation environments. To ascertain the wave-normal vectors, polarizations, and refractive indices of the plasma waves in a wide frequency band, the output signals detected by the MSC are fed into the two different wave receivers: one is the WaveForm Capture/Onboard Frequency Analyzer for waveform and spectrum observations in the frequency range from a few Hz up to 20 kHz, and the other is the High Frequency Analyzer for spectrum observations in the frequency range from 10 to 100 kHz. The noise equivalent magnetic induction of the MSC is $$20\,{\hbox {fT/Hz}}^{1/2}$$ at a frequency of 2 kHz, and the null depth of directionality is − 40 dB, which is equivalent to an angular error less than $$1^{\circ }$$ . The MSC on board the Arase satellite is the first experiment using a current-sensitive preamplifier for probing the plasma waves in the radiation belts.

Solar flare effects on D-region ionosphere using VLF measurements during low- and high-solar activity phases of solar cycle 24

Abstract: The effects of solar flares on the propagation of subionospheric VLF signals from NWC and NLK transmitter stations monitored at a low-latitude station, Suva (18.2°S, 178.4°E), Fiji, between December 2006 and December 2010 (an unprecedented solar minimum of solar cycles 23 and 24) and between January 2012 and December 2013 (moderate solar activity at the peak of solar cycle 24) have been analyzed to find solar flare time D-region changes. The amplitude and phase enhancements associated with solar flares were observed in the signals from both stations which are due to an increase in the electron density of the D-region as a result of extra ionization caused by the solar flares. The solar flare-induced perturbations in both the amplitude and phase of VLF signals were used to determine D-region ionospheric parameters: H′ (the ionospheric reflection height) and β (rate of increase in electron density with height) using Long Wave Propagation Capability (LWPC) version 2.1 modeling. A comparative analysis of the ionospheric D-region parameter changes carried out for this location shows a greater increase in β and decrease in H′ during low-solar activity period than during moderate-solar activity period, for the same class of flares. Our results also show greater differences in the values of β and H′ for strong flares in comparison with weak flares under both low- and moderate-solar activity conditions.

Polar motion prediction using the combination of SSA and Copula-based analysis

Abstract: The real-time estimation of polar motion (PM) is needed for the navigation of Earth satellite and interplanetary spacecraft. However, it is impossible to have real-time information due to the complexity of the measurement model and data processing. Various prediction methods have been developed. However, the accuracy of PM prediction is still not satisfactory even for a few days in the future. Therefore, new techniques or a combination of the existing methods need to be investigated for improving the accuracy of the predicted PM. There is a well-introduced method called Copula, and we want to combine it with singular spectrum analysis (SSA) method for PM prediction. In this study, first, we model the predominant trend of PM time series using SSA. Then, the difference between PM time series and its SSA estimation is modeled using Copula-based analysis. Multiple sets of PM predictions which range between 1 and 365 days have been performed based on an IERS 08 C04 time series to assess the capability of our hybrid model. Our results illustrate that the proposed method can efficiently predict PM. The improvement in PM prediction accuracy up to 365 days in the future is found to be around 40% on average and up to 65 and 46% in terms of success rate for the $${\hbox{PM}}_{x}$$ and $${\hbox{PM}}_{y}$$ , respectively.

High-energy electron experiments (HEP) aboard the ERG (Arase) satellite

Abstract: This paper reports the design, calibration, and operation of high-energy electron experiments (HEP) aboard the exploration of energization and radiation in geospace (ERG) satellite. HEP detects 70 keV–2 MeV electrons and generates a three-dimensional velocity distribution for these electrons in every period of the satellite’s rotation. Electrons are detected by two instruments, namely HEP-L and HEP-H, which differ in their geometric factor (G-factor) and range of energies they detect. HEP-L detects 70 keV–1 MeV electrons and its G-factor is 9.3 × 10−4 cm2 sr at maximum, while HEP-H observes 0.7–2 MeV electrons and its G-factor is 9.3 × 10−3 cm2 sr at maximum. The instruments utilize silicon strip detectors and application-specific integrated circuits to readout the incident charge signal from each strip. Before the launch, we calibrated the detectors by measuring the energy spectra of all strips using γ-ray sources. To evaluate the overall performance of the HEP instruments, we measured the energy spectra and angular responses with electron beams. After HEP was first put into operation, on February 2, 2017, it was demonstrated that the instruments performed normally. HEP began its exploratory observations with regard to energization and radiation in geospace in late March 2017. The initial results of the in-orbit observations are introduced briefly in this paper.

Fault model of the 2017 Jiuzhaigou Mw 6.5 earthquake estimated from coseismic deformation observed using Global Positioning System and Interferometric Synthetic Aperture Radar data

Abstract: On August 8, 2017, the Jiuzhaigou Mw 6.5 earthquake occurred in Sichuan province, southwestern China, along the eastern margin of the Tibetan Plateau. The epicenter is surrounded by the Minjiang, Huya, and Tazang Faults. As the seismic activity and tectonics are very complicated, there is controversy regarding the accurate location of the epicenter and the seismic fault of the Jiuzhaigou earthquake. To investigate these aspects, first, the coseismic deformation field was derived from Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) measurements. Second, the fault geometry, coseismic slip model, and Coulomb stress changes around the seismic region were calculated using a homogeneous elastic half-space model. The coseismic deformation field derived from InSAR measurements shows that this event was mainly dominated by a left-lateral strike-slip fault. The maximal and minimal displacements were approximately 0.15 m and − 0.21 m, respectively, along line-of-sight observation. The whole deformation field follows a northwest-trending direction and is mainly concentrated west of the fault. The coseismic slip is 28 km along the strike and 18 km along the dip. It is dominated by a left-lateral strike-slip fault. The average and maximal fault slip is 0.18 and 0.85 m, respectively. The rupture did not fully reach the ground surface. The focal mechanism derived from GPS and InSAR data is consistent with the kinematics and geometry of the Huya Fault. Therefore, we conclude that the northern section or the Shuzheng segment of the Huya Fault is the seismogenic fault. The maximal fault slip is located at 33.25°N and 103.82°E at a depth of ~ 11 km, and the release moment is approximately 6.635 × 1018 Nm, corresponding to a magnitude of Mw 6.49, which is consistent with results reported by the US Geological Survey, Global Centroid Moment Tensor, and other researchers. The coseismic Coulomb stress changes enhanced the stress on the northwest and southeast edges of the northern extension of the Huya Fault. Seismic risks cannot be ignored in the future although aftershocks are fewer in number in these regions.

Water sampling using a drone at Yugama crater lake, Kusatsu-Shirane volcano, Japan

Abstract: Remote sampling of water from Yugama crater lake at Kusatsu-Shirane volcano, Japan, was performed using a drone. Despite the high altitude of over 2000 m above sea level, our simple method was successful in retrieving a 250 mL sample of lake water. The procedure presented here is easy for any researcher to follow who operates a drone without additional special apparatus. We compare the lake water sampled by drone with that sampled by hand at a site where regular samplings have previously been carried out. Chemical concentrations and stable isotope ratios are largely consistent between the two techniques. As the drone can fly automatically with the aid of navigation by Global Navigation Satellite System (GNSS), it is possible to repeatedly sample lake water from the same location, even when entry to Yugama crater lake is restricted due to the risk of eruption.

InSAR analysis for detecting the route of hydrothermal fluid to the surface during the 2015 phreatic eruption of Hakone Volcano, Japan

Abstract: Although the 2015 Hakone Volcano eruption was a small-scale phreatic eruption with a discharged mass of only about 100 tons, interferometric synthetic aperture radar successfully detected surface deformations related to the eruption. Inversion model of the underground hydrothermal system based on measured ground displacements by ALOS-2/PALSAR-2 images showed that a crack opened at an elevation of about 530–830 m, probably at the time of the eruption. A geomorphological analysis detected several old NW–SE trending fissures, and the open crack was located just beneath one of the fissures. Thus, the crack that opened during the 2015 eruption could have been a preexisting crack that formed during a more voluminous hydrothermal eruption. In addition, the inversion model implies that a sill deflation occurred at an elevation of about 225 m, probably at the time of the eruption. The deflation of sill-like body represents a preexisting hydrothermal reservoir at an elevation of 100–400 m, which intruded fluid in the open crack prior to eruption. The volume changes of the open crack and the sill were calculated to be 1.14 × 105 m3 (inflation) and 0.49 × 105 m3 (deflation), respectively. A very local swelling (about 200 m in diameter) was also detected at the eruption center 2 months before the eruption. The local swelling, whose rate in satellite line-of-sight was 0.7–0.9 cm/day during May 2015 and declined in June, had been monitored until the time of the eruption, when its uplift halted. This was modeled as a point pressure source at an elevation of about 900 m (at a depth of about 80–90 m from the ground surface) and is considered to be a minor hydrothermal reservoir just beneath the fumarolic field. Our analysis shows that the northernmost tip of the open crack reached within 200 m of the surface. Thus, it is reasonable to assume that the hydrothermal fluid in the open crack found a way to the surface and formed the eruption.

Self-potential mapping using an autonomous underwater vehicle for the Sunrise deposit, Izu-Ogasawara arc, southern Japan

Abstract: We performed a simultaneous survey of self-potential and plume turbidity using an autonomous underwater vehicle (AUV) above the Sunrise deposit in the Myojin Knoll caldera of the Izu-Ogasawara arc. A 10-m-long electrode rod, on which five electrodes referenced with a common electrode were mounted, was connected at the tail of an AUV. The survey was conducted at a typical speed of 2 knots, covering the 1500 m × 1500 m area with a typical spacing of survey lines of 100 m. With AUV altitude of 100 m above the seafloor, a negative self-potential anomaly of a few millivolts was observed. The self-potential anomaly was found to spread 300 m × 300 m. The self-potential is probably attributable to the geo-battery mechanism: electric current is generated by redox reactions occurring around an ore body crossing a redox contrast. Assuming that the source of the self-potential is an electric current dipole, we can image a southward-dipping dipole with the moment of approximately 103 A m, approx. 30 m below the southern part of the ore deposit. Anomalies of turbidity, which are correlated to ambient temperature and which are signatures of discharged hydrothermal fluids, were distributed more broadly than the self-potential. Some turbidity anomalies were found without self-potential anomalies. They were probably transported by the ocean current. Spatial decoupling between the self-potential and turbidity anomalies suggests that the direct contribution of hydrothermal fluids to the self-potential anomalies is probably a secondary effect. The survey altitude of 100 m and the survey speed of 2 knots in the present study represent practical limitations for the self-potential survey when active hydrothermal fields are targeted. We have observed that the self-potential method responds exclusively to the presence of hydrothermal ore deposits. This behavior differs from other methods for exploring seafloor hydrothermal ore deposits: The geomagnetic method responds not only to ore deposits but also to volcanic bodies. The plume method can detect remote hydrothermal activities, but the source locations are not necessarily specified. The self-potential method is useful as an excellent exploration tool, particularly for initial surveys.

The extremely high-energy electron experiment (XEP) onboard the Arase (ERG) satellite

Abstract: The extremely high-energy electron experiment (XEP) onboard the Arase (ERG) satellite is designed to measure high-energy electrons in the Earth’s radiation belts. The XEP was developed by taking advantage of our technical heritage of high-energy particle detectors that are onboard Earth observation satellites of the Japan Aerospace Exploration Agency (JAXA) as the radiation monitor. The main target of the XEP is to precisely measure variations of relativistic electrons in the outer radiation belt even during magnetic storms. The measurement is scientifically required to address physical mechanisms of electron acceleration and loss. The XEP consists of five solid-state silicon detectors (SSDs) and a single-crystal inorganic scintillator of cerium-doped gadolinium orthosilicate (GSO) to measure electrons in the energy range of 0.4–20 MeV and has a 20° single field of view (FOV). It is also equipped with a plastic scintillator that surrounds the GSO scintillator to prevent particles from entering the detectors from outside the FOV. The XEP has started its observation of relativistic electrons and has successfully observed dynamic variations of relativistic electron fluxes in the outer radiation belt during magnetic storms. This paper describes the instrumentation of the XEP and presents an example of initial observation results.

Resistivity characterisation of Hakone volcano, Central Japan, by three-dimensional magnetotelluric inversion

Abstract: On 29 June 2015, a small phreatic eruption occurred at Hakone volcano, Central Japan, forming several vents in the Owakudani geothermal area on the northern slope of the central cones. Intense earthquake swarm activity and geodetic signals corresponding to the 2015 eruption were also observed within the Hakone caldera. To complement these observations and to characterise the shallow resistivity structure of Hakone caldera, we carried out a three-dimensional inversion of magnetotelluric measurement data acquired at 64 sites across the region. We utilised an unstructured tetrahedral mesh for the inversion code of the edge-based finite element method to account for the steep topography of the region during the inversion process. The main features of the best-fit three-dimensional model are a bell-shaped conductor, the bottom of which shows good agreement with the upper limit of seismicity, beneath the central cones and the Owakudani geothermal area, and several buried bowl-shaped conductive zones beneath the Gora and Kojiri areas. We infer that the main bell-shaped conductor represents a hydrothermally altered zone that acts as a cap or seal to resist the upwelling of volcanic fluids. Enhanced volcanic activity may cause volcanic fluids to pass through the resistive body surrounded by the altered zone and thus promote brittle failure within the resistive body. The overlapping locations of the bowl-shaped conductors, the buried caldera structures and the presence of sodium-chloride-rich hot springs indicate that the conductors represent porous media saturated by high-salinity hot spring waters. The linear clusters of earthquake swarms beneath the Kojiri area may indicate several weak zones that formed due to these structural contrasts.