Observation of the Kyucho in the Bungo Channel by HF radar

Observations of sea surface currents by HF radar were carried out in the Bungo Channel in summer 1992. The current ellipses of M₂ constituent obtained by the observational results agree quite well with those obtained by the ADCP observations, showing that the accuracy of the HF radar measurements is of the same level as ADCP. The results revealed the current structures and their change with the Kyucho in detail. The Kyucho is influenced by the complicated coastal geometry and does not propagate straightly into the Bungo Channel. It propagates further inward after charging the coastal bays with warm water. The current directions change largely, since the currents turn around the stagnant region in the bay filled with the warm water. The northward intrusion begins to be weakened in the southern part of the channel, while it still persists in the northern part. The northward current speeds of the observed Kyucho are about 50 cm/s and sometimes attain 60 to 70 cm/s.     

Developing a hybrid multi‐model for peak flood forecasting

The present study aims to develop a hybrid multi‐model using the soft computing approach. The model is a combination of a fuzzy logic, artificial neural network (ANN) and genetic algorithm (GA). While neural networks are low‐level computational structures that perform well dealing with raw data, fuzzy logic deal with reasoning on a higher level by using linguistic information acquired from domain experts. However, fuzzy systems lack the ability to learn and cannot adjust themselves to a new environment. Moreover, experts occasionally make mistakes and thus some rules used in a system may be false. A network type structure of the present hybrid model is a multi‐layer feed‐forward network, the main part is a fuzzy system based on the first‐order Sugeno fuzzy model with a fuzzification and a defuzzification processes. The consequent parameters are determined by least square method. The back‐propagation is applied to adjust weights of network. Then, the antecedent parameters of the membership function are updated accordingly by the gradient descent method. The GA was applied to select the fuzzy rule. The hybrid multi‐model was used to forecast the flood level at Chiang Mai (under the big flood 2005) and the Koriyama flood (2003) in Japan. The forecasting results are evaluated using standard global goodness of fit statistic, efficient index (EI), the root mean square error (RMSE) and the peak flood error. Moreover, the results are compared to the results of a neuro‐genetic model (NGO) and ANFIS model using the same input and output variables. It was found that the hybrid multi‐model can be used successfully with an efficiency index (EI) more than 0·95 (for Chiang Mai flood up to 12 h ahead forecasting) and more than 0·90 (for Koriyama flood up to 8 h ahead forecasting). In general, all of three models can predict the water level with satisfactory results. However, the hybrid model gave the best flood peak estimation among the three models. Therefore, the use of fuzzy rule base, which is selected by GA in the hybrid multi‐model helps to improve the accuracy of flood peak. Copyright © 2009 John Wiley & Sons, Ltd.     

Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

Evaluation of the resistance of coastal Pinus thunbergii Parlat. forests to the tsunami fluid force in Japan

Natural Hazards - Coastal forests can protect coastal areas from the effects of tsunamis. However, the coastal forests of Tohoku were destroyed by the Great East Japan Tsunami, which was triggered...     

A Method to Suppress Isobaric and Polyatomic Interferences for Measurements of Highly Siderophile Elements in Desilicified Geological Samples

Sample decomposition using inverse aqua regia at elevated temperatures and pressures (e.g., Carius tube or high‐pressure asher) is the most common method used to extract highly siderophile elements (HSEs: Ru, Rh, Pd, Re, Os, Ir, Pt and Au) from geological samples. Recently, it has been recognised that additional HF desilicification is necessary to better recover HSEs, potentially contained within silicate or oxide minerals in mafic samples, which cannot be dissolved solely by inverse aqua regia. However, the abundance of interfering elements tends to increase in the eluent when conventional ion‐exchange purification procedures are applied to desilicified samples. In this study, we developed an improved purification method to determine HSEs in desilicified samples. This method enables the reduction of the ratios of isobaric and polyatomic interferences, relative to the measured intensities of HSE isotope masses, to less than a few hundred parts per million. Furthermore, the total procedural blanks are either comparable to or lower than conventional methods. Thus, this method allows accurate and precise HSE measurements in mafic and ultramafic geological samples, without the need for interference corrections. Moreover, the problem of increased interfering elements, such as Zr for Pd and Cr for Ru, is circumvented for the desilicified samples.     

Interhemispheric observations of field line resonance frequencies as a continuous function of ground latitude in the auroral zones

In the austral summer of 2006–2007, the 48th Japanese Antarctic Research Expedition (JARE-48) installed two unmanned low-power magnetometers to form a closely spaced magnetometer network in combination with the permanent sites at Japan's Syowa Station in Antarctica. To identify field line resonances (FLRs), gradient methods are applied to the data from three adjacent sites in Antarctica and data from conjugate points in Antarctica and Iceland. By analyzing the data from the Antarctic and Icelandic sites individually, the structure of FLRs with high coherence is clearly identified. However, by analyzing the data from closely spaced Antarctic sites, it is more difficult to identify the signature of FLRs because of the inclusion of multiple signals related to the local geomagnetic pulsations over a broad frequency range. The frequency and resonance width of FLRs are determined by applying the amplitude phase gradient method (APGM) to the data from Antarctic sites. This yields the eigenfrequency as a continuous function of ground latitudes in the area surrounding Syowa Station. The mass density in the equatorial region at the L of the auroral zones is estimated from the obtained FLR frequency by numerically solving the standing Alfvén wave equation. The mass density thus obtained is consistent with observational results from previous in situ measurements by spacecraft. The results of the present study demonstrate that data from geomagnetic conjugate points are helpful in identifying FLR in cases in which the magnetometers are too close to each other to enable identification. Once FLR is identified, APGM can be applied to the identified FLR, yielding the FLR frequency as a continuous function of ground latitudes. Therefore, the magnetospheric equatorial mass density is readily estimated with high spatial resolution.     

Density Current Descending along Continental Slope and the Associated Deep Water Formation: Two-Dimensional Numerical Experiments with a Nonhydrostatic Model

Numerical experiments with a two-dimensional nonhydrostatic ocean model have been carried out to investigate the dynamical process of descending density current on a continental slope. The associated deep water formation has been also examined by tracking labeled particles. The descending flow along the continental slope occurs in the bottom Ekman layer. The net pressure gradient determining the volume transport consists of not only the pressure gradient due to density deviation but also the surface pressure gradient due to the depth-mean alongshore flow. Since these constituents have the opposite signs and strengthen each other, the oscillation with an alternation of intense up- and downslope flows appears around the shelf break. This temporal variation of the flow field causes the effective mixing on the slope between descending shelf and interior waters and forms the deep water as a mixture of them at a ratio of about 1:3. The present result is applied to the slope current around Antarctica, using velocity and density fields calculated by an ocean general circulation model. The Ekman volume transport is estimated at 0.97 Sv (1 Sv = 10⁶ m³s^–1) in the Weddell Sea, 0.35 Sv in the Ross Sea, and 1.8 Sv in total. About 70% of them is attributed to the depth-mean alongshore flow, such as the East Wind Drift and the Weddell Gyre driven by the wind. This suggests that the pressure gradient due to other factors than density deviation may play an important role in the deep and bottom water formation in the actual oceans.     

Geometry of the Nojima Fault at Nojima-Hirabayashi, Japan – I. A Simple Damage Structure Inferred from Borehole Core Permeability

The 1995 Kobe (Hyogo-ken Nanbu) earthquake, M = 7.2, ruptured the Nojima fault in southwest Japan. We have studied core samples taken from two scientific drillholes that crossed the fault zone SW of the epicentral region on Awaji Island. The shallower hole, drilled by the Geological Survey of Japan (GSJ), was started 75 m to the SE of the surface trace of the Nojima fault and crossed the fault at a depth of 624 m. A deeper hole, drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) was started 302 m to the SE of the fault and crossed fault strands below a depth of 1140 m. We have measured strength and matrix permeability of core samples taken from these two drillholes. We find a strong correlation between permeability and proximity to the fault zone shear axes. The half-width of the high permeability zone (approximately 15 to 25 m) is in good agreement with the fault zone width inferred from trapped seismic wave analysis and other evidence. The fault zone core or shear axis contains clays with permeabilities of approximately 0.1 to 1 microdarcy at 50 MPa effective confining pressure (10 to 30 microdarcy at in situ pressures). Within a few meters of the fault zone core, the rock is highly fractured but has sustained little net shear. Matrix permeability of this zone is approximately 30 to 60 microdarcy at 50 MPa effective confining pressure (300 to 1000 microdarcy at in situ pressures). Outside this damage zone, matrix permeability drops below 0.01 microdarcy. The clay-rich core material has the lowest strength with a coefficient of friction of approximately 0.55. Shear strength increases with distance from the shear axis. These permeability and strength observations reveal a simple fault zone structure with a relatively weak fine-grained core surrounded by a damage zone of fractured rock. In this case, the damage zone will act as a high-permeability conduit for vertical and horizontal flow in the plane of the fault. The fine-grained core region, however, will impede fluid flow across the fault.     

A preparation zone for volcanic explosions beneath Naka-dake crater,Aso volcano,as inferred from magnetotelluric surveys

The 1st crater of Naka-dake, Aso volcano, is one of the most active craters in Japan, and known to have a characteristic cycle of activity that consists of the formation of a crater lake, drying-up of the lake water, and finally a Strombolian-type eruption. Recent observations indicate an increase in eruptive activity including a decrease in the level of the lake water, mud eruptions, and red hot glows on the crater wall. Temporal variations in the geomagnetic field observed around the craters of Naka-dake also indicate that thermal demagnetization of the subsurface rocks has been occurring in shallow subsurface areas around the 1st crater. Volcanic explosions act to release the energy transferred from magma or volcanic fluids. Measurement of the subsurface electrical resistivity is a promising method in investigating the shallow structure of the volcanic edifices, where energy from various sources accumulates, and in investigating the behaviors of magma and volcanic fluids. We carried out audio-frequency magnetotelluric surveys around the craters of Naka-dake in 2004 and 2005 to determine the detailed electrical structure down to a depth of around 1 km. The main objective of this study is to identify the specific subsurface structure that acts to store energy as a preparation zone for volcanic eruption. Two-dimensional inversions were applied to four profiles across the craters, revealing a strongly conductive zone at several hundred meters depth beneath the 1st crater and surrounding area. In contrast, we found no such remarkable conductor at shallow depths beneath the 4th crater, which has been inactive for 70 years, finding instead a relatively resistive body. The distribution of the rotational invariant of the magnetotelluric impedance tensor is consistent with the inversion results. This unusual shallow structure probably reflects the existence of a supply path of high-temperature volcanic gases to the crater bottom. We propose that the upper part of the conductor identified beneath the 1st crater is mainly composed of hydrothermally altered zone that acts both as a cap to upwelling fluids supplied from deep-level magma and as a floor to infiltrating fluid from the crater lake. The relatively resistive body found beneath the 4th crater represents consolidated magma. These results suggest that the shallow conductor beneath the active crater is closely related to a component of the mechanism that controls volcanic activity within Naka-dake.