An improved transformation-elimination technique for the solution of perturbed Hamiltonian systems

A Lie transformation technique proposed by the author (1977) utilizing an approach due to Kolmogorov and Arnol'd is modified to eliminate the calculation of unproductive terms. Special features of the original method — the ability to eliminate simultaneously both short- and long-period terms, and [formal] quadratic convergence for non-degenerate system sand nearly quadratic convergence for degenerate systems—are retained in an even more efficient technique. Comparisons of the new method with other available Lie transform techniques is made.     

A search for natural or artificial objects located at the Earth-Moon libration points

Photographs in the vicinity of the Earth-Moon triangular libration points L4 and L5, and of the solar-synchronized positions in the associated halo orbits (A. A. Kamel, 1969, Ph.D. dissertation, Stanford University), were made during August–September 1979, using the 30-in Cassegrain telescope at Leuschner Observatory, Lafayette, California. An effective 2° square field was covered at each position. No discrete objects, either natural or artificial, were found. The detection limit was about 14th magnitude. The present work extends traditional SETI observations to include the search for interstellar probes (R. A. Freitas, Jr., 1980a, J. Brit. Interplanet. Soc.33, 95–100).     

Venus mesosphere and thermosphere temperature structure: II. Day-night variations

A two-dimensional nonlinear hydrodynamic model has been developed for studying the global scale winds, temperature, and compositional structure of the mesosphere and thermosphere of Venus. The model is driven by absorption of solar radiation. Ultraviolet radiation produces both heating and photodissociation. Infrared solar heating and thermal cooling are also included with an accurate NLTE treatment. The most crucial uncertainty in determining the solar drive is the efficiency by which $\text{λ < 1080}\text{A}\text{?}$ solar radiation is converted to heat. This question was analyzed in Part I, where it was concluded that essentially all hot atom and O(¹D) energy may be transferred to vibrational-rotational energy of CO₂ molecules. If this is so, the minimum possible euv heating occurs and is determined by the quenching of the resulting excess rotational energy. The hydrodynamic model is integrated with this minimum heating and neglecting any small-scale vertical eddy mixing. The results are compared with predictions of another model with the same physics except that it assumes that 30% of $\text{λ < 1080}\text{A}\text{?}$ radiation goes into heat and that the heating from longer-wavelength radiation includes the O(¹D) energy. For the low-efficiency model, exospheric temperatures are ?300°K on the dayside and drop to < 180°K at the antisolar point. For the higher-efficiency model, the day-to-night temperature variation is from ?600°K to ?250°K. Both versions of the model predict a wind of several hundred meters per second blowing across the terminator and abruptly weakening to small values on the nightside with the mass flow consequently going into a strong tongue of downward motion on the nightside of the terminator. The presence of this circulation could be tested observationally by seeing if its signature can be found in temperature measurements. Both versions of the model indicate that a self-consistent large-scale circulation would maintain oxygen concentrations with ?5% mixing ratios near the dayside F-1 ionospheric peak but ?40% at the antisolar point at the same pressure level.     

Are the high-latitude torsional oscillations of the sun real?

A numerical test is made to determine if the high-latitude torsional wave is generated from the low-latitude torsional pattern as a result of our reduction procedures. The results indicate that the high-latitude motions are not an artifact of the analysis, but are true solar features. We demonstrate also that the one-wave-per-hemisphere torsional oscillation does not result from the reduction procedure. These results place the observations in conflict with the predictions of - () models of the solar cycle.Now at Institute for Astronomy, University of Hawaii, Honolulu, Hawaii 96822, U.S.A.     

Sensitivity of glaciation to initial snow cover,CO2, snow albedo,and oceanic roughness in the NCAR CCM

We present results from numerical experiments made with a GCM, the NCAR CCM1, that were designed to estimate the annual balance between snow-fall accumulation and ablation for geographically important land regions for a variety of conditions. We also attempt to assess the reliability of these results by investigating model sensitivity to changes in prescribed physical parameters. Experiments were run with an initial imposition of 1 m of (midwinter) snowcover over all northern hemisphere land points. Over Alaska, western Canada, Siberia, and the Tibetan Plateau the model tended to retain this snow cover through the summer and in some cases increase its depth as well. We define these regions as glaciation sensitive and note some correspondence between them and source regions for the Pleistocene ice sheets. An experiment with greatly reduced CO₂ (100 ppm) showed a tendency towards spontaneous glaciation, i.e., the model remained snow-covered throughout the summer over the same geographic regions noted above. With 200 ppm CO₂ (roughly equal to values at the last glacial maximum), snow cover over these regions did not quite survive the summer on a consistent basis. Combining 200 ppm CO₂ and 1 m of initial northern hemisphere snow cover yielded glaciation-sensitive conditions, agreeing remarkably well with locations undergoing glaciation during the Pleistocene. To assess the reliability of these results, we have determined minimal model uncertainty by varying two of the empirical coefficients in the model within physically plausible ranges. In one case surface roughness of all ocean gridpoints was reduced by an order of magnitude, leading to local 10% reductions in precipitation (snowfall), a change hard to distinguish from inherent model variability. In the other case, the fraction of a land grid square assumed to be occupied by snow cover for albedo purposes was varied from one-half to unity. Large changes occurred in the degree of summer melting, and in some cases the sign of the net balance changed as fractional snow cover was changed. We conclude that the model may be able to reveal regions sensitive to glaciation, but that it cannot yield a reliable quantitative computation of the magnitude of the net snow accumulation that can be implicitly or explicitly integrated through time.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dilmenil     

Magmatism in the eastern Aleutian Arc: temporal characteristic of igneous activity on Akutan Island

Lavas from Akutan Island, located in the eastern Aleutian arc at the transition between continental and oceanic crust, show a gradual change in their petrologic and chemical characteristics over the last 4 million years. The oldest lavas exposed on the island, the Hot Springs Bay Volcanics (HSBV), range from magnesian basalt to dacite (45%–62% SiO₂). The most mafic basalts contain salitic clinopyroxene, Cr- and Al-rich spinel, and pargasitic amphibole suggesting that they were derived from relatively hydrous magmas at greater pressures than lavas from the younger Akutan Volcanics (AKV) and the modern volcano (MOD). AKV lavas also range between basalt and dacite (46%–63% SiO₂), but contain no hydrous phenocrysts and seem to have fractionated within a shallow level magma chamber. Lavas from the modern volcano are andesitic (52%–57% SiO₂) and have a mineral assemblage similar to that of AKV lavas of similar composition. With the exception of clinopyroxene and spinel in the most mafic lavas, the compositions of plagioclase (An_92?45), olivine (Fo_88?51), orthopyroxene (En_69?56), and titanomagnetite (15%–21% TiO₂) phenocrysts found in these lavas are within the range observed in lavas from other Aleutian volcanoes. Variations in the major element chemistry of the older lavas can be reproduced by fractional crystallization of the observed mineral assemblages, however closed system crystal fractionation models are inadequate to explain the trace element variations. During the last 4 million years, La/Yb ratios have decreased (6.5–3.3 for HSBV lavas and 2.9–1.9 for MOD lavas) whereas Ba/La ratios appear to have increased slightly (37–43 for HSBV and AKV, and 41–45 of MOD). The lower La/Yb ratios of MOD lavas correspond with lower total abundances of the REE and slightly higher Sr and Pb isotopic ratios. The increased⁸⁷Sr/⁸⁶Sr ratios and Pb isotopic ratios in the MOD lavas, the less enriched LREE, and the higher Ba/La ratios may result from partial melting of an arc source which has experienced previous melting events but has continued to be contaminated by a component from the subducting slab. It may also indicate a change in the degree of partial melting of the underlying mantle, which corresponds to a different percentage of a slab derived component being incorporated into the overlying mantle.     

Peak discharge of a Pleistocene lava-dam outburst flood in Grand Canyon, Arizona, USA

The failure of a lava dam 165,000 yr ago produced the largest known flood on the Colorado River in Grand Canyon. The Hyaloclastite Dam was up to 366 m high, and geochemical evidence linked this structure to outburst-flood deposits that occurred for 32 km downstream. Using the Hyaloclastite outburst-flood deposits as paleostage indicators, we used dam-failure and unsteady flow modeling to estimate a peak discharge and flow hydrograph. Failure of the Hyaloclastite Dam released a maximum 11 × 10⁹ m³ of water in 31 h. Peak discharges, estimated from uncertainty in channel geometry, dam height, and hydraulic characteristics, ranged from 2.3 to 5.3 × 10⁵ m³ s⁻¹ for the Hyaloclastite outburst flood. This discharge is an order of magnitude greater than the largest known discharge on the Colorado River (1.4 × 10⁴ m³ s⁻¹) and the largest peak discharge resulting from failure of a constructed dam in the USA (6.5 × 10⁴ m³ s⁻¹). Moreover, the Hyaloclastite outburst flood is the oldest documented Quaternary flood and one of the largest to have occurred in the continental USA. The peak discharge for this flood ranks in the top 30 floods (>10⁵ m³ s⁻¹) known worldwide and in the top ten largest floods in North America.