Simulation of the full two rigid body problem using polyhedral mutual potential and potential derivatives approach

Herein we investigate the coupled orbital and rotational dynamics of two rigid bodies modelled as polyhedra, under the influence of their mutual gravitational potential. The bodies may possess any arbitrary shape and mass distribution. A method of calculating the mutual potential’s derivatives with respect to relative position and attitude is derived. Relative equations of motion for the two body system are presented and an implementation of the equations of motion with the potential gradients approach is described. Results obtained with this dynamic simulation software package are presented for multiple cases to validate the approach and illustrate its utility. This simulation capability is useful both for addressing questions in dynamical astronomy and for enabling spacecraft missions to binary asteroid systems.     

Is there a black hole among the black-hole candidates?

None of the tentative black-hole identifications has remained without serious difficulties. As a class, the black-hole candidates differ in no obvious property from neutron-star binaries other than in their estimated high mass. We favour an interpretation of the black-hole candidates as binary neutron stars surrounded by a more or less massive accretion disc     

Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study

A Lagrangian particle method embedded within a 2-D finite element code, is used to study the transport and ocean–estuary exchange processes in the well-mixed Great Bay Estuarine System in New Hampshire, USA. The 2-D finite element model, driven by residual, semi-diurnal and diurnal tidal constituents, includes the effects of wetting and drying of estuarine mud flats through the use of a porous medium transport module. The particle method includes tidal advection, plus a random walk model in the horizontal that simulates sub-grid scale turbulent transport processes. Our approach involves instantaneous, massive [O(500,000)] particle releases that enable the quantification of ocean–estuary and inter-bay exchanges in a Markovian framework. The effects of the release time, spring–neap cycle, riverine discharge and diffusion strength on the intra-estuary and estuary–ocean exchange are also investigated.The results show a rather dynamic interaction between the ocean and the estuary with a fraction of the exiting particles being caught up in the Gulf of Maine Coastal Current and swept away. Three somewhat different estimates of estuarine residence time are calculated to provide complementary views of estuary flushing. Maps of residence time versus release location uncover a strong spatial dependency of residence time within the estuary that has very important ramifications for local water quality. Simulations with and without the turbulent random walk show that the combined effect of advective shear and turbulent diffusion is very effective at spreading particles throughout the estuary relatively quickly, even at low (1 m²/s) diffusivity. The results presented here show that a first-order Markov Chain approach has applicability and a high potential for improving our understanding of the mixing processes in estuaries.     

Libration effects for retrograde satellites in the restricted three-body problem

Numerical explorations of the restricted problem have shown that for stable large nonperiodic retrograde satellite orbits, the motion can be decomposed into a fast reference motion and a slow libration aroundB ₂ We study here this libration in the circular plane Hill's case, for which the reference motion is elliptic. We establish the equations of motion for the coordinates of the centre of this ellipse. We find two integrals of motion: the first is the semi-major axis of the ellipse; the second is essentially Jacobi's integral, translated into the new coordinates. We give a formula for the period of the libration and we find its limiting value for small libration amplitudes. A numerical verification gives very good agreement for all these results.     

The fusion crust of the Winchcombe meteorite: A preserved record of atmospheric entry processes

Fusion crusts form during the atmospheric entry heating of meteorites and preserve a record of the conditions that occurred during deceleration in the atmosphere. The fusion crust of the Winchcombe meteorite closely resembles that of other stony meteorites, and in particular CM2 chondrites, since it is dominated by olivine phenocrysts set in a glassy mesostasis with magnetite, and is highly vesicular. Dehydration cracks are unusually abundant in Winchcombe. Failure of this weak layer is an additional ablation mechanism to produce large numbers of particles during deceleration, consistent with the observation of pulses of plasma in videos of the Winchcombe fireball. Calving events might provide an observable phenomenon related to meteorites that are particularly susceptible to dehydration. Oscillatory zoning is observed within olivine phenocrysts in the fusion crust, in contrast to other meteorites, perhaps owing to temperature fluctuations resulting from calving events. Magnetite monolayers are found in the crust, and have also not been previously reported, and form discontinuous strata. These features grade into magnetite rims formed on the external surface of the crust and suggest the trapping of surface magnetite by collapse of melt. Magnetite monolayers may be a feature of meteorites that undergo significant degassing. Silicate warts with dendritic textures were observed and are suggested to be droplets ablated from another stone in the shower. They, therefore, represent the first evidence for intershower transfer of ablation materials and are consistent with the other evidence in the Winchcombe meteorite for unusually intense gas loss and ablation, despite its low entry velocity.     

A seasonal climate model of the atmospheres of the giant planets at the voyager encounter time: I. Saturn's stratosphere

A radiative seasonal model which incorporates a multilayer radiative transfer treatment at wave-lengths longward of 7 μm is presented and applied to Saturn's stratosphere. Opacities due to H₂-He, CH₄, C₂H₂, and C₂H₆ are included. Season-dependent insolation is shown to produce a strong hemispheric asymmetry decreasing with depth at the Voyager encounter times, and seasonal amplitudes of 30°K at the poles are predicted in the high stratosphere. The ring-modulated dependence of the insolation and the orbital eccentricity are shown to have a significant effect. Calculations agree closely with the Voyager 1 and 2 radio occultation ingress profiles recorded at 76°S and 36.5°S for CH₄/H₂ = 3.5 + 1.4/? 1.0 × 10^?3;the estimated errors include modeling systematic errors and uncertainties in the occultations profiles. The possible role of aerosols in the stratospheric heating is analyzed. The Voyager 2 egress profile recorded at 31°S cannot be reproduced by calculations. Some constraints on the C₂H₂ and C₂H₆ abundances are derived. The upper portion of the occultation profiles (p < 3mbar) can be matched for C₂H₂/H₂ = 1.0 + 1.3/?0.6 × 10^?7, C₂H₆/H₂ = 1.5 + 1.8/?0.9 × 10^?6 at 76°S and C₂H₂/H₂ = 4 + 6/?4 × 10^?8, C₂H₆/H₂ = 6 + 9/?6 × 10^?7 at 36.5°N. At the northern occultation latitude, the discrepancy with the concentrations derived from analysis of IRIS spectra by R. Courtin, D. Gautier, A. Marten, B. Bézard, and R. Hanel (1984, Astrophys. J.287) can be explained by a sharp variation of the mixing ratios of these gases with altitude in the upper stratosphere. Other interpretations are discussed.     

Estimation of phosphorus loads with sparse data for agricultural watersheds in the Czech Republic

Abstract

Agricultural watersheds in the Czech Republic are one of the primary sources of non-point-source phosphorus (P) loads in receiving waters. Since such non-point sources are generally located in headwater catchments, streamflow and P concentration data are sparse. We show how very short daily streamflow and P concentration records can be combined with nearby longer existing daily streamflow records to result in reliable estimates of daily and annual P concentrations and loads. Maintenance of variance streamflow record extension methods (MOVE) can be employed to extend short streamflow records. Constituent load regressions are used to predict daily P constituent loads from streamflow and other time varying characteristics. Annual P loads are then estimated for individual watersheds. Resulting annual P load estimates ranged from 0.21 to 95.4 kg year^-1 with a mean value of 11.77 kg year^-1. Similarly annual P yield estimates ranged from 0.01 to 0.3 kg ha^-1 year^-1 with an average yield of 0.07 kg ha^-1 year^-1. We document how short records of daily streamflow and P concentrations can be combined with a national network of daily streamflow records in the Czech Republic to arrive at meaningful and reliable estimates of annual P loads for small agricultural watersheds.

Citation Beránková, T., Vogel, R. M., Fiala, D. & Rosendorf, P. (2010) Estimation of phosphorus loads with sparse data for agricultural watersheds in the Czech Republic. Hydrol. Sci. J. 55(8), 1417–1426.     

Biodiversity in Hudson River shore zones: influence of shoreline type and physical structure

The shore zones of the Hudson River, like those of many developed waterways, are highly varied, containing a mix of seminatural and highly engineered shores. Our goal was to document the biodiversity supported by different kinds of shore zones in the Hudson. We chose six common types of shore zones, three of them ??natural?? (sand, unconsolidated rock, and bedrock), and three of them engineered (riprap, cribbing, and bulkheads). We measured selected physical characteristics (shore zone width, exposure, substrate roughness and grain size, shoreline complexity) of three examples of each of these shore types, and also sampled communities of terrestrial plants, fishes, and aquatic and terrestrial invertebrates. Community composition of most taxa differed across shore types, and frequently differed between wide, sheltered shores and narrow, exposed shores. Alien plant species were especially well represented along engineered shores. Nevertheless, a great deal of variation in biological communities was not explained by our six-class categorization of shore zones or the physical variables that we measured. No single shore type supported the highest values of all kinds of biodiversity, but engineered shore zones (especially cribbing and bulkheads) tended to have less desirable biodiversity characteristics than ??natural?? shore zones.     

Cytochrome P4501A biomarker indication of the timeline of chronic exposure of Barrow's goldeneyes to residual Exxon Valdez oil

We examined hepatic EROD activity, as an indicator of CYP1A induction, in Barrow’s goldeneyes captured in areas oiled during the 1989 Exxon Valdez spill and those from nearby unoiled areas. We found that average EROD activity differed between areas during 2005, although the magnitude of the difference was reduced relative to a previous study from 1996/1997, and we found that areas did not differ by 2009. Similarly, we found that the proportion of individuals captured from oiled areas with elevated EROD activity (?2 times unoiled average) declined from 41% in winter 1996/1997 to 10% in 2005 and 15% in 2009. This work adds to a body of literature describing the timelines over which vertebrates were exposed to residual Exxon Valdez oil and indicates that, for Barrow’s goldeneyes in Prince William Sound, exposure persisted for many years with evidence of substantially reduced exposure by 2 decades after the spill.     

Evaluating carbon storage on subalpine lake deltas

Mountainous regions are important contributors to the terrestrial organic carbon (OC) sink that affect global climate through the regulation of carbon‐based greenhouse gases. However, mountain OC dynamics are poorly quantified. We quantified OC storage in subalpine lake deltas in the Washington Central Cascades and Colorado Front Range with the objectives of determining the magnitude of transient carbon storage and understanding the differences in storage between the two ranges. We used field, laboratory, and GIS techniques to determine the magnitude of and controls on the subalpine lake delta OC pool in 26 subalpine lake deltas. Soil moisture, soil texture, mean basin slope, and delta valley confinement are significantly correlated with soil carbon on deltas. Average soil OC concentration on subalpine lake deltas ranges from 3 to 41%, and stocks range from 140 to 1256 Mg C/ha. Surprisingly, the carbon content of subalpine lake deltas is not significantly different between the two regions, despite stark contrasts in their climate, vegetation, and total ecosystem carbon stocks. We present a conceptual model that invokes geomorphic and biogeochemical processes to suggest that carbon is more likely to reach subalpine lake deltas from the upstream basin in the Colorado Front Range compared with the Washington Central Cascades, thus accounting for the similarity in OC storage between the two regions despite differences in total ecosystem carbon stocks and climate. This points to a complex interaction among carbon production, transport, and stability in each region, and supports the idea that geomorphic and biogeochemical processes determine the magnitude of transient OC storage more strongly than primary productivity or climate. Copyright © 2017 John Wiley & Sons, Ltd.