A model for numerical simulation of nonstationary sonar reverberation using linear spectral prediction

An innovative approach to the numerical generation of nonstationery reverberation time series is presented and demonstrated. The computer simulated reverberation time series are of high quality, in that they are accurate representations of those which would result from an actual sonar system (transmit/receive and horizontal/ vertical beampatterns; pulse type, shape, length, and power; frequency and sampling rate), platform (speed and depth), and environment (wind speed and direction, backscattering strengths, and propagation loss). Volume, surface, and/or bottom reverberation as seen by a multiple beam sonar on a moving platform is generated. The approach utilizes recent developments in linear spectral prediction research in which the spectra of stochastic processes are modeled as rational functions and algorithms are used to efficiently compute optimal estimates of coefficients which specify the spectra. A two-fold sequence is formulated; first, the expected reverberation spectra for all beams are predicted and, second, the stochastic time series are generated from the expected spectra. The expected spectra are predicted using a numerical implementation, referred to as the REVSPEC (reverberation spectrum) model, of a general formulation of Faure, Ol'shevskii, and Middleton. Given the spectra, the Levinson-Durbin method is used to solve the Yule-Walker equations of the autoregressive formulation of linear spectral prediction. The numerical implementation of the approach, referred to as the REVSIM (reverberation simulation) model, produces nonstationary coherent multiple-beam reverberation time series. The formulation of the REVSIM model is presented and typical results given. A comparison is made between the simulation outputs of the REVSIM model and those of the REVGEN (reverberation generator) model, a standard well-accepted time series simulation model, to demonstrate the validity of the new approach.     

Statistical Travel-time Tomography in Terms of Stacking Velocity

— Velocity evaluation is a key step in seismic analysis. The covariance of the true velocity field must be known when interpolating or simulating velocities from well measurements using geostatistical methods. In addition, inversion procedures often require information pertaining to this covariance. Traditionally it has been taken to be the covariance of stacking velocities. We present a simple example to show that this approximation can lead to significant errors. Better methods, such as those of Touati (1996) and Iooss (1998), use the variance of prestack picked travel times as a function of offset to infer that of the velocities. In this paper we extend their results on the estimation of the covariance of the reflected traveltimes, and obtain an explicit expression for the covariance of the square of the stacking slowness as a function of the covariance of the velocities. Although we are not able to invert the formula analytically to yield an explicit estimator for these parameters, the results obtained using it furnish a good and quick estimation of the velocity's covariance. This is illustrated with synthetic examples.     

Rate of Convergence of the Gibbs Sampler in the Gaussian Case

We show that the Gibbs Sampler in the Gaussian case is closely linked to linear fixed point iterations. In fact stochastic linear iterations converge toward a stationary distribution under the same conditions as the classical linear fixed point one. Furthermore the covariance matrices are shown to satisify a related fixed point iteration, and consequently the Gibbs Sampler in the gaussian case corresponds to the classical Gauss-Seidel iterations on the inverse of the covariance matrix, and the stochastic over-relaxed Gauss-Seidel has the same limiting distribution as the Gibbs Sampler. Then an efficient method to simulate a gaussian vector is proposed. Finally numerical investigations are performed to understand the effect of the different strategies such as the initial ordering, the blocking and the updating order for iterations. The results show that in a geostatistical context the rate of convergence can be improved significantly compared to the standard case.     

Mazzite,a new mineral,the natural counterpart of the synthetic zeolite Ω

Mont Semiol (also called Mont Semiouse), Montbrison, Loire, France, is the only place where the zeolite offrettite is known to occur. Now, a new mineral, Mazzite, has been found in the same locality. It occurs in the form of tiny needles, and its properties closely match those of the synthetic “molecular sieve” Ω. The chemical formula is Na_0.03K_1.91Ca_1.35Mg_1.99[Al_9.77Si_26.54O₇₂]·28.03 H₂O. Lattice dimensions are a=18.392 Å, c=7.646 Å, with P 6₃/mmc as most probable space group; D_calc.=2.108. Offrettite and mazzite occur in association with phillipsite, chabazite, calcite and siderite. The mineral has been so named in honour of Prof. Fiorenzo Mazzi.     

Stellerite from Villanova Monteleone,Sardinia

Stellerite (a zeolite) occurs as polycrystal aggregates filling cavities and fractures in the andesite near Villanova Monteleone (Sassari). Principal forms are (100), (010), (001) and (111). The indexed X-ray diffraction powder pattern is reported. The mineral is orthorhombic, with space group Fmmm (or F222 or Fmm2). The unit cell dimensions, chemical formula, density and optical properties are given.     

Radiative energy partition and cloud radiative forcing at a Po valley site

Values of downward and upward flux densities of solar and terrestrial radiation were continuously recorded between 1 December 2001 and 30 November 2002 using a four-components radiometer at S. Pietro Capofiume (SPC) in northern Italy (44°39′N, 11°37′E, alt. 11 m a.m.s.l.), which is characterized by a weakly-reflective surface. The aim of the study was to investigate the effects of clouds on surface radiation balance (SRB); the cloud fraction (N) has been retrieved through the inverted form of the parameterization proposed by Kasten and Czeplak [Solar Energy 24 (1980) 177] and cloud types estimated following the methodology of Duchon and O'Malley [J. Appl. Meteorol. 38 (1999) 132]. The cloud radiative forcing (CRF) was evaluated through the Bintanja and Van den Broeke [Int. J. Climatol. 16 (1996) 1281] formula and then associated with cloud type. Experimental results showed that during the measuring period the net shortwave (Sw) balance always decreased with increasing N, whereas the net longwave (Lw) balance increased in all seasons. The net radiation available at the surface decreased with increasing N in all seasons except in winter, where no significant dependency was detected.The analysis of the cloud radiative forcing indicates that all seasons were characterized by a net cooling of the surface except winter, where clouds seem to have no effects on the surface warming or cooling. Taking into account the dependence on solar radiation cycle, an intercomparison between the retrieved cloud types seems to indicate that the effect of stratus was a slight cooling whereas that of cumulus clouds was a stronger cooling of the surface. On the contrary, cirrus clouds seem to have slight warming effect on the surface.The annual trends of mean monthly values of shortwave and longwave radiation balances confirmed that the measurement site is characterized by a temperate climate. Moreover, in spite of the cooling effect of clouds, a monthly radiative energy surplus is available all year long for surface–atmosphere energy exchanges. The analysis is also instrumental for the detection of SRB variations.     

Comparing landslide inventory maps

Landslide inventory maps are effective and easily understandable products for both experts, such as geomorphologists, and for non experts, including decision-makers, planners, and civil defense managers. Landslide inventories are essential to understand the evolution of landscapes, and to ascertain landslide susceptibility and hazard. Despite landslide maps being compiled every year in the word at different scales, limited efforts are made to critically compare landslide maps prepared using different techniques or by different investigators. Based on the experience gained in 20 years of landslide mapping in Italy, and on the limited literature on landslide inventory assessment, we propose a general framework for the quantitative comparison of landslide inventory maps. To test the proposed framework we exploit three inventory maps. The first map is a reconnaissance landslide inventory prepared for the Umbria region, in central Italy. The second map is a detailed geomorphological landslide map, also prepared for the Umbria region. The third map is a multi-temporal landslide inventory compiled for the Collazzone area, in central Umbria. Results of the experiment allow for establishing how well the individual inventories describe the location, type and abundance of landslides, to what extent the landslide maps can be used to determine the frequency-area statistics of the slope failures, and the significance of the inventory maps as predictors of landslide susceptibility. We further use the results obtained in the Collazzone area to estimate the quality and completeness of the two regional landslide inventory maps, and to outline general advantages and limitations of the techniques used to complete the inventories.     

Quantitative Control of Migration: A Geostatistical Attempt

This paper is devoted to a geostatistical attempt at modeling migration errors when localizing a reflector in the ground. Starting with a probabilistic velocity model and choosing the simple geometrical optics background for the wave propagation in such media, we give the expression of the errors. This may be quantified provided the covariance of the velocity field is known. Variance of arrival times at constant offset is related to the covariance of the velocity field at hand. A practical application is given in the same paragraph. After that we give a typical schema for migration and uncertainty modeling: starting with seismic data, we make the weak seismic inversion. We then obtain the covariance of the velocity field that we use for simulating migration errors. The main issues of this methodology are discussed in the last paragraph.     

Testing geomorphic signal of active normal faulting: The case of the Cittanova Fault (Calabria,southern Italy)

Extracting tectonic signals from the landscape is an important challenge for constraining the style and rate of deformation associated with active faults, especially where their displacement history cannot be independently determined. Based on previous paleoseismological data coupled with new geomorphological field work and ¹⁴C dating of geomorphic markers, we analysed the geomorphic signal of the along‐strike differential throw of the Cittanova Fault in southern Calabria (Italy), the recent activity of which is already well documented and constrained. Through DTM‐derived stream power law parameters (SL and χ), we provide evidence of drainage network disequilibrium and reorganization in response to fault growth and deformation style. Furthermore, a methodological test of the reliability of the χ metric as a proxy for the differential throw along the strike of active normal faults provided good preliminary results, consistent with a strong inverse linear correlation with fault throw. Copyright © 2018 John Wiley & Sons, Ltd.     

Affluence drives the global displacement of land use

Increasing affluence is often postulated as a main driver for the human footprint on biologically productive areas, identified among the main causes of biodiversity loss, but causal relationships are obscured by international trade. Here, we trace the use of land and ocean area through international supply chains to final consumption, modeling agricultural, food, and forestry products on a high level of resolution while also including the land requirements of manufactured goods and services. In 2004, high-income countries required more biologically productive land per capita than low-income countries, but this connection could only be identified when land used to produce internationally traded products was taken into account, because higher-income countries tend to displace a larger fraction of land use. The equivalent land and ocean area footprint of nations increased by a third for each doubling of income, with all variables analyzed on a per capita basis. This increase came largely from imports, which increased proportionally to income. Export depended mostly on the capacity of countries to produce useful biomass, the biocapacity. Our analysis clearly shows that countries with a high biocapacity per capita tend to spare more land for nature. Biocapacity per capita can be increased through more intensive production or by reducing population density. The net displacement of land use from high-income to low-income countries amounted to 6% of the global land demand, even though high-income countries had more land available per capita than low-income countries. In particular, Europe and Japan placed high pressure on ecosystems in lower-income countries.