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

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

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

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

Analyzing and modeling gravity and magnetic anomalies using the SPHERE program and Magsat data

Computer codes were completed, tested, and documented for analyzing magnetic anomaly vector components by equivalent point dipole inversion. The codes are intended for use in inverting the magnetic anomaly due to a spherical prism in a horizontal geomagnetic field and for recomputing the anomaly in a vertical geomagnetic field. Modeling of potential fields at satellite elevations that are derived from three dimensional sources by program SPHERE was made significantly more efficient by improving the input routines. A preliminary model of the Andean subduction zone was used to compute the anomaly at satellite elevations using both actual geomagnetic parameters and vertical polarization. Program SPHERE is also being used to calculate satellite level magnetic and gravity anomalies from the Amazon River Aulacogen

Study of gravity and magnetic anomalies using MAGSAT data

The results of modeling satellite-elevation magnetic and gravity data using the constraints imposed by near surface data and seismic evidence shows that the magnetic minimum can be accounted for by either an intracrustal lithologic variation or by an upwarp of the Curie point isotherm. The long wavelength anomalies of the NOO's-vector magnetic survey of the conterminous U.S. were contoured and processed by various frequency filters to enhance particular characteristics. A preliminary inversion of the data was completed and the anomaly field calculated at 450 km from the equivalent magnet sources to compare with the POGO satellite data. Considerable progress was made in studing the satellite magnetic data of South America and adjacent marine areas. Preliminary versions of the 1 deg free-air gravity anomaly map (20 m gal contour interval) and the high cut (lambda approximately 8 deg) filtered anomaly maps are included

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