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).     

The impact of perturbations to ocean-model parameters on climate and climate change in a coupled model

Assessments of the impacts of uncertainties in parameters on mean climate and climate change in complex climate models have, to date, largely focussed on perturbations to parameters in the atmosphere component of the model. Here we expand on a previously published study which found the global impacts of perturbed ocean parameters on the rate of transient climate change to be small compared to perturbed atmosphere parameters. By separating the climate-change-induced ocean vertical heat transport in each perturbed member into components associated with the resolved flow and each parameterisation scheme, we show that variations in global mean heat uptake in different perturbed versions are an order of magnitude smaller than the average heat uptake. The lack of impact of the perturbations is attributed to (1) the relatively small impact of the perturbation on the direct vertical heat transport associated with the perturbed process and (2) a compensation between those direct changes and indirect changes in heat transport from other processes. Interactions between processes and changes appear to combine in complex ways to limit ensemble spread and uncertainty in the rate of warming. We also investigate regional impacts of the perturbations that may be important for climate change predictions. We find variations across the ensemble that are significant when measured against natural variability. In terms of the experimental set-up used here (models without flux adjustments) we conclude that perturbed physics ensembles with ocean parameter perturbations are an important component of any probabilistic estimate of future climate change, despite the low spread in global mean quantities. Hence, careful consideration should be given to assessing uncertainty in ocean processes in future probabilistic assessments of regional climate change.     

Spatial delineation of groundwater–surface water interactions through intensive in‐stream profiling

The dominance of ‘old’ pre‐event water in headwater storm runoff has been recorded in numerous upland catchment studies; however, the mechanisms by which this pre‐event water enters the stream channel are poorly understood. Understanding these processes is fundamental to determining the controls on surface water quality and associated impacts on stream ecology. Previous studies in the upland forested catchment of the Afon Hafren (River Severn) at Plynlimon, mid‐Wales, identified an active bedrock groundwater system that was discharging into the stream channel during storm response. Detailed analysis showed that these discharges were small and could not account for the majority of pre‐event storm water response identified at this site; pre‐event storm runoff had to be sourced predominantly from further upstream. An intensive stream survey was used to determine the spatial nature of groundwater–surface water (GW–SW) interactions in the Hafren Catchment. Detailed physico‐chemical in‐stream profiling identified a marked change in water quality indicating a significant discrete point of bedrock groundwater discharge upstream of the Hafren Transect study site. The in‐stream profiling showed the importance of high spatial resolution sampling as a key to understanding processes of GW–SW interaction and how quick and cost‐effective measurements of specific electrical conductance of stream waters could be used to highlight in‐stream heterogeneity. This approach is recommended for use in headwater catchments for initial characterisation of the stream channel in order to better locate instrumentation and to determine more effective targeted sampling protocols in upland catchment research. Copyright © 2012 John Wiley & Sons, Ltd.     

Merging chironomid training sets: implications for palaeoclimate reconstructions

Icelandic and Norwegian chironomid calibration or training sets were merged to investigate whether a larger combined training set would be useful to apply to subfossil chironomid data from Iceland for periods such as the early Holocene, the Holocene Thermal Maximum and the Little Ice Age, when temperatures can be expected to be outside the current temperature range of the Icelandic training set. Following taxonomic harmonisation, the Icelandic and Norwegian data sets were compared before being merged to form a combined Norwegian-Icelandic training set. Analyses showed that it was biologically and statistically valid to merge the two data sets. The resulting combined inference model for mean July air temperature had improved performance statistics (r²_jack = 0.87; RMSEP_jack = 1.13) when compared to the best performing Icelandic model (r²_jack = 0.61; RMSEP_jack = 0.83), due to the longer environmental gradient covered (Icelandic 6–11 °C; combined 3.5–16 °C), and to the increased number of samples (Icelandic = 53 lakes; combined = 207 lakes) and taxa (Icelandic = 47 taxa; combined = 133 taxa) present within the combined training set. The inference models were applied to an early Holocene chironomid sequence from Vatnamýri, north Iceland, and a 450-year recent record from Myfluguvatn, north-west Iceland, to compare the reconstructions produced. The various inference models produced similar trends and patterns of temperature reconstruction, but the inference model based on the combined training set produced a larger range of reconstructed temperatures than the Icelandic model. It was found that different inference models produced more variation in the reconstruction than when different training sets were used. A comparison of the Myfluguvatn reconstructions with meteorological observations showed that the combined Norwegian–Icelandic inference model produced more reliable results than the Icelandic or Norwegian inference models alone.     

MegaZ-LRG: a photometric redshift catalogue of one million SDSS luminous red galaxies

We describe the construction of MegaZ-LRG, a photometric redshift catalogue of over one million luminous red galaxies (LRGs) in the redshift range  0.4 < z < 0.7  with limiting magnitude   i < 20  . The catalogue is selected from the imaging data of the Sloan Digital Sky Survey (SDSS) Data Release 4. The 2dF-SDSS LRG and Quasar (2SLAQ) spectroscopic redshift catalogue of 13 000 intermediate-redshift LRGs provides a photometric redshift training set, allowing use of ann z, a neural network-based photometric-redshift estimator. The rms photometric redshift accuracy obtained for an evaluation set selected from the 2SLAQ sample is  σ _z = 0.049  averaged over all galaxies, and  σ _z = 0.040  for a brighter subsample  ( i < 19.0)  . The catalogue is expected to contain ∼5 per cent stellar contamination. The ann z code is used to compute a refined star/galaxy probability based on a range of photometric parameters; this allows the contamination fraction to be reduced to 2 per cent with negligible loss of genuine galaxies. The MegaZ-LRG catalogue is publicly available on the World Wide Web from http://www.2slaq.info .     

Evaluation of image-based multibeam sonar backscatter classification for benthic habitat discrimination and mapping at Stanton Banks,UK

In recent years, efforts have increased to develop quantitative, computer-directed methods for segmentation of multibeam (MBES) backscatter data. This study utilises MBES backscatter data acquired at Stanton Banks (UK) and subsequently processed through the QTC-Multiview software environment in a bid to evaluate the program's ability to perform unsupervised classification. Statistical comparison with ground-truth data (grab, stills and video) enabled cross validation of acoustic segmentation and biological assemblages observed at the site. 132 unspecified variables were extracted from user-specified rectangular patches of the backscatter image, reduced to three vectors by PCA, then clustered and classified by the software. Multivariate analyses of ground-truth data were conducted on 75 stills images and 51 grab samples. Video footage coincident with the stills was divided into 30 s segments and coded by dominant substrate and species. Cross tabulation determined the interrelationship between software classifications, multivariate analysis of the biological assemblages and coded video segments. Multiview optimally identified 19 classes using the automated clustering engine. These were revised to 6 habitats a posteriori, using combined analysis of ground-truth data and Multiview data products. These habitats broadly correspond to major physiographic provinces within the region. Multivariate statistical analysis reveals low levels of assemblage similarity (<35%) for samples occurring within Multiview classes, irrespective of the mode of acquisition. Coded video data is more spatially appropriate than the other methods of ground-truthing investigated, although it is less well suited to the extraction of truly quantitative data. Multivariate analysis indicates assemblages within physiographically distinct Multiview classes have a low degree of biological similarity, supporting the notion that abiotic proxies may be contraindicative of benthic assemblage variations. QTC-Multiview performs well as a mechanism for computer-assisted segmentation of MBES backscatter imagery into acoustic provinces; however a degree of caution is required prior to ascribing ecological significance to these classifications.     

Health profiles of Hamilton: Spatial characterisation of neighbourhoods for health investigations

This paper is part of a larger research program which employs a mixed-methods approach to study the determinants of health at the local level using specific neighborhoods in Hamilton, Ontario, Canada. In this paper, multivariate, spatial statistical techniques and geographic information systems are used to address questions about the characterization of neighbourhoods, based on socioeconomic determinants of health and risk factors such as smoking. While neighbourhood characterization has been a component of public health surveillance for some time, geostatistical techniques can now be used to derive more accurate representation of neighbourhoods for use in subsequent analysis. We utilize principal components analysis to reduce the data and extract the components that represent the underlying local processes. Principal components are also overlayed on comparative mortality figures to visualize where the socio-demographic determinants of health correspond spatially with mortality patterns. Predicted values from the components are then analysed for spatial clustering using local indicators of spatial association. The findings reveal a pattern of distinct neighbourhoods that will be used in subsequent quantitative and qualitative stages in the larger research programme. The results can also be used to inform public health policy and to target public health interventions.     

The braided Milk River, northern Montana, fails the Leopold–Wolman discharge-gradient test

The Milk River, the northernmost tributary to the Missouri–Mississippi River system, exhibits an anomalous sand-bed braiding reach in an otherwise meandering system. Shortly after leaving Alberta and entering Montana the river suddenly changes to braiding and maintains this pattern for 47 km before entering Fresno Reservoir. Measured stream gradient and bankfull discharge in the braiding reach severely fail the Leopold and Wolman [U.S. Geol. Surv. Prof. Pap. 282B (1957) 39] slope–discharge test for differentiating channel patterns. While channel slope has long been regarded as one of the primary variables associated with braiding, our data from the sand-bed Milk River do not support this relationship. Instead, the data show that the braiding reach has a lower channel slope (0.00047) than the meandering reach (0.00055). Coupled with a constant discharge the unit length stream power is comparable between the two reaches. At the morphologic transition between meandering and braiding, a dramatic reduction in channel bank strength occurs where the sampled silt–clay content declines from 65% in the meandering reach to 18% in the braiding. This enables channel widening which is reflected in a 60% reduction in unit area stream power in the braiding reach. Thus, sediment transport capacity declines and channel bars are deposited. During waning flows, these bars are dissected, producing a braiding morphology. We suggest that for sand-bed braiding rivers the silt–clay percentage in the channel banks may be more important than slope. A review of the original Leopold and Wolman [U.S. Geol. Surv. Prof. Pap. 282B (1957) 39] dataset, and many subsequent analyses, reveals that most braided rivers studied were gravel-bed. As a result, causal variables associated with braiding in sand-bed environments may need a thorough evaluation.     

Oxygen isotope geochemistry of the Mesozoic anorogenic complexes of Damaraland,northwest Namibia: evidence for crustal contamination and its effect on silica saturation

 Fifty-six new oxygen isotope analyses of minerals separated from nine Mesozoic anorogenic complexes of Damaraland in northwest Namibia have been used to estimate the δ¹⁸O values of the original magmas (δ_magma). These complexes range in composition from nepheline syenite to alkaline and peraluminous granites, often with a variety of rock types present at a single centre. The silica-undersaturated rock types show a relatively small spread of values of δ_magma from 6.0 to 7.0‰ (mean 6.6‰), which is consistent with their derivation from the mantle with little or no subsequent crustal contamination. The silica-oversaturated rocks show a wide range of δ_magma values from 4.9 to 12.0‰, with a mean value of 8.3‰. The high values of δ_magma in the silica-oversaturated rocks can only be explained with considerable involvement of the continental crust in their petrogenesis, and those rocks with δ_magma>10‰ are interpreted as essentially crustal melts. It is generally accepted that the Damaraland complexes were generated as a result of rifting across the Tristan plume, with the plume providing both mantle-derived magma and with it the heat required for crustal melting. In addition to their mantle-like oxygen isotope ratios, the undersaturated rocks of the Damaraland complexes possess the trace element characteristics (e.g. low Zr/Nb ratios) of ocean island basalt, suggesting that their parental magmas were produced from the plume itself. In contrast, the oversaturated complexes generally have higher Zr/Nb ratios that are consistent with a larger crustal input. The highest values of δ_magma in the Damaraland complexes are found in granitic rocks that intrude the central zone of the Pan-African Damara Orogen where presumably there is a substantial component of sedimentary origin in the lower to middle crust. Received: 3 April 1995/Accepted: 24 July 1995