169 research outputs found
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
Probabilistic 3-D time-lapse inversion of magnetotelluric data: application to an enhanced geothermal system
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
- âŠ