Automatic differentiation in geophysical inverse problems

Automatic differentiation (AD) is the technique whereby output variables of a computer code evaluating any complicated function (e.g. the solution to a differential equation) can be differentiated with respect to the input variables. Often AD tools take the form of source to source translators and produce computer code without the need for deriving and hand coding of explicit mathematical formulae by the user. The power of AD lies in the fact that it combines the generality of finite difference techniques and the accuracy and efficiency of analytical derivatives, while at the same time eliminating 'human' coding errors. It also provides the possibility of accurate, efficient derivative calculation from complex 'forward' codes where no analytical derivatives are possible and finite difference techniques are too cumbersome. AD is already having a major impact in areas such as optimization, meteorology and oceanography. Similarly it has considerable potential for use in non-linear inverse problems in geophysics where linearization is desirable, or for sensitivity analysis of large numerical simulation codes, for example, wave propagation and geodynamic modelling. At present, however, AD tools appear to be little used in the geosciences. Here we report on experiments using a state of the art AD tool to perform source to source code translation in a range of geoscience problems. These include calculating derivatives for Gibbs free energy minimization, seismic receiver function inversion, and seismic ray tracing. Issues of accuracy and efficiency are discussed

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/

Limits to reconstructing paleotopography from thermochronometer data

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95240/1/jgrf851.pd

The leading digit distribution of the worldwide Illicit Financial Flows

Benford's law states that in data sets from different phenomena leading digits tend to be distributed logarithmically such that the numbers beginning with smaller digits occur more often than those with larger ones. Particularly, the law is known to hold for different types of financial data. The Illicit Financial Flows (IFFs) exiting the developing countries are frequently discussed as hidden resources which could have been otherwise properly utilized for their development. We investigate here the distribution of the leading digits in the recent data on estimates of IFFs to look for the existence of a pattern as predicted by Benford's law and establish that the frequency of occurrence of the leading digits in these estimates does closely follow the law.Comment: 13 pages, 10 figures, 6 tables, additional data analyi

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

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

Statistics for Fission-Track Thermochronology

This chapter introduces statistical tools to extract geologically meaningful information from fission-track (FT) data using both the external detector and LA-ICP-MS methods. The spontaneous fission of 238U is a Poisson process resulting in large single-grain age uncertainties. To overcome this imprecision, it is nearly always necessary to analyse multiple grains per sample. The degree to which the analytical uncertainties can explain the observed scatter of the single-grain data can be visually assessed on a radial plot and objectively quantified by a chi-square test. For sufficiently low values of the chi-square statistic (or sufficiently high p values), the pooled age of all the grains gives a suitable description of the underlying ‘true’ age population. Samples may fail the chi-square test for several reasons. A first possibility is that the true age population does not consist of a single discrete age component, but is characterised by a continuous range of ages. In this case, a ‘random effects’ model can constrain the true age distribution using two parameters: the ‘central age’ and the ‘(over)dispersion’. A second reason why FT data sets might fail the chi-square test is if they are underlain by multimodal age distributions. Such distributions may consist of discrete age components, continuous age distributions, or a combination of the two. Formalised statistical tests such as chi-square can be useful in preventing overfitting of relatively small data sets. However, they should be used with caution when applied to large data sets (including length measurements) which generate sufficient statistical ‘power’ to reject any simple yet geologically plausible hypothesis

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

Thermochronology of mineral grains in the Red and Mekong Rivers, Vietnam: Provenance and exhumation implications for Southeast Asia

Sand samples from the mouths of the Red and Mekong Rivers were analyzed to determine the provenance and exhumation history of their source regions. U-Pb dating of detrital zircon grains shows that the main sources comprise crust formed within the Yangtze Craton and during the Triassic Indosinian Orogeny. Indosinian grains in the Mekong are younger (210-240 Ma) than those in the Red River (230-290 Ma), suggesting preferential erosion of the Qiangtang Block of Tibet into the Mekong. The Red River has a higher proportion of 700-800 Ma grains originally derived from the Yangtze Craton. 40Ar/ 39Ar dating of muscovite grains demonstrates that rocks cooled during the Indosinian Orogeny are dominant in both rivers, although the Mekong also shows a grain population cooling at 150-200 Ma that is not seen in the Red River and which is probably of original Qiangtang Block origin. Conversely, the Red River contains a significant mica population (350-500 Ma) eroded from the Yangtze Craton. High-grade metamorphic rocks exposed in the Cenozoic shear zones of southeast Tibet-Yunnan are minority sources to the rivers. However, apatite and zircon fission track ages show evidence for the dominant sources, especially in the Red River, only being exhumed through the shallowest 5-3 km of the crust since ̃25 Ma. The thermochronology data are consistent with erosion of recycled sediment from the inverted Simao and Chuxiong Basins, from gorges that incise the eastern flank of the plateau. Average Neogene exhumation rates are 104-191 m/Myr in the Red River basin, which is within error of the 178 ± 35 m/Myr estimated from Pleistocene sediment volumes. Sparse fission track data from the Mekong River support the Ar-Ar and U-Pb ages in favoring tectonically driven rock uplift and gorge incision as the dominant control on erosion, with precipitation being an important secondary influence. © 2006 by the American Geophysical Union