A complex ray-tracing tool for high-frequency mean-field flow interaction effects in jets

This paper presents a complex ray-tracing tool for the calculation of high-frequency Green’s functions in 3D mean field jet flows. For a generic problem, the ray solution suffers from three main deficiencies: multiplicity of solutions, singularities at caustics, and the determining of complex solutions. The purpose of this paper is to generalize, combine and apply existing stationary media methods to moving media scenarios. Multiplicities are dealt with using an equivalent two-point boundary-value problem, whilst non-uniformities at caustics are corrected using diffraction catastrophes. Complex rays are found using a combination of imaginary perturbations, an assumption of caustic stability, and analytic continuation of the receiver curve. To demonstrate this method, the ray tool is compared against a high-frequency modal solution of Lilley’s equation for an off-axis point source. This solution is representative of high-frequency source positions in real jets and is rich in caustic structures. A full utilization of the ray tool is shown to provide excellent results<br/

Non-Parametric Approximations for Anisotropy Estimation in Two-dimensional Differentiable Gaussian Random Fields

Spatially referenced data often have autocovariance functions with elliptical isolevel contours, a property known as geometric anisotropy. The anisotropy parameters include the tilt of the ellipse (orientation angle) with respect to a reference axis and the aspect ratio of the principal correlation lengths. Since these parameters are unknown a priori, sample estimates are needed to define suitable spatial models for the interpolation of incomplete data. The distribution of the anisotropy statistics is determined by a non-Gaussian sampling joint probability density. By means of analytical calculations, we derive an explicit expression for the joint probability density function of the anisotropy statistics for Gaussian, stationary and differentiable random fields. Based on this expression, we obtain an approximate joint density which we use to formulate a statistical test for isotropy. The approximate joint density is independent of the autocovariance function and provides conservative probability and confidence regions for the anisotropy parameters. We validate the theoretical analysis by means of simulations using synthetic data, and we illustrate the detection of anisotropy changes with a case study involving background radiation exposure data. The approximate joint density provides (i) a stand-alone approximate estimate of the anisotropy statistics distribution (ii) informed initial values for maximum likelihood estimation, and (iii) a useful prior for Bayesian anisotropy inference.Comment: 39 pages; 8 figure

Benford's law predicted digit distribution of aggregated income taxes: the surprising conformity of Italian cities and regions

The yearly aggregated tax income data of all, more than 8000, Italian municipalities are analyzed for a period of five years, from 2007 to 2011, to search for conformity or not with Benford's law, a counter-intuitive phenomenon observed in large tabulated data where the occurrence of numbers having smaller initial digits is more favored than those with larger digits. This is done in anticipation that large deviations from Benford's law will be found in view of tax evasion supposedly being widespread across Italy. Contrary to expectations, we show that the overall tax income data for all these years is in excellent agreement with Benford's law. Furthermore, we also analyze the data of Calabria, Campania and Sicily, the three Italian regions known for strong presence of mafia, to see if there are any marked deviations from Benford's law. Again, we find that all yearly data sets for Calabria and Sicily agree with Benford's law whereas only the 2007 and 2008 yearly data show departures from the law for Campania. These results are again surprising in view of underground and illegal nature of economic activities of mafia which significantly contribute to tax evasion. Some hypothesis for the found conformity is presented.Comment: 18 pages, 5 tables, 4 figures, 61 references, To appear in European Physical Journal

Application of Surface wave methods for seismic site characterization

Surface-wave dispersion analysis is widely used in geophysics to infer a shear wave velocity model of the subsoil for a wide variety of applications. A shear-wave velocity model is obtained from the solution of an inverse problem based on the surface wave dispersive propagation in vertically heterogeneous media. The analysis can be based either on active source measurements or on seismic noise recordings. This paper discusses the most typical choices for collection and interpretation of experimental data, providing a state of the art on the different steps involved in surface wave surveys. In particular, the different strategies for processing experimental data and to solve the inverse problem are presented, along with their advantages and disadvantages. Also, some issues related to the characteristics of passive surface wave data and their use in H/V spectral ratio technique are discussed as additional information to be used independently or in conjunction with dispersion analysis. Finally, some recommendations for the use of surface wave methods are presented, while also outlining future trends in the research of this topic

2D characterization of near-surface V P/V S: surface-wave dispersion inversion versus refraction tomography

International audienceThe joint study of pressure (P-) and shear (S-) wave velocities (Vp and Vs ), as well as their ratio (Vp /Vs), has been used for many years at large scales but remains marginal in near-surface applications. For these applications, and are generally retrieved with seismic refraction tomography combining P and SH (shear-horizontal) waves, thus requiring two separate acquisitions. Surface-wave prospecting methods are proposed here as an alternative to SH-wave tomography in order to retrieve pseudo-2D Vs sections from typical P-wave shot gathers and assess the applicability of combined P-wave refraction tomography and surface-wave dispersion analysis to estimate Vp/Vs ratio. We carried out a simultaneous P- and surface-wave survey on a well-characterized granite-micaschists contact at Ploemeur hydrological observatory (France), supplemented with an SH-wave acquisition along the same line in order to compare Vs results obtained from SH-wave refraction tomography and surface-wave profiling. Travel-time tomography was performed with P- and SH- wave first arrivals observed along the line to retrieve Vtomo p and Vtomo s models. Windowing and stacking techniques were then used to extract evenly spaced dispersion data from P-wave shot gathers along the line. Successive 1D Monte Carlo inversions of these dispersion data were performed using fixed Vp values extracted from Vtomo p the model and no lateral constraints between two adjacent 1D inversions. The resulting 1D Vsw s models were then assembled to create a pseudo-2D Vsw s section, which appears to be correctly matching the general features observed on the section. If the pseudo-section is characterized by strong velocity incertainties in the deepest layers, it provides a more detailed description of the lateral variations in the shallow layers. Theoretical dispersion curves were also computed along the line with both and models. While the dispersion curves computed from models provide results consistent with the coherent maxima observed on dispersion images, dispersion curves computed from models are generally not fitting the observed propagation modes at low frequency. Surface-wave analysis could therefore improve models both in terms of reliability and ability to describe lateral variations. Finally, we were able to compute / sections from both and models. The two sections present similar features, but the section obtained from shows a higher lateral resolution and is consistent with the features observed on electrical resistivity tomography, thus validating our approach for retrieving Vp/Vs ratio from combined P-wave tomography and surface-wave profiling

Asynchronous Antarctic and Greenland ice-volume contributions to the last interglacial sea-level highstand

The last interglacial (LIG; ~130 to ~118 thousand years ago, ka) was the last time global sea level rose well above the present level. Greenland Ice Sheet (GrIS) contributions were insufficient to explain the highstand, so that substantial Antarctic Ice Sheet (AIS) reduction is implied. However, the nature and drivers of GrIS and AIS reductions remain enigmatic, even though they may be critical for understanding future sea-level rise. Here we complement existing records with new data, and reveal that the LIG contained an AIS-derived highstand from ~129.5 to ~125 ka, a lowstand centred on 125–124 ka, and joint AIS + GrIS contributions from ~123.5 to ~118 ka. Moreover, a dual substructure within the first highstand suggests temporal variability in the AIS contributions. Implied rates of sea-level rise are high (up to several meters per century; m c−1), and lend credibility to high rates inferred by ice modelling under certain ice-shelf instability parameterisations

Post-orogenic shoshonitic magmas of the Yzerfontein pluton, South Africa: the 'smoking gun' of mantle melting and crustal growth during Cape granite genesis?

The post-orogenic Yzerfontein pluton, in the Saldania Belt of South Africa was constructed through numerous injections of shoshonitic magmas. Most magma compositions are adequately modelled as products of fractionation, but the monzogranites and syenogranites may have a separate origin. A separate high-Mg mafic series has a less radiogenic mantle source. Fine-grained magmatic enclaves in the intermediate shoshonitic rocks are autoliths. The pluton was emplaced between 533 ± 3 and 537 ± 3 Ma (LASF-ICP-MS U–Pb zircon), essentially synchronously with many granitic magmas of the Cape Granite Suite (CGS). Yzerfontein may represent a high-level expression of the mantle heat source that initiated partial melting of the local crust and produced the CGS granitic magmas, late in the Saldanian Orogeny. However, magma mixing is not evident at emplacement level and there are no magmatic kinships with the I-type granitic rocks of the CGS. The mantle wedge is inferred to have been enriched during subduction along the active continental margin. In the late- to post-orogenic phase, the enriched mantle partially melted to produce heterogeneous magma batches, exemplified by those that formed the Yzerfontein pluton, which was further hybridized through minor assimilation of crustal materials. Like Yzerfontein, the small volumes of mafic rocks associated with many batholiths, worldwide, are probably also lowvolume, high-level expressions of crustal growth through the emplacement of major amounts of mafic magma into the deep crust.IS