Search

Microbialites are organosedimentary structures that can be constructed by a variety of metabolically distinct taxa. Consequently, microbialite structures abound in the fossil record, although the exact nature of the biogeochemical processes that produced them is often unknown. One such class of ancient calcareous structures, Epiphyton and Girvanella, appear in great abundance during the Early Cambrian. Together with Archeocyathids, stromatolites and thrombolites, they formed major Cambrian reef belts. To a large extent, Middle to Late Cambrian reefs are similar to Precambrian reefs, with the exception that the latter, including terminal Proterozoic reefs, do not contain Epiphyton or Girvanella. Here we report the discovery in Pavilion Lake, British Columbia, Canada, of a distinctive assemblage of freshwater calcite microbialites, some of which display microstructures similar to the fabrics displayed by Epiphyton and Girvanella. The morphologies of the modern microbialites vary with depth, and dendritic microstructures of the deep water (> 30 m) mounds indicate that they may be modern analogues for the ancient calcareous structures. These microbialites thus provide an opportunity to study the biogeochemical interactions that produce fabrics similar to those of some enigmatic Early Cambrian reef structures.

We assemble data on the Pinatubo aerosol from space, air, and ground measurements, develop a composite picture, and assess the consistency and uncertainties of measurement and retrieval techniques. Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements agree with theoretical calculations in showing a dominant composition of H2SO4-H20 mixture, with H2SO4 weight fraction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are (1) volcanic ash, present at all heights initially and just above the tropopause until at least March 1992, and (2) much smaller H2SO4 fractions at the low temperatures of high-latitude winters and the tropical tropopause. Laboratory spectroscopy and calculations yield wavelength- and temperature-dependent refractive indices for the H2SO4-H20 droplets. These permit derivation of particle size information from measured optical depth spectra, for comparison to impactor and optical-counter measurements. All three techniques paint a generally consistent picture of the evolution of R(sub eff), the effective radius. In the first month after the eruption, although particle numbers increased greatly, R(sub eff) outside the tropical core was similar to preeruption values of approx. 0.1 to 0.2 microns, because numbers of both small (r < 0.2 microns) and large (r > 0.6 microns) particles increased. In the next 3-6 months, extracore R(sub eff) increased to approx. 0.5 microns, reflecting particle growth through condensation and coagulation. Most data show that R(sub eff) continued to increase for about 1 year after the eruption. R(sub eff) values up to 0.6 - 0.8 microns or more are consistent with 0.38 - 1 micron optical depth spectra in middle to late 1992 and even later. However, in this period, values from in situ measurements are somewhat less. The difference might reflect in situ undersampling of the very few largest particles, insensitivity of optical depth spectra to the smallest particles, or the inability of flat spectra to place an upper limit on particle size. Optical depth spectra extending to wavelengths lambda > 1 micron are required to better constrain R(sub eff), especially for R(sub eff) > 0.4 microns. Extinction spectra computed from in situ size distributions are consistent with optical depth measurements; both show initial spectra with lambda(sub max) <= 0.42 microns, thereafter increasing to 0.78 <= lambda(sub max) <= 1 micron. Not until 1993 do spectra begin to show a clear return to the preeruption signature of lambda(sub max) <= 0.42 microns. The twin signatures of large R(sub eff) (> 0.3 microns) and relatively flat extinction spectra (0.4 - 1 microns) are among the longest-lived indicators of Pinatubo volcanic influence. They persist for years after the peaks in number, mass, surface area, and optical depth at all wavelengths <= 1 microns. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between global maps of 0.5- and 1.0-micron optical depth derived from the Advanced Very High Resolution Radiometer (AVHRR) and Stratospheric Aerosol and Gas Experiment (SAGE) satellite sensors. However, there are important differences between the AVHRR and SAGE midvisible optical thickness products. We discuss possible reasons for these differences and how they might be resolved.

A new wavelength-dependent model of the single-scattering properties of the Martian dust is presented. The model encompasses the solar wavelengths (0.3 to 4.3 micrometers at 0.02 micrometer resolution) and does not assume a particular mineralogical composition of the particles. We use the particle size distribution, shape, and single-scattering properties at Viking Lander wavelengths presented by Pollack et al. [1995]. We expand the wavelength range of the aerosol model by assuming that the atmospheric dust complex index of refraction is the same as that of dust particles in the bright surface geologic units. The new wavelength-dependent model is compared to observations taken by the Viking Orbiter Infrared Thermal Mapper solar channel instrument during two dust storms. The model accurately matches afternoon observations and some morning observations. Some of the early morning observations are much brighter than the model results. The increased reflectance can be ascribed to the formation of a water ice shell around the dust particles, thus creating the water ice clouds which Colburn et al. [1989], among others, have predicted.

As the dense molecular cloud that was the precursor of our Solar System was collapsing to form a protosun and the surrounding solar-nebula accretion disk, infalling interstellar grains were heated much more effectively by radiation from the forming protosun than by radiation from the disk's accretion shock. Accordingly, we have estimated the temperatures experienced by these infalling grains using radiative diffusion calculations whose sole energy source is radiation from the protosun. Although the calculations are 1-dimensional, they make use of 2-D, cylindrically symmetric models of the density structure of a collapsing, rotating cloud. The temperature calculations also utilize recent models for the composition and radiative properties of interstellar grains (Pollack et al. 1994. Astrophys. J. 421, 615-639), thereby allowing us to estimate which grain species might have survived, intact, to the disk accretion shock and what accretion rates and molecular-cloud rotation rates aid that survival. Not surprisingly, we find that the large uncertainties in the free parameter values allow a wide range of grain-survival results: (1) For physically plausible high accretion rates or low rotation rates (which produce small accretion disks), all of the infalling grain species, even the refractory silicates and iron, will vaporize in the protosun's radiation field before reaching the disk accretion shock. (2) For equally plausible low accretion rates or high rotation rates (which produce large accretion disks), all non-ice species, even volatile organics, will survive intact to the disk accretion shock. These grain-survival conclusions are subject to several limitations which need to be addressed by future, more sophisticated radiative-transfer models. Nevertheless, our results can serve as useful inputs to models of the processing that interstellar grains undergo at the solar nebula's accretion shock, and thus help address the broader question of interstellar inheritance in the solar nebula and present Solar System. These results may also help constrain the size of the accretion disk: for example, if we require that the calculations produce partial survival of organic grains into the solar nebula, we infer that some material entered the disk intact at distances comparable to or greater than a few AU. Intriguingly, this is comparable to the heliocentric distance that separates the C-rich outer parts of the current Solar System from the C-poor inner regions.

High-resolution, near-infrared (1.09 to 2.5 micrometers) spectra of the night side of Venus have been obtained in 1990 and 1991 using the Fourier Transform Spectrometer at the 3.6-m Canada-France-Hawaii telescope. Absorptions due to H2O were detected in spectral windows near 2.3, 1.74, and 1.18 micrometers. Our analysis of these absorptions constrains the abundance of water vapor in three different altitude ranges located between the clouds and the surface: 30-40 km, 15-25 km and 0-15 km. A constant water vapor mixing ratio of 30 +/- 15 ppm below the clouds can fit the observations. These values are consistent with recent near-infrared studies of the night side of Venus at lower spectral resolution. The atmosphere of Venus appears to be dryer than originally suggested by the in-situ measurements made by the Pioneer Venus and Venera mass-spectrometers and gas-chromatographs.

We combine space, air, and ground measurements to develop a composite picture of the post-Pinatubo aerosol, and assess the consistency and uncertainties of various measurement and retrieval techniques. impactor and optical counter measurements, as well as retrievals from optical depth spectra, paint a generally consistent picture of the evolution of particle effective radii, R(sub eff). In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the preeruption value of 4.2 micrometers, because both small (r less than 0.25 micrometers) and large (r greater than 0.6 micrometers) particles increased in number, Over the next 3-6 months, R(sub eff) increased rapidly to about 0.5 micrometers. In general, R(sub eff) continued to increase for about a year after the eruption. The peak wavelength of optical depth spectra increased from initial values of less than 0.42 micrometers to values between 0.78 and 1 micrometer. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5 and 1.0-micrometer optical depth derived from the AVHRR and SAGE satellite measurements. It also sets a context for evaluating remaining uncertainties in each of these satellite data products. We also make consensus recommendations for particle composition, shape, and temperature- and wavelength-dependent refractive index, and show how the latter effect on backscatter spectra can influence particle sizes retrieved from multiwavelength lidar measurements.

We combine a variety of measurements to develop a composite picture of the post-Pinatubo aerosol and assess the consistency and uncertainties of the measurement and retrieval techniques. Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements are in accord with theoretical calculations in showing a dominant particle composition of H2SO4-H2O mixture, with H2SO4 weight fraction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are: (1) the presence of volcanic ash at all altitudes initially and in a layer just above the tropopause until at least March 1992, and (2) much smaller H2SO4 weight fractions at the low temperatures attained in high latitude winters and at the tropical tropopause, Laboratory spectroscopy and theoretical calculations yield wavelength- and temperature-dependent refractive indices for the dominant H2SO4-H2O droplets. These in turn permit derivation of particle size spectra from measured optical depth spectra, for comparison to direct measurements by impactors and optical counters. All three techniques paint a generally consistent picture of the evolution of R(sub eff), the effective, or area-weighted, particle radius. In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the preemption value of 0.1 to 0.2 microns, because both small (r less than 0.2 microns) and large (r greater than 0.6 micron particles increased in number. Over the next 3-6 months, R(sub eff) increased to about 0.5 microns reflecting particle growth through condensation and coagulation. In general, R(sub eff) continued to increase for about a year after the eruption. Extinction spectra computed from in situ size distribution measurements are consistent with optical depth measurements, which show spectra with maxima initially at wavelengths of 0.42 microns or less, and thereafter progressively increasing to between 0.78 and 1 micron. Not until 1993 do optical depth spectra begin to show a clear return to the preemption signature of maximizing at the shortest visible wavelengths or in the near UV. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5- 1.0- micron optical depth derived from the AVHRR and SAGE satellite measurements.

We combine a variety of measurements to develop a composite picture of the post-Pinatubo aerosol and assess the consistency and uncertainties of the measurement and retrieval techniques Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements are in accord with theoretical calculations In showing a dominant particle composition of H2SO4-H2O mature with H2SO4 weight traction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are: (1) the presence of volcanic ash at all attitudes initially and in a layer just above the tropopause until at least March 1992, and (2) much smaller H2SO4 weight fractions at the low temperatures attained In high latitude winters and at the tropical tropopause. Laboratory spectroscopy and theoretical calculations yield wavelength- and temperature-dependent refractive indices for the dominant H2SO4-H2O droplets. These in turn permit derivation of particle size spectra from measured optical depth spectra for comparison to direct measurements by impactors and optical counters, All three techniques paint a generally consistent picture of the evolution of R(sub eff) the effective, or area-weighted, particle radius. In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the pre-eruption value of 0.1 to 0.2 microns because both small (r less than 0.2 microns) and large (r greater than 0.6 microns) particles increased in number. Over the next 3-6 months, R(sub eff) increased to about 0.5 microns, reflecting particle growth through condensation and coagulation. In general, R(sub eff) continued to increase for about a year after the eruption. Extinction spectra computed from in situ size distribution measurements are consistent with optical depth measurements, which show spectra with maxima initially at wavelengths of 0.42 microns or less, and thereafter progressively increasing to between 0.78 and 1 micron. Not until 1993 does optical depth spectra begin to show a clear return to the preeruption signature of maximizing at the shortest visible wavelengths or in the near UV. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5- and 1.0-kilometer optical depth derived from the AVHRR and SAGE satellite measurements. It also sets a context for evaluating remaining uncertainties in each of these satellite data products. We also show how the effects of wavelength-dependent refractive index on backscatter spectra can influence particle sizes retrieved from multiwavelength lidar measurements.

We present moderate-resolution (lambda/delta lambda = 300 to 370) reflectance spectra of Mars from 2.04 to 2.44 micrometers that were obtained at UKIRT during the 1993 opposition. Seven narrow absorption features were detected and found to have a Mars origin. By comparison with solar and Mars atmospheric spectra, five of these features were attributed all or in part to Mars atmospheric CO2 or CO(2.052 +/- 0.003, 2.114 +/- 0.002, 2.150 +/- 0.003, 2.331 +/- 0.001, and 2.357 +/- 0.002 micrometers). Two of the bands (2.331 +/- 0.001 and 2.357 +/- 0.002 micrometers) appear to have widths and depths that are consistent with additional, nonatmospheric absorptions, although a solar contribution cannot be entirely ruled out. Two other weak bands centered at 2.278 +/- 0.002 and 2.296 +/- 0.002 micrometers may be at least partially mineralogic in origin. The data provide no conclusive identification of the mineralogy responsible for these absorption features. However, examination of terrestrial spectral libraries and previous moderate spectral resolution mineral studies indicates that the most likely origin of these features is either (bi)carbonate or (bi)sulfate anions in framework silicates or (Fe, Mg)-OH bonds in sheet silicates. If the bands are caused by phyllosilicate minerals, then an explanation must be found for the extremely narrow widths of the cation-OH features in the Mars spectra as compared to terrestrial minerals.

A coupled problem of diffusion and condensation is solved for the H2SO4-H2O system in Venus' cloud layer. The position of the lower cloud boundary and profiles of the H2O and H2SO4 vapor mixing ratios and of the H2O/H2SO4 ratio of sulfuric acid aerosol and its flux are calculated as functions of the column photochemical production rate of sulfuric acid, Phi(sub H2SO4). Variations of the lower cloud boundary are considered. Our basic model, which is constrained to yield f(sub H2O) (30 km) = 30 ppm (Pollack et al. 1993), predicts the position of the lower cloud boundary at 48.4 km coinciding with the mean Pioneer Venus value, the peak H2SO4 mixing ratio of 5.4 ppm, and the H2SO4 production rate Phi(sub H2SO4) = 2.2 x 10(exp 12)/sq cm/s. The sulfur to sulfuric acid mass flux ratio in the clouds is 1:27 in this model, and the mass loading ratio may be larger than this value if sulfur particles are smaller than those of sulfuric acid. The model suggests that the extinction coefficient of sulfuric acid particles with radius 3.7 micrometers (mode 3) is equal to 0.3/km in the middle cloud layer. The downward flux of CO is equal to 1.7 x 10(exp 12)/sq cm/s in this model. Our second model, which is constrained to yield f(sub H2O) = 10 ppm at the lower cloud boundary, close to the value measured by the Magellan radiooccultations, predicts the position of this boundary to be at 46.5 km, which agrees with the Magellan data; f(sub H2O) (30 km) = 90 ppm, close to the data of Moroz et al. (1983) at this altitude; Phi(sub H2SO4) = 6.4 x 10(exp 12)/sq cm/s; and Phi(sub co) = 4.2 x 10(exp 12)/sq cm/s. The S/H2SO4 flux mass ratio is 1:18, and the extinction coefficient of the mode 3 sulfuric acid particles is equal to 0.9 km in the middle cloud layer. A strong gradient of the H2SO4 vapor mixing ratio near the bottom of the cloud layer drives a large upward flux of H2SO4, which condenses and forms the excessive downward flux of liquid sulfuric acid, which is larger by a factor of 4-7 than the flux in the middle cloud layer. This is the mechanism of formation of the lower cloud layer. Variations of the lower cloud layer are discussed. Our modeling of the OCS and CO profiles in the lower atmosphere measured by Pollack et al. (1993) provides a reasonable explanation of these data and shows that the rate coefficient of the reaction SO3 + CO yields CO2 + SO2 is equal to 10(exp -11) exp(-(13,100 +/- 1000)/T) cc/s. The main channel of the reaction between SO3 and OCS is CO2 + (SO)2, and its rate coefficient is equal to 10(exp -11) exp(-(8900 +/- 500)/T) cc/s. In the conditions of Venus' lower atmosphere, (SO)2 is removed by the reaction (SO)2 + OCS yields CO + S2 + SO2. The model predicts an OCS mixing ratio of 28 ppm near the surface.

The climate of Venus is to a large degree controlled by the radiative properties of its massive atmosphere. In addition, outgassing due to volcanic activity, exospheric escape processes, and surface/atmosphere interactions may all be important in moderating the abundances of atmospheric CO2 and other volatiles. We have developed an evolutionary climate model for Venus using a systems approach that emphasizes feedbacks between elements in the climate system. Modules for atmospheric radiative transfer, surface/atmosphere interactions, tropospheric chemistry, and exospheric escape processes have so far been developed. Climate feedback loops result from interconnections between modules, in the form of the environmental parameters pressure, temperature, and atmospheric mixing ratios. The radiative transfer module has been implemented by using Rosseland mean opacities in a one dimensional grey radiative-convective model. The model has been solved for the static (time independent) case to determine climate equilibrium points. The dynamics of the model have also been explored by employing reaction/diffusion kinetics for possible surface atmosphere heterogeneous reactions over geologic timescales. It was found that under current conditions, the model predicts that the climate of Venus is at or near an unstable equilibrium point. The effects of constant rate volcanism and corresponding exsolution of volatiles on the stability of the climate model were also explored.

The rate at which the wind can redistribute sedimentary material is an important part of any planet's sedimentologic cycle, particularly for Mars, where the competing effects of other gradational processes are less than on Earth. The aeolian drift potential (DP) is a measure of the amount of material capable of being moved through a unit length by the wind for a given period of time. DP is a useful measure of the potential redistribution rate of windblown material on regional scales. The Martian aeolian DP was calculated from laboratory studies of sand movement conducted at Martian atmospheric densities and from surface stress, temperature, and pressure values for that region as determined from the Mars General (Atmospheric) Circulation Model (GCM) developed at the NASA/Ames Research Center. In our simulations for Mars, DP changes in both magnitude (as expected) and direction if the saltation threshold is altered.

A coupled problem of diffusion and condensation is solved for the H2SO4-H2O system in Venus' cloud layer. The position of the lower cloud boundary, profiles of the H2O and H2SO vapor mixing ratios, and of the H2O/H2SO4 ratio of sulfuric acid aerosol and its flux are calculated as functions of the column photochemical production rate of sulfuric acid. Variations of the lower cloud boundary are considered. Our basic model, which is constrained to yield f(sub H2O)(30km) = 30 ppm, predicts the position of the lower cloud boundary at 48.4 km coinciding with the mean Pioneer Venus value, the peak H2SO4 mixing ratio of 5.4 ppm, and the H2SO4 production rate phi(sub H2SO4) = 2.2 x 10(exp 12) cm(exp -2)s(exp -1).

The spectral image cubes obtained by the Near-Infrared Mapping Spectrometer (NIMS) on Galileo as it flew by Venus have been analyzed to constrain the vertical structure of the clouds, the nature of the aerosol particles, and the location and particle properties of the opacity variations responsible for high-contrast features observed in the near-infrared windows at 1.7 and 2.3 micrometers. A radiative transfer program was used to simulate mid-latitude curves of limb darkening at 3.7 micrometers. Best-fit models to these curves demonstrate that the upper clouds are dominated by mode 2 particles (r-bar = 1.0 micrometers), with a contribution of approximately 15% of opacity from mode 1 particles (r-bar = 0.3 micrometers). The low-latitude upper cloud is well represented by a dual scale-height model, with a particle scale height of approximately 1 km from an altitude of 61-63 km, and a scale height of approximately 6 km above this, up to the level where tau = 1 at approximately 71 km. This model also successfully simulates limb-darkening curves at 11.5 micrometers from the Pioneer Venus Orbiter Infrared Radiometer. Successful simulations of correlation plots of 1.7 vs 2.3 micrometers intensities reveal that mode 3 particles (r-bar = 3.65 micrometers) represent the dominant source of opacity in the lower and middle clouds, and that variation in total cloud opacity reflects chiefly the addition and removal of mode 3 particles near the cloud base. We find that the full spectrum of brightnesses at 1.7 and 2.3 micrometers implies that the total cloud optical depth varies from approximately 25 to approximately 40.

Using Venus nightside data obtained by the Galileo Near-Infrared Mapping Spectrometer (NIMS), we have studied the correlation of 1.74 and 2.30 micrometer radiation which is transmitted through the clouds. Since the scattering and absorption properties of the cloud particles are different at these two wavelengths, one can distinguish between abundance variations and variations in the properties of the cloud particles themselves. The correlation of intensities shows a clustering of data into five distinct branches. Using radiative transfer calculations, we interpret these branches as regions of distinct but different mixes of Mode 2' and 3 particles. The data and calculations indicate large differences in these modal ratios, the active cloud regions varying in content from nearly pure Mode 2' particles to almost wholly Mode 3. The spatial distribution of these branches shows large scale sizes and both hemispheric symmetries and asymmetries. High-latitude concentrations of large particles are seen in both hemispheres and there is banded structure of small particles seen in both the North and South which may be related. The mean particle size in the Northern Hemisphere is greater than found in the South. If these different branch regions are due to mixing of vertically stratified source regions (e.g. photochemical and condensation source mechanisms), then the mixing must be coherent over very large spatial scales.

The spectroscopic data of the Near-Infrared Mapping Spectrometer (NIMS), recorded during the Galileo flyby of Venus, are analysed to retrieve the water vapor abundance variations in the lower atmosphere of Venus at night. The 1.18 micrometer spectral window, which probes altitude levels below 20 km, is used for this purpose. Constraints on the CO2 continuum and far-wing opacity from existing ground-based high-resolution observations are included in the modelling of the NIMS spectra. The NIMS measurements can be fitted with a water vapor mixing ratio of 30 +/- 15 ppm, in agreement with analyses of ground-based nightside observations. The water vapor abundance shows no horizontal variations exceeding 20% over a wide latitude range (40 deg S, 50 deg N) on the nightside of Venus. Within the same selection of NIMS spectra, a large enhancement in the O2 fluorescence emission at 1.27 micrometer is observed at a latitude of 40 deg S, over a spatial area about 100 km wide.

Numerical calculations of the linear response of a 2D gaseous disk to the perturbations induced by a protoplanet and the corresponding torque are presented. When the pressure gradient is taken into account, torques are increased in disks with gradients in either surface density and sound speed, the effect of the latter being much greater for the same-sized gradient as measured by the power law index. The torques in turn may be used to calculate timescales for orbital migration of protoplanets.

Surface pressure measurement performed by the Viking landers show substantial variations in pressure on seasonal timescales that are characterized by two local minima and two local maxima. These variations have widely been attributed to the seasonal condensation and sublimation of CO2 in the two polar regions. It has been somewhat of a surprise that the amplitude of the minimum and maximum that is dominated by the CO2 cycle in the north was much weaker than the corresponding amplitude of the south-dominated extrema. Another surprise was that the seasonal pressure cycle during years 2 and 3 of the Viking mission was so similar to that for year 1, despite the occurrence of two global dust storms during year 1 and none during years 2 and 3. An energy balance model that incorporates dynamical factors from general circulation model (GCM) runs in which the atmospheric dust opacity and seasonal date were systematically varied was used to model the observed seasonal pressure variations. The energy balance takes account of the following processes in determining the rates of CO2 condensation and sublimation at each longitudinal and latitudinal grid point: solar radiation, infrared radiation from the atmosphere and surface, subsurface heat conduction, and atmospheric heat advection. Condensation rates are calculated both at the surface and in the atmosphere. In addition, the energy balance model also incorporates information from the GCM runs on seasonal redistribution of surface pressure across the globe. Estimates of surface temperature of the seasonal CO2 caps were used to define the infrared radiative losses from the seasonal polar caps. The seasonal pressure variations measured at the Viking lander sites were closely reproduced.

Quasistationary eddies are prominent in a large set of simulations of the Mars general circulation performed with the NASA-Ames GCM. Various spacecraft observations have at least hinted at the existence of such eddies in the Mars atmosphere. The GCM stationary eddies appear to be forced primarily by the large Mars topography, and (to a much lesser degree) by spatial variations in the surface albedo and thermal inertia. The stationary eddy circulations exhibit largest amplitudes at high altitudes (above 30-40 km) in the winter extratropical regions. In these regions they are of planetary scale, characterized largely by zonal wavenumbers 1 and 2. Southern Hemisphere winter appears to be dominated by a very strong wave 1 pattern, with both waves 1 and 2 being prominent in the Northern Hemisphere winter regime. This difference seems to be basically understandable in terms of differences in the topography in the two hemispheres. The stationary eddies in the northern winter extratropics are found to increase in amplitude with dust loading. This behavior appears to be at least partly associated with changes in the structure of the zonal-mean flow that favor a greater response to wave 1 topographic forcing. There are also strong stationary eddy circulations in the tropics and in the summer hemisphere. The eddies in the summer subtropics and extratropics arc substantially stronger in southern summer than in northern summer. The summer hemisphere stationary circulations are relatively shallow and are characterized by smaller zonal scales than those in the winter extratropics.

Using Venus nightside data obtained by the Galileo Near Infrared Mapping spectrometer, we hav studied the correclation of 1.74 and 2.30 (sub m) radiation which is transmitted through the clouds.

Most studies of active aeolian processes on Mars have focused on dust, i.e., particles approximately 1 micron in diameter that are transported in suspension by wind. The presence of sand dunes on Mars indicates that larger grains (approximately greater than 60 microns, transported primarily in saltation) are also present. Although indirect evidence suggests that some dunes may be active, definitive evidence is lacking. Nonetheless, numerous studies demonstrate that sand is substantially easier to transport by wind than dust, and it is reasonable to infer that sand transportation in saltation occurs under present Martian conditions. In order to assess potential source regions, transportation pathways, and sites of deposition for sand on Mars, an iterative sand transport algorithm was developed that is based on the Mars General Circulation Model of Pollack et al. The results of the dust transport model are then compared with observed surface features, such as dune field locations observed on images, and surficial deposits as inferred from Viking IRTM observations. Preliminary results suggest that the north polar dune fields in the vicinity of 270 degrees W, 70 degrees N originated from weathered polar layered plains centered at 280 degrees W, 85 degrees N, and that Thaumasia Fossae, southern Hellas Planitia, and the area west of Hellespontus Montes are sand depositional sites. Examples of transportation 'corridors' include a westward pathway in the latitudinal band 35 degrees N to 45 degrees N, and a pathway southward from Solis Planum to Thaumasia Fossae, among others.

One of the key radiative agents in the atmosphere of Mars is the suspended dust particles. A new analysis of two data sets of the Martian atmosphere is being carried out in order to better evaluate the radiative properties of the atmospheric dust particles. The properties of interest are the size distribution, optical constants, and other radiative properties, such as the single-scattering albedo and phase function. Of prime importance is the wavelength dependence of these radiative properties throughout the visible and near-infrared wavelengths. Understanding the wavelength dependence of absorption and scattering characteristics will provide a good definition of the influence that the atmospheric dust has on heating of the atmosphere.

The contribution of various factors to the thermal structure of Titan's past and present atmosphere are discussed. A one dimensional model of Titan's thermal structure is summarized. The greenhouse effect of Titan's atmosphere, caused primarily by pressure induced opacity of N2, CH4, and H2, is discussed together with the antigreenhouse effect dominated by the haze which absorbs incident sunlight. The implications for the atmosphere of the presence of an ocean on Titan are also discussed.

The presence of suspended dust in the Martian atmosphere, and its return to the planet's surface, is implicated in the formation of the polar layered terrain and the dichotomy in perennial CO2 polar cap retention in the two hemispheres. A three dimensional model was used to study Martian global dust storms. The model accounts for the interactive feedbacks between the atmospheric thermal and dynamical states and an evolving radiatively active suspended dust load. Results from dust storm experiments, as well as from simulations in which there is interest in identifying the conditions under which surface dust lifting occurs at various locations and times, indicate that dust transport due to atmospheric eddy motions is likely to be important in the arrival of suspended dust at polar latitudes. The layered terrain in both polar regions of Mars is interpreted as the reality of cyclical episodes of volatile (CO2, H2O) and dust deposition.

The interaction of three physical processes will determine the stability of the Martian polar caps as the surface pressure increases: the greenhouse effect, atmospheric heat transport, and the change in the CO2 frost point temperature. The contribution of each is readily determined in the Mars general circulation model (GCM). Therefore, we have initiated experiments with the GCM to determine how these processes interact, and how the atmosphere-polar cap system responds to increasing surface pressure. The experiments are carried out for northern winter solstice and generally assume the atmosphere to be free of dust. Each experiment starts from resting isothermal conditions and runs for 50 Mars days. Mars' current orbital parameters are used. The experiments are for surface pressures of 120, 480, and 960 mb, which represent 16, 64, and 128 times the current value. To date we have analyzed the 120 mb experiment and the results indicate the contrary to the simpler models, the polar caps actually advance instead of retreat. Other aspects of this investigation are presented.

Runoff channels and valley networks present on ancient, heavily cratered Martian terrain suggests that the climate of Mars was originally warm and wet. One explanation for the formation of these channels is that the surface was warmed by the greenhouse effect of a dense, CO2 atmosphere. However, recent work shows that this theory is not consistent for the early period of the solar system. One way to increase the surface temperature predicted is to assume that other greenhouse gases were present in Mars' atmosphere in addition to CO2 and H2O. This possible gas is ammonia, NH3. If ammonia was present in sufficient quantities, it could have raised the surface temperature to 273 K. An adequate source would have been volcanic outgassing if the NH3 produced was shielded from photolysis by an ultraviolet light absorber.

In time scales of a hundred thousand to a million years, the eccentricity of Mars orbit varies in a quasi-periodic manner between extremes as large as 0.14 and as small as 0 and the tilt of its axis of rotation with respect to the orbit normal also varies quasi-periodically between extremes as large as 35 deg and as small as 15 deg. In addition, the orientation of the axis precesses on comparable time scales. These astronomical variations are much more extreme than those experienced by the Earth. These variations are thought to have strongly modulated the seasonal cycles of dust, carbon dioxide, and water. One manifestation of the induced quasiperiodic climate changes may be the layered terrain of the polar regions, with individual layers perhaps recording variations in the absolute and/or relative deposition rates of dust and water in the polar regions, most likely in association with the winter time deposition of carbon dioxide ice. In an attempt to understand the manner in which atmospheric temperatures and winds respond to the astronomical forcings, we have initiated a series of numerical experiments with the NASA/Ames general circulation model of the Martian Atmosphere.

The 5-25 micron real and imaginary refraction indices are presented for palagonite and the silicates pyrophyllite, kaolinite, serpentine, montmorillonite, saponite, and orthopyroxene. Optical constants in the region of the H2O-bending fundamental near 6 microns are obtained for saponite, montmorillonite, and palagonite. It is established that, if a pellet of pure material can be polished to a mirror finish, the optical constants of such noncohesive materials as clays are easily derivable.

Recent airborne thermal infrared observations of Mars from the Kuiper Airborne Observatory (KAO) have provided evidence for the presence of carbonates, sulfates, and hydrates. Using the optical properties of calcite and anhydrite, it was estimated that CO3's and SO4's constituted about 1 to 3 and 10 to 15 wt. percent, repectively of the materials composing the atmospheric dust. Using the derived value as an estimate of total CO3 abundance, and making an assumption that the CO3's were uniformly distributed within the Martian regolith, it was estimated that such a CO3 reservoir could contain roughly 2 to 5 bars of CO2. While the results indicate that several volatile-bearing materials are present on Mars, the observations from the KAO are inherently limited in their ability to determine the spatial distributions of these materials. However, previous spacecraft observations of Mars provide both the spectral coverage necessary to identify these materials, as well as the potential for investigating their spatial variability. This has prompted us to pursue a reinvestigation of the Mariner 6 and 7 infrared spectrometer and Mariner 9 infrared interferometer spectrometer observations. The former data have been recently made available in digital format and calibration of wavelengths and intensities are almost complete. Additionally, we are pursuing the derivation of optical constants of more appropriate carbonates and sulfates.

The Cassini Mission is a joint undertaking of NASA and the European Space Agency (ESA) to explore the Saturnian System with a Saturn Orbiter and a Titan Probe. The launch vehicle and the Saturn Orbiter are the responsibility of NASA while the Huygens Probe (detachable Titan Probe) is the responsibility of ESA. The spacecraft will be launched in 1996 and the Huygens Probe will arrive at Titan in 2003. The Cassini Mission-Huygens Probe provides a unique opportunity to obtain detailed information about the atmosphere and, possibly, the surface of Titan. Titan possesses a substantial nitrogen atmosphere containing methane and many other organic compounds. Aerosols play an important role in the atmospheric processes on Titan. The Huygens Probe offers an opportunity to determine how organic particles are formed and grow which will clarify their role on Earth. A powerful analytical instrument, capable of addressing the above technology and other science questions, was recently proposed for the Huygens Probe. It is comprised of an aerosol and gas sampler and processor, and a gas chromatograph-ion mobility spectrometer. The instrument will be able to measure complex organics that make up the collected aerosols to the approximate 1 ppm level. Gases will be measured to approximately 10 ppb. Because the Titan atmosphere is expected to be quite complex, a gas chromatograph-ion mobility spectrometer is used to provide unequivocal identification of the components of the analytes. Further details of the science question to be investigated and the proposed instrument are described. The expected results and their implications are also addressed.

During the 1990 Galileo Venus flyby, the Near Infrared Mapping Spectrometer investigated the night-side atmosphere of Venus in the spectral range 0.7 to 5.2 micrometers. Multispectral images at high spatial resolution indicate substantial cloud opacity variations in the lower cloud levels, centered at 50 kilometers altitude. Zonal and meridional winds were derived for this level and are consistent with motion of the upper branch of a Hadley cell. Northern and southern hemisphere clouds appear to be markedly different. Spectral profiles were used to derive lower atmosphere abundances of water vapor and other species.

Near-IR images and spectra of the night side of Venus taken at the Anglo-Australian Telescope during February 1990 reveal four new thermal emission windows at 1.10, 1.18, 1.27, and 1.31 microns, in addition to the previously discovered windows at 1.74 and 2.3 microns. Images of the Venus night side show similar bright and dark markings in all windows, but their contrast is much lower at short wavelengths. The 1.27-micron window includes a bright, high-altitude O2 airglow feature in addition to a thermal contribution from the deep atmosphere. Simulations of the 1.27- and 2.3-micron spectra indicate water vapor mixing ratios near 40 + or - 20 ppm by volume between the surface and the cloud base.

A 1D quasi-spherical approximation is presently used under the assumption of hydrostatic equilibrium to model giant-planet formation processes that encompass the accretion and transport of angular momentum for both radiative and convective zones. These calculations have been conducted up to masses comparable to those of Saturn. Angular momentum transport is strong, leading to strong concentration of angular momentum in the outer convective zones which develop envelope masses greater than 30 earth masses. When the models are allowed to contract near the calculation's end, the ratio of centrifugal force to gravity at the outer radius rises sufficiently to validate the quasi-spherical approximation.

This study was initially conceived in order to aid in the interpretation of Martian surface and atmospheric aerosol mineralogy. As a result, the minerals included are biased toward samples which represent hydrated and hydroxylated silicates. Due to their physical particle size, clays and other materials, such as palagonite, cannot be prepared using typical preparation techniques. Yet in some cases, such as for Mars, these are the materials of perhaps the greatest interest. In order to obtain a suitable sample of these less cohesive materials for the laboratory measurements, a KBr pellet die was used and a pellet of the pure sample was prepared. For all clays and the palagonite, a pellet with highly reflective surfaces at visible wavelengths was produced. The reflectivities of all samples were determined, and to derive the optical constants of a materials as a function of wavelength, dispersion analysis was used which describes the real and imaginary indices of refraction as the contributions due to a sum of classical oscillators and relates them via Frensel's equations for non-normal incidence, the the measured near-normal reflectivity. The final values that were determined represent averages of several model fits to each data set using the same number of oscillators by varying the high frequency dielectric constant.

A baseline model of the composition and abundances of grains and gases in molecular cloud cores and accretion disks around young stars is defined by employing: a wide range of astronomical data and theory; the composition of primitive bodies in the solar system; and solar elemental abundances. It is proposed that in the coldest portions of these objects the major grain species include amorphous olivine, amorphous orthopyroxene, volatile and refractory organics, water ice, troilite, and metallic iron. Using a combination of laboratory measurements of optical constants and asymptotic theory, values of the real and imaginary indices of refraction of these grain species over a wavelength range that runs from the vacuum UV to the radio domain are derived. Auxiliary information for these grain species, such as their vaporization temperatures bulk densities, and their fractional abundances by mass. The above information on grain properties in molecular cloud cores and accretion disks is used to estimate the Rosseland mean opacity of the grains in both environments and the IR and microwave opacity of grains in accretion disks.

Thirteen spectra of Mars, in the 7.5 to 12.8 micron wavelength were obtained on 7 Dec. 1990 from the Infrared Telescope Facility (IRTF). For these observations, a grating with an ultimate resolving power of 120 to 250 was used and wavelengths were calibrated for each grating setting by comparison with the absorption spectrum of polystyrene measured prior to each set of observations. By sampling the Nyquist limit at the shortest wavelengths, an effective resolving power of about 120 over the entire wavelength range was achieved. A total of four grating settings were required to cover the entire wavelength region. A typical observing sequence consisted of: (1) positioning the grating in one of the intervals; (2) calibrating the wavelength of positions; and (3) obtaining spectra for a number of spots on Mars. Several observations of the nearby stellar standard star, alpha Tauri, were also acquired throughout the night. Each Mars spectrum represents an average of 4 to 6 measurements of the individual Mars spots. As a result of this observing sequence, the viewing geometry for a given location or spot on Mars does not change, but the actual location of the spot on Mars's surface varies somewhat between the different grating settings. Other aspects of the study are presented.

Recent reports concerning the mid-infrared reflectance properties of silicates coupled with recent observations of the Earth and other planets in the mid-infrared and the planned Thermal Emission Spectrometer scheduled as an instrument to be included on the Mars Observer all illustrate the increasing interest in the optical properties of materials in the mid-infrared and their direct application to remote sensing observations of other planetary surfaces. As the laboratory and observational data increase they will ultimately be modeled to aid in the understanding of the composition mineralogy, and distribution of the surface and atmospheric constituents on these bodies. In order to facilitate such quantitative analyses, knowledge regarding the optical constants (real (n) and imaginary (k) indices of refraction) of a wide variety of pertinent materials is required. Examples of the application of such quantitative analyses to the interpretation of Martian surface and atmospheric constituents, based on the optical constants of minerals, are presented. This study was initially conceived in order to aid in the interpretation of Martian surface and atmospheric aerosol mineralogy. As a result, the minerals included are biased toward samples which represent hydrated and hydroxylated silicates. These include the Al and Mg end members of the 1:1 layer lattice silicates, kaolinite and serpentine, respectively; an Al-bearing 2:1 layer lattice silicate, pyrophyllite; the Mg and Al smectite clays saponite and montmorillonite, respectively; and a palagonite, typically a poorly characterized alteration product of basaltic glass. The reflectivities of all samples were determined by placing them at the focus of a near-normal reflective attachment located in a Fourier transform spectrometer. To derive the optical constants of a material as a function of wavelength, researchers used the commonly employed technique of dispersion analysis.

The Martian southern hemisphere atmospheric water vapor column abundance measurements reported agree with Viking Orbiter atmospheric water detectors during early southern spring and southern autumnal equinox; profiles obtained in southern mid- and late summer, however, indicate the presence of twice as much water both in the southern hemisphere and planetwide. This discrepancy is accounted for by the high optical depths created by two global dust storms during the Viking year, while the present observations were obtained in the case of the relatively dust-free atmosphere of the 1988-1989 opposition.

The basic physics of the environmental perturbations caused by multiple nuclear detonations is explored, summarizing current knowledge of the possible physical, chemical, and biological impacts of nuclear war. Emphasis is given to the impact of the bomb-generated smoke (soot) particles. General classes of models that have been used to simulate nuclear winter are examined, using specific models as examples.

The midinfrared 2000-400/ cm (5-25 microns) optical constants (real (n) and imaginary (k) indices of refraction) are presented for: (1) pyrophyllite; (2) kaolinite; (3) serpentine; (4) montmorillonite; (5) saponite; (6) palagonite; and (7) orthopyroxene. Comparison of the values derived here with those previously presented for serpentine, montmorillonite, and palagonite is generally quite good and discrepancies between values are probably due to either chemical differences between the actual samples or different techniques used to derive the values. For montmorillonite, saponite, and palagonite we were able to derive optical constants in the region of the H2O-bending fundamental near 6 microns. We find that if a pellet of pure material can be produced with a mirror-like surface then the optical constants of clays and other noncohesive materials can be readily derived.

A program, called the Cassini mission, intended to study the Saturn system by utilizing a Saturn orbiter and a probe descending to the surface of Titan, is discussed. Winds are expected to cause perturbations to the probe local horizontal velocity, resulting in an anomalous drift in the probe location and a shift in the frequency of the probe telemetry, due to the Doppler effect. By using an iterative algorithm, in which the time variation of the probe telemetry frequency is monitored throughout the descent, and the probe trajectory is updated to reflect the effect of wind on the probe location, a highly accurate relative wind profile can be recovered. By adding a single wind velocity, measured by independent means, an absolute wind profile can be obtained. However, the accuracy of the zonal winds recovery is limited by errors in trajectory, and frequency.

The results of a comparison of laboratory spectral reflectance measurements of particulate mixtures of both hydrated silicates and palagonite with water ice, on the one hand, with two previously unpublished reflectance spectra of Callisto, yield direct support for the hypothesis that the measured reflectance of Callisto includes a substantial nonice component. Hapke's (1981) equations are used in a theoretical model for thorough analysis of telescopic data; a comparison of the calculation results thus obtained with measured reflectance data for Callisto indicate that three-component, ice/magnetite/serpentine mixtures are a better match for telescopic data than two-component ice mixtures.

Methane is the most abundant simple organic molecule in the outer solar system bodies. In addition to being a gaseous constituent of the atmospheres of the Jovian planets and Titan, it is present in the solid form as a constituent of icy surfaces such as those of Triton and Pluto, and as cloud condensate in the atmospheres of Titan, Uranus, and Neptune. It is expected in the liquid form as a constituent of the ocean of Titan. Cometary ices also contain solid methane. The optical constants for both solid and liquid phases of CH4 for a wide temperature range are needed for radiative transfer calculations, for studies of reflection from surfaces, and for modeling of emission in the far infrared and microwave regions. The astronomically important visual to near infrared measurements of solid methane optical constants are conspicuously absent from the literature. Preliminary results are presented on the optical constants of solid methane for the 0.4 to 2.6 micrometer region. Deposition onto a substrate at 10 K produces glassy (semi-amorphous) material. Annealing this material at approximately 33 K for approximately 1 hour results in a crystalline material as seen by sharper, more structured bands and negligible background extinction due to scattering. The constant k is reported for both the amorphous and the crystalline (annealed) states. Typical values (at absorption maxima) are in the .001 to .0001 range. Below lambda = 1.1 micrometers the bands are too weak to be detected by transmission through the films less than or equal to 215 micrometers in thickness, employed in the studies to date. Using previously measured values of the real part of the refractive index, n, of liquid methane at 110 K, n is computed for solid methane using the Lorentz-Lorenz relationship. Work is in progress to extend the measurements of optical constants n and k for liquid and solid to both shorter and longer wavelengths, eventually providing a complete optical constants database for condensed CH4.

One-dimensional, spherically symmetric hydrodynamic calculations of the effect of the impact of a 1-2 earth mass body on the gaseous envelope of Uranus, characterizing the impact as a partition of energy such that a fraction is deposited directly into the envelope, while the rest heats the core. A sharp transition is found between the cases of nearly-complete retention and dispersal, corresponding to energy-deposition values above or below a critical value. Within a realistic parameter range, the angular momentum of the Uranus system may be explained, together with the retention of a gaseous envelope of about 1 earth mass.

The satellites of Uranus, with densities between 1.3 and 1.7 g cm(-3) (from Voyager 2 observations) and the Pluto-Charon system, with a mean density of just above 1.8 g cm(-3) (from terrestrial observations of mutual eclipse events), are too dense to have a significant amount of methane ice in their interiors. However, the observed densities do not preclude contributions from such organic materials as the acid-insoluble residue in carbonaceous chondrites and laboratory-produced tholins, which have densities on the order of approximately 1.5 g cm(-3). These and other considerations have led researchers to investigate the carbon mass budget in the outer solar system, with an emphasis on understanding the contribution of organic materials. Modeling of the interiors of Pluto and Charon (being carried out by R. Reynolds and A. Summers of NASA/Ames), assuming rock and water ice as the only constituents, suggests a silicate mass fraction for this system on the order of 0.65 to 0.70. The present work includes the most recent estimates of the C/H enhancements and high z/low z ratios of the giant planets (Pollack and Bodenheimer, 1987), and involves a more careful estimation of the high z/low z mass ratio expected from solar abundances than was used in Pollack et al. (1986), including the influence of the fraction of C in CO on the amount of condensed water ice. These calculations indicate that for a particular fraction of C in CO and a given fraction of C-bearing planetesimals that dissolve in the envelope (most likely in the range 0.50 to 0.75), (1) Jupiter and Saturn require a larger fraction of C in condensed materials than Uranus and Neptune, but (2) the Jupiter and Saturn results are much less strongly constrained by the error bars on the observed C/H enhancements and high z/low z ratios than is the case for Uranus and Neptune. The clearest result is that in the region of the solar nebula near Uranus and Neptune, the minority of carbon that is not in gaseous CO (1) must include a nonzero amount of condensed material, but (2) is most likely not condensed material alone, i.e., there must be a third carbon-bearing component besides condensed material and gaseous CO. Given the implied dearth of methane ice, the condensed carbon is likely dominated by organic material, and the third component present in addition to CO and organics is assumed to be CH4 gas.

Differentially rotating disks of gases and solids occur in several astrophysical systems, in particular in the inner parts of protostellar nebulae, of which our own solar system is thought to be a relic. The objectives of this paper are to: (1) study localized turbulence in circumstances approximating those found in accretion disks using previously existing expertise in performing direct numerical simulation of turbulent, incompressible channel flows with low Reynolds numbers; (2) determine the limitations of such calculations; and (3) extend the type of numerical simulation (e.g., to include density and stratification and compressibility effects and to accommodate higher Reynolds numbers with sub-grid scale modeling) so that the relevant physical effects are realistically captured.

A reevaluation is presented of the 'nuclear winter' scenario of Turco et al. (1983). New pertinent data have emerged from laboratory studies, field experiments, and numerical models on the smoke-plume, mesoscale, and global scales. A full-scale nuclear exchange's probable soot injections lead, in three-dimensional climate simulations, to midsummer land temperature decreases averaging 10-20 C in northern midlatitudes, with local cooling of as much as 35 C. Anomalous circulation patterns due to solar heating of the soot could stabilize the upper atmosphere against overturning, thereby prolonging the soot's residence time in the atmosphere. Monsoon disruptions and severe ozone layer depletion are also foreseen.

Turbulent convection may have played a major role in determining the structure and evolution of the primordial solar nebula, but current, incomplete models of convection and turbulence give very different results and remain largely untested in the absence of detailed astronomical observations. Numerical simulations provide an 'experimental' database for comparison with these models, and, to this end, direct numerical simulations of turbulent convection were performed with modifications intended to mimic some of the unique physical features of thin accretion disks, such as the primordial solar nebula: internal heating, a gravitational acceleration that is linearly proportional to the distance from midplane of the nebula, and rapid rotation. Peclet numbers in the simulations are comparable to those in solar nebula models; Rossby numbers in the simulations are an order of magnitude larger than those in solar nebula models because of the unrealistically high Prandtl and low, Reynolds numbers required to resolve all scales of the convective flow. We find that, despite the loss of buoyancy at midplane, turbulent motions easily penetrate the midplane region with little loss of intensity, providing efficient transport of heat and turbulent kinetic energy throughout the interior. A simple mixing length model modified to include rotation is found to give convective heat fluxes for the interior flow in rough agreement with the numerical simulations. We discuss these preliminary results with regard to assumptions about heating distributions and convective heat fluxes made in standard solar nebula models. More definitive comparisons with solar nebula modelling will become possible when more realistic effects of shear, density, stratification, and compression are included.

Images of Neptune obtained by the narrow-angle camera of the Voyager 2 spacecraft reveal large-scale cloud features that persist for several months or longer. The features' periods of rotation about the planetary axis range from 15.8 to 18.4 hours. The atmosphere equatorward of -53 deg rotates with periods longer than the 16.05-hour period deduced from Voyager's planetary radio astronomy experiment (presumably the planet's internal rotation period). The wind speeds computed with respect to this radio period range from 20 meters per second eastward to 325 meters per second westward. Thus, the cloud-top wind speeds are roughly the same for all the planets ranging from Venus to Neptune, even though the solar energy inputs to the atmospheres vary by a factor of 1000.