Time-domain sparsity promoting least-squares reverse time migration with source estimation

Least-squares reverse-time migration is well known for its capability to generate artefact-free true-amplitude subsurface images through fitting observed data in the least-squares sense. However, when applied to realistic imaging problems, this approach is faced with issues related to overfitting and excessive computational costs induced by many wave-equation solves. The fact that the source function is unknown complicates this situation even further. Motivated by recent results in stochastic optimization and transform-domain sparsity promotion, we demonstrate that the computational costs of inversion can be reduced significantly while avoiding imaging artefacts and restoring amplitudes. While powerful, these new approaches do require accurate information on the source-time function, which is often lacking. Without this information, the imaging quality deteriorates rapidly. We address this issue by presenting an approach where the source-time function is estimated on the fly through a technique known as variable projection. Aside from introducing negligible computational overhead, the proposed method is shown to perform well on imaging problems with noisy data and problems that involve complex settings such as salt. In either case, the presented method produces high-resolution high-amplitude fidelity images including an estimate for the source-time function. In addition, due to its use of stochastic optimization, we arrive at these images at roughly one to two times the cost of conventional reverse-time migration involving all data.     

Slowness and azimuth determination for Bucovina array (BURAR) applying multiple signal techniques

The BURAR seismic array, located in Northern Romania (Bucovina region), is designed to monitor events located in an area poorly covered by other existing seismic stations. In order to use the BURAR array for single-station locations, it is crucial to calibrate the azimuth and slowness parameters, which are currently used in array techniques to locate earthquakes, blasts or nuclear explosions. The goal of this study is to apply “f–k” and plane wave fit techniques in order to constrain the slowness and azimuth parameters at BURAR for teleseismic, regional and local events. The analysis was carried out using P and S waves recorded for events occurred between 2004 and 2008 within a radius of 50° around BURAR. The azimuth values obtained applying both methods strongly deviated from the theoretical values for regions like Central Turkey, Bulgaria, Dodecanese Islands and other parts of Greece, while the ray parameter deviations with respect to a 1-D IASP91 reference model are less significant. For the local events, the anomalies are smaller, except the particular case of Vrancea intermediate-depth earthquakes for which strong azimuth deviations (33.5°), both positive and negative, are observed. We investigate how these systematic deviations in azimuth are explained by the structure lateral heterogeneities which characterize the study region.     

Gamma ray signatures of ultra high energy cosmic ray accelerators: electromagnetic cascade versus synchrotron radiation of secondary electrons

We discuss the possibility of observing ultra high energy cosmic ray sources in high energy gamma rays. Protons propagating away from their accelerators produce secondary electrons during interactions with cosmic microwave background photons. These electrons start an electromagnetic cascade that results in a broad band gamma ray emission. We show that in a magnetized Universe (B≳10⁻¹² G) such emission is likely to be too extended to be detected above the diffuse background. A more promising possibility comes from the detection of synchrotron photons from the extremely energetic secondary electrons. Although this emission is produced in a rather extended region of size ∼10 Mpc, it is expected to be point-like and detectable at GeV energies if the intergalactic magnetic field is at the nanogauss level.      

Size and shape characteristics of drumlins,derived from a large sample,and associated scaling laws

Ice sheets flowing across a sedimentary bed usually produce a landscape of blister-like landforms streamlined in the direction of the ice flow and with each bump of the order of 10² to 10³ m in length and 10¹ m in relief. Such landforms, known as drumlins, have mystified investigators for over a hundred years. A satisfactory explanation for their formation, and thus an appreciation of their glaciological significance, has remained elusive. A recent advance has been in numerical modelling of the land-forming process. In anticipation of future modelling endeavours, this paper is motivated by the requirement for robust data on drumlin size and shape for model testing.From a systematic programme of drumlin mapping from digital elevation models and satellite images of Britain and Ireland, we used a geographic information system to compile a range of statistics on length L, width W, and elongation ratio E (where E = L/W) for a large sample. Mean L, is found to be 629 m (n = 58,983), mean W is 209 m and mean E is 2.9 (n = 37,043). Most drumlins are between 250 and 1000 metres in length; between 120 and 300 metres in width; and between 1.7 and 4.1 times as long as they are wide. Analysis of such data and plots of drumlin width against length reveals some new insights. All frequency distributions are unimodal from which we infer that the geomorphological label of ‘drumlin’ is fair in that this is a true single population of landforms, rather than an amalgam of different landform types. Drumlin size shows a clear minimum bound of around 100 m (horizontal). Maybe drumlins are generated at many scales and this is the minimum, or this value may be an indication of the fundamental scale of bump generation (‘proto-drumlins’) prior to them growing and elongating. A relationship between drumlin width and length is found (with r² = 0.48) and that is approximately W = 7 L ^1/2 when measured in metres. A surprising and sharply-defined line bounds the data cloud plotted in E–W–L space, and records a scale-dependent maximum elongation limit (approximated by E_max = L^1/3, when L measured in metres). For a given length, for some reason as yet unknown, drumlins do not exceed the elongation ratio defined by this scaling law. We also report and compare our statistics to an amalgamated sample (25,907 drumlins) of measures derived from around 50 published investigations. Any theory must be able to explain the drumlin statistics and fundamental scaling properties reported herein and they thus provide powerful tests for drumlin modelling.     

The future of artificial satellite theories Hybrid ephemeris compression model

Since the time of Newton, astrodynamics has focused on the analytical solution of orbital problems. This was driven by the desire to obtain a theoretical understanding of the motion and the practical desire to be able to produce a computational result, Only with the advent of the computer did numerical integration become a practical consideration for solving dynamical problems. Although computer technology is not yet to the point of being able to provide numerical integration support for all satellite orbits, we are in a transition period which is being driven by the unprecedented increase in computational power, This transition will affect the future of analytical, semi-analytical and numerical artificial satellite theories in a dramatic way, In fact, the role for semi-analytical theories may disappear. During the time of transition, a central site may have the capacity to maintain the orbits using numerical integration, but the user may not have such a capacity or may need results in a more timely manner, One way to provide for this transition need is through the use of some type of satellite ephemeris compression. Through the combined use of a power series and a Fourier series, good quality ephemeris compression has been achieved for 7 day periods, The ephemeris compression requires less than 40 terms and is valid for all eccentricities and inclinations.     

Modulation of the high energy γ-ray flux from PSR B1259-63/SS2883 due to the orbital variation of the maximum energy of accelerated electrons

The inverse Compton (IC) scattering of ultrarelativistic electrons accelerated at the pulsar wind termination shock is generally believed to be responsible for TeV gamma-ray signal recently reported from the binary system PSR B1259-63/SS2883. In such a system the acceleration takes place in the presence of a dense radiation field provided by a companion Be2-type star. Thus it is natural to expect an orbital phase dependence of the acceleration efficiency in the system. The HESS collaboration reported the tendency of reduction of TeV γ-rays around the periastron. In this paper we study a possible explanation of this effect by the “early” (sub-TeV) cutoffs in the energy spectrum of accelerated electrons due to the enhanced rate of Compton losses close to the periastron.     

Interplay between scintillation and ionization in liquid xenon Dark Matter searches

We provide a new way of constraining the relative scintillation efficiency L_eff for liquid xenon. Using a simple estimate for the electronic and nuclear stopping powers together with an analysis of recombination processes we predict both the ionization and the scintillation yields. Using presently available data for the ionization yield, we can use the correlation between these two quantities to constrain L_eff from below. Moreover, we argue that more reliable data on the ionization yield would allow to verify our assumptions on the atomic cross sections and to predict the value of L_eff. We conclude that the relative scintillation efficiency should not decrease at low nuclear recoil energies, which has important consequences for the robustness of exclusion limits for low WIMP masses in liquid xenon Dark Matter searches.