Random noise attenuation via the randomized canonical polyadic decomposition

Tensor algebra provides a robust framework for multi-dimensional seismic data processing. A low-rank tensor can represent a noise-free seismic data volume. Additive random noise will increase the rank of the tensor. Hence, tensor rank-reduction techniques can be used to filter random noise. Our filtering method adopts the Candecomp/Parafac decomposition to approximates a N-dimensional seismic data volume via the superposition of rank-one tensors. Similar to the singular value decomposition for matrices, a low-rank Candecomp/Parafac decomposition can capture the signal and exclude random noise in situations where a low-rank tensor can represent the ideal noise-free seismic volume. The alternating least squares method is adopted to compute the Candecomp/Parafac decomposition with a provided target rank. This method involves solving a series of highly over-determined linear least-squares subproblems. To improve the efficiency of the alternating least squares algorithm, we uniformly randomly sample equations of the linear least-squares subproblems to reduce the size of the problem significantly. The computational overhead is further reduced by avoiding unfolding and folding large dense tensors. We investigate the applicability of the randomized Candecomp/Parafac decomposition for incoherent noise attenuation via experiments conducted on a synthetic dataset and field data seismic volumes. We also compare the proposed algorithm (randomized Candecomp/Parafac decomposition) against multi-dimensional singular spectrum analysis and classical prediction filtering. We conclude the proposed approach can achieve slightly better denoising performance in terms of signal-to-noise ratio enhancement than traditional methods, but with a less computational cost.     

Reweighting strategies in seismic deconvolution

Reweighting strategies have been widely used to diminish the influence of outliers in inverse problems. In a similar fashion, they can be used to design the regularization term that must be incorporated to solve an inverse problem successfully. Zero-order quadratic regularization, or damped least squares (pre-whitening) is a common procedure used to regularize the deconvolution problem. This procedure entails the definition of a constant damping term which is used to control the roughness of the deconvolved trace. In this paper I examine two different regularization criteria that lead to an algorithm where the damping term is adapted to successfully retrieve a broad-band reflectivity.
Synthetic and field data examples are used to illustrate the ability of the algorithm to deconvolve seismic traces.     

Hydrogeochemical and nitrate isotopic evolution of a semiarid mountainous basin aquifer of glacial-fluvial and paleolacustrine origin (Lake Titicaca,Bolivia): the effects of natural processes and anthropogenic activities

Hydrogeology Journal - A hydrogeochemical and stable isotopic ( $${\updelta }^{15}{\mathrm{N}}_{{\mathrm{NO}}_{3}}$$ and $${\updelta }^{18}{\mathrm{O}}_{{\mathrm{NO}}_{3}}$$ ) multitracer approach...     

The proximal marine record of Somma–Vesuvius volcanic activity in the Naples and Salerno bays,Eastern Tyrrhenian Sea,during the last 3 kyrs

The tephrostratigraphic analysis of nine gravity cores acquired on the continental shelf of the Naples and Salerno bays documents the proximal record of Somma–Vesuvius volcanic activity during the latest Holocene (last 3 kyrs). Five tephra layers from southern Naples Bay and three tephra layers from northern Salerno Bay were recognised in cores and their sedimentologic structures, textural parameters, stratigraphy and major and trace-element composition described.     

L-moments and C-moments

 It is well known that the computation of higher order statistics, like skewness and kurtosis, (which we call C-moments) is very dependent on sample size and is highly susceptible to the presence of outliers. To obviate these difficulties, Hosking (1990) has introduced related statistics called L-moments. We have investigated the relationship of these two measures in a number of different ways. Firstly, we show that probability density functions (pdf ) that are estimated from L-moments are superior estimates to those obtained using C-moments and the principle of maximum entropy. C-moments computed from these pdf's are not however, contrary to what one may have expected, better estimates than those estimated from sample statistics. L-moment derived distributions for field data examples appear to be more consistent sample to sample than pdf 's determined by conventional means. Our observations and conclusions have a significant impact on the use of the conventional maximum entropy procedure which typically uses C-moments from actual data sets to infer probabilities.     

Pre-Quaternary orogen-parallel extension in the Southern Apennine belt, Italy

The Southern Apennine fold and thrust belt differs from other parts of the peri-Tyrrhenian orogen. In most of the peri-Tyrrhenian belt, hinterland extension is oriented at a high-angle to the orogen axis and appears to be related to rifting and formation of oceanic crust within the Tyrrhenian basin. The Southern Apennines share the late-stage development of normal faults related to the opening of the Tyrrhenian Sea, but also experienced an episode of extension parallel to the strike of the tectonic belt. The orogen-parallel extension was apparently formed in response to the increase in length of the deformed belt during arcuation. Arcuation ostensibly was related to asymmetrical rifting in the hinterland, which was greater in the Southern Tyrrhenian Sea than in areas to the north, and proportionately greater shortening in the frontal parts of the southern belt as compared to regions in the north. During arcuation, extension was spatially concentrated within structural domains and was accomplished by displacement on low-angle detachment faults cutting through a previously imbricated thrust stack. During the Miocene-Pliocene, NNW-SSE extension in the interior of the Southern Appennine belt formed coveally with ENE-WSW shortening in the foreland. Longitudinal extension ceased in the Pleistocene, when younger high-angle normal faults formed in response to the easterly migration of Tyrrhenian Sea rifting and NE-SW extension associated with lithospheric stretching.     

New morpho-bathymetric data on marine hazard in the offshore of Gulf of Naples (Southern Italy)

Natural Hazards - Morpho-bathymetric data have been analyzed in order to provide new geological evidence on the Gulf of Naples, i.e., on the southern Ischia slope, on the Naples canyons and on the...