Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATMOSPHERES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

Model study of atmospheric transport using carbon 14 and strontium 90 as inert tracers

The observed excess carbon 14 in the atmosphere from 1963 to 1970 provides unique, but limited, data up to an altitude of about 35 km for testing the air motions calculated by 11 multidimensional atmospheric models. Strontium 90 measurements in the atmosphere from 1964 to mid-1967 provide data that have more latitude coverage than those of carbon 14 and are useful for testing combined models of air motions and aerosol settling. Model calculations for carbon 14 begin at October 1963, 9 months after the conclusion of the nuclear bomb tests; the initial conditions for the calculations are derived by three methods, each of which agrees fairly well with measured carbon 14 in October 1963 and each of which has widely different values in regions of the stratosphere where there were no carbon 14 measurements. The model results are compared to the stratospheric measurements, not as if the observed data were absolute standards, but in an effort to obtain new insight about the models and about the atmosphere. The measured carbon 14 vertical profiles at 31 deg N are qualitatively different from all of the models; the measured vertical profiles show a maximum mixing ratio in the altitude range of 20 to 25 km from October 1963 through July 1966, but all modeled profiles show mixing ratio maxima that increase in altitude from 20 km in October 1963 to greater than 40 km by April 1966. Both carbon 14 and strontium 90 data indicate that the models differ substantially among themselves with respect to stratosphere-troposphere exchange rate, but the modeled carbon 14 stratospheric residence times indicate that differences among the models are small with respect to transport rate between the middle stratosphere and the lower stratosphere. Strontium 90 data indicate that aerosol settling is important up to at least 35 km altitude. Relative to the measurements, about three quarters of the models transport carbon 14 from the lower stratosphere to the troposphere too rapidly, and all models with all three sets of initial conditions appear to sweep carbon 14 out of the midstratosphere (above 28 km) much more slowly than was observed from 1963 to 1970. Atmospheric dynamicists should seriously consider that the carbon 14 measurements imply with respect to advection in the midstratosphere.

Kinnison, D. E.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

Calculation of surface and top of atmosphere radiative fluxes from physical quantities based on ISCCP data sets. 1: Method and sensitivity to input data uncertainties

The largest uncertainty in upwelling shortwave (SW) fluxes (approximately equal 10-15 W/m(exp 2), regional daily mean) is caused by uncertainties in land surface albedo, whereas the largest uncertainty in downwelling SW at the surface (approximately equal 5-10 W/m(exp 2), regional daily mean) is related to cloud detection errors. The uncertainty of upwelling longwave (LW) fluxes (approximately 10-20 W/m(exp 2), regional daily mean) depends on the accuracy of the surface temperature for the surface LW fluxes and the atmospheric temperature for the top of atmosphere LW fluxes. The dominant source of uncertainty is downwelling LW fluxes at the surface (approximately equal 10-15 W/m(exp 2)) is uncertainty in atmospheric temperature and, secondarily, atmospheric humidity; clouds play little role except in the polar regions. The uncertainties of the individual flux components and the total net fluxes are largest over land (15-20 W/m(exp 2)) because of uncertainties in surface albedo (especially its spectral dependence) and surface temperature and emissivity (including its spectral dependence). Clouds are the most important modulator of the SW fluxes, but over land areas, uncertainties in net SW at the surface depend almost as much on uncertainties in surface albedo. Although atmospheric and surface temperature variations cause larger LW flux variations, the most notable feature of the net LW fluxes is the changing relative importance of clouds and water vapor with latitude. Uncertainty in individual flux values is dominated by sampling effects because of large natrual variations, but uncertainty in monthly mean fluxes is dominated by bias errors in the input quantities.

Zhang, Y.-C.↗

Energy-dependent effects of scattering atmospheres on X-ray pulsar pulse profiles

We propose that radiation-supported scattering atmospheres near accreting X-ray pulsars (XRPs) can explain energy-dependent features observed in the pulse profiles of 4U 1626-67, 4U 1538-52, 4U 1907+09 and Vela X-1. These atmospheres provide a physical model for the phenomenological annular emitting regions employed by Leahy to fit X-ray pulsar pulse profiles. We examine the effects of the scattering atmospheres under the assumptions that stable, optically thick atmospheres exist in a region where the optically thin resonant radiation force exceeds the force of gravity on ionized hydrogen. We predict that less complex pulse profiles will be observed at higher photon energies because the scattering atmospheres, which are supported by resonant Compton radiation pressure, become transparent to photons with energies greater than the cyclotron energy at the neutron star surface.

Sturner, Steven J.↗

Lead-induced atmospheric circulations

Leads provide a significant source of heat and moisture to the Arctic winter atmosphere, and plumes from wide leads have been observed to penetrate the Arctic inversion. We have developed a two-dimensional, high-resolution, deep anelastic numerical model to investigate the atmospheric convection from leads with widths ranging from 100 m to 10 km. A second-order turbulence closure scheme is used to parameterize the atmospheric turbulence in the horizontally inhomogeneous system. This study describes how the lead-induced circulations can enhance the vertical transport of heat into the atmospheric boundary layer. This model is compared with large-eddy simulation results and with lidar observations of a lead-induced ice crystal plume. The model is used to study the effect of varying lead widths and ambient atmospheric conditions on the resultant convection from leads, and some preliminary results are described.

Alam, Afshan↗

An analytic radiative transfer model for a coupled atmosphere and leaf canopy

A new analytical radiative transfer model of a leaf canopy is developed that approximates multiple-scattering radiance by a four-stream formulation. The canopy model is coupled to a homogeneous atmospheric model as well as a non-Lambertian lower boundary soil surface. The same four-stream formulation is also used for the calculation of multiple scattering in the atmosphere. Comparisons of radiance derived from the four-stream model with those calculated by an iterative numerical solution of the radiative transfer equation show that the analytic model has a very high accuracy, even with a turbid atmosphere and a very dense canopy in which multiple scattering dominates. Because the coupling of canopy and atmospheric models fully accommodates anisotropic surface reflectance and atmospheric scattering and its effect on directional radiance, the model is especially useful for application to directional radiance and measurements obtained by remote sensing. Retrieval of biophysical parameters using this model is under investigation.

Liang, Shunlin↗

Ocean color and atmospheric dimethyl sulfide: On their mesoscale variability

The mesoscale variability of dimethyl sulfide (DMS) and ocean color is explored to determine the feasibility of a predictive relationship. During NASA's Global Tropospheric Experiment/Chemical Instrumentation Test and Evaluation (GTE/CITE 3), simultaneous shipboard and aircraft studies were carried out in the North Atlantic, followed by aircraft studies in the South Atlantic. Surface concentrations of chlorophyll alpha were measured with an airborne spectroradiometer, the Ocean Data Acquisition System (ODAS), with simultaneous determinations of tropospheric DMS. Shipboard measurements of DMS in air and water as well as in situ chlorophyll alpha were taken in the North Atlantic. No relation was observed between shipboard aquatic DMS and chlorophyll alpha or primary productivity. Higher levels of aqueous DMS were not always reflected by atmospheric DMS, although shipboard and aircraft measurements of atmospheric DMS agreed very well. A significant relationship between atmospheric DMS and ocean color was seen once at low altitudes in both the North and South Atlantic only under clean air conditions. Atmospheric DMS levels during the North Atlantic experiment were probably lowered by the presence of mostly polluted air masses in the study area and were, overall, probably not representative of the in situ sea-to-air flux of DMS. Changes in concentration of aircraft-sensed chlorophyllous pigments were not reflected by atmospheric DMS. If a predictive algorithm is to be found, phytoplankton blooms should probably be the first place to study an ocean color-DMS relationship.

Matrai, Patricia A.↗

Viking Lander image analysis of Martian atmospheric dust

We have reanalyzed three sets of Viking Lander 1 and 2 (VL1 and VL2) images of the Martian atmosphere to better evaluate the radiative properties of the atmospheric dust particles. The properties of interest are the first two moments of the size distribution, the single-scattering albedo, the dust single-scattering phase function, and the imaginary index of refraction. These properties provide a good definition of the influence that the atmospheric dust has on heating of the atmosphere. Our analysis represents a significant improvement over past analyses (Pollack et al. 1977, 1979) by deriving more accurate brightness closer to the sun, by carrying out more precise analyses of the data to acquire the quantities of interest, and by using a better representation of scattering by nonspherical particles. The improvements allow us to better define the diffraction peak and hence the size distribution of the particles. For a lognormal particle size distribution, the first two moments of the size distribution, weighted by the geometric cross section, are found. The geometric cross-section weighted mean radius r(sub eff) is found to be 1.85 +/- 0.3 micrometers at VL2 during northern summer when dust loading was low and 1.52 +/- 0.3 micrometers at VL1 during the first dust storm. In both cases the best cross-section weighted mean variance nu(sub eff) of the size distribution is equal to 0.5 +/- 0.2 micrometers. The changes in size distribution, and thus radiative properties, do not represent a substantial change in solar energy deposition in the atmosphere over the Pollak et al. (1977, 1979) estimates.

Pollack, James B.↗

Diagnostic calculations of the circulation in the Martian atmosphere

The circulation of the Martian atmosphere during late southern summer is derived from atmospheric temperature and dust distributions retrieved from a subset of the Mariner 9 infrared interferometer spectrometer (IRIS) thermal emission spectra (LS = 343-348 deg) (Santee and Crisp, 1933). Zonal-mean zonal winds are calculated by assuming gradient wind balance and zero surface zonal wind. Both hemispheres have intense midlatitude westerly jets with velocities of 80-90 m/s near 50 km; in the southern tropics the winds are easterly with velocities of 40 m/s near 50 km. The net effect of the zonal-mean meridional circulation and large-scale waves can be approximated by the diabatic ciculation, which is defined from the atmospheric thermal structure and net radiative heating rates. The radiative transfer model described by Crisp (1990) and Santee (1993) is used to compute solar heating and thermal cooling rates from diurnal averages of the retrieved IRIS temperature and dust distributions. At pressures below 4 mbar, there are large net radiative heating rates (up to 5 K/d) in the equatorial region and large net radiative cooling rates (up to 12 K/d) in the polar regions. These net radiative heating rates are used in a diagnostic stream function model which solves for the meridonal and vertical components of the diabatic circulation simultaneously. We find a two cell circulation, with rising motion over the equator, poleward flow in both hemispheres, sinking motion over both polar regions, and return flow in the lowest atmospheric levels. The maximum poleward velocity is 3 m/s in the tropics at approximately 55 km altitude, and the maximum vertical velocity is 2.5 cm/s downward over the north pole at approximately 60 km altitude. If these large transport rates are sustained for an entire season, the Martian atmosphere above the 1-mbar level is overturned in about 38 days. This diabatic circulation is qualitatively similar to the terrestial diabatic circulation at the comparable season, but is more vigorous.

Santee, Michelle L.↗

Evolution of the martian atmosphere

Carbon dioxide, nitrogen, and the nonradiogenic and radiogenic noble gases are tracked from primordial inventories to their present states in a revised model of atmospheric evolution on Mars. Elemental and isotopic abundances evolve by hydrodynamic escape, impact erosion, outgassing, sputtering, photochemical escape of nitrogen, and carbonate formation and recycling. Atmospheric history is divided into early and late evolutionary periods, the first characterized by high CO2 pressures and a possible greenhouse and the second by a low pressure cap-regolith buffered system initiated by polar CO2 condensation approximately 3.7 Gyr ago. During early evolution the Xe isotopes are fractionated to their present composition by hydraulic escape, and CO2 pressure and isotopic history are dictated by the interplay of losses to erosion, sputtering, and carbonate precipitation, additions by outgassing and carbonate recycling, and perhaps also by feedback stabilization under greenhouse conditions. Atmospheric collapse near 3.7 Gyr leads to abrupt increases in the mixing ratios of preexisting Ar, Ne, and N2 at the exobase and their rapid removal by sputtering. Current abundances and isotopic compositions of these light species are therefore entirely determined by the action of sputtering and photochemical escape on gases supplied by outgassing during the late evolutionary epoch. The present atmospheric Kr inventory also derives almost completely from solar-like Kr degassed during this period. Consequently, among current observables, only the Xe isotopes and delta(C-13) survive as isotopic tracers of atmospheric history prior to its transition to low pressure.

Pepin, Robert O.↗

Quasi-periodic atmosphere-regolith-cap CO2 redistribution in the martian past

Our earlier Mars regolith-atmosphere-cap CO2 distribution model (Fanale et al., 1982, Icarus 50, 381-407) has been improved, revised, and extended back over Mars' mid to late history. The present model takes into account four new factors: (1) a more realistic long-term obliquity cycle, (2) thermal conduction as it affects the surface energy balance, (3) the changing solar constant, and (4) atmospheric erosion 3.5 byr ago to the present. Solar insolation and temperatures are computed for the full range of obliquities, latitudes, and epochs, and a CO2 adsorption relation is used, together with a conservation of mass constraint, to calculate atmospheric pressures and exchangeable CO2 mass as functions of obliquity and epoch for the regolith, atmosphere, and polar caps for two assumed thicknesses of a basalt regolith. It is found that the heat conduction term in the surface boundary condition has an important effect in reducing the range of atmospheric pressures over the obliquity cycle at all epochs.

Fanale, Fraser P.↗

Response of the upper atmosphere to variations in the solar soft x-ray irradiance

Terrestrial far ultraviolet (FUV) airglow emissions have been suggested as a means for remote sensing the structure of the upper atmosphere. The energy which leads to the excitation of FUV airglow emissions is solar irradiance at extreme ultraviolet (EUV) and soft x-ray wavelengths. Solar irradiance at these wavelengths is known to be highly variable; studies of nitric oxide (NO) in the lower thermosphere have suggested a variability of more than an order of magnitude in the solar soft x-ray irradiance. To properly interpret the FUV airflow, the magnitude of the solar energy deposition must be known. Previous analyses have used the electron impact excited Lyman-Birge-Hopfield (LBH) bands of N2 to infer the flux of photoelectrons in the atmosphere and thus to infer the magnitude of the solar irradiance. This dissertation presents the first simultaneous measurements of the FUV airglow, the major atmospheric constituent densities, and the solar EUV and soft x-ray irradiances. The measurements were made on three flights of an identical sounding rocket payload at different levels of solar activity. The linear response in brightness of the LBH bands to variations in solar irradiance is demonstrated. In addition to the N2 LBH bands, atomic oxygen lines at 135.6 and 130.4 nm are also studied. Unlike the LBH bands, these emissions undergo radiative transfer effects in the atmosphere. The OI emission at 135.6 nm is found to be well modeled using a radiative transfer calculation and the known excitation processes. Unfortunately, the assumed processes leading to OI 130.4 nm excitation are found to be insufficient to reproduce the observed variability of this emission. Production of NO in the atmosphere is examined; it is shown that a lower than previously reported variability in the solar soft x-ray irradiance is required to explain the variability of NO.

Bailey, Scott Martin↗

Research concerning the net flux of radiation in the atmosphere of Jupiter

The plan of the NFR (Net Flux of Radiation) team is for the data from the two solar channels (B and E) of NFR to be reduced with the goal of determining the solar heating rate. In order to determine the solar heating rate from the NFR measurements, effects due to the instrument's spatial and spectral response functions, to the temperature variation of the instrument (and associated drift of calibration), to the setting sun, and to the rotation of the probe (initially at a rate comparable to the NFR sampling frequency), all must be well modelled. In the past year, a forward modeling routine was created to simulate NFR data return in the B and E channels. The effects of varying parameters describing the atmospheric model (such as cloud location and thickness) and the descent profile (such as rotation rate) were investigated and an inversion routine was developed. For the forward modeling, existing radiative transfer codes were used to determine intensity fields within the Jovian atmosphere. A routine was developed to determine instantaneous instrument response by integrating the intensity field over the instrument response functions. A second routine was developed to determine the actual output of the NFR by integrating along an arbitrary descent trajectory. Near the top of the atmosphere, the upflux data alone are used to constrain the cloud structure of he atmosphere. To accomplish this, models are used to describe the variation in up flux between consecutive measurements in terms of variations of cloud opacity and variations in known parameters such as the solar zenith angle. This allows us to develop a zero-order model of cloud structure. Lower in the atmosphere, at levels where there is little or no azimuthal structure to the net flux measurements, both the up flux and net flux are used to derive layer transmission and reflection functions, which then determine layer opacity and single scattering albedo. A preliminary analysis of the data began in December 1995. In these data we could see the rapid oscillations expected at the beginning of the data due to probe rotation and the sun passing through the edge of the field of view. In addition, the time when this oscillation stopped was clearly visible. This sets the rough optical depth above the probe at this time.

Tomasko, M. G.↗

Composition/Structure/Dynamics of comet and planetary satellite atmospheres

This research program addresses two cases of tenuous planetary atmospheres: comets and Io. The comet atmospheric research seeks to analyze a set of spatial profiles of CN in comet Halley taken in a 7.4-day period in April 1986; to apply a new dust coma model to various observations; and to analyze observations of the inner hydrogen coma, which can be optically thick to the resonance scattering of Lyman-alpha radiation, with the newly developed approach that combines a spherical radiative transfer model with our Monte Carlo H coma model. The Io research seeks to understand the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and plasma given methods and algorithms developed for the study of neutral gas cometary atmospheres and the earth's polar wind and plasmasphere. Progress is reported on cometary Hydrogen Lyman-alpha studies; time-series analysis of cometary spatial profiles; model analysis of the dust comae of comets; and a global kinetic atmospheric model of Io.

Combi, Michael R.↗

Measuring the spectral emissivity of thermal protection materials during atmospheric reentry simulation

Hypersonic spacecraft reentering the earth's atmosphere encounter extreme heat due to atmospheric friction. Thermal Protection System (TPS) materials shield the craft from this searing heat, which can reach temperatures of 2900 F. Various thermophysical and optical properties of TPS materials are tested at the Johnson Space Center Atmospheric Reentry Materials and Structures Evaluation Facility, which has the capability to simulate critical environmental conditions associated with entry into the earth's atmosphere. Emissivity is an optical property that determines how well a material will reradiate incident heat back into the atmosphere upon reentry, thus protecting the spacecraft from the intense frictional heat. This report describes a method of measuring TPS emissivities using the SR5000 Scanning Spectroradiometer, and includes system characteristics, sample data, and operational procedures developed for arc-jet applications.

Marble, Elizabeth↗

Sulfur, Chlorine, and Flourine Degassing and Atmospheric Loading by the 1783 - 1784 AD Laki (Skaftar Fires) Eruption in Iceland

The 1783-1784 Laki tholeiitic basalt fissure eruption in Iceland was one of the greatest atmospheric pollution events of the past 250 years, with widespread effects in the northern hemisphere. The degassing history and volatile budget of this event are determined by measurements of pre-eruption and residual contents of sulfur, chlorine, and fluorine in the products of all phases of the eruption. In fissure eruptions such as Laki, degassing occurs in two stages: by explosive activity or lava fountaining at the vents, and from the lava as it flows away from the vents. Using the measured sulfur concentrations in glass inclusions in phenocrysts and in groundmass glasses of quenched eruption products, we calculate that the total accumulative atmospheric mass loading of sulfur dioxide was 122 Mt over a period of 8 months. This volatile release is sufficient to have generated approximately 250 Mt of H2SO4 aerosols, an amount which agrees with an independent estimate of the Laki aerosol yield based on atmospheric turbidity measurements. Most of this volatile mass (approximately 60 wt.%) was released during the first 1.5 months of activity. The measured chlorine and fluorine concentrations in the samples indicate that the atmospheric loading of hydrochloric acid and hydrofluoric acid was approximately 7.0 and 15.0 Mt, respectively. Furthermore, approximately 75% of the volatile mass dissolved by the Laki magma was released at the vents and carried by eruption columns to altitudes between 6 and 13 km. The high degree of degassing at the vents is attributed to development of a separated two-phase flow in the upper magma conduit, and implies that high-discharge basaltic eruptions such as Laki are able to loft huge quantities of gas to altitudes where the resulting aerosols can reside for months, or even 1-2 years. The atmospheric volatile contribution due to subsequent degassing of the Laki lava flow is only 18 wt.% of the total dissolved in the magma, and these emissions were confined to the lowest regions of the troposhere and therefore important only over Iceland. This study indicates that determination of the amount of sulfur degassed from the Laki magma batch by measurements of sulfur in the volcanic products (the petrologic method) yields a result which is sufficient to account for the mass of aerosols estimated by other methods.

Thordarson, T.↗

Viking Lander image analysis of Martian atmospheric dust

We have reanalyzed three sets of Viking Lander 1 and 2 (VL1 and VL2) images of the Martian atmosphere to better evaluate the radiative properties of the atmospheric dust particles. The properties of interest are the first two moments of the size distribution, the single-scattering albedo, the dust single-scattering phase function, and the imaginary index of refraction. These properties provide a good definition of the influence that the atmospheric dust has on heating of the atmosphere. Our analysis represents a significant improvement over past analyses (Pollack et al. 1977,1979) by deriving more accurate brightnesses closer to the sun, by carrying out more precise analyses of the data to acquire the quantities of interest, and by using a better representation of scattering by nonspherical particles. The improvements allow us to better define the diffraction peak and hence the size distribution of the particles. For a lognormal particle size distribution, the first two moments of the size distribution, weighted by the geometric cross section, are found. The geometric cross-section weighted mean radius (r(sub eff)) is found to be 1.85 +/- 0.3 microns at VL2 during northern summer when dust loading was low and 1.52 +/- 0.3 microns at VL1 during the first dust storm. In both cases the best cross-section weighted mean variance (nu(eff)) of the size distribution is equal to 0.5 +/- 0.2 microns. The changes in size distribution, and thus radiative properties, do not represent a substantial change in solar energy deposition in the atmosphere over the Pollack et al. (1977,1979) estimates.

Pollack, James B.↗

Diagnostic calculations of the circulation in the Martian atmosphere

The circulation of the Martian atmosphere during late southern summer is derived from atmospheric temperature and dust distributions retrieved from a subset of the Mariner 9 infrared interferometer spectrometer (IRIS) thermal emission spectra (L(sub s) = 343-348 deg). Zonal-mean zonal winds are calculated by assuming gradient wind balance and zero surface zonal wind. Both hemispheres have intense midlatitude westerly jets with velocities of 80-90 m/s near 50 km; in the southern tropics the winds are easterly with velocities of 40 m/s near 50 km. The net effect of the zonal mean meridional circulation and large-scale waves can be approximated by the diabatic circulation, which is defined from the atmospheric thermal structure and net radiative heating rates. The radiative transfer model described by Crisp (1990) and Santee (1993) is used to compute solar heating and thermal cooling rates from diurnal averages of the retrieved IRIS temperature and dust distributions. At pressures below 4 mbar, there are large net radiative heating rates (up to 5 K/d) in the equatorial region and large net radiative cooling rates (up to 12 K/d) in the polar regions. These net radiative heating rates are used in a diagnostic stream function model which solves for the meridional and vertical components of the diabatic circulation simultaneously. We find a two-cell circulation, with rising motion over the equator, poleward flow in both hemispheres, sinking motion over both polar regions, and return flow in the lowest atmospheric levels. The maximum poleward velocity is 3 m/s in the tropics at approx. 55 km altitude, and the maximum vertical velocity is 2.5 cm/s downward over the north pole at approx. 60 km altitude. If these large transport rates are sustained for an entire season, the Martian atmosphere above the 1-mbar level is overturned in about 38 days. This diabatic circulation is qualitatively similar to the terrestrial diabatic circulation at the comparable season, but is more vigorous.

Santee, Michelle L.↗