A Constrained Nonlinear Inversion Approach to Quantitative Interpretation of Self-potential Anomalies Caused by Cylinders,Spheres and Sheet-like Structures

An interpretative method based on a nonlinearly mathematical optimization concept has been developed in this paper, in order to interpret self-potential anomalies (SP) due to horizontal cylinder, vertical cylinder, sphere and sheet-like structures. This interpretative method comprises three main steps. The first step is to formulate mathematically a nonlinearly constrained minimization problem (NCMP) to describe the geophysical problem related to the studied structure. The second one is to suggest an interior penalty function in order to convert the nonlinearly constrained minimization problem (NCMP) into a nonlinearly unconstrained minimization one (NUMP). The third step is to solve the converted nonlinearly unconstrained minimization problem (NUMP) by the well-known Hooke and Jeeves direct search algorithm in order to estimate the geophysical parameters of the studied structure, i.e., depth, polarization angle, electric dipole moment (magnitude of polarization) and geometric shape factor. The Hooke and Jeeves direct search algorithm is purposely chosen for being robust and its application to SP data allows a rapid convergence towards the optimal estimate of parameters. This interpretative method was first tested on theoretical synthetic models with different random noise, where a very close agreement was obtained between assumed and evaluated parameters.The validity of the proposed interpretative method is also tested on practical field examples taken from Turkey, India and Germany, where available SP data existed and was previously analyzed by different interpretative methods. The agreement between the results obtained by the developed method and those obtained by other published methods is good.Acknowledgment Authors would like to thank Dr. I. Othman Director General of the Atomic Energy Commission of Syria for his interest and continuous encouragement to achieve this work. Special thanks to the reviewers for their constructive suggestions aimed at enhancing the quality of this paper.     

A Robust Nonlinear Inversion for the Interpretation of Magnetic Anomalies Caused by Faults,Thin Dikes and Spheres Like Structure Using Stochastic Algorithms

A geophysical interpretative method is proposed to depth, amplitude coefficient (effective magnetization intensity), and index parameter (effective magnetization inclination) determination of a buried structure from magnetic field data anomaly due to a fault, a thin dike or a sphere-like structure. The method is based on the nonlinearly constrained mathematical modelling and also on the stochastic optimization approaches. The proposed interpretative method was first tested on a theoretical synthetic model with different random errors, where a very close agreement was obtained between the assumed and the evaluated parameters. The validity of this method was also tested on practical field data taken from United States, Australia, India and Brazil, where available magnetic data existed and were previously analyzed by different interpretative methods. The agreement between the results obtained by our developed method and those obtained by the other geophysical methods is good.     

A Versatile Nonlinear Inversion to Interpret Gravity Anomaly Caused by a Simple Geometrical Structure

A geophysical interpretative method is proposed to depth, amplitude coefficient and geometrical shape factor determination of a buried structure from an observed gravity anomaly related to a cylinder or a sphere-like structure.The method is based on nonlinearly constrained mathematical modelling and also on stochastic optimization approaches. The proposed interpretative method first has been tested on theoretical synthetic models with different random errors at a certain depth, where a very close agreement has been observed between assumed and evaluated parameters. Subsequent field data have been considered for which the interpreted results by other methods are available for comparison. The agreement between the obtained results by the proposed technique and by other geophysical methods is good. A statistical analysis has been also carried out to demonstrate the accuracy and the precision of the suggested interpretative method.     

Source parameters of some large earthquakes in the eastern Mediterranean region based on an itertive maximum entropy technique

TeleseismicP-waves of some large earthquakes that occurred in the eastern Mediterranean region have been analysed by using an iterative maximum entropy technique in order to obtain the independent spectral parameters, the long-period spectral level ₀ and the corner frequencyf ₀ of the far-field displacement spectra.Based on these parameters, the seismic source parameters seismic momentM ₀, source dimensionr, fault lengthl, average displacement u, shear stress drop , radiated energyE _s and apparent stressn are calculated for the considered earthquakes by using Brune's and Madariaga's models.The striking feature of the source parameters obtained in this study is the low stress drop value which varies between 5 and 15 bars. If Madariaga's model had been used, higher stress drop values would have been obtained.The low stress drop earthquakes in the eastern Mediterranean region might be interpreted either by the possible presence of low strength material near the source or by the partial stress drop model.     

Hydrochemical and microbiological characteristics of the Gusino-Ubukunskaya group of water bodies

Results of hydrochemical and microbiological studies of five largest water bodies belonging to the Gusino-Ubukunskaya group are presented. Present-day environmental state of the ecosystem subject to anthropogenic load is described. The possibility of using water quality parameters of the investigated lakes as indicators of the degree of contamination of waters flowing into the Baikal Lake is shown.Translated from Vodnye Resursy, Vol. 32, No. 1, 2005, pp. 79–84.Original Russian Text Copyright © 2005 by Khakhinov, Namsaraev, Ulzetueva, Barkhutova, Abidueva, Banzaraktsaeva.     

Least squares inversion of one-dimensional magnetotelluric data: An assessment of procedures employed by Brown University

The pirmary objective of this paper is to review inversion algorithms employed to interpret magnetotelluric data by the geophysical/electromagnetics group at Brown University. Examples are presented from a variety of tectonic regimes — the Rio Grande rift, the Colorado Plateau, and Iceland. Following a brief tutorial in which the basic relations for the inverse prrolem are developed, we discuss the detectability and resolvability of various elements of the models. The correlation between certain model parameters is stressed, a phenomena well-known to interpreters. In this paper we emphasize the importance of identifying the precise range of values for which two or more parameters are coupled — these ranges are defined by what wer term saturation limits. Generalized inverse theory is employed to arrive at optimal models which are evaluated in terms of resolution and information matrices. The degree to which models fit a given dat base can be assessed using several quantitative approaches. We have found that it is impotant not only to determine the best fitting model, but also to identify those marginal solutions which also fall within the range of uncertainty of the data. The results from our direct invverse method agree favorably with those from Monte Carlo methods.     

Electrical conductivity of the earth's crust and upper mantle

This review summarizes recent results of electrical resistivity studies of the earth's crust and upper mantle. Where available, the data are discussed in the context of further regional geophysical information. Electrical resistivity is very sensitive to a wide range of petrological and physical parameters, e.g., to carbon, fluids, volatiles and enhanced temperatures, making electrical resistivity methods a powerful tool in crust and upper mantle investigations. Yet, the general increase in resistivity data of the crust and mantle has not ended the battle of explanations for anomalous crustal conductivities.     

Interpretation of magnetic anomalies using a genetic algorithm

A genetic algorithm (GA) is an artificial intelligence method used for optimization. We applied a GA to the inversion of magnetic anomalies over a thick dike. Inversion of nonlinear geophysical problems using a GA has advantages because it does not require model gradients or well-defined initial model parameters. The evolution process consists of selection, crossover, and mutation genetic operators that look for the best fit to the observed data and a solution consisting of plausible compact sources. The efficiency of a GA on both synthetic and real magnetic anomalies of dikes by estimating model parameters, such as depth to the top of the dike (H), the half-width of the dike (B), the distance from the origin to the reference point (D), the dip of the thick dike (δ), and the susceptibility contrast (k), has been shown. For the synthetic anomaly case, it has been considered for both noise-free and noisy magnetic data. In the real case, the vertical magnetic anomaly from the Pima copper mine in Arizona, USA, and the vertical magnetic anomaly in the Bayburt–Sar?han skarn zone in northeastern Turkey have been inverted and interpreted. We compared the estimated parameters with the results of conventional inversion methods used in previous studies. We can conclude that the GA method used in this study is a useful tool for evaluating magnetic anomalies for dike models.     

Does the interplanetary magnetic field and the solar activity contribute to the tropospheric circulation?

A significant tendency is shown for both Etesians and (+, –) sector boundaries of the interplanetary magnetic field (IMF) to occur on the same solar rotation days, during the main period of the Etesians effect (July–August). In addition, the solar activity seems to control the Etesians distribution within the IMF sector structure. In the epoch of maximum there is a significant tendency of Etesians to occur during toward IMF days. In contrast, in the epoch of minimum Etesians occur mainly during away IMF days. Finally, in the epoch of intermediate the Etesians are uniformly distributed in away and toward IMF days. Since these conclusions are statistically significant at high confidence levels, it is fair to assume that IMF and solar activity seem to contribute, to some extent, to the Etesians occurrences, as well as to their distribution within the solar rotation and the IMF sector structure; that is, some solar contribution to the tropospheric circulation is implied.     

Perturbation of eigenfrequencies of a cylinder due to a ridge on the surface

Summary An infinite cylinder is assumed to have a small ridge on the surface. The effect of the ridge on the eigenfrequencies of flute mode of oscillation is obtained in the linear approximation. The usefulness of these results in seismological problems is indicated.     

Thermomagnetic properties of some Late Mesozoic diabase dikes of South Spitsbergen

Summary Laboratory investigations of thermomagnetic properties of specimens from seven diabase dikes of South Spitsbergen have been performed. Studies of temperature variation of isothermal saturation remanence combined with Hopkinson effect measurements and investigations of low temperature remanence give information on the mineral composition of the magnetic phases present in the magnetic fraction of the diabases. The experiments show that the clean natural remanence of these dikes has been acquired at the time of their, consolidation in the Late Mesozoic, though three of them may have been partially remagnetized at the time of Tertiary orogeny. Since the dikes reveal several normal and reversed directions of clean natural remanence, it is concluded that they were formed at different times in the Late Mesozoic, even though some of them may have been subjected to secondary influences at a later date.     

A brief outline of the magnetic susceptibility anisotropy of various rock types

Summary In the first part of this paper, the main geological and geophysical applications of magnetic susceptibility anisotropy are briefly described. In the second part, the data on the magnetic anisotropy of various rock types are summarized and some conclusions concerning its significance are deduced. For this purpose, all accessible data were used.Dedicated to 90th Birthday of Professor Frantiek Fiala     

Numerical modeling of the high-latitude F-layer anomalies

The results of a detailed numerical study of the behaviour of the convective polar ionosphere are presented. The developed theoretical model produces three-dimensional distribution of electron density, electron and ion temperature. The effects of auroral particle precipitation on the density and temperature structures are studied for winter and low solar activity conditions. The high-latitude ionospheric features, such as the tongue of ionization, the main trough the polar ionization peak, the auroral ionization peak, the high-latitude ionization hole, the tongue of electron temperature, the high latitude minima of electron, temperature, and the ion temperature hot spot are obtained from calculations.Numerically obtained results are used for determining the HF propagation paths in the polar ionosphere. The effects of ionospheric irregularities on high frequency ionospheric radio waves are investigated by using a three-dimensional ray-tracing computer program. Ray-path trajectories are presented for different values of the elevation angle of transmission. From our study, it was found that large-scale irregularity structures of the high-latitude ionosphere, in the presence of the earth's magnetic field, significantly affect high frequency radio wave propagation.     

Statistical methods for river runoff prediction

Methods used to analyze one type of nonstationary stochastic processes—the periodically correlated process—are considered. Two methods of one-step-forward prediction of periodically correlated time series are examined. One-step-forward predictions made in accordance with an autoregression model and a model of an artificial neural network with one latent neuron layer and with an adaptation mechanism of network parameters in a moving time window were compared in terms of efficiency. The comparison showed that, in the case of prediction for one time step for time series of mean monthly water discharge, the simpler autoregression model is more efficient.Translated from Vodnye Resursy, Vol. 32, No. 2, 2005, pp. 133–145.Original Russian Text Copyright © 2005 by Pisarenko, Lyubushin, Bolgov, Rukavishnikova, Kanyu, Kanevskii, Saveleva, Demyanov, Zalyapin.     

Complex image method for calculating electric and magnetic fields produced by an auroral electrojet of finite length

The electromagnetic field due to ionospheric currents has to be known when evaluating space weather effects at the earths surface. Forecasting methods of these effects, which include geomagnetically induced currents in technological systems, are being developed. Such applications are time-critical, so the calculation techniques of the electromagnetic field have to be fast but still accurate. The contribution of secondary sources induced within the earth leads to complicated integral formulas for the field at the earths surface with a time-consuming computation. An approximate method of calculation based on replacing the earth contribution by an image source having mathematically a complex location results in closed-form expressions and in a much faster computation. In this paper we extend the complex image method (CIM) to the case of a more realistic electrojet system consisting of a horizontal line current filament with vertical currents at its ends above a layered earth. To be able to utilize previous CIM results, we prove that the current system can be replaced by a purely horizontal current distribution which is equivalent regarding the total (= primary + induced) magnetic field and the total horizontal electric field at the earths surface. The latter result is new. Numerical calculations demonstrate that CIM is very accurate and several magnitudes faster than the exact conventional approach.     

Seismic velocity-ratio in the crust and uppermost mantle in southeast Wellington province,New Zealand

The regional variation of the seismic velocity-ratio () over a 200 km long traverse has been studied by means of microearthquake surveys. The Wadati-plot method is used with a minimum of four P and S arrivals for each of 49 earthquakes. The area as a whole is found to be characterized by a value of 1.74–1.76 for earthquakes of depth 12–40 km, except in a 50 km long section near Wellington, where is low at 1.60. This low has been attributed to the fault zones in the region. A small change of is observed between the upper crust (5 km) and lower crust (12 km), but there is no change of between the lower crust and uppermost mantle.     

Space Plasma Exploration by Active Radar (SPEAR): an overview of a future radar facility

SPEAR is a new polar cap HF radar facility which is to be deployed on Svalbard. The principal capabilities of SPEAR will include the generation of artificial plasma irregularities, operation as an all-sky HF radar, the excitation of ULF waves, and remote sounding of the magnetosphere. Operation of SPEAR in conjunction with the multitude of other instruments on Svalbard, including the EISCAT Svalbard radar, and the overlap of its extensive field-of-view with that of several of the HF radars in the SuperDARN network, will enable in-depth diagnosis of many geophysical and plasma phenomena associated with the cusp region and the substorm expansion phase. Moreover, its ability to produce artificial radar aurora will provide a means for the other instruments to undertake polar cap plasma physics experiments in a controlled manner. Another potential use of the facility is in field-line tagging experiments, for coordinated ground-satellite experiments. Here the scientific objectives of SPEAR are detailed, along with the proposed technical specifications of the system.     

A Local Least-Squares Modification of Stokes’ Formula

The combination of Stokes formula and an Earth Gravity Model (EGM) for geoid determination has become a standard procedure. However, the way of modifying Stokes formula vary from author to author, and numerous methods of modification exist. Most methods are deterministic, with the primary goal of reducing the truncation bias committed by limiting the area of Stokes integration around the computation point, but there are also some stochastic methods with the explicit goal to reduce the global mean square error of the geoid height estimator stemming from the truncation bias as well as the random errors of the EGM and the gravity data. The latter estimators are thus, at least from a theoretical point of view, optimal in a global mean sense, but in a local sense they may be far from optimality.Here we take advantage of the error variance-covariance matrices of the EGM and the terrestrial gravity data to derive the modification parameters of Stokes kernel in a local least-squares sense. The solution is given for the unbiased type of modification of Stokes formula of Sjöberg (1991).     

Paleomagnetic data analysis and the continents' drift

Summary The results from the representation of the Earth's magnetic field of the different geological epochs (since the Silurian epoch till contemporary times) with the field of an average magnetic dipole show that the eccentric dipole is displaced in the direction of the Pacific Ocean. The difference between the poles' migration curves of the different continents can be explained by the eccentricity of this dipole.     

An electromagnetic prospecting instrument for subsurface conducting zones

Summary A portable electromagnetic prospecting unit comprising a transmitter, receiver and resolver, designed and constructed to operate with a vertical setup of coils at a single frequency of 1250 cps, measures the components of the magnetic field at the receiver station resolved in phase and at quadrature with field near the transmitter. This may be used to locate vertical or steeply dipping conductors down to a depth of about 0.6 times the maximum workable transmitter receiver separation of 120 meters, the accuracy of measurements of the resolved components being better than 2%. The variation of resolved components over a buried conductor indicates its position and depth.This paper is published by the kind permission of the Director General, Geological Survey of India.