Wave dispersion and optimal mass modelling for one-dimensional periodic structures

Discrete analysis methods are frequently used for the study of the structure and soil. However, the assumption of the displacement interpolation function makes the waves dispersive, which means the numerical dispersion. The wave dispersion induced by the discretization depends on the mass modelling. Also, the existence of added lumped masses makes waves dispersive even for the continuum modelling. In order to examine these wave dispersions, a one-dimensional periodic structure is adopted as an analysis model and the dynamic transfer matrix method is applied. A wave solution and a finite element solution are used for the evaluation of the transfer matrix. The phase and group velocities in the structure are explicitly represented. These values are compared among the continuum modelling and the discretization modelling in which several consistent mass ratios are adopted. The optimal consistent mass ratio, which makes the wave velocity of the discrete model the same as that of the continuum model, is newly developed here. The validity of this mass modelling technique is presented by examining the frequency response function and impulse response function.     

Uniformity of Longitudinal Distribution of Coronal Mass Ejections

It was found that the heliographic longitudinal distribution of the occurrence of coronal mass ejections (CMEs) was not uniform around the solar minimum of cycle 21. This is confirmed by using the ² test for goodness of fit. Sunspot groups and flares show similar non-uniform distributions during the same period. This suggests that there was a longitudinal range where both large-scale and small-scale solar magnetic activities were very low around the solar minimum of cycle 21. The longitudinal range with low activity was more than 150 degrees and continued for approximately four rotations. As solar activity increased, CMEs occurred almost uniformly in heliographic longitude.     

Degradation ofn-paraffin mixture by marine microorganisms in enriched seawater medium

The ability to degraden-paraffin mixture of two bacterial strains,Caulobacter sp. andFlavobacterium sp., isolated from sea water of Tokyo Bay was studied experimentally in the enriched seawater (ESW) medium. These bacteria degraded actively the mixture ofn-tridecane,n-tetradecane,n-pentadecane andn-hexadecane. The maximum rate of degradation was observed after a lag period of 2 to 8 day and these bacteria were found to degrade then-paraffin mixture at rates calculated to be in a range from 3.3×10^?12 to 3.4×10^?11 mg-oil cell^?1 h^?1 at 20°C. The maximum degradation rate,r _m mg-oil l^?1 h^?1, was correlated with the amount of the initial totaln-paraffin,S mg-oil l^?1, as expressed by the following equation: $$rm = (rm)\max \left( {\frac{S}{{S + Km}}} \right)$$ where (r _m )_max denotes the largest value ofr _m whenn-paraffin exists in large excess andK _m is a constant and represents the amount ofn-paraffin at which the degradation rate,r _m , reaches 1/2 of its largest value, (r _m )_max. The values of (r _m )_max andK _m were calculated to be as follows: In the case ofCaulobacter sp. (strain KM-1), (r _m )_max=6.0 mg-oil l^?1 h^?1 andK _m =191 mg-oillesw ^?1; in the case ofFlavobacterium sp., (r _m )_max=5.47 mg-oil l^?1 h^?1 andK _m =152 mg-oillesw ^?1.     

Soft X-Ray and White-Light Coronal Observations of the Post-Flare Arch on 2 November 1991

A giant post-flare arch observed on 2-3 November 1991 was analyzed using the soft X-ray telescope (SXT) on board Yohkoh and the Mark III (MK3) K-coronameter at the High Altitude Observatory/Mauna Loa Solar Observatory. The rising arch was observed in both soft X-ray and K-corona observations. The estimated rising speed from the MK3 observation was approximately 4 km s^-1. A V-shaped depression area was observed on the south side of the giant arch. Change in the K-corona observations was faint while the arch was rising. According to the solar wind observations by the Pioneer Venus Orbiter and the Interplanetary Cometary Explorer, this giant arch event may have been associated with an interplanetary shock.     

Organohalogen and organotin compounds in killer whales mass-stranded in the Shiretoko Peninsula, Hokkaido, Japan

Blubber and liver samples were obtained for analysis of wide ranges of contaminants from killer whales (Orcinus orca) which were locked away in drifting sea ice on the coast of Rausu, the Shiretoko Peninsula in Eastern Hokkaido, Japan in February 2005. Among the organohalogen compounds analyzed, DDTs were the predominant contaminants with concentrations ranging from 28 to 220 microg/g on a lipid-weight basis followed by PCBs and other organochlorine pesticides. PBDEs levels were two or three orders of magnitude lower than those of PCBs and DDTs. 2,3,7,8-Tetrachlorodibenzo-p-dioxin toxic equivalents (TEQs) derived by WHO mammal-TEF in killer whales were in the range of 110-440 pgTEQ/g. Mono-ortho coplanar PCBs contributed to 75-98% of total TEQs, indicating coplanar PCBs are significant contaminants for risk assessment in this species. The fact that hepatic residue levels of butyltins (from 13 to 770 ng/g wet weight) were much higher than those of phenyltins may be reflecting extensive use of tributyltin as antifouling paint.     

航空激光雷达用于森林测量的数据处理方法研究

航空激光雷达是最近十几年开发的最新测量技术.利用激光扫描器所获得的高度数据可以建立精密的地形模型,同时可以获得森林生产力的信息.无论是地形测量还是植被估算,在数据处理时都需要将来自地表和植被的信号分离开,即滤波.数据滤波的方法很多,但是都不成熟,而且处于技术保密状态.本文介绍目前主要的信号分类方法,并开发了动态阈值分离法.该方法以森林为解析对象,通过反复运算逐渐逼近,完成地表和植被的分离,从而实现了地形与森林测量的双重目的.通过引入地形因子把原地形转换为理想平地,大大减少了过度滤波或者滤波不充分的情况.本方法基本实现自动处理,容易操作.利用本方法计算的高度误差(均方差,RMS)为0.4-0.5m.本方法还适用于城市地区高程模型的建立.     

Simulation of soil liquefaction distribution in downtown Mashiki during 2016 Kumamoto earthquake using nonlinear site response

Bulletin of Earthquake Engineering - Several sites located between Road No.28 and Akitsu River in downtown Mashiki were liquefied during the mainshock of the 2016 Kumamoto earthquake. According to...     

耦合融雪的分布式流域"降雨-径流"数值模型

针对中国对耦合融雪的"降雨-径流"过程数值计算模型方面的研究较少,本文联合应用能量平衡方程及运动波理论的基础方程式,结合GIS技术构建分布式流域耦合融雪的"降雨-径流"数值计算模型,并应用流域实测资料进行了有效性验证。研究结果表明:在不考虑融雪条件下,实测流量与计算结果之间存在着较大的误差;在考虑融雪条件下,实测流量与计算结果之间的误差在基准允许范围之内(小于3%);在积雪与融雪区,对流域"降雨-径流"过程的数值计算须考虑融雪。该研究为积雪融雪区提供了一种耦合融雪的"降雨-径流"过程数值计算的新方法。     

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.     

Magnetotelluric resistivity modeling for 3D characterization of geothermal reservoirs in the Western side of Mt. Aso, SW Japan

Geothermal reservoirs are usually located at a depth range of 2 to 5 km, so to efficiently utilize such resources an advanced prospecting method is needed to detect these deep geologic structures. This study aimed to three-dimensionally characterize geothermal reservoirs by a combination of Magnetotelluric (MT) survey, inversion analysis of apparent resistivity, and interpolation of the resistivity data obtained. The western side of Mt. Aso crater, southwest Japan, was chosen as the case study area. Three hot springs exist there and a fault is assumed to go in the direction connecting them. A MT survey was carried out at 26 sites and the data processed by a remote reference method to reduce artificial noises. Based on skewness and Mohr circle analyses of the impedance tensor, the local geologic structure at each site could be approximated as horizontally layered and therefore, a one-dimensional inversion analysis was applied to the MT raw data. The resultant resistivity column data were then interpolated by the three-dimensional optimization principle method. The resistivity distributions obtained clarified continuous conductors with especially low resistivity (less than 10 Ω·m) at the hot springs along the fault. Because the resistivity decreases largely with an abundance of clay minerals, the conductors could be considered to correspond with the cap rocks. Thus, two geothermal reservoirs, whose shapes were estimated to be pillar, were detected under the cap rocks in an elevation range of − 1000 to − 3000 m. By comparing the resistivity distributions with the temperature distributions based on fluid-flow calculations at a steady state using FEM, the validity of the location and dimension of the estimated reservoirs were confirmed.