Application of bivariate frequency analysis to the derivation of rainfall–frequency curves

Bivariate distributions have been recently employed in hydrologic frequency analysis to analyze the joint probabilistic characteristics of multivariate storm events. This study aims to derive practical solutions of application for the bivariate distribution to estimate design rainfalls corresponding to the desired return periods. Using the Gumbel mixed model, this study constructed rainfall–frequency curves at sample stations in Korea which provide joint relationships between amount, duration, and frequency of storm events. Based on comparisons and analyses of the rainfall–frequency curves derived from univariate and bivariate storm frequency analyses, this study found that conditional frequency analysis provides more appropriate estimates of design rainfalls as it more accurately represents the natural relationship between storm properties than the conventional univariate storm frequency analysis.     

A Method for Installing Piezometers in Large Cobble Bed Rivers

An impact drive point method is described for emplacing piezometers in a cobble river bottom where this has previously been difficult without the use of drilling rigs. To force the drive point piezometers through coble, the vibrational impact of an air-powered hammer was carried directly to the drive point by the use of an internal drive rod. After insertion to depth, the drive rod was removed from the lower portion of the piezometer and a standpipe was added to extend the piezometer above the river level. Piezometers installed in this way have permitted water quality analysis and dynamic measurement of vertical potentials in cobble sediments ranging in size from 2.5 to >30 cm and the method has been successfully used in the Columbia River, USA, and Töss River, Switzerland. This innovative method provides information on the hydrodynamics of pore water in highly permeable, cobble deposits that are common in high-energy river and lake bottoms. Piezometers installed using the internal drive rod method facilitate the assessment of the temporal and spatial dynamics of recharge and discharge at the ground water/surface water interface and analyses of the ecological connectivity between the hyporheic zone and surface water of rivers and streams. This information will lead to improved management decisions related to our nation's ground water and surface water supplies.     

Evaluation of riverbed magnetic susceptibility for mapping biogeochemical hot spots in groundwater-impacted rivers

Redox hot spots occurring as metal-rich anoxic groundwater discharges through oxic wetland and river sediments commonly result in the formation of iron (Fe) oxide precipitates. These redox-sensitive precipitates influence the release of nutrients and metals to surface water and can act as ‘contaminant sponges’ by absorbing toxic compounds. We explore the feasibility of a non-invasive, high-resolution magnetic susceptibility (MS) technique to efficiently map the spatial variations of magnetic Fe oxide precipitates in the shallow bed of three rivers impacted by anoxic groundwater discharge. Laboratory analyses on Mashpee River (MA, USA) sediments demonstrate the sensitivity of MS to sediment Fe concentrations. Field surveys in the Mashpee and Quashnet rivers (MA, USA) reveal several discrete high MS zones, which are associated with likely anoxic groundwater discharge as evaluated by riverbed temperature, vertical head gradient, and groundwater chemistry measurements. In the East River (CO, USA), widespread cobbles/rocks exhibit high background MS from geological ferrimagnetic minerals, thereby obscuring the relatively small enhancement of MS from groundwater induced Fe oxide precipitates. Our study suggests that, in settings with low geological sources of magnetic minerals such as lowland rivers and wetlands, MS may serve as a complementary tool to temperature methods for efficiently mapping Fe oxide accumulation zones due to anoxic groundwater discharges that may function as biogeochemical hot spots and water quality control points in gaining systems.     

High precision glacial–interglacial benthic foraminiferal Sr/Ca records from the eastern equatorial Atlantic Ocean and Caribbean Sea

Glacial–interglacial variation in the marine Sr/Ca ratio has important implications for coral Sr thermometry [J.W. Beck et al., Science 257 (1992) 644–647]. A possible variation of 1–3% was proposed based on ocean models [H.M. Stoll and D.P. Schrag, Geochim. Cosmochim. Acta 62 (1998) 1107–1118]. Subsequently, studies have used fossil foraminifera to test this prediction [P.A. Martin et al., Geochem. Geophys. Geosyst. 1 (1999); H.M. Stoll et al., Geochim. Cosmochim. Acta 63 (1999) 3535–3547; H. Elderfield et al., Geochem. Geophys. Geosyst. 1 (2000)]. But whether some component of foraminiferal Sr/Ca variation can be uniquely ascribed to seawater Sr variation is still not clear. To address this question, we developed cleaning and analysis techniques and measured Sr/Ca ratios on individual shells of the modern benthic foraminifer Cibicidoides wuellerstorfi. We showed that different size shells have different Sr/Ca ratios; however, samples with shell sizes of 355–500 μm appear to have normally distributed Sr/Ca ratios (1σ=1.8%). For multi-shell measurements (with estimated errors of 0.12–0.39%), the ratio varied by as much as 7.2±0.5% during the last glaciation for two Caribbean records at the same site and by 3.7±0.5% over the past 40,000 yr for one record from the Sierra Leone Rise in the eastern equatorial Atlantic. The two Caribbean records are very similar indicating that the behavior of shell Sr uptake was identical locally and that the shell Sr/Ca ratio faithfully reflects the local environment. The Atlantic record differs from the Caribbean records by as much as several percent. Thus, the foraminiferal Sr/Ca changes cannot be solely due to changes in seawater Sr/Ca unless the glacial deep ocean had spatial variation in Sr/Ca well in excess of the modern ocean. Certain similarities between the three records do exist. Notably, the rate of change of Sr/Ca is similar between 9 and 0 ka (−0.25%/kyr) and between 25 and 16 ka (+0.16%/kyr). This suggests that during these intervals, benthic foraminiferal Sr/Ca was affected by similar large-scale variables. One of these variables may be the average marine Sr/Ca ratio; however, comparison with model predictions [H.M. Stoll and D.P. Schrag, Geochim. Cosmochim. Acta 62 (1998) 1107–1118] suggests other factors must also be considered. The discrepancies between the two sites may be related to the different water mass histories for the Caribbean and eastern Atlantic. Our results suggest that variation of the seawater Sr budget only partially contributed to C. wuellerstorfi Sr/Ca records, while other significant factors still need to be quantified. At present we cannot confidently determine past seawater Sr/Ca variation from our foraminiferal records.     

Neuro‐control of seismically excited steel structure through sensitivity evaluation scheme

The neuro‐controller training algorithm based on cost function is applied to a multi‐degree‐of‐freedom system; and a sensitivity evaluation algorithm replacing the emulator neural network is proposed. In conventional methods, the emulator neural network is used to evaluate the sensitivity of structural response to the control signal. To use the emulator, it should be trained to predict the dynamic response of the structure. Much of the time is usually spent on training of the emulator. In the proposed algorithm, however, it takes only one sampling time to obtain the sensitivity. Therefore, training time for the emulator is eliminated. As a result, only one neural network is used for the neuro‐control system. In the numerical example, the three‐storey building structure with linear and non‐linear stiffness is controlled by the trained neural network. The actuator dynamics and control time delay are considered in the simulation. Numerical examples show that the proposed control algorithm is valid in structural control. Copyright © 2001 John Wiley & Sons, Ltd.     

Principal axes of m-DOF structures PartⅠ:Static loading

This paper is the first in a two-part series that discusses the principal axes of M-DOF structures subjected to static and dynamic loads. The primary purpose of this series is to understand the magnitude of the dynamic response of structures to enable better design of structures and control modification devices/systems. Under idealized design conditions, the structural responses are obtained by using single direction input ground motions in the direction of the intended control devices/systems,and by assuming that the responses of the structure is decoupleable in three mutually perpendicular directions. This standard practice has been applied to both new and retrofitted structures using various seismic protective systems. Very limited information is available on the effects of neglecting the impact of directional couplings (cross effects - of which torsion is a component) of the dynamic response of structures. In order to quantify such effects, it is necessary to examine the principal axes of structures under both static and dynamic loading.This first paper deals with quantitative definitions of principal axes and "cross effects" of three-dimensional structures under static load by using linear algebra. It shows theoretically that, for three-dimensional structures, such principal axes rarely exist. Under static loading conditions, the cross effect is typically small and negligible from the viewpoint of engineering applications. However, it provides the theoretical base for subsequent quantification of the response couplings under dynamic loads, which is reported in part Ⅱ of this series.     

Understanding of the Common Land Model performance for water and energy fluxes in a farmland during the growing season in Korea

The Common Land Model (CLM) is one of the most widely used land surface models (LSMs) due to the practicality of its simple parameterization scheme and its versatility in embracing a variety of field datasets. The improved assessment of land surface water and energy fluxes using CLM can be an alternative approach for understanding the complex land–atmosphere interactions in data‐limited regions. The understanding of water and energy cycles in a farmland is crucial because it is a dominant land feature in Korea and Asia. However, the applications of CLM to farmland in Korea are in paucity. The simulations of water and energy fluxes by CLM were conducted against those from the tower‐based measurements during the growing season of 2006 at the Haenam site (a farmland site) in Korea without optimization. According to the International Geosphere–Biosphere Programme (IGBP) land cover classification, a homogeneous cropland was selected initially for this study. Although the simulated soil moisture had a similar pattern to that of the observed, the former was relatively drier (at 0·1 m³ m^?3) than the latter. The simulated net radiation showed good agreement with the observed, with a root mean squared error (RMSE) of 41 W m^?2, whereas relatively large discrepancies between the simulation and observation were found in sensible (RMSE of 66 W m^?2) and latent (RMSE of 60 W m^?2) heat fluxes. On the basis of the sensitivity analysis, soil moisture was more receptive to land cover and soil texture parameterizations when compared to soil temperature and turbulent fluxes. Despite the uncertainty in the predictive capability of CLM employed without optimization, the initial performance of CLM suggests usefulness in a data‐limited heterogeneous farmland in Korea. Further studies are required to identify the controls on water and energy fluxes with an improved parameterization. Copyright © 2010 John Wiley & Sons, Ltd.     

Thermal evolution of the Tertiary Shimanto Belt, Muroto Peninsula, Shikoku, Japan

Abstract The Shimanto accretionary complex on the Muroto Peninsula of Shikoku comprises two major units of Tertiary strata: the Murotohanto Sub-belt (Eocene-Oligocene) and the Nabae Sub-belt (Oligocene-Miocene). Both sub-belts have been affected by thermal overprints following the peak of accretion-related deformation. Palaeotemperatures for the entire Tertiary section range from ～ 140 to 315°C, based upon mean vitrinite reflectance values of 0.9–5.0%R_m. Values of illite crystallinity index are consistent with conditions of advanced diagenesis and anchimetamorphism. Illite/mica b₀ lattice dimensions indicate that burial pressures were probably no greater than 2.5kbar. In general, levels of thermal maturity are higher for the Murotohanto Sub-belt than for the Nabae Sub-belt. The Eocene-Oligocene strata also display a spatial decrease in thermal maturity from south to north and this pattern probably was caused by regional-scale differential uplift following peak heating. Conversely, the palaeothermal structure within the Nabae Sub-belt is fairly uniform, except for the local effects of mafic intrusions at the tip of Cape Muroto. There is a paleotemperature difference of ～ 90°C across the boundary between the Murotohanto and Nabae Sub-belts (Shiina-Narashi fault), and this contrast is consistent with approximately 1200 m of post-metamorphic vertical offset. Subduction prior to Middle Miocene probably involved the Kula or fused Kula-Pacific plate and the background geothermal gradient during the Eocene-Oligocene phase of accretion was ～ 30–35°C/km. Rapid heating of the Shimanto Belt evidently occurred immediately after a Middle Miocene reorganization of the subduction boundary. Hot oceanic lithosphere from the Shikoku Basin first entered the subduction zone at ～ 15 Ma; this event also coincided with the opening of the Sea of Japan and the rapid clockwise rotation of southwest Japan. The background geothermal gradient at that time was ～ 70°C/km. Whether or not all portions of the inherited (Eocene-Oligocene) palaeothermal structure were overprinted during the Middle Miocene remains controversial.     

cDNA cloning and expression of a cytochrome P450 1A (CYP1A) gene from the hermaphroditic fish Rivulus marmoratus

We previously described the genomic structure of the cytochrome P450 1A (CYP1A) gene from the hermaphroditic fish Rivulus marmoratus [Kim, I.-C., Kim, Y.J., Yoon, Y.-D, Kawamura, S., Lee, Y.-S., Lee, J.-S., 2004a. Cloning of cytochrome P450 1A (CYP1A) genes from the hermaphroditic fish Rivulus marmoratus and the Japanese medaka Oryzias latipes. Mar. Environ. Res. 58, 125–129]. To further characterize R. marmoratus CYP1A, we cloned the cDNA sequence of a CYP1A gene from this species and also expressed its recombinant protein in an E. coli system. We exposed R. marmoratus to 4-nonylphenol, and found a small induction of CYP1A mRNA in the treated animals. In this paper, we discuss the characteristics of R. marmoratus CYP1A gene as well as its potential use in a biomonitoring assay.