Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical harmonic method”

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

On determining the spectrum of primordial inhomogeneity from the COBE DMR sky maps: Method

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.

Gorski, Krzysztof M.↗
NASA NTRS1985

Separation of core and crustal magnetic field sources

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.

Shure, L.↗
NASA NTRS2001

Satellite Remote Sensing of Tropical Precipitation and Ice Clouds for GCM Verification

This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Evans, K. Franklin↗
NASA NTRS1984

Continental and oceanic crustal magnetization modelling

Inversion of magnetic data from the MAGSAT satellite, to arrive at intensities of magnetization of the Earth's crust, was performed by two different methods. The first method uses a spherical harmonic model of the magnetic field. The coefficients believed to represent sources in the Earth's crust can then be inverted to arrive at vertical dipole moments per unit area at the Earth's surface. The spherical harmonic models contain coefficients of degrees of harmonics up to 23. The dipole moment per unit area for a surface element can then be determined by summing the contribution for each individual degree of harmonic. The magnetic moments were calculated for continental and oceanic areas separately as well as over certain latitudinal segments. Of primary concern was to determine whether there are any differences between continental and oceanic areas. The second analysis with magnetization intensities was made using narrower ranges of degrees of harmonics, assuming that higher degrees are present in the core field signal.

Harrison, C. G. A.↗
NASA NTRS1998

Computation of the Effects of Inhomogeneous Clouds on Retrieval of Remotely Sensed Properties

Current and future earth observation programs depend on satellite measurements of radiance to retrieve the properties of clouds on a global basis. At present, this retrieval is made assuming that the clouds in the instrument field of view are plane parallel and independent of adjacent pixels. While this assumption is known to be false except in very limited cases, its impact can be evaluated, and if possible corrected, based on emerging theoretical techniques. In this study, the Spherical Harmonic Discrete Ordinate Method (SHDOM, Evans, 1996) has been used to assess the sensitivity of the retrieval to a variety of cloud parameters. SHDOM allows the plane parallel assumption to be relaxed and makes 2D and even 3D radiative solutions practical. A previous study (Chambers et al., 1996) assessed the effect of horizontal inhomogeneity in 45 LANDSAT scenes of boundary layer clouds over ocean. The four scenes studied here represent overcast, broken, scattered and strongly thermally forced cloud fields and are used to perform sensitivity studies to a wider variety of parameters. Comparisons are made at three solar zenith angles (theta (sub 0) = 0, 49, and 63 degrees) to avoid ambiguity in the results due to solar zenith angle.

Chambers, Lin H.↗
NASA NTRS1996

Independent Pixel and Two Dimensional Estimates of LANDSAT-Derived Cloud Field Albedo

A theoretical study has been conducted on the effects of cloud horizontal inhomogeneity on cloud albedo bias. A two-dimensional (2D) version of the Spherical Harmonic Discrete Ordinate Method (SHDOM) is used to estimate the albedo bias of the plane parallel (PP-IPA) and independent pixel (IPA-2D) approximations for a wide range of 2D cloud fields obtained from LANDSAT. They include single layer trade cumulus, open and closed cell broken stratocumulus, and solid stratocumulus boundary layer cloud fields over ocean. Findings are presented on a variety of averaging scales and are summarized as a function of cloud fraction, mean cloud optical depth, cloud aspect ratio, standard deviation of optical depth, and the gamma function parameter Y (a measure of the width of the optical depth distribution). Biases are found to be small for small cloud fraction or mean optical depth, where the cloud fields under study behave linearly. They are large (up to 0.20 for PP-IPA bias, -0.12 for IPA-2D bias) for large v. On a scene average basis PP-IPA bias can reach 0.30, while IPA-2D bias reaches its largest magnitude at -0.07. Biases due to horizontal transport (IPA-2D) are much smaller than PP-IPA biases but account for 20% RMS of the bias overall. Limitations of this work include the particular cloud field set used, assumptions of conservative scattering, constant cloud droplet size, no gas absorption or surface reflectance, and restriction to 2D radiative transport. The LANDSAT data used may also be affected by radiative smoothing.

Chambers, L. H.↗
NASA NTRS2024

Sensitivity studies of nighttime top-of-atmosphere radiances from artificial light sources using a 3-D radiative transfer model for nighttime aerosol retrievals

By accounting for surface-based light source emissions and top-of-atmosphere (TOA) downward lunar fluxes, we adapted the spherical harmonics discrete ordinate method (SHDOM) 3-dimensional (3-D) radiative transfer model (RTM) to simulate nighttime 3-D TOA radiances as observed from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day/Night Band (DNB) on board the Suomi-NPP satellite platform. Used previously for daytime 3-D applications, these new SHDOM enhancements allow for the study of the impacts of various observing conditions and aerosol properties on simulated VIIRS-DNB TOA radiances. Observations over Dakar, Senegal, selected for its bright city lights and a large range of aerosol optical depth (AOD), were investigated for potential applications and opportunities for using observed radiances containing VIIRS-DNB “bright pixels” from artificial light sources to conduct aerosol retrievals. We found that using the standard deviation (SD) of such bright pixels provided a more stable quantity for nighttime AOD retrievals than direct retrievals from TOA radiances. Further, both the mean TOA radiance and SD of TOA radiances over artificial sources are significantly impacted by satellite viewing angles. Light domes, the enhanced radiances adjacent to artificial light sources, are strong functions of aerosol properties and especially aerosol vertical distribution, which may be further utilized for retrieving aerosol layer height in future studies. Through inter-comparison with both day- and nighttime Aerosol Robotic Network (AERONET) data, the feasibility of retrieving nighttime AODs using 3-D RTM SHDOM over artificial light sources was demonstrated. Our study shows strong potential for using artificial light sources for nighttime AOD retrievals, while also highlighting larger uncertainties in quantifying surface light source emissions. This study underscores the need for surface light emission source characterizations as a key boundary condition, which is a complex task that requires enhanced input data and further research. We demonstrate how quality-controlled nighttime light data from NASA’s Black Marble product suite could serve as a primary input into estimations of surface light source emissions for nighttime aerosol retrievals.

Jianglong Zhang↗
NASA NTRS2021

GRACE Fast Mascons from Spherical Harmonics and a Regularization Design Trade Study

Mass concentration (mascon) solutions have become a prominent medium for investigating time-variable gravity recovered by GRACE and GRACE Follow-On. While GRACE Level-2 spherical harmonic products require various post-processing techniques to eliminate correlated noise, mascon formulations employ spatial regularization strategies during the estimation step to improve signal recovery. However, mascon estimation has traditionally required large computing resources and GRACE Level-1B processing capabilities. In this study, we show that a typical mascon estimation system can be reformulated to allow for the estimation of regularized mascons from Level-2 spherical harmonics. Provided that spherical harmonic solution covariances are available, the computed mascons will be mathematically equivalent to similar mascons estimated from Level-1B observations. This method is computationally efficient, better leverages GRACE spherical harmonics than past methods, and matches the performance of typical mascon solutions without locking scientists into pre-determined regularization designs. We develop a proof-of-concept solution using ITSG-Grace2018 and compare results with traditional mascons from JPL and NASA GSFC. We then assess the effects of spherical harmonic truncation and use of regularization correlations on basin signal recovery. We find that spherical expansions to degree and order 60 provide the minimum expansion necessary to study most basins, while larger expansions help further localize signals. We also find that diagonal regularizations (i.e., regularizations that do not contain inter-mascon correlations) are adversely affected by leakage, especially across boundaries such as coastlines where signals are not highly correlated, whereas including inter-mascon correlations and regional boundaries in the regularization greatly improves signal recovery.

M. J. Croteau↗
NASA NTRS1971

Lunar gravity derived from long-period satellite motion, a proposed method

A method was devised to determine the spherical harmonic coefficients of the lunar gravity field. The method consists of a two-step data reduction and estimation process. Pseudo-Doppler data were generated simulating two different lunar orbits. The analysis included the perturbing effects of the L1 lunar gravity field, the earth, the sun, and solar radiation pressure. Orbit determinations were performed on these data and long-period orbital elements were obtained. The Kepler element rates from these solutions were used to recover L1 lunar gravity coefficients. Overall results of the experiment show that lunar gravity coefficients can be accurately determined and that the method is dynamically consistent with long-period perturbation theory.

Ferrari, A. J.↗
NASA NTRS1971

An Empirical Method for Determining the Lunar Gravity Field

A method has been devised to determine the spherical harmonic coefficients of the lunar gravity field. This method consists of a two-step data reduction and estimation process. In the first step, a weighted least-squares empirical orbit determination scheme is applied to Doppler tracking data from lunar orbits to estimate long-period Kepler elements and rates. Each of the Kepler elements is represented by an independent function of time. The long-period perturbing effects of the earth, sun, and solar radiation are explicitly modeled in this scheme. Kepler element variations estimated by this empirical processor are ascribed to the non-central lunar gravitation features. Doppler data are reduced in this manner for as many orbits as are available. In the second step, the Kepler element rates are used as input to a second least-squares processor that estimates lunar gravity coefficients using the long-period Lagrange perturbation equations.

Ferrari, A. J.↗
NASA NTRS1973

Lunar gravity derived from long-period satellite motion - A proposed method.

A new method has been devised to determine the spherical harmonic coefficients of the lunar gravity field. This method consists of a two-step data reduction and estimation process. In the first step, a weighted least-squares empirical orbit determination scheme is applied to Doppler tracking data from lunar orbits to estimate long-period Kepler elements and rates. Each of the Kepler elements is represented by an independent function of time. The long-period perturbing effects of the earth, sun, and solar radiation are explicitly modeled in this scheme. Kepler element variations estimated by this empirical processor are then ascribed to the non-central lunar gravitation features. Doppler data are reduced in this manner for as many orbits as are available. In the second step, the Kepler element rates are used as input to a second least-squares processor that estimates lunar gravity coefficients using the long-period Lagrange perturbation equations.

Ferrari, A. J.↗
NASA NTRS1985

Geophysical interpretation of the magnetic anomalies of the Earth derived from MAGSAT data

The ambiguities about geophysical implications based on the correlation of scalar magnetic anomalies and geological features were investigated. A method was developed to convert scalar magnetic anomalies into a map of the lateral variations of magnetic susceptibility of the lithosphere. This map is directly correlated with the causative sources. The method is based on spherical harmonic analysis of lateral variations seen on the scalar magnetic anomaly map and those of the lithospheric magnetic susceptibility. The harmonic coefficients are related through the fundamental causality relationship governing a magnetized body and its associated scalar magnetic anomaly. The main features of the resulting magnetic susceptibility anomalies are outlined.

Arkani-Hamed, J.↗
NASA NTRS1991

A simple method for correcting spatially resolved solar intensity oscillation observations for variations in scattered light

A measurement of the intensity distribution in an image of the solar disk will be corrupted by a spatial redistribution of the light that is caused by the earth's atmosphere and the observing instrument. A simple correction method is introduced here that is applicable for solar p-mode intensity observations obtained over a period of time in which there is a significant change in the scattering component of the point spread function. The method circumvents the problems incurred with an accurate determination of the spatial point spread function and its subsequent deconvolution from the observations. The method only corrects the spherical harmonic coefficients that represent the spatial frequencies present in the image and does not correct the image itself.

Jefferies, S. M.↗
NASA NTRS2000

Geoid Recovery Using Geophysical Inverse Theory Applied to Satellite to Satellite Tracking Data

This report describes a new method for determination of the geopotential, or the equivalent geoid. It is based on Satellite-to-Satellite Tracking (SST) of two co-orbiting low earth satellites separated by a few hundred kilometers. The analysis is aimed at the GRACE Mission, though it is generally applicable to any SST data. It is proposed that the SST be viewed as a mapping mission. That is, the result will be maps of the geoid or gravity, as contrasted with determination of spherical harmonics or Fourier coefficients. A method has been developed, based on Geophysical Inverse Theory (GIT), that can provide maps at a prescribed (desired) resolution and the corresponding error map from the SST data. This computation can be done area by area avoiding simultaneous recovery of all the geopotential information. The necessary elements of potential theory, celestial mechanics, and Geophysical Inverse Theory are described, a computation architecture is described, and the results of several simulations presented. Centimeter accuracy geoids with 50 to 100 km resolution can be recovered with a 30 to 60 day mission.

Gaposchkin, E. M.↗
NASA NTRS2006

Terrestrial Water Mass Load Changes from Gravity Recovery and Climate Experiment (GRACE)

Recent studies show that data from the Gravity Recovery and Climate Experiment (GRACE) is promising for basin- to global-scale water cycle research. This study provides varied assessments of errors associated with GRACE water storage estimates. Thirteen monthly GRACE gravity solutions from August 2002 to December 2004 are examined, along with synthesized GRACE gravity fields for the same period that incorporate simulated errors. The synthetic GRACE fields are calculated using numerical climate models and GRACE internal error estimates. We consider the influence of measurement noise, spatial leakage error, and atmospheric and ocean dealiasing (AOD) model error as the major contributors to the error budget. Leakage error arises from the limited range of GRACE spherical harmonics not corrupted by noise. AOD model error is due to imperfect correction for atmosphere and ocean mass redistribution applied during GRACE processing. Four methods of forming water storage estimates from GRACE spherical harmonics (four different basin filters) are applied to both GRACE and synthetic data. Two basin filters use Gaussian smoothing, and the other two are dynamic basin filters which use knowledge of geographical locations where water storage variations are expected. Global maps of measurement noise, leakage error, and AOD model errors are estimated for each basin filter. Dynamic basin filters yield the smallest errors and highest signal-to-noise ratio. Within 12 selected basins, GRACE and synthetic data show similar amplitudes of water storage change. Using 53 river basins, covering most of Earth's land surface excluding Antarctica and Greenland, we document how error changes with basin size, latitude, and shape. Leakage error is most affected by basin size and latitude, and AOD model error is most dependent on basin latitude.

Seo, K.-W.↗
NASA NTRS1981

A stochastic-dynamic model for global atmospheric mass-field statistics

Global atmospheric mass field error correlations based on satellite observations and on numerical forecasts show strong and systematic latitude dependence. A model for the latitude dependent spatial correlation structure of mass field forecast errors is derived from dynamical considerations. Three methods of solution were tested. In the first method, the equation was solved by expansion in spherical harmonics, and the correlation function was computed analytically using the expansion coefficients. In the second method, the finite difference equivalent of the equation was solved using a fast poisson solver. The correlation function was computed using stratified sampling of the individual realizations. In the third method, a higher order equation was derived, and solved directly in finite differences by two successive applications of the fast poisson solver. The three methods were compared for accuracy and efficiency, and the third method was chosen as clearly superior.

Ghil, M.↗
NASA NTRS1992

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗