Crustal structure using receiver function in the east part of A＇nyemaqen suture belt

Twenty broadband seismographs were deployed along Hongyuan, Sichuan to Wuwei, Gansu. 81 teleseismic events were recorded in one year. We computed receiver functions from teleseismic waveform data and obtained S wave velocity structure beneath each station along the profile by using receiver function inversion method. The results revealed that the crustal structure is very complex and crustal average S wave velocity is to be on the low side. Low velocity structure generally exists in the depth range of 10～40 km in the crust between Aba arc fault and northern edge fault of Qinling earth＇s axis and it is a tectonic feature of complex geological process such as ancient A＇nyemaqen Tethys ocean from closing and side colliding to subducted plate exhumed or thrust rock slice lifted. The Moho is about 50 km depth along the profile and is slightly deeper in the south than in the north.     

Three-dimensional thermal structure of the Chinese continental crust and upper mantle

We invert S-wave velocities for the 3D upper-mantle temperatures, in which the position with a temperature crossing the 1300℃ adiabat is corresponding to the top of the seismic low velocity zone. The temperatures down to the depth of 80 km are then calculated by solving steady-state thermal conduction equation with the constraints of the inverted upper-mantle temperatures and the surface temperatures, and then surface heat flows are calculated from the crustal temperatures. The misfit between the calculated and observed surface heat flow is smaller than 20% for most regions. The result shows that, at a depth of 25 km, the crustal temperature of eastern China (500―600℃) is higher than that of western China (<500℃). At a depth of 100 km, temperatures beneath eastern and southeastern China are higher than the adiabatic temperature of 1300℃, while that beneath west China is lower. The Tarim craton and the Sichuan basin show generally low temperature. At a depth of 150 km, temperatures beneath south China, eastern Yangtze craton, North China craton and around the Qiangtang terrane are higher than the adiabatic temperature of 1300℃, but is the lowest beneath the Sichuan basin and the regions near the Indian-Eurasian collision zone. At a depth of 200 km, very low temperature occurs beneath the Qinghai-Tibet Plateau and the south to the Tarim craton.     

Automatic Interpretation of Magnetic Data Using Euler Deconvolution with Nonlinear Background

The voluminous gravity and magnetic data sets demand automatic interpretation techniques like Naudy, Euler and Werner deconvolution. Of these techniques, the Euler deconvolution has become a popular choice because the method assumes no particular geological model. However, the conventional approach to solving Euler equation requires tentative values of the structural index preventing it from being fully automatic and assumes a constant background that can be easily violated if the singular points are close to each other. We propose a possible solution to these problems by simultaneously estimating the source location, depth and structural index assuming nonlinear background. The Euler equation is solved in a nonlinear fashion using the optimization technique like conjugate gradient. This technique is applied to a published synthetic data set where the magnetic anomalies were modeled for a complex assemblage of simple magnetic bodies. The results for close by singular points are superior to those obtained by assuming linear background. We also applied the technique to a magnetic data set collected along the western continental margin of India. The results are in agreement with the regional magnetic interpretation and the bathymetric expressions.     

Effect of surface stress on magneto-elastic surface waves in finitely conducting media

This paper deals with the investigation of the effect of surface stress and conductivity on the propagation of surface wave in isotropic, homogeneous, elastic media under the action of a primary magnetic field. Formulation of the general surface wave propagation problem has been made, and the corresponding frequency equation has been derived. Frequency equations for Rayleigh wave, surface shear wave and Stoneley wave have been deduced from that of general surface wave as special cases. The effects of surface stress and magnetic field on the wave velocities and attenuation factors of Rayleigh wave and surface shear wave are shown by numerical calculations and graphs. Some important wave velocity equations, as obtained by other authors, have been deduced as special cases from the wave velocity equation for Stoneley wave. It is found that the combined effect of surface stress and magnetic field modulates the wave velocity ratios and attenuation factors of Rayleigh wave and surface shear wave to a considerable extent.     

Use of variograms for field magnetometry analysis in Upper Silesia Industrial Region

Measurements of topsoil magnetic susceptibility are often used for quick assessment of soil contamination of anthropogenic origin, with heavy metals or other pollutants. However, because of complicated correlations between low-field magnetic susceptibility (shortened to magnetic susceptibility) of topsoil and soil pollution, the outcome of a field magnetometry survey can not be related directly to soil pollution. For each case study, the results should be interpreted on their own taking into account not only the type of pollution but also pedogenic, biogenic and environmental factors. In practice, it is very difficult to measure and consider all these factors. Here we illustrate the merit of geostatistical methods, which are focused on the spatial variability of a phenomenon, in the interpretation of soil magnetometry results. This article presents the analysis of spatial variability of top soil layers magnetic susceptibility-within the Upper Silesia Industrial Region (USIR)-using semivariance analysis. It also explains how to adjust the sampling density of field magnetometry measurements to spatial variability of the soil pollution as well as to the spatial scale of the investigated area. For this purpose, the values of magnetic susceptibility have been measured by using various sampling densities at areas of different size located within USIR. This enabled to determine the main scales of magnetic susceptibility spatial variability of soils within USIR using semivariance. A few distinct scales of variability were found from the site scale to a more regional scale. Variability ranges of 30 km, 12 km, and 5 km refer to the large regional scale, whereas smaller ranges of few hundreds down to a few tens of meters, can be attributed to the local (site) scale. In addition, the precision of the measuring campaigns, performed within USIR with different sampling densities, was compared through the analysis of the spatial variability of the soil magnetic susceptibility signal by using ordinary kriging. jarek97@yahoo.com, piotr.fabijanczyk@is.pw.edu.pl     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.     

Modeling doline populations with logistic growth functions

Initiation and evolution of a large karst depression population consisting of parent and daughter dolines were spatially and temporally modeled using logistic growth functions. Logistic growth models are well suited for analyzing doline population initiation and evolution because they reflect the density‐dependent growth mechanisms present in the evolution of karst depressions. Seven assumptions based upon previous studies were refined into mathematical statements and tested using more than 2000 dolines from a subpopulation of both parent and daughter karst depressions on the Western Highland Rim and Pennyroyal Plain of Tennessee and Kentucky. Logistic growth models quantify the initiation and evolution of doline populations and interface well with recent models describing the evolution of three‐dimensional conduit systems. Logistic growth models should apply to modeling other doline populations in karst terranes with hydraulically efficient, three‐dimensional conduit systems. Copyright © 2006 John Wiley & Sons, Ltd.     

Earthquake‐resistant structural design through energy demand and capacity

An energy‐based earthquake‐resistant structural design method is proposed. The proposed method uses specific input energy spectra, modal or time‐history analyses, and energy distribution among structural members. For a given member strength and stiffness, a relationship between the energy attributable to damage absorbed by a member and its cumulative ductility demand can be determined. Member strength, stiffness and energy capacity are design parameters which are simultaneously used in the design. The method can avoid soft‐storey design. The damage is measured based on a cumulative basis considering earthquake magnitude, frequency, and duration. Tests have been carried out to determine energy absorbing capacities of various structural components. More efforts are needed to make the energy‐based earthquake‐resistant structural design practical, but ssimple formulations for this method are possible. Copyright © 2007 John Wiley & Sons, Ltd.     

Stochastic assessment of GPS carrier phase measurements for precise static relative positioning

 Global positioning system (GPS) carrier phase measurements are used in all precise static relative positioning applications. The GPS carrier phase measurements are generally processed using the least-squares method, for which both functional and stochastic models need to be carefully defined. Whilst the functional model for precise GPS positioning is well documented in the literature, realistic stochastic modelling for the GPS carrier phase measurements is still both a controversial topic and a difficult task to accomplish in practice. The common practice of assuming that the raw GPS measurements are statistically independent in space and time, and have the same accuracy, is certainly not realistic. Any mis-specification in the stochastic model will inevitably lead to unreliable positioning results. A stochastic assessment procedure has been developed to take into account the heteroscedastic, space- and time-correlated error structure of the GPS measurements. Test results indicate that the reliability of the estimated positioning results is improved by applying the developed stochastic assessment procedure. In addition, the quality of ambiguity resolution can be more realistically evaluated. Received: 13 February 2001 / Accepted: 3 September 2001     

GPS network monitors the Western Alps' deformation over a five-year period: 1993–1998

  The Western Alps are among the best studied collisional belts with both detailed structural mapping and also crustal geophysical investigations such as the ECORS and EGT seismic profile. By contrast, the present-day kinematics of the belt is still largely unknown due to small relative motions and the insufficient accuracy of the triangulation data. As a consequence, several tectonic problems still remain to be solved, such as the amount of N–S convergence in the Occidental Alps, the repartition of the deformation between the Alpine tectonic units, and the relation between deformation and rotation across the Alpine arc. In order to address these problems, the GPS ALPES group, made up of French, Swiss and Italian research organizations, has achieved the first large-scale GPS surveys of the Western Alps. More than 60 sites were surveyed in 1993 and 1998 with a minimum observation of 3 days at each site. GPS data processing has been done by three independent teams using different software. The different solutions have horizontal repeatabilities (N–E) of 4–7 mm in 1993 and 2–3 mm in 1998 and compare at the 3–5-mm level in position and 2-mm/yr level in velocity. A comparison of 1993 and 1998 coordinates shows that residual velocities of the GPS marks are generally smaller than 2 mm/yr, precluding a detailed tectonic interpretation of the differential motions. However, these data seem to suggest that the N–S compression of the Western Alps is quite mild (less than 2 mm/yr) compared to the global convergence between the African and Eurasian plate (6 mm/yr). This implies that the shortening must be accomodated elsewhere by the deformation of the Maghrebids and/or by rotations of Mediterranean microplates. Also, E–W velocity components analysis supports the idea that E–W extension exists, as already suggested by recent structural and seismotectonic data interpretation. Received: 27 November 2000 / Accepted: 17 September 2001