GOCE Gravity Field Processing Strategy

A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.     

The organisation of fish assemblages in the regulated Lima basin, Northern Portugal

In order to understand the structure of fish assemblages in the modified Lima basin (Northern Portugal), two distinct datasets concerning the presence and abundance of fish species were subjected to multivariate analysis. On the River Lima two types of flow modification are present within kilometres of one another: (a) a reduced and constant flow due to hypolimnetic release; and (b) an intense and irregular flow. A comparison of their influence on fish assemblages revealed a gradient of assemblage types from tributaries to main river sites. The latter were characterised by a strong dominance of cyprinids, particularly Iberian barbel (Barbus bocagei). The former harboured two kinds of fish assemblages: those closer to the river mouth were dominated by the cyprinids Iberian chub (Squalius carolitertii) and Iberian nase (Chondrostoma polylepis), which were also frequently present in the main river; while in those further upstream the predominant species was the brown trout (Salmo trutta). Although explanatory variables such as distance from source, altitude, substrate coarseness and width were the primary correlates of fish assemblage composition, dam construction and flow regulation also had a significant effect upon assemblage structure, particularly by: i) reducing the importance of migratory species; ii) constraining the presence of trout in the regulated segments; and iii) simplifying the community, especially in the case of the constant and reduced flow regime.     

Controls on modern alluvial fan processes in the central Alps,northern Italy

Alluvial fan development in Alpine areas is often affected by catastrophic sedimentary processes associated with extreme ?oods events, causing serious risks for people living on the fans. Hazard assessment in these areas depends on proper identi?cation of the dominant sedimentary processes on the fans. Data from a set of 209 alluvial fans from the central Alps of Italy are presented in this paper and analysed with the help of various statistical techniques (linear regression, principal components analysis, cluster analysis, discriminant analysis and logistic regression). First, we used modern sedimentary facies and historical records (?ood events since 15th century), to distinguish between the two dominant sedimentary processes on alluvial fans: debris ?ows and stream?ows. Then, in order to analyse the main controls on past and present fan processes, 36 morphological, geological and land‐use variables were analysed. As with observations for arid‐environment fans, catchment morphology is the most in?uential factor in the study area, whereas geology and land use are minor controls. The role of climatic change and landsliding within the catchments also seems to be very important and is discussed. Statistical techniques also help in differentiating groups of alluvial fans by sets of controlling factors, including stage and type of evolution. Finally, by using discriminant analysis and logistic regression, we classi?ed alluvial fans according to the dominant sedimentary process, with a success rate ranging between 75 and 92 per cent. Copyright © 2004 John Wiley & Sons, Ltd.     

Estimating the seismic displacement of friction pendulum isolators based on non‐linear response history analysis

A procedure for developing equations that estimate the isolator displacement due to strong ground motion is applied to buildings isolated with the friction pendulum system. The resulting design equations, based on rigorous non‐linear analysis, offer an alternative to the iterative equivalent‐linear methods used by current U.S. building codes. The governing equations of the system are reduced to a form such that the median normalized displacement of the system due to an ensemble of ground motions is found to depend on only the isolation period—a function of the curvature of the isolator—and the friction force at incipient slip normalized by peak ground velocity. The normalization is effective in minimizing the dispersion of the normalized displacement for an ensemble of ground motions, implying that the median normalized displacement is a reliable estimate of response. The design equations reflect the significant (20 to 38%) increase in displacement when the excitation includes two lateral components of ground motion instead of just one component. Equivalent‐linear methods are shown to underestimate by up to 30% the exact median displacement determined by non‐linear response history analysis for one component of ground motion, and building codes include at most a 4.4% increase for a second component. Copyright © 2003 John Wiley & Sons, Ltd.     

Strain field measurements of rubber by image analysis and design criteria for laminated rubber bearings (LRB)

Although seismic isolation rubber bearings in bridges and buildings have proven to be a very effective passive method for reducing earthquake‐induced forces, a detailed mechanical modeling of the rubber that is used in bearings under large strains has not been established. Therefore, a 3D model of failure behavior and the design criteria for the safety evaluation of seismic isolation bearings have not yet been developed. This paper presents: (1) correlation‐based template‐matching algorithms to measure large strain fields of continua; (2) a failure criterion for rubber; and (3) the design criteria for the safety evaluation of laminated algorithms, data‐validation algorithms were developed and implemented to eliminate possible unrealistic displacement vectors present in the measured displacement field. The algorithms were successfully employed in the strain field measurement of LRB and rubber materials that are subjected to failure. The measured local strains for rubber material at failure were used to develop a failure criterion for rubber. The validity of the proposed criterion was evaluated by applying it to the LRB; the criterion was introduced into a 3D finite element model of LRB, compared with the experimental results of bearings failure, and verified. Finally, design criteria are proposed for LRB for the safety evaluation. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic response of a cable‐stayed bridge deck under multi‐component non‐stationary random ground motion

A Markov method of analysis is presented for obtaining the seismic response of cable‐stayed bridges to non‐stationary random ground motion. A uniformly modulated non‐stationary model of the random ground motion is assumed which is specified by the evolutionary r.m.s. ground acceleration. Both vertical and horizontal components of the motion are considered to act simultaneously at the bridge supports. The analysis duly takes into account the angle of incidence of the earthquake, the spatial correlation of ground motion and the quasi‐static excitation. A cable‐stayed bridge is analysed under a set of parametric variations in order to study the non‐stationary response of the bridge. The results of the numerical study indicate that (i) frequency domain spectral analysis with peak r.m.s. acceleration as input could provide more r.m.s. response than the peak r.m.s. response obtained by the non‐stationary analysis; (ii) the longitudinal component of the ground motion significantly influences the vertical vibration of the bridge; and (iii) the angle of incidence of the earthquake has considerable influence on the deck response. Copyright © 2003 John Wiley & Sons, Ltd.     

Frequency domain modal analysis of earthquake input energy to highly damped passive control structures

A new complex modal analysis‐based method is developed in the frequency domain for efficient computation of the earthquake input energy to a highly damped linear elastic passive control structure. The input energy to the structure during an earthquake is an important measure of seismic demand. Because of generality and applicability to non‐linear structures, the earthquake input energy has usually been computed in the time domain. It is shown here that the formulation of the earthquake input energy in the frequency domain is essential for deriving a bound on the earthquake input energy for a class of ground motions and for understanding the robustness of passively controlled structures to disturbances with various frequency contents. From the viewpoint of computational efficiency, a modal analysis‐based method is developed. The importance of overdamped modes in the energy computation of specific non‐proportionally damped models is demonstrated by comparing the energy transfer functions and the displacement transfer functions. Through numerical examinations for four recorded ground motions, it is shown that the modal analysis‐based method in the frequency domain is very efficient in the computation of the earthquake input energy. Copyright © 2004 John Wiley & Sons, Ltd.     

A modal pushover analysis procedure to estimate seismic demands for unsymmetric‐plan buildings

An Erratum has been published for this article in Earthquake Engng. Struct. Dyn. 2004; 33:1429. Based on structural dynamics theory, the modal pushover analysis (MPA) procedure retains the conceptual simplicity of current procedures with invariant force distribution, now common in structural engineering practice. The MPA procedure for estimating seismic demands is extended to unsymmetric‐plan buildings. In the MPA procedure, the seismic demand due to individual terms in the modal expansion of the effective earthquake forces is determined by non‐linear static analysis using the inertia force distribution for each mode, which for unsymmetric buildings includes two lateral forces and torque at each floor level. These ‘modal’ demands due to the first few terms of the modal expansion are then combined by the CQC rule to obtain an estimate of the total seismic demand for inelastic systems. When applied to elastic systems, the MPA procedure is equivalent to standard response spectrum analysis (RSA). The MPA estimates of seismic demand for torsionally‐stiff and torsionally‐flexible unsymmetric systems are shown to be similarly accurate as they are for the symmetric building; however, the results deteriorate for a torsionally‐similarly‐stiff unsymmetric‐plan system and the ground motion considered because (a) elastic modes are strongly coupled, and (b) roof displacement is underestimated by the CQC modal combination rule (which would also limit accuracy of RSA for linearly elastic systems). Copyright © 2004 John Wiley & Sons, Ltd.     

Approximate modal decomposition of inelastic dynamic responses of wall buildings

Two approximate methods for decomposing complicated inelastic dynamic responses of wall buildings into simple modal responses are presented. Both methods are based on the equivalent linear concept, where a non‐linear structure is represented by a set of equivalent linear models. One linear model is used for representing only one vibration mode of the non‐linear structure, and its equivalent linear parameters are identified from the inelastic response time histories by using a numerical optimizer. Several theoretical relations essential for the modal decomposition are derived under the framework of complex modal analysis. Various numerical examinations have been carried out to check the validity of the proposed modal decomposition methods, and the results are quite satisfactory in all cases. Fluctuating bending moment and shear at any location along the wall height contributed by each individual vibration mode can be obtained. Modal contributions to shear and flexural strength demands, as well as the corresponding modal properties, under various seismic loading conditions can also be identified and examined in detail. Furthermore, the effects of higher vibration modes on seismic demands of wall buildings are investigated by using the modal decomposition methods. Several new insights into the complicated inelastic dynamics of multi‐story wall buildings are presented. Copyright © 2004 John Wiley & Sons, Ltd.