Field calibration of a formula for entrance mixing of river inflows to lakes: Lake Taupo,North Island,New Zealand

Field measurements were used to validate predictions for the initial dilution of negatively buoyant, cold‐water inflows to Lake Taupo, as part of a study to quantify mixing processes associated with the two largest inflows to the lake. The predictions were made using a formulation originally derived for positively buoyant, warm‐water inflows to cooling ponds. The formulation predicts the total dilution of an inflow during its inertia‐dominated phase between its entrance to the lake and the point where buoyancy forces are great enough to cause the inflow to plunge and form a submerged density current. In one of the measured inflows, the inflowing jet was free to entrain from both sides; in the other, entrainment was restricted on one side by attachment of the inflowing jet to the shoreline of a bay just upstream of the plunge point. In the former example, the unmodified coefficients from the cooling pond formulation provided an excellent prediction of initial dilution. In the latter example, entrainment was reduced and different coefficients were derived. In both examples the inflows remained attached to the lake bed throughout their course until their liftoff at depths of 45–55 m to form interflows. The difference between coefficients for the two inflows indicates that the coefficient values should be considered site‐specific. The formulation is not valid for inflows that separate from the bottom of the inflow channel before plunging. The entrance mixing formulation was incorporated in a more general model of lake stratification, DYRESM, which already includes a well‐documented routine for routing underflows down submerged channels on the bed of a lake. Application of the model to the inflows measured in Lake Taupo gave good results for two model outputs that were not involved in the calibration of the entrance mixing formulation, but that are affected by the result of the initial dilution calculation—the temperatures in the river plume after it has plunged, and the insertion depth.     

Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

On the computation of the electrical potential inside a horizontally‐layered half‐space

A new formulation is proposed for the electrical potential developed inside a horizontally‐layered half‐space for a direct current point‐source at the surface. The recursion formula for the kernel coefficient in the potential integral is simpler than the generally used two‐coefficient recursion. The numerical difficulties that may occur during the computation of the integrals and near the source axis are examined and solutions are proposed. The set of equations permits a stable and accurate computation of the tabular potential everywhere in the medium.     

Man-induced hydrological changes, metazooplankton communities and invasive species in the Berre Lagoon (Mediterranean Sea, France)

The Berre Lagoon has been under strong anthropogenic pressure since the early 1950s. The opening of the hydroelectric EDF power plant in 1966 led to large salinity drops. The zooplankton community was mainly composed of two common brackish species: Acartia tonsa and Brachionus plicatilis. Since 2006, European litigation has strongly constrained the input of freshwater, maintaining the salinity above 15. A study was performed between 2008 and 2010 to evaluate how these modifications have impacted the zooplankton community. Our results show that the community is more diverse and contains several coastal marine species (i.e., Centropages typicus, Paracalanus parvus and Acartia clausi). A. tonsa is still present but is less abundant, whereas B. plicatilis has completely disappeared. Strong predatory marine species, such as chaetognaths, the large conspicuous autochtonous jellyfish Aurelia aurita and the invasive ctenophore Mnemiopsis leidyi, are now very common as either seasonal or permanent features of the lagoon.     

Nonlinear earthquake analysis of high arch dam–water–foundation rock systems

A nonlinear finite element model for earthquake response analysis of arch dam–water–foundation rock systems is proposed in this paper. The model includes dynamic dam–water and dam–foundation rock interactions, the opening of contraction joints, the radiation damping of semi‐unbounded foundation rock, the compressibility of impounded water, and the upstream energy propagating along the semi‐unbounded reservoir. Meanwhile, a new equivalent force scheme is suggested to achieve free‐field input in the model. The effects of the earthquake input mechanism, joint opening, water compressibility, and radiation damping on the earthquake response of the Ertan arch dam (240 m high) in China are investigated using the proposed model. The results show that these factors significantly affect the earthquake response of the Ertan arch dam. Such factors should therefore be considered in the earthquake response analysis and earthquake safety evaluation of high arch dams. Copyright © 2011 John Wiley & Sons, Ltd.     

Seismic isolation for low‐to‐medium‐rise buildings using granulated rubber–soil mixtures: numerical study

This paper presents the preliminary research works on a potential seismic isolation method that makes use of scrap rubber tires for the protection of low‐to‐medium‐rise buildings. The method involves mixing shredded rubber tire particles with soil materials and placing the mixtures around building foundations, which provides a function similar to that of a cushion. Meanwhile, the stockpiling of scrap tires is a significant threat to our environment, and the engineering community has been looking for long‐term viable solutions to the recycling and reuse of rubber. A finite element program has been developed for modeling the time‐domain dynamic responses of soil–foundation–structure system, by which the effectiveness and robustness of the proposed method have been evaluated. In general, the structural responses, in terms of acceleration and inter‐story drift, can be reduced by 40–60%. Copyright © 2012 John Wiley & Sons, Ltd.     

Fragility analysis of low‐rise masonry in‐filled reinforced concrete buildings by a coefficient‐based spectral acceleration method

This study presents a seismic fragility analysis of low‐rise masonry in‐filled (MI) reinforced concrete (RC) buildings using a proposed coefficient‐based spectral acceleration method. The coefficient‐based method, without requiring any complicated finite element analysis, is a simplified procedure for assessing the spectral acceleration demand (or capacity) of buildings subjected to earthquakes. This paper begins with a calibration of the proposed coefficient‐based method for low‐rise MI RC buildings using published experimental results obtained from shaking table tests. Spectral acceleration‐based fragility curves for low‐rise MI RC buildings under various inter‐story drift limits are then constructed using the calibrated coefficient‐based method. A comparison of the experimental and estimated results indicates that the simplified coefficient‐based method can provide good approximations of the spectral accelerations at peak loads of low‐rise MI RC buildings, if a proper set of drift‐related factors and initial fundamental periods of structures are used. Moreover, the fragility curves constructed using the coefficient‐based method can provide a satisfactory vulnerability evaluation for low‐rise MI RC buildings under a given performance level. Copyright © 2011 John Wiley & Sons, Ltd.     

Inelastic displacement demands in steel structures and their relationship with earthquake frequency content parameters

This paper deals with the estimation of peak inelastic displacements of SDOF systems, representative of typical steel structures, under constant relative strength scenarios. Mean inelastic deformation demands on bilinear systems (simulating moment resisting frames) are considered as the basis for comparative purposes. Additional SDOF models representing partially‐restrained and concentrically‐braced (CB) frames are introduced and employed to assess the influence of different force‐displacement relationships on peak inelastic displacement ratios. The studies presented in this paper illustrate that the ratio between the overall yield strength and the strength during pinching intervals is the main factor governing the inelastic deformations of partially‐restrained models and leading to significant differences when compared with predictions based on bilinear structures, especially in the short‐period range. It is also shown that the response of CB systems can differ significantly from other pinching models when subjected to low or moderate levels of seismic demand, highlighting the necessity of employing dedicated models for studying the response of CB structures. Particular attention is also given to the influence of a number of scalar parameters that characterise the frequency content of the ground motion on the estimated peak displacement ratios. The relative merits of using the average spectral period T_aver, mean period T_m, predominant period T_g, characteristic period T_c and smoothed spectral predominant period T_o of the earthquake ground motion, are assessed. This paper demonstrates that the predominant period, defined as the period at which the input energy is maximum throughout the period range, is the most suitable frequency content scalar parameter for reducing the variability in displacement estimations. Finally, noniterative equivalent linearisation expressions based on the secant period and equivalent damping ratios are presented and verified for the prediction of peak deformation demands in steel structures. Copyright © 2011 John Wiley & Sons, Ltd.     

Characteristics of concentrated flow hydraulics for rangeland ecosystems: implications for hydrologic modeling

Concentrated flow is often the dominant source of water erosion following disturbance on rangelands. Because of the lack of studies that explain the hydraulics of concentrated flow on rangelands, cropland‐based equations have typically been used for rangeland hydrology and erosion modeling, leading to less accurate predictions due to different soil and vegetation cover characteristics. This study investigates the hydraulics of concentrated flow using unconfined field experimental data over diverse rangeland landscapes within the Great Basin Region, United States. The results imply that the overall hydraulics of concentrated flow on rangelands differ significantly from those of cropland rills. Concentrated flow hydraulics on rangelands are largely controlled by the amount of cover or bare soil and hillslope angle. New predictive equations for concentrated flow velocity (R² = 0·47), hydraulic friction (R² = 0·52), and width (R² = 0·4) representing a diverse set of rangeland environments were developed. The resulting equations are applicable across a wide span of ecological sites, soils, slopes, and vegetation and ground cover conditions and can be used by physically‐based rangeland hydrology and erosion models to estimate rangeland concentrated flow hydraulic parameters. Published in 2011. This article is a US Government work and is in the public domain in the USA.