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.     

Crustal structure of the continental margin of Korea in the East Sea (Japan Sea) from deep seismic sounding data: evidence for rifting affected by the hotter than normal mantle

Despite the various opening models of the southwestern part of the East Sea (Japan Sea) between the Korean Peninsula and the Japan Arc, the continental margin of the Korean Peninsula remains unknown in crustal structure. As a result, continental rifting and subsequent seafloor spreading processes to explain the opening of the East Sea have not been adequately addressed. We investigated crustal and sedimentary velocity structures across the Korean margin into the adjacent Ulleung Basin from multichannel seismic (MCS) reflection and ocean bottom seismometer (OBS) data. The Ulleung Basin shows crustal velocity structure typical of oceanic although its crustal thickness of about 10 km is greater than normal. The continental margin documents rapid transition from continental to oceanic crust, exhibiting a remarkable decrease in crustal thickness accompanied by shallowing of Moho over a distance of about 50 km. The crustal model of the margin is characterized by a high-velocity (up to 7.4 km/s) lower crustal (HVLC) layer that is thicker than 10 km under the slope base and pinches out seawards. The HVLC layer is interpreted as magmatic underplating emplaced during continental rifting in response to high upper mantle temperature. The acoustic basement of the slope base shows an igneous stratigraphy developed by massive volcanic eruption. These features suggest that the evolution of the Korean margin can be explained by the processes occurring at volcanic rifted margins. Global earthquake tomography supports our interpretation by defining the abnormally hot upper mantle across the Korean margin and in the Ulleung Basin.     

Theory of satellite ground-track crossovers

The fundamental geometry of satellite ground tracks and their crossover problem are investigated. For idealized nominal ground tracks, the geometry is governed by a few constant parameters whose variations lead to qualitative changes in the crossover solutions. On the basis that the theory to locate crossovers has not been studied in sufficient detail, such changes are described in regard to the number of crossover solutions in conjunction with their bifurcations. Employing the spinor algebra as a tool for establishing the ground-track crossing condition, numerical methodologies to locate crossovers appearing in general dual-satellite ground-track configurations are also presented. The methodologies are applied to precisely determined orbital ephemerides of the GEOSAT, ERS-1, and TOPEX/POSEIDON altimeter satellites. Received: 19 November 1996 / Accepted: 12 May 1997     

Application of an optimum design technique for determining the coefficient of consolidation by using piezocone test data

For normally consolidated clay, several researchers have developed a number of theoretical time factors to determine the coefficient of consolidation from piezocone test results. However, depending on assumptions and analytical techniques, it could vary considerably, even for a specific degree of consolidation. In this paper a method is proposed to determine a consistent coefficient of consolidation by applying the concept of an optimum design technique over all ranges of the degree of consolidation. Initial excess pore pressure distribution is assumed to be capable of being obtained by the successive spherical cavity expansion theory. The dissipation of pore pressure is simulated by means of a two-dimensional linear-uncoupled axi-symmetric consolidation analysis. The minimization of differences between measured and predicted excess pore pressure was carried out by the BFGS unconstrained optimum design algorithm with a one-dimensional golden section search technique. By analyzing numerical examples and in-situ test results, it was found that the adopted optimum design technique gives consistent and convergent results.     

Regime interplay: the case of biodiversity and climate change

In this paper, the interplay between the biodiversity and climate change regimes is examined in the context of a case study. An understanding of how the regimes interact is crucial to their successful operation, since it provides insight into how policy options might be drawn up to create synergy rather than conflict among the respective regimes’ objectives. Despite the increasing importance of the study of regime interplay, due to its potential influence on regime effectiveness, there is limited understanding of why regimes interact and how their interplay may affect their successful operation, since few empirical studies have been conducted in this area. By analysing the causes of the regime interplay and the potential effects thereof, this paper draws policy implications regarding how the interacting regimes might be better co-ordinated to create synergy.     

Influence of the total atmospheric optical depth and cloud cover on solar irradiance components

Broadband solar irradiance data obtained in the spectral range 400–940 nm at Kwangju, South Korea from 1999–2000 have been analyzed to investigate the effects of cloud cover and atmospheric optical depth on solar radiation components. Results from measurements indicate that the percentage of direct and diffuse horizontal components of solar irradiance depend largely on total optical depth (TOD) and cloud cover. During summer and spring, the percentages of diffuse solar irradiance relative to the global irradiance were 5.0% and 4.9% as compared to 2.2% and 3.0% during winter and autumn. The diffuse solar irradiance is higher than the direct in spring and summer by 24.2%, and 40.6%, respectively, which may largely be attributed to the attenuation (scattering) of radiation by heavy dust pollution and large cloud amount. In cloud-free conditions with cloud cover ≤2/10, the fraction of the direct and diffuse components were 66.0% and 34.0%, respectively, with a mean daily global irradiance value of 7.92±2.91 MJ m⁻² day⁻¹. However, under cloudy conditions (with cloud cover ≥8/10), the diffuse and direct fractions were 97.9% and 2.2% of the global component, respectively. The annual mean TOD under cloudless conditions (cloud cover≤2/10) yields 0.74±0.33 and increased to as much as 3.15±0.67 under cloudy conditions with cloud amount ≥8/10. An empirical formula is derived for estimating the diffuse and direct components of horizontal solar irradiance by considering the total atmospheric optical depth (TOD). Results from statistical models are shown for the estimation of solar irradiance components as a function of TOD with sufficient accuracy as indicated by low standard error for each solar zenith angle (SZA).     

Similarities between strike-slip faults at different scales and a simple age determining method for active faults

Abstract Several differently scaled strike‐slip faults were examined. The faults shared many geometric features, such as secondary fractures and linkage structures (damage zones). Differences in fault style were not related to specific scale ranges. However, it was recognized that differences in style may occur in different tectonic settings (e.g. dilational/contractional relays or wall/linkage/tip zones), different locations along the master fault or different fault evolution stages. Fractal dimensions were compared for two faults (Gozo and San Andreas), which supports the idea of self‐similarity. Fractal dimensions for traces of faults and fractures of damage zones were higher (D ～1.35) than for the main fault traces (D ～1.005) because of increased complexity due to secondary faults and fractures. Based on the statistical analysis of another fault evolution study, single event movements in earthquake faults typically have a maximum earthquake slip : rupture length ratio of approximately 10^?4, although this has only been established for large earthquake faults because of limited data. Most geological faults have a much higher maximum cumulative displacement : fault length ratio; that is, approximately 10^?2 to 10^?1 (e.g. Gozo, ～10^?2; San Andreas, ～10^?1). The final cumulative displacement on a fault is produced by accumulation of slip along ruptures. Hence, using the available information from earthquake faults, such as earthquake slip, recurrence interval, maximum cumulative displacement and fault length, the approximate age of active faults can be estimated. The lower limit of estimated active fault age is expressed with maximum cumulative displacement, earthquake slip and recurrence interval as T ? (d_max /u) · I(M).     

On the Effects of Non-planar Geometry for Blind Thrust Faults on Strong Ground Motion

— We quantify the effects of complex fault geometry on low-frequency (<1 Hz) strong ground motion using numerical modeling of dynamic rupture. Our tests include the computation of synthetic seismograms for several simple rupture scenarios with planar and curved fault approximations of the 1994 Northridge earthquake. We use the boundary integral equation method (BIEM) to compute the dynamic rupture process, which includes the normal stress effects along the curved fault geometries. The wave propagation and computation of synthetic seismograms are modeled using a fourth-order finite-difference method (FDM). The near-field ground motion is significantly affected by the acceleration, deceleration and arrest of rupture due to the curvature of the faults, as well as the variation in directivity of the rupture. For example, a 6-km-long hanging-wall or footwall splay with a maximum offset of 1 km can change 1-Hz peak velocities by up to a factor of 2-3 near the fault. Our tests suggest that the differences in waveform are larger on the hanging wall compared to those on the footwall, although the differences in amplitude are larger in the forward rupture direction (footwall). The results imply that kinematic ground motion estimates may be biased by the omission of dynamic rupture effects and even relatively gentle variation in fault geometry, and even for long-period waves.     

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.