A detrital apatite fission‐track study of the CIROS‐2 sedimentary record: Tracing ice pathways in the Ross Sea area over the last 5 Ma

The Ross Sea is a crucial area to investigate pathways of ice during the Cenozoic as it records the evolution of both the East and West Antarctic Ice Sheets. This work is based on detrital apatite fission track (AFT) data extracted from the sedimentary record of well CIROS‐2, which spans the last 5 Ma. The AFT data show a large range of ages, and most of the grains fit well with two main components that fall between 24 and 42 Ma and between 43 and 70 Ma, whereas the other components are not regularly distributed through the well, thus indicating a mixture of provenance from different areas along the Transantarctic Mountains. As a whole, our work suggests glacial expansion over the McMurdo Sound during the Pliocene, and ice periodically invading and retreating in Pleistocene.     

Optimal low-thrust trajectories to asteroids through an algorithm based on differential dynamic programming

In this paper an optimisation algorithm based on Differential dynamic programming is applied to the design of rendezvous and fly-by trajectories to near Earth objects. Differential dynamic programming is a successive approximation technique that computes a feedback control law in correspondence of a fixed number of decision times. In this way the high dimensional problem characteristic of low-thrust optimisation is reduced into a series of small dimensional problems. The proposed method exploits the stage-wise approach to incorporate an adaptive refinement of the discretisation mesh within the optimisation process. A particular interpolation technique was used to preserve the feedback nature of the control law, thus improving robustness against some approximation errors introduced during the adaptation process. The algorithm implements global variations of the control law, which ensure a further increase in robustness. The results presented show how the proposed approach is capable of fully exploiting the multi-body dynamics of the problem; in fact, in one of the study cases, a fly-by of the Earth is scheduled, which was not included in the first guess solution.     

Detection of Flow Distortions and Systematic Errors in Sonic Anemometry Using the Planar Fit Method

We describe the coordinate transformations that can be used to convert the velocity components measured by a set of sonic anemometers with time-dependent tilt fluctuations into a single, time-independent coordinate system. By applying the planar fit method (PFM) to each anemometer dataset, it is possible, for planar flows, to locate the flow plane at each measurement point and compare its orientation with the topography. Installation on a ship is also considered. An application of this method to intercomparison data has led to the detection of an instrument error due to a misalignment between the assembly of the sonic transducers and the anemometer pedestal. If this error occurs, pedestal levelling does not guarantee that measurements are unbiased. A correction method is proposed and the results of two experiments are shown. Flow planarity at different levels and flow distortion caused by the mast are highlighted. The influence of the error on the evaluation of the Reynolds stresses using PFM or the double rotation method and the triple rotation method is discussed and the tilt corrected stresses calculated using the three methods compared.     

Changes of the susceptibility to lava flow invasion induced by morphological modifications of an active volcano: the case of Mount Etna,Italy

Changes of the susceptibility to lava flow invasion at Mount Etna are quantified by using lava flow simulations on four Digital Elevation Models documenting the morphostructural modifications of the volcano in the time interval 1986–2007. The probabilistic code DOWNFLOW is used to derive the areas invaded by several thousands of lava flows obtaining, for each DEM, maps of the susceptibility to lava flow invasion and of the lava flow hazard. These maps show, for the first time, the evolution of these surficial properties with time, and render a quantitative image of the effects of topographic changes on the preferential lava flow drainage paths. The results illustrate how the emplacement of new lava flows and the growth of scoria cones affect the probability of inundation by lava flows. We conclude that the persistent activity of this volcano requires a frequent updating of the topography for a reliable lava flow hazard assessment.     

Land uplift due to subsurface fluid injection

The subsurface injection of fluid (water, gas, vapour) occurs worldwide for a variety of purposes, e.g. to enhance oil production (EOR), store gas in depleted gas/oil fields, recharge overdrafted aquifer systems (ASR), and mitigate anthropogenic land subsidence. Irrespective of the injection target, some areas have experienced an observed land uplift ranging from a few millimetres to tens of centimetres over a time period of a few months to several years depending on the quantity and spatial distribution of the fluid used, pore pressure increase, geological setting (depth, thickness, and area extent), and hydro-geomechanical properties of the injected formation. The present paper reviews the fundamental geomechanical processes that govern land upheaval due to fluid injection in the subsurface and presents a survey of some interesting examples of anthropogenic uplift measured in the past by the traditional levelling technique and in recent times with the aid of satellite technology. The examples addressed include Long Beach, Santa Clara Valley, and Santa Ana basin, California; Las Vegas Valley, Nevada; Cold Lake and other similar sites, Canada; Tokyo and Osaka, Japan; Taipei, Taiwan; Krechba, Algeria; Upper Palatinate, Germany; Chioggia and Ravenna, Italy.     

Hazard assessment at Mount Etna using a hybrid lava flow inundation model and satellite-based land classification

Using a lava flow emplacement model and a satellite-based land cover classification, we produce a map to allow assessment of the type and quantity of natural, agricultural and urban land cover at risk from lava flow invasion. The first step is to produce lava effusion rate contours, i.e., lines linking distances down a volcano??s flank that a lava flow will likely extend if fed at a given effusion rate from a predetermined vent zone. This involves first identifying a vent mask and then running a downhill flow path model from the edge of every pixel around the vent mask perimeter to the edge of the DEM. To do this, we run a stochastic model whereby the flow path is projected 1,000 times from every pixel around the vent mask perimeter with random noise being added to the DEM with each run so that a slightly different flow path is generated with each run. The FLOWGO lava flow model is then run down each path, at a series of effusion rates, to determine likely run-out distance for channel-fed flow extending down each path. These results are used to plot effusion rate contours. Finally, effusion rate contours are projected onto a land classification map (produced from an ASTER image of Etna) to assess the type and amount of each land cover class falling within each contour. The resulting maps are designed to provide a quick look-up capability to assess the type of land at risk from lava extending from any location at a range of likely effusion rates. For our first (2,000 m) vent zone case used for Etna, we find a total of area of ~680 km² is at risk from flows fed at 40 m³ s^?1, of which ~6 km² is urban, ~150 km² is agriculture and ~270 km² is grass/woodland. The model can also be run for specific cases, where we find that Etna??s 1669 vent location, if active today, would likely inundate almost 11 km² of urban land, as well as 15.6 km² of agricultural land, including 9.5 km² of olive groves and 5.2 km² of vineyards and fruit/nut orchards.     

Sun–Earth L_{1} and L_{2} to Moon transfers exploiting natural dynamics

This paper examines the design of transfers from the Sun–Earth libration orbits, at the \(L_{1}\) and \(L_{2}\) points, towards the Moon using natural dynamics in order to assess the feasibility of future disposal or lifetime extension operations. With an eye to the probably small quantity of propellant left when its operational life has ended, the spacecraft leaves the libration point orbit on an unstable invariant manifold to bring itself closer to the Earth and Moon. The total trajectory is modeled in the coupled circular restricted three-body problem, and some preliminary study of the use of solar radiation pressure is also provided. The concept of survivability and event maps is introduced to obtain suitable conditions that can be targeted such that the spacecraft impacts, or is weakly captured by, the Moon. Weak capture at the Moon is studied by method of these maps. Some results for planar Lyapunov orbits at \(L_{1}\) and \(L_{2}\) are given, as well as some results for the operational orbit of SOHO.     

Improved Leap—Frog Symplectic Integrators for Orbits of Small Eccentricity in the Perturbed Kepler Problem

I have improved the precision of the leap–frog symplectic integrators for perturbed Kepler problems at small eccentricities, without significant loss of CPU time. The integration scheme proposed is competitive, in some situations, with the so-called mixed variable integrators.     

Analysis of the Chaotic Behaviour of Orbits Diffusing along the Arnold Web

In a previous work [Guzzo et al. DCDS B 5, 687–698 (2005)] we have provided numerical evidence of global diffusion occurring in slightly perturbed integrable Hamiltonian systems and symplectic maps. We have shown that even if a system is sufficiently close to be integrable, global diffusion occurs on a set with peculiar topology, the so-called Arnold web, and is qualitatively different from Chirikov diffusion, occurring in more perturbed systems. In the present work we study in more detail the chaotic behaviour of a set of 90 orbits which diffuse on the Arnold web. We find that the largest Lyapunov exponent does not seem to converge for the individual orbits while the mean Lyapunov exponent on the set of 90 orbits does converge. In other words, a kind of average mixing characterizes the diffusion. Moreover, the Local Lyapunov Characteristic Numbers (LLCNs), on individual orbits appear to reflect the different zones of the Arnold web revealed by the Fast Lyapunov Indicator. Finally, using the LLCNs we study the ergodicity of the chaotic part of the Arnold web.     

Is the universe homogeneous? (On large scales)

I critically discuss in a pedagogical and phenomenological way a few crucial tests challenging the recent claims by Pietronero and collaborators that there is no evidence from available galaxy catalogues that the Universe is actually homogeneous above a certain scale. In a series of papers, these authors assert that observations are consistent with a fractal distribution of objects extending to the limit of the present data. I show that while galaxies are indeed clustered in a scale-free (fractal) way on small and intermediate scales, this behaviour does not continue indefinitely. Although the specific wavelength at which the galaxy distribution apparently turns to homogeneity is dangerously close to the size of the largest samples presently available, there are serious hints suggesting that this turnover is real and that its effects are detected in the behaviour of statistical estimators. The most recent claims of a continuing fractal hierarchy up to scales of several hundreds Megaparsecs seem to be ascribable to the use of incomplete samples or to an improper treatment of otherwise high-quality data sets. The fractal perspective, nevertheless, represents a fruitful way to look at the clustering properties of galaxies, when properly coupled to the traditional gravitational instability scenario. In the last part of this paper I will try to clarify, at a very simple level, some of the confusion existing on the actual scaling properties of the galaxy distribution, and discuss how these can provide hints on the evolution of the large-scale structure of the Universe.