Engineering PapersSearch

Engineering topics

Rubincam, David Parry

Publications and source records attributed to Rubincam, David Parry.

At least 19 records

Possible Lack of Low-Mass Meteoroids in the Earth's Meteoroid Flux Due to Space Erosion?

The Earth's cumulative meteoroid flux, as found by Halliday et al. (1996), may have a shallower slope for meteoroid masses in the range 0.1-2.5 kg compared to those with masses greater than 2.5 kg when plotted on a log flux vs. log mass graph. This would indicate a lack of low-mass objects. While others such as Ceplecha (1992) find no shallow slope, there may be a reason for a lack of 0.1-2.5 kg meteoroids which supports Halliday et al.'s finding. Simple models show that a few centimeters of space erosion in stony meteoroids can reproduce the bend in Halliday et al.'s curve at approximately 2.5 kg and give the shallower slope.

Rubincam, David Parry

Goulds Belt, Interstellar Clouds, and the Eocene Oligocene Helium-3 Enhancement

Drag from hydrogen in the interstellar cloud which formed Gould's Belt may have sent interplanetary dust particle (IDPs) and small meteoroids with embedded helium to the Earth, perhaps explaining part the helium-3 flux increase seen in the sedimentary record near the Eocene-Oligocene transition. Assuming the Solar System passed through part of the cloud, IDPs in the inner Solar System may have been dragged to Earth, while dust and small meteoroids in the asteroid belt up to centimeter size may have been dragged to the resonances, where their orbital eccentricities were pumped up into Earth-crossing orbits; however, this hypotheses does not explain the Popigai and Chesapeake Bay impacts.

Rubincam, David Parry

Goulds Belt, Interstellar Clouds, and the Eocene-Oligocene Helium-3 Spike

Drag from hydrogen in the interstellar cloud which formed Gould's Belt may have sent small meteoroids with embedded helium to the Earth, perhaps explaining part or all of the (sup 3) He spike seen in the sedimentary record at the Eocene-Oligocene transition. Assuming the Solar System passed through part of the cloud, meteoroids in the asteroid belt up to centimeter size may have been dragged to the resonances, where their orbital eccentricities were pumped up into Earth-crossing orbits.

Gould's Belt

Cosmic Ray Exposure Ages of Stony Meteorites: Space Erosion or Yarkovsky?

Space erosion from dust impacts may set upper limits on the cosmic ray exposure (CRE) ages of stony meteorites. A meteoroid orbiting within the asteroid belt is bombarded by both cosmic rays and interplanetary dust particles. Galactic cosmic rays penetrate only the first few meters of the meteoroid; deeper regions are shielded. The dust particle impacts create tiny craters on the meteoroid's surface, wearing it away by space erosion (abrasion) at a particular rate. Hence a particular point inside a meteoroid accumulates cosmic ray products only until that point wears away, limiting CRE ages. The results would apply to other regolith-free surfaces in the solar system as well, so that abrasion may set upper CRE age limits which depend on the dusty environment. Calculations based on N. Divine's dust populations and on micrometeoroid cratering indicate that stony meteoroids in circular ecliptic orbits at 2 AU will record 21Ne CRE ages of approx.176 x 10(exp 6) years if dust masses are in the range 10(exp -21) - 10(exp -3) kg. This is in broad agreement with the maximum observed CRE ages of approx. 100 x 10(exp 6) years for stones. High erosion rates in the inner solar system may limit the CRE ages of Near-Earth Asteroids (NEAs) to approx. 120 x 10(exp 6) years. If abrasion should prove to be approx. 6 times quicker than found here, then space erosion may be responsible for many of the measured CRE ages of main belt stony meteorites. In that case the CRE ages may not measure the drift time to the resonances due to the Yarkovsky effects as in the standard scenario, and that for some reason Yarkovsky is ineffective.

Yarkovsky

Shape Shifting Satellites in Binary Near-Earth Asteroids: Do Meteoroid Impacts Play a Role in BYORP Orbital Evolution?

Less than catastrophic meteoroid impacts over 10(exp 5) years may change the shape of small rubble-pile satellites in binary NEAs, lengthening the average BYORP (binary Yarkovsky-Radzievskii-Paddack) rate of orbital evolution. An estimate of shape-shifting meteoroid fluxes give numbers close enough to causing random walks in the semimajor axis of binary systems to warrant further investigation

Rubincam, David Parry

The Precession Index, A Nonlinear Energy Balance Model, And Seversmith Psychroterms

An important component of Milankovitch's astronomical theory of climate change is the precession index. The precession index, along with the Earth's tilt and orbital eccentricity, are believed to be the major controlling factors of climate change in the last few million years. The precession index is e sin omega(sub s) where e is the Earth's orbital eccentricity and omega(sub s) measures how close the Sun is to the Earth at midsummer. When omega(sub s) = 90deg the Sun is close to the Earth during northern summer, and at 270deg it is far from the Earth during northern summer. The precession index varies with time, because both the eccentricity e and the parameter omega(sub s) are constantly changing due to disturbances in the Earth's orbit by other planets, and due to the precession of the Earth, The change is largely periodic, with a period of about 23,000 years.

Rubincam, David Parry

Saturn's Rings, the Yarkovsky Effects, and the Ring of Fire

The dimensions of Saturn's A and B rings may be determined by the seasonal Yarkovsky effect and the Yarkovsky-Schach effect; the two effects confine the rings between approximately 1.68 and approximately 2.23 Saturn radii, in reasonable agreement with the observed values of 1.525 and 2.267. The C ring may be sparsely populated because its particles are transients on their way to Saturn; the infall may create a luminous Ring of Fire around Saturn's equator. The ring system may be young: in the past heat flow from Saturn's interior much above its present value would not permit rings to exist.

Rubincam, David Parry

Is Asteroid 951 Gaspra in a Resonant State with Its Spin Increasing Due to YORP?

Asteroid 951 Gaspra appears to be in an obliquity resonance with its spin increasing due to the YORP effect. Gaspra, an asteroid 5.8 km in radius, is a prograde rotator with a rotation period of 7.03 hours. A three million year integration indicates its orbit is stable over at least this time span. From its known shape and spin axis orientation and assuming a uniform density, Gaspra's axial precession period turns out to be nearly commensurate with its orbital precession period, which leads to a resonance condition with consequent huge variations in its obliquity. At the same time its shape is such that the Yarkovsky-O'Keefe-Radzievskii-Paddack effect (YORP effect for short) is increasing its spin rate. The YORP cycle normally leads from spin-up to spin-down and then repeating the cycle; however, it appears possible that resonance trapping can at least temporarily interrupt the YORP cycle, causing spin-up until the resonance is exited. This behavior may partially explain why there is an excess of fast rotators among small asteroids. YORP may also be a reason for small asteroids entering resonances in the first place.

Rubincam, David Parry

Polar Wander on Triton and Pluto Due to Volatile Migration

Polar wander may occur on Triton and Pluto because of volatile migration. Triton, with its low obliquity, can theoretically sublimate volatiles (mostly nitrogen) at the rate of approximately 10(exp 14) kilograms per year from the equatorial regions and deposit them at the poles. Assuming Triton to be rigid on the sublimation timescale, after approximately 10(exp 5) years the polar caps would become large enough to cancel the rotational flattening, with a total mass equivalent to a global layer approximately 120-250 m in depth. At this point the pole wanders about the tidal bulge axis, which is the line joining Triton and Neptune. Rotation about the bulge axis might be expected to disturb the leading side/trailing side cratering statistics. Because no such disturbance is observed, it may be that Triton's mantle viscosity is too high but its surface volatile inventory is too low to permit wander. On the other hand, its mantle viscosity might be low, so that any uncompensated cap load might be expected to wander toward the tidal bulge axis. In this case, the axis of wander passes through the equator from the leading side to the trailing side; rotation about this wander axis would not disturb the cratering statistics. Low-viscosity polar wander may explain the bright southern hemisphere: this is the pole which is wandering toward the equator. In any case the permanent polar caps may be geologically very young. Polar wander may possibly take place on Pluto, due to its obliquity oscillations and perihelion-pole geometry. However, Pluto is probably not experiencing any wander at present. The Sun has been shining strongly on the poles over the last half of the obliquity cycle, so that volatiles should migrate to the equator, stabilizing the planet against wander. Spacecraft missions to Triton and Pluto which measure the dynamical flattening could give information about the accumulation of volatiles at the poles. Such information is best obtained by measuring gravity and topography from orbiters, as was done for Mars with the highly successful Mars Global Surveyor.

Rubincam, David Parry

Gravitational Core-Mantle Coupling and the Acceleration of the Earth

Gravitational core-mantle coupling may be the cause of the observed variable acceleration of the Earth's rotation on the 1000 year timescale. The idea is that density inhomogeneities which randomly come and go in the liquid outer core gravitationally attract density inhomogeneities in the mantle and crust, torquing the mantle and changing its rotation state. The corresponding torque by the mantle on the core may also explain the westward drift of the magnetic field of 0.2 deg per year. Gravitational core-mantle coupling would stochastically affect the rate of change of the Earth's obliquity by just a few per cent. Its contribution to polar wander would only be about 0.5% the presently observed rate. Tidal friction is slowing down the rotation of the Earth, overwhelming a smaller positive acceleration from postglacial rebound. Coupling between the liquid outer core of the Earth and the mantle has long been a suspected reason for changes in the length-of-day. The present investigation focuses on the gravitational coupling between the density anomalies in the convecting liquid outer core and those in the mantle and crust as a possible cause for the observed nonsecular acceleration on the millenial timescale. The basic idea is as follows. There are density inhomogeneities caused by blobs circulating in the outer core like the blobs in a lava lamp; thus the outer core's gravitational field is not featureless. Moreover, these blobs will form and dissipate somewhat randomly. Thus there will be a time variability to the fields. These density inhomogeneities will gravitationally attract the density anomalies in the mantle.

Rubincam, David Parry

A Geometrical Approach to Bell's Theorem

Bell's theorem can be proved through simple geometrical reasoning, without the need for the Psi function, probability distributions, or calculus. The proof is based on N. David Mermin's explication of the Einstein-Podolsky-Rosen-Bohm experiment, which involves Stern-Gerlach detectors which flash red or green lights when detecting spin-up or spin-down. The statistics of local hidden variable theories for this experiment can be arranged in colored strips from which simple inequalities can be deduced. These inequalities lead to a demonstration of Bell's theorem. Moreover, all local hidden variable theories can be graphed in such a way as to enclose their statistics in a pyramid, with the quantum-mechanical result lying a finite distance beneath the base of the pyramid.

Rubincam, David Parry

Insolation and the Precession Index

Simple nonlinear climate models yield a precession index-like term in the temperature. Despite its importance in the geologic record, the precession index e sin omega, where e is the Earth's orbital eccentricity and omega is the Sun's perigee in the geocentric frame, is not present in the insolation at the top of the atmosphere. Hence there is no one-for-one mapping of 23,000 and 19,000 year periodicities from the insolation to the paleoclimate record; a nonlinear climate model is needed to produce these periods. Two such models, a grey body and an energy balance climate model with an added quadratic term, produce e sin omega terms in temperature. These terms, which without feedback mechanisms achieve extreme values of about plus or minus 0.48 K for the grey body and plus or minus 0.64 K for the energy balance model, simultaneously cool one hemisphere while they warm the other. Moreover, they produce long-term cooling in the northern hemisphere when the Sun's perigee is near northern solstice and long-term warming in the northern hemisphere when the perigee is near southern solstice. Thus this seemingly paradoxical mechanism works against the standard model which requires cool northern summers (Sun far from Earth in northern summer) to build up northern ice sheets, so that if the standard model is correct it may be more efficient than previously thought. Alternatively, the new mechanism could possibly be dominant and indicate southern hemisphere control of the northern ice sheets, wherein the southern oceans undergo a long-term cooling when the Sun is close to the Earth during southern summer. The cold water eventually flows north, cooling the northern hemisphere. This might explain why the northern oceans lag the southern ones when it comes to orbital forcing.

Rubincam, David Parry

Pluto and Charon: A Case of Precession-Orbit Resonance?

Pluto may be the only known case of precession-orbit resonance in the solar system. The Pluto-Charon system orbits the Sun with a period of 1 Plutonian year, which is 250.8 Earth years. The observed parameters of the system are such that Charon may cause Pluto to precess with a period near 250.8 Earth years. This gives rise to two possible resonances, heretofore unrecognized. The first is due to Pluto's orbit being highly eccentric, giving solar torques on Charon with a period of 1 Plutonian year. Charon in turn drives Pluto near its precession period. Volatiles, which are expected to shuttle across Pluto's surface between equator and pole as Pluto's obliquity oscillates, might change the planet's dynamical flattening enough so that Pluto crosses the nearby resonance, forcing the planet's equatorial plane to depart from Charon's orbital plane. The mutual tilt can reach as much as 2 deg after integrating over 5.6 x 10(exp 6) years, depending upon how close Pluto is to the resonance and the supply of volatiles. The second resonance is due to the Sun's traveling above and below Charon's orbital plane; it has a period half that of the eccentricity resonance. Reaching this half-Plutonian year resonance requires a much larger but still theoretically possible amount of volatiles. In this case the departure of Charon from an equatorial orbit is about 1 deg after integrating for 5.6 x 10(exp 6) years. The calculations ignore libration and tidal friction. It is not presently known how large the mutual tilt can grow over the age of the solar system, but if it remains only a few degrees, then observing such small angles from a Pluto flyby mission would be difficult. It is not clear why the parameters of the Pluto-Charon system are so close to the eccentricity resonance.

Rubincam, David Parry

America's foremost early astronomer

The life of 18th century astronomer, craftsman, and partriot David Rittenhouse is detailed. As a craftsman, he distinguished himself as one of the foremost builders of clocks. He also built magnetic compasses and surveying instruments. The finest examples of his craftsmanship are considered two orreries, mechanical solar systems. In terms of astronomical observations, his best-known contribution was his observation of the transit of Venus in 1769. Rittenhouse constructed the first diffraction grating. Working as Treasurer of Pennsylvania throughout the Revolution, he became the first director of the Mint in 1792. Astronomical observations in later life included charting the position of Uranus after its discovery.

Rubincam, David Parry

Asteroid orbit evolution due to thermal drag

Thermal drag, a variant of the Yarkovsky effect, may act on small asteroids with sizes from a few meters to a few tens of meters. Yarkovsky thermal drag comes from an asteroid's absorbing sunlight in the visible and reradiating it in the infrared. Since the infrared photons have momentum, by action-reaction, they kick the asteroid when they leave its surface. The reradiation, which is asymmetric in latitude over the asteroid, gives a net force along the asteroid's pole. Due to the asteroid's thermal inertia, averaging this force over one orbital period produces a net drag if the spin axis has a component in the orbital plane. Thermal drag tends to circularize orbits. It can increase or decrease orbital inclinations. An object whose spin axis points in random directions over its lifetime displays little change in orbital inclination. Thermal drag appears to have little to do with the delivery of chondrites from the asteroid belt; the thermal drag timescale (10(exp 8) years for meter-sizzed objects) is long compared with their cosmic ray exposure ages, and aphelia in the asteroid belt are not expected for mature thermal drag orbits. However, Yarkovsky thermal drag may act on the recently discovered near-Earth asteroids, which have radii of 10-30 m. Asteroid 1992 DA, for instance, might have its orbit shrunk by 0.1 AU in 3 x 10(exp 7) years, removing it from an Earth-crossing orbit. The near-Earth asteroids also tend to have small to moderate orbital eccentricities, as expected for highly evolved thermal drag objects. However, the time needed to bring them in from the asteroid belt (about 10(exp 9) years) is long compared with the collisonal and dynamical lifetimes (both about 10(exp 8) years) for Earth-crossing objects, arguing against their emplacement by thermal drag.

Rubincam, David Parry

Variations of Mars gravitational field and rotation due to seasonal CO2 exchange

About a quarter of the Martian atmospheric mass is exchanged between the atmosphere and the polar caps in the course of a Martian year: CO2 condenses to form (or add to) the polar caps in winter and sublimes into the atmosphere in summer. This paper studies the effect of this CO2 mass redistribution on Martian rotation and gravitational field. Two mechanisms are examined: (1) the waxing and waning of solid CO2 in the polar caps and (2) the geographical distribution of gaseous CO2 in the atmosphere. In particular, the net peak-to-peak changes in J2 and J3 over a Martian year are both found to be as much as about 6 x 10 to the -9th. A simulation suggests that these changes may be detected by the upcoming Mars Observer under favorable but realistic conditions.

Chao, B. Fong

Mars - Change in axial tilt due to climate?

The average tilt of Mars' equator with respect to its orbital plane may have increased significantly over the age of the solar system. Obliquity oscillations might have induced changes in the climate, which altered the mass distribution and hence the solar torque on the planet. Viscous deformation attributable to loading by the large polar caps expected at low obliquity may have induced secular changes in the axial tilt. Earth-like effective viscosities can account for virtually the entire present obliquity of 24.4 degrees. Thus the present average tilt of Mars may not be primordial.

Rubincam, David Parry

Drag on the Lageos satellite

The drag data on the Lageos satellite collected between 1976 and 1987 are analyzed, and the cause of fluctuations observed at times when the orbit intersects the earth's shadow is investigated. The average drag on the Lageos consists of three types of drag: the Yarkovsky thermal drag, which accounts for about 70 percent of the observed drag; the neutral particle drug, accounting for about 14 percent of the drag for a hydrogen number density of 5 x 10 to the 9th/cu m; and the charged particle drag, which accounts for 12 percent of average drag at a number density of 3 x 10 to the 9th and a temperature of 5000 K in the plasmasphere for a satellite potential of -1 V. The observed drag agrees with the particle environment at the Lageos altitude and the drag models. However, the fluctuations in Lageos's along-track acceleration when the orbit intersects the earth's shadow are still not well understood.

Rubincam, David Parry