Scientific theories of how subduction and plate tectonics began on Earth--and what the tectonic structure of Earth was before this--remain enigmatic and contentious. Understanding viable scenarios ...for the onset of subduction and plate tectonics is hampered by the fact that subduction initiation processes must have been markedly different before the onset of global plate tectonics because most present-day subduction initiation mechanisms require acting plate forces and existing zones of lithospheric weakness, which are both consequences of plate tectonics. However, plume-induced subduction initiation could have started the first subduction zone without the help of plate tectonics. Here, we test this mechanism using high-resolution three-dimensional numerical thermomechanical modelling. We demonstrate that three key physical factors combine to trigger self-sustained subduction: (1) a strong, negatively buoyant oceanic lithosphere; (2) focused magmatic weakening and thinning of lithosphere above the plume; and (3) lubrication of the slab interface by hydrated crust. We also show that plume-induced subduction could only have been feasible in the hotter early Earth for old oceanic plates. In contrast, younger plates favoured episodic lithospheric drips rather than self-sustained subduction and global plate tectonics.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, KISLJ, NUK, PILJ, PNG, SAZU, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
Aims. The DustPedia project is capitalising on the legacy of the Herschel Space Observatory, using cutting-edge modelling techniques to study dust in the 875 DustPedia galaxies – representing the ...vast majority of extended galaxies within 3000 km s-1 that were observed by Herschel. This work requires a database of multiwavelength imagery and photometry that greatly exceeds the scope (in terms of wavelength coverage and number of galaxies) of any previous local-Universe survey. Methods. We constructed a database containing our own custom Herschel reductions, along with standardised archival observations from GALEX, SDSS, DSS, 2MASS, WISE, Spitzer, and Planck. Using these data, we performed consistent aperture-matched photometry, which we combined with external supplementary photometry from IRAS and Planck. Results. We present our multiwavelength imagery and photometry across 42 UV-microwave bands for the 875 DustPedia galaxies. Our aperture-matched photometry, combined with the external supplementary photometry, represents a total of 21 857 photometric measurements. A typical DustPedia galaxy has multiwavelength photometry spanning 25 bands. We also present the Comprehensive & Adaptable Aperture Photometry Routine (CAAPR), the pipeline we developed to carry out our aperture-matched photometry. CAAPR is designed to produce consistent photometry for the enormous range of galaxy and observation types in our data. In particular, CAAPR is able to determine robust cross-compatible uncertainties, thanks to a novel method for reliably extrapolating the aperture noise for observations that cover a very limited amount of background. Our rich database of imagery and photometry is being made available to the community.
Aims. The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the ...DustPedia sample. Methods. This study is performed by exploiting the multi-wavelength DustPedia database, from ultraviolet (UV) to sub-millimeter bands, in addition to molecular (12CO) and atomic (Hi) gas maps and metallicity abundance information available in the literature. We fitted the surface-brightness profiles of the tracers of dust and stars, the mass surface-density profiles of dust, stars, molecular gas, and total gas, and the SFR surface-density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (αCO) per galaxy by using dust- and gas-mass profiles. Results. Although each galaxy has its own peculiar behavior, we identified a common trend of the exponential scale-lengths versus wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec2 radius decrease from UV to 70 μm, from 0.4 to 0.2, and then increase back up to ~0.3 at 500 microns. The main result is that, on average, the dust-mass surface-density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 μm surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 μm scale-length tend to increase from earlier to later types, the scale-length at 70 μm tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 μm that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types. Our αCO determinations are in the range (0.3−9) M⊙ pc-2 (K km s-1)-1, almost invariant by using a fixed dust-to-gas ratio mass (DGR) or a DGR depending on metallicity gradient.
In high density environments, the gas content of galaxies is stripped, leading to a rapid quenching of their star formation activity. This dramatic environmental effect, which is not related to ...typical passive evolution, is generally not taken into account in the star formation histories (SFHs) usually assumed to perform spectral energy distribution (SED) fitting of these galaxies, yielding a poor fit of their stellar emission and, consequently, biased estimate of the star formation rate (SFR). In this work, we aim at reproducing this rapid quenching using a truncated delayed SFH that we implemented in the SED fitting code CIGALE. We show that the ratio between the instantaneous SFR and the SFR just before the quenching (rSFR) is well constrained as long as rest-frame UV data are available. This SED modeling is applied to the Herschel Reference Survey (HRS) containing isolated galaxies and sources falling in the dense environment of the Virgo cluster. The latter are Hi-deficient because of ram pressure stripping. We show that the truncated delayed SFH successfully reproduces their SED, while typical SFH assumptions fail. A good correlation is found between rSFR and Hi−def, the parameter that quantifies the gas deficiency of cluster galaxies, meaning that SED fitting results can be used to provide a tentative estimate of the gas deficiency of galaxies for which Hi observations are not available. The HRS galaxies are placed on the SFR-M∗ diagram showing that the Hi-deficient sources lie in the quiescent region, thus confirming previous studies. Using the rSFR parameter, we derive the SFR of these sources before quenching and show that they were previously on the main sequence relation. We show that the rSFR parameter is also recovered well for deeply obscured high redshift sources, as well as in the absence of IR data. SED fitting is thus a powerful tool for identifying galaxies that underwent a rapid star formation quenching.
Aims.
The purpose of this work is the characterization of the main scaling relations between all of the interstellar medium (ISM) components, namely dust, atomic, molecular, and total gas, and ...gas-phase metallicity, as well as other galaxy properties, such as stellar mass (
M
star
) and galaxy morphology, for late-type galaxies in the Local Universe.
Methods.
This study was performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage spanning from
T
= 1−10,
M
star
from 6 × 10
7
to 3 × 10
11
M
⊙
, star formation rate from 6 × 10
−4
to 60
M
⊙
yr
−1
, and oxygen abundance from 12 + log(O/H) = 8−9.5. Molecular and atomic gas data were collected from the literature and properly homogenized. All the masses involved in our analysis refer to the values within the optical disks of galaxies. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent CO-to-H
2
conversion factor (
X
CO
). The analysis was performed by means of the survival analysis technique, in order to properly take into account the presence of both detection and nondetection in the data.
Results.
We confirm that the dust mass correlates very well with the total gas mass, and find –for the first time– that the dust mass correlates better with the atomic gas mass than with the molecular one. We characterize important mass ratios such as the gas fraction, the molecular-to-atomic gas mass ratio, the dust-to-total gas mass ratio (DGR), and the dust-to-stellar mass ratio, and study how they relate to each other, to galaxy morphology, and to gas-phase metallicity. Only the assumption of a metallicity-dependent
X
CO
reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. The DGR, the gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for galaxies rich in molecular gas this trend breaks down. To our knowledge, this trend has never been found before, and provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies. We discuss several scenarios related to this finding.
Conclusions.
The DustPedia database of late-type galaxies is an extraordinary tool for the study of the ISM scaling relations, thanks to its homogeneous collection of data for the different ISM components. The database is made publicly available to the whole community.
Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a ...sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10 000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction >60%, metallicity 12 + log(O∕H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction >80%) than the typical dust-to-metal ratio (Md∕MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass.
ABSTRACT
There is a severe tension between the observed star formation rate (SFR)–stellar mass (M⋆) relations reported by different authors at z = 1–4. In addition, the observations have not been ...successfully reproduced by state-of-the-art cosmological simulations that tend to predict a factor of 2–4 smaller SFRs at a fixed M⋆. We examine the evolution of the SFR–M⋆ relation of z = 1–4 galaxies using the skirt simulated spectral energy distributions of galaxies sampled from the Evolution and Assembly of GaLaxies and their Environments simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR–M⋆ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \, \mu m}$, MIPS – 24, 70, and 160 ${\rm \, \mu m}$ or SPIRE – 250, 350, and 500 ${\rm \, \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by −0.5 dex at z ∼ 4 but overpredict the measurements by 0.3 dex at z ∼ 1. Relations relying on dust-corrected rest-frame ultraviolet luminosities, are flatter since they overpredict/underpredict SFRs for low/high star-forming objects and yield deviations from the intrinsic relation from 0.10 to −0.13 dex at z ∼ 4. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.
Context. A crucial ingredient for numerically solving the three-dimensional radiative transfer problem is the choice of the grid that discretizes the transfer medium. Many modern radiative transfer ...codes, whether using Monte Carlo or ray tracing techniques, are equipped with hierarchical octree-based grids to accommodate a wide dynamic range in densities. Aims. We critically investigate two different aspects of octree grids in the framework of Monte Carlo dust radiative transfer. Inspired by their common use in computer graphics applications, we test hierarchical k-d tree grids as an alternative for octree grids. On the other hand, we investigate which node subdivision-stopping criteria are optimal for constructing of hierarchical grids. Methods. We implemented a k-d tree grid in the 3D radiative transfer code SKIRT and compared it with the previously implemented octree grid. We also considered three different node subdivision-stopping criteria (based on mass, optical depth, and density gradient thresholds). Based on a small suite of test models, we compared the efficiency and accuracy of the different grids, according to various quality metrics. Results. For a given set of requirements, the k-d tree grids only require half the number of cells of the corresponding octree. Moreover, for the same number of grid cells, the k-d tree is characterized by higher discretization accuracy. Concerning the subdivision stopping criteria, we find that an optical depth criterion is not a useful alternative to the more standard mass threshold, since the resulting grids show a poor accuracy. Both criteria can be combined; however, in the optimal combination, for which we provide a simple approximate recipe, this can lead to a 20% reduction in the number of cells needed to reach a certain grid quality. An additional density gradient threshold criterion can be added that solves the problem of poorly resolving sharp edges and strong density gradients. Conclusions. We advocate the use of k-d trees and the proposed combination of criteria to set up hierarchical grids for 3D radiative transfer. These recipes are straightforward for implementing and should help to develop faster and more accurate 3D radiative transfer codes.
Here, we study the 3-D subduction initiation process induced by the interaction between a hot thermo-chemical mantle plume and oceanic lithosphere using thermo-mechanical viscoplastic finite ...difference marker-in-cell models. Our numerical modeling results show that self-sustaining subduction is induced by plume–lithosphere interaction when the plume is sufficiently buoyant, the oceanic lithosphere is sufficiently old and the plate is weak enough to allow the buoyant plume to pass through it. Subduction initiation occurs following penetration of the lithosphere by the hot plume and the downward displacement of broken, nearly circular segments of lithosphere (proto-slabs) as a result of partially molten plume rocks overriding the proto-slabs. Our experiments show four different deformation regimes in response to plume–lithosphere interaction: a) self-sustaining subduction initiation, in which subduction becomes self-sustaining; b) frozen subduction initiation, in which subduction stops at shallow depths; c) slab break-off, in which the subducting circular slab breaks off soon after formation; and d) plume underplating, in which the plume does not pass through the lithosphere and instead spreads beneath it (i.e., failed subduction initiation). These regimes depend on several parameters, such as the size, composition, and temperature of the plume, the brittle/plastic strength and age of the oceanic lithosphere, and the presence/absence of lithospheric heterogeneities. The results show that subduction initiates and becomes self-sustaining when the lithosphere is older than 10 Myr and the non-dimensional ratio of the plume buoyancy force and lithospheric strength above the plume is higher than approximately 2. The outcomes of our numerical experiments are applicable for subduction initiation in the modern and Precambrian Earth and for the origin of plume-related corona structures on Venus.
•Four different responses result from plume–lithosphere interaction.•Responses depend on plume's buoyancy, age and strength of the oceanic lithosphere.•Self-sustained subduction initiates if overlaying oceanic lithosphere is >10 Myr.•Results can apply for subduction initiation in modern and Precambrian Earth.
Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering ...of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy.