N-body integrators with individual time steps from Hierarchical splitting

We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping, these integrators are momentum conserving and show excellent energy conservation in conjunction with a symmetrized time step criterion. We use an explicit but approximate formula for the time symmetrization that is compatible with the use of individual time steps. No iterative scheme is necessary. We implement these ideas in the HUAYNO¹ code and present tests of the integrators and show that the presented integration schemes shows good energy conservation, with little or no systematic drift, while conserving momentum and angular momentum to machine precision for long term integrations.     

First ever in situ observations of Venus’ polar upper atmosphere density using the tracking data of the Venus Express Atmospheric Drag Experiment (VExADE)

On its highly elliptical 24 h orbit around Venus, the Venus Express (VEX) spacecraft briefly reaches a periapsis altitude of nominally 250 km. Recently, however, dedicated and intense radio tracking campaigns have taken place in August 2008, October 2009, February and April 2010, for which the periapsis altitude was lowered to the 186–176 km altitude range in order to be able to probe the upper atmosphere of Venus above the North Pole for the first time ever in situ. As the spacecraft experiences atmospheric drag, its trajectory is measurably perturbed during the periapsis pass, allowing us to infer total atmospheric mass density at the periapsis altitude. A Precise Orbit Determination (POD) of the VEX motion is performed through an iterative least-squares fitting process to the Doppler tracking data, acquired by the VEX radioscience experiment (VeRa). The drag acceleration is modelled using an initial atmospheric density model (VTS3 model, Hedin, A.E., Niemann, H.B., Kasprzak, W.T., Seiff, A. [1983]. J. Geophys. Res. 88, 73–83). A scale factor of the drag acceleration is estimated for each periapsis pass, which scales Hedin’s density model in order to best fit the radio tracking data. Reliable density scale factors have been obtained for 10 passes mainly from the second (October 2009) and third (April 2010) VExADE campaigns, which indicate a lower density by a factor of about 1.8 than Hedin’s model predicts. These first ever in situ polar density measurements at solar minimum have allowed us to construct a diffusive equilibrium density model for Venus’ thermosphere, constrained in the lower thermosphere primarily by SPICAV-SOIR measurements and above 175 km by the VExADE drag measurements (Müller-Wodarg et al., in preparation). The preliminary results of the VExADE campaigns show that it is possible to obtain with the POD technique reliable estimates of Venus’ upper atmosphere densities at an altitude of around 175 km. Future VExADE campaigns will benefit from the planned further lowering of VEX pericenter altitude to below 170 km.     

Analysis of energetic electron drop-outs in the upper atmosphere of Titan during flybys in the dayside magnetosphere of Saturn

We discuss the high energy electron absorption signatures at Titan during the Cassini dayside magnetospheric encounters. We use the electron measurements of the Low Energy Measurement System of the Magnetospheric Imaging Instrument. We also examine the mass loading boundary based on the ion data of the Ion Mass Spectrometer sensor of the Cassini Plasma Spectrometer. The dynamic motion of the Kronian magnetopause and the periodic charged particle flux and magnetic field variations – associated with the magnetodisk of Saturn – of the subcorotating magnetospheric plasma creates a unique and complex environment at Titan. Most of the analysed flybys (like T25–T33 and T35–T51) cluster at similar Saturn Local Time positions. However the instantaneous direction of the incoming magnetospheric particles may change significantly from flyby to flyby due to the very different magnetospheric field conditions which are found upstream of Titan within the sets of encounters.The energetic magnetospheric electrons gyrate along the magnetic field lines of Saturn, and at the same time bounce between the mirror points of the magnetosphere. This motion is combined with the drift of the magnetic field lines. When these flux tubes interact with the upper atmosphere of Titan, their content is depleted over approximately an electron bounce period. These depletion signatures are observed as sudden drop-outs of the electron fluxes. We examined the altitude distribution of these drop-outs and concluded that these mostly detected in the exo-ionosphere of Titan and sometimes within the ionosphere.However there is a relatively significant scatter in the orbit to orbit data, which can be attributed to the which can be attributed to the variability of the plasma environment and as a consequence, the induced magnetosphere of Titan. A weak trend between the incoming electron fluxes and the measured drop-out altitudes has also been observed.     

Biological and Physical Clogging in Infiltration Wells: Effects of Well Diameter and Gravel Pack

Gravity-driven infiltration into the shallow subsurface via small-diameter wells (SDWs), i.e., wells with an inner diameter smaller than 7.5 cm (3 inches) and no gravel pack) has proven to be a cost-efficient and flexible tool for managed aquifer recharge (MAR), as it provides relatively high recharge rates with minimal construction effort. SDWs have a significantly smaller open filter area than larger diameter wells with gravel pack, making the infiltration of low-quality waters through these wells more at risk clogging. To investigate their susceptibility for biological and physical clogging, 24 physical models with different well setups were evaluated by infiltrating either nutrient-poor but turbid water or nutrient-rich but clear water. The experiments showed that smaller diameters and the lack of a gravel pack increase the well's susceptibility to both kinds of clogging. However, this effect was observed to be much more pronounced for physical than for biological clogging. Our conclusion is that SDWs show severe disadvantages with respect to the infiltration of highly turbid waters in comparison to large diameter wells with a gravel pack. Nevertheless, this disadvantage is much less severe when it comes to the infiltration of clear but nutrient-rich waters (e.g., treated wastewater). Depending on the economic and geological circumstances of a MAR-project, this disadvantage could be outweighed by the significantly lower construction costs of SDWs.     

Sediment cascades and the entangled relationship between human impact and natural dynamics at the pre-pottery Neolithic site of Göbekli Tepe,Anatolia

This study presents a meta-analysis of radiocarbon ages for the environs of Göbekli Tepe – one of the oldest monumental structures worldwide – using cumulative probability functions to diachronically assess phases of geomorphodynamic activity as controlled by natural or anthropogenic drivers. We employ sediment cascades as a heuristic framework to study the complex responses of the geomorphological system to various triggers at local to supra-regional scales. Possible triggers include climatic variability as documented by supra-regional hydroclimatic proxy data, regional demographic trends, and local to regional socioeconomic developments such as the emergence of sedentism or the introduction and dispersal of livestock herding. Our results show that phases of intensified geomorphodynamic activity occurred between ca. 7.4–7.0 and 5.8–3.3 ka BP. These phases roughly coincide with phases of population growth in southern Turkey and climatic variations in Turkey and the Levant. The phase between ca. 5.8–3.3 ka BP also corresponds to the time when organized agriculture and the seeder plough were introduced. Also, the identified phases are in agreement with the general trend of varying geomorphodynamic activity in the Eastern Mediterranean as driven by human impact and climatic change. However, neither the Younger Dryas–Holocene transition nor the development of herding during the Pre-Pottery Neolithic left a clear signature. We demonstrate how the different depositional environments in the studied landscape compartments vary with respect to their spatiotemporal coverage and discuss challenges when trying to understand processes that once shaped landscapes of past societies. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Hydrological modelling for assessing spatio-temporal groundwater recharge variations in the water-stressed Amathole Water Supply System,Eastern Cape,South Africa

To increase the resilience of regional water supply systems in South Africa in the face of anticipated climatic changes and a constant increase in water demand, water supply sources require diversification. Many water-stressed metropolitan regions in South Africa depend largely on surface water to cover their water demand. While climatic and river discharge data is widely available in these regions, information on groundwater resources – which could support supply source diversification – is scarce. Groundwater recharge is a key parameter that is used to estimate groundwater amounts that can be sustainably exploited at a sub-watershed level. Therefore, the objective of this study was to develop a reliable hydrological modelling routine that enables the assessment of regional spatio-temporal variations of groundwater recharge to discern the most promising areas for groundwater development. Accordingly, we present a semi-distributed hydrological modelling approach that incorporates water balance routines coupled with baseflow modelling techniques to yield spatio-temporal variations of groundwater recharge on a regional level. The approach is demonstrated for the actively managed catchment areas of the Amathole Water Supply System situated in a semi-arid part of the Eastern Cape of South Africa. In the investigated study area, annual groundwater recharge exhibits a high spatio-temporal heterogeneity and is estimated to vary between ~0.5% and 8% of annual rainfall. Despite some uncertainties induced by limited data availability, calibration and validation of the model were found to be satisfactory and yielded model results similar to (point) data of annual groundwater recharge reported in earlier studies. Our approach is therefore found to derive crucial information for efficiently targeting more detailed groundwater exploration studies and could work as a blueprint for orientating groundwater potential exploration in similar environments.     

A catalogue of M51 type galaxy associations

A catalog of 232 apparently interacting galaxy pairs of the M51 class is presented. Catalog members were identified from visual inspection of multi-band images in the IRSA archive. The major findings in the compilation of this catalog are (1) A surprisingly low number of the main galaxies in M51 systems are early type spirals and barred spirals. (2) Over 70% of the main galaxies in M51 systems are 2-armed spirals. (3) Some systems that were classified as M51 types in previous studies are not M51 types as defined in this catalog. There were a number of systems previously classified as M51 systems for which the “companion” is identified as an HII region within the main galaxy or foreground star within the Milky Way. (4) It was found that only 18% of the M51 type companions have redshift measurements in the literature. There is a significant need for spectroscopic study of the companions in order to improve the value of the catalog as a sample for studying the effects of M51 type interaction on galaxy dynamics, morphology, and star formation. Further spectroscopy will also help constrain the statistics of possible chance projections between foreground and background galaxies in the catalog. The catalog also contains over 430 additional systems which are classified as “possible M51” systems. The reasons for classifying certain systems as possible M51 systems are discussed.     

The integrated Sachs–Wolfe effect in cosmologies with coupled dark matter and dark energy

The subject of this paper is the derivation of the integrated Sachs–Wolfe (iSW) effect in cosmologies with coupled dark matter and dark energy fluids. These couplings influence the iSW effect in three ways: the Hubble function assumes a different scaling, the structure growth rate shows a different time evolution and, in addition, the Poisson equation, which relates the density perturbations to fluctuations in the gravitational potential, is changed, due to the violation of the scaling  ρ∝ a ⁻³  of the matter density ρ with scalefactor a . Exemplarily, I derive the iSW spectra for a model in which dark matter decays into dark energy, investigate the influence of the dark matter decay rate and the dark energy equation of state on the iSW signal, and discuss the analogies for gravitational lensing. Quite generally, iSW measurements should reach similar accuracy in determining the dark energy equation of state parameter and the coupling constant.