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.     

Metamorphism and metamorphic K–Ar ages of the Mesozoic accretionary complex in Northland, New Zealand

Abstract A southwest dipping Mesozoic accretionary complex, which consists of tectonically imbricated turbiditic mudstone and sandstone, hemipelagic siliceous mudstone, and bedded cherts and basaltic rocks of pelagic origin, is exposed in northern North Island, New Zealand. Interpillow limestone is sometimes contained in the basaltic rocks. The grade of subduction‐related metamorphism increases from northeast to southwest, indicating an inverted metamorphic gradient dip. Three metamorphic facies are recognized largely on the basis of mineral parageneses in sedimentary and basaltic rocks: zeolite, prehnite‐pumpellyite and pumpellyite‐actinolite. From the apparent interplanar spacing d₀₀₂ data for carbonaceous material, which range from 3.642 to 3.564 Å, the highest grade of metamorphism is considered to have attained only the lowermost grade of the pumpellyite‐actinolite facies for which the highest temperature may be approximately 300°C. Metamorphic white mica K–Ar ages are reported for magnetic separates and <2 µm hydraulic elutriation separates from 27 pelitic and semipelitic samples. The age data obtained from elutriation separates are approximately 8 m.y. younger, on average, than those from magnetic separates. The age difference is attributed to the possible admixture of nonequilibrated detrital white mica in the magnetic separates, and the age of the elutriation separates is considered to be the age of metamorphism. If the concept, based on fossil evidence, of the subdivision of the Northland accretionary complex into north and south units is accepted, then the peak age of metamorphism in the north unit is likely to be 180–130 Ma; that is, earliest Middle Jurassic to early Early Cretaceous, whereas that in the south unit is 150–130 Ma; that is, late Late Jurassic to early Early Cretaceous. The age cluster for the north unit correlates with that of the Chrystalls Beach–Taieri Mouth section (uncertain terrane), while the age cluster for the south unit is older than that of the Younger Torlesse Subterrane in the Wellington area, and may be comparable with that of the Nelson and Marlborough areas (Caples and Waipapa terranes).     

Heat influx and exhumation of the Shimanto accretionary complex: Miocene fission track ages from the Kii Peninsula, southwest Japan

Abstract   To determine how local geological events contributed to the evolution of accretionary complexes and eventual exposure of rocks with different structural levels, geochronological mapping was carried out using fission track (FT) analysis at the Kii Peninsula, southwest Japan. At this site, the original zonal structure of Cretaceous accretionary complexes parallel to the subduction zone is disturbed by the northward projection of the Shimanto accretionary complex. Twenty-six zircon FT ages were obtained from an area of ∼12 km in an east–west direction and ∼15 km in a north–south direction, and classified into three groups: (i) ages ∼15 Ma (range ∼10–20 Ma), which are distributed along the northwest–southeast valley; (ii) ages of ∼50 Ma in the northwest of the study area; and (iii) ages older than those in Groups 1 and 2. Based on results from eight zircon FT length distributions, the Miocene ages appear to be the result of spatial variations in heat influx and cooling after the regional exhumation of the area, as recorded by FT ages of ∼50 Ma.     

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.     

Integrated spatial sampling modeling of geospatial data

Sampling is to, by efficient selection of samples, acquire the accurate information about the population (the research object) at less cost. Spatial sampling is a kind of sampling toward geospatial objects or features with spatial correlation. The differences between effi-cient sampling and completely universal survey lie in quality, time and cost. Sampling provides a kind of economical, prompt and accurate survey[13]. Efficient spatial sampling can be regarded as the optimization of the sampl…     

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.