Constraining dark matter properties with SPI

Using the high-resolution spectrometer SPI on board the International Gamma-Ray Astrophysics Laboratory ( INTEGRAL ), we search for a spectral line produced by a dark matter (DM) particle with a mass in the range  40 keV < M _DM < 14 MeV  , decaying in the DM halo of the Milky Way. To distinguish the DM decay line from numerous instrumental lines found in the SPI background spectrum, we study the dependence of the intensity of the line signal on the offset of the SPI pointing from the direction toward the Galactic Centre. After a critical analysis of the uncertainties of the DM density profile in the inner Galaxy, we find that the intensity of the DM decay line should decrease by at least a factor of 3 when the offset from the Galactic Centre increases from 0° to 180°. We find that such a pronounced variation of the line flux across the sky is not observed for any line, detected with a significance higher than 3σ in the SPI background spectrum. Possible DM decay origin is not ruled out only for the unidentified spectral lines, having low (∼3σ) significance or coinciding in position with the instrumental ones. In the energy interval from 20 keV to 7 MeV, we derive restrictions on the DM decay line flux, implied by the (non-)detection of the DM decay line. For a particular DM candidate, the sterile neutrino of mass M _DM, we derive a bound on the mixing angle.     

Constraints on decaying dark matter from XMM–Newton observations of M31

We derive constraints on the parameters of the radiatively decaying dark matter (DM) particle, using the XMM–Newton EPIC spectra of the Andromeda galaxy (M31). Using the observations of the outer (5–13 arcmin) parts of M31, we improve the existing constraints. For the case of sterile neutrino DM, combining our constraints with the latest computation of abundances of sterile neutrinos in the Dodelson–Widrow (DW) scenario, we obtain the lower mass limit   m _s < 4 keV  , which is stronger than the previous one   m _s < 6 keV  , obtained recently by Asaka, Laine & Shaposhnikov. Comparing this limit with the most recent results on Lyman α forest analysis of Viel et al.  ( m _s > 5.6 keV  ), we argue that the scenario in which all the DM is produced via the DW mechanism is ruled out. We discuss, however, other production mechanisms and note that the sterile neutrino remains a viable candidate for DM, either warm or cold.     

Seismic Source Characteristics of Nuclear and Chemical Explosions in Granite from Hydrodynamic Simulations

Seismic source characteristics of low-yield (0.5–5 kt) underground explosions are inferred from hydrodynamic simulations using a granite material model on high-performance (parallel) computers. We use a non-linear rheological model for granite calibrated to historical near-field nuclear test data. Equivalent elastic P-wave source spectra are derived from the simulated hydrodynamic response using reduced velocity potentials. Source spectra and parameters are compared with the models of Mueller and Murphy (Bull Seism Soc Am 61:1675–1692, 1971, hereafter MM71) and Denny and Johnson (Explosion source phenomenology, pp 1–24, 1991, hereafter DJ91). The source spectra inferred from the simulations of different yields at normal scaled depth-of-burial (SDOB) match the MM71 spectra reasonably well. For normally buried nuclear explosions, seismic moments are larger for the hydrodynamic simulations than MM71 (by 25 %) and for DJ91 (by over a factor of 2), however, the scaling of moment with yield across this low-yield range is consistent for our calculations and the two models. Spectra from our simulations show higher corner frequencies at the lower end of the 0.5–5.0 kt yield range and stronger variation with yield than the MM71 and DJ91 models predict. The spectra from our simulations have additional energy above the corner frequency, probably related to non-linear near-source effects, but at high frequencies the spectral slopes agree with the f ^?2 predictions of MM71. Simulations of nuclear explosions for a range of SDOB from 0.5 to 3.9 show stronger variations in the seismic moment than predicted by the MM71 and DJ91 models. Chemical explosions are found to generate higher moments by a factor of about two compared to nuclear explosions of the same yield in granite and at normal depth-of-burial, broadly consistent with comparisons of nuclear and chemical shots at the US Nevada Test Site (Denny, Proceeding of symposium on the non-proliferation experiment, Rockville, Maryland, 1994). For all buried explosions, the region of permanent deformation and material damage is not spherical but extends along the free surface above and away from the source. The effect of damage induced by a normally buried nuclear explosion on seismic radiation is explored by comparing the motions from hydrodynamic simulations with those for point-source elastic Green’s functions. Results show that radiation emerging at downward takeoff angles appears to be dominated by the expected isotropic source contribution, while at shallower angles the motions are complicated by near-surface damage and cannot be represented with the addition of a simple secondary compensated linear vector dipole point source above the shot point. The agreement and differences of simulated source spectra with the MM71 and DJ91 models motivates the use of numerical simulations to understand observed motions and investigate seismic source features for underground explosions in various emplacement media and conditions, including non-linear rheological effects such as material strength and porosity.     

Atmospheric boundary layer over steep surface waves

Turbulent air-sea interactions coupled with the surface wave dynamics remain a challenging problem. The needs to include this kind of interaction into the coupled environmental, weather and climate models motivate the development of a simplified approximation of the complex and strongly nonlinear interaction processes. This study proposes a quasi-linear model of wind-wave coupling. It formulates the approach and derives the model equations. The model is verified through a set of laboratory (direct measurements of an airflow by the particle image velocimetry (PIV) technique) and numerical (a direct numerical simulation (DNS) technique) experiments. The experiments support the central model assumption that the flow velocity field averaged over an ensemble of turbulent fluctuations is smooth and does not demonstrate flow separation from the crests of the waves. The proposed quasi-linear model correctly recovers the measured characteristics of the turbulent boundary layer over the waved water surface.     

Applied Mathematics in EM Studies with Special Emphasis on an Uncertainty Quantification and 3-D Integral Equation Modelling

Despite impressive progress in the development and application of electromagnetic (EM) deterministic inverse schemes to map the 3-D distribution of electrical conductivity within the Earth, there is one question which remains poorly addressed—uncertainty quantification of the recovered conductivity models. Apparently, only an inversion based on a statistical approach provides a systematic framework to quantify such uncertainties. The Metropolis–Hastings (M–H) algorithm is the most popular technique for sampling the posterior probability distribution that describes the solution of the statistical inverse problem. However, all statistical inverse schemes require an enormous amount of forward simulations and thus appear to be extremely demanding computationally, if not prohibitive, if a 3-D set up is invoked. This urges development of fast and scalable 3-D modelling codes which can run large-scale 3-D models of practical interest for fractions of a second on high-performance multi-core platforms. But, even with these codes, the challenge for M–H methods is to construct proposal functions that simultaneously provide a good approximation of the target density function while being inexpensive to be sampled. In this paper we address both of these issues. First we introduce a variant of the M–H method which uses information about the local gradient and Hessian of the penalty function. This, in particular, allows us to exploit adjoint-based machinery that has been instrumental for the fast solution of deterministic inverse problems. We explain why this modification of M–H significantly accelerates sampling of the posterior probability distribution. In addition we show how Hessian handling (inverse, square root) can be made practicable by a low-rank approximation using the Lanczos algorithm. Ultimately we discuss uncertainty analysis based on stochastic inversion results. In addition, we demonstrate how this analysis can be performed within a deterministic approach. In the second part, we summarize modern trends in the development of efficient 3-D EM forward modelling schemes with special emphasis on recent advances in the integral equation approach.     

Karibibite,a new FeAs mineral from South West Africa

Karibibite (ideally, Fe₂As₄O₉) occurs in vugs in massive loellingite of the Karibib pegmatite area, South West Africa. It is brownish yellow and finely fibrous. The thickness of the soft, single fibers is less than 1 micron, unsuitable for single-crystal X-ray study. Electron diffraction and X-ray powder pattern indicate that the mineral is orthorhombic, with $\text{a}{}_0\text{= 27.91}\text{A}\text{?}\text{, b}{}_0\text{= 6.53}\text{A}\text{?}\text{and c}{}_0\text{(}\text{fiber axis}\text{) = 7.20}\text{A}\text{?}$. The space group cannot be given. The mineral is paramagnetic with yellow fluorescence and is pleochroic with γ > 2.10, α = 1.96, 2V_α large, d = 4.07. It is soluble in acids and alkali hydroxide. Decomposition starts around 320 °C. The infra-red absorption spectrum indicates absence of AsO₄ groups. The mineral is classified tentatively as an oxide or arsenite.     

Evaluation of formation pore pressure behind the casing using borehole gravity data

Reliable estimates of the fluid pressure in the pore space of rocks are critical for different aspects of petroleum exploration and production including injection operations and scenarios of water flooding. Numerous approaches are available for formation pore pressure evaluation, however, these measurements become a challenge inside a cased borehole, and a list of possible options is short: either the casing is to be perforated, or the production tubing needs to be disconnected to perform the pressure tests. We present a method for through-casing evaluation of formation pore pressure without shutting down production. We suggest equipping an observation well with a borehole gravimeter and acquiring time variations of the vertical component of the gravity field. Changes in gravity occur during gas production and are related to time variations of formation pore pressure. Gravity changes obtained in the observation well are supposed to be inverted for time-dependent formation pore pressure variations beyond the casing. Our results and recommendations are based on numerical modeling of pore pressure spatial distribution during gas field exploitation and relevant changes in borehole gravity. Benchmark models were elaborated in order to consider a dynamic process of pressure changes in time and space under conditions similar to those in the Medvezhye gas field (Russia). Different modeling scenarios are considered for early and late stages of gas field exploitation. The sensitivity analysis was performed to estimate quantitatively a sensitivity of borehole temporal gravity changes to variations in formation pore pressure behind the casing. Based on resolution analysis we justify the possibility to extract the gravity measurements directly related to changes in pore pressure from the total changes in the gravity field due to reservoir exploitation. The impact of pore pressure on the gravity field measured in boreholes during the water flooding is also evaluated, and obtained results are discussed.