J2 Invariant Relative Orbits for Spacecraft Formations

An analytic method is presented to establish J ₂ invariant relative orbits. Working with mean orbit elements, the secular drift of the longitude of the ascending node and the sum of the argument of perigee and mean anomaly are set equal between two neighboring orbits. By having both orbits drift at equal angular rates on the average, they will not separate over time due to the J₂ influence. Two first order conditions are established between the differences in momenta elements (semi-major axis, eccentricity and inclination angle) that guarantee that the drift rates of two neighboring orbits are equal on the average. Differences in the longitude of the ascending node, argument of perigee and initial mean anomaly can be set at will, as long as they are setup in mean element space. For near polar orbits, enforcing both momenta element constraints may result in impractically large relative orbits. It this case it is shown that dropping the equal ascending node rate requirement still avoids considerable relative orbit drift and provides substantial fuel savings.This revised version was published online in October 2005 with corrections to the Cover Date.     

Analytical short-term orbit predictions withJ 3 andJ 4 in terms of KS elements

Analytical theory for short-term orbit motion of satellite orbits with Earth's zonal harmonicsJ ₃ andJ ₄ is developed in terms of KS elements. Due to symmetry in KS element equations, only two of the nine equations are integrated analytically. The series expansions include terms of third power in the eccentricity. Numerical studies with two test cases reveal that orbital elements obtained from the analytical expressions match quite well with numerically integrated values during a revolution. Typically for an orbit with perigee height, eccentricity and inclination of 421.9 km, 0.17524 and 30 degrees, respectively, maximum differences of 27 and 25 cm in semimajor axis computation are noted withJ ₃ andJ ₄ term during a revolution. For application purposes, the analytical solutions can be used for accurate onboard computation of state vector in navigation and guidance packages.     

Quasi-critical orbits for artificial lunar satellites

We study the problem of critical inclination orbits for artificial lunar satellites, when in the lunar potential we include, besides the Keplerian term, the J ₂ and C ₂₂ terms and lunar rotation. We show that, at the fixed points of the 1-D averaged Hamiltonian, the inclination and the argument of pericenter do not remain both constant at the same time, as is the case when only the J ₂ term is taken into account. Instead, there exist quasi-critical solutions, for which the argument of pericenter librates around a constant value. These solutions are represented by smooth curves in phase space, which determine the dependence of the quasi-critical inclination on the initial nodal phase. The amplitude of libration of both argument of pericenter and inclination would be quite large for a non-rotating Moon, but is reduced to <0°.1 for both quantities, when a uniform rotation of the Moon is taken into account. The values of J ₂, C ₂₂ and the rotation rate strongly affect the quasi-critical inclination and the libration amplitude of the argument of pericenter. Examples for other celestial bodies are given, showing the dependence of the results on J ₂, C ₂₂ and rotation rate.     

Refinement of the geopotential scale factorR 0 on

Altimetric measurements of the GEOSAT satellite were used for the determination of geopotential scale factorR ₀. The geopotential valueW ₀ on the geoid surface was then computed (W ₀ =GM/R ₀).The GEOSAT Geophysical Data Records (GDR's) covering an initial period of the Exact Repeat Mission (ERM) were filtered and processed. The necessary corrections were made in order to allow a precise detection of the sea surface. Gravitational geopotential, rotation and permanent tides were taken into account and the equipotential surface which is the best approximation of the sea surface was found.The determination of the potential valueW ₀ on the mean geoid surface in this way is very promising. An associated value withW ₀ - the geopotential scale factorR ₀ - seems to be a very good Earth dimension defining quantity. Moreover, there are many possible applications ofW ₀ (R ₀) in modern geophysics.The incorporation of one of these parameters - we now recommendR ₀ - into the set of the Primary Geodetic Parameters (PGP) is discussed and suggested.     

A procedure of selection of meteors from major streams for determination of mean orbits

A procedure of selection of meteoroids from major streams is suggested and applied to the IAU Lund photographic database modified by a check for internal consistency among orbital elements (3411 orbits). Limits for choice of stream members were defined by break points on the plots of the cumulative numberN _C vs. the Southworth-HawkinsD discriminant. For the break points were considered the points from which the dependenceN _C vs.D changes to a quasi-linear one, and with the increasingD, N _C changes only moderately. Except for the Taurids which desire a separate analysis, theN _C vs.D diagrams are presented for the following major meteoroid streams: Quadrantids, Lyrids, Aquarids, Capricornids, N and S Aquarids, Perseids, Orionids, Leonids and Geminids. The mean orbits, velocities and radiants of the streams are derived and compared with the osculating orbits of their parent bodies. The limitingD _B was found to be a function of the number of the stream membersN _CB. Omitting the exceptionally concentrated Geminids, the relation is in the formD _B = 0.058 *ln(N _CB) – 0.04.     

The Non-Linear Stability of L4 in the Restricted Three-Body Problem when the Bigger Primary is a Triaxial Rigid Body

The non-linear stability of L ₄ in the restricted three-body problem has been studied when the bigger primary is a triaxial rigid body with its equatorial plane coincident with the plane of motion. It is found that L ₄ is stable in the range of linear stability except for three mass ratios:

Motion close to the Hopf Bifurcation Of the vertical family of periodic orbits of L4

The paper deals with different kinds of invariant motions (periodic orbits, 2D and 3D invariant tori and invariant manifolds of periodic orbits) in order to analyze the Hamiltonian direct Hopf bifurcation that takes place close to the Lyapunov vertical family of periodic orbits of the triangular equilibrium point L₄ in the 3D restricted three-body problem (RTBP) for the mass parameter, μ greater than (and close to) μ_R (Routh’s mass parameter). Consequences of such bifurcation, concerning the confinement of the motion close to the hyperbolic orbits and the 3D nearby tori are also described.     

Solutions and periodicity of satellite relative motion under even zonal harmonics perturbations

Finding relative satellite orbits that guarantee long-term bounded relative motion is important for cluster flight, wherein a group of satellites remain within bounded distances while applying very few formationkeeping maneuvers. However, most existing astrodynamical approaches utilize mean orbital elements for detecting bounded relative orbits, and therefore cannot guarantee long-term boundedness under realistic gravitational models. The main purpose of the present paper is to develop analytical methods for designing long-term bounded relative orbits under high-order gravitational perturbations. The key underlying observation is that in the presence of arbitrarily high-order even zonal harmonics perturbations, the dynamics are superintegrable for equatorial orbits. When only J ₂ is considered, the current paper offers a closed-form solution for the relative motion in the equatorial plane using elliptic integrals. Moreover, necessary and sufficient periodicity conditions for the relative motion are determined. The proposed methodology for the J ₂-perturbed relative motion is then extended to non-equatorial orbits and to the case of any high-order even zonal harmonics (J _2n , n ≥ 1). Numerical simulations show how the suggested methodology can be implemented for designing bounded relative quasiperiodic orbits in the presence of the complete zonal part of the gravitational potential.     

Infinite Feigenbaum sequences and spirals in the vicinity of the Lagrangian periodic solutions

We studied systematically cases of the families of non-symmetric periodic orbits in the planar restricted three-body problem. We took interesting information about the evolution, stability and termination of bifurcating families of various multiplicities. We found that the main families of simple non-symmetric periodic orbits present a similar dynamical structure and bifurcation pattern. As the Jacobi constant changes each branch of the characteristic of a main family spirals around a focal point-terminating point in x- at which the Jacobi constant is C _∞ = 3 and their periodic orbits terminate at the corotation (at the Lagrangian point L₄ or L₅). As the family approaches asymptotically its termination point infinite changes of stability to instability and vice versa occur along its characteristic. Thus, infinite bifurcation points appear and each one of them produces infinite inverse Feigenbaum sequences. That is, every bifurcating family of a Feigenbaum sequence produces the same phenomenon and so on. Therefore, infinite spiral characteristics appear and each one of them generates infinite new inner spirals and so on. Each member of these infinite sets of the spirals reproduces a basic bifurcation pattern. Therefore, we have in general large unstable regions that generate large chaotic regions near the corotation points L₄, L₅, which are unstable. As C varies along the spiral characteristic of every bifurcating family, which approaches its focal point, infinite loops, one inside the other, surrounding the unstable triangular points L₄ or L₅ are formed on their orbits. So, each terminating point corresponds to an asymptotic non-symmetric periodic orbit that spirals into the corotation points L₄, L₅ with infinite period. This is a new mechanism that produces very large degree of stochasticity. These conclusions help us to comprehend better the motions around the points L₄ and L₅ of Lagrange.     

Bifurcation points and intersections of families of periodic orbits in the three-dimensional restricted three-body problem

Intersections of families of three-dimensional periodic orbits which define bifurcation points are studied. The existence conditions for bifurcation points are discussed and an algorithm for the numerical continuation of such points is developed. Two sequences of bifurcation points are given concerning the family of periodic orbits which starts and terminates at the triangular equilibrium pointsL ₄,L ₅. On these sequences two trifurcation points are identified forµ = 0.124214 andµ = 0.399335. The caseµ = 0.5 is studied in particular and it is found that the space families originating at the equilibrium pointsL ₂,L ₃,L ₄,L ₅ terminate on the same planar orbitm _1v of the familym.     

Frozen orbits for satellites close to an Earth-like planet

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.     

Bifurcations of plane with three-dimensional periodic orbits in the elliptic restricted problem

We show that the procedure employed in the circular restricted problem, of tracing families of three-dimensional periodic orbits from vertical self-resonant orbits belonging to plane families, can also be applied in the elliptic problem. A method of determining series of vertical bifurcation orbits in the planar elliptic restricted problem is described, and one such series consisting of vertical-critical orbits (a _v=+1) is given for the entire range (0,1/2) of the mass parameter . The initial segments of the families of three-dimensional orbits which bifurcate from two of the orbits belonging to this series are also given.     

Astrophysical properties of red giants in three open clusters older than the hyades

Results fromCMT ₁T₂ 1T2 broad-band and DDO intermediate-band photometry are presented for G and K giants in the old open clusters NGC. 2482, NGC 3680, and IC 4651. Two independent photometric criteria have been used to separate red field stars from the physical members of the clusters. Recent calibrations of the DDO andCMT ₁T₂ systems have been used to derive reddening, distance moduli, metallicities, effective temperatures, and surface gravities. Rough estimates of masses have also been made. The giants of NGC 2482 and IC 4651 have CN strengths nearly identical to the Hyades giants, while those of NGC 3680 are slightly richer in CN than the nearby K giants.CMT1T2 abundance analysis in NGC 2482 and NGC 3680 yield [Fe/H]_MT = - 0.1 ± 0.1 as derived from the iron lines, while abundances derived from the CNO - contaminated (C - M) index are 0.4 dex higher. BothCMT ₁T₂ and DDO data support the conclusion that 1C 4651, with [Fe/H] = + 0.2 ± 0.1, is on the metalrich side of the distribution of intermediate and old open clusters. Finally, the mass results suggest that the clump stars in NGC 3680 and. IC 4651 could have undergone mass loss before reaching their helium core burning phase of evolution. Supported in part by the Consejo Nacional de Investigaciones Cientificas y Tecnicas (CONICET) of Argentina. Visiting astronomer of Cerro Tololo Inter-American Observatory supported by the National Science Foundation under contract No. AST 74-04128.     

Tolerance of allogromiid Foraminifera to severely elevated carbon dioxide concentrations: Implications to future ecosystem functioning and paleoceanographic interpretations

Increases in the partial pressure of carbon dioxide (pCO₂) in the atmosphere will significantly affect a wide variety of terrestrial fauna and flora. Because of tight atmospheric–oceanic coupling, shallow-water marine species are also expected to be affected by increases in atmospheric carbon dioxide concentrations. One proposed way to slow increases in atmospheric pCO₂ is to sequester CO₂ in the deep sea. Thus, over the next few centuries marine species will be exposed to changing seawater chemistry caused by ocean–atmospheric exchange and/or deep-ocean sequestration. This initial case study on one allogromiid foraminiferal species (Allogromia laticollaris) was conducted to begin to ascertain the effect of elevated pCO₂ on benthic Foraminifera, which are a major meiofaunal constituent of shallow- and deep-water marine communities. Cultures of this thecate foraminiferan protist were used for 10–14-day experiments. Experimental treatments were executed in an incubator that controlled CO₂ (15 000; 30 000; 60 000; 90 000; 200 000 ppm), temperature and humidity; atmospheric controls (i.e., ~ 375 ppm CO₂) were executed simultaneously. Although the experimental elevated pCO₂ values are far above foreseeable surface water pCO₂, they were selected to represent the spectrum of conditions expected for the benthos if deep-sea CO₂ sequestration becomes a reality. Survival was assessed in two independent ways: pseudopodial presence/absence and measurement of adenosine triphosphate (ATP), which is an indicator of cellular energy. Substantial proportions of A. laticollaris populations survived 200 000 ppm CO₂ although the mean of the median [ATP] of survivors was statistically lower for this treatment than for that of atmospheric control specimens. After individuals that had been incubated in 200 000 ppm CO₂ for 12 days were transferred to atmospheric conditions for ~ 24 h, the [ATP] of live specimens (survivors) approximated those of the comparable atmospheric control treatment. Incubation in 200 000 ppm CO₂ also resulted in reproduction by some individuals. Results suggest that certain Foraminifera are able to tolerate deep-sea CO₂ sequestration and perhaps thrive as a result of elevated pCO₂ that is predicted for the next few centuries, in a high-pCO₂ world. Thus, allogromiid foraminiferal “blooms” may result from climate change. Furthermore, because allogromiids consume a variety of prey, it is likely that they will be major players in ecosystem dynamics of future coastal sedimentary environments.     

On Sitnikov-like motions generating new kinds of 3D periodic orbits in the R3BP with prolate primaries

The existence of new equilibrium points is established in the restricted three-body problem with equal prolate primaries. These are located on the Z-axis above and below the inner Eulerian equilibrium point L ₁ and give rise to a new type of straight-line periodic oscillations, different from the well known Sitnikov motions. Using the stability properties of these oscillations, bifurcation points are found at which new types of families of 3D periodic orbits branch out of the Z-axis consisting of orbits located entirely above or below the orbital plane of the primaries. Several of the bifurcating families are continued numerically and typical member orbits are illustrated.     

Rosette Central Configurations, Degenerate Central Configurations and Bifurcations

In this paper we find a class of new degenerate central configurations and bifurcations in the Newtonian n-body problem. In particular we analyze the Rosette central configurations, namely a coplanar configuration where n particles of mass m₁ lie at the vertices of a regular n-gon, n particles of mass m₂ lie at the vertices of another n-gon concentric with the first, but rotated of an angle π /n, and an additional particle of mass m₀ lies at the center of mass of the system. This system admits two mass parameters μ = m₀/m₁ and ε = m₂/m₁. We show that, as μ varies, if n > 3, there is a degenerate central configuration and a bifurcation for every ε > 0, while if n = 3 there is a bifurcation only for some values of ε.     

On separation of major meteoroid streams from the sporadic background

The paper deals with a search for chosen photographic meteoroid streams compiled from the IAU Meteor Data Center Lund catalogue from which less than 2% of the orbits had to be removed due to internal inconsistency among the orbit parameters. Additional 35 orbits were removed due to extremely high hyperbolic velocities. The final set consists of 3411 orbits. Members of the Quadrantids, Lyrids, Perseids and Geminids were searched for, firstly, by a stream-search procedure utilizing the Southworth-HawkinsD-criterion. This choice, as a rule, represents the most abundant filament of the stream. Secondly, rate distribution histograms ofD were divided into region of shower meteors and region of sporadic background meteors. The searched database with a relatively low abundance of sporadic meteors in the analyzed periods simplified this choice, and followingly, fitting the obtained values by means of power and exponential functions, the limitingD _s for particular showers were derived. The derivedD _s appears as the optimum value, as for higherD, the number of sporadic meteors included in the stream sample increases more rapidly than the number of additional shower meteors, and for smallerD, the number of shower meteors decreases quicker than the number of omitted sporadic meteors. The following counts of shower meteorsN and limitingD _s were found: Quadrantids (39, 0.22), Lyrids (11, 0.15), Perseids (595, 0.53) and Geminids (224, 0.32). Efficiency of the procedure was tested comparing the number of sporadic meteors in the region of radiant area and the neighbouring regions of the same size.     

The Phase Space Structure Around L4 in the Restricted Three-Body Problem

The phase space structure around L ₄ in the restricted three-body problem is investigated. The connection between the long period family emanating from L ₄ and the very complex structure of the stability region is shown by using the method of Poincarés surface of section. The zero initial velocity stability region around L ₄ is determined by using a method based on the calculation of finite-time Lyapunov characteristic numbers. It is shown that the boundary of the stability region in the configuration space is formed by orbits suffering slow chaotic diffusion.     

Nonlinear dynamical analysis for displaced orbits above a planet

Nonlinear dynamical analysis and the control problem for a displaced orbit above a planet are discussed. It is indicated that there are two equilibria for the system, one hyperbolic (saddle) and one elliptic (center), except for the degenerate h _z ^max, a saddle-node bifurcation point. Motions near the equilibria for the nonresonance case are investigated by means of the Birkhoff normal form and dynamical system techniques. The Kolmogorov–Arnold–Moser (KAM) torus filled with quasiperiodic trajectories is measured in the τ ₁ and τ ₂ directions, and a rough algorithm for calculating τ ₁ and τ ₂ is proposed. A general iterative algorithm to generate periodic Lyapunov orbits is also presented. Transitions in the neck region are demonstrated, respectively, in the nonresonance, resonance, and degradation cases. One of the important contributions of the paper is to derive necessary and sufficiency conditions for stability of the motion near the equilibria. Another contribution is to demonstrate numerically that the critical KAM torus of nontransition is filled with the (1,1)-homoclinic orbits of the Lyapunov orbit.     

Integration theory for the restricted problem of three axisymmetric bodies

In this paper the circular planar restricted problem of three axisymmetric ellipsoids S _i (i = 1, 2, 3), such that their equatorial planes coincide with the orbital plane of the three centres of masses, be considered. The equations of motion of infinitesimal body S ₃ be obtained in the polar coordinates. Using iteration approach we have given an approximation for another integral, which independent of the Jacobian integral, in the case of P-type orbits (near circular orbits surrounding both primaries).