S U M M A R YWe determine finite-frequency sensitivity kernels for seismic interferometry based upon noise cross-correlation measurements. Under the assumptions that noise is spatially uncorrelated but non-uniform, we determine ensemble-averaged cross correlations between synthetic seismograms at two geographically distinct locations. By minimizing a measure of the difference between observed and simulated ensemble cross correlations-subject to the constraint that the simulated wavefield satisfies the seismic wave equation-we obtain ensemble sensitivity kernels. These ensemble kernels reflect the sensitivity of ensemble cross-correlation measurements to variations in model parameters, for example, mass density, shear and compressional wave speeds and the spatial distribution of noise. Ensemble kernels are calculated based upon the interaction between two wavefields: an ensemble forward wavefield and an ensemble adjoint wavefield. To obtain the ensemble forward wavefield, one first calculates a generating wavefield obtained by inserting a signal determined by the characteristics of the noise at the location of the first receiver, saving the results of this calculation at locations where noise is generated, that is, typically on (a portion of) the Earth's surface. Next, one uses this generating wavefield as the source of the ensemble forward wavefield associated with the first receiver. The ensemble adjoint wavefield is obtained by using measurements between simulated and observed ensemble cross correlations as a seismic source located at the second receiver. The interaction between ensemble forward and adjoint wavefields 'paints' ensemble sensitivity kernels. We illustrate the construction of ensemble kernels and their nature in two and three dimensions using a spectral-element method. In addition to a 'banana-doughnut' feature connecting the two receivers, as in traditional finite-frequency earthquake tomography, some noise cross-correlation sensitivity kernels exhibit hyperbolic 'jets' protruding from each receiver in a direction away from the other receiver. Ensemble sensitivity kernels for long-period (T > ∼50 s) non-uniform noise in global models exhibit sensitivity along the minor and major arcs. These kernels reflect the fact that measurements typically involve long time-series that include multi-orbit surface waves. Like free oscillations, such measurements are sensitive to structure along the great circle through the two receivers. From the perspective of noise cross-correlation tomography, we discuss strategies for inversions in terrestrial and helioseismology.

International audienceWe present forward and adjoint spectral-element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite-frequency Fr'echet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite-frequency kernels generated by Version 2.0 'Sesame' of our widely used open source spectral-element package SPECFEM3D

No abstract

We have drawn connections between imaging in exploration seismology, adjoint methods, and emerging finite-frequency tomography. All of these techniques rely on spatial and temporal constructive interference between observed and simulated waveforms to map locations of structural anomalies. Modern numerical methods and computers have facilitated the accurate and efficient simulation of 3D acoustic, (an)elastic, and poroelastic wave propagation. Using a 2D cross section of the SEG/EAGE salt model, we have determined how such waveform simulations might be harnessed to improve onshore and offshore seismic imaging strategies and capabilities. We have found that the density sensitivity kernel in adjoint tomography is related closely to the imaging principle in exploration seismology, and that in elastic modeling the impedance kernel actually is a better diagnostic tool for reflector identification. The shear- and compressional-wave speed sensitivity kernels in adjoint tomography are related closely to finite-frequency banana-doughnut kernels, and these kernels are well suited for mapping larger-scale structure, i.e., for transmission tomography. These ideas have been substantiated by addressing problems in subsalt time-lapse migration.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ∼ 10 −6 pc and rapidly expands outwards.(2) A central core forms, with a mass of ∼ 1 M , supported by gas pressure gradients and rotation. (3) Growing gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ∼ 5 × 10 37 − 5 × 10 38 erg s −1 , of the order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ∼ 10 −4 − 10 −3 pc, and ∼ 0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, because the radiation is capable of escaping due to anisotropy in the optical depth and associated gradients.

The neutrino-cooled accretion disk, which was proposed to work as the central engine of gamma-ray bursts, encounters difficulty in interpreting the X-ray flares after the prompt gamma-ray emission. In this paper, the magnetic coupling between the inner disk and the central black hole is taken into consideration. For mass accretion rates around 0.001 ∼ 0.1 M ⊙ s −1 , our results show that the luminosity of neutrino annihilation can be significantly enhanced due to the coupling effects. As a consequence, after the gamma-ray emission, a remnant disk with mass M disk 0.5 M ⊙ may power most of the observed X-ray flares with the rest frame duration less than 100 seconds. In addition, a comparison between the magnetic coupling process and the Blandford-Znajek mechanism is shown on the extraction of black hole rotational energy.

No abstract

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ∼ 10 M and a highly variable luminosity of 10 38 − 10 39 erg s −1 , comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ∼ 10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10 −3 pc, which is ∼100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.