Recovering magnetic susceptibility from electromagnetic data over a one-dimensional earth

While the inversion of electromagnetic data to recover electrical conductivity has received much attention, the inversion of those data to recover magnetic susceptibility has not been fully studied. In this paper we invert frequency-domain electromagnetic (EM) data from a horizontal coplanar system to recover a 1-D distribution of magnetic susceptibility under the assumption that the electrical conductivity is known. The inversion is carried out by dividing the earth into layers of constant susceptibility and minimizing an objective function of the susceptibility subject to fitting the data. An adjoint Green's function solution is used in the calculation of sensitivities, and it is apparent that the sensitivity problem is driven by three sources. One of the sources is the scaled electric field in the layer of interest, and the other two, related to effective magnetic charges, are located at the upper and lower boundaries of the layer. These charges give rise to a frequency-independent term in the sensitivities. Because different frequencies penetrate to different depths in the earth, the EM data contain inherent information about the depth distribution of susceptibility. This contrasts with static field measurements, which can be reproduced by a surface layer of magnetization. We illustrate the effectiveness of the inversion algorithm on synthetic and field data and show also the importance of knowing the background conductivity. In practical circumstances, where there is no a priori information about conductivity distribution, a simultaneous inversion of EM data to recover both electrical conductivity and susceptibility will be required.     

Palaeosecular variation over the past 5 Myr based on a new generalized database

A generalized database of global palaeomagnetic data from 3719 lava flows and thin dykes of age 0–5 Ma has been constructed for use with a relational database management system. The database includes all data whose virtual geomagnetic poles (VGP) lie within 45 of the spin axis and can be used for studies of palaeosecular variation and for geomagnetic field modelling. Because many of these data were collected and processed more than 15–20 years ago, each result has been characterized according to the demagnetization procedures carried out. Analysis of these data in terms of the latitude variation of the angular dispersion of VGPs (palaeosecular variation from lavas) strongly suggests that careful data selection is required and that many of the older studies may need to be redone using more modern methods. Differences between the angular dispersions for separate normal- and reverse-polarity data sets confirm that many older studies have not been adequately cleaned magnetically. Therefore, the use of the database for geomagnetic field modelling should be carried out with some caution. Using a VGP cut-off angle that varies with latitude, the best data set consists of 2636 results that show a smooth increase of VGP angular dispersion with latitude. Model G for palaeosecular variation, which is based on modelling of the antisymmetric (dipole) and symmetric (quadrupole) dynamo families, provides a good fit to these results.     

Climatic Conditions for Wine Grape Growing

Vineyards in New Zealand are pushing further south than anywhere else in the Southern Hemisphere, this despite poor understanding of the latitudinal limits of grape growing. European experience shows a limit of 50^oN. Valuable pointers and cautionary tales from Germany encourage a view that the southern limit in New Zealand may be 45^oS. This is broadly confirmed by comparative climate data from Germany and New Zealand.     

Distribution and morphology of sinkholes triggered by flooding following Tropical Storm Alberto at Albany, Georgia, USA

Flooding of the Flint River in July 1994 triggered the collapse of at least 312 sinkholes in the karstic Dougherty Plain at Albany, Georgia. We examined the distribution and morphology of these new sinkholes to evaluate the mode of formation, to characterize early stages of the evolution of sinkhole form, and to estimate the lowering of the surface associated with the development of new sinkholes.Eighty-eight percent of sinkholes occur inside the limits of flooding, especially in areas of sandy overburden, and they often follow joint-controlled linear trends. Sinkhole dimensions are log-normally distributed with median values of circumference = 5.7 m, length = 1.8 m, width = 1.6 m, and depth = 0.7 m; asymmetry (L:W) = 1.2. Cross-sectional forms range from narrow cylinders to large bowls, with many sinkholes having undercut sides.Flooding triggered the formation of sinkholes by saturating and liquefying overburden, which caused soil arches to collapse and flow into cavities in bedrock. The prevalence of sinkholes near the periphery of flooding suggests that drainage and loss of buoyant support as flood waters subsided may also have contributed to failure. A volume ratio index is used to quantify the three dimensional geometric form of sinkholes. Initially, small cylindrical shafts open over a bedrock joint, followed by progressive slumping that leads to widening and increases in volume to a final bowl form. Estimates of the aggregate volume of overburden transported underground in flooded areas range from 7,990 to 11,130 m3. Averaged over flooded areas, this accounts for 0.26 to 0.37 mm/km² lowering of the surface. Based on a 500 year recurrence interval for the flood event, values for lowering of the surface range from 0.52 to 0.74 mm per 1,000 years. These values are an order of magnitude less than estimates of carbonate dissolution and suggest that transport of overburden underground is limited by triggering events.