Groundmass pargasite in the 1991–1995 dacite of Unzen volcano: phase stability experiments and volcanological implications

Pargasite commonly occurs in the dacitic groundmass of the 1991–1995 eruption products of Unzen volcano. We described the occurrence and chemical compositions of amphibole in the dacite, and also carried out melting experiments to determine the low-pressure stability limit of amphibole in the dacite. The 1991–1995 ejecta of the Unzen volcano show petrographic evidence of magma mixing, such as reverse compositional zoning of plagioclase and amphibole phenocrysts, and we used a groundmass separate as a starting material for the experiments. Reversed experiments show that the maximum temperature for the crystallization of amphibole is 930°C at 196 MPa, 900°C at 98 MPa, and 820°C at 49 MPa. Compared with the experimental results on the Mount St. Helens dacite, present experiments on the Unzen dacitic groundmass show that amphibole is stable to pressures ca. 50 MPa lower at 850°C. Available Fe–Ti oxide thermometry indicates the crystallization temperature of the groundmass of the Unzen dacite to be 880±30°C, suggesting that the groundmass pargasite crystallized at >70 MPa, corresponding to a depth of more than 3 km in the conduit. The chlorine content of the groundmass pargasite is much lower than that of phenocrystic magnesiohornblende in the 1991–1995 dacite of Unzen volcano, indicating that vesiculation/degassing of magma took place before the crystallization of the groundmass pargasite. The present study shows that the magma was water oversaturated and that the degassing of magma along with magma mixing caused crystallization of the groundmass amphibole at depths of more than 3 km in the conduit.     

Mg-Fe partitioning between olivine and ultramafic melts at high pressures

Precise determination of the partitioning of Mg and Fe²⁺ between olivine and ultramafic melt has been made at pressures from 5 to 13 GPa using a MA-8 type multi-anvil high-pressure apparatus (PREM) installed at Earthquake Research Institute, University of Tokyo. A very short rhenium capsule (<100 μm sample thickness) was adopted to minimize temperature variation within the sample container. Synthetic gels with the composition of the upper mantle peridotite were used as starting materials to promote the homogeneity. Analyses of quenched melts and coexisting olivines were made with an electron probe microanalyzer. The obtained partition coefficient, K_D [=(FeO/MgO)^ol/(FeO/MgO)^melt], decreases from 0.35 to 0.25 with increasing pressure from 5 to 13 GPa, suggesting a negative correlation between pressure and K_D above 5 GPa. Our result is consistent with a parabolic relationship between K_D and degree of polymerization (NBO/T) of melts reported by previous studies at lower pressures. The negative correlation between pressure and K_D suggests that olivine crystallizing in a magma ocean becomes more Mg-rich with depth and that primary magmas generated in the upper mantle become more Fe-rich with depth than previously estimated.     

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.     

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     

Seasonal barotropic sea surface height fluctuation in relation to regional ocean mass variation

A consistency between seasonal fluctuation of actual sea surface height (SSH) and those caused by mass and density variations in gyre-scale regions is examined. The SSH obtained from satellite altimetry (altimetric SSH) is adopted as the actual SSH. SSH caused by mass variation (mass-related SSH) is simulated using a barotropic global ocean model forced by water flux, wind stress and surface pressure. SSH caused by density variation (steric SSH) is calculated from water density profile, i.e. temperature and salinity profiles. The model SSH well represents mass-related SSH for gyre-scale regional means, and seasonal fluctuation of the altimetric SSH corrected for the model SSH is similar to that of steric SSH above a pressure level larger than 300 dbar. The results indicate that the mass-related SSH does not much respond to the baroclinic adjustment to the seasonally varying wind stress curl. The mass-related SSH forced by wind stress and surface pressure should be accounted for regional evaluation, though it is not necessary for global mean evaluation. Detection of steric SSH from altimetric SSH would be useful for assimilation approaches in which the altimetric SSH is treated as the variable reflecting subsurface temperature and salinity.     

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.     

Slip Distribution and Seismic Moment of the 2010 and 1960 Chilean Earthquakes Inferred from Tsunami Waveforms and Coastal Geodetic Data

The slip distribution and seismic moment of the 2010 and 1960 Chilean earthquakes were estimated from tsunami and coastal geodetic data. These two earthquakes generated transoceanic tsunamis, and the waveforms were recorded around the Pacific Ocean. In addition, coseismic coastal uplift and subsidence were measured around the source areas. For the 27 February 2010 Maule earthquake, inversion of the tsunami waveforms recorded at nearby coastal tide gauge and Deep Ocean Assessment and Reporting of Tsunamis (DART) stations combined with coastal geodetic data suggest two asperities: a northern one beneath the coast of Constitucion and a southern one around the Arauco Peninsula. The total fault length is approximately 400 km with seismic moment of 1.7 × 10²² Nm (Mw 8.8). The offshore DART tsunami waveforms require fault slips beneath the coasts, but the exact locations are better estimated by coastal geodetic data. The 22 May 1960 earthquake produced very large, ~30 m, slip off Valdivia. Joint inversion of tsunami waveforms, at tide gauge stations in South America, with coastal geodetic and leveling data shows total fault length of ~800 km and seismic moment of 7.2 × 10²² Nm (Mw 9.2). The seismic moment estimated from tsunami or joint inversion is similar to previous estimates from geodetic data, but much smaller than the results from seismic data analysis.     

Alteration processes recorded by back-arc mantle peridotites from oceanic core complexes,Shikoku Basin,Philippine Sea

We determined the mineralogical and petrological characteristics of ultramafic rocks dredged from two oceanic core complexes: the Mado Megamullion and 23°30′N non-transform offset massif, which are located within the Shikoku back-arc basin in the Philippine Sea. The ultramafic rocks are strongly serpentinized, but can be classified as harzburgite/lherzolite or dunite, based on relict primary minerals and their pseudomorphs. Strongly elongated pyroxene porphyroclasts with undulatory extinction indicate high-temperature (≥700 °C) strain localization on a detachment fault within the upper mantle at depths below the brittle–viscous transition. During exhumation, the peridotites underwent impregnation by magmatic or hydrothermal fluids, lizardite/chrysotile serpentinization at ≤300 °C, antigorite crystallization, and silica metasomatism that formed talc. These features indicate that the detachment fault zones formed a fluid pathway and facilitated a range of fluid–peridotite interactions.     

Phase-lagged warming and the disruption of climatic rhythms during the Matuyama–Brunhes magnetic polarity transition

The influence of the geomagnetic field on climate, a long hotly disputed issue, was examined for the Matuyama–Brunhes (MB) magnetic polarity transition using palynological and paleomagnetic data sets from a high accumulation rate (ca. 50 cm/kyr) sediment core from Osaka Bay. During the period from marine oxygen isotope stages (MIS) 20 to 18, climate change is well correlated with global ice volume variation in the precession cycle, with the exception of the early half of MIS 19. The postglacial warming after substage 20.2 was interrupted by cooling that began just before the sea-level highstand correlated with substage 19.3 and persisted until about a mid-point between 19.3 and 19.2, followed by a rapid warming. The thermal maximum clearly postdates the highest sea-level highstand by 6–7 kyr, and the connection between orbital forcing and climate was disrupted. The cooling event coincided with the center of the paleointensity low during the MB transition. This unusual climate cooling across a sea-level peak is very likely related to the field intensity decrease. The data from Osaka Bay may suggest an instance where the geomagnetic field has influenced climate in the past.     

Availability of an equation to evaluate error by optical path discretization in radiative transfer computation based on the successive order of scattering method

The availability of an equation to evaluate the influence of multiple scattering in the single scattering process corresponding to a layer of arbitrary optical thickness was established. In order to confirm the validity of this equation, the radiance distribution in this layer was computed using a plane–parallel layer model based on the successive order of scattering method. The relative errors in a radiance distribution computed were evaluated as a function of optical thickness by the derived equation. It was shown that this equation provides a theoretical background for determining layer thickness using the plane–parallel layer model.