A development of the discontinuous deformation analysis for rock fall analysis

Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd.     

Energy and spectral analysis of confined solar flares from radio and X-ray observations

The energy and spectral shape of radio bursts may help us understand the generation mechanism of solar eruptions,including solar flares,coronal mass ejections,e...     

An upper limit of visible emission from the plume generated by the impact of the K fragment of SL9

We monitored Jupiter with a video movie camera attached to a 1.3 m reflector during the impact event of the K fragment of D/Shoemaker-Levy 9. The camera was sensitive to visible radiation. A plume, which was generated by the impact and prominent in infrared wavelength, was not detected. Jovian images with hypothetical plumes were synthesized to examine how luminous the plume should have been to be recognized. It has been revealed that the plume was no more than twice as bright as Io in visible wavelength.     

Abundances of 23 field RR Lyrae stars

We present stellar parameters and abundances of 15 elements （Na, Mg, A1, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Ni, Cu, Y and Ba） for 23 field RR Lyrae variables based on high-resolution （R ~ 60 000） and high signal-to-noise （S/N~ 200） spectra obtained using the High Dispersion Spectrograph on the Subaru Telescope. Six stars in the sample have more than one spectrum observed at different pulsation phases. The derived abundance ratios of [X/Fe] for 14 elements （except for Ba） do not vary during the pulsation cycle. An interesting curve of [Ba/Fe] versus phase is detected for the first time and it shows decreasing [Ba/Fe] with increasing temperature at a given metallicity. Combining with data in the literature, abundances of most RR Lyrae stars as a function of [Fe/H] follow the same trends as those of dwarf stars, but [Sc/Fe] and [Y/Fe] ratios of RR Lyrae stars at solar metallicity are lower than those of dwarf stars. The kinematics of RR Lyrae stars indicate that three comparatively metal-rich RR Lyrae stars might originate from the thick disk and they show higher [a/Fe] ratios than RR Lyrae stars with thin disk kinematics. Among 23 RR Lyrae stars, two special objects are found with abnormal abundances; TV Lib has high [a/Fe], [Sc/Fe], [Y/Fe] and [Ba/Fe] ratios while TW Her has solar [a/Fe] but significantly lower [Sc/Fe], [Y/Fe] and [Ba/Fe] ratios as compared with other RR Lyrae stars.     

Quaternary Volcanic-Related Acid Hydrothermal Alteration in the Ugusu Silica Stone Deposit, Western Izu Peninsula, Central Japan: Geology and Alteration

Acid alteration areas accompanying Quaternary volcanoes are widespread in the western Izu Peninsula, central Japan. The Ugusu alteration area is the largest among them and is mined for “silica stone” at the silica body in the core of the alteration area. Silica zone defined by previous studies is subdivided into highly leached brecciated silica zone and residual silica zone, which extend along a NNW‐SSE–NNE‐SSW direction of fractures/faults. Fe‐rich, alunite, advanced argillic alteration, and intermediate argillic alteration zones occur toward the outside surrounding the two silica zones. The ascent of acid hydrothermal fluid would be responsible for the silica zones and surrounding alteration zones at an earlier stage, while the hydrothermal brecciation and silica veins were caused by a limited supply of silica‐saturated fluids at later stages. Based on the available mineral stability relations and fluid inclusion thermometries, the formation temperatures are estimated to be: >300°C for the residual silica zone; >290°C for the diaspore association in the advanced argillic alteration zone; and <260°C for the kaolinite association in the intermediate argillic alteration zone. The later stage quartz druses have been formed at 200–260°C. The Ugusu–Fukata acid hydrothermal systems were active at 1.5–1.2 Ma, which were temporally related to the Tanaba Andesite volcanism. Hydrothermal system at the Seikoshi gold–silver deposit survived until 0.7 Ma after the volcanism. In the western half of the Izu Peninsula, subduction of the Philippine Sea plate underneath the Suruga Trough caused nearly N–S‐trending maximum horizontal compressive stress (σ_Hmax) and the resultant formation of similarly trending alteration areas and Au‐Ag vein‐type deposits in the Ugusu‐Toi‐Seikoshi area. From a practical viewpoint, at the Ugusu silica stone deposit, the fracture‐controlled vertical morphology of the silica body provides an important guide for exploration. Because the alteration zones occur both in the lower and upper sides of the silica bodies, it is important to make sure to which side the alteration zones correspond.     

Detection and accuracy of landslide movement by InSAR analysis using PALSAR-2 data

Interferometric synthetic aperture radar (InSAR) analysis is a radar technique for generating large-area maps of ground deformation using differences in the phase of microwaves returning to a satellite. In recent years, high-resolution SAR sensors have been developed that enable small-scale slope deformation to be detected, such as the partial block movement of a landslide. The L-band SAR (PALSAR-2) is mounted on Advanced Land Observing Satellite-2 (ALOS-2), which was launched on 24 Mar. 2014. Its main improvements compared with ALOS are enhanced resolution of as high as 3 m with a high-frequency recurrence period (14 days). Owing to its high resolution and the use of the L-band, PALSAR-2 can obtain reflective data passing through a tree canopy surface, unlike the other synthetic aperture radars. Therefore, the coherence of InSAR in mountainous forest areas is less likely to decrease, making it advantageous for the extraction of slope movement. In this study, to verify the accuracy of InSAR analysis using PALSAR-2 data, we compared the results of InSAR analysis and the measurement of the displacement in a landslide by global navigation satellite system (GNSS) observation. It was found that the average difference between the displacements obtained by InSAR analysis and the field measurements by GNSS was only 15.1 mm in the slant range direction, indicating the high accuracy of InSAR analysis. Many of the areas detected by InSAR analysis corresponded to the locations of surface changes due to landslide activity. Additionally, in the areas detected by InSAR analysis using multiple datasets, the ground changes due to landslide movement were confirmed by site investigation.     

Magmatic processes forming replacement textures with fluorite alignments in feldspars in an evolved trachyte from Oki-Dogo Island,Sea of Japan

Internal microtextures of ternary alkali feldspars in sanidine trachyte from Oki-Dogo Island were examined using an electron microprobe analyzer, a scanning electron microscope, a transmission electron microscope and cathodoluminescence instruments, to develop the understanding of volcanic processes of alkaline magmas related to feldspar crystallization. The examined trachyte is an evolved rock of the Oki-Dogo Pliocene trachyte group. Its phenocryst feldspars are commonly associated with lamellar-wavy-domain textures with scales approximately from 100 nm up to several hundreds of μm that show complex and gradual variations in composition: however, anti-rapakivi zoning textures common in other Oki-Dogo alkaline rocks are almost completely absent in the trachyte. These textures are produced by extensive magmatic ion-exchange replacement reactions progressively advanced in the evolved magma. Characteristic braided fluorite alignments are developed consistently with lamellar-wavy-domain textures in phenocryst feldspars, and similar braided alignments are also present in groundmass feldspars with complicated microtextures. Most of fluorite grains are <100 nm in diameter, and the patterns of braided fluorite alignments vary greatly in individual feldspars. The whole occurrence of the feldspar microtextures represents an extreme example of diffusion-controlled replacement reactions, progressively advanced in the dry (relatively anhydrous) trachyte magma. The genetic processes forming fluorite alignments in feldspars are related to magma compositions, especially F and P contents, and the crystallization of F-bearing minerals, especially of fluorapatite.     

Modeling the interaction between fluid flow and particle movement with discontinuous deformation analysis (DDA) method

Coupled hydro‐mechanical processes in granular media represent the interaction between solid particles movement and fluid flow during external and/or internal loading. The processes attract attention in geotechnical engineering since they cause many disasters such as landslide, slope collapse, boiling/quicksand and soil liquefaction, etc. This study presents a numerical method combining discontinuous deformation analysis (DDA) for mechanical calculation and finite element method for fluid flow simulation to model the interaction between solid particles' movement and fluid flow from microscopic point of view. The term of hydraulic pressure was newly formulated and introduced to the original DDA. The simplified examples are given to verify the new method, and the computational results correlate well with the theoretical calculations, further development is also considered. Copyright © 2010 John Wiley & Sons, Ltd.     

Laser-Dye Imaging of the pH Field in a Laboratory Experiment

We established a precise DeLIF (Dual Emission Laser Light Fluorescence) system for the instant measurement of the pH field in laboratory experiments. This method is very important in studies of the biological impacts of CO₂ sequestration, i.e., effects of pH on organisms. The technique is demonstrated by applying it to the pH measurement around a liquid CO₂ droplet, which has a large gradient of CO₂ concentration near its interface, the concentration value on the boundary being known as the solubility of CO₂. In the Laser Induced Fluorescence (LIF) method using quinine and an ultraviolet laser light source, the relatively unstable pulse of the laser source decreased the measurement reliability. We therefore investigated many kinds of fluorescent dye and demonstrated the DeLIF measurement, identifying a process for cancelling fluctuations of laser light intensity. Quinine and sulforhodamine were used as fluorescence dyes. The molar absorption coefficient and emission intensity of quinine depend on the pH of aqueous solutions, while those of sulforhodamine are not as much affected by pH. The ratio of these fluorescence intensities represents pH at each point of the image. As a result, the pH field was successfully measured with an uncertainty of only approximately ±0.005, for 3.0 < pH. This revised version was published online in July 2006 with corrections to the Cover Date.     

Accurate ocean tide modeling in southeast Alaska and large tidal dissipation around Glacier Bay

An accurate prediction of ocean tides in southeast Alaska is developed using a regional, barotropic ocean model with a finite difference scheme. The model skill is verified by the observational tidal harmonics in southeast Alaska including Glacier Bay. The result is particularly improved in Glacier Bay compared to the previous model described by Foreman et al. (2000). The model bathymetry dominates the model skill. We re-estimate tidal energy dissipation in the Alaska Panhandle and suggest a value for tidal energy dissipation of 3.4 GW associated with the M₂ constituent which is 1.5 times the estimation of Foreman et al. (2000). A large portion of the M₂ energy budget entering through Chatham Strait is dissipated in the vicinity of Glacier Bay. Moreover, it is shown that the developed model has the potential to correct the ocean tide loading effect in geodetic data more efficiently than the model of Foreman et al. (2000), especially around Glacier Bay.