Thick and thin models of the evolution of carbon dioxide on Mars

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.     

Seasonal variations in prey selection by estuarine black-headed gulls (Larus ridibundus)

The number of black-headed gulls (Larus ridibundus) in the Clyde Estuary is large. In summer the average density has reached 1350 gulls km^?2 and in winter 180 gulls km^?2. This paper compares prey selection and feeding efficiency in gulls during summer and winter on tidal flats, and considers how seasonal differences may be adaptations to cope with seasonal changes in prey availability.Gross and net rates of energy intake were highest in summer because gulls captured more of the polychaete N. diversicolor than the amphipod C. volutator. In winter, gulls selected for C. volutator and therefore an energetically less profitable diet. Throughout the year gulls selected more C. volutator relative to N. diversicolor than expected on energetic grounds and so apparently did not maximize potential net rate of energy intake.Gulls used three techniques to capture prey and made most intensive use of the ‘crouch’ technique. Crouching gulls attained a much higher net rate of energy intake than ‘upright’ or ‘paddling’ gulls.A log-linear model showed that (a) season, water depth and gull density determined feeding technique and (b) feeding technique and season independently determined foraging success and prey selection. Thus gull density and water depth acted on prey selection through imposed variations in feeding technique.Reasons for gulls selecting energetically unprofitable C. volutator and for the use of several distinct feeding techniques are discussed.     

Screening and ranking framework for geologic CO<Subscript>2</Subscript> storage site selection on the basis of health,safety, and environmental risk

A screening and ranking framework (SRF) has been developed to evaluate potential geologic carbon dioxide (CO₂) storage sites on the basis of health, safety, and environmental (HSE) risk arising from CO₂ leakage. The approach is based on the assumption that CO₂ leakage risk is dependent on three basic characteristics of a geologic CO₂ storage site: (1) the potential for primary containment by the target formation; (2) the potential for secondary containment if the primary formation leaks; and (3) the potential for attenuation and dispersion of leaking CO₂ if the primary formation leaks and secondary containment fails. The framework is implemented in a spreadsheet in which users enter numerical scores representing expert opinions or published information along with estimates of uncertainty. Applications to three sites in California demonstrate the approach. Refinements and extensions are possible through the use of more detailed data or model results in place of property proxies.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

Characterization and cycling of atmospheric mercury along the central US Gulf Coast

Concentrations of atmospheric Hg species, elemental Hg (Hg^°), reactive gaseous Hg (RGM), and fine particulate Hg (Hg-PM_2.5) were measured at a coastal site near Weeks Bay, Alabama from April to August, 2005 and January to May, 2006. Mean concentrations of the species were 1.6 ± 0.3 ng m⁻³, 4.0 ± 7.5 pg m⁻³ and 2.7 ± 3.4 pg m⁻³, respectively. A strong diel pattern was observed for RGM (midday maximum concentrations were up to 92.7 pg m⁻³), but not for Hg^° or Hg-PM_2.5. Elevated RGM concentrations (>25 pg m⁻³) in April and May of 2005 correlated with elevated average daytime O₃ concentrations (>55 ppbv) and high light intensity (>500 W m⁻²). These conditions generally corresponded with mixed continental-Gulf and exclusively continental air mass trajectories. Generally lower, but still elevated, RGM peaks observed in August, 2005 and January–March, 2006 correlated significantly (p < 0.05) with peaks in SO₂ concentration and corresponded to periods of high light intensity and lower average daytime O₃ concentrations. During these times air masses were dominated by trajectories that originated over the continent. Elevated RGM concentrations likely resulted from photochemical oxidation of Hg^° by atmospheric oxidants. This process may have been enhanced in and by the near-shore environment relative to inland sites. The marine boundary layer itself was not found to be a significant source of RGM.     

Lagrangian diffusion equation and its application to oceanic dispersion

The Lagrangian diffusion equation appropriate for the dispersion of current followers (e. g., floats, drogues, drifters) is proposed. The analytical solution to the equation is obtained for a uniform deformation field, characterized by Lagrangian deformations and anisotropic eddy diffusivities both varying with time. Expressions are derived for the patch area and its elongation and rotation. For small values of elapsed time after the initial release the patch area can be accounted for by the exponential of the cumulative value of the horizontal divergence; the relative rate of change of the patch area can be accounted for by the horizontal divergence.     

Instrument calibration of ocean bottom seismographs

To increase the accuracy of measuring sea floor motion with ocean bottom seismometers, we calibrate the seismometer system on the ocean floor. Data from the sea floor calibration, augmented with electronic and land calibration data, enables us to find the OBS transfer function to an accuracy of 0.5% in the frequency range of 0.1 to 32 Hz. We are able to distinguish between temperature, instrument and OBS ground coupling effects, all of which alter the transfer function. This paper reviews our method of calibration and discusses the effects of temperature and some of the instrument design features on the vertical seismometer transfer function.