MAGSAT scalar anomalies

To facilitate processing large data arrays, elements of spherical Earth analysis programs NVERTSM, SMFLD, NVERTG and GLFD were implemented and tested on the LARS IBM 4341 computer. Currently, the problem of inverting 2 deg MAGSAT scalar anomalies for the region (80 W, 60 E) longitude and (40 S, 70 N) latitude is being implemented on the LARS-computer for quantitative comparison with free air gravity anomaly, geothermal and tectonic data. Gravity and MAGSAT anomalies from a subset of this region (30 W, 60 E), (40 S, 70 N) were already processed for a paper on satellite magnetic anomalies of Africa and Europe

Satellite elevation magnetic anomaly maps

The problem of inverting 2 deg average MAGSAT scalar anomalies for the region 80 W, 60 E longitude and 40 S, 70 N latitude was attempted on the LARS computer; however, the effort was aborted due to insufficient allocation of CPU-time. This problem is currently being resubmitted and should be implemented shortly for quantitative comparison with free-air gravity anomaly, geothermal, and tectonic data

The south-central United States magnetic anomaly

The South-Central United States Magnetic Anomaly is the most prominent positive feature in the MAGSAT scalar magnetic field over North America. The anomaly correlates with increased crustal thickness, above average crustal velocity, negative free air gravity anomalies and an extensive zone of Middle Proterozoic anorogenic felsic basement rocks. Spherical dipole source inversion of the MAGSAT scalar data and subsequent calculation of reduced to pole and derivative maps provide constraints for a crustal magnetic model which corresponds geographically to the extensive Middle Proterozoic felsic rocks trending northeasterly across the United States. These felsic rocks contain insufficient magnetization or volume to produce the anomaly, but are rather indicative of a crustal zone which was disturbed during a Middle Proterozoic thermal event which enriched magnetic material deep in the crust

Regularization-robust preconditioners for time-dependent PDE constrained optimization problems

In this article, we motivate, derive and test �effective preconditioners to be used with the Minres algorithm for solving a number of saddle point systems, which arise in PDE constrained optimization problems. We consider the distributed control problem involving the heat equation with two diff�erent functionals, and the Neumann boundary control problem involving Poisson's equation and the heat equation. Crucial to the eff�ectiveness of our preconditioners in each case is an eff�ective approximation of the Schur complement of the matrix system. In each case, we state the problem being solved, propose the preconditioning approach, prove relevant eigenvalue bounds, and provide numerical results which demonstrate that our solvers are eff�ective for a wide range of regularization parameter values, as well as mesh sizes and time-steps

Fast iterative solution of reaction-diffusion control problems arising from chemical processes

PDE-constrained optimization problems, and the development of preconditioned iterative methods for the efficient solution of the arising matrix system, is a field of numerical analysis that has recently been attracting much attention. In this paper, we analyze and develop preconditioners for matrix systems that arise from the optimal control of reaction-diffusion equations, which themselves result from chemical processes. Important aspects in our solvers are saddle point theory, mass matrix representation and effective Schur complement approximation, as well as the outer (Newton) iteration to take account of the nonlinearity of the underlying PDEs

Utility of correlation techniques in gravity and magnetic interpretation

Internal correspondence uses Poisson's Theorem in a moving-window linear regression analysis between the anomalous first vertical derivative of gravity and total magnetic field reduced to the pole. The regression parameters provide critical information on source characteristics. The correlation coefficient indicates the strength of the relation between magnetics and gravity. Slope value gives delta j/delta sigma estimates of the anomalous source. The intercept furnishes information on anomaly interference. Cluster analysis consists of the classification of subsets of data into groups of similarity based on correlation of selected characteristics of the anomalies. Model studies are used to illustrate implementation and interpretation procedures of these methods, particularly internal correspondence. Analysis of the results of applying these methods to data from the midcontinent and a transcontinental profile shows they can be useful in identifying crustal provinces, providing information on horizontal and vertical variations of physical properties over province size zones, validating long wavelength anomalies, and isolating geomagnetic field removal problems

Utility of correlation techniques in gravity and magnetic interpretation

Two methods of quantitative combined analysis, internal correspondence and clustering, are presented. Model studies are used to illustrate implementation and interpretation procedures of these methods, particularly internal correspondence. Analysis of the results of applying these methods to data from the midcontinent and a transcontinental profile show they can be useful in identifying crustal provinces, providing information on horizontal and vertical variations of physical properties over province size zones, validating long wave-length anomalies, and isolating geomagnetic field removal problems. Thus, these techniques are useful in considering regional data acquired by satellites

Spherical Earth analysis and modeling of lithospheric gravity and magnetic anomalies

A comprehensive approach to the lithospheric analysis of potential field anomalies in the spherical domain is provided. It has widespread application in the analysis and design of satellite gravity and magnetic surveys for geological investigation

Gravity and magnetic anomaly modeling and correlation using the SPHERE program and Magsat data

The spherical Earth inversion, modeling, and contouring software were tested and modified for processing data in the Southern Hemisphere. Preliminary geologic/tectonic maps and selected cross sections for South and Central America and the Caribbean region are being compiled and as well as gravity and magnetic models for the major geological features of the area. A preliminary gravity model of the Andeas Beniff Zone was constructed so that the density columns east and west of the subducted plates are in approximate isostatic equilibrium. The magnetic anomaly for the corresponding magnetic model of the zone is being computed with the SPHERE program. A test tape containing global magnetic measurements was converted to a tape compatible with Purdue's CDC system. NOO data were screened for periods of high diurnal activity and reduced to anomaly form using the IGS-75 model. Magnetic intensity anomaly profiles were plotted on the conterminous U.S. map using the track lines as the anomaly base level. The transcontinental magnetic high seen in POGO and MAGSAT data is also represented in the NOO data

Gravity and magnetic anomaly data analysis

Progress on the analysis MAGSAT data is reported. The MAGSAT data from 40 deg S to 70 deg N latitude and 30 deg W to 60 E longitude was reduced to radial polarization. In addition, gravity anomaly data from this area were processed and a variety of filtered maps were prepared for combined interpretation of the gravity and magnetic data in conjunction with structural and tectonic maps of the area. The VERSATEC listings and cross-reference maps of variable and array names for the spherical Earth analysis programs NVERTSM, SMFLD, NVERTG, and GFLD were also prepared