Families among high-inclination asteroids

We present a new classification of families identified among the population of high-inclination asteroids. We computed synthetic proper elements for a sample of 18,560 numbered and multi-opposition objects having sine of proper inclination greater than 0.295. We considered three zones at different heliocentric distances (inner, intermediate and outer region) and used the standard approach based on the Hierarchical Clustering Method (HCM) to identify families in each zone. In doing so, we used slightly different approach with respect to previously published methodologies, to achieve a more reliable and robust classification. We also used available SDSS color data to improve membership and identify likely family interlopers. We found a total of 38 families, as well as a significant number of clumps and clusters deserving further investigation.     

Neural network interpolation of the magnetic field for the LISA Pathfinder Diagnostics Subsystem

LISA Pathfinder is a science and technology demonstrator of the European Space Agency within the framework of its LISA mission, which aims to be the first space-borne gravitational wave observatory. The payload of LISA Pathfinder is the so-called LISA Technology Package, which is designed to measure relative accelerations between two test masses in nominal free fall. Its disturbances are monitored and dealt by the diagnostics subsystem. This subsystem consists of several modules, and one of these is the magnetic diagnostics system, which includes a set of four tri-axial fluxgate magnetometers, intended to measure with high precision the magnetic field at the positions of the test masses. However, since the magnetometers are located far from the positions of the test masses, the magnetic field at their positions must be interpolated. It has been recently shown that because there are not enough magnetic channels, classical interpolation methods fail to derive reliable measurements at the positions of the test masses, while neural network interpolation can provide the required measurements at the desired accuracy. In this paper we expand these studies and we assess the reliability and robustness of the neural network interpolation scheme for variations of the locations and possible offsets of the magnetometers, as well as for changes in environmental conditions. We find that neural networks are robust enough to derive accurate measurements of the magnetic field at the positions of the test masses in most circumstances.     

Searching for curvature pion radiation from protons in strongly magnetized pulsars

By using rather conservative estimates based on the simplest polar cap model, we search the ATNF Pulsar Catalogue for strongly magnetized stars that could accelerate relativistic protons up to the curvature pion production threshold. The best candidate turns out to be the 16 ms pulsar J0537-6910, but the corresponding characteristic parameter χ=a/m _p is yet too small to give origin to observable signals. We show that, for pulsars with period P≈1 ms, a surface polar magnetic field B≈10¹² G is required in order to induce detectable curvature pion radiation from accelerated protons in the magnetosphere. Some other emission processes are also considered.     

Latent Exurban Development: City Expansion Along the Rural-To-Urban Gradient in Growing and Declining Regions of Southern Europe

This paper quantifies the extent of exurban development in Mediterranean Europe. The assessment was carried out by studying changes in the urban-to-rural population density gradient between the years 1950 and 2010. Three of the six urban regions in this study have experienced population growth and moderate urban concentration, while two regions appear to be shifting toward population decline and urban de-concentration after having experienced compact expansion. A phase of recent re-urbanization has been observed in one region. Altogether, these findings indicate a common path of urban expansion among representative Mediterranean regions between 1950 and 1980 while, in the following period, the cities experienced distinct development phases. From this study, we conclude that exurban development is mainly the product of a shift from compact and dense to semi-compact and intermediate-density settlements. [Key words: semi-dense urban growth, density-distance curve, Mediterranean Europe].     

Interactions between the rotifer Euchlanis dilatata and the cladocerans Alona glabra and Macrothrix triserialis in relation to diet type

We studied, separately, the effect of the cladocerans Alona glabra and Macrothrix triserialis on the population growth of the rotifer Euchlanis dilatata using three different food types (Chlorella and Scenedesmus, live or dead, offered on a comparable dry weight basis, 3.16 μg C ml⁻¹). Regardless of the presence or absence of competition, E. dilatata cultured using Chlorella (live or dead) showed a higher population growth as compared to Scenedesmus. Compared to heat-killed C. vulgaris, E. dilatata cultured alone and fed on live algae reached higher abundances. Both cladoceran species had an adverse effect on the population growth of Euchlanis, but did not eliminate the rotifers. Of the two cladocerans, M. triserialis had a far more adverse effect on the population densities of E. dilatata as compared to A. glabra. E. dilatata also influenced the population growth of A. glabra and M. triserialis, the former being more adversely affected than the latter. When the egg ratios of E. dilatata were plotted as a function of population density, a significantly inverse relation was obtained only in treatments containing Chlorella but not for other treatments. Statistically the rate of population increase (r) per day of E. dilatata was adversely affected by the diet type and the presence of competing species. The interaction of diet type X competition was also significant. For A. glabra, the presence of E. dilatata and the diet type, but not their interaction, significantly influenced the r. On the other hand for M. triserialis, only the diet type was significant. Our results suggest that the competitive outcome between rotifers and cladocerans was dependent on diet type.     

Free‐field rotations during earthquakes: Relevance on buildings

Soil rotations around horizontal axes, during an earthquake, are studied through records collected by closely spaced arrays of strong motion accelerometers. The cross power spectrum of accelerations at nearby stations has been generally utilized to describe the spatial distribution of the motion. A number of cross spectra have been obtained during the training of these arrays. To take profit of these elaborations, a mathematical relation is established between the cross power spectrum and the power spectrum of the rotation. Rotation data presented by Liu et al, concerning 52 earthquake records collected at a single station in Taiwan, are compared with rotation data computed according to our procedure. The two series of data are suitably normalized to the peak horizontal acceleration. The data are shown in function of the distance from epicentre. The same ratio, computed according to our procedure, is in good agreement with the average value of these data. Direct measurements and the present approach have lead to evaluations of rotation higher than those predicted by mathematical investigations on the basis of the wave propagation theory, for comparable circumstances. The relevance of this input motion for relatively tall structures is examined, with reference to the structural effects that the horizontal motion concurrently provides. Meaningful will be ranked those effects of the order of magnitude of 20% or higher than those implied by the horizontal excitation. For understanding the relevance on building structures, the procedure has two areas of concern: 1) the coherence implicit in the cross power spectra, which depends on the interpolation process of the original records, collected in the arrays of instruments, and 2) the relative importance of the vertical to the horizontal input motion. As to the second item, the relevance of the rotation component on structures largely depends on the relative importance of the vertical to the horizontal input motion. When the records in an area a few km from the epicenter are considered, the response spectrum of vertical motion is comparable and in some records even higher, than that of horizontal motion, over the entire range of frequencies. This has been observed as well for the 2009 earthquake event of L'Aquila, Italy, and that at the Christchurch (New Zealand) 2011. When the response spectrum of vertical motion is comparable to that of horizontal motion, the effects of rotational motions on most engineering structures can be meaningful. Copyright © 2011 John Wiley & Sons, Ltd.     

Assessment of rockfall susceptibility by integrating statistical and physically-based approaches

In Val di Fassa (Dolomites, Eastern Italian Alps) rockfalls constitute the most significant gravity-induced natural disaster that threatens both the inhabitants of the valley, who are few, and the thousands of tourists who populate the area in summer and winter.To assess rockfall susceptibility, we developed an integrated statistical and physically-based approach that aimed to predict both the susceptibility to onset and the probability that rockfalls will attain specific reaches. Through field checks and multi-temporal aerial photo-interpretation, we prepared a detailed inventory of both rockfall source areas and associated scree-slope deposits. Using an innovative technique based on GIS tools and a 3D rockfall simulation code, grid cells pertaining to the rockfall source-area polygons were classified as active or inactive, based on the state of activity of the associated scree-slope deposits. The simulation code allows one to link each source grid cell with scree deposit polygons by calculating the trajectory of each simulated launch of blocks. By means of discriminant analysis, we then identified the mix of environmental variables that best identifies grid cells with low or high susceptibility to rockfalls. Among these variables, structural setting, land use, and morphology were the most important factors that led to the initiation of rockfalls.We developed 3D simulation models of the runout distance, intensity and frequency of rockfalls, whose source grid cells corresponded either to the geomorphologically-defined source polygons (geomorphological scenario) or to study area grid cells with slope angle greater than an empirically-defined value of 37° (empirical scenario). For each scenario, we assigned to the source grid cells an either fixed or variable onset susceptibility; the latter was derived from the discriminant model group (active/inactive) membership probabilities.Comparison of these four models indicates that the geomorphological scenario with variable onset susceptibility appears to be the most realistic model. Nevertheless, political and legal issues seem to guide local administrators, who tend to select the more conservative empirically-based scenario as a land-planning tool.     

Inverse analysis of stochastic moment equations for transient flow in randomly heterogeneous media

We present a nonlinear stochastic inverse algorithm that allows conditioning estimates of transient hydraulic heads, fluxes and their associated uncertainty on information about hydraulic conductivity (K) and hydraulic head (h  ) data collected in a randomly heterogeneous confined aquifer. Our algorithm is based on Laplace-transformed recursive finite-element approximations of exact nonlocal first and second conditional stochastic moment equations of transient flow. It makes it possible to estimate jointly spatial variations in natural log-conductivity (Y=lnK)$(Y=\ln K)$, the parameters of its underlying variogram, and the variance–covariance of these estimates. Log-conductivity is parameterized geostatistically based on measured values at discrete locations and unknown values at discrete “pilot points”. Whereas prior values of Y at pilot point are obtained by generalized kriging, posterior estimates at pilot points are obtained through a maximum likelihood fit of computed and measured transient heads. These posterior estimates are then projected onto the computational grid by kriging. Optionally, the maximum likelihood function may include a regularization term reflecting prior information about Y. The relative weight assigned to this term is evaluated separately from other model parameters to avoid bias and instability. We illustrate and explore our algorithm by means of a synthetic example involving a pumping well. We find that whereas Y and h can be reproduced quite well with parameters estimated on the basis of zero-order mean flow equations, all model quality criteria identify the second-order results as being superior to zero-order results. Identifying the weight of the regularization term and variogram parameters can be done with much lesser ambiguity based on second- than on zero-order results. A second-order model is required to compute predictive error variances of hydraulic head (and flux) a posteriori. Conditioning the inversion jointly on conductivity and hydraulic head data results in lesser predictive uncertainty than conditioning on conductivity or head data alone.     

Application of a mixing-ratios based formulation to model mixing-driven dissolution experiments

We address the question of how one can combine theoretical and numerical modeling approaches with limited measurements from laboratory flow cell experiments to realistically quantify salient features of complex mixing-driven multicomponent reactive transport problems in porous media. Flow cells are commonly used to examine processes affecting reactive transport through porous media, under controlled conditions. An advantage of flow cells is their suitability for relatively fast and reliable experiments, although measuring spatial distributions of a state variable within the cell is often difficult. In general, fluid is sampled only at the flow cell outlet, and concentration measurements are usually interpreted in terms of integrated reaction rates. In reactive transport problems, however, the spatial distribution of the reaction rates within the cell might be more important than the bulk integrated value. Recent advances in theoretical and numerical modeling of complex reactive transport problems [De Simoni M, Carrera J, Sanchez-Vila X, Guadagnini A. A procedure for the solution of multicomponent reactive transport problems. Water Resour Res 2005;41:W11410. doi: 10.1029/2005WR004056, De Simoni M, Sanchez-Vila X, Carrera J, Saaltink MW. A mixing ratios-based formulation for multicomponent reactive transport. Water Resour Res 2007;43:W07419. doi: 10.1029/2006WR005256] result in a methodology conducive to a simple exact expression for the space–time distribution of reaction rates in the presence of homogeneous or heterogeneous reactions in chemical equilibrium. The key points of the methodology are that a general reactive transport problem, involving a relatively high number of chemical species, can be formulated in terms of a set of decoupled partial differential equations, and the amount of reactants evolving into products depends on the rate at which solutions mix. The main objective of the current study is to show how this methodology can be used in conjunction with laboratory experiments to properly describe the key processes that occur in a complex, geochemically-active system under chemical equilibrium conditions. We model three CaCO₃ dissolution experiments reported in Singurindy et al. [Singurindy O, Berkowitz B, Lowell RP. Carbonate dissolution and precipitation in coastal environments: Laboratory analysis and theoretical consideration. Water Resour Res 2004;40:W04401. doi: 10.1029/2003WR002651, Singurindy O, Berkowitz B, Lowell RP. Correction to Carbonate dissolution and precipitation in coastal environments: laboratory analysis and theoretical consideration. Water Resour Res 2005;41:W11701. doi: 10.1029/2005WR004433], in which saltwater and freshwater were mixed in different proportions. The integrated reaction rate within the cell estimated from the experiments are modeled independently by means of (a) a state-of-the-art reactive transport code, and (b) the uncoupled methodology of [12, 13], both of which use dispersivity as a single, adjustable parameter. The good agreement between the results from both methodologies demonstrates the feasibility of using simple solutions to design and analyze laboratory experiments involving complex geochemical problems.