Ephemeral rollover points and clinothem evolution in the modern Po Delta based on repeated bathymetric surveys

Reconstructions of ancient delta systems rely typically on a two-dimensional (2D) view of prograding clinothems but may miss their three-dimensional (3D) stratigraphic complexity which can, instead, be best documented on modern delta systems by integrating high-resolution geophysical data, historical cartography, core data and geomorphological reconstructions offshore. We quantitatively compare three precisely positioned, high-resolution multi-beam bathymetry maps in the delta front and pro delta sectors (0.3 to 10 m water depth) of Po di Pila, the most active of the modern Po Delta five branches. By investigating the detailed morphology of the prograding modern Po Delta, we shed new light on the mechanisms that control the topset to foreset transition in clinothems and show the temporal and spatial complexity of a delta and its pro delta slope, under the impact of oceanographic processes. This study documents the ephemeral nature of the rollover point at the transition between sandy topset (fluvial, delta plain to mouth-bar) and muddy seaward-dipping foreset deposits advancing, in this case, in >20 m of water depth. Three multibeam surveys, acquired between 2013 and 2016, document the complexity in space and time of the topset and foreset regions and their related morphology, a diagnostic feature that could not be appreciated using solely 2D, even very high-resolution, seismic profiles. In addition, the comparison of bathymetric surveys gathered with one-year lapses shows the migration of subaqueous sand dunes on the clinothem topset, the formation of ephemeral cut-and-fill features at the rollover point (few m below mean sea level), the presence of collapse depressions derived by sagging of sediments and fluid expulsion (possibly induced by storm waves) on the foreset, and splays of sand likely reflecting gravity flows on the lower foreset. Though the modern Po Delta is anthropogenic in many respects, its subaqueous clinothem can be studied as a scale model for ancient clinothems that are less resolved geometrically and far less constrained chronologically.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

Doublet craters on Venus

Of the impact craters on Earth larger than 20 km in diameter, 10-15% (3 out of 28) are doublets, having been formed by the simultaneous impact of two well-separated projectiles. The most likely scenario for their formation is the impact of well-separated binary asteroids. If a population of binary asteroids is capable of striking the Earth, it should also be able to hit the other terrestrial planets as well. Venus is a promising planet to search for doublet craters because its surface is young, erosion is nearly nonexistent, and its crater population is significantly larger than the Earth's. After a detailed investigation of single craters separated by less than 150 km and “multiple” craters having diameters greater than 10 km, we found that the proportion of doublet craters on Venus is at most 2.2%, significantly smaller than Earth's, although several nearly incontrovertible doublets were recognized. We believe this apparent deficit relative to the Earth's doublet population is a consequence of atmospheric screening of small projectiles on Venus rather than a real difference in the population of impacting bodies. We also examined “splotches,” circular radar reflectance features in the Magellan data. Projectiles that are too small to form craters probably formed these features. After a careful study of these patterns, we believe that the proportion of doublet splotches on Venus (14%) is comparable to the proportion of doublet craters found on Earth (10-15%). Thus, given the uncertainties of interpretation and the statistics of small numbers, it appears that the doublet crater population on Venus is consistent with that of the Earth.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

Error estimates with L1 solutions

A method is given, based on the pseudoinverse of the equations of condition, to obtain error estimates for the solution in the normL ₁ of an over-determined linear system. The computational labor to obtain the errors, while not trivial, is less than that for various competing methods, particularly if there are many more equations of condition than unknowns. The error estimates for anL ₁ solution are substantially larger than those for a least squares solution of the some data. It is suggested that a complete discussion of a linear system include at least bothL ₁ and least squares solutions with their respective errors and the condition number of the linear system.     

Recovering the initial conditions of our local Universe from NOG and PSCz catalogues

We apply the ztrace algorithm to the optical NOG and infrared PSC z galaxy catalogues to reconstruct the pattern of primordial fluctuations that have generated our local Universe. We check that the density fields traced by the two catalogues are well correlated, and consistent with a linear relation [either in δ or in  log (1 +δ)  ] with relative bias (of NOG with respect to PSC z )   b _rel= 1.1 ± 0.1  . The relative bias relation is used to fill the optical zone of avoidance at  | b | < 20°  using the PSC z galaxy density field.
We perform extensive testing on simulated galaxy catalogues to optimize the reconstruction. The quality of the reconstruction is predicted to be good at large scales, up to a limiting wavenumber   k _lim≃ 0.4 h Mpc⁻¹  beyond which all information is lost. We find that the improvement arising from the denser sampling of the optical catalogue is compensated by the uncertainties connected to the larger zone of avoidance.
The initial conditions reconstructed from the NOG catalogue are found (analogously to those from the PSC z ) to be consistent with a Gaussian paradigm. We use the reconstructions to produce sets of initial conditions ready to be used for constrained simulations of our local Universe.     

Longevity of fluorine-bearing tremolite on Venus

There is a general belief that hydrous minerals cannot exist on Venus under current surface conditions. This view was challenged when Johnson and Fegley (2000, Icarus 146, 301-306) showed that tremolite (Ca₂Mg₅Si₈O₂₂(OH)₂), a hydrous mineral, is stable against thermal decomposition at current Venus surface temperatures, e.g., 50% decomposition in 4 Ga at 740 K. To further explore hydrous mineral thermal stability on Venus, we experimentally determined the thermal decomposition kinetics of fluorine-bearing tremolite. Fluor-tremolite is thermodynamically more stable than OH-tremolite and should decompose more slowly. However how much slower was unknown. We measured the decomposition rate of fluorine-bearing tremolite and show that its decomposition is several times to greater than ten times slower than that of OH-tremolite. We also show that F-bearing tremolite is depleted in fluorine after decomposition and that fluorine is lost as a volatile species such as HF gas. If tremolite ever formed on Venus, it would probably also contain fluorine. The exceptional stability of F-bearing tremolite strengthens our conclusions that if hydrous minerals ever formed on Venus, they could still be there today.     

Dynamics of Two Planets in the 3/2 Mean-motion Resonance: Application to the Planetary System of the Pulsar PSR B1257+12

This paper considers the dynamics of two planets, as the planets B and C of the pulsar PSR B1257+12, near a 3/2 mean-motion resonance. A two-degrees-of-freedom model, in the framework of the general three-body planar problem, is used and the solutions are analyzed through surfaces of section and Fourier techniques in the full phase space of the system.On the possibility of existence of a fourth planet in distant orbit, see Wolszczan et al., 2000     

An estimate of the contribution of high nergy cosmic-ray protons to the absorbed dose inventory of a cometary nucleus

We discuss the contribution of cosmic-ray protons at all energies above 1 MeV to the absorbed doses of the surface layers of a comet. Since there exists no calculation which takes into account proton energy losses by means of losses to electrons and nuclear collisions (in a cascade process), and losses due to the low energy end-products of the cascade, we have made a rough estimate of all of these contributions. An analytical formula is proposed that allows a rapid estimate of ranges and the dose absorbed at any depth. We give dose-depth curves for two extreme values of the energy at which nuclear collisions begin to dominate the slowing-down process, and for an intermediate value we display the dose-depth curve down to 20 m from the surface. The relevance of these findings to dosimetry in comets and some alterations of cometary material are considered. The need for improving the analytical expression proposed is stressed.