Contribution to the theory of the perturbed Kepler motion of a material point in the presence of a small tangential force
Equations of motion for geopotential determination of central gravitational field of earth influenced by tangential force
SEARCH · Search NASA
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.
Equations of motion for geopotential determination of central gravitational field of earth influenced by tangential force
The mathematical formulation required for predicting secular variation in the geopotential is developed for the case of a spherically symmetric, self-gravitating, viscoelastic earth model and an arbitrary surface load which can include a gravitational self-consistent ocean loading component. The theory is specifically applied to predict the present-day secular variation in the zonal harmonics of the geopotenial arising from the surface mass loading associated with the late Pleistocene glacial cycles. A procedure is outlined in which predictions of the present-day geopotential signal due to the late Pleistocene glacial cycles may be used to derive bounds on the net present-day mass flux from the Antarctic and Greenland ice sheets to the local oceans.
Two methods are used to investigate the vertical tilts of planetary waves as functions of zonal wavenumber and frequency. The vertical tilts are computed by cross-spectral analysis of the geopotential heights at different pressures. In the midlatitude troposphere, the eastward-moving waves had a westward tilt with height, as expected, but the westward-moving waves with frequencies higher than 0.2/d showed statistically significant eastward vertical tilts. For a free Rossby wave, this implies that the Eliassen-Palm flux is downward along with its energy propagation. A downward energy propagation suggests an upper-level source of these waves. It is proposed that the eastward-tilting waves were forced by the nonlinear interaction of stationary waves and baroclinically unstable cyclone-scale waves. The predicted vertical tilt and phase speed were consistent with the observations. In addition, simulations of a general circulation model were analyzed. In the control run, eastward-tilting waves disappeared when the sources of stationary waves were removed. This is consistent with the present theory.
Linear estimation theory, along with a new technique to compute relative data weights, was applied to the determination of the Earth's geopotential field and other geophysical model parameters using a combination of satellite ground-based tracking data, satellite altimetry data, and the surface gravimetry data. The relative data weights for the inhomogeneous data sets are estimated simultaneously with the gravity field and other geophysical and orbit parameters in a least squares approach to produce the University of Texas gravity field models. New techniques to perform calibration of the formal covariance matrix for the geopotential solution were developed to obtain a reliable gravity field error estimate. Different techniques, which include orbit residual analysis, surface gravity anomaly residual analysis, subset gravity solution comparisons and consider covariance analysis, were applied to investigate the reliability of the calibration.
A theory is presented that points out the existence of several long period and 'cross effects' in the coefficients in the expansion of the geopotential and in the motion of satellites. The tidal potential, defined as small periodic variations in the geopotential, was calculated. The influence of these geopotential variations on satellite perturbation is examined. Spherical harmonics were employed.
The rotation of Earth breaks time-reversal and reflection symmetries in an opposite sense north and south of the equator, leading to a topological origin for certain atmospheric and oceanic equatorial waves. Away from the equator, the rotating shallow-water and stably stratified primitive equations exhibit Poincaré inertia–gravity waves that have nontrivial topology as evidenced by their strict superinertial time scale and a phase singularity in frequency–wavevector space. This nontrivial topology then predicts, via the principle of bulk-interface correspondence, the existence of two equatorial waves along the equatorial interface, the Kelvin and Yanai waves. To directly test the nontrivial topology of Poincaré-gravity waves in observations, we examine ERA5 data and study cross correlations between the wind velocity and geopotential height of the midlatitude stratosphere at the 50 hPa height. We find the predicted vortex and antivortex in the relative phase of the geopotential height and velocity at the high frequencies of the waves. By contrast, lower-frequency planetary waves are found to have trivial topology also as expected from theory. These results demonstrate a new way to understand stratospheric waves and provide a new qualitative tool to investigate waves in other components of the climate system.
In recent years, observational evidence has been obtained for westward traveling planetary waves in the middle atmosphere with the aid of global data from satellites. There is no doubt that the fair portion of the observed traveling waves can be understood as the manifestation of the normal mode Rossby waves which are theoretically derived from the tidal theory. Some observational aspects of the structure and behavior of the normal model Rossby waves in the upper stratosphere are reported. The data used are the global stratospheric geopotential thickness and height analyses which are derived mainly from the Stratospheric Sounding Units (SSUs) on board TIROS-N and NOAA satellites. A clear example of the influence of the normal mode Rossby wave on the mean flow is reported. The mechanism considered is interference between the normal mode Rossby wave and the quasi-stationary wave.
Precision mean element (PME) satellite theories play a key role in orbit dynamics analyses. These theories employ: nonsingular orbital elements comprehensive force models Generalized Method of Averaging Numerical interpolation concepts The Draper Semianalytical Satellite Theory (DSST) (Refs. 1 - 6), whose development was led by the author, and the independently-developed Universal Semianalytical Method (USM) (Ref. 7) are examples of such theories. These theories provide the capability to tailor the force modeling to meet the desired computational speed vs. accuracy trade-off. The flexibility of such theories is demonstrated by their ability to include complicated atmosphere density models and spacecraft models in the perturbation theory context. The value of high speed satellite theories, in this era of computational plenty, is that they allow new ways of looking at astrodynamical problems such as orbit design (Refs. 8, 9) and atmosphere density updating (Refs. 10, 11). In the mid to late-1980 s, the geodynamics community led the development of very precise geopotential models such as GEM T2 and GEM T3 (Ref. 12), and with the subsequent analysis of the TOPEX flight data, JGM-2 and JGM-3 (Ref. 13). These were high degree and order geopotentials, at least 50 x 50. In 1993, the DSST implementation in the GTDS program was extended to include the 50 x 50 geopotential models (Ref. 14). The 50 x 50 geopotential, J2000 integration coordinate system, and solid Earth tide capabilities were integrated in GTDS by Scott Carter (Ref. 15). This capability demonstrated 1 m accuracy versus the TOPEX Precise Orbit Ephemerides. Subsequently the DSST Standalone program was also extended to include high degree and order geopotential models (Ref. 5). More recently GTDS has been hosted in the Linux PC environment. However, all of these efforts have been limited to modeling the motion of an artificial Earth satellite. They did not consider the additional complexities associated with lunar, planetary, or other natural satellite orbiters. Such complexities include: additional coordinate systems (associated with the direction of the north pole of rotation and the prime meridian of the new central bodies) (Ref. 16) normalized gravity model coefficients (desirable for high degree and order fields) (Ref. 17) indirect oblateness
Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.
The differential equations of the tidal oscillations of the earth were established under the assumption that the interior of the earth is laterally inhomogeneous. The theory was developed using vectorial and dyadic symbolism to shorten the exposition and to reduce the differential equations to a symmetric form convenient for programming and for numerical integration. The formation of tidal buldges on the surfaces of discontinuity and the changes in the internal density produce small periodic variations in the exterior geopotential which are reflected in the motion of artificial satellites. The analoques of Love elastic parameters in the expansion of exterior tidal potential reflect the asymmetric and inhomogeneous structure of the interior of the earth.
Upon looking at the various colossal interferometers, etc., discussed at this conference to test gravitational theory, one cannot avoid feeling that easier approaches exist. The use of low velocity, neutral atom matter waves in place of electromagnetic waves in sensitive inertial interferometer configurations is proposed. For applications, spacecraft experiments to sense a drag-free condition, to measure the Lense-Thirring precession, to measure the gravitomagnetic effect and/or the earth's geopotential (depending on altitude), and to detect long period gravitational waves are considered. Also, a terrestrial precision test of the equivalence principle on spin polarized atoms, capable of detecting effects of the 5th force is considered. While the ideas described herein are preliminary, the orders of magnitude are sufficiently tantalizing to warrant further study. Although existing proposed designs may be adequate for some of these experiments, the use of matter-wave interferometry offers reduced complexity and cost, and an absence of cryogenics.
A theory for the long-term rotational motion of the quasi-rigid Earth was constructed by numerical integration. The theory spans 72,000 years centered about 1968 A.D., and provides accurate rotational and positional data for the Earth in the recent past and the near future. The physical model is termed dynamically consistent because developments for the active forces and torques are truncated based solely on their magnitudes regardless of their origin. The model includes all appropriate forces and torques due to the geopotential and tidal effects as well as lunisolar and planetary contributions. The elastic and inelastic deformations due to tidal action were too small to affect the mass properties of the Earth at the truncation level of the model. However, long-term dissipative effects of the tidal forces and torques were not negligible. These considerations gave the model its quasi-rigid characterization. The numerical output provided both rotational and orbital-element data. The data were fitted throughout the 72,000-year range using Chebyshev polynomial series.
The equations of motion of an artificial satellite are given in nonsingular variables. Any term in the geopotential is considered as well as luni-solar perturbations up to an arbitrary power of r/r', r' being the geocentric distance of the disturbing body. Resonances with tesseral harmonics and with the moon or sun are also considered. By neglecting the shadow effect, the disturbing function for solar radiation is also developed in nonsingular variables for the long periodic perturbations. Formulas are developed for implementation of the theory in actual computations.
The equations of motion of an artificial satellite are given in nonsingular variables. Any term in the geopotential is considered as well as luni-solar perturbations up to an arbitrary power of r/r prime; r prime being the geocentric distance of the disturbing body. Resonances with tesseral harmonics and with the moon or sun are also considered. By neglecting the shadow effect, the disturbing function for solar radiation is also developed in nonsingular variables for the long periodic perturbations. Formulas are developed for implementation of the theory in actual computations.
The disturbing functions due to the geopotential and Luni-solar attractions are linear and bilinear forms in spherical harmonics. Making use of recurrence relations for the solid spherical harmonics and their derivatives, recurrence formulas are obtained for high degree terms as function of lower degree for any term of those disturbing functions and their derivative with respect to any element. The equations obtained are effective when a numerical integration of the equations of motion is appropriate. In analytical theories, they provide a fast way of obtaining high degree terms starting from initial very simple functions.
A number of observational, theoretical, and modeling studies have been conducted on the remote influence of the tropics on the midlatitudes. Teleconnection in the geopotential height field between the tropics and the extratropics associated with the Southern Oscillation in the so-called Pacific-North America (PNA) pattern, has often been cited as an example of such remote influence. Simmons et al. (1983) have shown that teleconnections such as PNA may arise both from barotropic instability of the extratropical zonal flow and from tropical forcing. The present investigation is only concerned with the dynamics of the part of the atmospheric teleconnection which arises from tropical forcing. Some fundamentals of Rossby wave theory in unsheared flow are discussed along with the linear response to equatorial forcing for horizontally, sheared zonal mean flow. A global shallow-water spectral model is employed in conjunction with a theoretical analysis to show the teleconnection patterns which develop under different large-scale mean circulation conditions.
A physical model of Mars is presented on the basis of light-scattering observations of the Martian atmosphere and surface and interior data obtained from observations of the geopotential field. A general description of the atmosphere is presented, with attention given to the circulation and the various cloud types, and data and questions on the blue haze-clearing effect and the seasonal darkening wave are summarized and the Mie scattering model developed to explain these observations is presented. The appearance of the planet from earth and spacecraft through Mariner 9 is considered, and attention is given to the preparation of topographical contour maps, the canal problem and large-scale lineaments observed from Mariner 9, the gravity field and shape of the planet and the application of Runcorn's geoid/convection theory to Mars. Finally, a summary of Viking results is presented and their application to the understanding of Martian geophysics is discussed.
Attention is given to the low-frequency variability of outgoing longwave radiation (OLR) fluctuations, their possible correlations over different parts of the globe, and their relationships with teleconnections obtained from other meteorological parameters, for example, geopotential and temperature fields. Simultaneous relationships with respect to the Southern Oscillation (Namais, 1978; Barnett, 1981) signal and the reference OLR fluctuation over the equatorial central Pacific are investigated. Emphasis is placed on the relative importance of the Southern Oscillation (SO) signal over preferred regions. Using lag cross-correlation statistics, possible lagged relationships between the tropics and midlatitudes and their relationships with the SO are then investigated. Only features that are consistent with present knowledge of the dynamics of the system are emphasized. Certain features which may not meet rigorous statistical significance tests but yet are either expected a priori from independent observations or are predicted from dynamical theories are also explored.