System identification of a building from multiple seismic records

This paper describes the identification of finite dimensional, linear, time‐invariant models of a 4‐story building in the state space representation using multiple data sets of earthquake response. The building, instrumented with 31 accelerometers, is located on the University of California, Irvine campus. Multiple data sets, recorded during the 2005 Yucaipa, 2005 San Clemente, 2008 Chino Hills and 2009 Inglewood earthquakes, are used for identification and validation. Considering the response of the building as the output and the ground motion as the input, the state space models that represent the underlying dynamics of the building in the discrete‐time domain corresponding to each data set are identified. The time‐domain Eigensystem Realization Algorithm with the Observer/Kalman filter identification procedure are adopted in this paper, and the modal parameters of the identified models are consistently determined by constructing stabilization diagrams. The four state space models identified demonstrate that the response of the building is amplitude dependent with the response frequency and damping, being dependent on the magnitude of ground excitation. The practical application of this finding is that the consistency of this building response to future earthquakes can be quickly assessed, within the range of ground excitations considered (0.005g–0.074g), for consistency with prior response—this assessment of consistent response is discussed and demonstrated with reference to the four earthquake events considered in this study. Inclusion of data sets relating to future earthquakes will enable the findings to be extended to a wider range of ground excitation magnitudes. Copyright © 2010 John Wiley & Sons, Ltd.     

Resolution studies for Magnetic Resonance Sounding (MRS) using the singular value decomposition

We present a new approach to analyse the subsurface water content distribution obtained by inversion of MRS data in terms of resolution and penetration depth. It is based on a singular value decomposition (SVD) of the MRS forward operator to derive the model resolution matrix including regularisation parameters, i.e. including noise conditions. The approach takes loop size, subsurface resistivity distribution and noise conditions as input parameters affecting MRS into account and allows an assessment on resolution and penetration.The application of the new approach shows that the loop diameter must be carefully chosen depending on the investigation site to obtain optimal resolution and penetration. Using the introduced resolution measures the quality and reliability of the estimated model can be assessed.     

Use of USLE/GIS technology integrated with geostatistics to assess soil erosion risk in different land uses of Indagi Mountain Pass—Çankırı, Turkey

The universal soil loss equation (USLE) is an erosion model to estimate average soil loss that would generally result from splash, sheet, and rill erosion from agricultural plots. Recently, use of USLE has been extended as a useful tool predicting soil losses and planning control practices by the effective integration of the GIS-based procedures to estimate the factor values on a grid cell basis. This study was performed for five different lands uses of Indağı Mountain Pass, Cankırı to predict the soil erosion risk by the USLE/GIS methodology for planning conservation measures in the site. Of the USLE factors, rainfall-runoff erosivity factor (USLE-R) and topographic factor (USLE-LS) were greatly involved in GIS. These were surfaced by correcting USLE-R site-specifically using DEM and climatic data and by evaluating USLE-LS by the flow accumulation tool using DEM and watershed delineation tool to consider the topographical and hydrological effects on the soil loss. The study assessed the soil erodibility factor (USLE-K) by randomly sampled field properties by geostatistical analysis. Crop management factor for different land-use/land cover type and land use (USLE-C) was assigned to the numerical values from crop and flora type, canopy and density of five different land uses, which are plantation, recreational land, cropland, forest and grassland, by means of reclassifying digital land use map available for the site. Support practice factor (USLE-P) was taken as a unit assuming no erosion control practices. USLE/GIS technology together with the geostatistics combined these major erosion factors to predict average soil loss per unit area per unit time. Resulting soil loss map revealed that spatial average soil loss in terms of the land uses were 1.99, 1.29, 1.21, 1.20, 0.89 t ha⁻¹ year⁻¹ for the cropland, grassland, recreation, plantation and forest, respectively. Since the rate of soil formation was expected to be so slow in Central Anatolia of Turkey and any soil loss of more than 1 ton ha⁻¹ year⁻¹ over 50–100 years was considered as irreversible for this region, soil erosion in the Indağı Mountain Pass, to the great extent, attained the irreversible state, and these findings should be very useful to take mitigation measures in the site.     

Assessment of the reliability of 2D inversion of apparent resistivity data

The reliability of inversion of apparent resistivity pseudosection data to determine accurately the true resistivity distribution over 2D structures has been investigated, using a common inversion scheme based on a smoothness‐constrained non‐linear least‐squares optimization, for the Wenner array. This involved calculation of synthetic apparent resistivity pseudosection data, which were then inverted and the model estimated from the inversion was compared with the original 2D model. The models examined include (i) horizontal layering, (ii) a vertical fault, (iii) a low‐resistivity fill within a high‐resistivity basement, and (iv) an upfaulted basement block beneath a conductive overburden. Over vertical structures, the resistivity models obtained from inversion are usually much sharper than the measured data. However, the inverted resistivities can be smaller than the lowest, or greater than the highest, true model resistivity. The substantial reduction generally recorded in the data misfit during the least‐squares inversion of 2D apparent resistivity data is not always accompanied by any noticeable reduction in the model misfit. Conversely, the model misfit may, for all practical purposes, remain invariant for successive iterations. It can also increase with the iteration number, especially where the resistivity contrast at the bedrock interface exceeds a factor of about 10; in such instances, the optimum model estimated from inversion is attained at a very low iteration number. The largest model misfit is encountered in the zone adjacent to a contact where there is a large change in the resistivity contrast. It is concluded that smooth inversion can provide only an approximate guide to the true geometry and true formation resistivity.     

Numerical analysis of seismoelectric wave propagation in spatially confined geological units

Using numerical modelling, we investigate the evolution of seismoelectric effects induced by seismic excitation in spatially confined lithological units. Typical geometries represent clay lenses embedded in an aquifer or petroleum deposits in a host rock. In fluid‐saturated rocks, seismic waves can generate electromagnetic fields due to electrokinetic coupling mechanisms associated with such processes in the vicinity of the fluid‐mineral interface. Two seismoelectric phenomena are investigated: (1) the co‐seismic field associated with the seismic displacement at each point in a subsurface and (2) the interface response generated at layer boundaries. Our modelling uses a simplified time‐domain formulation of the coupled problem and an efficient 2D finite‐element implementation. To gain insight into the morphogenetic field behaviour of the seismoelectric effects, several numerical simulations for various target geometries were treated. Accordingly, we varied both the thickness of the confined units and the value of the electrical bulk conductivity in porous media. Analysis of these effects shows differences between interface responses for electrically conductive versus resistive units. So the pertinent contrast in electrical bulk conductivity controls the shape and structure of these seismoelectric conversion patterns. Moreover, the seismoelectric interface response captures both the petrophysical and geometrical characteristics of the geological unit. These models demonstrate the value of using seismoelectric interface response for reservoir characterization in either hydrogeological or hydrocarbon exploration studies.     

ESR and Th/U dating of speleothems from Aladağlar Mountain Range (AMR) in Turkey

Electron spin resonance (ESR) and ²³⁰Th/²³⁴U ages of speleothem samples collected from karstic caves located around 3000 m elevation in the Alada?lar Mountain Range (AMR), south-central Turkey, were determined in order to provide new insight and information regarding late Pleistocene climate. ESR ages were validated with the ²³⁰Th/²³⁴U ages of test samples. The ESR ages of 21 different layers of six speleothem samples were found to range mostly between about 59 and 4 ka, which cover the Marine Oxygen Isotope Stages (MIS) MIS 3 to MIS 1. Among all, only six layers appear to have deposited during MIS 8 and 5. Most of the samples dated were deposited during the late glacial stage (MIS 2). It appears that a cooler climate with a perennial and steady recharge was more conducive to speleothem development rather than a warmer climate with seasonal recharge in the AMR during the late Quaternary. This argument supports previous findings that suggest a two -fold increase in last glacial maximum mean precipitation in Turkey with respect to the present value.     

Shadow of Kerr-Taub-NUT black hole

The shadow of a rotating black hole with nonvanishing gravitomagnetic charge has been studied. It was shown that in addition to the angular momentum of black hole the gravitomagnetic charge term deforms the shape of the black hole shadow. From the numerical results we have obtained that for a given value of the rotation parameter, the presence of a gravitomagnetic charge enlarges the shadow and reduces its deformation with respect to the spacetime without gravitomagnetic charge. Finally we have studied the capture cross section for massive particles by black hole with the nonvanishing gravitomagnetic charge.     

Improvements in Inversion of Magnetic Resonance Exploration-Water Content, Decay Time, and Resistivity

In this review article, we present recent developments and Improvements in magnetic resonance sounding (MRS), a newly established geophysical exploration method that provides unique information about hydrogeophysical properties due to its direct sensitivity to hydrogen protons and proton dynamics. Starting with the most sophisticated and complete MRS formulation, we give a detailed view on how to solve the equation, i.e., inverting exactly for all model parameters: water content, decay time, and resistivity. Giving a short review of general inversion schemes used in geophysics, the special properties of MRS inversion are evaluated and the development of MRS inversion over recent years is shown. We present the extension of MRS to magnetic resonance tomography (MRT), i.e., the extension to two-dimensional investigations and appropriate inversions. Finally, we address restrictions, limitations, and inconsistencies as well as future developments.