N-body simulations of cohesion in dense planetary rings: A study of cohesion parameters

We present results from a large suite of simulations of Saturn’s dense A and B rings using a new model of particle sticking in local simulations (Perrine, R.P., Richardson, D.C., Scheeres, D.J. [2011]. Icarus 212, 719–735). In this model, colliding particles can be incorporated into or help fragment rigid aggregations on the basis of certain user-specified parameters that can represent van der Waals forces or interlocking surface frost layers.Our investigation is motivated by laboratory results that show that interpenetration of surface layers can allow impacting frost-covered ice spheres to stick together. In these experiments, cohesion only occurs below specific impact speeds, which happen to be characteristic of impact speeds in Saturn’s rings. Our goal is to determine if weak bonding is consistent with ring observations, to constrain cohesion parameters in light of existing ring observations, to make predictions about particle populations throughout the rings, and to discover other diagnostics that may constrain bonding parameters.We considered the effects of five parameters on the equilibrium characteristics of our ring simulations: speed-based merge and fragmentation limits, bond strength, ring surface density, and patch orbital distance (i.e., the A or B ring), some with both monodisperse and polydisperse comparison cases. In total, we present data from 95 simulations.We find that weak cohesion is consistent with observations of the A and B rings (e.g., French, R.G., Nicholson, P.D. [2000]. Icarus 145, 502–523), and we present a range of simulation parameters that reproduce the observed size distribution and maximum particle size. It turns out that the parameters that match observations differ between the A and B rings, and we discuss the potential implications of this result. We also comment on other observable consequences of cohesion for the rings, such as optical depth and scale height effects, and discuss whether very large objects (e.g., “propeller” source objects) are grown bottom-up from cohesion of smaller ring particles.     

Numerical simulations of granular dynamics II: Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties.We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration.We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.     

Mineral chemical provenance of Neolithic pitchstone artefacts from Ballygalley,County Antrim,Northern Ireland

A large number of pitchstone fragments and artefacts have recently been discovered at a Neolithic settlement site in Ballygalley, Co. Antrim, Northern Ireland. They consist predominantly of flakes and un‐reworked lumps and cores, with only one complete tool being found. Since no sources of workable pitchstone exist in Ireland, the source must have been the abundant pitchstone volcanic rocks found on the Hebridean Islands of northwest Scotland. The composition of the glass from a number of artefacts is highly siliceous, indicating that they were derived from pitchstones on the Island of Arran; pitchstones from all other Scottish locations are less silica‐rich. In addition, analysis of pyroxene and amphibole microcrystallites within the pitchstone suggests that the Corriegills area of Arran is the most likely source of the Ballygalley artefacts, although the precise outcrop has proved elusive. These finds, and others across Ireland, show that raw materials were being transported and probably traded over considerable distances despite there being suitable alternative sources of material for making tools (flint, etc.) available in the local area. This suggests that the pitchstone had a very specialist use. © 2002 Wiley Periodicals, Inc.     

Macroseepage of methane and light alkanes at the La Brea tar pits in Los Angeles

Geologic seepage of methane and light (C₂-C₅) alkanes was measured at the La Brea Tar Pits in Los Angeles. Samples were collected using flux chambers with stainless steel canisters and analyzed using gas chromatography. Average seepage rates from individual seepage sites were 970 ± 330 mg/h of methane, 14.0 ± 5.5 mg/h of ethane, 9.1 ± 3.7 mg/h of propane, 3.7 ± 1.6 mg/h of i-butane, 0.33 ± 0.16 mg/h of n-butane, 260 ± 120 μg/h of i-pentane, and 5.3 ± 1.9 μg/h of n-pentane, while maximum seepage rates exceeded 17 g/h of methane, 270 mg/h of ethane, 190 mg/h of propane, 95 mg/h of i-butane, 10 mg/h of n-butane, 7 mg/h of i-pentane, and 0.1 mg/h of n-pentane. These absolute fluxes have an additional unknown amount of error associated with them due to sampling methodology, and should be taken as the lower limit of emissions. Samples collected revealed generally dry gas, with high methane emissions relative to the light alkanes. Overall emissions from the tar pits were found to come not only from the active geologic seepage, but also from the outgassing of the standing asphalt at the site. Using the gas ratios, which are negligibly affected by errors introduced by sampling methodology, observed in this study, daily emissions of C₂ – C₅ alkanes from the La Brea area were estimated to be 4.7 ± 1.6 Mg, which represents 2–3 % of total emissions in the entire Los Angeles region.     

Co-management and the co-production of knowledge: Learning to adapt in Canada's Arctic

Co-management institutional arrangements have an important role in creating conditions for social learning and adaptation in a rapidly changing Arctic environment, although how that works in practice has not been clearly articulated. This paper draws on three co-management cases from the Canadian Arctic to examine the role of knowledge co-production as an institutional trigger or mechanism to enable learning and adapting. Experience with knowledge co-production across the three cases is variable but outcomes illustrate how co-management actors are learning to learn through uncertainty and environmental change, or learning to be adaptive. Policy implications of this analysis are highlighted and include the importance of a long-term commitment to institution building, an enabling policy environment to sustain difficult social processes associated with knowledge co-production, and the value of diverse modes of communication, deliberation and social interaction.     

Geochemistry of ice-covered, meromictic Lake A in the Canadian High Arctic

The geochemical processes occurring within meromictic Lake A (maxdepth 120 m) on northern Ellesmere Island, Canada, were investigated to determinethe history of the lake and to provide a baseline for future studies. The lake contained seawaterdiluted by freshwater input that had been mixed prior to the lake's isolation from tidalaction. Input of freshwater after isolation of the lake created vertical stratification resultingin the creation of distinct oxic, suboxic and anoxic zones. Dissolved oxygen was present to 13 m,and sulphide beneath 32 m. Manganese and iron cycling dominated the redox chemistrybetween these depths. Total manganese concentrations reached 176 M, higher thanin most other natural stratified lake or marine environments.     

Analysis of Interpolation Methods for Kinematic DGPS Control in Aerial Photogrammetry

Kinematic differential Global Positioning System (DGPS) positioning is routinely used in industry for directly observing an aircraft's position at each instant of photographic exposure during a photogammetric survey. A critical aspect of the subsequent data processing is estimation of the aircraft position at the exact time of exposure. GPS measurements are acquired at a uniform sampling rate, typically 1 Hz. The exposure times, however, do not generally coincide with these times. As a result, the exposure station positions must be interpolated from the adjacent GPS positions. This is typically done using a low-order polynomial, expressed as a function of time, for each coordinate dimension. However, trajectory perturbations induced by atmospheric turbulence can render such interpolation methods ineffective. This article will convey the results of an investigation into the use of several different interpolation models with airborne GPS data during straight, level flight. The fundamental task of time series reconstruction will first be addressed, in which several possible interpolation models are described. Two 10-Hz, airborne GPS data sets were collected to test the accuracy of each model. The error properties resulting from the application of each model to these data will be presented and analyzed in terms of time-domain statistics and frequency-domain characteristids. It will be demonstrated that interpolation error can be significantly reduced, especially in the height dimension, through judicious choice of an interpolator. ? 2000 John Wiley & Sons, Inc.     

Rates of cave and landform development in the Yorkshire Dales from speleothem age data

Uranium-series ages have been obtained for 87 speleothems collected from nine major cave systems in the Craven district of northern England. Large systems such as Ease Gill Caverns, the West Kingsdale caves, and Gaping Gill-Ingleborough Cave, which contain relict high-level tunnels, are found to be older than 350,000 years (the limit of the ²³⁰Th/²³⁴U dating method). There is little evidence to indicate a significant enlargement of these passages since this time. Estimates of the age of Victoria Cave from ²³⁴U/²³⁸U isotopic ratios suggest that the cave has been fully relict for more than 500,000 years. Ages of in situ speleothems immediately adjacent to local base level cave streams show that mean maximum downcutting rates in limestone channels are about 2 to 5 cm 1000 yr. These rates are significantly lower than those obtained from direct measurements on limestone bedrock in stream channels but are comparable to areal denudation rates based on solute budgets. Using the present elevation of caves with respect to adjacent valley floor levels, mean maximum valley entrenchment rates are found to range between 5 cm/ky and < 20cm/ky, which corresponds to 6 m to < 24 m of lowering per glacial/interglacial cycle. These rates suggest that upper beds of the limestone were incised to form the Yorkshire Dales between 1 and 2 million years ago. The results indicate that the erosional effects of individual glaciations are not as severe as previously proposed in the literature.     

Stable isotope studies of fluid inclusions in speleothems and their paleoclimatic significance

Fluid inclusions found trapped in speleothems (cave deposited travertine) are interpreted as samples of seepage water from which enclosing calcium carbonate was deposited. The inclusions are assumed to have preserved their D/H ratios since the time of deposition. Initial ¹⁸O/¹⁶O ratios can be inferred from δD because rain- and snow-derived seepage waters fall on the meteoric water line (δD = 8δ¹⁸O + 10). Estimates of temperature of deposition of the carbonate can be calculated from inclusion D/H ratios and δ¹⁸O of enclosing calcite in Pleistocene speleothems. For most speleothems investigated (0–200,000 yr old) δ¹⁸O of calcite appears to have decreased with increasing temperature of deposition indicating that the dominant cause of climate-dependent change in δ¹⁸O of calcite was the change in K_cw, the isotope fractionation equilibrium constant, with temperature; δ¹⁸O of meteoric precipitation generally increased with increasing temperature, but not sufficiently to compensate for the decrease in K_cw.     

The Zn abundance and isotopic composition of diatom frustules,a proxy for Zn availability in ocean surface seawater

We have developed cleaning methods for extracting diatomopal from bulk marine sediment samples, for measurement of both zinc (Zn) abundance and isotope composition. This cleaning technique was then applied to a set of Holocene core-top samples from the Southern Ocean. The measured δ⁶⁶Zn (reported relative to the JMC_Lyon standard) and Zn/Si ratios from the Southern Ocean diatom_opal samples range from 0.7 to 1.5‰, and from 14 to 0.9 μmol/mol, respectively. The Zn abundance and isotope composition data show a clear correlation with opal burial rates and other oceanographic parameters. In common with previous work, we interpret the systematic changes in the Zn/Si ratio to be linked to the variability in the concentrations of bioavailable Zn in the ambient surface seawater where the diatom opal is formed. This variability is likely to be primarily controlled by the degree to which Zn is taken up into phytoplankton biomass. The observed systematic pattern in the δ⁶⁶Zn compositions of the diatom_opal core-top samples is, similarly, likely to reflect changes in the δ⁶⁶Zn composition of the ambient Zn in the surface waters above the core-top sites, which is progressively driven towards isotopically heavier values by preferential incorporation of the lighter isotopes into phytoplankton organic material. Thus, the systematic relationship between Zn isotopes and abundance observed in the core-top diatom_opal samples suggests a potential tool for investigating the biogeochemical cycling of Zn in the past surface ocean for down-core diatom_opal material. In this respect, it may be possible to test hypotheses that attribute variations in atmospheric CO₂ on glacial–interglacial timescales to the degree to which trace metals limited primary productivity in HNLC zones.