Integral and spherical-harmonic analyses of the geomagnetic field for 1955.0, part i
Analysis of geomagnetic field for 1955 by spherical harmonics
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Analysis of geomagnetic field for 1955 by spherical harmonics
Use of Doppler tracking data from the Lunar Orbiter series of spacecraft in an analysis of the spherical harmonic coefficients of the lunar gravitational field through thirteenth degree and order. The value obtained for the mass of the moon, GM = 4902.84 cu km/sec/sec, is in good agreement with previous results and with results obtained by alternate procedures. Acceleration contour plots, derived from the gravitational coefficients, show correlations with surface features on the near side of the moon, but are of questionable validity for the far side because of the lack of direct tracking data on the far side. Based on the most recent gravitational field data, the current estimate for the polar moment of inertia of the moon is C/Ma squared = 0.4019 super + 0.004 sub - 0.002. This value indicates that the interior of the moon can be homogeneous, but some results presented strongly suggest that the moon is differentiated, with an excess of mass in the direction toward the earth.
The effect of a standing wave field in a rectangular chamber on the shape and location of an acoustically positioned drop or bubble is calculated. The sample deformation and equilibrium position are obtained from an analysis of the spherical harmonic projections of the total surface stress tensor. The method of calculation relies on the assumed condition that the sample is only slightly distorted from a spherical form. The equilibrium location of a levitated drop is combined with a formula introduced by Hasegawa (1979) to calcualte the ka dependence of the radiation force function. The present theory is valid for large as well as small ka values. Calculations in the small ka limit agree with previous theories and experimental results. Examples are presented for nonplane-wave modes as well as plane-wave rectangular modes.
A sixteenth-degree and sixteenth-order spherical harmonic lunar gravity field has been derived from the long-term Keplerian variations in the orbits of the Apollo subsatellites and Lunar Orbiter 5. This model resolves the major mascon gravity anomalies of the lunar near side and is in very good agreement with line-of-sight acceleration results. The far-side map shows the major ringed basins to be strong localized negative anomalies located in broad regions of positive gravity which correspond closely to the highlands. The rms pressure levels calculated from equivalent-surface height variations show that the moon and earth support nearly equal pressures, whereas Mars is appreciably stronger. The moon appears to support larger loads than earth owing to its weaker central gravity field and perhaps a colder upper lithosphere. Significant differences between the low-degree gravity and topography spectra indicate that the longer-wavelength topographic features are isostatically compensated.
The atomic hydrogen distribution at 250 km during December 1974 solstice was inferred, considering charge exchange equilibrium, from Atmosphere Explorer-C measurements of n(H(+)), n(O(+) and N(O). An empirical model, derived from the observations by least suqare analysis in terms of spherical harmonics, has the following characteristics: (1) n(H) increases by as much as a factor of two between the summer and winter hemispheres, (2) the n(H) diurnal variation is largest at the equator and (3) the diurnal variation is larger in the winter hemisphere than in the summer. Similar analysis of the gas temperature derived from n(N2) measurements reveals that all n(H) and Tg spherical harmonic coefficients are anticorrelated. Both the diurnal and latitudinal (annual) n(H) and Tg amplitudes are in substantial agreement with the zero flux condition, in which exospheric flow dominates the hydrogen distribution. The observed diurnal phase lag of n(H) with respect to Tg is about one hour, agreeing with theory.
Fluid motions in the electrically conducting core and magnetized crustal rocks are the two major sources of the magnetic field observed on or slightly above the Earth's surface. The exact separation of these two contributions is not possible without imposing a priori assumptions about the internal source distribution. Nonetheless models like these were developed for hundreds of years Gauss' method, least squares analysis with a truncated spherical harmonic expansion was the method of choice for more than 100 years although he did not address separation of core and crustal sources, but rather internal versus external ones. Using some arbitrary criterion for appropriate truncation level, we now extrapolate downward core field models through the (approximately) insulating mantle. Unfortunately our view can change dramatically depending on the degree of truncation for describing core sources.
The computation is described of a geopotential model to deg 360, a sea surface topography model to deg 10/15, and adjusted Geosat orbits for the first year of the exact repeat mission (ERM). This study started from the GEM-T2 potential coefficient model and it's error covariance matrix and Geosat orbits (for 22 ERMs) computed by Haines et al. using the GEM-T2 model. The first step followed the general procedures which use a radial orbit error theory originally developed by English. The Geosat data was processed to find corrections to the a priori geopotential model, corrections to a radial orbit error model for 76 Geosat arcs, and coefficients of a harmonic representation of the sea surface topography. The second stage of the analysis took place by doing a combination of the GEM-T2 coefficients with 30 deg gravity data derived from surface gravity data and anomalies obtained from altimeter data. The analysis has shown how a high degree spherical harmonic model can be determined combining the best aspects of two different analysis techniques. The error analysis was described that has led to the accuracy estimates for all the coefficients to deg 360. Significant work is needed to improve the modeling effort.
A full understanding of LISA (Laser Interferometer Space Antenna)'s science capability will require accurate models of incident waveform signals and the instrumental response. While Fisher matrix analysis is useful for some estimates, a full Bayesian treatment is needed for important cases at the limit of LISA's capability. We will apply fast analysis algorithms enabling accurate treatment with EOB (effective-one-body) waveforms and the full-featured LISA response to study the significance of higher spherical harmonics and mergers in LISA analysis.
We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.
Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.
Detailed analysis of the Automatic Rendezvous and Capture problem indicate a need for three different regions of mathematical description for the GN&C algorithms: (1) multi-vehicle orbital mechanics to the rendezvous interface point, i.e., within 100 nm; (2) relative motion solutions (such as Clohessy-Wiltshire type) from the far-field to the near-field interface, i.e., within 1 nm and; (3) close proximity motion - the near-field motion where the relative differences in the gravitational and orbit inertial accelerations can be neglected from the equations of motion. Limit boundaries to these regions can be precisely defined by further analysis and will be functions of the tracking measurement accuracies and the computer resources available for the solution of the algorithms. This paper analyzes the relative motion in Regions 2 and 3 above and present the derivation and discussion of the general case of non-spherical gravitational perturbed relative motion. Mathematical deviations from the numerically integrated spherical gravity case and solutions from the Clohessy-Wiltshire equations are presented in the analysis. Based upon this preliminary analysis, it is recommended that further efforts be used to assess the relative position and velocity differences in Region 2 due to non-spherical gravity harmonics and that viable GN&C algorithms be developed to include these gravity perturbations (especially the effects of the first gravity harmonic, J2).
The TOPEX/Poseidon satellite data can be used to improve the knowledge of the earth's gravitational field. The GPS data are especially useful for improving the gravity field over the world's oceans, where the current tracking data are sparse. Using realistic scenario for processing 10 days of GPS data, a covariance analysis is performed to obtain the expected improvement to the GEM-T2 gravity field. The large amount of GPS data and the large number of parameters (1979 parameters for the gravity field, plus carrier-phase biases, etc.) required special filtering techniques for efficient solution. The gravity-bin technique is used to compute the covariance matrix associated with the spherical harmonic gravity field. The covariance analysis shows that the GPS data from one 10-day arc of TOPEX/Poseidon with no a priori constraints can resolve medium degree and order (3-26) parameters with sigmas (standard deviations) that are an order of magnitude smaller than the corresponding sigmas of GEM-T2. When the information from GEM-T2 is combined with the TOPEX/Poseidon GPS measurements, an order-of-magnitude improvement is observed in low- and medium-degree terms with significant improvements spread over a wide range of degree and order.
An analysis of the magnetic field of Saturn is presented which takes into account both the Voyager 1 and 2 vector magnetic field observations. The analysis is based on the traditional spherical harmonic expansion of a scale potential to derive the magnetic field within 8 Saturn radii. A third-order zonal harmonic model fitted to Voyager 1 and 2 observations is found to be capable of predicting the magnetic field characteristics at one encounter based on those observed at another, unlike models including dipole and quadrupole terms only. The third-order model is noted to lead to significantly enhanced polar surface field intensities with respect to dipole models, and probably represents the axisymmetric part of a complex dynamo field.
The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.
The natural approach to a spectral analysis of data distributed on the sky employs spherical harmonic decomposition. A common problem encountered in practical astronomy is the lack of full sky coverage in the available data. For example, the removal of the Galactic plane data from the Cosmic Background Explorer (COBE) Differential Microwave Radiometer (DMR) sky maps compromises Fourier analysis of the cosmic microwave background (CMB) temperature distribution due to the loss of orthogonality of the spherical harmonics. An explicit method for constructing orthonormal functions on an incomplete (e.g., Galaxy-cut) sphere is presented. These functions should be used in the proper Fourier analysis of the COBE DMR sky maps to provide the correct input for the determination of the spectrum of primordial inhomogeneity. The results of such an analysis are presented in an accompanying Letter. A similar algebraic construction of appropriate functions can be devised for other astronomical applications.
The Zimm plot has been widely used to characterize the molecular dimensions of polymers from small-angle scattering experiments, where the reciprocal intensity is analyzed as a function of the square of the magnitude of the scattering wavevector Q. Here, this work explores the benefits of analyzing the reciprocal scattering intensity from deformed polymers, extending the original Zimm plot to anisotropic materials. In the small-angle limit, a tensorial extension of the Guinier law is found for the gyration tensor and the reciprocal single-chain structure factor. In the high-Q limit, application of the spherical harmonic expansion technique to the reciprocal structure factor permits direct model-independent analysis of spatially dependent molecular deformation of polymers. Additionally, the contributions from high-order spherical harmonics become insignificant in the reciprocal-intensity representation. The proposed generalized Zimm plot approach is demonstrated computationally with the affine deformation model and the Rouse model, and experimentally with small-angle neutron scattering measurements of deformed polystyrene melts.
New spherical harmonic models of the gravity and topography of Mars place important constraints on the structure and dynamics of the interior. The gravity and topography models are significantly phase coherent for harmonic degrees n less than 30 (wavelengths greater than 700 km). Loss of coherence below that wavelength is presumably due to inadequacies of the models, rather than a change in behavior of the planet. The gravity/topography admittance reveals two very different spectral domains: for n greater than 4, a simple Airy compensation model, with mean depth of 100 km, faithfully represents the observed pattern; for degrees 2 and 3, the effective compensation depths are 1400 and 550 km, respectively, strongly arguing for dynamic compensation at those wavelengths. The gravity model has been derived from a reanalysis of the tracking data for Mariner 9 and the Viking Orbiters, The topography model was derived by harmonic analysis of the USGS digital elevation model of Mars. Before comparing gravity and topography for internal structure inferences, we must ensure that both are consistently referenced to a hydrostatic datum. For the gravity, this involves removal of hydrostatic components of the even degree zonal coefficients. For the topography, it involves adding the degree 4 equipotential reference surface, to get spherically referenced values, and then subtracting the full degree 50 equipotential. Variance spectra and phase coherence of orthometric heights and gravity anomalies are addressed.
An improved theoretical model of Venusian global gravity has been obtained by fitting a tenth degree spherical harmonic series to 78 orbital arcs of Doppler tracking data from the Pioneer Venus Orbiter. Maps of the free-air anomaly and its formal error are presented. Isostatic anomaly and 'Geoid' maps are also presented, and their geophysical implications are discussed in details. Comparison with equivalent resolution topographic models reveals a strong correlation between long wavelength gravity and the topography of Venus. Analysis of the second degree harmonics showed two aspects of the orientation of the inertial axes of Venus: (1) a significant (about three degrees) departure of the axis of greatest inertia from the rotational axis; and (2) a near alignment of the axis of least inertia with the location of the subterrestrial point at the time of the next inferior conjunction with earth (December 16, 2101). A series of contour maps of the Venusian free-air anomalies is provided.