Doublet craters on Venus

Of the impact craters on Earth larger than 20 km in diameter, 10-15% (3 out of 28) are doublets, having been formed by the simultaneous impact of two well-separated projectiles. The most likely scenario for their formation is the impact of well-separated binary asteroids. If a population of binary asteroids is capable of striking the Earth, it should also be able to hit the other terrestrial planets as well. Venus is a promising planet to search for doublet craters because its surface is young, erosion is nearly nonexistent, and its crater population is significantly larger than the Earth's. After a detailed investigation of single craters separated by less than 150 km and “multiple” craters having diameters greater than 10 km, we found that the proportion of doublet craters on Venus is at most 2.2%, significantly smaller than Earth's, although several nearly incontrovertible doublets were recognized. We believe this apparent deficit relative to the Earth's doublet population is a consequence of atmospheric screening of small projectiles on Venus rather than a real difference in the population of impacting bodies. We also examined “splotches,” circular radar reflectance features in the Magellan data. Projectiles that are too small to form craters probably formed these features. After a careful study of these patterns, we believe that the proportion of doublet splotches on Venus (14%) is comparable to the proportion of doublet craters found on Earth (10-15%). Thus, given the uncertainties of interpretation and the statistics of small numbers, it appears that the doublet crater population on Venus is consistent with that of the Earth.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

Error estimates with L1 solutions

A method is given, based on the pseudoinverse of the equations of condition, to obtain error estimates for the solution in the normL ₁ of an over-determined linear system. The computational labor to obtain the errors, while not trivial, is less than that for various competing methods, particularly if there are many more equations of condition than unknowns. The error estimates for anL ₁ solution are substantially larger than those for a least squares solution of the some data. It is suggested that a complete discussion of a linear system include at least bothL ₁ and least squares solutions with their respective errors and the condition number of the linear system.     

Longevity of fluorine-bearing tremolite on Venus

There is a general belief that hydrous minerals cannot exist on Venus under current surface conditions. This view was challenged when Johnson and Fegley (2000, Icarus 146, 301-306) showed that tremolite (Ca₂Mg₅Si₈O₂₂(OH)₂), a hydrous mineral, is stable against thermal decomposition at current Venus surface temperatures, e.g., 50% decomposition in 4 Ga at 740 K. To further explore hydrous mineral thermal stability on Venus, we experimentally determined the thermal decomposition kinetics of fluorine-bearing tremolite. Fluor-tremolite is thermodynamically more stable than OH-tremolite and should decompose more slowly. However how much slower was unknown. We measured the decomposition rate of fluorine-bearing tremolite and show that its decomposition is several times to greater than ten times slower than that of OH-tremolite. We also show that F-bearing tremolite is depleted in fluorine after decomposition and that fluorine is lost as a volatile species such as HF gas. If tremolite ever formed on Venus, it would probably also contain fluorine. The exceptional stability of F-bearing tremolite strengthens our conclusions that if hydrous minerals ever formed on Venus, they could still be there today.     

Dynamics of Two Planets in the 3/2 Mean-motion Resonance: Application to the Planetary System of the Pulsar PSR B1257+12

This paper considers the dynamics of two planets, as the planets B and C of the pulsar PSR B1257+12, near a 3/2 mean-motion resonance. A two-degrees-of-freedom model, in the framework of the general three-body planar problem, is used and the solutions are analyzed through surfaces of section and Fourier techniques in the full phase space of the system.On the possibility of existence of a fourth planet in distant orbit, see Wolszczan et al., 2000     

The Common Origin of the High Inclination TNO's

Numerical integrations of the four major planets orbits inside a primordialplanetesimals disk show that a fraction of Neptune primordial scatteredobjects are deposited into the classical Kuiper Belt at Solar System age. Theseobjects exhibit inclinations as high as 40° and can account forpresent high inclinations population in the classical Kuiper Belt. The samemechanism can also originate high perihelion scattered objects like 2000 CR₁₀₅. The process that in the end produced such objects can be divided into two phases, a migration phase where nonconservative dynamics acted to producesome stable objects already at 10⁸ years and a nonmigrating phase that helped to establish some other objects as stable TNO's. Low inclination CKBO's have inprinciple an origin through the resonance sweeping process, although someresults from numerical integrations at least suggest a possible origin also fromthe primordial Neptune scattered population.     

An Example of Island-Arc Petrogenesis: Geochemistry and Petrology of the Southern Luzon Arc, Philippines

Volcanism throughout the Luzon arc is associated with eastwardsubduction of the South China Sea floor along the Manila Trench.The southern section of the arc, the focus of this study, extendsfrom the Lingayen-Dingalan fault to the small islands just southof Luzon. Two segments appear to exist along this section ofthe arc the northern Bataan and southern Mindoro segments whichare separated by the Macolod Corridor. The volcanic rocks have typical arc phenocryst mineralogies:olivine, clinopyroxene, plagioclase, and titanomagnetite inthe most mafic rocks and clinopyroxene, plagioclase, orthopyroxene,titanomagnetite, ? amphibole in the more felsic samples. Complexzoning, sieve textures, and decoupling of incompatible traceelements suggest that processes such as assimilation have takenplace. The rocks from the study area range from basalts to rhyolitesand show typical calc-alkaline features. The rocks of the MacolodCorridor and Mindoro segment are particularly enriched in largeion lithophile elements (LILE), light rare earth elements (LREE),and radiogenic Sr compared with the Bataan segment. The datafall within the mantle array on Sr-Nd isotopic diagrams andgrade toward higher Sr and lower Nd isotopic values from northto south. A likely source for the volcanics of this study is either amid-ocean ridge basalt (MORB)-type mantle that undergoes higherdegrees of partial melting than regions involved in MORB generationor a previously depleted source. We suggest that the high fieldstrength element (HFSE) anomalies have been derived throughdifferential element partitioning during fluid transport fromthe subducted lithosphere to the mantle wedge. Continental crustal material seems to play a significant roleparticularly in the Macolod Corridor and the Mindoro segment,based on the high LILE, La/Sm ratios, radiogenic Sr isotopes,and ¹⁸O values. The Macolod Corridor and the Mindoro segmenthave undergone source contamination by crustal material fromthe North Palawan-Mindoro crustal block either during the collisionof this block with the Manila Trench or by subduction of sedimentsrich in this crustal material. A similar component has alsobeen detected in the Bataan segment but in minor amounts. Thetrace element and isotopic differences between the northernand southern sections of the arc are interpreted in terms ofvariable composition (i. e., variable amounts of a crustal componentintroduced from the Palawan-Mindoro crustal terrain) of themetasomatic fluids released into the source.     

An estimate of the contribution of high nergy cosmic-ray protons to the absorbed dose inventory of a cometary nucleus

We discuss the contribution of cosmic-ray protons at all energies above 1 MeV to the absorbed doses of the surface layers of a comet. Since there exists no calculation which takes into account proton energy losses by means of losses to electrons and nuclear collisions (in a cascade process), and losses due to the low energy end-products of the cascade, we have made a rough estimate of all of these contributions. An analytical formula is proposed that allows a rapid estimate of ranges and the dose absorbed at any depth. We give dose-depth curves for two extreme values of the energy at which nuclear collisions begin to dominate the slowing-down process, and for an intermediate value we display the dose-depth curve down to 20 m from the surface. The relevance of these findings to dosimetry in comets and some alterations of cometary material are considered. The need for improving the analytical expression proposed is stressed.