Measurement of bubble size distributions in vesiculated rocks with implications for quantitative estimation of eruption processes

This paper outlines methods for determining a bubble size distribution (BSD) and the moments of the BSD function in vesiculated clasts produced by volcanic eruptions. It reports the results of applications of the methods to 11 natural samples and discusses the implications for quantitative estimates of eruption processes. The analysis is based on a quantitative morphological (stereological) method for 2-dimensional imaging of cross-sections of samples. One method determines, with some assumptions, the complete shape of the BSD function from the chord lengths cut by bubbles. The other determines the 1st, 2nd and 3rd moments of distribution functions by measurement of the number of bubbles per unit area, the surface area per unit volume, and the volume fraction of bubbles. Comparison of procedures and results of these two distinct methods shows that the latter yields rather more reliable results than the former, though the results coincide in absolute and relative magnitudes. Results of the analysis for vesiculated rocks from eleven subPlinian to Plinian eruptions show some interesting systematic correlations both between moments of the BSD and between a moment and the eruption column height or the SiO₂ content of magma. These correlations are successfully interpreted in terms of the nucleation and growth processes of bubbles in ascending magmas. This suggests that bubble coalescence does not predominate in sub-Plinian to Plinian explosive eruptions. The moment-moment correlations put constraints on the style of the nucleation and growth process of bubbles. The scaling argument suggests that a single nucleation event and subsequent growth with any kind of bubble interaction under continuous depressurization, which leads to an intermediate growth law between the diffusional growth ( ) at a constant depressurization rate and the Ostwald ripening ( ) under a constant pressure, where R_m and t are the mean radius of bubble and the effective time of diffusion respectively, occurred in the eruptions. It is emphasized that the BSD in vesiculated rocks from terrestrial volcanoes can be used to estimate quantitatively eruption processes such as the initial saturation pressure and magma ascent velocity in a volcanic conduit.     

Multiple stages of phengite formation in Sanbagawa schists

Summary Basic and pelitic schists of the garnet and biotite zones in the Sanbagawa belt in Shikoku, Japan, commonly contain phengite developed at different stages of metamorphism. Textures of the phengite and associated minerals show that Al-rich phengite is a prograde product and Al-poor phengite is a retrograde product.

---

Mehrfache Stadien der Phengit-Bildung in den Sanbagawa-Schiefern

Zusammenfassung Basische und pelitische Schiefer der Granat- und Biotitzonen im Sanbagawa Gürtel in Shikoku, Japan, führen Phengite, die während verschiedener Stadien der Metamorphose entstanden sind. Texturen der Phengite und assoziierter Minerale zeigen, daß Al-reicher Phengit das Produkt prograder Metamorphose ist, während Al-armer Phengit auf retrograde Metamorphose zurückgeht.

---

---

With 5 Figures

---

On leave from the Geological Survey of Slovakia, Spissa, Spisska, Nova Ves, Slovakia     

Strain and tilt changes measured during a water injection experiment at the Nojima Fault zone, Japan

Abstract In order to make geophysical and geological investigations of the Nojima Fault on Awaji Island, Japan, three boreholes measuring 1800 m, 800 m and 500 m deep were drilled into the fault zone. The fault is one of the seismic source faults of the 1995 Hyogo-ken Nanbu earthquake of M 7.2. A new multicomponent borehole instrument was installed at the bottom of the 800 m borehole and continuous observations of crustal strain and tilt have been made using this instrument since May 1996. A high-pressure water injection experiment within the 1800 m borehole was done in February and March 1997 to study the geophysical response, behavior, permeability, and other aspects of the fault zone. The injection site was located approximately 140 m horizontally and 800 m vertically from the instrument. Associated with the water injection, contraction of approximately 0.7 × 10⁻⁷ str (almost parallel to the fault) and tilt of approximately 1 × 10^-7 rad in the sense of upheaval toward the injection site were observed. In addition to these controlled experiments, the strainmeter and tiltmeter also recorded daily variations. We interpret strain and tilt changes to be related to groundwater discharge and increased ultra-micro seismicity induced by the injected water.     

Mesozoic radiolarians and radiolarian-bearing sequences in the circum-Pacific regions: A report of the symposium-'Radiolarians and Orogenic Belts'

Abstract This paper contains extended abstracts of the seven papers presented at the symposium 'Radiolarians and Orogenic Belts' held at the seventh meeting of the International Association of Radiolarian Paleontologists (INTERRAD). Important results of the symposium include the following: (1) Upper Paleozoic and Mesozoic cherts are widely distributed within accretionary complexes in the circum-Pacific orogenic belt. Radiolarian dating reveals that long durations of chert sedimentation in a pelagic environment are recorded on both sides of Pacific-rim accretionary complexes (e.g. New Zealand, Japan, Russian Far East, Canadian Cordillera). (2) Triassic radiolarian faunas from New Zealand and the Omolon Massif, northeast Siberia are similar in composition and are characterized by the absence of typical Tethyan elements. This suggests that radiolarian faunal provincialism may have been established as early as the Triassic. High-latitude radiolarian taxa exhibit a bi-polar distribution pattern. (3) The Lower Triassic interval in chert dominant pelagic sequences is mechanically weaker than other levels and acted as a décollement in accretionary events. This lithologic. contrast in physical property is considered to reflect radiolarian evolution, such as the end-Permian mass extinction.     

Behaviour of subaqueous tunnels during earthquakes

Three-dimensional models of a subaqueous tunnel were built on a shaking table and vibrated for the purpose of investigating the dynamic behaviour of the tunnel. The material of the soft surface layer was gelatin gel and that of the tunnel was silicon rubber. In addition, earthquake observations at the site of the real subaqueous tunnel have been carried out using accelerometers and strain gauges set on the walls of the tunnel. At the site, micro-tremors were also measured. In these studies, it was noted that non-uniform displacements of the ground along the tunnel axis were caused by variations of the dynamic behaviour of the ground, and that the tunnel deforms correspondingly to the varying deformations of the ground. From these experimental results, a mathematical model of the subaqueous tunnel was formulated for anti-seismic design studies.     

Solar Wind Speed and Expansion Rate of the Coronal Magnetic Field in Solar Maximum and Minimum Phases

Relationships between solar wind speed and expansion rate of the coronal magnetic field have been studied mainly by in-ecliptic observations of artificial satellites and some off-ecliptic data by Ulysses. In this paper, we use the solar wind speed estimated by interplanetary scintillation (IPS) observations in the whole heliosphere. Two synoptic maps of SWS estimated by IPS observations are constructed for two Carrington rotations CR 1830 and 1901; CR 1830 starting on the 11th of June, 1990 is in the maximum phase of solar activity cycle and CR 1901 starting on the 29th of September, 1995 is in the minimum phase. Each of the maps consist of 64800 (360×180) data points. Similar synoptic maps of expansion rate of the coronal magnetic field (RBR) calculated by the so-called potential model are also constructed under a radial field assumption for CR 1830 and CR1901. Highly significant correlation (r=–0.66) is found between the SWS and the RBR during CR1901 in the solar minimum phase; that is, high-speed winds emanate from photospheric areas corresponding to low expansion rate of the coronal magnetic field and low speed winds emanate from photospheric areas of high expansion rate. A similar result is found during CR 1830 in solar maximum phase, though the correlation is relatively low (r=–0.29). The correlation is improved when both the data during CR 1830 and CR 1901 are used together; the correlation coefficient becomes –0.67 in this case. These results suggest that the correlation analysis between the SWS and the RBR can be applied to estimate the solar wind speed from the expansion rate of the coronal magnetic field, though the correlation between them may depend on the solar activity cycle. We need further study of correlation analysis for the entire solar cycle to get an accurate empirical equation for the estimation of solar wind speed. If the solar wind speed is estimated successfully by an empirical equation, it can be used as an initial condition of a solar wind model for space weather forecasts.     

On-going orogeny in the outer-arc of the Timor–Tanimbar region, eastern Indonesia

The Timor–Tanimbar islands of eastern Indonesia form a non-volcanic arc in front of a 7 km deep fore-arc basin that separates it from a volcanic inner arc. The Timor–Tanimbar Islands expose one of the youngest high P/T metamorphic belts in the world, providing us with an excellent opportunity to study the inception of orogenic processes, undisturbed by later tectonic events.Structural and petrological studies of the high P/T metamorphic belt show that both deformation and metamorphic grade increase towards the centre of the 1 km thick crystalline belt. Kinematic indicators exhibit top-to-the-north sense of shear along the subhorizontal upper boundaries and top-to-the-south sense in the bottom boundaries of the high P/T metamorphic belt. Overall configuration suggests that the high P/T metamorphic rocks extruded as a thin sheet into a space between overlying ophiolites and underlying continental shelf sediments. Petrological study further illustrates that the central crystalline unit underwent a Barrovian-type overprint of the original high P/T metamorphic assemblages during wedge extrusion, and the metamorphic grade ranged from pumpellyite-actinolite to upper amphibolite facies.Quaternary uplift, marked by elevation of recent reefs, was estimated to be about 1260 m in Timor in the west and decreases toward Tanimbar in the east. In contrast, radiometric ages for the high P/T metamorphic rocks suggest that the exhumation of the high P/T metamorphic belt started in western Timor in Late Miocene time and migrated toward the east. Thus, the tectonic evolution of this region is diachronous and youngs to the east. We conclude that the deep-seated high P/T metamorphic belt extrudes into shallow crustal levels as a first step, followed by doming at a later stage. The so-called ‘mountain building’ process is restricted to the second stage. We attribute this Quaternary rapid uplift to rebound of the subducting Australian continental crust beneath Timor after it achieved positive buoyancy, due to break-off of the oceanic slab fringing the continental crust. In contrast, Tanimbar in the east has not yet been affected by later doming. A wide spectrum of processes, starting from extrusion of the high P/T metamorphic rocks and ending with the later doming due to slab break-off, can be observed in the Timor–Tanimbar region.     

Inclusion-trail geometry of albite porphyroblasts in a fold structure in the Sambagawa Belt, central Shikoku, Japan

The orientation of straight inclusion trails within albite porphyroblasts from basic schists has been measured around a north-closure fold, in the Besshi district of the Sambagawa Belt, central Shikoku, Japan. The porphyroblasts are aligned with their longest dimension parallel to both the subhorizontal, east–west-directed mineral lineation and to the fold axis. There is a systematic variation in inclusion-trail geometry between the upper (northern) and lower (southern) fold limbs. The shear sense deduced from quartz c-axis fabrics is top-to-the-west in the upper limb and top-to-the-east in the lower limb. Based on observed variations in porphyroblast inclusion trails, the structural history can be modelled as follows: (i) shear flow caused east–west stretching and folding of the metamorphic zonation; (ii) east–west ductile shear resulted in opposing senses of shear in the upper and lower limbs as the eclogite body situated in the core of the fold was extruded to the east.     

Enhanced chlorophyll associated with island-induced cyclonic eddies in the eastern channel of the Tsushima Straits

The relationship between island-induced cyclonic eddies and chlorophyll a (chl-a) was investigated using field data and satellite images in the eastern channel of the Tsushima Straits. The maximum chl-a concentration around the leeward eddy of the Tsushima Islands was two or three times greater than that of outside the eddy. Two different mechanisms of chl-a enhancement associated with island-induced cyclonic eddies were found in the post-bloom periods. In summer, when nutrients were depleted in the surface layer, eddy pumping increased the nutrient supply in the euphotic zone, resulting in enhanced chl-a around the shallow thermocline near the eddy core. In late autumn, when the mixed layer deepened over the euphotic zone, the mixed layer depth became shallow due to the doming effect of the cyclonic eddy, therefore, the improved irradiance condition led to an increase in the chl-a concentration in the surface mixed layer. Nighttime satellite visible images showed a number of fishing vessels in the lee region of the Tsushima Islands, implying that the enhanced phytoplankton biomass may have resulted in good feeding conditions for fishes and squids in the Tsushima Straits.