Climate sensitivity to changes in land surface characteristics

Using a recently developed global vegetation distribution, topography, and shorelines for the Early Eocene in conjunction with the Genesis version 2.0 climate model, we investigate the influences that these new boundary conditions have on global climate. Global mean climate changes little in response to the subtle changes we made; differences in mean annual and seasonal surface temperatures over northern and southern hemispheric land, respectively, are on the order of 0.5°C. In contrast, and perhaps more importantly, continental scale climate exhibits significant responses. Increased peak elevations and topographic detail result in larger amplitude planetary 4 mm/day and decreases by 7–9 mm/day in the proto Himalayan region. Surface temperatures change by up to 18°C as a direct result of elevation modifications. Increased leaf area index (LAI), as a result of altered vegetation distributions, reduces temperatures by up to 6°C. Decreasing the size of the Mississippi embayment decreases inland precipitation by 1–2 mm/day. These climate responses to increased accuracy in boundary conditions indicate that “improved” boundary conditions may play an important role in producing modeled paleoclimates that approach the proxy data more closely.     

Constraints on Mercury’s Na exosphere: Combined MESSENGER and ground-based data

We have used observations of sodium emission obtained with the McMath-Pierce solar telescope and MESSENGER’s Mercury Atmospheric and Surface Composition Spectrometer (MASCS) to constrain models of Mercury’s sodium exosphere. The distribution of sodium in Mercury’s exosphere during the period January 12-15, 2008, was mapped using the McMath-Pierce solar telescope with the 5″ × 5″ image slicer to observe the D-line emission. On January 14, 2008, the Ultraviolet and Visible Spectrometer (UVVS) channel on MASCS sampled the sodium in Mercury’s anti-sunward tail region. We find that the bound exosphere has an equivalent temperature of 900-1200 K, and that this temperature can be achieved if the sodium is ejected either by photon-stimulated desorption (PSD) with a 1200 K Maxwellian velocity distribution, or by thermal accommodation of a hotter source. We were not able to discriminate between the two assumed velocity distributions of the ejected particles for the PSD, but the velocity distributions require different values of the thermal accommodation coefficient and result in different upper limits on impact vaporization. We were able to place a strong constraint on the impact vaporization rate that results in the release of neutral Na atoms with an upper limit of 2.1 × 10⁶ cm⁻² s⁻¹. The variability of the week-long ground-based observations can be explained by variations in the sources, including both PSD and ion-enhanced PSD, as well as possible temporal enhancements in meteoroid vaporization. Knowledge of both dayside and anti-sunward tail morphologies and radiances are necessary to correctly deduce the exospheric source rates, processes, velocity distribution, and surface interaction.     

The effect of magnetic fields on the formation of circumstellar discs around young stars

We present first results of our simulations of magnetic fields in the formation of single and binary stars using a recently developed method for incorporating Magnetohydrodynamics (MHD) into the Smoothed Particle Hydrodynamics (SPH) method. An overview of the method is presented before discussing the effect of magnetic fields on the formation of circumstellar discs around young stars. We find that the presence of magnetic fields during the disc formation process can lead to significantly smaller and less massive discs which are much less prone to gravitational instability. Similarly in the case of binary star formation we find that magnetic fields, overall, suppress fragmentation. However these effects are found to be largely driven by magnetic pressure. The relative importance of magnetic tension is dependent on the orientation of the field with respect to the rotation axis, but can, with the right orientation, lead to a dilution of the magnetic pressure-driven suppression of fragmentation.     

Resonant Compton upscattering in anomalous X-ray pulsars

A significant new development in the study of Anomalous X-ray Pulsars (AXPs) has been the recent discovery by INTEGRAL and RXTE of flat, hard X-ray components in three AXPs. These non-thermal spectral components differ dramatically from the steeper quasi-power-law tails seen in the classic X-ray band in these sources. A prime candidate mechanism for generating this new component is resonant, magnetic Compton upscattering. This process is very efficient in the strong magnetic fields present in AXPs. Here an introductory exploration of an inner magnetospheric model for upscattering of surface thermal X-rays in AXPs is offered, preparing the way for an investigation of whether such resonant upscattering can explain the 20–150 keV spectra seen by INTEGRAL. Characteristically flat emission spectra produced by non-thermal electrons injected in the emission region are computed using collision integrals. A relativistic QED scattering cross section is employed so that Klein–Nishina reductions are influential in determining the photon spectra and fluxes. Spectral results depend strongly on the magnetospheric locale of the scattering and the observer’s orientation, which couple directly to the angular distributions of photons sampled.     

Lyman α radiative transfer in a multiphase medium

Hydrogen Lyman α (Lyα) is our primary emission-line window into high-redshift galaxies. Despite an extensive literature, Lyα radiative transfer in the most realistic case of a dusty, multiphase medium has received surprisingly little detailed theoretical attention. We investigate Lyα resonant scattering through an ensemble of dusty, moving, optically thick gas clumps. We treat each clump as a scattering particle and use Monte Carlo simulations of surface scattering to quantify continuum and Lyα surface scattering angles, absorption probabilities, and frequency redistribution, as a function of the gas dust content. This atomistic approach speeds up the simulations by many orders of magnitude, making possible calculations which are otherwise intractable. Our fitting formulae can be readily adapted for fast radiative transfer in numerical simulations. With these surface scattering results, we develop an analytic framework for estimating escape fractions and line widths as a function of gas geometry, motion, and dust content. Our simple analytic model shows good agreement with full Monte Carlo simulations. We show that the key geometric parameter is the average number of surface scatters for escape in the absence of absorption,     , and we provide fitting formulae for several geometries of astrophysical interest. We consider the following two interesting applications. (i) Equivalent widths ( EWs ). Lyα can preferentially escape from a dusty multiphase interstellar medium if most of the dust lies in cold neutral clouds, which Lyα photons cannot penetrate. This might explain the anomalously high EWs sometimes seen in high-redshift/submillimetre sources. (ii) Multiphase galactic outflows . We show the characteristic profile is asymmetric with a broad red tail, and relate the profile features to the outflow speed and gas geometry. Many future applications are envisaged.