ALBERT

Description: NASA/GRO grant NAG 5-2081, at the University of Chicago, has provided support for a broad program of theoretical research in nuclear astrophysics and related areas, with regard to gamma-ray and hard X-ray emission from classical nova explosions. This research emphasized the possible detection of 22Na gamma-ray line emission from nearby novae involving ONeMg white dwarfs, the detailed examination of 26Al production in novae, and the possible detection of the predicted early gamma ray emission from novae that arises from the decay of the short lived, positron emitting isotopes of CNO elements. Studies of nova related problems have consumed an increasing fraction of the Principal Investigator's research efforts over the past decade. Current research addresses problems associated with the standard model for the outbursts of the classical novae: the occurrence of thermonuclear runaways (TNR) in the accreted hydrogen rich envelopes on white dwarfs in close binary systems (see, e.g., the reviews by Truran 1982; and Shara 1989). Research in progress and planned for the next three years has three main objectives: (1) to gain an improved understanding of the early evolution of the light curves of, particularly, the fastest novae; (2) to gain an improved understanding of the relative importance of the various possible mechanisms of envelope hydrogen depletion (e.g. winds, common envelope driven mass loss, and nuclear burning) to the long term evolution of novae in outburst; and (3) to seek to provide a somewhat more definitive statement of the role of classical novae in nucleosynthesis. Our proposed 2-D studies of convection during the early phases of the TNR and our systematic attempt to incorporate an improved treatment of radiation hydrodynamics into the hydrodynamic code utilized in our calculations, are particularly relevant to the first of these objectives. Further 2-D studies of the effects of common envelope evolution are intended to provide more realistic constraints on the mass depletion mechanisms. Finally, detailed calculations of the thermonuclear history of the matter ejected in novae will be carried out for representative nova configurations involving both carbon-oxygen (CO) and oxygen-neon-magnesium (ONeMg) white dwarfs.

Description: The funds from this grant were used to support observations and analysis with the International Ultraviolet Explorer (IUE) satellite telescope. The main area of scientific research concerned the variability analyses of ultraviolet spectra of Active Galactic Nuclei, primarily quasars, Seyfert galaxies, and BL Lacertae objects. The Colorado group included, at various times, the P.I. (J.M. Shull), Research Associate Dr. Rick Edelson, and graduate students Jon Saken, Elise Sachs, and Steve Penton. A portion of the work was also performed by CU undergraduate student Cheong-ming Fu. A major product of the effort was a database of all IUE spectra of active galactic nuclei. This database is being analyzed to obtain spectral indices, line fluxes, and continuum fluxes for over 500 AGN. As a by-product of this project, we implemented a new, improved technique of spectral extraction of IUE spectra, which has been used in several AGN-WATCH campaigns (on the Seyfert galaxy NGC 4151 and on the BL Lac object PKS 2155-304).

Description: The progress made during 1995 on the Monte-Carlo gamma-ray spectrum simulation program BSIMUL is discussed. Several features have been added, including the ability to model shields that are tapered cylinders. Several simulations were made on the Near Earth Asteroid Rendezvous detector.

Description: This paper reports a new balance for the measurement of three components of force - lift, drag and pitching moment - in impulsively starting flows which have a duration of about one millisecond. The basics of the design of the balance are presented and results of tests on a 15 deg semi-angle cone set at incidence in the T4 shock tunnel are compared with predictions. These results indicate that the prototype balance performs well for a 1.9 kg, 220 mm long model. Also presented are results from initial bench tests of another application of the deconvolution force balance to the measurement of thrust produced by a 2D scramjet nozzle.

Description: The second guide star catalog (GSC-2) is a project that aims to create a complete catalog of stars and galaxies to about the 18th magnitude, and would contain colors, magnitudes, positions and proper motions. The catalog would provide an object list for the construction of an input catalog for possible future astrometric satellites, such as the global astrometric interferometer for astrophysics (GAIA) satellite. With a schedule that is compatible with the projected timeframe of GAIA, the GSC-2 could be available in time to support the preparatory astrophysical observations for an input catalog. The availability of the digitalized images for the preparation of finding charts would improve the efficiency of this task.

Description: We present radio observations of the gravitational lens PKS 1830-211 at 8.4 and 15 GHz acquired using the Very Large Array. The observations were made over a 13 month period. Significant flux density changes over this period provide strong constraints on the time delay between the two lensed images and suffest a value of 44 +/- 9 days. This offers new direct evidence that this source is indeed a gravitational lens. The lens distance is dependent upon the model chosen, but reasonable limits on the mass of the lensing galaxy suggest that it is unlikely to be at a redshift less than a few tenths, and may well be significantly more distant.

Description: The power laws are approximately f(exp -1.9), f(exp -1.9), and f(exp -2.1) respectively for the Grigg-Shjellerup (GS), Giacobini-Zinner (GZ), and Halley (H) comets. Other than similarities in the power spectra, the magnetic field turbulence is considerably different at the three comets. Phase steepening is demonstrated to occur at the trailing edges of the GS waves. This is probably due to nonlinear steepening plus dispersion of the left-hand mode components, i.e., the turbulence is whistler-mode. This too can be explained by nonlinear steepening plus dispersion of the magnetosonic waves. At the level of GS and GZ turbulence development when the spacecraft measurements were made, classical three-wave processes, such as the decay or modulation instabilities do not appear to play important roles. It is most likely that the nonlinear steepening and dispersive time scales are more rapid than three-wave processes, and the latter had not had time to develop for the relatively new turbulence. The wave turbulence at Halley is linearly polarized. The exact nature of this turbulence is still not well understood. Several possibilities are suggested, based on a preliminary analyses.

Description: A preliminary comparison of the GEOS-1 (Goddard Earth Observing System) data assimilation system convective cloud mass fluxes with fluxes from a cloud-resolving model (the Goddard Cumulus Ensemble Model, GCE) is reported. A squall line case study (10-11 June 1985 Oklahoma PRESTORM episode) is the basis of the comparison. Regional (central U. S.) monthly total convective mass flux for June 1985 from GEOS-1 compares favorably with estimates from a statistical/dynamical approach using GCE simulations and satellite-derived cloud observations. The GEOS-1 convective mass fluxes produce reasonable estimates of monthly-averaged regional convective venting of CO from the boundary layer at least in an urban-influenced continental region, suggesting that they can be used in tracer transport simulations.

Description: A study of X-ray emission from five short-period Algol-type binaries based on observations with Advanced Satellite for Cosmology and Astrophysics (ASCA) and ROSAT is presented. We have observed RZ Cas with both satellites, and beta Per, U Cep, delta Lib, and TW Dra with ROSAT. Significant intensity variations are seen in the X-ray emission from RZ Cas, U Cep, TW Dra, and delta Lib. These variations seem unrelated to the eclipsing behavior of these systems and are probably due to either rotational modulation of compact active regions on the surfaces of the chromospherically active secondary components or to flaring activity in the systems. The spectra of all but one of the systems require the presence of at least two discrete plasma components with different temperatures (0.6 - 0.7 keV, and approximately 2 keV) and the abundances of the medium-Z elements 20% - 50% of the solar photospheric values. The high resolving power and signal-to-noise ratio of the ASCA spectra allow us to individually constrain the coronal abundances of O, Ne, Mg, Si, S, and Fe in RZ Cas. We demonstrate that, if we use the elemental abundances and temperatures obtained from the analysis of their ASCA spectra as (fixed) inputs, to fit the ROSAT PSPC spectra well requires the presence of a third component (kT approximately 0.2 - 0.3 keV) in RZ Cas and beta Per. A continuous emission measure model of the power-law type (EM(T) variesas (T/T(sub max)(sup alpha)) generally gives a poor fit to the ASCA and ROSAT data on most sources. Circumstellar or circumbinary absorbing matter seems to be present in some of these systems, as indicated by the variable total column density needed to fit their X-ray spectra.

Description: The U.S. upper Midwest was subjected to severe flooding during the summer of 1993. Heavy rainfall in the Mississippi River basin from April through July caused flooding on many Midwest rivers, including the Mississippi, Illinois, Missouri, and Kansas Rivers. The flood crest of 15.1 m at St. Louis, Missouri, on 1 August 1993 was the highest ever measured, surpassing the previous record of 13.2 m set on 28 April 1973. Damage estimates include at least 47 flood-related deaths and a total damage cost of $12 billion. Remotely sensed imagery of severe flooding in the U.S. Midwest was obtained under cloud-free skies on 29 July 1993 by the MODIS (Moderate Resolution Imaging Spectroradiometer) Airborne Simulator (MAS). The MAS is a newly developed scanning spectrometer with 50 spectral bands in the wavelength range 0.55-14.3 micrometers. Estimation of the total flooded area in the MAS scene acquired near St. Louis was accomplished by comparing the MAS scene to a Landsat-5 thematic mapper (TM) scene of the same area acquired on 14 April 1984 in nonflood conditions. For comparison, the MAS band centered at 0.94 micrometers and the TM band centered at 1.65 micrometers were selected because of the high contrast seen in these bands between land and water-covered surfaces. An estimate of the area covered by water in the MAS and TM scenes was obtained by developing land/water brightness thresholds from histograms of the MAS and TM digital image data. Afetr applying the thresholds, the difference between the area covered by water in the MAS and TM scenes, and hence the flooded area in the MAS scene, was found to be about 396 sq km, or about 153 square miles.

Description: Shortwave radiative fluxes that reach the earth's surface are key factors that influence atmospheric and oceanic circulations as well as surface climate. Yet, information on these fluxes is meager. Surface site data are generally available from only a limited number of observing stations over land. Much less is known about the large-scale variability of the shortwave radiative fluxes over the oceans, which cover most of the globe. Recognizing the need to produce global-scale fields of such fluxes for use in climate research, the World Climate Research Program has initiated activities that led to the establishment of the Surface Radiation Budget Climatology Project with the ultimate goal to determine various components of the surface radiation budget from satellite data. In this paper, the first global products that resulted from this activity are described. Monthly and daily data on a 280-km grid scale are available. Samples of climate parameters obtainable from the dataset are presented. Emphasis is given to validation and limitations of the results. For most of the globe, satellite estimates have bias values between +/- 20 W/sq m and root mean square (rms) values are around 25 W/sq m. There are specific regions with much larger uncertainties however.

Description: The deployment of a space-based Doppler lidar would provide information that is fundamental to advancing the understanding and prediction of weather and climate. This paper reviews the concepts of wind measurement by Doppler lidar, highlights the results of some observing system simulation experiments with lidar winds, and discusses the important advances in earth system science anticipated with lidar winds. Observing system simulation experiments, conducted using two different general circulation models, have shown (1) that there is a significant improvement in the forecast accuracy over the Southern Hemisphere and tropical oceans resulting from the assimilation of simulated satellite wind data, and (2) that wind data are significantly more effective than temperature or moisture data in controlling analysis error. Because accurate wind observations are currently almost entirely unavailable for the vast majority of tropical cyclones worldwide, lidar winds have the potential to substan- tially improve tropical cyclone forecasts. Similarly, to improve water vapor flux divergence calculations, a direct measure of the ageostrophic wind is needed since the present level of uncer- tainty cannot be reduced with better temperature and moisture soundings alone.

Description: The 183-GHz water vapor line was tentatively detected on Mars in January 1991, with the IRAM 30-m millimeter antenna, under extremely dry atmospheric conditions. The measurement refers to the whole disk. The spectral line, although marginally detected, can be fit with a constant H2O mixing ratio of 1.0 x 10(exp -5), which corresponds to a water abundance of 1 pr-microns; in any case, an upper limit of 3 pr-microns is inferred. This value is comparable to the very small abundances measured by Clancy (1992) 5 weeks before our observation and seems to imply both seasonal and long-term variations in the martian water cycle.

Description: Titan's large free eccentricity results in significant tidal dissipation. This can be used to constrain the existence and depth of hydrocarbon oceans. A hydrodynamical-numerical 2 deg ocean tide model has been constructed to investigate this connection. The model allows some simple land configurations. The results indicate that existence of such an ocean over the age of the solar system is hard to explain, as is the existence of the large eccentricity itself. If such an ocean exists, it is likely to be more than 500 m deep, ignoring the influence of land masses.

Description: The 'satellite-gauge model' (SGM) technique is described for combining precipitation estimates from microwave satellite data, infrared satellite data, rain gauge analyses, and numerical weather prediction models into improved estimates of global precipitation. Throughout, monthly estimates on a 2.5 degrees x 2.5 degrees lat-long grid are employed. First, a multisatellite product is developed using a combination of low-orbit microwave and geosynchronous-orbit infrared data in the latitude range 40 degrees N - 40 degrees S (the adjusted geosynchronous precipitation index) and low-orbit microwave data alone at higher latitudes. Then the rain gauge analysis is brougth in, weighting each field by its inverse relative error variance to produce a nearly global, observationally based precipitation estimate. To produce a complete global estimate, the numerical model results are used to fill data voids in the combined satellite-gauge estimate. Our sequential approach to combining estimates allows a user to select the multisatellite estimate, the satellite-gauge estimate, or the full SGM estimate (observationally based estimates plus the model information). The primary limitation in the method is imperfections in the estimation of relative error for the individual fields. The SGM results for one year of data (July 1987 to June 1988) show important differences from the individual estimates, including model estimates as well as climatological estimates. In general, the SGM results are drier in the subtropics than the model and climatological results, reflecting the relatively dry microwave estimates that dominate the SGM in oceanic regions.

Description: A number of planetary objects exhibit large radar reflectivity and polarization ratios, and more recently, a similar behavior has been observed over a vast portion of the Earth's surface: the percolation facies of the Greenland Ice Sheet. Surface-based ranging radar data and snow stratigraphy studies demonstrated that the radar properties of that portion of Greenland are caused by enhanced scattering from massive, large, solid-ice bodies buried in the top few meters of the dry, cold, clean snowy surface of the ice sheet and created by seasonal melting and refreezing events. Here, we model the icy inclusions as randomly oriented, discrete, noninteracting, dielectric cylinders embedded in a transparent snow medium. An exact analytical solution is used to compute the scattered field from the cylinders. Using this model, we correctly predict the polarimetric radar observations gathered by an airborne imaging system at three wavelengths (5.6, 24, and 68 cm), between 19 deg and 65 deg incidence angle. The diameter and number density of the cylinders that are inferred from the radar data using the model are consistent with in situ observations of the icy inclusions. The large radar reflectivity and polarization ratios are interpreted as arising internal reflections of the radar signals in the icy inclusions that first-order external reflection models fail to predict. The results compare favorably with predictions from the coherent backscatter or weak localization theory and may provide a complementary framework for interpreting exotic radar echoes from other planetary objects.

Description: 24 years of Lunar Laser Ranging (LLR) observations and 16 years of Very Long Baseline Interferometry (VLBI) observations are combined in a global analysis to yield improved estimates of the Earth's precession and nutation. The correction to the International Astronomical Union (IAU) (1976) precession constant inferred from this joint VLBI/LLR analysis is -3.00 +/- 0.20 milliarcsec/yr (mas/yr). A significant obliquity rate correction of -0.20 +/- 0.08 mas/yr is also found. In all, 32 forced nutation coefficients are estimated. These coefficients confirm that the IAU (1980) nutation theory is in error by several mas. The estimated nutation coeficients are found to vary by as much as several tenths of mas, depending on the a priori nutation model used to analyze the VLBI and LLR data. Forced circular nutations derived from this analysis agree with the ZMOA-1990-2 nutation theory at the 0.2 mas level for the 18.6 yr terms, and at the 0.05 mas level for the other terms (periods less than or = 1 yr). A retrograde free core nutation with an amplitude of 0.20 mas is also detected. Its phase is found to be very sensitive to the precise value of the free core nutation period used in the solution. Separate analyses of four independent subsets of the LVBI data indicate no significant variations of the free core nutation since 1988. The pre-1988 estimates of the free core nutation are consistent with the post-1988 estimates but are not accurate enough to rule out possible variations of the free core nutation at these earlier epochs.

Description: The Ulysses spacecraft has gathered data from within flows from the Sun's southern polar coronal hole, the first in situ measurement of this region. We present a brief analysis of the heliospheric magnetic field data from this region, using a fractal method. As is the case near the ecliptic, estimated spectral exponents are near 5/3 on spacecraft scales of seconds to minutes. On longer time scales, however, there appears to be a significantly different population in polar flows, which is similar to that found by the Helios spacecraft in fast solar wind flows at 0.3 AU.

Description: The lower stratospheric variability of equatorial water vapor, measured by the Microwave Limb Sounder (MLS), follows an annual cycle modulated by the quasi-biennial oscillation. At levels higher in the stratosphere, water vapor measurements exhibit a semi-annual oscillatory signal with the largest amplitudes at 2.2 and 1hPa. Zonal-mean cross sections of MLS water vapor are consistent with previous satellite measurements from the limb infrared monitor of the stratosphere (LIMS) and the stratospheric Aerosol and Gas Experiment 2 (SAGE 2) instruments in that they show water vapor increasing upwards and the polewards from a well defined minimum in the tropics. The minimum values vary in height between the retrieved 46 and 22hPa pressure levels.

Description: A calculation, employing a detailed model of neutral oxygen, is carried out to give fluorescent line intensities expected in a long-proposed photoexcitation by accidental resonance (PAR) process in which hydrogen Lyman-beta photoexcites the oxygen spectrum. The results pertain to the optically thin case but provide an upper limit to the fluorescent intensities which can be attained. They are applied to analyze line ratios involving the strong 8446 A line observed in classical novae during the diffusion-enhanced and Orion phases. Operation of the PAR process in the novae is verified. It is found that photoexcitation rates in the ejecta reach values greater than 0.1/sec, corresponding to hydrogen Lyman-beta radiation field intensities greater than 1250 ergs/cm/sec/sr.

Description: We have investigated the properties of the OMC-2 and OMC-3 cores in the Orion giant molecular cloud using high spatial spectral resolution observations of several transitions of the (13)CO, C(18)O, C(S-32) and C(S-34) molecules taken with the SEST telescope. The OMC-2 core consists of one clump (22 solar mass) with a radius of 0.11 pc surrounded by a cluster of 11 discrete infrared sources. The H2 column density and volume density in the center of this clump are 2 x 10(exp 22)/sq cm and 9 x 10(exp 5)/cu cm respectively. From a comparison between physical parameters derived from C(18)O and C(S-32) observations we conclude that the molecular envelope around the core has been completely removed by these sources and that only the very dense gas is left. OMC-3 shows a more complex elongated structure in C(18)O and CS than OMC-2. The C(S-32) and C(S-34) maps show that the denser region can be separated into at least sub-cores of roughly equal sizes (radius approximately equals 0.13 pc), with n(H2) = 6 x 10(exp 5)/cu cm, and a mass of 10 solar mass (from C(S-32)). The very different masses obtained for the central core from C(18)O and C(S-32) (55 and 12 solar mass respectively) indicate that a massive envelope is still present around the very dense sub-cores. We report the first detection of several molecular outflows in OMC-3. The presence of an IRAS source and the first detection of these outflows confirm that star formation is going on in OMC-3. Based on the different physical properties of these regions compared with OMC-1, OMC-2 appears to be in an intermediate evolutionary stage between OMC-1 and OMC-3.

Description: The depletion of condensable elements onto grains in gaseous nebulae can provide evidence that dust is well mixed with the ionized gas. Al and CA are two of the most depleted elements in the general interstellar medium, and it is therefore important to measure their abundances within the ionized region of nebulae. We compute a large grid of photoionization models and identify sets of line ratios which are relatively insensitive to stellar and nebular parameters, and are thus excellent diagnostics for determining relative abundances. Based on the absence of the ((Ca II) lambda lambda 291, 7324 doublet and the detection of Al II) lambda lambda 2660, 2669 in the ultraviolet, we determine the extent of aluminum and calcium depletion onto grains in NGC 7027 and the Orion Nebula. Our results show a approximately 0.3 dex depletion for Al, but a depletion of more than two and a half orders of magnitude for Ca. A similar calculation based on Mg II lambda 2798 yields roughly a 0.8 dex depletion for Mg. This reaffirms the discrepancy between depletion determined from high and low ionization Mg lines. We also find evidence for a 'depletion gradient' in Ca in NGC 7027, since the calcium depletion we infer for the outer, more neutral regions using (Ca II) is somewhat higher than that inferred for the inner high-ionization region, using (Ca v). This gradient can test current models of the survival of grains within hot ionized gas.

Description: The absolute rate constant for the reaction O((3)P) + HOBr has been measured between T = 233K and 423K using the discharge-flow kinetic technique coupled to mass spectrometric detection. The value of the rate coefficient at room temperature is (2.5 +/- 0.6) x 10(exp -11)cu cm/molecule/s and the derived Arrhenius expression is (1.4 +/- 0.5) x 10(exp -10) exp((-430 +/- 260)/T)cu cm/molecule/s. From these rate data the atmospheric lifetime of HOBr with respect to reaction with O((3)P) is about 0.6h at z = 25 km which is comparable to the photolysis lifetime based on recent measurements of the UV cross section for HOBr. Implications for HOBr loss in the stratosphere have been tested using a 1D photochemical box model. With the inclusion of the rate parameters and products for the O + HOBr reaction, calculated concentration profiles of BrO increase by up to 33% around z = 35 km. This result indicates that the inclusion of the O + HOBr reaction in global atmospheric chemistry models may have an impact on bromine partitioning in the middle atmosphere.

Description: Synthetic-aperture radar (SAR) interferograms produced from ESA's ERS-1 satellite, provide the first synoptic view of ice flow dynamics of the western sector of the Greenland Ice Sheet. Glacial motion is detected in the radar ranging direction at millimetric scales, across a complete sequence of snow accumulation and melting regimes, despite significant varations in their radar scattering properties. Ice flow evolves from a slow, regular motion at the higher elevations. At lower elevations, motion is strongly convoluted by meter-scale undulations in surface topography, which have a unique interferometric signature that enables a novel approach for retrieving flow direction. Inferred flow directions, combined with surface displacements in the radar ranging direction, yield ice velocity estimates that are within 6% of in-situ measurements gathered along a 40 km survey line. Application of repeat-pass SAR interferometry to the entire Greenland Ice Sheet should enable precise mapping of its ice flow dynamics at an unprecedented level of spatial detail.

Description: The modern scenario of evolution of massive binary systems predicts the existence of a subclass of binary radio pulsars (PSRs) with black holes (BHs). Their Galactic number was evaluated as approximately 1 per 1000 single pulsars (Lipunov et al. 1994b). Distinctive properties of such binaries would be (1) mass of the unseen companion M(sub c) greater than 3-4 solar mass and (2) absence of eclipses of the pulsar radiation with no distinctive variance of the dispersion measure along the pulsar orbit. The pulsars themselves must be similar to standard isolated ones. The recently discovered binary 1 s pulsar PSR B0042-73 = PSR J0045-7319 in the Small Magellanic Cloud (SMC) with a massive companion in a highly elongated (eccentricity e =0.8) 51 day orbit (Kaspi et al. 1994) may be the first such pulsar with a BH. The paradoxical fact that the first pulsar discovered in the SMC proved to be in a binary system can be naturally understood if its companion actually is a 10-30 solar mass black hole. We illustrate this fact by the numerical calculation of evolution of radio pulsars after a star formation burst.

Description: We calculate the structure of an effectively optically thin and geometrically thin accretion disk in the Kerr geometry, including electron-positron pairs. We show that the properties of the disk solutions are strongly dependent on the angular momentum of the central black hole. We find that close to a rapidly rotating hole there can be an appreciable pair density even for modest accretion rates. Pair critical accretion rates recently discovered in Newtonian disk models are also shown to be present in the general relativistic models, and we show that the geometrically thin disk approximations easily break down for rapidly rotating holes.

Description: We examine the hydrodynamic origin of relativistic outflows in active galactic nuclei (AGN). Specifically, we propose that the presence of a population of relativistic hadrons in the AGN 'central engine' and the associated neutron production suffices to produce outflows which under rather general conditions could be relativistic. The main such condition is that the size of the neutron production region be larger than the neutron flight path tau(sub n) approximately 3 x 10(exp 13) cm. This condition guarantees that the mean energy per particle in the proton fluid, resulting from the decay of the neutrons outside their production region, be greater than the proton rest mass. The expansion of this fluid can then lead naturally to a relativistic outflow by conversion of its internal energy to directed motion. We follow the development of such flows by solving the mass, energy as well as the kinetic equation for the proton gas in steady state, taking into account the source terms due to compute accurately the adiabatic index of the expanding gas, and in conjunction with Bernoulli's equation the detailed evolution of the bulk Lorentz factor. We further examine the role of large-scale magnetic fields in confining these outflows to produce the jets observed at larger scales.

Description: A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Description: Analysis of the version 16 Halogen Occultation Experiment (HALOE) CH4 data shows that this long-lived trace gas is well correlated with potential vorticity (PV) computed from National Meteorological Center balanced winds. Analyzing late September and October 1992 data, we show that very low CH4 values are confined to the interior of a vortex edge defined by the maximum gradient in PV. The CH4 and HF time tendency is used to estimate the descent rate in the Antarctic vortex. After removing a component of the trend correlated with the HALOE sampling pattern, we compute the lower stratosphere vertical descent rates and net heaing rates in the spring Antarctic vortex. Our computations of the spring Antarctic vortex heating rates give -0.5 to -0.1 K/day. Over the winter season, the overall lower stratospheric descent rate averages about 1.8-1.5 km/month. These computations are in line with radiative transfer estimates of the heating and descent rate. The HALOE data thus appear to be consistent with the picture of an isolated lower stratospheric Antarctic vortex.

Description: Using delta C-13 measurements in atmospheric CO2 from a cooperative global air sampling network, we determined the partitioning of the net uptake of CO2 between ocean and land as a function of latitude and time. The majority of delta C-13 measurements were made at the Institute of Arctic and Alpine Research (INSTAAR) of the University of Colorado. We perform an inverse deconvolution of both CO2 and delta C-13 observations, using a two-dimensional model of atmospheric transport. Also, the discrimination against C-13 by plant photosynthesis, as a function of latitude and time, is calculated from global runs of the simple biosphere (SiB) model. Uncertainty due to the longitudinal structure of the data, which is not represented by the model, is studied through a bootstrap analysis by adding and omitting measurement sites. The resulting error estimates for our inferred sources and sinks are of the order of 1 GTC (1 GTC = 10(exp 15) gC). Such error bars do not reflect potential systematic errors arising from our estimates of the isotopic disequilibria between the atmosphere and the oceans and biosphere, which are estimated in a separate sensitivity analysis. With respect to global totals for 1992 we found that 3.2 GTC of carbon dissolved into the ocean and that 1.5 GTC were sequestered by land ecosystems. Northern hemisphere ocean gyres north of 15 deg N absorbed 2.7 GTC. The equatorial oceans between 10 deg S and 10 deg N were a net source to the atmosphere of 0.9 GTC. We obtained a sink of 1.6 GTC in southern ocean gyres south of 20 deg S, although the deconvolution is poorly constrained by sparse data coverage at high southern latitudes. The seasonal uptake of CO2 in the northern gyres appears to be correlated with a bloom of phytoplankton in surface waters. On land, northern temperate and boreal ecosystems between 35 deg N and 65 deg N were found to be a major sink of CO2 in 1992, as large as 3.5 GTC. Northern tropical ecosystems (equator-30 deg N) appear to be a net source to the source to the atmosphere of 2 GTC which could reflect biomass burning. A small sink, 0.3 GTC, was inferred for southern tropical ecosystems (30 deg S-equator).

Description: Global sets of surface radiation budget (SRB) have been obtained from satellite programs. These satellite-based estimates need validation with ground-truth observations. This study validates the estimates of monthly mean surface insolation contained in two satellite-based SRB datasets with the surface measurements made at worldwide radiation stations from the Global Energy Balance Archive (GEBA). One dataset was developed from the Earth Radiation Budget Experiment (ERBE) using the algorithm of Li et al. (ERBE/SRB), and the other from the International Satellite Cloud Climatology Project (ISCCP) using the algorithm of Pinker and Laszlo and that of Staylor (GEWEX/SRB). Since the ERBE/SRB data contain the surface net solar radiation only, the values of surface insolation were derived by making use of the surface albedo data contained GEWEX/SRB product. The resulting surface insolation has a bias error near zero and a root-mean-square error (RMSE) between 8 and 28 W/sq m. The RMSE is mainly associated with poor representation of surface observations within a grid cell. When the number of surface observations are sufficient, the random error is estimated to be about 5 W/sq m with present satellite-based estimates. In addition to demonstrating the strength of the retrieving method, the small random error demonstrates how well the ERBE derives from the monthly mean fluxes at the top of the atmosphere (TOA). A larger scatter is found for the comparison of transmissivity than for that of insolation. Month to month comparison of insolation reveals a weak seasonal trend in bias error with an amplitude of about 3 W/sq m. As for the insolation data from the GEWEX/SRB, larger bias errors of 5-10 W/sq m are evident with stronger seasonal trends and almost identical RMSEs.

Description: The topology of the magnetic field in the heliosheath is illustrated using plots of the field lines. It is shown that the Archimedean spiral inside the terminal shock is rotated back in the heliosheath into nested spirals that are advected in the direction of the interstellar wind. The 22-year solar magnetic cycle is imprinted onto these field lines in the form of unipolar magnetic envelopes surrounded by volumes of strongly mixed polarity. Each envelope is defined by the changing tilt of the heliospheric current sheet, which is in turn defined by the boundary of unipolar high-latitude regions on the Sun that shrink to the pole at solar maximum and expand to the equator at solar minimum. The detailed shape of the envelopes is regulated by the solar wind velocity structure in the heliosheath.

Description: Participants in this workshop, which convened in Venice, Italy, 6-8 May 1993, met to consider the current state of climate monitoring programs and instrumentation for the purpose of climatological prediction on short-term (seasonal to interannual) timescales. Data quality and coverage requirements for definition of oceanographic heat and momentum fluxes, scales of inter- and intra-annual variability, and land-ocean-atmosphere exchange processes were examined. Advantages and disadvantages of earth-based and spaceborne monitoring systems were considered, as were the structures for future monitoring networks, research programs, and modeling studies.

Description: The work described here makes it possible to identify anomalous wind behavior such as the nighttime meridional wind abatements that occur at F-region heights. A new analysis technique uses a simple empirical wind model to simulate measurements of 'normal' winds (as measured by the Neutral Atmosphere and Temperature Experiment (NATE) that flew on the Atmosphere Explorer-E (AE-E)) to highlight anomalous wind measurements made by the satellite while in circular orbits at 270-290 km altitude. Our approach is based on the recognition that the 'in orbit' wind variation must show the combined effects of the diurnal wind variation as seen from the ground with the latitude variation of the satellite orbit. For the data period 77250-78035 examined thus far, the wind abatement always occurred with a corresponding pressure or temperature maximum, and was detected on 12 out of the 36 nights with data. This study has revealed that the wind abatement occur only during or shortly after increases in solar extreme ultraviolet (EUV) flux, as indicated by daily radio flux measurements. In the past, nighttime wind reversals at mid-latitudes have been associated with increased geomagnetic activity. This study indicates that intensified solar EUV heating may be responsible for anomalous thermospheric nighttime winds at mid-latitudes.

Description: SN 1987A hard X-ray continuum spectra obtained on 1987 October 29, 1988 April 9-10, and 1988 November 11 from balloon-flight measurments are presented. The spectra, spanning the energy range from 25 keV to 300 keV, have been analyzed using a detector response matrix inversion technique that converts the spectra form counts/s/sq cm keV to photons/s/sq cm keV allowing direct comparison with theoretical models. The results indicate that the bulk of the (56)Co is mixed moderately through the inner regions of the supernova envelope but they do not preclude the mixing of a small amount of the (56)Co farther out into the envelope necessary to account for the observed (56)Co line fluxes. The effect of the ratio (57)Co to (56)Co on the 1988 November 11 continuum spectrum is discussed.

Description: The UV line profile structure of high-ionization resonance lines found with the International Ultraviolet Explorer (IUE) in the brightest of four multiply imaged sources (images-A) in the candidate gravitational lens UM 425 = QSO 1120+019 indicates broad absorption line (BAL) structure. The deep-broad trough associated with the O IV line extends to velocities approiximately -12,000 km/s, and contains disrete features that suggest multicomponent velocity structure. This structure may include contributions from C IV absorption from the early-type galaxy that is believed to lens UM 425. A strong absorption feature in the blue wing of the Lyman-alpha lambda 1216 emission line may be a Lyman alpha absorption system at a Z(sub Ly alpha) = 1.437 +/- 0.003, or it may be formed by the superposition of the broad N V lambda lambda 1238, 1242 absorption trough on the extended blue emission wing of the QSO Lyman-alpha line. We obtained a redshift of Z(sub QSO) = 1.471 +/- 0.003 from Lyman-alpha lambda 1215, consistent with the redshift found by Meylan and Djorgovski in the optical. The Lyman-alpha line appears unusally weak due to the presence of N V lambda 1240 BAL absorption. A Lyman-limit absorption system at lambda 912 was not observed in the QSO rest frame. The detection of BAL structure in the other weaker ground-state resonance lines of N II (l) and S IV (l) was not found, suggesting these lines are formed in a region that is distinct from the BAL component. Detection of BAL structure in the other fainter images in this system with Hubble Space Telescope (HST) instrumentation, similar to structure observed here in image A, could provide evidence that UM 425 is a gravitational lens.

Description: The concept of the well-known Langley plot technique, used for the calibration of ground-based instruments, has been generalized for application to satellite instruments. In polar regions, near summer solstice, the solar backscattered ultraviolet (SBUV) instrument on the Nimbus 7 satellite samples the same ozone field at widely different solar zenith angles. These measurements are compared to assess the long-term drift in the instrument calibration. Although the technique provides only a relative wavelength-to-wavelength calibration, it can be combined with existing techniques to determine the drift of the instrument at any wavelength. Using this technique, we have generated a 12-year data set of ozone vertical profiles from SBUV with an estimated accuracy of +/- 5% at 1 mbar and +/- 2% at 10 mbar (95% confidence) over 12 years. Since the method is insensitive to true changes in the atmospheric ozone profile, it can also be used to compare the calibrations of similar SBUV instruments launched without temporal overlap.

Description: We investigate the physical basis for the timescale of impulsive-phase, redshifted Lyman-alpha emission in stellar flares on the assumption that it is determined by energy losses in a nonthermal proton beam that is penetrating the chromosphere from above. The temporal evolution of ionization and heating in representative model chromospheres subjected to such beams is calculated. The treatment of 'stopping' of beam protons takes into account their interactions with (1) electrons bound in neutral hydrogen, (2) nuclei of neutral hydrogen, (3) free electrons, and (4) ambient thermal protons. We find that, for constant incident beam flux, the system attains an equilibrium with the beam energy input to the chromosphere balanced by radiative losses. In equilibrium, the beam penetration depth is constant, and erosion of the chromosphere ceases. If the redshifted, impulsive-phase stellar flare Lyman-alpha emission is produced by downstreaming hydrogen formed through charge exchange between beam protons and ambient hydrogen, then the emission should end when the beam no longer reaches neutral hydrogen. The durations of representative emission events calculated on this assumption range from 0.1 to 14 s. The stronger the beam, the shorter the timescale over which the redshifted Lyman-alpha emission can be observed.

Description: An automated scheme to characterize precipitation echoes within small windows in the radar field is presented and applied to previously subjectively classified tropical rain cloud systems near Darwin, Australia. The classification parameters are (a) E(sub e), effective efficiency, as determined by cloud-top and cloud-base water vapor saturation mixing ratios; (b) BBF, brightband fraction, as determined by the fraction of the radar echo area in which the maximal reflectivity occurs within +/- 1.5 km of the 0 C isotherm level; and (c) del(sub r) Z, radial reflectivity gradients (dB/km). These classification criteria were applied to tropical rain cloud systems near Darwin, Australia, and to winter convective rain cloud systems in Israel. Both sets of measurements were made with nearly identical networks of C-band radars and rain gauge networks. The results of the application of these objective classification criteria to several independently predetermined rain regimes in Darwin have shown that better organized rain systems have smaller del(sub r) Z and larger BBF. Similarly, smaller del(sub r)Z and larger BBF were also observed from maritime rain cloud systems, as compared to continental rain cloud systems with the same degree of organization. Continental rain cloud system, regardless of their degree of organization, have larger depths, as expressed by E(sub e). The rainfall analyses presented in this study are based exclusively on rain gauge measurements, while radar information was used only to classify the individual gauge measurements.

Description: Through analysis of available optical spectrophotometric data and radio flux density measurements in the literature, it is demonstrated that a good correlation exists between the radio power and bolometric luminosity of the optically-selected OSOs in the Bright Quasar Sample (BOS) of Schmidt and Green (1983). We have recently used VLBI measurements of a sample of ultraluminous infrared galaxies to infer the likely existence of radio-quiet Active Galactic Nuclei (AGNs) deeply enshrouded in dust within their nuclei (Lonsdale, Smith, and Lonsdale 1993). We employ the radio-bolometric luminosity correlation for the BQS quasars to test whether these hypothetical buried AGNs can be energetically responsible for the observed far-infrared luminosities of the ultraluminous infrared galaxies. The ultraluminous infrared galaxies are shown to follow the same relation between radio core power and bolometric luminosity as the radio-quiet QSOs, suggesting that buried AGNs can account for essentially all the observed infrared luminosity, and raising the possibility that any starburst which may be in progress may not be energetically dominant. The broader implications of the radio-optical correlation in quasars for AGNs and luminous infrared galaxy models and the use of radio astronomy as a probe of the central powerhouse in radio quiet AGNs and luminous infrared galaxies are briefly discussed.

Description: The proposed design and construction of the Fizeau astrometric mapping explorer (FAME), a small astrometric instrument for use on an artificial satellite, is reported on. The instrument and spacecraft are designed to slowly spin and will repeatedly scan great circles on the sky so that, over a period of time, it will cover the complete sky and repeat in a manner similar to that of the Hipparcos satellite. The instrument will use the two fixed dilute aperture telescopes to measure a fixed angle between stars and detect the positions, magnitude and color of all stars crossing its field of view to a visual magnitude of approximately 15 mag. The aim of the instrument is to obtain a catalog of positions, proper motions and parallaxes of all stars down to about 15 mag, with a magnitude dependent accuracy of positions of 20 to 800 micro arcsec, proper motions of 20 to 800 micro arcsec per year and parallaxes of 20 to 800 micro arcsec.

Description: The concepts related to the operation and design of the global astrometric interferometer for astrophysics (GAIA) bring together solutions chosen for the astrometry satellite and interferometric techniques. Like the Hipparcos satellite, GAIA is a continuously scanning instrument for which the integration time on any observed object is limited by the field of view of the detector. If a final astrometric accuracy of 10 microarcsec is aimed at, a field of 1 deg in diameter is needed. A design is presented for the proposed 2.6 m baseline Fizeau interferometer with two 40 cm apertures and overall dimensions compatible with the size of the Ariane 5 payload shroud. It has a 0.9 deg diffraction limited field of view. The response of the optical system to small perturbations on each optical element is given in terms of the fringe visibility, which is shown to be dependent on the sub-aperture spot separation. The robustness of the design to thermal, mechanical and manufacturing errors is discussed. The unavoidable distortion present in wide field optical systems is analyzed in terms of displacement of the interference fringes.

Description: A direct link to an extragalactic reference system is considered as being a principle aim of the global astrometric interferometer for astrophysics (GAIA) mission. The data available from an extragalactic data base and a quasi stellar object (QSO) catalog were used to obtain an estimation of the number of QSO link candidates. The quality of presently available data and the expected accuracy of the extragalactic link are discussed. It is concluded that at least 150 QSO's must be observed by GAIA in order to guarantee an accuracy of better than 1 microarcsec/year for the link. New observations will be needed before the GAIA launch in order to reduce uncertainties in the positions, magnitudes and redshifts for some known quasars. The variability of QSO's with magnitudes near the GAIA observation limit can raise a potential problem. The motions of nearby QSO's are expected to be much smaller than 2 microarcsec/year, and therefore, will not affect the accuracy of the link in the proposed GAIA mission.

Description: We present interferometric observations of ice-sheet motion in western Greenland based on pairs of ERS-1 synthetic aperture radar (SAR) images. These observations provide the first detailed regional view of ice motion associated with dynamically supported topography near the margin of an ice sheet. The interferograms of this area are much more complicated than other interferograms of ice sheets presented to date. We devote the largest part of this paper to explaining the source of the complexity in these interferograms. A synthetic interferogram based on a simple model helps to illustrate the effects of different components of the ice velocity field in interferometric data and suggests a method for estimating the large-scale ice velocity field from such data.

Description: We report on large-scale ab initio multiconfiguration Hartree-Fock calculations for the UV0.01 multiplet, 2s(sup 2)2p(sup 2)P(sub J) - 2s2p(sup 2 4)P(sub J prime), in N III. The resulting transition probabilities agree very well with recent semiempirical calculations, and the lifetimes for two of the three upper levels agree with experiments. The deviation for the third level is discussed. Comparisons made with the highest quality IUE echelle spectra available -- those of RR Tel and V1016 Cyg (both photoionized sources with electron densities below 10(exp 8)/cu cm) -- show that computed branching ratios of lines sharing a common upper level are in agreement with observations to within uncertainties of +/- 10%. High-quality solar limb data or stellar data from the Hubble Space Telescope (HST) could, in principle, be used to determine whether the theoretical or measured lifetimes for the discrepant level are in error. Unfortunately, stellar data for high-density plasmas (N(sub e) greater than 10(exp 11)/cu cm are needed) do not yet exist, and existing solar data lack the photometric precision to address this problem.

Description: 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.

Description: High-contrast peaks in the cosmic microwave background (CMB) anisotropy can appear as unresolved sources to observers. We fit simluated CMB maps generated with a cold dark matter model to a set of unresolved features at instrumental resolution 0.5 deg-1.5 deg to derive the integral number density per steradian n (greater than absolute value of T) of features brighter than threshold temperature absolute value of T and compare the results to recent experiments. A typical medium-scale experiment observing 0.001 sr at 0.5 deg resolution would expect to observe one feature brighter than 85 micro-K after convolution with the beam profile, with less than 5% probability to observe a source brighter than 150 micro-K. Increasing the power-law index of primordial density perturbations n from 1 to 1.5 raises these temperature limits absolute value of T by a factor of 2. The MSAM features are in agreement with standard cold dark matter models and are not necessarily evidence for processes beyond the standard model.

Description: We have simulated full-sky maps of the cosmic microwave background (CMB) anisotropy expected from cold dark matter (CDM) models at 0.5 deg and 1.0 deg angular resolution. Statistical properties of the maps are presented as a function of sky coverage, angular resolution, and instrument noise, and the implications of these results for observability of the Doppler peak are discussed. The rms fluctuations in a map are not a particularly robust probe of the existence of a Doppler peak; however, a full correlation analysis can provide reasonable sensitivity. We find that sensitivity to the Doppler peak depends primarily on the fraction of sky covered, and only secondarily on the angular resolution and noise level. Color plates of the simulated maps are presented to illustrate the anisotropies.

Description: Results of a simulation study on the effects of optically thick stratospheric sulfate aerosol layers on the backscattered ultraviolet radiation (buv) in the range 256-340 nm are presented. In general, the increased Mie scattering produced by the aerosols results in an enhancement of the buv radiation. The increase is approximately linear with optical depth and strongly depends on solar zenith angle and aerosol layer altitude in relation to the ozone maximum. The effect is greatest at those wavelenghts whose contribution functions peak in the vicinity of the densest part of the aerosol layer. The aerosol induced perturbation of the buv field affects the ozone profile retrieval from space measurements by the Solar Backscatter Ultraviolet Experiment (SBUV) experiment. In the tropical stratosphere, the retrieved ozone between 25 and 45 km is underestimated as a result of increased Mie scattering. On the other hand, an algorithm related effect causes the retrieved ozone below 25 km to be overestimated by an amount similar to the stratospheric deficit.

Description: Soft gamma repeaters are characterized by recurrent activity consisting of short-duration outbursts of high-energy emission that is typically of temperature less than 40 keV. One recent model of repeaters is that they originate in the environs of neutron stars with superstrong magnetic fields, perhaps greater than 10(exp 14) G. In such fields, the exotic process of magnetic photon splitting gamma yields gamma gamma acts very effectively to reprocess gamma-ray radiation down to hard X-ray energies. In this Letter, the action of photon splitting is considered in some detail, via the solution of photon kinetic equations, determining how it limits the hardness of emission in strongly magnetized repeaters, and thereby obtaining observational constraints to the field in SGR 1806-20.

Description: The COMPTEL sky map of the 1.8 MeV line emission from Al-26 shows an extended diffuse distribution along the Galactic plane with a peculiar large scale asymmetry about the Galactic center (GC) and a clumpy structure with several noticeable hot spots. The most prominent hot spot at the GC appears shifted to positive longitude by about 2 deg. Nearby supernova remnants or Wolf-Rayet stars are plausible explanations for individual hot spots such as in the Vela region. We show that the global asymmetry and most hot spots can be understood by a more general model: the Al-26 sources are confined in the spiral arms of our Galaxy, and coagent star formation is responsible for additional clumpiness. The off-centered GC hot spot is probably due to the observed asymmetrically distributed circumnuclear molecular gas interacting with a central stellar bar, whose existence has been inferred from IR observations.

Description: Application of the window probability matching method to radar and rain gauge data that have been objectively classified into different rain types resulted in distinctly different Z(sub e)-R relationships for the various classifications. These classification parameters, in addition to the range from the radar, are (a) the horizontal radial reflectivity gradients (dB/km); (b) the cloud depth, as scaled by the effective efficiency; (c) the brightband fraction within the radar field window; and (d) the height of the freezing level. Combining physical parameters to identify the type of precipitation and statistical relations most appropriate to the precipitation types results in considerable improvement of both point and areal rainfall measurements. A limiting factor in the assessment of the improved accuracy is the inherent variance between the true rain intensity at the radar measured volume and the rain intensity at the mouth of the rain guage. Therefore, a very dense rain gauge network is required to validate most of the suggested realized improvement. A rather small sample size is required to achieve a stable Z(sub e)-R relationship (standard deviation of 15% of R for a given Z(sub e)) -- about 200 mm of rainfall accumulated in all guages combined for each classification.

Description: Two dimensional spectral types for each of the stars observed in the global astrometric interferometer for astrophysics (GAIA) mission would provide additional information for the galactic structure and stellar evolution studies, as well as helping in the identification of unusual objects and populations. The classification of the large quantity generated spectra requires that automated techniques are implemented. Approaches for the automatic classification are reviewed, and a metric-distance method is discussed. In tests, the metric-distance method produced spectral types with mean errors comparable to those of human classifiers working at similar resolution. Data and equipment requirements for an automated classification survey, are discussed. A program of auxiliary observations is proposed to yield spectral types and radial velocities for the GAIA-observed stars.

Description: Low temperature Zn diffusion in GaSb, where the minimum temperature was 450 C, was studied. The pseudo-closed box (PCB) method was used for Zn diffusion into GaAs, AlGaAs, InP, InGaAs and InGaAsP. The PCB method avoids the inconvenience of sealed ampoules and proved to be simple and reproducible. The special design of the boat for Zn diffusion ensured the uniformality of Zn vapor pressure across the wafer surface, and thus the uniformity of the p-GaSb layer depth. The p-GaSb layers were studied using Raman scattering spectroscopy and the x-ray rocking curve method. As for the postdiffusion processing, an anodic oxidation was used for a precise thinning of the diffused GaSb layers. The results show the applicability of the PCB method for the large-scale production of the GaSb structures for solar cells.

Description: Indium phosphide is considered to be a strong contender for many photovoltaic space applications because of its radiation resistance and its potential for high efficiency. An overview of recent progress is presented, and possible future research directions for indium phosphide space solar cells are discussed. The topics considered include radiation damage studies and space flight experiments.

Description: Understanding solar cell response to pulsed laser outputs is important for the evaluation of power beaming applications. The time response of high efficiency GaAs and silicon solar cells to a 25 nS monochromatic pulse input is described. The PC-1D computer code is used to analyze the cell current during and after the pulse for various conditions.

Description: Four-terminal mechanically stacked solar cells were developed for advanced space arrays with line-focus reflective concentrators. The top cells are based on AlGaAs/GaAs multilayer heterostructures prepared by low temperature liquid phase epitaxy. The bottom cells are based on heteroepitaxial InP/InGaAs liquid phase epitaxy or on homo-junction GaSb, Zn-diffused structures. The sum of the highest reached efficiencies of the top and bottom cells is 29.4 percent. The best four-terminal tandems have an efficiency of 27 to 28 percent. Solar cells were irradiated with 1 MeV electrons and their performances were determined as a function of fluence up to 10(exp 16) cm(exp-2). It was shown that the radiation resistance of developed tandem cells is similar to the most radiative stable AlGaAs/GaAs cells with a thin p-GaAs photoactive layer.

Description: An investigation of solar cells based on AlGaAs/GaAs heterostructures with an internal Bragg reflector as the back-surface reflector is presented. The Bragg reflector is grown by low pressure metalorganic chemical vapor deposition on n-GaAs substrates in a horizontal resistively heated reactor. The Bragg reflector with its maximum reflectance centered at a wavelength of 860 nm consists of 12 pairs of AlAs/GaAs layers. The resulting Bragg reflector has a thickness of 0.072 micrometers for AlAs and 0.059 micrometers for GaAs. The multi-layered quasi-dielectric stack selectively reflects weakly absorbed photons with energies near to the GaAs band gap for a second pass through the photoactive region, thus increasing the photocurrent. The use of the Bragg reflector allows the external quantum efficiency to be increased in the long wavelength of the spectrum. The use of the Bragg reflector and an antireflective coating and prismatic cover allowed an efficiency of 23.4 percent to be obtained.

Description: The high efficiency solar cells based on multilayer AlGaAs/GaAs heterostructures, prepared by low temperature liquid phase epitaxy (LPE), were developed and tested. An investigation of the low temperature LPE process for the crystallization of AlGaAs heterostructures of as high as 24.0 to 24.7 percent under AMO conditions at concentration ratios of 20 to 100x, were reached. Developed solar cells show substantial radiation resistance to the damage induced by 3.75 MeV electrons.

Description: We consider the evolution and heating of dust embedded in the hot interstellar medium of isolated elliptical galaxies. We first construct a new set of galaxy models spanning a decade in luminosity which we use to study the evolution of dust as it is ejected by stellar sources, merges and flows with the interstellar gas, and is sputtered away due to the presence of the hot gas. We find that although grains can flow a considerable distance from the parent star in some cases before being sputtering away, the grain size distribution at a given location is accurately determined by assuming in situ sputtering of dust grains as they are ejected by the parent star. We find that dust heating is dominated by absorption of ambient starlight. Heating due to collisions with energetic electrons in the hot gas and absorption of thermal X-rays is smaller by orders of magnitude. We also find that for the largest galaxy considered (L(sub B) approximately 10 (exp 11) solar luminosity) the energy lost from the hot gas due to electronic collisions with the dust is about an order of magnitude less than that lost due to thermal emission in X-rays. The gas in smaller galaxies loses even less of its energy via heating of grains. In deriving this result, we find that the inclusion of grain sputtering is a crucial ingredient of the calculation. The presence of dust grains is therefore energetically unimportant for elliptical galaxy cooling flow models.

Description: We have studied the implications of recent trends in the annual mean and the amplitude of the annual harmonic of ozone in the upper stratosphere from the 15 years of the combined data from the Nimbus-7 SBUV and the NOAA-11 SBUV/2 instruments. This was done in the context of the GSFC 2D model predictions of these trends which are based on plausible scenarios of anthropogenic Cly increase in the atmosphere. The comparison of the observed and model-estimated annual mean ozone trends show some similarity in their latitude and altitude characteristics. Both the model and data show a maximum ozone decrease of -6 to -10 % per decade at high latitudes in the upper stratosphere. However, there are also significant differences between the observed and computed trends which may be related to both the model uncertainty and the uncertainty in correcting for the long term instrument drift. The observations also suggest a decrease of 10-25 % per decade in the annual amplitude of ozone at 2 mb between 40 deg - 60 deg in both hemispheres, with a relatively larger interannual variability in the northern hemisphere. These values are in general agreement with the model predictions and thus provide additional support in favor of the chlorine induced changes in ozone in the upper stratosphere.

Description: BrO measured from the NASA ER-2 during Airborne Arctic Stratospheric Expedition (AASE) II exhibited a mean value (for 20-minute averages) of 5.4 +/- 0.3 pptv, with a standard deviation of 3.1 pptv. Ratios of BrO to available inorganic bromine (Br(sub y)) show only slight increases in polar regions relative to midlatitudes. A comparison between observed latitudinal and diurnal variations of this same ratio and that calculated by photochemical models shows reasonable agreement in behavior, but significant discrepancies in magnitude. It is unclear whether this difference is due to errors in measurements, models or both.

Description: Measurements of stratospheric hydrofluoric acid (HF) have been made by the JPL MkIV interferometer during high-altitude balloon flights. Infrared solar absorption spectra were acquired near 35 deg N at altitudes between local tropopause and 38 km. Volume mixing ratio profiles of HF derived from 4 flights (1990-93), in conjunction with simultaneously observed N2O profiles, indicate an average rate of HF increase of (5.5 +/- 0.3)% per year, in agreement with time-dependent, two-dimensional model simulations (6% per year) and ATMOS measurements.

Description: Near-infrared images of the Galactic bulge at 1.25, 2.2, 3.5, and 4.9 microns obtained by the Diffuse Infrared Background Experiment (DIRBE) onboard the Cosmic Background Explorer (COBE) satellite are used to characterize its morphology and to determine its infrared luminosity and mass. Earlier analysis of the DIRBE observations (Weiland et al. 1994) provided supporting evidence for the claim made by Blitz & Spergel (1991) that the bulge is bar-shaped with its near end in the first Galactic quadrant. Adopting various triaxial analytical functions to represent the volume emissivity of the source, we confirm the barlike nature of the bulge and show that triaxial Gaussian-type functions provide a better fit to the data than other classes of functions, including an axisymmetric spheroid. The introduction of a `boxy' geometry, such as the one used by Kent, Dame, & Fazio (1991) improves the fit to the data. Our results show that the bar is rotated in the plane with its near side in the first Galactic quadrant creating an angle of 20 deg +/- 10 deg between its major axis and the line of sight to the Galactic center. Typical axis ratios of the bar are (1:0.33 +/- 0.11:0.23 +/- 0.08), resembling the geometry of prolate spheroids. There is no statistically significant evidence for an out-of-plane tilt of the bar at 2.2 microns, and marginal evidence for a tilt of approximately equal 2 deg at 4.9 microns. The introduction of a roll around the intrinsic major axis of the bulge improves the `boxy' appearance of some functions. A simple integration of the observed projected intensity of the bulge gives a bulge luminosity of 1.2 x 10(exp 9), 4.1 x 10(exp 8), 2.3 x 10(exp 8), and 4.3 x 10(exp 7) solar luminosity, respectively, at 1.25, 2.2, 3.5, and 4.9 microns wavelength for a Galactocentric distance of 8.5 kpc. The 2.2 microns luminosity function of the bulge population in the direction of Baade's window yields a bolometric luminosity of L(sub bol) = 5.3 x 10(exp 9) solar luminosity. Stellar evolutionary models relate this luminosity to the number of main-sequence progenitor stars that currently populate the red giant branch. Combined with the recent determination of the main-sequence turnoff mass for the bulge by the Hubble Space Telescope (Holtzman et al. 1993) we derive a photometrically determined bulge mass of approximately equal to 1.3 x 10(exp 10) solar mass for a Salpeter initial mass function extended down to 0.1 solar mass.

Description: We report the results from observations of National Radio Astronomy Observatory (NRAO) 140 carried out in the X-ray band using ROSAT and Advanced Satellite for Cosmology and Astrophysics (ASCA) in 1992 and 1994. We find the source to be absorbed by an effective hydrogen column density of approximately 3 x 10(exp 21) atoms/sq cm during both epochs, consistent with the combined atomic and molecular column inferred from radio measurements along this line of sight through the Perseus cloud complex. We compare these results with previous X-ray absorption measurements and briefly comment upon the origin of the excess absorption which has been seen toward this source. We find the ASCA spectrum of NRAO 140 is well described by a power law of energy index alpha = 0.73 +/- 0.03 and also yields the tightest constraint to date on Fe K-shell emission, with 90% confidence upper limits of 38 and 31 eV for a narrow line at a rest frame energy of 6.4 and 6.7 keV respectively. This, along with a lack of hardening towards higher energies, suggests that either NRAO 140 is devoid of cold reprocessing material, the reprocessing material has a geometry in which the imprinted features are weak, and/or the X-ray emission is relativistically beamed towards us. We also report the detection of a serendipitous source in both the Advanced Satellite for Cosmology and Astrophysics Gas Imaging Spectrometer (ASCA GIS) and ROSAT fields of view. We identify this source as the binary star system IX Per and find its spectrum to be well fitted by a two-temperature Raymond-Smith plasma.

Description: We have produced two-dimensional maps of the intensity ratio, Q(sub 60), of 60 micron infrared to 20 cm radio continuum emission, for a set of 25 nearby galaxies, mostly spirals. The ratio maps were obtained from infrared images made using IRAS data with the maximum correlation method, and radio images made using VLA data. Before taking the ratio, the radio images were processed so as to have the same resolution properties as the infrared images; the final spatial resolution in all cases is approximately 1 min, corresponding to 1 - 2 kpc for most galaxies. This resolution represents a significant improvement over previous studies. Our new high-resolution maps confirm the slow decrease of Q(sub 60) with increasing radial distance from the nucleus, but show additional structure which is probably associated with separate sites of active star formation in the spiral arms. The maps show Q(sub 60) to be more closely related to infrared surface brightness than to the radial distance r in the galaxy disk. We note also that the Q(sub 60) gradients are absent (or at least reduced) for the edge-on galaxies, a property which can be attributed to the dilution of contrast due to the averaging of the additional structure along the line of sight. The results are all in qualitative agreement with the suggestion that the radio image represents a smeared version of the infrared image, as would be expected on the basis of current models in which the infrared-radio correlation is driven by the formation of massive stars, and the intensity distribution of radio emission is smeared as a result of the propagation of energetic electrons accelerated during the supernova phase.

Description: Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS) measurements of lower stratospheric ClO during 1992-93 and 1993-94 Arctic winters are presented. Enhanced ClO in the 1992-93 winter was first observed in early December, and extensively during February when temperatures were continually low enough for polar stratospheric cloud (PSCs). Sporadic episodes of enhanced ClO were observed for most of the 1993-94 winter as minimum temperatures hovered near the PSC threshold, with largest ClO amounts occurring in early March after a sudden deep cooling in late February.

Description: The U-B color has been suggested as a predictor of the presence of water of hydration on asteroids. Photometry from the Eight-Color Asteroid Survey (ECAS) was used to test this concept. An overlap in U-B color prevents this magnitude difference from distinguishing between surface material that was thermally processed at higher temperatures and surface material that was aqueously altered. Two tests of the presence of water of hydration using visible spectral region photometry failed to flag those few higher albedo M- and E-class asteroids having photometry that shows a 3.0-micrometers water of hydration absorption. These asteroids probably contain little or no oxidized iron in their surface material.

Description: The thermal-infrared (longwave) emission from a vegetated terrain is generally anisotropic, i.e., the emission temperature varies with the view direction. If a directional measurement of temperature is considered to be equal to the effective temperature of the hemispheric emission, then the estimate of the latter can be significantly in error. The view-direction (zenith angle theta(sub eq) at which the emission equivalence does hold is determined in our modeling study. In a two-temperature field-of-view (soil and plants), theta(sub eq) falls in a narrow range depending on plant density and canopy architecture. Theta(sub eq) does not depend on soil and (uniform) plant temperatures nor on their ratio, even though the pattern of emission vs. the view direction depends crucially on this ratio. For a sparse canopy represented as thin, vertical cylindrical stalks (or vertical blades uniformly distributed in azimuth) with horizontal facets, theta(sub eq) ranges from 48 to 53 deg depending on the optical density of the vertical elements alone. When plant elements are modeled as small spheres, theta(sub eq) lies between 53 to 57 deg (for the same values of the canopy optical density). Only for horizontal leaves (a truly planophile canopy) is the temperature measured from any direction equal to the temperature of the hemispheric emission. When the emission temperature changes with optical depth within the canopy at a specified rate, theta(sub eq) depends to some extent on that rate. For practically any sparsely vegetated surface, a directional measurement at the zenith angle of 50 deg offers an appropriate evaluation of the hemispheric emission, since the error in the estimate will, at most, only slightly exceed 1% (around 4 W/sq m). Estimates of the hemispheric emission through a nadir measurement, on the other hand, can be in error in some cases by about 10%, i.e., on the order of 40 W/sq m.

Description: The Monte Carlo technique of simulating diffusive particle acceleration at shocks has made spectral predictions that compare extremely well with particle distributions observed at the quasi-parallel region of the earth's bow shock. The current extension of this work to compare simulation predictions with particle spectra at oblique interplanetary shocks has required the inclusion of significant cross-field diffusion (strong scattering) in the simulation technique, since oblique shocks are intrinsically inefficient in the limit of weak scattering. In this paper, we present results from the method we have developed for the inclusion of cross-field diffusion in our simulations, namely model predictions of particle spectra downstream of oblique subluminal shocks. While the high-energy spectral index is independent of the shock obliquity and the strength of the scattering, the latter is observed to profoundly influence the efficiency of injection of cosmic rays into the acceleration process.

Description: Numerical simulation experiments were conducted to delineate the influence of in situ deforestation data on episodic rainfall by comparing two ensembles of five 5-day integrations performed with a recent version of the Goddard Laboratory for Atmospheres General Circulation Model (GCM) that has a simple biosphere model (SiB). The first set, called control cases, used the standard SiB vegetation cover (comprising 12 biomes) and assumed a fully forested Amazonia, while the second set, called deforestation cases, distinguished the partially deforested regions of Amazonia as savanna. Except for this difference, all other initial and prescribed boundary conditions were kept identical in both sets of integrations. The differential analyses of these five cases show the following local effects of deforestation. (1) A discernible decrease in evapotranspiration of about 0.80 mm/d (roughly 18%) that is quite robust in the averages for 1-, 2-, and 5-day forecasts. (2) A decrease in precipitation of about 1.18 mm/d (roughly 8%) that begins to emerge even in 1-2 day averages and exhibits complex evolution that extends downstream with the winds. (3) A significant decrease in the surface drag force (as a consequence of reduced surface roughness of deforested regions) that, in turn, affects the dynamical structure of moisture convergence and circulation. The surface winds increase significantly during the first day, and thereafter the increase is well maintained even in the 2- and 5-day averages.

Description: A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Description: Four spectral emission features of the N(sub KK) = 4(sub 04) -3(sub 13) rotational transition of methylene (CH2) have been detected at signal levels 5-7 sigma above noise toward the hot core of the Orion-KL nebula and the molecular cloud in proximity to the continuum source W51 M. Specifically, in both sources we have resolved the F = 6-5, 5-4, and 4-3 hyperfine transitions of the J = 5-4 fine-structure levels and detected the blended hyperfine structure of the J = 4-3 fine structure levels. At the J = 3-2 fine-structure levels, we have observed new transitions of NS, a known interstellar molecule, which severely contaminates the search for CH2 hyperfine transitions. These new sensitive observations finally confirm the existence of interstellar CH2 which was tentatively reported by us some years ago.

Description: The star, Altair (A7 IV-V), is clearly shown to have Lyman-alpha emission of chromospheric origin, while no evidence is found for the Mg II emission reported in previous investigations. We present non-Local Thermodymanic Equilibrium (non-LTE) semiempirical models incorporating partial redistribution of the chromosphere of Altair that reproduce the observed Lyman-alpha emission and the Mg II resonance absorption at 2800 A. We unambiguously establihed that chromospheres exist at spectral types as early as A7 on the main sequence, and we also demonstrate that it very unlikely that the observed emission originates in a corotating expanding wind. This result represents a new challenge for chromospheric heating theories. It may indicate that both differential rotation and convection layers, at least near the equator, exist in this fast rotating (v sin i = 220 km/s) star.

Description: A steady-state scheme for data assimilation in the context of a single, short period (relative to a day), sun-synchronous, polar-orbiting satellite is examined. If the satellite takes observations continuously, the gains, which are the weights for blending observations and predictions together, are steady in time. For a linear system forced by random noise, the optimal steady-state gains (Wiener gains) are equivalent to those of a Kalman filter. Computing the Kalman gains increases the computational cost of the model by a large factor, but computing the Wiener gains does not. The latter are computed by iteration using prior estimates of the gains to assimilate simulated observations of one run of the model, termed 'truth' into another run termed 'prediction'. At each stage, the prediction errors form the basis for the next estimate of the gains. Steady state is achieved after three or four iterations. Further simplification is achieved by making the gains depend on longitudinal distance from the observation point, not on absolute longitude. For a single-layer primitive equation model, the scheme works well even if only the mass field is observed but not the velocity field. Although the scheme was developed for Mars Observer, it should be applicable to data retrieved from Earth atmosphere satellites, for example, UARS.

Description: Brightness temperature difference (BTD) values are calculated for selected Geostationary Operational Environmental Satellite (GOES-6) channels (3.9, 12.7 micrometer) and Advanced Very High Resolution Radiometer channels (3.7, 12.0 micrometer). Daytime and nighttime discrimination of particle size information is possible given the infrared cloud extinction optical depth and the BTD value. BTD values are presented and compared for cirrus clouds composed of equivalent ice spheres (volume, surface area) versus randomly oriented hexagonal ice crystals. The effect of the hexagonal ice crystals is to increase the magnitude of the BTD values calculated relative to equivalent ice sphere (volume, surface area) BTDs. Equivalent spheres (volume or surface area) do not do a very good job of modeling hexagonal ice crystal effects on BTDs; however, the use of composite spheres improves the simulation and offers interesting prospects. Careful consideration of the number of Legendre polynomial coefficients used to fit the scattering phase functions is crucial to realistic modeling of cirrus BTDs. Surface and view-angle effects are incorporated to provide more realistic simulation.

Description: We use the two-point correlation function of the extrema points (peaks and valleys) in the Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) 2 year sky maps as a test for non-Gaussian temperature distribution in the cosmic microwave background anisotropy. A maximum-likelihood analysis compares the DMR data to n = 1 toy models whose random-phase spherical harmonic components a(sub lm) are drawn from either Gaussian, chi-square, or log-normal parent populations. The likelihood of the 53 GHz (A+B)/2 data is greatest for the exact Gaussian model. There is less than 10% chance that the non-Gaussian models tested describe the DMR data, limited primarily by type II errors in the statistical inference. The extrema correlation function is a stronger test for this class of non-Gaussian models than topological statistics such as the genus.

Description: An examination and analysis of video images of lightning, captured by the payload bay TV cameras of the space shuttle, provided a variety of examples of lightning in the stratosphere above thunderstorms. These images were obtained on several recent shuttle flights while conducting the Mesoscale Lightning Experiment (MLE). The images of stratospheric lightning illustrate the variety of filamentary and broad vertical discharges in the stratosphere that may accompany a lightning flash. A typical event is imaged as a single or multiple filament extending 30 to 40 km above a thunderstorm that is illuminated by a series of lightning strokes. Examples are found in temperate and tropical areas, over the oceans, and over the land.

Description: During its full-sky survey, the Energetic Gamma Ray Experiment Telescope (EGRET) high-energy instrument aboard the Compton Gamma-Ray Observatory detected gamma-ray emission in the energy range above 30 MeV from a source identified as the blazar PKS 0420-014. This object was observed during two separate viewing periods in 1992 February/March and May/June. An intensity decrease above 100 MeV of a factor of at least 1.5 from a maximum of (5.0 +/- 1.4) x 10(exp -7) photons 1 sq cm/s was observed during that time interval indicating extensive variability. The photon spectrum in the range between 30 and 10,000 MeV at the time of the maximum intensity is well represented by a power law with an exponent of -1.9 +/- 0.3. Some similarities with other EGRET detected blazars are briefly discussed.

Description: Get-Away Special (GAS) G-301, named the Flying Falcon and scheduled for launch on the STS-77 Space Shuttle in April, 1996, is being prepared to perform an experiment designed by the Department of Geology, Bowling Green State University (BGSU). The experiment will employ a new type of infrared imager designed and built by a consortium of Teltron Technologies Inc., Hudson Research Inc., and BGSU that is an uncooled, quantrum ferro-electric, infrared return beam vidicon (IRBV) camera capable of detecting thermal infrared radiation throughout the 2.0-50.0 micron wavelength region, and to which an integral, unable Fabry-Perot filter and a telescopic lens have been added. The primary objectives in the experiment include the mapping of methane plumes from solid waste landfills and wetlands in the midwestern U.S., the mapping of methane plumes offshore in the Gulf of Mexico and in the Middle East, brief monitoring for precursors of volcanoes or earthquakes in the South China sea and the East Pacific Rise (about 300 km west of Easter Island), and the mapping of silica content in exposed outcrops and residual soils of the southwestern U.S. and Middle East.

Description: The Thermal Energy Storage-1 (TES-1) is a flight experiment that flew on the Space Shuttle Columbia (STS-62), in March 1994, as part of the OAST-2 mission. TES-1 is the first experiment in a four experiment suite designed to provide data for understanding the long duration microgravity behavior of thermal energy storage fluoride salts that undergo repeated melting and freezing. Such data have never been obtained before and have direct application for the development of space-based solar dynamic (SD) power systems. These power systems will store solar energy in a thermal energy salt such as lithium fluoride or calcium fluoride. The stored energy is extracted during the shade portion of the orbit. This enables the solar dynamic power system to provide constant electrical power over the entire orbit. Analytical computer codes have been developed for predicting performance of a spaced-based solar dynamic power system. Experimental verification of the analytical predictions is needed prior to using the analytical results for future space power design applications. The four TES flight experiments will be used to obtain the needed experimental data. This paper will focus on the flight results from the first experiment, TES-1, in comparison to the predicted results from the Thermal Energy Storage Simulation (TESSIM) analytical computer code. The TES-1 conceptual development, hardware design, final development, and system verification testing were accomplished at the NASA lewis Research Center (LeRC). TES-1 was developed under the In-Space Technology Experiment Program (IN-STEP), which sponsors NASA, industry, and university flight experiments designed to enable and enhance space flight technology. The IN-STEP Program is sponsored by the Office of Space Access and Technology (OSAT).

Description: A prototype of the Josephson-effect spectrum analyzer developed for the millimeter wave band is described. The measurement results for spectra obtained in the frequency band from 50 to 250 GHz are presented.

Description: Monitoring spatial and temporal changes of soil moisture are of importance to hydrology, meteorology, and agriculture. This paper reports a result on study of using L-band SAR imagery to estimate soil moisture and surface roughness for bare fields. Due to limitations of the Small Perturbation Model, it is difficult to apply this model on estimation of soil moisture and surface roughness directly. In this study, we show a simplified model derived from the Integral Equation Model for estimation of soil moisture and surface roughness. We show a test of this model using JPL L-band AIRSAR data.

Description: In the tropical rain forests of Manu, in Peru, where forest biomass ranges from 4 kg/sq m in young forest succession up to 100 kg/sq m in old, undisturbed floodplain stands, the P-band polarimetric radar data gathered in June of 1993 by the AIRSAR (Airborne Synthetic Aperture Radar) instrument separate most major vegetation formations and also perform better than expected in estimating woody biomass. The worldwide need for large scale, updated biomass estimates, achieved with a uniformly applied method, as well as reliable maps of land cover, justifies a more in-depth exploration of long wavelength imaging radar applications for tropical forests inventories.

Description: Radar backscatter intensity as measured by calibrated synthetic aperture radar (SAR) systems is primarily controlled by three factors: local incidence angle, wavelength-scale roughness, and dielectric permittivity of surface materials. In order to make adequate use of radar observations for geological investigations of surface type, the relationships between lithology and the above characteristics must be adequately understood. In arid terrains weathering signatures (e.g. fracturing, debris grain size and shape, slope characteristics) are controlled to some extent by lithologic characteristics of the parent bedrock. These textural features of outcrops and their associated debris control radar backscatter to varying degrees. The quad-polarization JPL AIRSAR system allows sampling of textures at three distinct wavelength scales: C-band (5.66 cm), L-band (23.98 cm), and P-band (68.13 cm). This paper presents a discussion of AIRSAR data using recent field observations of weathered felsic and basaltic volcanic rock units exposed in the southern part of the Lunar Crater Volcanic Field, in the Pancake Range of central Nevada. The focus is on the relationship of radar backscatter at multiple wavelengths to weathering style and parent bedrock lithology.

Description: Silicic (silica-rich) lava flows, such as rhyolite, rhyodacite, and dacite, possess unique physical properties primarily because of the relatively high viscosity of the molten lava. Silicic flows tend to be thicker than basaltic flows, and the resulting large-scale morphology is typically a steep-sided dome or flow lobe, with aspect ratios (height/length) sometimes approaching unity. The upper surfaces of silicic domes and flows are normally emplaced as relatively cool, brittle slabs that fracture as they are extruded from the central vent areas, and are then rafted away toward the flow margin as a brittle carapace above a more ductile interior layer. This mode of emplacement results in a surface with unique roughness characteristics, which can be well-characterized by multiparameter synthetic aperture radar (SAR) observations. In this paper, we examine the scattering properties of several silicic domes in the Inyo volcanic chain in the Eastern Sierra of California, using AIRSAR and TOPSAR data. Field measurements of intermediate-scale (cm to tens of m) surface topography and block size are used to assess the mechanisms of the scattering process, and to quantify the unique roughness characteristics of the flow surfaces.

Description: The replacement of semidesert grassland by woody shrubland is a widespread form of desertification. This change in physiognomy and species composition tends to sharply reduce the productivity of the land for grazing by domestic livestock, increase soil erosion and reduce soil fertility, and greatly alter many other aspects of ecosystem structure and functioning. Remote sensing methods are needed to assess and monitor shrubland encroachment. Detection of woody shrubs at low density would provide a particularly useful baseline on which to access changes, because an initially low shrub density often tends to increase even after cessation of the disturbance (e.g., overgrazing, drought, or fire suppression) responsible for triggering the initial stages of the invasion (Grover and Musick, 1990). Limited success has been achieved using optical remote sensing. In contrast to other forms of desertification, biomass does not consistently decrease with a shift from grassland to shrubland. Estimation of green vegetation amount (e.g., by NDVI) is thus of limited utility, unless the shrubs and herbaceous plants differ consistently in phenology and the area can be viewed during a season when only one of these is green. The objective of this study was to determine if the potential sensitivity of active microwave remote sensing to vegetation structure could be used to assess the degree of shrub invasion of grassland. Polarimetric Airborne Synthetic Aperture Radar (AIRSAR) data were acquired for a semiarid site containing varied mixtures of shrubs and herbaceous vegetation and compared with ground observations of vegetation type and other landsurface characteristics. In this preliminary report we examine the response of radar backscatter intensity to shrub density. The response of other multipolarization parameters will be examined in future work.

Description: Conventional representations of polarization response are referred to a horizontally and vertically polarized basis. Recent studies by Freeman and Durden, van Zyl, and others suggest that alternative polarimetric features which more easily resolve the contributions of simple scattering mechanisms such as odd-bounce, even-bounce, and diffuse scattering could offer several advantages in terrain classification. The circular polarization covariance matrix is a potential source of such features. In this paper, we derive its relationship to the Stokes matrix, describe some of its properties, and compare the utility of linear and circular polarimetric features in classifying an AIRSAR scene containing urban, park, and ocean terrain.

Description: The relationship between habitat area, spatial dynamics of the landscape, and species diversity is an important theme in population and conservation biology. Of particular interest is how populations of various species are affected by increasing habitat edges due to fragmentation. Over the last decade, assumptions regarding the effects of habitat edges on biodiversity have fluctuated wildly, from the belief that they have a positive effect to the belief that they have a clearly negative effect. This change in viewpoint has been brought about by an increasing recognition of the importance of geographic scale and a reinterpretation of natural history observations. In this preliminary report from an ongoing project, we explore the use of remote sensing technology and geographic information systems to further our understanding of how species diversity and population density are affected by habitat heterogeneity and landscape composition. A primary feature of this study is the investigation of SAR for making more rigorous investigations of habitat structure by exploiting the interaction between radar backscatter and vegetation structure and biomass. A major emphasis will be on the use of SAR data to define relative structural types based on measures of structural consolidation using the vegetation surface area to volume ratio (SA/V). Past research has shown that SAR may be sensitive to this form of structural expression which may affect biodiversity.

Description: Landforms in arid regions record the interplay between tectonic forces and climate. Alluvial fans are a common landform in desert regions where the rate of uplift is greater than weathering or sedimentation. Changes in uplift rate or climatic conditions can lead to isolation of the currently forming fan surface through entrenchment and construction of another fan either further from the mountain front (decreased uplift or increased runoff) or closer to the mountain front (increased uplift or decreased runoff). Thus, many alluvial fans are made up of a mosaic of fan units of different age, some older than 1 million years. For this reason, determination of the stages of fan evolution can lead to a history of uplift and runoff. In an attempt to separate the effects of tectonic (uplift) and climatic (weathering, runoff, sedimentation) processes on the shapes of alluvial fan units, a modified conic equation developed by Troeh (1965) was fitted to TOPSAR digital topographic data for the Trail Canyon alluvial fan in Death Valley, California. This allows parameters for the apex position, slope, and radial curvature to be compared with unit age.

Description: During recent years signature analysis, classification, and modeling of Synthetic Aperture Radar (SAR) data as well as estimation of geophysical parameters from SAR data have received a great deal of interest. An important requirement for the quantitative use of SAR data is the accurate estimation of the backscattering coefficient sigma(exp 0). In terrain with relief variations radar signals are distorted due to the projection of the scene topography into the slant range-Doppler plane. The effect of these variations is to change the physical size of the scattering area, leading to errors in the radar backscatter values and incidence angle. For this reason the local incidence angle, derived from sensor position and Digital Elevation Model (DEM) data must always be considered. Especially in the airborne case, the antenna gain pattern can be an additional source of radiometric error, because the radar look angle is not known precisely as a result of the the aircraft motions and the local surface topography. Consequently, radiometric distortions due to the antenna gain pattern must also be corrected for each resolution cell, by taking into account aircraft displacements (position and attitude) and position of the backscatter element, defined by the DEM data. In this paper, a method to derive an accurate estimation of the backscattering coefficient using NASA/JPL AIRSAR data is presented. The results are evaluated in terms of geometric accuracy, radiometric variations of sigma(exp 0), and precision of the estimated forest biomass.

Description: An empirical model was developed to infer soil moisture and surface roughness from radar data. The accuracy of the inversion technique is assessed by comparing soil moisture obtained with the inversion technique to in situ measurements. The effect of vegetation on the inversion is studied and a method to eliminate the areas where vegetation impairs the algorithm is described.

Description: In 1993 NASA's Jet Propulsion Laboratory Airborne Synthetic Aperture Radar system (AIRSAR) was deployed to South America to collect multi-parameter radar data over pre-selected targets. Among the sites targeted was a series of wind streaks located in the Altiplano of Bolivia. The objective of this investigation is to study the effect of wavelength, polarization, and incidence angle on the visibility of wind streaks in radar data. Because this is a preliminary evaluation of the recently acquired data we will focus on one scene and, thus, only on the effects of wavelength and polarization. Wind streaks provide information on the near-surface prevailing winds and on the abundance of winderodible material, such as sand. The potential for a free-flyer radar system that could provide global radar images in multiple wavelengths, polarizations, and incidence angles requires definition of system parameters for mission planning. Furthermore, thousands of wind streaks were mapped from Magellan radar images of Venus; their interpretation requires an understanding of the interaction of radar with wind streaks and the surrounding terrain. Our experiment was conducted on wind streaks in the Altiplano of Bolivia to address these issues.

Description: The Thermal Infrared Multispectral Scanner (TIMS) was flown on the NASA C-130 aircraft for a series of 12 flights during HAPEX-Sahel at altitudes ranging from 0.25 to 6 km (0.6 to 15 m resolution). TIMS provides coverage of the 8 to 12 micrometer thermal infrared band in 6 contiguous channels. Thus it is possible to observe the spectral behavior of the surface emissivity over this wavelength interval.

Description: Thermal Infrared Multispectral Scanner (TIMS) data were acquired over the McDowell Mountains northeast of Scottsdale, Arizona during August 1994. The raw data were processed to emphasize lithologic differences using a decorrelation stretch and assigning bands 5, 3, and 1 to red, green, and blue, respectively. Processed data of alluvium flanking the mountains exhibit moderate color variation. The objective of this study was to determine, using a quantitative approach, what environmental variable(s), in the absence of bedrock, is/are responsible for influencing the spectral properties of the desert alluvial surface.

Description: This summary describes the spatial, spectral, and radiometric calibration of Thermal Infrared Multispectral Scanner (TIMS) performed at the Jet Propulsion Laboratory (JPL) Thermal Infrared Calibration Facility (TIRCAL) between May and August, 1994. The 1994 calibration of TIMS was the first to make use of the new EXABYTE (8mm helical-scan tape) recording system. With the new recorder, the TIMS data tapes may be read directly on any computer system that has an EXABYTE tape drive. We analyzed the calibration data sets using image processing procedures written in Interactive Data Language.

Description: In mineral exploration, the ability to distinguish and map petrochemical variations of magmatic rocks can be a useful reconnaissance tool. Alkalinity is one such petrochemical parameter and is used in the characterization of granitoid rocks. In quartz normative plutonic rocks, alkalinity is related to the composition and abundance of feldspars. Together with quartz abundance, knowledge of feldspar modes allows the classification of these igneous rocks according to the Streckeisen diagram. Alternative classification schemes rely on whole rock geochemistry instead of mineral identifications. The relative ease of obtaining whole rock analyses means that geochemical classifications tend to be favored in exploration geology. But the technique of thermal infrared spectroscopy of rocks yields information on mineralogy and is one that can be applied remotely. The goal of the current work then is to establish whether data from TIMS can be used to distinguish the mineralogical variations that relate to alkalinity. An ideal opportunity to test this thesis arises from the work presented in a paper by Dewitt (1989). This paper contains the results of mapping and analysis of Proterozoic plutonic rocks in north-central Arizona. The map resulting from this work delineates plutons according to alkalinity in an effort to establish a trend or polarity in the regional magmatism. Also contained within this paper are brief descriptions of the mineralogy of half of the region's plutons. This combination of mineralogical and geochemical information was the rationale behind choosing this area as a site for TIMS over flights. A portion of the region centered on the northern Bradshaw Mountains was selected because it contains plutons of all three alkalinity classifications (alkali-calcic, calc-alkalic, and calic) present on DeWitt's map within a relatively small area. The site was flown in August of 1994 and the data received a few days before the writing of this manuscript. Most of this paper is devoted to the description of laboratory based spectroscopy and spectral simulations. These are required to gain insight into the correct procedures for enhancing the relatively small differences in the low spectral resolution TIMS data.

Description: This summary describes the 9 March 1994 Thermal Infrared Multispectral Scanner (TIMS) airborne calibration experiment conducted at Castaic Lake, California. This experiment was a collaborative effort between the TIMS and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) science teams at the Jet Propulsion Laboratory (JPL). TIMS was flown on the NASA/Ames Research Center C130 with the new retractable air fence installed in the TIMS instrument bay. The purpose of this experiment was to determine if the fence would reduce the air turbulence in the TIMS instrument bay, thereby reducing the errors in calibration caused by wind-blast cooling of the blackbody reference sources internal to TIMS. Previous experiments have indicated that the wind blast effect could cause TIMS to over-estimate surface temperatures by more than 10 C. We have examined the TIMS data from twelve lines flown over Castaic Lake. Four of the lines were flown at an altitude of approximately 2.5 km (MSL), four at an altitude of approximately 6.7 km, and four at approximately 8.3 km. At each altitude there were flights with northern and southern headings, with the aircraft level and at a positive pitch (nose-up attitude). The suite of twelve flights was designed to subject the TIMS/air fence system to different wind conditions and air temperatures. The TIMS flights were supported by a ground-truth team, who measured lake surface temperatures from a boat, and an atmosphere characterization team, who launched an airsonde and measured solar irradiance with a Reagan Sun Photometer. The Reagan measurements were used to construct a time-series of estimates of the total abundance of water vapor in the atmospheric column. These estimates were used to constrain modifications of the airsonde water vapor profile measurements made when processing the TIMS data with a customized version of the MODTRAN radiative transfer code.

Description: The sedimentary geology of the western Grand Canyon consists of gently northeast dipping sandstones, shales, and carbonates. However, due to facies changes within the units, the geomorphology varies from that seen by visitors at the National Park Headquarters. There, the cliff and slope expression of the rocks is replaced in the west by a series of mesas, ridges, and horizontal platforms. The largest of these occurs on the Esplanade Sandstone within the Supai Formation. The Esplanade is formed by slope retreat of the overlying units and resistance to erosion by the underlying limestones. It is onto this platform that the lavas of the Uinkaret Plateau were emplaced. The Uinkaret lava field lies 120 km south of St. George, Utah and is tectonically defined by two major normal faults -- the Hurricane to the west and the Toroweap to the east. The purpose of this investigation was to collect visible, near and thermal infrared data at different periods of the day and year. It is expected that these data will provide the ability to retrieve water temperatures; monitor sediment loads; map and examine any changes in the near shore vegetation communities and understand some of the intricacies of the geology. This paper will serve, to some degree, as a progress report on the Grand Canyon study, since only a fraction of the data has been received and processed thus far. Data from the Thermal Infrared Multispectral Scanner (TIMS) and the Landsat Thematic Mapper simulator (NS001) were acquired simultaneously on April 4, 1994. A second data acquisition occurred on August 27, 1994. Initial analysis of the TIMS data indicates a remarkably noise-free data set with minimal atmospheric attenuation. environments is evident.

Description: This paper describes an empirical method to correct TIMS (Thermal Infrared Multispectral Scanner) data for atmospheric effects by transferring calibration from a laboratory thermal emission spectrometer to the TIMS multispectral image. The method does so by comparing the laboratory spectra of samples gathered in the field with TIMS 6-point spectra for pixels at the location of field sampling sites. The transference of calibration also makes it possible to use spectra from the laboratory as endmembers in unmixing studies of TIMS data.