New fully kinetic model for the study of electric potential, plasma, and dust above lunar landscapes

We have developed a new fully kinetic electrostatic simulation, HYBes, to study how the lunar landscape affects the electric potential and plasma distributions near the surface and the properties of lifted dust. The model embodies new techniques that can be used in various types of physical environments and situations. We demonstrate the applicability of the new model in a situation involving three charged particle species, which are solar wind electrons and protons, and lunar photoelectrons. Properties of dust are studied with test particle simulations by using the electric fields derived from the HYBes model. Simulations show the high importance of the plasma and the electric potential near the surface. For comparison, the electric potential gradients near the landscapes with feature sizes of the order of the Debye length are much larger than those near a flat surface at different solar zenith angles. Furthermore, dust test particle simulations indicate that the landscape relief influences the dust location over the surface. The study suggests that the local landscape has to be taken into account when the distributions of plasma and dust above lunar surface are studied. The HYBes model can be applied not only at the Moon but also on a wide range of airless planetary objects such as Mercury, other planetary moons, asteroids, and nonactive comets.

Is Vesta an intact and pristine protoplanet?

It is difficult to find a Vesta model of iron core, pyroxene and olivine-rich mantle, and HED crust that can match the joint constraints of (a) Vesta's density and core size as reported by the Dawn spacecraft team; (b) the chemical trends of the HED meteorites, including the depletion of sodium, the FeO abundance, and the trace element enrichments; and (c) the absence of exposed mantle material on Vesta's surface, among Vestoid asteroids, or in our collection of basaltic meteorites. These conclusions are based entirely on mass-balance and density arguments, independent of any particular formation scenario for the HED meteorites themselves. We suggest that Vesta either formed from source material with non-chondritic composition or underwent after its formation a radical physical alteration, possibly caused by collisional processes, that affected its global composition and interior structure. (C) 2015 Elsevier Inc. All rights reserved.

SPH calculations of asteroid disruptions: The role of pressure dependent failure models

We present recent improvements of the modeling of the disruption of strength dominated bodies using the Smooth Particle Hydrodynamics (SPH) technique. The improvements include an updated strength model and a friction model, which are successfully tested by a comparison with laboratory experiments. In the modeling of catastrophic disruptions of asteroids, a comparison between old and new strength models shows no significant deviation in the case of targets which are initially non-porous, fully intact and have a homogeneous structure (such as the targets used in the study by Benz and Asphaug, 1999). However, for many cases (e.g. initially partly or fully damaged targets and rubble-pile structures) we find that it is crucial that friction is taken into account and the material has a pressure dependent shear strength. Our investigations of the catastrophic disruption threshold (27, as a function of target properties and target sizes up to a few 100 km show that a fully damaged target modeled without friction has a Q(D)*:, which is significantly (5-10 times) smaller than in the case where friction is included. When the effect of the energy dissipation due to compaction (pore crushing) is taken into account as well, the targets become even stronger (Q(D)*; is increased by a factor of 2-3). On the other hand, cohesion is found to have an negligible effect at large scales and is only important at scales less than or similar to 1 km. Our results show the relative effects of strength, friction and porosity on the outcome of collisions among small (less than or similar to 1000 km) bodies. These results will be used in a future study to improve existing scaling laws for the outcome of collisions (e.g. Leinhardt and Stewart, 2012). (C) 2014 Elsevier Ltd. All rights reserved.

Modeling asteroid collisions and impact processes

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.

(Table 3) Stable isotope composition of sof t tissues of specimens from the Sea of Okhotsk

We discovered and investigated several cold-seep sites in four depth zones of the Sea of Okhotsk off Northeast Sakhalin: outer shelf (160-250 m), upper slope (250-450 m), intermediate slope (450-800 m), and Derugin Basin (1450-1600 m). Active seepage of free methane or methane-rich fluids was detected in each zone. However, seabed photography and sampling revealed that the number of chemoautotrophic species decreases dramatically with decreasing water depth. At greatest depths in the Derugin Basin, the seeps were inhabited by bacterial mats and bivalves of the families Vesicomyidae (Calyptogena aff. pacifica, C. rectimargo, Archivesica sp.), Solemyidae (Acharax sp.) and Thyasiridae (Conchocele bisecta). In addition, pogonophoran tubeworms of the family Sclerolinidae were found in barite edifices. At the shallowest sites, on the shelf at 160 m, the seeps lack chemoautotrophic macrofauna; their locations were indicated only by the patchy occurrence of bacterial mats. Typical seep-endemic metazoans with chemosynthetic symbionts were confined to seep sites at depths below 370 m. A comparative analysis of the structure of seep and background communities suggests that differences in predation pressure may be an important determinant of this pattern. The abundance of predators such as carnivorous brachyurans and asteroids, which can invade seeps from adjacent habitats and efficiently prey on sessile seep bivalves, decreased very pronouncedly with depth. We conclude from the obvious correlation with the conspicuous pattern in the distribution of seep assemblages that, on the shelf and at the upper slope, predator pressure may be high enough to effectively impede any successful settlement of viable populations of seep-endemic metazoans. However, there was also evidence that other depth-related factors, such as bottom-water current, sedimentary regimes, oxygen concentrations and the supply of suitable settling substrates, may additionally regulate the distribution of seep fauna in the area. As a consequence of the pronounced pattern in the distribution of seep communities, their ecological significance as food sources of surrounding background fauna increased with water depth. Isotopic analyses suggest that in the Derugin Basin seep colonists feed on chemoautotrophic seep organisms, either directly or by preying on metazoans with chemosynthetic symbionts. In contrast, seep organisms apparently do not contribute to the nutrition of the adjacent background fauna on the shelf and at the slope. In this area, elevated epifaunal abundances at seep sites were caused primarily by the availability of suitable settling substrates rather than by an enrichment of food supply.

He isotopes in Late Cretaceous sediments

Late Cretaceous (100-73 Ma) pelagic limestones were measured for helium concentration and isotopic composition to characterize the interplanetary dust flux using 3He as a tracer. In the Bottaccione section near Gubbio, Italy, three intervals of elevated 3He concentration were detected: K1 in the Campanian stage at ~79 Ma, K2 in the Santonian stage at ~ 85 Ma, and K3 in the Turonian stage at ~91 Ma. All three of these episodes are associated with high 3He/4He and 3He/non-carbonate ratios, consistent with their derivation from an enhanced extraterrestrial 3He flux rather than decreased carbonate sedimentation or dissolution. While K2 is modest in magnitude and duration and thus is of limited significance, K1 and K3 are each identified by a few myr interval with an ~4-fold enhancement in mean 3He flux compared with pre-event levels. Samples from ODP Hole 762C in the Indian Ocean spanning both K2 and K3 (93-83 Ma) confirm the presence of a peak in the Turonian stage, suggesting that K3 is a global event. The K1 and K3 3He events are similar in most respects to the two peaks previously detected in the Cenozoic, suggesting a similar origin. These have been attributed to a major asteroid collision in the Late Miocene and to a shower of either comets or asteroids in the Late Eocene. Based on the age and temporal evolution of K1, we suggest that it most likely records the collision which produced the Baptistina asteroid family independently dated at ~80 Ma. The K3 event is less easily explained. It is characterized by an unusually spiky and erratic temporal progression, suggesting an unusual abundance of very 3He rich particles not previously seen in the sedimentary 3He record. We suggest this episode arises either from a comet shower or from an asteroid shower possibly associated with dust-producing lunar impacts.

Soft-bottom macrobenthic faunal associations in the southern Chilean glacial fjord complex

Macrobenthic associations were investigated at 29 sampling stations with a semi-quantitative Agassiz trawl, ranging from the South Patagonian Icefield to the Straits of Magellan in the South Chilean fjord system. A total of 1,895 individuals belonging to 131 species were collected. 19 species belong to colonial organisms, mainly Bryozoa (17 species) and Octocorallia (2 species). The phylum Echinodermata was the most diverse in species number (47 species), with asteroids (25 species) and ophiuroids (13 species) being the best represented within this taxon. Polychaeta was the second dominant group in terms of species richness (46 species). Multidimensional scaling ordination (MDS) separated two station groups, one related to fjords and channels off the South Patagonian Icefield and the second one to stations surrounding the Straits of Magellan. 45 species account for 90% of the dissimilarity between these two groups. These differences can mainly be explained by the influence of local environmental conditions determined by processes closely related to the pres- ence/absence of glaciers. Abiotic parameters such as water depth, type of sediment and chemical features of the superficial sediment were not correlated with the numbers of individuals caught by the Agassiz trawl in each group of sampling stations.

Global Scale Impacts

Global scale impacts modify the physical or thermal state of a substantial fraction of a target asteroid. Specific effects include accretion, family formation, reshaping, mixing and layering, shock and frictional heating, fragmentation, material compaction, dilatation, stripping of mantle and crust, and seismic degradation. Deciphering the complicated record of global scale impacts, in asteroids and meteorites, will lead us to understand the original planet-forming process and its resultant populations, and their evolution in time as collisions became faster and fewer. We provide a brief overview of these ideas, and an introduction to models.

Collisional Formation and Modeling of Asteroid Families

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.

The dust environment of comet 67P/Churyumov-Gerasimenko from Rosetta OSIRIS and VLT observations in the 4.5 to 2.9 AU heliocentric distance range inbound

Context. The ESA Rosetta spacecraft, currently orbiting around cornet 67P/Churyumov-Gerasimenko, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. Aims. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. Methods. A Monte Carlo dust tail code, which has already been used to characterise the dust environments of several comets and active asteroids, has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from. Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. Results. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is alpha = -3, instead of alpha = -2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 AU, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 AU. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

A new approach on the long term dynamics of neo's under Yarkorsky effect

To calculate the force associated with the Yarkovsky effect the temperature distribution on the surface of the asteroid should be determined; it depends on the asteroid orbit, size and shape, spin axis orientation and period, mass, density of surface layers, albedo, thermal conductivity, capacity and IR emissivity of the material. The uncertainty of many of these parameters invites to develop simplified methods to calculate the influence of the Yarkovsky effect on long term dynamics of asteroids. In this paper we present one of this method based in a special perturbation procedure developed in our group.

The ion beam shepherd: A new concept for asteroid deflection

A novel slow push asteroid deflection strategy has been recently proposed in which an Earth threatening asteroid can be deflected by exploiting the momentum transmitted by a collimated beam of quasi-neutral plasma impinging against the asteroid surface. The beam can be generated with state-of-the art ion engines from a hovering spacecraft with no need for physical attachment or gravitational interaction with the celestial body. The spacecraft, placed at a distance of a few asteroid diameters, would need an ion thruster pointed at the asteroid surface as well as a second propulsion system to compensate for the ion engine reaction and keep the distance between the asteroid and the shepherd satellite constant throughout the deflection phase. A comparison in terms of required spacecraft mass per total imparted deflection impulse shows that the method outperforms the gravity tractor concept by more than one order of magnitude for asteroids up to about 200 m diameter. The two methods would yield comparable performance for asteroids larger than about 2 km

High Fidelity Models for Near-Earth Object Dynamics

Motivado por los últimos hallazgos realizados gracias a los recientes avances tecnológicos y misiones espaciales, el estudio de los asteroides ha despertado el interés de la comunidad científica. Tal es así que las misiones a asteroides han proliferado en los últimos años (Hayabusa, Dawn, OSIRIX-REx, ARM, AIMS-DART, ...) incentivadas por su enorme interés científico. Los asteroides son constituyentes fundamentales en la evolución del Sistema Solar, son además grandes concentraciones de valiosos recursos naturales, y también pueden considerarse como objectivos estratégicos para la futura exploración espacial. Desde hace tiempo se viene especulando con la posibilidad de capturar objetos próximos a la Tierra (NEOs en su acrónimo anglosajón) y acercarlos a nuestro planeta, permitiendo así un acceso asequible a los mismos para estudiarlos in-situ, explotar sus recursos u otras finalidades. Por otro lado, las asteroides se consideran con frecuencia como posibles peligros de magnitud planetaria, ya que impactos de estos objetos con la Tierra suceden constantemente, y un asteroide suficientemente grande podría desencadenar eventos catastróficos. Pese a la gravedad de tales acontecimientos, lo cierto es que son ciertamente difíciles de predecir. De hecho, los ricos aspectos dinámicos de los asteroides, su modelado complejo y las incertidumbres observaciones hacen que predecir su posición futura con la precisión necesaria sea todo un reto. Este hecho se hace más relevante cuando los asteroides sufren encuentros próximos con la Tierra, y más aún cuando estos son recurrentes. En tales situaciones en las cuales fuera necesario tomar medidas para mitigar este tipo de riesgos, saber estimar con precisión sus trayectorias y probabilidades de colisión es de una importancia vital. Por ello, se necesitan herramientas avanzadas para modelar su dinámica y predecir sus órbitas con precisión, y son también necesarios nuevos conceptos tecnológicos para manipular sus órbitas llegado el caso. El objetivo de esta Tesis es proporcionar nuevos métodos, técnicas y soluciones para abordar estos retos. Las contribuciones de esta Tesis se engloban en dos áreas: una dedicada a la propagación numérica de asteroides, y otra a conceptos de deflexión y captura de asteroides. Por lo tanto, la primera parte de este documento presenta novedosos avances de apliación a la propagación dinámica de alta precisión de NEOs empleando métodos de regularización y perturbaciones, con especial énfasis en el método DROMO, mientras que la segunda parte expone ideas innovadoras para la captura de asteroides y comenta el uso del “ion beam shepherd” (IBS) como tecnología para deflectarlos. Abstract Driven by the latest discoveries enabled by recent technological advances and space missions, the study of asteroids has awakened the interest of the scientific community. In fact, asteroid missions have become very popular in the recent years (Hayabusa, Dawn, OSIRIX-REx, ARM, AIMS-DART, ...) motivated by their outstanding scientific interest. Asteroids are fundamental constituents in the evolution of the Solar System, can be seen as vast concentrations of valuable natural resources, and are also considered as strategic targets for the future of space exploration. For long it has been hypothesized with the possibility of capturing small near-Earth asteroids and delivering them to the vicinity of the Earth in order to allow an affordable access to them for in-situ science, resource utilization and other purposes. On the other side of the balance, asteroids are often seen as potential planetary hazards, since impacts with the Earth happen all the time, and eventually an asteroid large enough could trigger catastrophic events. In spite of the severity of such occurrences, they are also utterly hard to predict. In fact, the rich dynamical aspects of asteroids, their complex modeling and observational uncertainties make exceptionally challenging to predict their future position accurately enough. This becomes particularly relevant when asteroids exhibit close encounters with the Earth, and more so when these happen recurrently. In such situations, where mitigation measures may need to be taken, it is of paramount importance to be able to accurately estimate their trajectories and collision probabilities. As a consequence, advanced tools are needed to model their dynamics and accurately predict their orbits, as well as new technological concepts to manipulate their orbits if necessary. The goal of this Thesis is to provide new methods, techniques and solutions to address these challenges. The contributions of this Thesis fall into two areas: one devoted to the numerical propagation of asteroids, and another to asteroid deflection and capture concepts. Hence, the first part of the dissertation presents novel advances applicable to the high accuracy dynamical propagation of near-Earth asteroids using regularization and perturbations techniques, with a special emphasis in the DROMO method, whereas the second part exposes pioneering ideas for asteroid retrieval missions and discusses the use of an “ion beam shepherd” (IBS) for asteroid deflection purposes.

Singularities in DROMO formulation. Analysis of deep flybys

The singularities in Dromo are characterized in this paper, both from an analytical and a numerical perspective. When the angular momentum vanishes, Dromo may encounter a singularity in the evolution equations. The cancellation of the angular momentum occurs in very speci?c situations and may be caused by the action of strong perturbations. The gravitational attraction of a perturbing planet may lead to rapid changes in the angular momentum of the particle. In practice, this situation may be encountered during deep planetocentric ?ybys. The performance of Dromo is evaluated in di?erent scenarios. First, Dromo is validated for integrating the orbit of Near Earth Asteroids. Resulting errors are of the order of the diameter of the asteroid. Second, a set of theoretical ?ybys are designed for analyzing the performance of the formulation in the vicinity of the singularity. New sets of Dromo variables are proposed in order to minimize the dependency of Dromo on the angular momentum. A slower time scale is introduced, leading to a more stable description of the ?yby phase. Improvements in the overall performance of the algorithm are observed when integrating orbits close to the singularity.

First EURONEAR NEA discoveries from La Palma using the INT

Since 2006, the European Near Earth Asteroids Research (EURONEAR) project has been contributing to the research of near-Earth asteroids (NEAs) within a European network. One of the main aims is the amelioration of the orbits of NEAs, and starting in 2014 February we focus on the recovery of one-opposition NEAs using the Isaac Newton Telescope (INT) in La Palma in override mode. Part of this NEA recovery project, since 2014 June EURONEAR serendipitously started to discover and secure the first NEAs from La Palma and using the INT, thanks to the teamwork including amateurs and students who promptly reduce the data, report discoveries and secure new objects recovered with the INT and few other telescopes from the EURONEAR network. Five NEAs were discovered with the INT, including 2014 LU14, 2014 NL52 (one very fast rotator), 2014 OL339 (the fourth known Earth quasi-satellite), 2014 SG143 (a quite large NEA), and 2014 VP. Another very fast moving NEA was discovered but was unfortunately lost due to lack of follow-up time. Additionally, another 14 NEA candidates were identified based on two models, all being rapidly followed-up using the INT and another 11 telescopes within the EURONEAR network. They include one object discovered by Pan-STARRS, two Mars crossers, two Hungarias, one Jupiter trojan, and other few inner main belt asteroids (MBAs). Using the INT and Sierra Nevada 1.5 m for photometry, then the Gran Telescopio de Canarias for spectroscopy, we derived the very rapid rotation of 2014 NL52, then its albedo, magnitude, size, and its spectral class. Based on the total sky coverage in dark conditions, we evaluate the actual survey discovery rate using 2-m class telescopes. One NEA is possible to be discovered randomly within minimum 2.8 deg2 and maximum 5.5 deg2. These findings update our past statistics, being based on double sky coverage and taking into account the recent increase in discovery.