Results are presented from a 3-D Pluto general circulation model (GCM) that includes conductive heating and cooling, non-local thermodynamic equilibrium (non-LTE) heating by methane at 2.3 and 3.3 microns, non-LTE cooling by cooling by methane at 7.6 microns, and LTE CO rotational line cooling. The GCM also includes a treatment of the subsurface temperature and surface-atmosphere mass exchange. An initially 1 m thick layer of surface nitrogen frost was assumed such that it was large enough to act as a large heat sink (compared with the solar heating term) but small enough that the water ice subsurface properties were also significant. Structure was found in all three directions of the 3-D wind field (with a maximum magnitude of order 10 m/s in the horizontal directions and 10$^-5$ microbar/s in the vertical direction). Prograde jets were found at several altitudes. The direction of flow over the poles was found to very with altitude. Broad regions of up-welling and down-welling were also found. Predictions of vertical temperature profiles are provided for the Alice and REX instruments on New Horizons, while predictions of light curves are provided for ground-based stellar occultation observations. With this model methane concentrations of 0.2% and 1.0% and 8 and 24 microbar surface pressures are distinguishable. For ground-based stellar occultations, a detectable difference exists between light curves with the different methane concentrations, but not for different initial global mean surface pressures.

We present a 3D general circulation model of Pluto and Triton's atmospheres, which uses radiative-conductive-convective forcing. In both the Pluto and Triton models, an easterly (prograde) jet is present at the equator with a maximum magnitude of 10-12 m/s and 4 m/s, respectively. Neither atmosphere shows any significant overturning circulation in the meridional and vertical directions. Rather, it is horizontal motions (mean circulation and transient waves) that transport heat meridionally at a magnitude of 1 and 3 x 10^7 W at Pluto's autumn equinox and winter solstice, respectively (seasons referenced to the Northern Hemisphere). The meridional and dayside-nightside temperature contrast is small (<5 K). We find that the lack of vertical motion can be explained on Pluto by the strong temperature inversion in the lower atmosphere. The height of the Voyager 2 plumes on Triton can be explained by the dynamical properties of the lower atmosphere alone (i.e., strong wind shear) and does not require a thermally defined troposphere (i.e., temperature decreasing with height at the surface underlying a region of temperature increasing with height). The model results are compared with Pluto stellar occultation light curve data from 1988, 2002, 2006, and 2007 and Triton light curve data from 1997.

We use the Massachusetts Institute of Technology general circulation model (GCM) dynamical core, in conjunction with a Newtonian relaxation scheme that relaxes to a gray, analytical solution of the radiative transfer equation, to simulate a tidally locked, synchronously orbiting super-Earth exoplanet. This hypothetical exoplanet is simulated under the following main assumptions: (1) the size, mass, and orbital characteristics of GJ 1214b (Charbonneau et al., 2009), (2) a greenhouse-gas dominated atmosphere, (3), the gas properties of water vapor, and (4) a surface. We have performed a parameter sweep over global mean surface pressure (0.1, 1, 10, and 100 bar) and global mean surface albedo (0.1, 0.4, and 0.7). Given assumption (1) above, the period of rotation of this exoplanet is 1.58 Earth-days, which we classify as the rapidly rotating regime. Our parameter sweep differs from Heng and Vogt (2011), who performed their study in the slowly rotating regime and using Held and Suarez (1994) thermal forcing. This type of thermal forcing is a prescribed function, not related to any radiative transfer, used to benchmark Earth's atmosphere. An equatorial, westerly, superrotating jet is a robust feature in our GCM results. This equatorial jet is westerly at all longitudes. At high latitudes, the flow is easterly. The zonal winds do show a change with global mean surface pressure. As global mean surface pressure increases, the speed of the equatorial jet decreases between 9 and 15 hours local time (substellar point is located at 12 hours local time). The latitudinal extent of the equatorial jet increases on the nightside. Furthermore, the zonal wind speed in the equatorial and midlatitude jets decreases with increasing surface albedo. Also, the latitudinal width of the equatorial jet decreases as surface albedo increases.

We observed a stellar occultation by Titan on 2003 November 14 from La Palma Observatory using ULTRACAM with three Sloan filters: u', g', and i' (358, 487, and 758 nm, respectively). The occultation probed latitudes 2 degrees S and 1 degrees N during immersion and emersion, respectively. A prominent central flash was present in only the i' filter, indicating wavelength-dependent atmospheric extinction. We inverted the light curves to obtain six lower-limit temperature profiles between 335 and 485 km (0.04 and 0.003 mb) altitude. The i' profiles agreed with the temperature measured by the Huygens Atmospheric Structure Instrument [Fulchignoni, M. et al., 2005. Nature 438, 785-791] above 415 km (0.01 mb). The profiles obtained from different wavelength filters systematically diverge as altitude decreases, which implies significant extinction in the light curves. Applying an extinction model [Elliot, J.L., Young, L.A., 1992. Astron. J. 103, 991-1015] gave the altitudes of line of sight optical depth equal to unity: 396 +/- 7 km and 401 +/- 20 km (u' immersion and emersion); 354 +/- 7 km and 387 +/- 7 km (g' immersion and emersion); and 336 +/- 5 km and 318 +/- 4 km (i' immersion and emersion). Further analysis showed that the optical depth follows a power law in wavelength with index 1.3 +/- 0.2. We present a new method for determining temperature from scintillation spikes in the occulting body's atmosphere. Temperatures derived with this method are equal to or warmer than those measured by the Huygens Atmospheric Structure Instrument. Using the highly structured, three-peaked central flash, we confirmed the shape of Titan's middle atmosphere using a model originally derived for a previous Titan occultation [Hubbard, W.B. et al., 1993. Astron. Astrophys. 269, 541-563].

We present millimeter interferometry images of the CO J=1-0 line emission arising in the circumstellar envelope of HD 56126 (a.k.a. IRAS 07134+1005), which is one of the best studied 21-micron proto-planetary nebulae (PPNs). The CO emission extends from 1.2" to 7" in radius from the central star and appears consistent with a simple expanding envelope, as expected for a post-AGB star. We quantitatively model the molecular envelope using a radiative transfer code that we have modified for detached shells. Our best fit model reveals that two sequential winds created the circumstellar envelope of HD 56126: an AGB wind that lasted 6500 years with a mass-loss rate of 5.1x10^{-6} M_odot yr^{-1} and a more intense superwind that lasted 840 years with a mass-loss rate of 3x10^{-5} M_odot yr^{-1} and that ended the star's life on the AGB 1240 years ago. The total mass of this envelope is 0.059 M_odot which indicates a lower limit progenitor mass for the system of 0.66 M_odot, quite reasonable for this low-metallicity star which probably resides in the thick disk of the Galaxy. Comparison with images of the dust emission reveal a similar structure with the gas in the inner regions. Using 2-DUST, we model the dust emission of this source so that the model is consistent with the CO emission model and find an average gas-to-dust mass ratio of 75.