Cyclic magma storages and transfers at Piton de La Fournaise volcano (La Réunion hotspot) inferred from deformation and geochemical data     

Aline Peltier  Vincent Famin  Patrick Bachlery  Valrie Cayol  Yo Fukushima  Thomas Staudacher 《Earth and Planetary Science Letters》2008,270(3-4):180-188

Piton de La Fournaise is in a period of intense volcanic activity since 1998. To constrain the magma dynamics responsible for this activity, we combined GPS ground deformation monitoring interpreted through numerical modelling and geochemistry. Two cycles of continuous volcano inflation are evidenced for the May 2004–December 2005 period, with a rest from March to October 2005. These inflations are consistent with two cycles of compatible major element enrichment in the emitted lavas. Numerical models indicate that the pressurization of a single magma reservoir may be responsible for the observed pre-eruptive inflations of the volcano. The reservoir, located at 2300 m depth, has a radius of  500 m. At the beginning of each cycle, dykes propagate from the roof of the reservoir and yield eruptions of differentiated basalt near the summit. At the end of the cycle, dykes propagate from the eastern sidewall of the reservoir and yield distal eruptions of primitive magmas away from the summit. The volumes of magma emitted during the primitive eruptions seem too large to explain the surface deformations and therefore suggest some refill of the reservoir by deeper magmas. Our results may be used to predict the location and lava volume of future eruptions at Piton de La Fournaise volcano, depending on the timing of these eruptions within a cycle of volcanic activity.  相似文献   

Characteristics of isolated hybrid tremor (HBT) during a calm activity period at Aso Volcano     

Takehiko Mori  Yasuaki Sudo  Tomoki Tsutsui  Shin Yoshikawa 《Bulletin of Volcanology》2008,70(9):1031-1042

Isolated-type tremors having two events with different dominant frequencies are characteristic seismological phenomena observed during the fumarolic activity stage at Aso Volcano. These isolated tremors are called hybrid tremors (HBT) and comprise two parts: an initial part named the “HF-part” with a dominant frequency in the high-frequency region (approximately 10 Hz) and the following part named the “LF-part” with a dominant frequency in the low-frequency region (approximately 2 Hz). The LF-part is observed after the HF-part, and the HBT is accompanied by a long-period tremor (LPT). Hypocenters and source parameters are estimated using seismograms recorded at 64 stations around Nakadake crater. The amplitude distributions of all HF-parts have almost similar trends. Similarly, the amplitude distributions of all LF-parts have almost similar trends. However, the amplitude distributions of HF- and LF-parts are not similar. From these results, we proposed that the hypocenters and source parameters of HF- and LF-parts are not common, but each of them have common hypocenters and source parameters. The hypocenter region of HF-parts was estimated to be just beneath the fumarole region south of the 1st crater: the volume fluctuation is the major source factor. The hypocenter region of LF-parts is estimated to be at a depth of approximately 300 m beneath the first crater: the strike–slip component is the major source parameter. The hypocentral depth of LF-parts is located at the upper end of the crack estimated to be the source of the LPTs. The LPTs and HBTs are observed almost simultaneously. We consider that volcanic fluid is involved in the source mechanisms of both HBT and LPT.  相似文献   

Fast transform from geocentric to geodetic coordinates   总被引：3，自引：0，他引：3  

T. Fukushima 《Journal of Geodesy》1999,73(11):603-610

 A new iterative procedure to transform geocentric rectangular coordinates to geodetic coordinates is derived. The procedure solves a modification of Borkowski's quartic equation by the Newton method from a set of stable starters. The new method runs a little faster than the single application of Bowring's formula, which has been known as the most efficient procedure. The new method is sufficiently precise because the resulting relative error is less than 10⁻¹⁵, and this method is stable in the sense that the iteration converges for all coordinates including the near-geocenter region where Bowring's iterative method diverges and the near-polar axis region where Borkowski's non-iterative method suffers a loss of precision. Received: 13 November 1998 / Accepted: 27 August 1999  相似文献   

Application of spectral decomposition algorithm for mapping water quality in a turbid lake (Lake Kasumigaura,Japan) from Landsat TM data     

Youichi Oyama  Bunkei Matsushita  Takehiko Fukushima  Kazuo Matsushige  Akio Imai 《ISPRS Journal of Photogrammetry and Remote Sensing》2009,64(1):73-85

The remote sensing of Case 2 water has been far less successful than that of Case 1 water, due mainly to the complex interactions among optically active substances (e.g., phytoplankton, suspended sediments, colored dissolved organic matter, and water) in the former. To address this problem, we developed a spectral decomposition algorithm (SDA), based on a spectral linear mixture modeling approach. Through a tank experiment, we found that the SDA-based models were superior to conventional empirical models (e.g. using single band, band ratio, or arithmetic calculation of band) for accurate estimates of water quality parameters. In this paper, we develop a method for applying the SDA to Landsat-5 TM data on Lake Kasumigaura, a eutrophic lake in Japan characterized by high concentrations of suspended sediment, for mapping chlorophyll-a (Chl-a) and non-phytoplankton suspended sediment (NPSS) distributions. The results show that the SDA-based estimation model can be obtained by a tank experiment. Moreover, by combining this estimation model with satellite-SRSs (standard reflectance spectra: i.e., spectral end-members) derived from bio-optical modeling, we can directly apply the model to a satellite image. The same SDA-based estimation model for Chl-a concentration was applied to two Landsat-5 TM images, one acquired in April 1994 and the other in February 2006. The average Chl-a estimation error between the two was 9.9%, a result that indicates the potential robustness of the SDA-based estimation model. The average estimation error of NPSS concentration from the 2006 Landsat-5 TM image was 15.9%. The key point for successfully applying the SDA-based estimation model to satellite data is the method used to obtain a suitable satellite-SRS for each end-member.  相似文献   

Subsidence associated with the LUSI mud eruption,East Java,investigated by SAR interferometry     

Yo Fukushima  Jim Mori  Manabu Hashimoto  Yasuyuki Kano 《Marine and Petroleum Geology》2009

A mud volcano LUSI initiated its eruption on 29 May 2006, adjacent to a hydrocarbon exploration well in East Java. Ground subsidence in the vicinity of the LUSI eruptive vent was well recorded by a Synthetic Aperture Radar (SAR) PALSAR onboard the Japanese ALOS satellite. We apply an Interferometric SAR (InSAR) technique on ten PALSAR data scenes, acquired between 19 May 2006 and 21 May 2007, in order to obtain continuous maps of ground displacements around LUSI. Although the displacements in the area closest to the eruptive vent (spatial extension of about 1.5 km) are not detectable because of the erupted mud, all the processed interferograms indicate subsidence in an ellipsoidal area of approximately 4 km (north–south) × 3 km (east–west), centered at the main eruptive vent. In particular, interferograms spanning the first four months until 4 Oct. 2006 and the subsequent 46 days between 4 Oct. 2006 and 19 Nov. 2006 show at least about 70 cm and 80 cm of displacements away from the satellite, respectively. Possible causes of the subsidence, i.e., 1) loading effect of the erupted mud, 2) creation of a cylindrical mud conduit, and 3) pressure decrease and depletion of materials at depth, are investigated. The effects of the first two causes are found to be insufficient to explain the total amount of subsidence observed in the first six months. The third possibility is quantitatively examined using a boundary element approach by modeling the source of deformation as a deflating oblate spheroid. The spheroid is estimated to lie at depths of a few hundred to a thousand meters. The estimated depths are significantly shallower than determined from analyses of erupted mud samples; the difference is explained by presence of significant amount of inelastic deformation including compaction and downward transfer of material.  相似文献   

A model of river flow forecasting for a small forested mountain catchment     

Yoshihiro Fukushima 《水文研究》1988,2(2):167-185

A model for quantitatively expressing the hydrological cycle in a forested mountain catchment is proposed as a HYCYMODEL. HYCYMODEL is able to predict both short- and long-term hydrographs because the model parameters remain independent of time. It shows a good applicability for ten years of continuous data at both hourly and daily intervals for the Kiryu catchment—a forested mountain basin. Since HYCYMODEL does not need hydrograph separation between storm flow and base flow, it is a particularly attractive model.  相似文献   

The Fermi coordinate system in the post-Newtonian framework     

Toshio Fukushima 《Celestial Mechanics and Dynamical Astronomy》1988,44(1-2):61-75

The Fermi coordinate system (or the non-rotating proper reference frame) is studied within the framework of the post-Newtonian approximation of general relativity. First its general functional form is found to consist of three parts; (1) the motion of the space origin of the Fermi coordinate system relative to the background one, (2) the Fermi-Walker tetrad transported along the world line of the space origin, and (3) the spacelike geodesic starting from the space origin. Next, the post-Newtonian expressions of the latter two are obtained under the condition that the first is given. Then the full coordinate transformation formula connecting the Fermi coordinate system to the background one is derived explicitly. The effectiveness of the Fermi coordinate system is discussed and the effective region is found to be a cylinder with the radius of about 0.5 kpc for the Fermi coordinate system comoving with the Earth. The mathematical way to derive the generalized Fermi coordinate system which Ashby and Bertotti defined is also shown.  相似文献   

Seasonal dependences of geomagnetic variations in the polar region in connection with large-amplitude annual Z-variation at the geomagnetic pole     

Tsugunobu Nagai  Naoshi Fukushima 《Planetary and Space Science》1979,27(12):1513-1522

A noticeable depression of the vertical component Z of the geomagnetic field is observed in the polar cap in summer. From the statistical analysis of the equivalent overhead current patterns for daily geomagnetic variations in the summer and winter polar regions for three different conditions of IMF (interplanetary magnetic field), it was concluded that the annual variation of geomagnetic Z in the vicinity of the geomagnetic pole is attributed to the relative spatial shift of the twin-vortex current patterns over the polar cap from summer to winter. In winterthe clockwise current vortex in the dawn sector extends almost over the entire polar cap (except for the periods when the B_z component of IMF has a large positive value), and this will result in the positive variation of the Z-value at the geomagnetic pole. In summer the counter-clockwise current vortex in the dusk sector always extends over the whole polar cap even when B_z of IMF is positiveso that the variation of Z becomes negative. The persistent existence of current vortex in the dawn sector is important for the further study of magnetospheric convection when B_z is positive.  相似文献   

A physical‐model study of the statistics of seismic waveform fluctuations in random heterogeneous media     

C. Sivaji  O. Nishizawa  G. Kitagawa  Y. Fukushima 《Geophysical Journal International》2002,148(3):575-595

  相似文献   

Characterizing the flush of stream chemical runoff from forested watersheds     

Zhao Zhang  Fulu Tao  Peijun Shi  Wei Xu  Yu Sun  Takehiko Fukushima  Yuichi Onda 《水文研究》2010,24(20):2960-2970

A flush can be defined as stream chemical exhibiting higher concentrations during the prophase of a storm event at an event scale, or exhibiting progressively lower concentrations during several successive storms at a seasonal scale. Investigating the flush characteristics of chemical runoff from forested watersheds is important and helpful to understand the chemical dynamics as well as to design a sampling schedule strategy during storm events. Here, three parameters describing the flush characteristic are quantified and the flush characteristics of chemicals from four Japanese forested watersheds (Mie, Kochi, Nagano and Tokyo) were investigated at both event and seasonal levels. We found that the characteristics of the flush were complicated, and depended on the constituents of the hydrochemistry, climate and runoff quality. Generally, the flush occurs more readily for particulate components than for those in solution; the flush on nitrate‐nitrogen is weaker in regions of nitrogen saturation, such as Nagano and Tokyo, than in Mie and Kochi. Rainfall feature was the main factor controlling the flush of particulate components. However, the source available in a watershed plays a main role on the flushes of dissolve chemicals. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献