RIDGE REGRESSION INVERSION APPLIED TO CRUSTAL RESISTIVITY SOUNDING DATA FROM SOUTH AFRICA

Geophysics ◽  
1977 ◽  
Vol 42 (5) ◽  
pp. 995-1005 ◽  
Author(s):  
W. R. Petrick ◽  
W. H. Pelton ◽  
S. H. Ward

Ridge regression inversion has been used to test the applicability of various one‐dimensional crustal models to the interpretation of deep Schlumberger sounding data from southern Africa (Van Zijl and Joubert, 1975). Four main models were investigated: a simple three‐layered earth, a layered earth with a transition zone exhibiting a linear decrease in log resistivity with depth, a similar earth with the transition zone determined by cubic splines, and a model having exponential resistivity behavior at depth. The last model corresponds to temperature‐dependent semiconduction through solid mineral grains (Brace, 1971). It was found that all of these models are capable of fitting the sounding data from southwestern Africa, while all except the semiconduction model fit the data from southeastern Africa. One is, thereby, immediately alerted to the problem of lack of resolution in Schlumberger sounding data where geologic control is not available. A major with the inversion of Schlumberger data alone is that accurate information is obtainable only for the resistivity‐thickness product of the resistive portion of the crust. On the other hand, magnetotelluric data, when available, tends to provide information on the thickness, but very little information on the true resistivity of the section. In order to resolve both resistivity and thickness it is possible to invert simultaneously Schlumberger and magnetotelluric (MT) data. Results obtained from the combined inversion of the African resistivity data and hypothetical MT data show that a considerable improvement in model resolution can be achieved using MT amplitude data even of poor accuracy from a relatively limited frequency range (0.1 to 100 Hz), whereas inclusion of MT phase information is of negligible additional benefit. Unfortunately, no significant test can be made, from data available at the time of our analysis, of the applicability of one‐dimensional inversion in a geologic circumstance which probably demands more dimension.

1987 ◽  
Vol 35 (3) ◽  
pp. 295-313
Author(s):  
K. Harmsen

A one-dimensional model for the movement of zinc sulphate through calcium-saturated soil is presented. Processes considered include mass flow, ion exchange, precipitation and dissolution. Precipitation occurs when the solubility product of gypsum is exceeded. In the presence of gypsum, ion exchange takes place at two separate interfaces, which move with different velocities through the soil. At the first interface precipitation of gypsum takes place in conjunction with ion exchange, and at the second interface the gypsum dissolves again and ion exchange proceeds until equilibrium is reached with the leaching solution. The composition of the transition zone between the two interfaces is calculated from the conditions of mass conservation and electroneutrality, the solubility product of gypsum and assuming a linear ion exchange equation. It is shown that the concentration of sulphate in the transition zone is higher than in the leaching solution, due to dissolution of gypsum at the second interface. In the presence of gypsum, zinc penetrates deeper into the soil than in its absence, but the fraction of the adsorption complex saturated with zinc is smaller. (Abstract retrieved from CAB Abstracts by CABI’s permission)


Geophysics ◽  
1976 ◽  
Vol 41 (4) ◽  
pp. 766-770 ◽  
Author(s):  
F. E. M. Lilley

Observed magnetotelluric data are often transformed to the frequency domain and expressed as the relationship [Formula: see text]where [Formula: see text] [Formula: see text] and [Formula: see text] [Formula: see text] represent electric and magnetic components measured along two orthogonal axes (in this paper, for simplicity, to be north and east, respectively). The elements [Formula: see text] comprise the magnetotelluric impedance tensor, and they are generally complex due to phase differences between the electric and magnetic fields. All quantities in equation (1) are frequency dependent. For the special case of “two‐dimensional” geology (where structure can be described as having a certain strike direction along which it does not vary), [Formula: see text] with [Formula: see text]. For the special case of “one‐dimensional” geology (where structure varies with depth only, as if horizontally layered), [Formula: see text] and [Formula: see text].


The one-dimensional inverse electromagnetic scattering problem for the inversion of amplitude data of either linear polarization state is investigated. The method exploits the complex structure of the field scattered from a class of inhomogeneous dielectrics and enables the analytic signal to be reconstructed from measurements of the amplitude alone. The method is demonstrated and exemplified with experimental data in both transverse electric and transverse magnetic polarization states. The implications of the method as a means for regularization of scattered data are briefly discussed.


1983 ◽  
Vol 88 (B3) ◽  
pp. 2407 ◽  
Author(s):  
Shimon Coen ◽  
Franchesca Quercia ◽  
Maria Mackiewicz

Geophysics ◽  
1985 ◽  
Vol 50 (10) ◽  
pp. 1618-1627 ◽  
Author(s):  
A. P. Raiche ◽  
D. L. B. Jupp ◽  
H. Rutter ◽  
K. Vozoff

One‐dimensional earth models consisting of uniform horizontal layers are useful both as actual representations of earth structures and as host models for more complex structures. However, there are often inherent difficulties in establishing layer thicknesses and resistivities from one type of measurement alone. For example, the dc resistivity method is sensitive to both conductive and resistive layers, but as these layers become thin, nonuniqueness becomes a severe problem. Electromagnetic (EM) methods are good for establishing the parameters of conductive layers, but they are quite insensitive to resistive layers. The use of both coincident loop transient EM (TEM) and Schlumberger methods, together with a joint inverse computer program, can vastly improve interpretation of layered‐earth parameters. The final model is less dependent upon starting guesses, error bounds are much improved, and nonuniqueness is much less of a problem. These advantages are illustrated by interpretation of real field data as well as by a theoretical study of four different types of earth models.


Sign in / Sign up

Export Citation Format

Share Document