System identification of the Vincent Thomas suspension bridge using earthquake records

The Vincent Thomas Bridge in the Los Angeles metropolitan area, is a critical artery for commercial traffic flow in and out of the Los Angeles Harbor, and is at risk in the seismically active Southern California region, particularly because it straddles the Palos Verdes fault zone. A combination of linear and non‐linear system identification techniques is employed to obtain a complete reduced‐order, multi‐input–multi‐output (MIMO) dynamic model of the Vincent Thomas Bridge based on the dynamic response of the structure to the 1987 Whittier and 1994 Northridge earthquakes. Starting with the available acceleration measurements (which consists of 15 accelerometers on the bridge structure and 10 accelerometers at various locations on its base), an efficient least‐squares‐based time‐domain identification procedure is applied to the data set to develop a reduced‐order, equivalent linear, multi‐degree‐of‐freedom model. Although not the main focus of this study, the linear system identification method is also combined with a non‐parametric identification technique, to generate a reduced‐order non‐linear mathematical model suitable for use in subsequent studies to predict, with good fidelity, the total response of the bridge under arbitrary dynamic environments. Results of this study yield measurements of the equivalent linear modal properties (frequencies, mode shapes and non‐proportional damping) as well as quantitative measures of the extent and nature of non‐linear interaction forces arising from strong ground shaking. It is shown that, for the particular subset of observations used in the identification procedure, the apparent non‐linearities in the system restoring forces are quite significant, and they contribute substantially to the improved fidelity of the model. Also shown is the potential of the identification technique under discussion to detect slight changes in the structure's influence coefficients, which may be indicators of damage and degradation in the structure being monitored. Difficulties associated with accurately estimating damping for lightly damped long‐span structures from their earthquake response are discussed. The technical issues raised in this paper indicate the need for added spatial resolution in sensor instrumentation to obtain identified mathematical models of structural systems with the broadest range of validity. Copyright © 2003 John Wiley & Sons, Ltd.     

Data Assimilation (4D-VAR) for Shallow-Water Flow: The Case of the Chicoutimi River

This is a web presentation of the work presented at the 10th Annual Conference of the CFD Society of Canada, “CFD 2002”, at the University of Windsor on June 9-11, 2002. This discussion paper presents the four-dimensional variational data assimilation (4D-VAR) technique as a tool to forecast floods. This discussion will be limited to hydrological forecast. We assume that the weather, here a large rainstorm, had already been forecasted by the meteorological services. In the 4D-VAR technique, we need to minimize, in the sense of Lagrange, a cost function which measures the difference between the forecast and the observations. The physical equations acts as a set of constraints. Here, the model is the shallow-water equations modified to include sediment transport. The minimum was found by using the steepest descent algorithm. This is made possible because the gradient of the cost function can be calculated analytically by using the adjoint equations of the model. To illustrate the 4D-VAR technique, the bypass of a simple theoretical dam as well as the more complex overflowing of the Chicoutimi River at the Chute-Garneau dam (during the 1996 flood) are investigated.     

Osmium-strontium-neodymium-lead isotopic covariations in mid-ocean ridge basalt glasses and the heterogeneity of the upper mantle

Osmium, strontium, neodymium, and lead isotopic data have been obtained for 30 hand picked samples of basaltic glass from the Pacific, Atlantic and Indian mid-oceanic ridges. Large variations in Os isotopic ratios exist in the glasses, from abyssal peridotite-like values to radiogenic compositions similar to oceanic island basalts (¹⁸⁷Os/¹⁸⁶Os and ¹⁸⁷Os/¹⁸⁸Os ratios range from 1.06 to 1.36 and from 0.128 to 0.163, respectively). Os isotopic and elemental data suggest the existence of mixing correlations. This relationship might be ascribed to secondary contamination processes; however, such a hypothesis cannot account for the negative correlation observed between Os and Nd isotopes and the existence of complementary covariations between Os and SrPb isotopes. In this case, OsSrNdPb isotopic variations are unrelated to late post-eruption or shallow level contamination. These relationships provide strong evidence that the Os isotopic composition of the samples are derived from the mantle and thus implies a global chemical heterogeneity of the oceanic upper mantle. The results are consistent with the presence of recycled oceanic crust in the mantle sources of mid-ocean ridge basalts, and indicate that the unique composition of the upper mantle below the Indian ocean results from its contamination by a mantle component characterized by radiogenic Os and particularly unradiogenic Nd and Pb isotopic compositions.     

Crystal morphologies in pillow basalts: implications for mid-ocean ridge processes

We present the results of a detailed petrological study of a sparsely phyric basalt (MAPCO CH98-DR11) dredged along the Mid-Atlantic Ridge (30°41′N). The sample contains microphenocrysts of olivine that display four different rapid-growth morphologies. Comparison of these morphologies with those obtained in dynamic crystallization experiments allows us to constrain the thermal history of the sample. The dendritic morphology (swallowtail, chain and lattice olivine) is directly related to the final quenching during magma–seawater interaction. In contrast, the three other morphologies, namely the complex polyhedral crystal, the closed hopper and the complex swallowtail morphology result from several cycles of cooling–heating (corresponding to a maximum degree of undercooling of 20–25°C) during crystal growth. These thermal variations occurred before eruption and are interpreted to be the result of turbulent convection in a small magmatic body beneath the ridge. The results suggest that the Mid-Atlantic Ridge is underlain by a mush zone that releases batches of liquid during tectonic segregation. Aphyric basalts are emitted during eruptions controlled by the tectonic activity, whereas phyric basalts correspond to small fractions of magma from the mush zone mobilized by reinjections of primitive magmas.     

Pedogenic factors affecting magnetic susceptibility of the last interglacial palaeosol S1 in the Chinese Loess Plateau

The magnetic susceptibility has been used as a quantitative or semi‐quantitative proxy for reconstructing the summer monsoon intensity in the Chinese Loess Plateau based on extensive studies on climatic or/and environmental mechanisms producing the magnetic susceptibility signatures. However, the precise nature of the link between past climates and the susceptibility signatures has remained uncertain primarily due to lack of our understanding in the ?nalizing and preserving processes of the signatures. This paper attempts to examine the reliability or acceptability of this summer monsoon proxy from non‐magnetic perspectives of soil‐forming processes. We chose nine sections along two transects: one across the western part of the Chinese Loess Plateau and another across the eastern part. Several conclusions can be drawn from our analytical data. First, clay translocation within the S1 palaeosol pro?les, as indicated by ?eld‐observed clay coatings on ped faces in Bt and Bk horizons and demonstrated by laboratory‐analysed clay contents, must have moved some of the magnetic minerals downward so that the susceptibility re?ects only the post‐translocation distribution of the magnetic‐susceptibility‐producing minerals. Second, the best‐developed palaeosol S1S3 at most of the sections studied is not expressed by the magnetic susceptibility because this palaeosol developed in underlying coarse loess (L2) and coarse textures tend to lower the susceptibility. Third, carbonate concentration is normally negatively correlated with the magnetic susceptibility or simply suppresses the magnetic susceptibility peak when the susceptibility enhancement exceeds the carbonate dilution effect. To conclude, extreme caution must be observed when using magnetic susceptibility signatures to retrieve high‐resolution records of the last interglacial palaeoclimate in the Chinese Loess Plateau. Copyright © 2004 John Wiley & Sons, Ltd.     

Optimal design of viscoelastic dampers using eigenvalue assignment

In this study a procedure for determining the optimum size and location of viscoelastic dampers is proposed using the eigenvalue assignment technique. Natural frequencies and modal damping ratios, required to realize a given target response, are determined first by the convex model. Then the desired dynamic structural properties are realized by optimally distributing the damping and stiffness coefficients of viscoelastic dampers using non‐linear programming based on the gradient of eigenvalues. This optimization method provides information on the optimal location as well as the magnitude of the damper parameters. The proposed procedure is applied to the retrofit of a 10‐story shear frame, and to a three‐dimensional structure with an asymmetric plan. The analysis results confirm that the responses of model structures retrofitted by the proposed method correspond well with the given target response. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic structural control using an electric servomotor active mass driver system

This study investigates an electric‐type active mass driver (AMD) system for structural vibration control. Composed primarily of an electric servomotor and a ball screw, the electrical AMD system is free from noise problems, oil leakage, and labor‐intensive maintenance that commonly are associated with hydraulic AMD systems. The desired stroke amplification of the mass and the power demand of the servomotor can be adjusted via the ball screw pitch, which in turn affects the effectiveness and efficiency of the system. Meanwhile, an instantaneous optimal direct output feedback control algorithm is adopted. Numerical simulation is performed using a five‐story steel frame as the object structure under the conditions of the 1940 El Centro earthquake. The AMD system proves to be effective and efficient within a certain range of the ball screw pitch. The reductions of the peak responses can reach as high as 70% if properly designed. Requiring only the velocity measurement of the top floor for on‐line feedback control, the proposed control algorithm is recommended for practical implementation. Copyright © 2004 John Wiley & Sons, Ltd.     

Magnetic fabrics and rock magnetism of Proterozoic dike swarm from the southern São Francisco Craton, Minas Gerais State, Brazil

Magnetic fabric and rock magnetism studies were performed on 32 mafic dikes of a Proterozoic dike swarm from the southern São Francisco Craton (SFC; Minas Gerais State, SE Brazil). Magnetic anisotropies were determined by applying anisotropy of low-field magnetic susceptibility (AMS) and anisotropy of remanent magnetization (ARM). The latter was performed imposing both anhysteretic (total (AAR) and partial pAAR)) and isothermal remanence magnetizations (AIRM). Partial anhysteretic remanence anisotropy was performed based on remanent coercivity spectra from a pilot specimen of each site. In most sites, AMS is dominantly carried by ferromagnetic minerals, however, in some sites, the paramagnetic contribution exceeds 70% of bulk susceptibility. Rock magnetism and thin section analysis allow classifying the dikes as non-hydrothermalized and hydrothermalized. Magnetic measurement shows that the mean magnetic susceptibility is usually lower than 5×10⁻³ (SI). Ti-poor titanomagnetites up to pure magnetite pseudo-single-domain (PSD) grain sizes carry the majority of magnetic fabrics for non-hydrothermalized dikes whereas coarse to fine grained Ti-poor titanomagnetites carry the majority of magnetic fabrics for hydrothermalized dikes.Three primary AMS fabrics are recognized which are coaxial with ARM fabric, except for two dikes, from both non-hydrothermalized and hydrothermalized dikes. Normal AMS fabric surprisingly is not dominant (31%). The parallelism between AMS, pAAR_0–30, pAAR_30–60 and pAAR_60–90 fabrics in the hydrothermalized dikes indicates that magnetic grains formed due to late-stage crystallization or to remobilization of iron oxides due to hydrothermal alteration after dike emplacement have acquired a mimetic fabric coaxial with the primary fabric given by coarse-grained early crystallized Ti-poor titanomagnetites. This fabric is interpreted as magma flow in which the analysis of K_max inclination permitted the inference that the dikes were fed by horizontal or subhorizontal fluxes (K_max<30°). Intermediate AMS fabric is the most important (41%) in the investigated swarm. It is interpreted as due to vertical compaction of a static magma column with the minimum stress along the dike strike. ARM determinations for these sites also remained intermediate except for two dikes. In one of them, AIRM fabric resulted in normal AMS fabric while for the other AAR fabric resulted in inverse AMS fabric. A combination of AMS and ARM fabrics suggest that magmatic fabric for both dikes were overprinted by some late local event, probably related to Brasiliano orogenic processes after dike emplacement. InverseInverse AMS fabric is a minority (four dikes). ARM determinations also remained inverse suggesting a primary origin for inverse AMS fabric.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.