Laboratory sandbox validation of pollutant source location methods

Inverse methods can be used to recover the pollutant source location from concentration data. In this paper, the relative effectiveness of two proposed methods, simultaneous release function and source location identification (SRSI) and backward probability model based on adjoint state method (BPM-ASM) are evaluated using real data collected by using experimental equipment. The device is a sandbox that reproduces an unconfined aquifer in which all the variables are controlled. A numerical model was calibrated using experimental observations. The SRSI is a stochastic procedure which finds the source location and the release history by means of a Bayesian geostatistical approach (GA). The BPM-ASM provides the backward probability location of the pollutant detected at a monitoring point by means of a reverse transport simulation. The results show that both methods perform well. While the simultaneous release function and SRSI method requires a preliminary delineation of a probable source area and some weak hypotheses about the statistical structure of the unknown release function, the backward probability model requires some hypothesis about the contaminant release time. A case study was performed using two observation points only, and despite the scarcity of data, both methodologies were able to accurately reconstruct the true source location. The GA has the advantage to recover the release history function too, whilst the backward probability model works well with fewer data. If there are many observations, both methodologies may be computationally heavy. A transfer function approach has been adopted for the numerical definition of the sensitivity matrix in the SRSI method. The reliability of the experimental equipment was tested in previous laboratory works, conducted under several different conditions.     

Application of wavelet transform analysis to landslide generated waves

The aim of this paper is to illustrate the results of a preliminary study on aerial landslide-generated waves, which has been mainly intended to establish a method for analyzing water surface records. Some simple physical experiments, reproducing the Scott Russell's wave generator, were carried out in a small two-dimensional wave flume; the Wavelet Transform (WT) is applied to analyze wave measurements and it is shown that useful information can be obtained by means of this technique. The celerity of impulsively generated waves, reflection by an overflow structure and seiching phenomena of the flume are studied. A discussion of the results along with some remarks about ongoing research is also given.     

Remote sensing of Mercury's magnetospheric plasma environment via energetic neutral atoms imaging

The potentiality of exploring the plasma dynamical behavior of Mercury via energetic neutral atom (ENA) analysis is briefly discussed. The analysis of the Hermean ENA would give clues to the global magnetospheric configuration, clarifying crucial issues like the existence of radiation belts. In order to demonstrate the feasibility of ENA imaging in the Hermean environment, an empirical model of the equatorial proton distributions, originally developed for the Earth's magnetosphere, is here modified in order to be applicable to the expected conditions, according to the present knowledge of the planet. Simulations of the ENA radiation are then presented and discussed. It is shown that the simulated radiation is significantly high, especially when pointing the field of view towards the Hermean surface. The way of determining the proton flux distributions and the cross-tail electric field characteristics by using the ENA signal is discussed. Finally, a preliminary design of an ENA detector suitable for the ESA mission to Mercury ‘BepiColombo’ is proposed and discussed.     

Catastrophic fragmentation and formation of families: Preliminary results from a new numerical model

Preliminary results of an improved version of the semiempirical model for catastrophic break up processes developed by Paolicchi et al., (1989) are presented. Among the several changes with respect to the old version, the most important seem to be related to the new treatment of gravitational effects, including self-compression and reaccumulation of fragments. In particular, the new model is able to analyze processes involving both cm-sized objects, like those studied by means of laboratory experiments, as well as much larger bodies, for which self-gravitational effects are dominant; moreover, in this latter case the model seems in principle adequate to describe with the same physics very different phenomena, like the formation of plausible asteroid families and the creation of single, rapidly spinning, objects. This fact, if confirmed by refined analyses, may be of high importance for our general understanding of asteroid collisional evolution.     

Symmetries and Rank Deficiency in the Orbit Determination Around Another Planet

We study the possible degeneracies for the normal matrix of the observations from the Earth of the motion of a satellite around a planet, and give the possible solutions to the loss of precision in the orbit determination caused by the rank deficiency. Finally we discuss the methods available to control the instability in the orbit determination resulting from the degeneracy.     

The collapse of primordial stars: importance of radiation

We simulate the collapse of a primordial protostellar cloud by means of a 1D hydrodynamics code accounting for chemical evolution, radiative transfer and radiation pressure. We find that the role of radiation pressure is negligible throughout the whole simulations, i.e. Until shortly after the formation of a central hydrostatic core. We also estimate the luminosity and the spectrum of such collapsing clouds. The luminosity is initially due to a number of H₂ lines and is of the order of 10^33-34 erg s^-1. It then grows to values ≳10³⁶ erg s^-1 by the time the core forms, and results from both HH lines and continuum radiation. This revised version was published online in September 2006 with corrections to the Cover Date.     

Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21 cm line from the earliest phases of star and galaxy formation in the Universe. This 21 cm signal provides a new and unique window both on the time of the formation of the first stars and accreting black holes and the subsequent period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.     

A Minimal Spanning Tree algorithm for source detection in γ-ray images

We developed a source detection algorithm based on the Minimal Spanning Tree (MST), that is a graph-theoretical method useful for finding clusters in a given set of points. This algorithm is applied to γ-ray bi-dimensional images where the points correspond to the arrival direction of photons, and the possible sources are associated with the regions where they clusterize. Some filters to select these clusters and to reduce the spurious detections are introduced. An empirical study of the statistical properties of MST on random fields is carried out in order to derive some criteria to estimate the best filter values. We also introduce two parameters useful to verify the goodness of candidate sources. To show how the MST algorithm works in practice, we present an application to an EGRET observation of the Virgo field, at high Galactic latitude and with a low and rather uniform background, in which several sources are detected.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.