An assessment of inversion methods for AEM data applied to environmental studies

Airborne electromagnetic (AEM) surveys are currently being flown over populated areas and applied to detailed problems using high flight line densities. Interpretation information is supplied through a model of the subsurface resistivity distribution. Theoretical and survey data are used here to study the character and reliability of such models. Although the survey data were obtained using a fixed-wing system, the corresponding associations with helicopter, towed-bird systems are discussed. Both Fraser half-space and 1D inversion techniques are considered in relation to their ability to distinguish geological, cultural and environmental influences on the survey data. Fraser half-space modelling provides the dual interpretation parameters of apparent resistivity and apparent depth at each operational frequency. The apparent resistivity was found to be a remarkably stable parameter and appears robust to the presence of a variety of at-surface cultural features. Such features provide both incorrect altitude data and multidimensional influences. Their influences are observed most strongly in the joint estimate of apparent depth and this accounts for the stability of the apparent resistivity. Positive apparent depths, in the example data, result from underestimated altitude measurements. It is demonstrated that increasingly negative apparent depths are associated with increasing misfits between a 1D model and the data. Centroid depth calculations, which are a transform of the Fraser half-space parameters, provide an example of the detection of non-1D influences on data obtained above a populated area. 1D inversion of both theoretical and survey data is examined. The simplest use of the 1D inversion method is in providing an estimate of a half-space resistivity. This can be undertaken prior to multilayer inversion as an initial assessment. Underestimated altitude measurements also enter the problem and, in keeping with the Fraser pseudo-layer concept, an at-surface highly resistive layer of variable thickness can be usefully introduced as a constrained parameter. It is clearly difficult to ascribe levels of significance to a ‘measure’ of misfit contained in a negative apparent depth with the dimensions of metres. The reliability of 1D models is better assessed using a formal misfit parameter. With the misfit parameter in place, the example data suggest that the 1D inversion methods provide reliable apparent resistivity values with a higher resolution than the equivalent information from the Fraser half-space estimates.     

Theoretical study of slope effects in resistivity surveys and applications

When electric soundings are made over an irregular terrain, topographic effects can influence the values of apparent resistivity and lead to erroneous 1D interpretation. A 3D finite-element method has been applied to study the topographical effect of a slope on Schlumberger soundings parallel to the strike. When the resistivity survey is performed at the top of the slope, the apparent resistivity values can be two times higher than in the flat-earth case, depending on the angle (α) and height (H) of the slope, and on the distance (X) between the sounding and the slope top. The results are presented as nondimensional curves which can be used for evaluating topographic anomalies for any value of the parameters α, H and X. It is numerically shown that the topographic effects can be removed from measurements on horizontally layered structures with an irregular earth surface. Real measurements were performed in different geological conditions over an irregular terrain. The correction method based on the nondimensional curves has been applied to the data and has enabled the determination of the correct layered ground configuration using 1D interpretation.     

AIRBORNE RESISTIVITY MAPPING USING A MULTIFREQUENCY ELECTROMAGNETIC SYSTEM*

Airborne electromagnetic methods are most commonly used in mineral exploration. However, new developments, such as multifrequency capability and digital on-board field recording, as well as improvements in instrumentation resulting in high signal-to-noise ratios in recorded data, have made their application in geological mapping possible. A three-frequency airborne EM survey carried out over an area northwest of Timmins, Ontario, was interpreted in terms of thickness and resistivity of the layers of a two-layer earth section. Since both in-phase and quadrature components are measured, this provides six independent parameters at each point in space. Based on prior geological information and a preliminary interpretation of the field records, two two-layer models of the subsurface seemed to be appropriate for most of the survey area. An automatic computerized interpretation procedure was devised to interpret the field data at each point in terms of thickness and resistivity parameters of those two models. When the geology is more complex, the data do not fit the models and no interpretations are made. Two maps illustrating the variation of the resistivity and the thicknesses of the layers were constructed from the interpreted data. These maps agree with the known geological information about the distribution of glacial clay in the area. Areas where the layered models do not fit are known to be areas where the geology is complex with a large number of dykes and other lateral inhomogeneities. The study shows that multifrequency airborne EM surveys can be very useful in geological mapping over inaccessible terrain and can significantly help the mapping geologist where outcrops are scarce.     

Geoelectric structure of the Ein Gedi sinkhole occurrence site at the Dead Sea shore in Israel

Hazardous sinkholes started to appear in alluvial fans and unconsolidated sediments along the western Dead Sea coast in 1990. Since then hundreds of sinkholes have appeared from north to south along the shoreline. The Electrical Resistivity Tomography (ERT) method was used to achieve a better understanding of the subsurface geoelectric structure at the sinkhole development sites, taking into account that electric parameters (such as resistivity or conductivity) are very sensitive to formation properties and their variations in time. Fifteen image lines were surveyed at the Ein Gedi area during a period of active sinkhole development (in 2001–2002) over an area of 300 × 550 m². Resistivity cross-sections and maps were constructed from 2-D linear surveys. The process of sinkhole formation in the surveyed area is located in a strip 50–70 m wide and 300–500 m long, extending approximately in a north–south direction. The sinkholes are arranged along a tortuous line within this strip. On resistivity maps and sections this U-shaped zone appears as an alternation of high resistivity anomalies of 350–1000 Ωm (at sinkhole group locations) with narrow background resistivity zones of 50–100 Ωm. The large size of resistivity anomalies (250 × 300 m²), which are considerably greater than those of the sinkholes, form one of the features of the sinkhole sites in the Ein Gedi area. The anomalies continue down to the water table or even deeper (maximum of 25–35 m depth). A low resistivity layer of 1–8 Ωm underlies them. The combined analysis of the image results and other geophysical data shows that high resistivity anomalies are associated with the decompaction of the soil mass at the sinkhole development sites and surrounding areas. Recent studies have shown that sinkholes in the Ein Gedi area are developing along the salt western edge located at a depth of 50 m. The subsurface high resistivity anomaly conforms to the sinkhole line (and salt boundary). They are presumably located above the great dissolution caverns at the salt edge. The heterogeneity of the resistivity structure within the high resistivity anomaly (seen in both lateral and vertical planes) confirms that a disintegration of internal formation structure takes place. Away from the sinkhole sites the subsurface resistivity distribution is homogeneous.     

The dynamics of a rigid foundation on the surface of an elastic half-space

The response of a rigid rectangular foundation block resting on an elastic half-space has been determined by considering first the displacement functions for any position on the surface of an unloaded half-space due to a harmonic point force. The influence of the foundation has been taken into account by assuming a relaxed condition at the interface, i.e. a uniform displacement under the foundation and that the sum of the point forces must equal the total applied force. The three motions of vertical, horizontal and rocking have been considered and numerical values for the in-phase and the quadrature components of the displacement functions are presented for a Poisson's ratio of 0.25. The effect of the mass and inertia of the foundation can be allowed for by an impedance matching technique. Response curves and non-dimensional resonant frequency curves are given for a square and a rectangular foundation for different mass and inertia ratios and for several values of Poisson's ratio. These curves are for design purposes and are an addition to similar published curves for circular and infinitely long rectangular foundations. Some of the calculated results have been verified by a laboratory experiment.     

Induction arrows for a buried conducting plate

A scaled electromagnetic analogue model is used to study the behaviour of single station induction arrows over a conducting plate embedded in a homogeneous poorly conducting medium. For a profile crossing the edge of the conducting plate, the in-phase arrow, for all locations, points towards the bulk of the good conductor while the quadrature arrow reverses direction directly over the edge, consistently pointing towards the edge. The results also show that both the in-phase and quadrature arrow lengths are reduced with increasing conductor thickness and depth of burial. The model results are applied to examples of a uniform depth ocean and a conducting sill embedded in a poorly conducting Earth.     

Auswertung von Consolfunkfeuerpeilungen

Summary A table for the North Sea is given containing the Consol position lines radiating from Bushmills and Stavanger. Taking the coordinates of two points of the bearing from the tables according to the number of dots or dashes observed, one can immediately plot a part of the position line near the DR position.

---

---

Hierzu Tafel 4–7 mit Tabelle 1–4     

Reducing discrepancies in ground and satellite-observed eruption heights

The plume height represents a crucial piece of evidence about an eruption, feeding later assessment of its size, character, and potential impact, and feeding real-time warnings for aviation and ground-based populations. There have been many observed discrepancies between different observations of maximum plume height for the same eruption. A comparison of maximum daily height estimates of volcanic clouds over Indonesia and Papua New Guinea during 1982–2005 shows marked differences between ground and satellite estimates, and a general tendency towards lower height estimates from the ground. Without improvements in the quality of these estimates, reconciled among all available methods, warning systems will be less effective than they should be and the world's record of global volcanism will remain hard to quantify. Examination of particular cases suggests many possible reasons for the discrepancies. Consideration of the satellite and radar cloud observations for the 1991 Pinatubo eruptions shows that marked differences can exist even with apparently good observations. The problem can be understood largely as a sampling issue, as the most widely reported parameter, the maximum cloud height, is highly sensitive to the frequency of observation. Satellite and radar cloud heights also show a pronounced clumping near the height of the tropopause and relative lack of eruptions reaching only the mid-troposphere, reinforcing the importance of the tropopause in determining the eruption height in convectively unstable environments. To reduce the discrepancies between ground and satellite estimates, a number of formal collaboration measures between vulcanological, meteorological and aviation agencies are suggested.     

Pyroclastic current dynamic pressure from aerodynamics of tree or pole blow-down

The common occurrence of tree and pole blow-down from pyroclastic currents provides an opportunity to estimate properties of the currents. Blow-down may occur by uprooting (root zone rupture), or flexure or shear at some point on the object. If trees are delimbed before blow-down, each tree or pole can be simulated by a cylinder perpendicular to the current. The force acting on a cylinder is a function of flow dynamic pressure, cylinder geometry, and drag coefficient. Treated as a cantilever of circular cross-section, the strength for the appropriate failure mode (rupture, uprooting or flexure) can then be used to estimate the minimum necessary current dynamic pressure. In some cases, larger or stronger standing objects can provide upper bounds on the dynamic pressure. This analysis was treated in two ways: (1) assuming that the current properties are vertically constant; and (2) allowing current velocity and density to vary vertically according to established models for turbulent boundary layers and stratified flow. The two methods produced similar results for dynamic pressure. The second, along with a method to approximate average whole-current density, offers a means to estimate average velocity and density over the height of the failed objects. The method is applied to several example cases, including Unzen, Mount St. Helens, Lamington, and Merapi volcanoes. Our results compare reasonably well with independent estimates. For several cases, we found that it is possible to use the dynamic pressure equations developed for vertically uniform flow, along with the average cloud density multiplied by a factor of 2–5, to determine average velocity over the height of the failed object.     

An analogue model study of EM induction in the North China-Korea coastal region

The behaviour of the magnetic field variations over the North China-Korea coastal region is studied with the aid of a scaled laboratory analogue model. The model source frequencies simulate naturally occurring geomagnetic variations of 3–60 min periods. In-phase and quadrature magnetic B_x, B_y and B_z field measurements for the modelled region are presented for E- and B-polarizations. Large anomalous in-phase and quadrature model magnetic fields are observed over the Korea-Japan strait for E-polarization and over the Bohai strait for B-polarization due to current channelling through the straits. Large responses of the peninsulas in the shallow coastal areas occur at short periods but decrease abruptly with increasing period. In general, the effects of peninsulas, straits, bays and irregularities in the coastlines play an important role in the magnetic field responses both on-shore and off-shore for this complex North China-Korea coastal region. Model and field site induction arrows are compared for three sites west of Bohai Bay.     

Topographic effects on telluric field measurements

Using the Schwartz-Christofel transformation and numerical integration, the effect of a sloping topographic irregularity on the telluric field measurements in a sedimentary basin is estimated. Results show that in the vicinity of the topographic feature, the distortion introduced would increase with the angle of inclination of the sloping feature. It is noticed that, for moderate inclinations (20°–50°), the telluric field measured near the topographic feature is within 10% of its undisturbed value for distance greater than 0.1H–0.3H from the topographic feature, whereH is the thickness of the sedimentary column overlying the resistive basement. Suitable charts are prepared to aid as means to arrive at the estimates of errors for various angles of inclinations of such topographic feature and also to help in formulating approximate rules of thumb for selection of station sites in a field survey to minimize such topographic effects.     

On the interpretation of vertical gravity gradients produced by magmatic intrusions

A method for the computation and interpretation of gravity and height changes and vertical gravity gradients produced by magmatic intrusions in a layered elastic–gravitational medium is presented. The methodology assumes a planar medium geometry, which consists of welded layers overlying a half-space. The medium is elastic and gravitating. The intrusion is treated as a point source and can be located at any depth inside the medium. The theoretical elastic–gravitational model allows the computation of the so-called geometric and orthometric vertical displacements as well as gravity changes of different types. The corresponding vertical gravity gradients can also be computed. We present several examples of theoretical computations and study different types of these geometric and orthometric gravity gradients and the information we can get from the application of the methodology described. The results presented show that the use of these gradients is a useful tool to obtain information on the dynamics of the injection processes, including the detection of new magma recharge. We show that using the elastic–gravitational deformation model we can explain non-linear gravity–height relationships that appear in volcanic areas. We also present the application of the methodology to Mayon volcano, Philippines, delineating the intrusion of new magma, consistent with the produced eruptions after the observation period.     

Regional study of the anomalous change in apparent resistivity before the Tangshan earthquake (M=7.8, 1976) in China

Electrical resistivity measurements have been conducted as a possible means for obtaining precursory earthquake information. Before five great earthquakes (M>7,h<25 km) in China, the apparent resistivity _a showed systematic variations within a region 200 km from the epicenters. In particular, 9 stations in the Tangshan-Tianjin-Beijing region prior to the Tangshan earthquake (M=7.8,h=11 km, 27 July 1976) showed a consistent decrease of apparent resistivity around the epicenter, with a maximum resistivity change of 6% and a period of variation of 2–3 years. Simultaneous water table observations in this region showed a declining water table, and ground surface observations indicated a slight (5 mm) uplift in the epicenter region relative to its surroundings.In order to develop an explanation for the observed change of apparent resistivity associated with these great earthquakes, we have used Archie's Law to explore the effects of changes in rock porosity, water content and electrolyte resistivity on measured resistivity.Tentative conclusions of this study are as follows: (1) the apparent resistivity change is opposite to the effect expected from the simultaneous water table trend; (2) the dilatancy needed to give such resistivity variations (assuming Archie's Law holds) is much larger than that needed to explain the observed uplift (by 2–3 orders of magnitude); (3) salinity change in the pore electrolyte is a possible explanation for the variation in resistivity: an increase in the salinity would cause a proportional decrease in resistivity; the data needed to test this hypothesis, however, are lacking; and (4) the effect of changing geometry of rock pores or cracks due to pressure solution may provide an explanation for the decrease in apparent resistivity; it is different in nature from the effect of a volume change in response to stress although the geometry change is also closely related to the stress change.     

矿井瞬变电磁法全空间视电阻率解释方法研究

从全空间中心回线源响应公式出发,推导了全空间晚期和全区视电阻率公式;通过与半空间公式比较和对实测资料的分析,研究了全空间视电阻率公式的应用效果.结果显示,全空间视电阻率为半空间视电阻率的(5/2)^2/3倍,其晚期视电阻率与真值的相对误差小于0.56%,而相应的半空间视电阻率与真值误差均大于42%.与常规解释方法相比,采用全空间全区视电阻率换算方法能更好地圈定低阻异常范围,且电阻率差异大.因此,该公式更适用于矿井瞬变电磁法资料处理与解释,并具有更好的效果.     

Model‐based multi‐metric control of uniaxial shake tables

Shake tables provide a direct means by which to evaluate structural performance under earthquake excitation. Because the entire structure is mounted on the base plate and subjected to the ground motion in real time, dynamic effects and rate‐dependent behavior can be accurately represented. Shake table control is not straightforward as the desired signal is an acceleration record, while most actuators operate in displacement feedback for stability. At the same time, the payload is typically large relative to the capacity of the actuator, leading to pronounced control‐structure interaction. Through this interaction, the dynamics of the specimen influence the dynamics of the shake table, which can be problematic when specimens change behavior because of damage or other nonlinearities. Moreover, shake tables are themselves inherently nonlinear, making it difficult to accurately recreate a desired acceleration record over a broad range of amplitudes and frequencies. A model‐based multi‐metric shake table control strategy is proposed to improve tracking of the desired acceleration of a uniaxial shake table, remaining robust to nonlinearities including changes in specimen condition. The proposed strategy is verified for the shake table testing of both linear and nonlinear structures. Copyright © 2013 John Wiley & Sons, Ltd.     

Amplification effects from earthquakes and ambient noise in the Dead Sea rift (Israel)

Seven sites were instrumented in the Parsa area located in the seismically active Dead Sea rift system. Moderate and weak motions generated by earthquakes and ambient noise were used to identify amplifications due to geological and topographic effects.Three observation methods were applied to estimate site effects: (1) conventional soil–bedrock station-pair spectral ratios for earthquake motions and microtremors; (2) horizontal-to-vertical component spectral ratios for shear-waves observed simultaneously at a site (receiver function estimates) and (3) horizontal-to-vertical spectral ratios of microtremor measurements (Nakamura estimate). The site response spectra of soil sites exhibited significant peaks between 1 and 3 Hz with amplification factors typically within the range of 2.5–4.0. A bedrock site on the high plateau near the escarpment top showed a peak between 2 and 3 Hz, mainly due to an EW oscillation of the NS topographic feature. Our observations indicated that seismograms recorded in the tunnel were either enriched or depleted at certain frequencies owing to interference of incident and surface-reflected waves.     

Are topographic data sets fractal?

The scale invariant properties of fractal sets make them attractive models for topographic profiles because those profiles are the end product of a complex system of physical processes operating over many spatial scales. If topographic data sets are fractal, their power spectra will be well represented by lines in log-log space with slopess such that –3s<–1. The power spectra from a Digital Elevation Model (30 meter sample spacing) of the Sierra Nevada Batholith and from Seabeam center beam depths (425 meter sample spacing) along a flowline in the South Atlantic are curved. Straight sections in the spectra can be identified but the slopes of those sections are strongly dependent upon the particulars of the data analysis. Fractal geometry must be used with caution in the discussion of topographic data sets.     

Effects of polarization calibration on aerosol optical depth retrieval: An ocean case sensitivity analysis

Reflectance measurements of both the visible and infrared bands of passive remote sensing sensors are widely used to retrieve aerosol optical depth(AOD) information. This is performed commonly for data obtained over both ocean and land, and these measurements allow for the off line development of a lookup table using radiative transfer models. Owing to molecular and aerosol effects, the reflected light received by the sensor is usually highly polarized. The linear polarization effect may be up to 100%, and the polarization factor of a sensor optical system will change the total intensity as well as the polarization status of the signal reaching the detector. The detector response will be different when the incident light polarization status changes, even if the total intensity remains constant. However, if the polarization calibration is neglected, it will cause obvious errors in the aerosol data retrieval. This is especially true for aerosol optical depth retrieval over an ocean. This measurement relies directly on the reflectance output of the sensor. Cases involving land surfaces are not discussed herein because the inhomogeneous properties conceal the error due to polarization. Taking the 550 and 860 nm bands as examples, the difference between the real top-of-atmosphere(TOA) reflectance and the reflectance reaching the detector is calculated using three different sensor polarization standards according to the Sea-viewing Wide Field-of-view Sensor(Sea Wi FS) and Moderate Resolution Imaging Spectroradiometer(MODIS) standards. The differences in AOD retrieval are also demonstrated using the lookup table developed previously from a vector radiative transfer code. The results reveal that under a normal situation in which the AOD is 0.15, the maximum AOD retrieval error could reach 0.04 in 550 nm but only 0.02 in 860 nm for the dust aerosol model. For the soot aerosol model, the maximum AOD retrieval error is 0.1 in 550 nm and 0.12 in 860 nm, indicating that the lack of polarization calibration will lead to large errors in aerosol retrieval over an ocean.     

Geology of the quaternary volcanic centres of the east Anatolia

Following the collision along the Bitlis–Zagros suture, a north–south convergence between the Arabian Platform and Laurasia has continued uninterrupted until the present. As a result, the continental crust has been shortened, thickened and consequently elevated to form the Turkish–Iranian high plateau. On the high plateau volcanic activity began during the Neogene, intensified during the late Miocene–Pliocene and continued until historical times. Large volcanic centres have been developed during the Quaternary which form significant peaks above the Turkish–Iranian high plateau. Among the Quaternary volcanoes, the major volcanic centres are Ararat, Tendürek, Suphan and Nemrut. Ararat (Ağri Daği) is the largest volcanic center and is a compound stratovolcano, consisting of Greater Ararat and lesser Ararat. The former represents the highest elevation of Anatolia reaching over 5000 m in height. Tendürek is a double-peaked shield volcano, which produced a voluminous amount of basalt lava as extensive pahoehoe, and aa flows. It has an ill-defined semi-caldera. Suphan is an isolated stratovolcano, capped by silicic dome. It represents the second highest topographic elevation in Anatolia, with a height of over 4000 m. A cluster of subsidiary cones and small domes surrounds the volcano. Nemrut is the largest member of a group of volcanoes, which trend north–south. It is a stratovolcano, having a well-defined collapse caldera and a caldera lake. Various volcanic ejecta have been extruded from these volcanic centres over the last 1 to 2 million years. The Quaternary volcanic centres, although temporally and spatially closely associated, display a wide range of lavas from basalt to rhyolite. The volcanoes have diverse compositional trends; Ararat is distinctly subalkaline, Suphan is mildly subalkaline, Nemrut is mildly alkaline and Tendürek is strongly alkaline. The major and trace element compositions together with the isotope ratios indicate that their magmas were generated from a heterogeneous mantle source. Each of the volcanic centres has undergone a partly different magmatic evolution.     

Geoelectromagnetic and geothermic investigations in the Ihlara Valley geothermal field

The Ihlara Valley is situated within a volcanic arc that is formed by the collision of the eastern Mediterranean plate system with the Anatolian plate. In this study we will present data from a reservoir monitoring project over the Ihlara-Ziga geothermal field, located 22 km east of Aksaray, in central Anatolia.Although identified geothermal resources in the Ihlara Valley are modest, substantial undiscovered fields have been inferred primarily from the volcanic and tectonic setting but also from the high regional heat flow (150–200 mWm⁻²) on the Kir ehir Massif.In 1988 and 1990, geoelectromagnetic surveys were undertaken by MTA-Ankara to confirm the presence of a relatively shallow (≈ 0.5–1 km), hydrothermally caused conductive layer or zone. CSAMT and Schlumberger resistivity data show good correspondence with each other, and 2-D geoelectric models are also in harmony with geologic data and gravity anomalies.The depth of the resistive basement, which is interpreted as Paleozoic limestone, is 200–250 m in the western part and increases eastward (≈ 600–750 m). This may imply N-S-oriented normal faulting within the survey area. The parameters of the top layer are a resistivity of 25 to 95 ohm m and a thickness of between 100 and 250 m. The thickness of the conductive tuffs between the top layer and the basement, whose resistivity is about 4–5 o hmm, also increases eastward (from 100 to 450 m). The apparent resistivity maps for the frequencies between 32 and 2 Hz reveal a localized low resistivity anomaly to the east of Belisirma.