Dynamical and compositional assessment of near‐Earth object mission targets

Abstract— Using an H‐plot analysis, we identify 234 currently known near‐Earth objects that are accessible for rendezvous with a “best case” delta‐V of less than 7 km/s. We provide a preliminary compositional interpretation and assessment of these potential targets by summarizing the taxonomic properties for 44 objects. Results for one‐half (22) of this sample are based on new spectroscopic measurements presented here. Our approach provides an easy‐to‐update method for giving guidelines to both observers and mission analysts for focusing on objects for which actual mission opportunities are most likely to be found. Observing prospects are presented for categorizing the taxonomic properties of the most accessible targets that are not yet measured.     

Ray tracing in elliptical anisotropic media using the linear traveltime interpolation (LTI) method applied to traveltime seismic tomography

The linear traveltime interpolation (LTI) method is a suitable ray‐tracing technique for modelling first‐arrival times in isotropic media. LTI is extended to elliptical anisotropic media and applied to a tomographic inversion procedure. A theoretical formulation is first derived and then LTI implementation is discussed in terms of source–receiver arrays and cell size. The method is then combined with the tomographic inversion procedure adopted. The matching of the ray tracing with inversion in elliptical anisotropic media posed a double non‐linear problem. Thus two assumptions were made: the velocity in each cell is uniform and the main directions of anisotropy are known. To take into account the geometrical characteristics of the area under investigation (depth and velocity of the weathering, and thickness of the inner media), cells of varying size were considered. No hypothesis was made on anisotropy weakness.
The algorithm was first tested on synthetic models and then applied to a field survey. On comparing the results of the synthetic models and the field survey with those obtained with a linear raypath approximation, it was found that there were fewer data misfits.     

Geoelectrical survey of structures closed in depth by two dipping discontinuities

Alfano's integral equation method, previously used to calculate VES theoretical curves in parallel or diverging dipping planes including an indefinite conductive medium (open structures), was applied to calculate VES theoretical curves in closed structures, i.e. when the conductive medium is limited in depth. The theoretical curves, as calculated for various cases, show slopes steeper than the corresponding two-layers curves and resistivity asymptotic values equal to those of the resistive layer; their trends can be used to evaluate the dips of contacts.Two surveys, carried out on similar superficial situations in the Appennines area, have been interpreted by means of the new theoretical curves, and have shown different structural situations in depth.     

Background simulations for the Large Area Detector onboard LOFT

The Large Observatory For X-ray Timing (LOFT), currently in an assessment phase in the framework the ESA M3 Cosmic Vision programme, is an innovative medium-class mission specifically designed to answer fundamental questions about the behaviour of matter, in the very strong gravitational and magnetic fields around compact objects and in supranuclear density conditions. Having an effective area of ～10 m² at 8 keV, LOFT will be able to measure with high sensitivity very fast variability in the X-ray fluxes and spectra. A good knowledge of the in-orbit background environment is essential to assess the scientific performance of the mission and optimize the design of its main instrument, the Large Area Detector (LAD). In this paper the results of an extensive Geant-4 simulation of the instrumentwillbe discussed, showing the main contributions to the background and the design solutions for its reduction and control. Our results show that the current LOFT/LAD design is expected to meet its scientific requirement of a background rate equivalent to 10 mCrab in 2?30 keV, achieving about 5 mCrab in the most important 2–10 keV energy band. Moreover, simulations show an anticipated modulation of the background rate as small as 10 % over the orbital timescale. The intrinsic photonic origin of the largest background component also allows for an efficient modelling, supported by an in-flight active monitoring, allowing to predict systematic residuals significantly better than the requirement of 1 %, and actually meeting the 0.25 % science goal.     

A large class of non-zonal oceanic flows satisfying the arnold-blumen sufficient condition for stability

Abstract

Recently Andrews has discussed an example of a topographically-forced non-zonal now satisfying the Arnold-Blumen sufficient condition for stability. At large distances from the topographic centre this flow becomes purely zonal and westward. After underlining the richness of solutions of the Andrews model, the present paper goes on to show that Andrews' technique can be applied successfully to a model where the vertical profile of static stability resembles those found in the ocean. In particular we obtain a large class of hydrodynamic stable flows, forced by the bottom topography, for continuously stratified fluids (two layers each with uniform Brunt-Väisälä frequency).     

Discrete mechanics: Some remarks

A detailed investigation about Discrete Mechanics (a numerical method developed by D. Greenspan and R. A. LaBudde) in Celestial Mechanics applications has been carried out. The method shows a remarkable precession of the argument of pericentre, whose value depends mainly on the time step scaling used and on the length of the integration step. It has been found that this behaviour is due to the structure of Discrete Mechanics and that the precession, taken as a function of the length of the integration step, follows with good agreement an inverse square law. The existence of this relationship allowed the development of a modified version of the method: its consistence with respect to the original one is discussed both from an analytical and a computational point of view. Furthermore the Modified Discrete Mechanics allowed the extension of the method to the case of the potential due to an oblate spheroid.     

Novel spaceways for reaching the Moon: an assessment for exploration

The renewed interest of the major space agencies for the exploration of the Moon has made a review of the present/near future scenario and the related accessible mission profiles desirable. In particular the application of the dynamical systems approach to spaceflight dynamics could bring a significant contribution. A simple method for evaluating the efficiency of these novel spaceways for reaching the Moon if compared to more traditional mission profiles is presented and some general considerations on their utilization for automatic precursor missions as well as for the setting up and the maintenance of a Moon Base are discussed.     

Integrated geophysical methods for the characterisation of an archaeological site (Massenzio Basilica — Roman forum, Rome, Italy)

A geophysical study that involved different techniques was carried out with the aim to improve the knowledge of the archaeological site where the Basilica of Maxentius was founded and to discern individual covered structures (foundations).Vertical Electrical Sounding (VES), seismic refraction and Electrical Resistivity Tomography (ERT) studies were performed in the archaeological site. VES and seismic refraction allowed to characterise the main geological formations of the hill where the Basilica was built and to distinguish the concrete floor and backfilling. Electrical data were processed using different algorithms; their results were compared to appraise the inverted models' robustness.ERT inversion algorithms were used to delineate shape and size of a much more complex structure, that were originally expected from archaeological excavation plan. The results of the commercial program were used as a posteriori information to include them in the algorithm proposed by the authors; the sequential use of the programs defined a processing procedure.The integrated use of different geophysical techniques reduced a great deal the intrinsic ambiguities of each method. Direct explorations (boreholes and archaeological excavations) confirmed the geophysical results.     

2D data modelling by electrical resistivity tomography for complex subsurface geology

A new tool for two‐dimensional apparent‐resistivity data modelling and inversion is presented. The study is developed according to the idea that the best way to deal with ill‐posedness of geoelectrical inverse problems lies in constructing algorithms which allow a flexible control of the physical and mathematical elements involved in the resolution. The forward problem is solved through a finite‐difference algorithm, whose main features are a versatile user‐defined discretization of the domain and a new approach to the solution of the inverse Fourier transform. The inversion procedure is based on an iterative smoothness‐constrained least‐squares algorithm. As mentioned, the code is constructed to ensure flexibility in resolution. This is first achieved by starting the inversion from an arbitrarily defined model. In our approach, a Jacobian matrix is calculated at each iteration, using a generalization of Cohn's network sensitivity theorem. Another versatile feature is the issue of introducing a priori information about the solution. Regions of the domain can be constrained to vary between two limits (the lower and upper bounds) by using inequality constraints. A second possibility is to include the starting model in the objective function used to determine an improved estimate of the unknown parameters and to constrain the solution to the above model. Furthermore, the possibility either of defining a discretization of the domain that exactly fits the underground structures or of refining the mesh of the grid certainly leads to more accurate solutions. Control on the mathematical elements in the inversion algorithm is also allowed. The smoothness matrix can be modified in order to penalize roughness in any one direction. An empirical way of assigning the regularization parameter (damping) is defined, but the user can also decide to assign it manually at each iteration. An appropriate tool was constructed with the purpose of handling the inversion results, for example to correct reconstructed models and to check the effects of such changes on the calculated apparent resistivity. Tests on synthetic and real data, in particular in handling indeterminate cases, show that the flexible approach is a good way to build a detailed picture of the prospected area.