Planetary orbital equations in externally-perturbed systems: position and velocity-dependent forces

The increasing number and variety of extrasolar planets illustrates the importance of characterizing planetary perturbations. Planetary orbits are typically described by physically intuitive orbital elements. Here, we explicitly express the equations of motion of the unaveraged perturbed two-body problem in terms of planetary orbital elements by using a generalized form of Gauss’ equations. We consider a varied set of position and velocity-dependent perturbations, and also derive relevant specific cases of the equations: when they are averaged over fast variables (the “adiabatic” approximation), and in the prograde and retrograde planar cases. In each instance, we delineate the properties of the equations. As brief demonstrations of potential applications, we consider the effect of Galactic tides. We measure the effect on the widest-known exoplanet orbit, Sedna-like objects, and distant scattered disk objects, particularly with regard to where the adiabatic approximation breaks down. The Mathematica code which can help derive the equations of motion for a user-defined perturbation is freely available upon request.     

Trophic modeling of the Northern Humboldt Current Ecosystem, Part II: Elucidating ecosystem dynamics from 1995 to 2004 with a focus on the impact of ENSO

The Northern Humboldt Current Ecosystem is one of the most productive in the world in terms of fish production. Its location near to the equator permits strong upwelling under relatively low winds, thus creating optimal conditions for the development of plankton communities. These communities ultimately support abundant populations of grazing fish such as the Peruvian anchoveta, Engraulis ringens. The ecosystem is also subject to strong inter-annual environmental variability associated with the El Niño Southern Oscillation (ENSO), which has major effects on nutrient structure, primary production, and higher trophic levels. Here our objective is to model the contributions of several external drivers (i.e. reconstructed phytoplankton changes, fish immigration, and fishing rate) and internal control mechanisms (i.e. predator-prey) to ecosystem dynamics over an ENSO cycle. Steady-state models and time-series data from the Instituto del Mar del Perú (IMARPE) from 1995 to 2004 provide the base data for simulations conducted with the program Ecopath with Ecosim. In simulations all three external drivers contribute to ecosystem dynamics. Changes in phytoplankton quantity and composition (i.e. contribution of diatoms and dino- and silicoflagellates), as affected by upwelling intensity, were important in dynamics of the El Niño of 1997–98 and the subsequent 3 years. The expansion and immigration of mesopelagic fish populations during El Niño was important for dynamics in following years. Fishing rate changes were the most important of the three external drivers tested, helping to explain observed dynamics throughout the modeled period, and particularly during the post-El Niño period. Internal control settings show a mix of predator–prey control settings; however a “wasp-waist” control of the ecosystem by small pelagic fish is not supported.     

Direct terrestrial-marine correlation demonstrates surprisingly late onset of the last interglacial in central Europe

An interdisciplinary study of a small sedimentary basin at Neumark Nord 2 (NN2), Germany, has yielded a high-resolution record of the palaeomagnetic Blake Event, which we are able to place at the early part of the last interglacial pollen sequence documented from the same section. We use this data to calculate the duration of this stratigraphically important event at 3400 ± 350 yr. More importantly, the Neumark Nord 2 data enables precise terrestrial-marine correlation for the Eemian stage in central Europe. This shows a remarkably large time lag of ca. 5000 yr between the MIS 5e ‘peak’ in the marine record and the start of the last interglacial in this region.     

From subduction to transform motion: a seabeam survey of the Hellenic trench system

Preliminary results of a multi-narrow beam survey of the Hellenic trench system, in the Eastern Mediterranean, are presented. The southwestern Ionian branch is divided in small basins, partly filled with Pleistocene sediments. The morphology suggests that the basins are deformed by a compressional stress acting roughly perpendicularly to the trench along N50°E. This direction is the direction of the regional slip vector of the shallow thrust-type earthquakes. The structure of the southeastern Pliny-Strabo branch is quite different. Narrow en-e´chelon slots, oriented N40°E, have been mapped within the main troughs oriented N60°E. The regional earthquake slip vector is also oriented along N40°E. We conclude that the Hellenic trench system is an active subduction system, dominated by thrust along the Ionian branch and by transform motion along the Pliny-Strabo branch.     

A surface complexation model for sulfate and selenate on iron oxides consistent with spectroscopic and theoretical molecular evidence

Sulfate and selenate adsorption on iron oxides are important reactions in natural systems under a very wide range of pH values, ionic strengths, and electrolyte compositions. Under such conditions, spectroscopic and theoretical calculations have demonstrated the potential importance of a variety of surface species. Understanding the variations in the surface speciation of these oxyanions is fundamental to prediction of their partitioning between minerals and aqueous solutions. In the present study, published experimental spectroscopic and theoretical molecular evidence of the identity of sulfate/selenate surface species are integrated with a surface complexation model consistent with a wide variety of experimental adsorption, surface titration, and proton coadsorption data to define the surface speciation of sulfate and selenate on iron oxides under a wide range of conditions. The analysis was carried out with the extended triple layer model (ETLM) taking into account the electrostatics of water dipole desorption during ligand exchange reactions. On seven out of eight goethites studied, sulfate and selenate surface reactions can be represented by the formation of a monodentate-mononuclear inner-sphere and a bidentate-binuclear outer-sphere (or H-bonded) species according to

     

On a modal damping identification model for non-classically damped linear building structures subjected to earthquakes

A simple modal damping identification model developed by the present authors for classically damped linear building frames is extended here to the non-classically damped case. The modal damping values are obtained with the aid of the frequency domain modulus of the roof-to-basement transfer function and the resonant frequencies of the structure (peaks of the transfer function) as well as the modal participation factors and mode shapes of the undamped structure. The assumption is made that the modulus of the transfer function of the non-classically damped structure matches the one of the classically damped structure in a discrete manner, i.e., at the resonant frequencies of that function modulus. This proposed approximate identification method is applied to a number of plane building frames with and without pronounced non-classical damping under different with respect to their frequency content earthquakes and its limitations and range of applicability are assessed with respect to the accuracy of both the identified damping ratios and that of the seismic structural response obtained by classical mode superposition and use of those identified modal damping ratios.     

Microstructure and hydraulic conductivity of a compacted lime-treated soil

Under a given compaction energy and procedure, it is known that maximum dry density of a soil is lowered due to lime addition. This modification of maximum dry density could alter the hydraulic conductivity of the soil. The main object of this study was to assess the impact of lime-stabilization on a silt soil microstructure and then on saturated hydraulic conductivity. An investigation at the microscopic level with mercury intrusion porosimetry showed that lime treatment induced the formation of a new small class, with a diameter lower than 3 × 10³ Å in the compacted soil. This class is responsible for the difference in dry density between the treated and the untreated sample after compaction. It is shown that this small pores class was not altered by the compaction water content, the compaction procedure or the dry density. As in untreated soils, only the larger pores were modified by the compaction water content and the compaction procedure in the lime treated samples. The hydraulic conductivity appeared to be only related to the largest pores volume of the tested silt, regardless of lime treatment. Therefore, this study demonstrated that even if addition of lime resulted in a dramatic change of the maximum dry density of the tested silty soil, its effect on hydraulic conductivity is limited.