On the nature of “overshoot” in a collisionless shock

In this paper, by comparing experimental data on bow shock with MHD-relationships on a flat shock discontinuity, allowing for the presence behind the front of turbulent electrostatic oscillations and of an ion beam, an analysis is made of the nature of the “overshoot” of magnetic field (density) behind the front of a collisionless shock wave. It is shown that the large value of plasma compression in the overshoot region ($\text{n}{}_2{}^{\mathrm{ov}}\text{n}{}_1\text{) ～ 6}$, in excess of the maximum allowable value of density jump ($\text{n}{}_2\text{n}{}_1\text{)|}{}_{\mathrm{<mtext>max</mtext>}}\text{= (γ +}\text{1}\text{γ ? 1}\text{)|}{}_{\mathrm{\gamma }}\text{=}{}_{\mathrm{<mtext>5</mtext><mtext>3</mtext>}}\text{= 4}\text{at a Mach number}\text{M → ∞}$, is attributable to the presence in the “overshoot” of a high level of lowerhybrid electrostatic oscillations with an energy density W ? nT.     

Love numbers of the giant planets

We calculate the Love numbers k_n for n = 2 to 10, and determine the “gravitational noise” from tides. The new values k₂ for Jupiter, Saturn, and Uranus yield new estimates for the planetary dissipation functions: $\text{Q}{}_{\mathrm{<mtext>J</mtext>}}\text{? 2.5 × 10}{}^4$, $\text{Q}{}_{\mathrm{<mtext>S</mtext>}}\text{? 1.4 × 10}{}^4\text{, Q}{}_{\mathrm{<mtext>U</mtext>}}\text{? 5 × 10}{}^3$.     

Cyclotron side-band emissions from ring-current electrons

VLF-emissions with subharmonic cyclotron frequency from magnetospheric electrons have been detected by the S³-A satellite (Explorer 45) whose orbit is close to the magnetic equatorial plane where the wave-particle interaction is most efficient. These emissions are observed during the main phase of a geomagnetic storm in the nightside of the magnetosphere outside of the plasmasphere around L = 3–5. The emissions consist essentially of two frequency regimes, one below the equatorial electron gyro-frequency, , and the other above . The emissions below are whistler mode and there is a sharp band of “missing emissions” along . The emissions above are electrostatic mode and the frequency ranges up to . It is concluded that these emissions are generated by the enhanced relativity low energy (1–5 keV) ring current electrons, penetrating into the nightside magnetosphere during the main phase of a magneto storm. Although the high energy (50–350 keV) electrons showed remarkable changes of pitch angle distribution, their associations with VLF-emissions are not so significant as those of low energy electrons.     

Changes in the concentration of mesospheric ozone during the total solar eclipse

Measurements of dayglow radiance of $\text{O}{}_2\text{(}{}^1\text{Δ}{}_{\mathrm{g}}\text{)}$ and OH(7,2) bands are reported. Ground based photometers were used to monitor zenith radiance of 1270 and 694 nm emissions during the total solar eclipse of 16 February 1980. Altitude distribution of 1270 nm intensity was derived from ground based observations. A set of altitude distributions of $\text{O}{}_2\text{(}{}^1\text{Δ}{}_{\mathrm{g}}\text{)}$ were thus obtained throughout the eclipse. These altitude distributions were converted into ozone distributions using the rate equations for formation and loss of ozone and $\text{O}{}_2\text{(}{}^1\text{Δ}{}_{\mathrm{g}}\text{)}$ molecules. Results indicate an increase in the ozone concentration at mid-eclipse. OH(7,2) emission did not show enhancement during totality. This may mean that there was no increase in OH concentration during the eclipse.     

Geopotential coefficients of order 29, from analysis of the orbit of satellite 1967-104B

The orbit of the satellite 1967-104B has been analysed as it passed through 29:2 resonance with the Earth's gravitational field between January 1977 and September 1978. From the changes in inclination and eccentricity the following lumped 29th-order geopotential harmonic coefficients were obtained: $\text{10}{}^9\text{C}\text{?}{}_{29}{}^{0.2}\text{= 4.1 ± 0.8}$, $\text{10}{}^9\text{S}\text{?}{}_{29}{}^{0.2}\text{= 10.3 ± 2.4}$, $\text{10}{}^9\text{C}\text{?}{}_{29}{}^{1.1}\text{= ? 160 ± 19}$, $\text{10}{}^9\text{S}\text{?}{}_{29}{}^{1.1}\text{= 79 ± 10}$, $\text{10}{}^9\text{C}\text{?}{}_{29}{}^{\mathrm{?1.3}}\text{= 38 ± 14}$, $\text{10}{}^9\text{S}\text{?}{}_{29}{}^{\mathrm{?1.3}}\text{= 19 ± 5}$. These values have been compared with existing comprehensive geopotential models: the best agreement is with the model of Rapp (1981).     

Analysis of HI content in galaxies

We have collected data on 241 galaxies from 13 sources and made a statistical analysis after reduction to a uniform system. We found that the Hubble sequence is one of increasing $\text{〈}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{M}{}_{\mathrm{<mtext>T</mtext>}}\text{〉}$ and $\text{〈}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{L}{}_{\mathrm{<mtext>B</mtext>}}\text{〉}$, these mean values increasing monotonically from .0016 and .024 at E to .084 and .83 at Im, but the dispersion is large.The HI content in barred spiral is greater than that in ordinary spirals, and this is consistent with their statistics of angular momentum and colour.The HI content is related to colour; it is greater in bluer systems. The large dispersion suggests that it also depends on some other factors, but these are smoothed out when averaged over each type, resulting in a linear relation between $\text{〈log(}\text{M}{}_{\mathrm{<mtext>M</mtext>}}\text{M}{}_{\mathrm{<mtext>T</mtext>}}\text{〉}$ and 〈(B ? V^O_T)〉. Unlike the colour-colour diagram, the large dispersion on the $\text{log}\text{(}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{L}{}_{\mathrm{<mtext>B</mtext>}}\text{) ? (B ? V}{}^0{}_{\mathrm{<mtext>T</mtext>}}\text{)}$ is not related to peculiar galaxies.     

The determination and analysis of the orbit of NIMBUS 1 ROCKET, 1964-52B: Part 2: Variations in orbital eccentricity

The influence of aerodynamic drag and the geopotential on the motion of the satellite 1964-52B is considered. A model of the atmosphere is adopted that allows for oblateness, and in which the density behaviour approximates to the observed diurnal variation. A differential equation governing the variation of the eccentricity, e, combining the effects of air drag with those of the Earth's gravitational field is given. This is solved numerically using as initial conditions 310 computed orbits of 1964-52B.The observed values of eccentricity are modified by the removal of perturbations due to luni-solar attraction, solid Earth and ocean tides, solar radiation pressure and low-order long-periodic tesseral harmonic perturbations. The method of removal of these effects is given in some detail. The behaviour of the orbital eccentricity predicted by the numerical solution is compared with the modified observed eccentricity to obtain values of atmospheric parameters at heights between 310 and 430 km. The daytime maximum of air density is found to be at 14.5 hours local time. Analysis of the eccentricity near 15th order resonance with the geopotential yielded values of four lumped geopotential harmonics of order 15, namely: $\text{10}{}^9\text{C}{}^{\mathrm{1,0}}{}_{15}\text{= ?78.8 ± 7.0}$, $\text{10}{}^9\text{S}{}^{\mathrm{1,0}}{}_{15}\text{= ?69.4 ± 5.3}$, $\text{10}{}^9\text{C}{}^{\mathrm{?1,2}}{}_{15}\text{= ?41.6 ± 3.5}$$\text{10}{}^9\text{S}{}^{\mathrm{?1,2}}{}_{15}\text{= ?26.1 ± 8.9}$, at inclination 98.68°.     

A possible mechanism of formation of faculae

We propose a new heating mechanism of faculae. We think that the formation of faculae is a result of the Joule dissipation of the Hall current generated by the interaction of the convection field of granules in an active region and the inter-granular magnetic field. For a region to generate effectively Hall current, its characteristic length must be such that the magnetic Reynolds number is less than 1. The equation of energy balance in the facula region is

$\text{16σT}{}^3{}_{\mathrm{p}}\text{(T}{}_{\mathrm{l}}\text{? T}{}_{\mathrm{p}}\text{)}{}^{\mathrm{n}}\text{HP}{}_{\mathrm{s}}{}^{\mathrm{a}}\text{H}{}^{\mathrm{?}}\text{=}\text{Qn}{}_{\mathrm{s}}\text{m}{}_{\mathrm{i}}\text{u}{}_{\mathrm{x}}{}^2\text{(τ}{}^2{}_{\mathrm{in}}\text{ω}{}_{\mathrm{i}}\text{)}$

.For five observational models of faculae, we calculated the corresponding velocity fields, and the results are in basic agreement with the observed fields. The present mechanism explains the dependence of the facula brightness on the magnetic and velocity fields, the apparent distribution of the faculae on the solar disk and suggest a possible interpretation of the five structures of faculae.     

Energy level perturbation of hydrogen in the Robertson-Walker metric

We calculated the perturbations in the energy level of the hydrogen atom in the four states $\text{1S}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{, 2S}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{2P}{}_{\mathrm{<mtext>3</mtext><mtext>2</mtext>}}\text{and}\text{2S}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, caused by the background curvature of the Robertson-Walker metric. We found that it is only when the radius of curvature is less than 10^?7 cm that an energy level split of the size of the Lamb shift can occur between $\text{2S}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{and}\text{2P}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$.     

An estimator of the underlying size distribution of overlapping impact-craters

The stochastic model of lunar type impact-crater formation which assumes (a) random impacts, (b) circular craters, each obliterating any portions of earlier craters lying within, and (c) a probability P_i(t) that a newly formed crater (primary or secondary) has an area a_i is analyzed to develop a method of estimating P_i from the final overlapping pattern. It is found that if each crater is weighted by the fraction of the rim which is visible and which lies in an observation area A, then the expected value of the weighted sum $\text{Ω}{}_{\mathrm{i}}$ of craters of area a_i is simply proportional to P_i for any degree of coverage under several conditions, including (a) constant P_i for all i, and (b) P_i stepping from a constant early value to zero (for some i's) with otherwise arbitrary bombardment. Furthermore, in the general case, the expected value of the contribution $\text{ΔΩ}{}_{\mathrm{i}}\text{(t}{}_0\text{)}\text{to}\text{Ω}{}_{\mathrm{i}}$ produced during t₀ ± Δt/2 is found to be proportional to P_i(t₀). Thus measurement of $\text{Ω}{}_{\mathrm{i}}$ in the first two cases, or of $\text{ΔΩ}{}_{\mathrm{i}}$ if crater age data is available in the last case, provides an estimate of the desired P_i. Therefore the $\text{Ω}{}_{\mathrm{i}}$ introduce the correct weighting factors that just compensate for the effect of overlap.Expressions for the variances of $\text{Ω}{}_{\mathrm{i}}\text{and}\text{Ω = Σ}{}_{\mathrm{i}}\text{Ω}{}_{\mathrm{i}}$ are derived from which it is shown that under the above conditions, $\text{Ω}{}_{\mathrm{i}}\text{/Ω}\text{or}\text{ΔΩ}{}_{\mathrm{i}}\text{/ΔΩ}$ are consistent estimators of P_i. Formal evaluation of the variances is carried out in the special case of constant P_i and no secondary cratering. A criterion for the degree of coverage is given; in particular it is shown that the expectation of $\text{σ = Σ}{}_{\mathrm{i}}\text{a}{}_{\mathrm{i}}\text{Ω}{}_{\mathrm{i}}$ at saturation is just A.     

Deuteron whistler and trans-equatorial propagation of the ion cyclotron whistler

New ion cyclotron whistlers which have the asymptotic frequency of one half the local proton gyrofrequency, $\text{G}{}_{\mathrm{p}}\text{2}$, and the minimum (or equatorial) proton gyrofrequency, G_pm, along the geomagnetic field line passing through the satellite have been found in the low-latitude topside ionosphere from the spectrum analysis of ISIS VLF electric field data received at Kashima, Japan. Ion cyclotron whistlers with asymptotic frequency of G_pm or $\text{G}{}_{\mathrm{pm}}\text{2}$ are observed only in the region of $\text{B}{}_{\mathrm{m}}\text{>}\text{B}\text{2}$ or rarely $\text{B}{}_{\mathrm{m}}\text{>}\text{B}\text{4}$, where B is the local magnetic field and B_m is the mini magnetic field along the geomagnetic field line passing through the satellite.The particles with one half the proton gyrofrequency may be the deuteron or alpha particle. Theoretical spectrograms of the electron whistlers (R-mode) and the ion cyclotron whistlers (L-mode) propagating along the geomagnetic field lines are computed for the appropriate distributions of the electron density and the ionic composition, and compared with the observed spectrograms.The result shows that the ion cyclotron whistler with the asymptotic frequency of $\text{G}{}_{\mathrm{p}}\text{2}$ is the deuteron whistler, and that the ion cyclotron whistlers with the asymptotic frequency of G_pm or $\text{G}{}_{\mathrm{pm}}\text{2}$ are caused by the trans-equatorial propagation of the proton or deuteron whistler from the other hemisphere.     

Latitudinal structures of discrete arcs resulting from viscous interaction between sheared plasma flows

Latitudinal structures of discrete arcs are modelled as a consequence of the quasi-steady magnetosphere-ionosphere coupling involving viscous interaction between sunward and anti-sunward plasma flows in the magnetosphere. The quasi-steady state in the magnetosphere and ionosphere coupling is described by the magnetospheric and ionospheric current conservation and the field-aligned currentpotential relation assuming adiabatic electron motion along field lines. The upward and downward fieldaligned currents are assumed to be stably maintained by vorticity-induced space charges in the region of plasma flow reversal, where divergence of the magnetospheric electric field E_⊥ is negative and positive, respectively. By introducing the effective conductance Σ_dc arising from the anomalous viscosity, a specific relation between the dc field-aligned current density J_∥ and the magnetospheric electric field E_⊥ is derived as $\text{J}{}_{\mathrm{\parallel }}\text{=−Σ}{}_{\mathrm{dc}}\text{div}\text{E}{}_{\mathrm{\perp }}$. Sufficiently large potential drops to accelerate auroral electrons are shown to exist along the auroral field lines originating from the flow reversal region with div $\text{E}{}_{\mathrm{\perp }}\text{< 0}$. It is shown that the latitudinal structure of a discrete arc is primarily determined by the magnetospheric potential structure and the characteristic width is on the order of 10 km at the ionospheric altitude.     

Titan aerosols: Optical properties and vertical distribution

An analysis of Titan's solar phase variation as a function of wavelength together with the continuum geometric albedo makes it possible to set limits on the real part of the refractive index and on the average particle size of the aerosol component of Titan's atmosphere: $\text{1.5 ?n}{}_{\mathrm{r}}\text{< 2.0}\text{and}\text{0.20 μm <}\text{r}\text{?0.35 μm}$. If n_ris known $\text{r}$ can be determined to within a few percent, and varies inversely with n_r. Using this information in a two-layer model of a methane-aerosol atmosphere and comparing the result with Titan's visible and near-infrared methane spectrum leads to the conclusion that the top layer of Titan's atmosphere contains 0.01 km atm of methane and 2.5 extinction optical depths of aerosol, while the data are consistent with a bottom layer containing 2.2 km atm of methane and about 7.5 aerosol optical depths for $\text{n}{}_{\mathrm{<mtext>r</mtext>}}\text{= 1.7,}\text{r}\text{= 0.25 μ}\text{m}$.     

New high-speed photometry of EH Lib and a binary interpretation of its period change

Six times of maxima of the ultrashort-period cepheid variable EH Librae were measured in 1980 May to June and in 1981 January, with a three-channel photocounting high-speed photoelectric photometer. These, together with all the photoelectric times of maxima over the past 30 years, are used to re-examine the nature of the change of the period. We found that we can fix the times of maxima by the following formula

$\text{T}{}_{\max }\text{= T}{}_0\text{+P}{}_0\text{E+}\text{1}\text{2}\text{βE}{}^2\text{+A}\text{sin}\text{2π}\text{EP}{}_0\text{E}{}_0$

where T₀ = HJD 2433438.6088 and P₀ = 0.0884132445 d represent the initial maximum epoch and the pulsation period, β = ?2.8 × 10^?8/yr; A = 0.0015 d, P₀ = 6251 d = 17.1 yr are the semi-amplitude and the period of the sine curve, and E is the number of periods elapsed since T₀, and (E₀ = 70700).If we interpret this 17.1 year periodicity as a modulation of the phase of maximum by binary motion, then the semi-amplitude of the orbital radial velocity variation is K = 2πasini/E₀ = 0.45 km/s and the mass function is

$\text{f(m)=}\text{m}{}^3{}_2\text{sin}{}^3\text{i}\text{(m}{}_1\text{m}{}_2\text{)}{}^2\text{=}\text{(a}\text{sin}\text{i)}{}^3\text{E}{}^2{}_0\text{=6 x 10}{}^{\mathrm{?5}}\text{M}{}_{\mathrm{\odot }}$

     

The influence of ethane and acetylene upon the thermal structure of the Jovian atmosphere

Ethane and acetylene, both of which possess more efficient emission bands than methane, have been incorporated into a thermal structure model for the atmosphere of Jupiter. Choosing for illustrative purposes the mixing ratios $\text{[}\text{C}{}_2\text{H}{}_6\text{]}\text{[}\text{H}{}_2\text{]}\text{= 10}{}^{\mathrm{?5}}$ and $\text{[}\text{C}{}_2\text{H}{}_2\text{]}\text{[}\text{H}{}_2\text{]}\text{= 5 × 10}{}^{\mathrm{?7}}$, it is found that these hydrocarbon gases lower the atmospheric temperature within the thermal inversion region by as much as 20 K, subsequently reducing the emission intensity of the 7.7 μm CH₄ band below the observed result. It is qualitatively shown, however, that this cooling by C₂H₆ and C₂H₂ could be compensated by aerosol heating resulting from a uniformily mixed aerosol which absorbs 15% of the incident solar radiation. Such aerosol heating has been suggested by uv albedo observations.     

Electron runaway in turbulent astrophysical plasmas

An astrophysical electron acceleration process is described which involves turbulent plasma effects: the acceleration mechanism will operate in ‘collision free’ magnetoactive astrophysical plasmas when ion-acoustic turbulence is generated by an electric field which acts parallel to the ambient magnetic lines of force. The role of ‘anomalous’ (ion-sound) resistivity is crucial in maintaining the parallel electric field. It is shown that, in spite of the turbulence, a small fraction of the electron population can accelerate freely, i.e. runaway, in the high parallel electric potential. The number density n^(B) of the runaway electron component is of order $\text{n}{}^{\mathrm{(B)}}\text{?}\text{n}\text{2}\text{(}\text{c}{}_{\mathrm{s}}\text{U}{}_{\mathrm{\parallel }}{}^{\mathrm{?}}\text{)}{}^2$, where n = background electron number density, c_s = ion-sound speed and U_∥^? = relative drift velocity between the electron and ion populations. The runaway mechanism and the number density n^(B) do not depend critically on the details of the non-linear saturation of the ion-sound instability.     

On fourier-analysis of geomagnetic pulsations pc-1, “two-fold” and “three-fold” pulsations pc-1 and fundamental frequencies of the magnetospheric cavity—I

The paper gives the results of detailed studies of the frequency spectra $\text{S}{}_{\mathrm{s}}\text{(?)}$ of the chain of the wave packets F_s(t) of geomagnetic pulsations PC-1 recorded at the Novolazarevskaya station. The bulk of the energy of F_s(t) is concentrated in the vicinity of the central frequencies $\text{?}{}_{\mathrm{s0}}$ of spectra—the carrier frequencies of the signals. The velocity $\text{V}{}_0\text{≌ 6.10}{}^3\text{km s}{}^{\mathrm{?1}}$ of the flux of protons generating these signals correspond to them. The spectra of the signals have oscillations—“satellites” irregularly distributed in frequency. These satellites, as the authors believe, testify to the presence of the individual groups of protons of low concentration whose velocities vary within 10³–10⁴ km s^?1.Their energy is only of the order of 10^?2–10^?3 of the energy of the main proton flux. Clearly pronounced maxima on double and triple frequencies $\text{? = 2?}{}_{\mathrm{s0}}\text{and}\text{3?}{}_{\mathrm{s0}}$ are detected. They show that the generation of pulsations PC-1 is accompanied by the generation on the overtones of wave packets called in this paper “two-fold” and “three-fold” pulsations PC-1. Intensive symmetrical satellites of a modulation character have been discovered on frequencies $\text{?}{}^{\mathrm{\pm }}{}_{\mathrm{sK}}$. Frequency differences $\text{Δ?}{}_{\mathrm{sK}}{}^{\mathrm{\pm }}\text{=}\text{¦?}{}_{\mathrm{s0}}\text{? ?}{}_{\mathrm{sK}}{}^{\mathrm{\pm }}\text{¦}\text{= (0.011,0.022}\text{and}\text{0.035)}\text{Hz}$ correspond to them. The authors believe that the values of Δ?^±_sK are resonance frequencies of the magnetospheric cavity in which geomagnetic pulsations PC-1 are generated. It is established that the values of Δ?^±_sK coincide closely with the carrier frequencies of geomagnetic pulsations PC-3 and PC-4 generated in the magnetosphere. This leads to the conclusion that the resonance oscillations of the magnetospheric cavity are their source. Thus, the generation of geomagnetic pulsations of different types and resonance oscillations in the magnetosphere are integrated into a unified process. The importance of the results obtained and the necessity to check further their trustworthiness and universality, using experimental data gathered in different conditions, is stressed.     

The locating of hydromagnetic whistler sources and determination of their generating proton spectra

Results are given of the calculations of the group delay time propagating τ(ω, φ₀) of hydromagnetic whistlers, using outer ionospheric models closely resembling actual conditions. The τ(ω, φ₀) dependencies were compared with the experimental data of τ_exp(ω, φ₀) obtained from sonagrams. The sonagrams were recorded in the frequency range $\text{? ? (0.5?2.5) Hz}$ at observation points located at geomagnetic latitudes φ₀ = (53?66)° and in the vicinity of the geomagnetic poles. This investigation has led us to new and important conclusions.The wave packets (W.P.) forming hydromagnetic whistlers (H.W.) are mainly generated in the plasma regions at L = 3.5?4.0. This is not consistent with ideas already expressed in the literature that their generation region is $\text{L ? 3?10}$. The overwhelming majority of the τ_exp values differ considerably from the times at which wave packets would, in theory, propagate along the magnetic field lines corresponding to those of the geomagnetic latitudes φ₀ of the observation points. The second important fact is that the W.P. frequency ω is less than $\text{Ω}{}_{\mathrm{H}}$ everywhere along its propagation trajectory, including the apogee of the magnetic force line ($\text{Ω}{}_{\mathrm{H}}$ is the proton gyrofrequency). Proton flux spectra $\text{E ? (30?120)}$ keV, responsible for H.W. generation, were determined. Comparison of the Explorer-45 and OGO-3 measurements published in the literature, with our data, showed that the proton flux density energy responsible for the H.W. excitation $\text{N}{}_{\mathrm{p}}\text{(}\text{MV}{}_6{}^2\text{2}\text{) ? (5 × 10}{}^{\mathrm{?3}}\text{?10}{}^{\mathrm{?1}}\text{)}\text{H}{}_{\mathrm{a}}{}^2\text{8π}$ where H_a is the magnetic field force in the generation region of these W.P. The electron concentration is $\text{N}{}_{\mathrm{a}}\text{? (10}{}^2\text{?10}{}^3\text{) cm}{}^{\mathrm{?3}}$. The values given in the literature are $\text{N}{}_{\mathrm{a}}\text{? (10}{}^{\mathrm{?10}}\text{?10}{}^3\text{) cm}{}^{\mathrm{?3}}$. The e data considered also leads to the conclusion that the generating mechanism of the W.P. studied probably always co-exists with the mechanism of their amplification.     

15th order resonance terms using the decaying orbit of TETR-3

The orbit of TETR-3 (1971-83B), inclination: 33°, passed through resonance with 15th order geopotential terms in February 1972. The resonance caused the orbit inclination to increase by 0.015°. Analysis of 48 sets of mean Kepler elements for this satellite in 1971–1972 (across the resonance) has established the following strong constraint for high degree, 15th order gravitational terms (normalized):

$\text{10}{}^9\text{(}\text{C, S}\text{)}{}_{15}\text{= (28.3 ± 3.0, 7.4 ± 3.0) =}\text{0.001(}\text{C, S}\text{)}{}_{\mathrm{15,15}}\text{?0.015(}\text{C, S}\text{)}{}_{\mathrm{17,15}}\text{+0.073(}\text{C, S}\text{)}{}_{\mathrm{19,15}}\text{?0.219(}\text{C, S}\text{)}{}_{\mathrm{21,15}}\text{+0.477(}\text{C, S}\text{)}{}_{\mathrm{23,15}}\text{?0.781(}\text{C, S}\text{)}{}_{\mathrm{25,15}}\text{+1.000(}\text{C, S}\text{)}{}_{\mathrm{27,15}}\text{?0.0963(}\text{C, S}\text{)}{}_{\mathrm{29,15}}\text{+0.622(}\text{C, S}\text{)}{}_{\mathrm{31,15}}\text{?0.119(}\text{C, S}\text{)}{}_{\mathrm{33,15}}\text{?0.290(}\text{C, S}\text{)}{}_{\mathrm{35,15}}\text{+0.403(}\text{C, S}\text{)}{}_{\mathrm{37,15}}\text{?0.223(}\text{C, S}\text{)}{}_{\mathrm{39,15}}\text{?0.058(}\text{C, S}\text{)}{}_{\mathrm{41,15}}\text{+…}$

This result combined with previous results on high inclination 15th order and other resonant orbits suggests that the coefficients of the gravity field beyond the 15th degree are smaller than Kaula's rule ($\text{10}{}^{\mathrm{?5}}\text{l}{}^2$).