The Combined Atmospheric Photochemistry and Ion Tracing code: Reproducing the Viking Lander results and initial outflow results

A Combined Atmospheric Photochemistry and Ion Tracing code (CAPIT) has been developed to explore ion loss into space at Mars. The CAPIT code includes the major photochemical reactions of Mars’ ionosphere, ion tracing in the presence of magnetic fields, and plasma wave heating of ions. In particular, we examine whether O⁺ escape from the day-side ionosphere is limited by ion production (UV input) or by external energy input to the ions. To verify the code, it is demonstrated that the CAPIT solutions reproduce the Viking 1 Lander’s ion density and temperature profiles. Using Viking 1 Lander conditions as a baseline, ion outflow rates are examined as function of solar wind energy input via plasma waves and UV ionization rates. The O⁺ outflow rates predicted by the simulation are comparable to the outflow rates estimated by observation. The results indicate that plasma waves are a viable source of energy to O⁺ ions and suggest that present-day O⁺ outflow rates at Mars are source limited by photochemical production (UV input) during periods of strong energy input (plasma wave activity), but otherwise regulated by both UV input and energy input. These results imply that ion heating by plasma waves can influence the present-day loss of O⁺.     

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).     

UVIS observations of the FUV OI and CO 4P Venus dayglow during the Cassini flyby

We analyze FUV spatially-resolved dayglow spectra obtained at 0.37 nm resolution by the UVIS instrument during the Cassini flyby of Venus. We use a least-squares fit method to determine the brightness of the OI emissions at 130.4 and OI 135.6 nm, and of the bands of the CO fourth positive system which are dominated by fluorescence scattering. We compare the brightness observed along the UVIS foot track of the two OI multiplets with that deduced from a model of the excitation of these emissions by photoelectron impact on O atoms and resonance scattering of the solar 130.4 nm emission. The large optical thickness 130.4 nm emission is accounted for using a radiative transfer model. The airglow intensities are calculated along the foot track and found to agree with the observed 130.4 nm brightness within ∼10%. The modeled OI 135.6 nm brightness is also well reproduced by the model. The oxygen density profile of the VTS3 model is found to be consistent with the observations. We find that self-absorption of the (0, v″) bands of the fourth positive emission of CO is important and we derive a CO vertical column of about 6.4 × 10¹⁵ cm⁻² in close agreement with the value provided by the VTS3 empirical atmospheric model.     

Tsunamis on Mars: Earth analogues of projected Martian sediment

Bolide impacts on Mars, within the proposed ocean boundaries (“contacts 1 and 2”) in the northern lowlands, would certainly have generated ultra high energy waves similar to tsunamis on Earth. Impacts into putative Noachian and Hesperian seas of variable areal extents and depths would have experienced high-energy inundations (transgressions), which would have left an imprint in the stack of deposits adjacent to the proposed shorelines. On Earth, the principal influencing factors for tsunami-wave energy are the character of shoreline topography and coastal water depth, which control wave compression and shoreline friction. Shorelines with narrow embayments and steep offshore gradients produce wave compression and increased collision of grains within the carried load contrasted with linear shorelines and shallow offshore gradients that dissipate energy. Steep offshore gradients produce concentrated major wave friction with the bed engendering high kinetic energy in the wave during emplacement of tsunami-generated sediment, which differs from shallow offshore beds that produce lower frictional effects over a wider area and drawdown of wave energy. Thus, overprinting of transported quartz grains on Earth is greatest where wave energy is highest, attenuated down to minor or nil overprinting where wave energy is less. Such grain overprinting in the form of energy-induced microtextures would also be observed in other grain types such as olivine and plagioclase, as such mineralogies are expected to dominate the Martian landscape based on orbital and local field (lander and rover) perspectives. Kinetic energy variation in tsunamis is controlled more by the square of velocity than mass, the resulting collisional effects of which produce swarms of v-shaped percussion microfeatures on quartz and other silicate mineral surfaces when velocity and compression are highest. This work indicates that a valid test for the ocean hypothesis is targeting “coastal” areas adjacent to narrow embayments where offshore depths are known to be highest, as possible tsunami-emplaced sediments, especially those that have been protected from atmospheric conditions through relatively rapid burial, may reveal a high frequency of percussion cracks, features of which appear to be unique to such terrestrial environments.     

New intrinsic-colour calibration for uvby-β photometry

A new intrinsic-colour calibration $((b-y{)}_o-\beta )$ is presented for the uvby-β photometric system, making use of re-calibrated Hipparcos parallaxes and published reddening maps. This new calibration for $(b-y{)}_o-\beta$, our Eq. (1), has been based upon stars with $d_{\mathit{Hip}}<70\text{pc}$ in the photometric catalogues of Schuster and Nissen (1988), Schuster et al., 2004, Schuster et al., 2006, provides a small dispersion, $\pm 0.009$, and has a positive “standard” $+2.239\mathrm{\Delta }\beta$ coefficient, which is not too different from the coefficients of Crawford (1975a, +1.11) and of Olsen (1988, +1.34). For 61 stars with spectra from CASPEC, UVES/VLT, and FIES/NOT databases, without detectable Na I lines, the average reddening value $〈E(b-y)〉=-0.001\pm 0.002$ shows that any zero-point correction to our intrinsic-colour equation must be minuscule.     

Simulated impacts of artificial groundwater recharge and discharge on the source area and source volume of an Atlantic Coastal Plain stream, Delaware, USA

A numerical groundwater-flow model was used to characterize the source area and volume of Phillips Branch, a baseflow-dominated stream incising a highly permeable unconfined aquifer on the low relief Delmarva Peninsula, USA. Particle-tracking analyses indicate that the source area (5.51 km²) is ～20% smaller than the topographically defined watershed (6.85 km²), and recharge entering ～37% of the surface watershed does not discharge to Phillips Branch. Groundwater residence time within the source volume ranges from a few days to almost 100 years, with 95% of the volume “flushing” within 50 years. Artificial discharge from groundwater pumping alters the shape of the source area and reduces baseflow due to the interception of stream flow paths, but has limited impacts on the residence time of groundwater discharged as baseflow. In contrast, artificial recharge from land-based wastewater disposal substantially reduces the source area, lowers the range in residence time due to the elimination of older flow paths to the stream, and leads to increased discharge to adjacent surface-water bodies. This research suggests that, in this and similar hydrogeologic settings, the “watershed” approach to water-resource management may be limited, particularly where anthropogenic stresses alter the transport of soluble contaminants through highly permeable unconfined aquifers.     

Multiscalarity of the Nutrient–Chlorophyll Relationship in Coastal Phytoplankton

The relationship between nitrogen (N) and phytoplankton chlorophyll a (Chl) establishes a basis for understanding eutrophication in coastal marine ecosystems. A substantial literature exists on cross-ecosystem analysis of this relationship, but there is little information on cross-scale patterns. A collection of observational records in Bedford Basin (Canada) was used to construct the N–Chl relationship at four time scales: intra-day, intra-annual, interannual, and interdecadal. Additionally, a dataset of contingent observations from 16 biogeochemical ocean provinces was used to construct the N–Chl relationship at large spatial scale. In Bedford Basin, N statistically predicts Chl at time scales that are short (intra-day, intra-annual) and long (interdecadal) but not intermediate (interannual). There is an apparent stoichiometric regularity in the dependence of Chl on N that crosses time scales. Further, an apparent similitude exists between the local pattern at long time scale and the global pattern at large space scale.     

Regional patterns of U.S. household carbon emissions

Market-based policies to address fossil fuel-related externalities including climate change typically operate by raising the price of those fuels. Increases in energy prices have important consequences for a typical U.S. household that spent almost $4,000 per year on electricity, fuel oil, natural gas, and gasoline in 2005. A key question for policymakers is how these consequences vary over different regions and subpopulations across the country—especially as adjustment and compensation programs are designed to protect more vulnerable regions. To answer this question, we use non-publicly available data from the U.S. Consumer Expenditure Survey over the period 1984–2000 to estimate long-run geographic variation in household use of electricity, fuel oil, natural gas, and gasoline, as well as the associated incidence of a $10 per ton tax on carbon dioxide (ignoring behavioral response). We find substantial variation: incidence from the tax range from $97 dollars per year per household in New York County, New York to $235 per year per household in Tensas Parish, Louisiana. This variation can be explained by differences in energy use, carbon intensity of electricity generation, and electricity regulation.