A heating process of Kuchi-erabu-jima volcano,Japan, as inferred from geomagnetic field variations and electrical structure

Since August 2000, we have recorded the total intensity of the geomagnetic field at the summit area of Kuchi-erabu-jima volcano, where phreatic eruptions have repeatedly occurred. A time series analysis has shown that the variations in the geomagnetic field since 2001 have a strong relationship to an increase in volcanic activity. These variations indicate thermal demagnetization of the subsurface around the presently active crater. The demagnetization source for the early variations, until summer 2002, was estimated at about 200 m below sea level. For the variations since 2003, the source was modeled on the basis of the expansion of a uniformly magnetized ellipsoid. The modeling result showed that the source is located at 300 m above sea level beneath the crater. We carried out an audio-frequency magnetotelluric survey with the aim of obtaining a relation between the demagnetization source and the shallow structure of the volcano. A two-dimensional inversion applied to the data detected two good conductors, a shallow thin one which is restricted to a region around the summit area, while the other extends over the edifice at depths between 200 and 800 m. These conductors are regarded as clay-rich layers with low permeability, which were assumed to be generated through hydrothermal alteration. The demagnetization source for the early variations was possibly located at the lower part of the deep conductor and the source after 2003 lies between the two conductors, where groundwater is considered to be abundant. Based on these results, as well as on seismological, geodetic, and geochemical information, we propose a heating process of the Kuchi-erabu-jima volcano. In the initial stage, high-temperature volcanic gases supplied from the deep-seated magma remained temporarily at the level around the lower part of the less permeable deep conductor since the ascent path had not yet been established. Then, when the pathway developed as a result of repeated earthquakes, it became possible for a massive flux of volcanic gases to ascend through the conductor. The high temperature gases reached the aquifer located above the conductor and the heat was efficiently transported to the surrounding rocks through the groundwater. As a consequence, an abrupt increase of the gas flux and diffusion of the heat through the aquifer occurred and the high-temperature zone expanded. Since the high-temperature zone is located beneath another conductor, which acts as caprock, we assume that the energy of the phreatic explosion is accumulated there.     

Changes in the lower limit of mountain permafrost between 1973 and 2004 in the Khumbu Himal, the Nepal Himalayas

Because the Khumbu Himal of the Nepal Himalayas lacks long-term climate records from weather stations, mountain permafrost degradation serves as an important indicator of climate warming. In 1973, the permafrost lower limit was estimated to be 5200–5300 m above sea level (ASL) on southern-aspect slopes in this region. Using ground-temperature measurements, we examined the mountain permafrost lower limit on slopes with the same aspect in 2004. The results indicate that the permafrost lower limit was 5400–5500 m ASL in 2004. The permafrost lower limit was estimated to be 5400 to 5500 m on slopes with a southern aspect in the Khumbu Himal in 1991 using seismic reflection soundings. Thus, it is possible that the permafrost lower limit has risen 100–300 m between 1973 and 1991, followed by a stable limit of 5400 to 5500 m over the last decade. An increase in mean annual air temperature of approximately 0.2 to 0.4 °C from the 1970s to the 1990s has indicated a rise in the permafrost lower limit of 40 to 80 m at the Tibetan Plateau. The rise in the mountain permafrost lower limit in the Khumbu Himal exceeds that of the Tibetan Plateau, suggesting the possibility of greater climate warming in the Khumbu Himal.     

Time variations of hard X-rays from Sco X-1

Simultaneous hard X-ray and optical observations of Sco X-1 were carried out on 1971 May 1 at Hyderabad, India, when Sco X-1 was optically bright. The X-ray intensity observed by balloon-borne counter telescopes increased in coincidence with optical enhancements, while the plasma temperature derived by fitting the X-ray spectrum in the energy range 20–40 keV to the thermal bremsstrahlung spectrum did not appreciably change over the whole period of observation.     

Nuclear field shift effect as a possible cause of Te isotopic anomalies in the early solar system—An alternative explanation of Fehr et al. (2006 and 2009)

Abstract— We explore the possibility that Te isotopic anomalies measured in Ca‐Al‐rich inclusions (Fehr et al. 2009) and in leachates of carbonaceous chondrites (Fehr et al. 2006) may be due to mass‐independent effects controlled by nuclear field shift rather than to nucleosynthetic processes. Fehr et al.'s spectrum of mass‐independent anomalies of Te isotopes shows a smooth correlation with mass number and nuclear charge distribution. Ratios of even to odd isotopes, as the ¹²⁵Te/¹²⁶Te ratio used by these authors for normalization are particularly prone to nuclear field shift effects. We show that the alternative normalization of isotopic ratios to ¹³⁰Te/¹²⁶Te strongly reduces the trend of isotopic fractionation with mass number, leaving only ¹²⁵Te as truly anomalous. For both normalizations (¹²⁵Te/¹²⁶Te and ¹³⁰Te/¹²⁶Te), Fehr et al.'s results fit the theory of Bigeleisen (1996), which suggests that the nuclear field shift effect can potentially account for the observed Te isotope abundances, as an alternative to nucleosynthetic processes. We propose that these mass‐independent effects may be acquired during accretion of sulfides from the solar nebula.     

Compaction and fragmentation of porous gypsum targets from low-velocity impacts

We performed low-velocity impact experiments of gypsum spheres with porosity ranging from 0 to 61% and diameter ranging from 25 to 83 mm. The impact velocity was from 0.2 to 22 m/s. The target was an iron plate. The outcome of gypsum spheres with porosity 31-61% was different from those of non-porous ice [Higa M., Arakawa, M., Maeno, N., 1996. Planet. Space Sci. 44, 917-925; Higa M., Arakawa, M., Maeno, N., 1998. Icarus 133, 310-320] and non-porous gypsum. In between the intact and fragmentation modes, the outcome of the non-porous ice and gypsum was crack growth at the impact point. However, the outcome of the porous gypsum was compaction. We found that the restitution coefficients of the porous gypsum spheres were all in a similar range, in spite of the difference of the porosity and size at impact velocities up to about 10 m/s where they begin to be fragmented in pieces. Moreover, there is not a large difference between the restitution coefficient of porous and non-porous gypsum. These results collectively indicate that restitution coefficient of gypsum spheres of cm-size is not strongly dependent upon the porosity and compaction process.     

Ion upflow and downflow at the topside ionosphere observed by the EISCAT VHF radar

We have determined the MLT distribution and KpK_P dependence of the ion upflow and downflow of the thermal bulk oxygen ion population based on a data analysis using the EISCAT VHF radar CP-7 data obtained at Tromsø during the period between 1990 and 1996: (1) both ion upflow and downflow events can be observed at any local time (MLT), irrespective of dayside and nightside, and under any magnetic disturbance level, irrespective of quiet and disturbed levels; (2) these upflow and downflow events are more frequently observed in the nightside than in the dayside; (3) the upflow events are more frequently observed than the downflow events at any local time except midnight and at any K_P level and the difference of the occurrence frequencies between the upflow and downflow events is smaller around midnight; and (4) the occurrence frequencies of both the ion upflow and downflow events appear to increase with increasing KP level, while the occurrence frequency of the downflow appears to stop increasing at some KP level     

Radiocarbon dates of elevated shorelines,Indonesia and Malaysia. Part 1

Four new radiocarbon dates of elevated strandlines in tectonically active areas of eastern Indonesia and East Malaysia indicate average rates of uplift that range between 4.5 and 9 mm annually during the past 24,000 yr. These values are at least three times higher than former estimates from eastern Indonesia. Another radiocarbon date from the south arm of Sulawesi—also tectonically mobile—indicates a rate of uplift of 1.4–2.5 mm per year which corresponds with earlier determinations. This particular case, however, suggests that the sample was probably located close to a north-south axis about which southern Sulawesi was tilted during the Quaternary.     

Propagating tsunami wave and subsequent resonant response signals detected by HF radar in the Kii Channel,Japan

Signals from the tsunami waves induced by the March 11, 2011 moment magnitude (M_w) 9.0 Tohoku-Oki earthquake and from subsequent resonances were detected as radial velocity variability by a high-frequency ocean surface radar (HF radar) installed on the eastern coast of the Kii Channel, at a range of about 1000 km from the epicenter along the eastern to southern coasts of Honshu Island. A time–distance diagram of band-passed (9–200 min) radial velocity along the beam reveals that the tsunami waves propagated from the continental shelf slope to the inner channel as progressive waves for the first three waves, and then natural oscillations were excited by the waves; and that the direction of the tsunami wave propagation and the axis of the natural oscillations differed from that of the radar beam. In addition, spectral analyses of the radial velocities and sea surface heights obtained in the channel and on the continental shelf slope suggest complex natural oscillation modes excited by the tsunami waves.     

Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture

The hybrid finite-discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general-purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU-parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi-static and dynamic loading conditions are conducted using the GPGPU-parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed-up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU-parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU-parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.     

Estimating the impact parameters of cosmic dust particles using a piezoelectric lead zirconate titanate detector

A cosmic dust detector for installation on a satellite is currently being developed using piezoelectric lead zirconate titanate (PZT), which can possess both functions of the collector and the transducer. The characteristics of the PZT detector have been studied by bombarding it with hypervelocity particles supplied by a Van de Graaff accelerator. The front surface of the detector used in this study was covered with a white paint to reduce any increase in the temperature due to the solar radiation. There was a linear relationship between the rise time of the signal produced by the detector and the particle's velocities, which were above 10 km/s on impact. This implies that individual particle velocities on impact can be inferred through the empirical formula derived from the data obtained from the PZT detector.