Engineering topics

Voorhies, Coerte V.

Publications and source records attributed to Voorhies, Coerte V..

18 records

NASA NTRS2006

Thickness of the Magnetic Crust of Mars from Magneto-Spectral Analysis

Previous analysis of the magnetic spectrum of Mars showed only a crustal source field. The observational spectrum was fairly well fitted by the spectrum expected from random dipolar sources scattered on a spherical shell about 46 plus or minus 10 km below Mars' 3389.5 km mean radius. This de-correlation depth overestimates the typical depth of extended magnetized structures, and so was judged closer to mean source layer thickness than twice its value. To better estimate the thickness of the magnetic crust of Mars, six different magnetic spectra were fitted with the theoretical spectrum expected from a novel, bimodal distribution of magnetic sources. This theoretical spectrum represents both compact and extended, laterally correlated sources, so source shell depth is doubled to obtain layer thickness. The typical magnetic crustal thickness is put at 47.8 plus or minus 8.2 km. The extended sources are enormous, typically 650 km across, and account for over half the magnetic energy at low degrees. How did such vast regions form?

Voorhies, Coerte V.↗

NASA NTRS2004

Narrow Scale Flow and a Weak Field by the Top of Earth's Core: Evidence from Orsted, Magsat and Secular Variation

As Earth's main magnetic field weakens, our magnetic shield against the onslaught of the solar wind thins. And the field strength needed to fend off battering by solar coronal mass ejections is decreasing, just when the delicate complexity of modem, vulnerable, electro-technological systems is increasing at an unprecedented rate. Recently, a working group of distinguished scientist from across the nation has asked NASA's Solid Earth and Natural Hazards program a key question: What are the dynamics of Earth s magnetic field and its interactions with the Earth system? Paleomagnetic studies of crustal rocks magnetized in the geologic past reveal that polarity reversals have occurred many times during Earth s history. Networked super-computer simulations of core field and flow, including effects of gravitational, pressure, rotational Coriolis, magnetic and viscous forces, suggest how this might happen in detail. And space-based measurements of the real, time-varying magnetic field help constrain estimates of the speed and direction of fluid iron flowing near the top of the core and enable tests of some hypotheses about such flow. Now scientists at NASA s Goddard Space Flight Center have developed and applied methods to test the hypotheses of narrow scale flow and of a dynamically weak magnetic field near the top of Earth s core. Using two completely different methods, C. V. Voorhies has shown these hypotheses lead to specific theoretical forms for the "spectrum" of Earth s main magnetic field and the spectrum of its rate of change. Much as solar physicists use a prism to separate sunlight into its spectrum, from long wavelength red to short wavelength blue light, geophysicists use a digital prism, spherical harmonic analysis, to separate the measured geomagnetic field into its spectrum, from long to short wavelength fields. They do this for the rate of change of the field as well.

Voorhies, Coerte V.↗

NASA NTRS2004

Crustal Magnetic Spectra from Correlated Sources on Mars

The spectral method for distinguishing crustal from core-source magnetic fields has been re-examined, modified and applied to both a comprehensive geomagnetic field model and an altitude normalized magnetic map of Mars. These observational spectra are fairly fitted by theoretical forms expected from certain elementary classes of magnetic sources. For Earth we found fields from a core radius 3512 +/- 64 km and a crust represented by a shell of random dipolar sources at radius 6367 +/- 14 km. For Mars we found only a field from a crust represented in the same way, but 46 +/- 10 km below the planetary menu radius. More realistic theoretical spectra, allowing for crustal thickness, oblateness and magnetization by a planet centered dipole, were derived and discussed, as were spectral effects of laterally correlated sources. The main effect of laterally correlated sources is to soften the spectrum at high degrees. We tend to over-estimate source shell depth when this is omitted. To include this effect simply size and magnetization distribution functions for extended sources are recast as a characteristic diameter and mean square magnetization amplitude for an ensemble of vertically magnetized spherical caps on a shell. For small caps, and at moderated degrees, the practical derivatives of the log-theoretical spectrum with respect to amplitude, shell radius, and cap diameter are approximately proportional to l, n, and -n **2, respectively. Separation of diameter from amplitude from layer thickness. Results from applications to observational spectra are discussed, noting that there are now several fine field models for Mars; moreover, the terrestrial magnetic spectrum at high degrees, as revised, updated with high precision Oerstad data, and upgrade with high resolution Champ data, appears softer than before.

Voorhies, Coerte V.↗

NASA NTRS2004

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothesis passes the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale range. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length-scale 1/k, if one presumes kinetic energy is converted in either eddy-overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small-scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magnetogeostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as a barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

NASA NTRS2004

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynmcal hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothes pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at lentgh-scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto-geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core- mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order l20 years is pseudo-scale-invarient. Magnetic spectra of other planets may differ; however, if a transition to non-conductmg fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

NASA NTRS2004

A Geomagnetic Estimate of Mean Paleointensity

To test a statistical hypothesis about Earth's magnetic field against paleomagnetism, the present field is used to estimate time averaged paleointensity. The estimate used the modern magnetic multipole spectrum R(n), which gives the mean square induction represented by spherical harmonics of degree n averaged over the sphere of radius a = 6371.2 km. The hypothesis asserts that low degree multi-pole powers of the coresource field are distributed as chi-squared with 2n+1 degrees of freedom and expectation values, where c is the 3480 km radius of the Earth's core. (This is compatible with a usually mainly geocentric axial dipolar field). Amplitude K is estimated by fitting theoretical to observational spectra through degree 12. The resulting calibrated expectation spectrum is summed through degree 12 to estimate expected square intensity F(exp 2). The sum also estimates F(exp 2) averaged over geologic time, in so far as the present magnetic spectrum is a fair sample of that generated in the past by core geodynamic processes. Additional information is included in the original extended abstract.

Voorhies, Coerte V.↗

NASA NTRS2000

On Geomagnetism and Paleomagnetism I

A partial description of Earth's broad scale, core-source magnetic field has been developed and tested three ways. The description features an expected, or mean, spatial magnetic power spectrum that is approximately inversely proportional to horizontal wavenumber atop Earth's core. This multipole spectrum describes a magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Temporal variations of core multipole powers about mean values are to be expected and are described statistically, via trial probability distribution functions, instead of deterministically, via trial solution of closed transport equations. The distributions considered here are closed and neither require nor prohibit magnetic isotropy. The description is therefore applicable to, and tested against, both dipole and low degree non-dipole fields. In Part 1, a physical basis for an expectation spectrum is developed and checked. The description is then combined with main field models of twentieth century satellite and surface geomagnetic field measurements to make testable predictions of the radius of Earth's core. The predicted core radius is 0.7% above the 3480 km seismological value. Partial descriptions of other planetary dipole fields are noted.

Voorhies, Coerte V.↗

NASA NTRS1998

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aid core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile, akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core; (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core relative to the mantle is calculated to be at most 1.5 deg./yr.

Voorhies, Coerte V.↗

NASA NTRS1998

On Geomagnetism and Paleomagnetism

A statistical description of Earth's broad scale, core-source magnetic field has been developed and tested. The description features an expected, or mean, spatial magnetic power spectrum that is neither "flat" nor "while" at any depth, but is akin to spectra advanced by Stevenson and McLeod. This multipole spectrum describes the magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Natural variations of core multipole powers about their mean values are to be expected over geologic time and are described via trial probability distribution functions that neither require nor prohibit magnetic isotropy. The description is thus applicable to core-source dipole and low degree non-dipole fields despite axial dipole anisotropy. The description is combined with main field models of modem satellite and surface geomagnetic measurements to make testable predictions of: (1) the radius of Earth's core, (2) mean paleomagnetic field intensity, and (3) the mean rates and durations of both dipole power excursions and durable axial dipole reversals. The predicted core radius is 0.7% above the 3480 km seismologic value. The predicted root mean square paleointensity (35.6 mu T) and mean Virtual Axial Dipole Moment (about 6.2 lx 1022 Am(exp 2)) are within the range of various mean paleointensity estimates. The predicted mean rate of dipole power excursions, as defined by an absolute dipole moment <20% of the 1980 value, is 9.04/Myr and 14% less than obtained by analysis of a 4 Myr paleointensity record. The predicted mean rate of durable axial dipole reversals (2.26/Myr) is 2.3% more than established by the polarity time-scale for the past 84 Myr. The predicted mean duration of axial dipole reversals (5533 yr) is indistinguishable from an observational value. The accuracy of these predictions demonstrates the power and utility of the description, which is thought to merit further development and testing. It is suggested that strong stable stratification of Earth's uppermost outer core leads to a geologically long interval of no dipole reversals and a very nearly axisymmetric field outside the core. Statistical descriptions of other planetary magnetic fields are outlined.

Voorhies, Coerte V.↗

NASA NTRS1996

Accurate Predictions of Mean Geomagnetic Dipole Excursion and Reversal Frequencies, Mean Paleomagnetic Field Intensity, and the Radius of Earth's Core Using McLeod's Rule

The geomagnetic spatial power spectrum R(sub n)(r) is the mean square magnetic induction represented by degree n spherical harmonic coefficients of the internal scalar potential averaged over the geocentric sphere of radius r. McLeod's Rule for the magnetic field generated by Earth's core geodynamo says that the expected core surface power spectrum (R(sub nc)(c)) is inversely proportional to (2n + 1) for 1 less than n less than or equal to N(sub E). McLeod's Rule is verified by locating Earth's core with main field models of Magsat data; the estimated core radius of 3485 kn is close to the seismologic value for c of 3480 km. McLeod's Rule and similar forms are then calibrated with the model values of R(sub n) for 3 less than or = n less than or = 12. Extrapolation to the degree 1 dipole predicts the expectation value of Earth's dipole moment to be about 5.89 x 10(exp 22) Am(exp 2)rms (74.5% of the 1980 value) and the expected geomagnetic intensity to be about 35.6 (mu)T rms at Earth's surface. Archeo- and paleomagnetic field intensity data show these and related predictions to be reasonably accurate. The probability distribution chi(exp 2) with 2n+1 degrees of freedom is assigned to (2n + 1)R(sub nc)/(R(sub nc). Extending this to the dipole implies that an exceptionally weak absolute dipole moment (less than or = 20% of the 1980 value) will exist during 2.5% of geologic time. The mean duration for such major geomagnetic dipole power excursions, one quarter of which feature durable axial dipole reversal, is estimated from the modern dipole power time-scale and the statistical model of excursions. The resulting mean excursion duration of 2767 years forces us to predict an average of 9.04 excursions per million years, 2.26 axial dipole reversals per million years, and a mean reversal duration of 5533 years. Paleomagnetic data show these predictions to be quite accurate. McLeod's Rule led to accurate predictions of Earth's core radius, mean paleomagnetic field intensity, and mean geomagnetic dipole power excursion and axial dipole reversal frequencies. We conclude that McLeod's Rule helps unify geo-paleomagnetism, correctly relates theoretically predictable statistical properties of the core geodynamo to magnetic observation, and provides a priori information required for stochastic inversion of paleo-, archeo-, and/or historical geomagnetic measurements.

Voorhies, Coerte V.↗

NASA NTRS1993

Steady induction effects in geomagnetism. Part 1C: Geomagnetic estimation of steady surficial core motions: Application to the definitive geomagnetic reference field models

In the source-free mantle/frozen-flux core magnetic earth model, the non-linear inverse steady motional induction problem was solved using the method presented in Part 1B. How that method was applied to estimate steady, broad-scale fluid velocity fields near the top of Earth's core that induce the secular change indicated by the Definitive Geomagnetic Reference Field (DGRF) models from 1945 to 1980 are described. Special attention is given to the derivation of weight matrices for the DGRF models because the weights determine the apparent significance of the residual secular change. The derived weight matrices also enable estimation of the secular change signal-to-noise ratio characterizing the DGRF models. Two types of weights were derived in 1987-88: radial field weights for fitting the evolution of the broad-scale portion of the radial geomagnetic field component at Earth's surface implied by the DGRF's, and general weights for fitting the evolution of the broad-scale portion of the scalar potential specified by these models. The difference is non-trivial because not all the geomagnetic data represented by the DGRF's constrain the radial field component. For radial field weights (or general weights), a quantitatively acceptable explication of broad-scale secular change relative to the 1980 Magsat epoch must account for 99.94271 percent (or 99.98784 percent) of the total weighted variance accumulated therein. Tolerable normalized root-mean-square weighted residuals of 2.394 percent (or 1.103 percent) are less than the 7 percent errors expected in the source-free mantle/frozen-flux core approximation.

Voorhies, Coerte V.↗

NASA NTRS1993

Simultaneous solution for core magnetic field and fluid flow beneath an electrically conducting mantle

The effects of laterally homogeneous mantle electrical conductivity were included in steady, frozen-flux core surface flow estimation along with refinements in method and weighting. The refined method allows simultaneous solution for both the initial radial geomagnetic field component at the core-mantle boundary (CMB) and the sub-adjacent fluid motion; it also features Gauss' method for solving the non-linear inverse problem associated with steady motional induction. The tradeoff between spatial complexity of the derived flows and misfit to the weighted Definitive Geomagnetic Reference Field models (DGRF's) is studied for various mantle conductivity profiles. For simple flow and a fixed initial geomagnetic condition, a fairly high deep-mantle conductivity performs better than either insulating or weakly conducting profiles; however, a thin, very high conductivity layer at the base of the mantle performs almost as well. Simultaneous solution for both initial geomagnetic field and flow reduces the misfit per degree of freedom even more than does changing the mantle conductivity profile. Moreover, when both core field and flow are estimated, the performance of the solutions and the derived flows become insensitive to the conductivity profile.

Voorhies, Coerte V.↗

NASA NTRS1993

Steady induction effects in geomagnetism. Part 1B: Geomagnetic estimation of steady surficial core motions: A non-linear inverse problem

The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation (SV) indicated by models of the observed geomagnetic field is examined in the source-free mantle/frozen-flux core (SFI/VFFC) approximation. This inverse problem is non-linear because solutions of the forward problem are deterministically chaotic. The SFM/FFC approximation is inexact, and neither the models nor the observations they represent are either complete or perfect. A method is developed for solving the non-linear inverse motional induction problem posed by the hypothesis of (piecewise, statistically) steady core surface flow and the supposition of a complete initial geomagnetic condition. The method features iterative solution of the weighted, linearized least-squares problem and admits optional biases favoring surficially geostrophic flow and/or spatially simple flow. Two types of weights are advanced radial field weights for fitting the evolution of the broad-scale portion of the radial field component near Earth's surface implied by the models, and generalized weights for fitting the evolution of the broad-scale portion of the scalar potential specified by the models.

Voorhies, Coerte V.↗

NASA NTRS1992

Steady induction effects in geomagnetism. Part 1A: Steady motional induction of geomagnetic chaos

Geomagnetic effects of magnetic induction by hypothetically steady fluid motion and steady magnetic flux diffusion near the top of Earth's core are investigated using electromagnetic theory, simple magnetic earth models, and numerical experiments with geomagnetic field models. The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation indicated by broad-scale models of the observed geomagnetic field is examined and solved. In Part 1, the steady surficial core flow estimation problem is solved in the context of the source-free mantle/frozen-flux core model. In the first paper (IA), the theory underlying such estimates is reviewed and some consequences of various kinematic and dynamic flow hypotheses are derived. For a frozen-flux core, fluid downwelling is required to change the mean square normal magnetic flux density averaged over the core-mantle boundary. For surficially geostrophic flow, downwelling implies poleward flow. The solution of the forward steady motional induction problem at the surface of a frozen-flux core is derived and found to be a fine, easily visualized example of deterministic chaos. Geomagnetic effects of statistically steady core surface flow may well dominate secular variation over several decades. Indeed, effects of persistent, if not steady, surficially geostrophic core flow are described which may help explain certain features of the present broad-scale geomagnetic field and perhaps paleomagnetic secular variation.

Voorhies, Coerte V.↗

NASA NTRS1991

On the joint inversion of geophysical data for models of the coupled core-mantle system

Joint inversion of magnetic, earth rotation, geoid, and seismic data for a unified model of the coupled core-mantle system is proposed and shown to be possible. A sample objective function is offered and simplified by targeting results from independent inversions and summary travel time residuals instead of original observations. These data are parameterized in terms of a very simple, closed model of the topographically coupled core-mantle system. Minimization of the simplified objective function leads to a nonlinear inverse problem; an iterative method for solution is presented. Parameterization and method are emphasized; numerical results are not presented.

Voorhies, Coerte V.↗

NASA NTRS1987

The time-varying geomagnetic field

Information (obtained from Magsat and other sources) published by the U.S. authors between 1983 and 1986 on the large-scale geomagnetic field of deep internal origin and its secular variation (SV) is reviewed. Results on the main field modeling, including reference fields, and on the separaton of core and crustal fields are discussed together with the advances made in geomagnetic field theory and applications of novel methods. Consideration is also given to global SV, short-term global SV, and regional CV analyses; nonsecular impulses, featuring the sudden geomagnetic jerk of 1969; and the frozen-flux core approximation, including downward continuation to the core-mantle boundary constraints, westward drift, and estimates of the fluid velocity at the top of the core.

Voorhies, Coerte V.↗

NASA NTRS1986

Steady surficial core motions - An alternate method

A new method for deriving steady surficial core motions from geomagnetic field models is presented and applied. The original method determines the steady velocity field at the top of a frozen-flux core which best fits, in both the spatial and temporal linear least squares sense, a model of the large scale geomagnetic secular variation (SV) at the base of a source-free mantle (Voorhies, 1986). New method is based upon the same physical assumptions, but fits SV at earth's surface instead of the core-mantle boundary (CMB). The newly derived flow differs somewhat from prior solutions, and shows but slight westward drift; however, other key global properties are virtually unchanged. Over a 15 year interval the new method provides a better fit to the time-varying geomagnetic field at earth's surface than does the old.

Voorhies, Coerte V.↗

NASA NTRS1986

Steady flows at the top of earth's core derived from geomagnetic field models

Select models of the main geomagnetic field and its secular variation are extrapolated to the base of an insulating mantle and used to estimate the adjacent fluid motion of a perfectly conducting outer core. The assumption of steady motion provides formally unique solutions and is tested along with that of no upwelling. The hypothesis of no upwelling is found to be substantially worse than that of steady motion. Although the actual motion is not thought to be steady, the large-scale secular variation at the top of the core can be adequately described by a large-scale, combined toroidal-poloidal circulation which is steady for intervals of at least a decade or two. The derived flows include a bulk westward drift but are complicated by superimposed jets, gyres, and surface divergence indicative of vigorous vertical motion at depth. The circulation pattern and key global properties including rms speed, upwelling, and westward drift are found to be fairly insensitive to variations in modeling parameters.

Voorhies, Coerte V.↗