Geostatistical interpolation of space–time rainfall on Tamshui River basin,Taiwan

Taiwan suffers from heavy storm rainfall during the typhoon season. This usually causes large river runoff, overland flow, erosion, landslides, debris flows, loss of power, etc. In order to evaluate storm impacts on the downstream basin, a real‐time hydrological modelling is used to estimate potential hazard areas. This can be used as a decision‐support system for the Emergency Response Center, National Fire Agency Ministry, to make ‘real‐time’ responses and minimize possible damage to human life and property. This study used 34 observed events from 14 telemetered rain‐gauges in the Tamshui River basin, Taiwan, to study the spatial–temporal characteristics of typhoon rainfall. In the study, regionalized theory and cross‐semi‐variograms were used to identify the spatial‐temporal structure of typhoon rainfall. The power form and parameters of the cross‐semi‐variogram were derived through analysis of the observed data. In the end, cross‐validation was used to evaluate the performance of the interpolated rainfall on the river basin. The results show the derived rainfall interpolator represents the observed events well, which indicates the rainfall interpolator can be used as a spatial‐temporal rainfall input for real‐time hydrological modelling. Copyright © 2007 John Wiley & Sons, Ltd.     

Physical factors influencing immature-fish communities in the surf zones of sandy beaches in northwestern Kyushu Island,Japan

We aim to understand the relationships between physical conditions and characteristics of the immature-fish community in surf zones of sandy beaches. Therefore, we obtained fish samples between March 2007 and February 2008 and analyzed certain physical conditions in the surf zones of 21 sandy beaches on the coastline of the northwestern Kyushu Island, Japan. We collected a total of 83 species and 6458 immature individuals. In a BIO-ENV analysis, the highest correlation was observed between fish assemblage and S20 (i.e., the slope from the shoreline to the sites where the depth was 20 m) and current velocity (CV) values. Stepwise multiple linear regression analyses revealed that the number of species and individuals decrease with an increase in the S20 and CV values. These results show that species richness and the abundance of immature-fish increase under shelving and calm conditions. Thus, immature-fish assemblages are strongly influenced by the prevailing physical conditions. Moreover, in six of the 10 dominant species, a negative correlation was observed between CV and abundance. On the other hand, S20 was found to be the explanatory variable only in the case of the most dominant species, i.e., Gymnogobius breunigii. Furthermore, a positive correlation was observed between S1 (i.e., the slope from the shoreline to the sites where the depth was 1.0 m at the mean tidal level) and median particle size (i.e., MPS of the sediments) and the abundances of Sillago japonica and Favonigobius gymnauchen, respectively, and a negative correlation with salinity, in the case of Acanthogobius lactipes. We conclude that the characteristics of the fish community in surf zones on sandy beaches are determined by not only the shelving and calm conditions, which influence fish assemblages and abundances, but also the habitat diversity, which influences the diversity of fish species.     

Comparison of the structural,spectroscopic and phenanthrene binding characteristics of humic acids from soils and lake sediments

The structural, spectroscopic and phenanthrene binding characteristics were compared for humic acids (HA) extracted from two different sources: eight soils and six lake sediments. The elemental analysis revealed that HA from sediments had higher H/C, N/C, and (N + O)/C ratios compared to HA from soils, reflecting a lower degree of humification and more autochthonous organic input upon the formation of the HA for lake sediments versus soil environments. HA from soils exhibited a higher content of aromatic carbon structures than the sediment HA based on ¹³C NMR results. The source of HA was easily distinguished by comparing the synchronous fluorescence spectra of each HA group. The presence of a protein like fluorescence was prominent for the HA from sediment while it was minor for the HA from soil. Irrespective of the HA source, however, humification index (HIX) exhibited a common positive correlation with the aromatic content, and a negative correlation with O-alkyl carbon structures of the HA. The correlations were consistent with the general structural trends of humification processes, suggesting that HIX may serve as a source independent predictor to describe the structural information and humification degree of terrigenous HA. Aromatic carbon structures did contribute to enhancing the phenanthrene binding for both sources of HA. However, the primary structures associated with non-ideal phenanthrene binding (i.e., non-linear sorption isotherm) appear to differ by the HA source because the opposed correlations were obtained between aromaticity and the isotherm linearity for the two HA groups. Our results suggest that the HA structural function associated with specific non-linear sorption for hydrophobic organic contaminants (HOCs) may be more strongly governed by the HA source than by the apparent physico-chemical properties.     

Study of the initiation and propagation of excavation damaged zones around openings in argillaceous rock

The study of the creation and evolution of the excavation disturbed zone (EDZ) in argillaceous rocks is a major issue for the safety of nuclear wastes underground repositories. In this context, the argillaceous Tournemire site has provided a unique opportunity to study the evolution of the EDZ with time thanks to the existence of three openings of different ages. A thorough characterization of the EDZ has been conducted by different means such as visual observation, analysis of samples extracted from drilled boreholes, EDZ permeability measurements, etc. On the basis of these measurements, a conceptual model of the EDZ initiation and propagation at the Tournemire site has been proposed. In order to validate this model, numerical simulations of increasing complexity have been carried out. In a first attempt, the response of the rock mass to the excavation phase, followed by seasonal cyclic variations of temperature and relative humidity inside the opening, has been simulated by means of a purely mechanical analysis, using a simple elastic material model. The EDZ has been estimated by post-processing the calculated stress states, using a Mohr–Coulomb failure criterion. The results obtained show that no EDZ could be predicted unless adopting a low cohesion value for the rock mass. Moreover, the deferred nature of the EDZ formation in Tournemire could not be reproduced. These limitations have then been suppressed by using a coupled viscoplastic-damaging mechanical model, the parameters of which have been identified from different laboratory experiments. With this model, a time evolution of the EDZ could be predicted, but the EDZ pattern could not match the one observed in situ. Finally, in view of the importance of the hydraulic couplings, unsaturated hydro-mechanical calculations have been carried out to investigate the effect of the numerous seasonal variations cycles and the resulting shrinkage.     

利用单频双程波动方程计算初至走时及其振幅

通过在频率域双程波动方程模拟算法中加入一个复数频率（实部表示频率,虚部表示衰减因子）压制地震波初至走时之后的能量,从而把初至走时及其振幅的计算问题转换为单一频率波场中最大能量走时和振幅的拾取问题,然后利用单一频率域波场的相位项和振幅项分别计算初至走时及其振幅.本文还提出利用参数分析方法求取最优的复数频率,并给出数值计算例子,将本方法的计算结果与有限差分程函方程初至走时和最大能量走时振幅进行比较,结果表明,该方法具有适应于任意复杂介质和多炮多接收点走时和振幅的计算.     

Two‐dimensional waveform inversion of multi‐component data in acoustic‐elastic coupled media

In order to account for the effects of elastic wave propagation in marine seismic data, we develop a waveform inversion algorithm for acoustic‐elastic media based on a frequency‐domain finite‐element modelling technique. In our algorithm we minimize residuals using the conjugate gradient method, which back‐propagates the errors using reverse time migration without directly computing the partial derivative wavefields. Unlike a purely acoustic or purely elastic inversion algorithm, the Green's function matrix for our acoustic‐elastic algorithm is asymmetric. We are nonetheless able to achieve computational efficiency using modern numerical methods. Numerical examples show that our coupled inversion algorithm produces better velocity models than a purely acoustic inversion algorithm in a wide variety of cases, including both single‐ and multi‐component data and low‐cut filtered data. We also show that our algorithm performs at least equally well on real field data gathered in the Korean continental shelf.     

Comparison of scaling methods for waveform inversion

Waveform inversion can lead to faint images for later times due to geometrical spreading. The proper scaling of the steepest-descent direction can enhance faint images in waveform inversion results. We compare the effects of different scaling techniques in waveform inversion algorithms using the steepest-descent method. For the scaling method we use the diagonal of the pseudo-Hessian matrix, which can be applied in two different ways. One is to scale the steepest-descent direction at each frequency independently. The other is to scale the steepest-descent direction summed over the entire frequency band. The first method equalizes the steepest-descent directions at different frequencies and minimizes the effects of the band-limited source spectrum in waveform inversion. In the second method, since the steepest-descent direction summed over the entire frequency band is divided by the diagonal of the pseudo-Hessian matrix summed over the entire frequency band, the band-limited property of the source wavelet spectrum still remains in the scaled steepest-descent directions. The two scaling methods were applied to both standard and logarithmic waveform inversion. For standard waveform inversion, the method that scales the steepest-descent direction at every frequency step gives better results than the second method. On the other hand, logarithmic waveform inversion is not sensitive to the scaling method, because taking the logarithm of wavefields automatically means that results for the steepest-descent direction at each frequency are commensurate with each other. If once the steepest-descent directions are equalized by taking the logarithm of wavefields in logarithmic waveform inversion, the additional equalizing effects by the scaling method are not as great as in conventional waveform inversion.     

Organohalogen contaminants in striped dolphins (Stenella coeruleoalba) from Japan: Present contamination status, body distribution and temporal trends (1978-2003)

Organohalogen contaminants including PCBs, DDTs, CHLs, HCHs, HCB, PBDEs and HBCDs were determined in striped dolphins (Stenella coeruleoalba) found stranded at Gogo-shima (n = 6, 2003) and collected from Taiji (n = 15, 1978-1992) in Japan. All target compounds were significantly detected in all the specimens, indicating ubiquitous contamination of oceanic cetaceans in northwest Pacific Ocean. Examination of body distribution of organohalogens in the six specimens from Gogo-shima showed no significant difference in concentrations among the analyzed tissues, except for brain, which had lower levels possibly due to the existence of blood-brain barrier. For evaluating temporal trends, archived blubber samples of adult male stripped dolphins collected in 1978, 1979, 1986 and 1992 were analyzed. Concentrations of PCBs, DDTs and HCHs did not change significantly during 1978-2003. In contrast, remarkable increasing trends of PBDEs and HBCDs were observed, suggesting growing consumption in Japan and surrounding countries in recent years.     

Field observations on correlation of fatty acid profiles between suspended particulate matter and green-lipped mussels in subtropical waters of Hong Kong

By analysis of the fatty acid profiles in mussel tissues and suspended particulate matter (SPM) in water, the present study showed a significant relationship of the trophic linkage between mussels and the SPM. At seven locations from inner to outer areas along the eutrophic Tolo Harbour and Tolo Channel, Hong Kong, the composition (as percentage of total fatty acids) of both monoenoic and polyunsaturated fatty acids in the tissues of green-lipped mussels Perna viridis and SPM in water had significant correlation (p<0.01). In particular, the composition of docosahexaenoic acid (DHA, 22:6n3) (as percentage of total fatty acids) in mussel tissues was statistically correlated with that in SPM (p<0.01), implying that mussels incorporate DHA, the biomarker of dinoflagellates. Principal component analysis further demonstrated that the fatty acid profiles of SPM were different among locations in the harbour, mid and outer channel of the study area, so were the mussel tissue fatty acid profiles. Cluster analysis of phytoplankton data also revealed the dominance of diatoms in the inner harbour and channel areas, whereas dinoflagellates were abundant in the outer channel waters. The possible implications of using benthic suspension feeders such as green-lipped mussels P. viridis for controlling phytoplankton abundance in coastal waters are discussed.     

Spatio-temporal changes of marine macrobenthic community in sub-tropical waters upon recovery from eutrophication. II. Life-history traits and feeding guilds of polychaete community

A two-year study was conducted in the vicinity of a harbour in sub-tropical Hong Kong, to examine the progress of recovery of macrobenthic community, based on analyses of both life-history traits and trophic guilds of polychaetes, upon cessation of organic pollution caused by sewage discharge. Seventy seven out of 83 species collected were classified under four ecological groups based on the life-history traits and sensitivity to organic gradients. The mean ATZI marine biotic index (AMBI) derived from these ecological groups showed spatial difference among the five sampling locations. In particular, the presence of different percentages of polychaete species from Groups III (unbalanced community) and IV (polluted community) suggested the presence of pollution stress in certain degree at all sampling locations. However, no significant temporal changes were noted over the study period. From all polychaete species identified, they were classified into 13 feeding guilds. The mean diversity of these feeding guilds at most of the sampling locations was significantly higher than that at one of the inside-harbour locations. The composition of feeding guilds was also significantly different spatially. At one of the inside-harbour locations, the dominant feeding guilds were motile/discretely motile surface deposit feeders with tentaculates or unarmed pharynx, and motile omnivores with jawed pharynx in the first year of study, but were replaced by motile, jawed carnivores in the second year of study. The increased proportion of carnivores over the study period can be seen as a sign of recovery in the community structure since abundance of predators is commonly higher in habitats with better environmental quality. The implications of using life-history traits and feeding guild analyses for benthic community are discussed.