ABSTRACT
This work presents AutoLens, the first entirely automated modeling suite for the analysis of galaxy-scale strong gravitational lenses. AutoLens simultaneously models the lens galaxy’s light ...and mass whilst reconstructing the extended source galaxy on an adaptive pixel-grid. The method’s approach to source-plane discretization is amorphous, adapting its clustering and regularization to the intrinsic properties of the lensed source. The lens’s light is fitted using a superposition of Sersic functions, allowing AutoLens to cleanly deblend its light from the source. Single-component mass models representing the lens’s total mass density profile are demonstrated, which in conjunction with light modeling can detect central images using a centrally cored profile. Decomposed mass modeling is also shown, which can fully decouple a lens’s light and dark matter and determine whether the two components are geometrically aligned. The complexity of the light and mass models is automatically chosen via Bayesian model comparison. These steps form AutoLens’s automated analysis pipeline, such that all results in this work are generated without any user intervention. This is rigorously tested on a large suite of simulated images, assessing its performance on a broad range of lens profiles, source morphologies, and lensing geometries. The method’s performance is excellent, with accurate light, mass, and source profiles inferred for data sets representative of both existing Hubble imaging and future Euclid wide-field observations.
Chemical and physical Earth models agree little as to the radioactive power of the planet. Each predicts a range of radioactive powers, overlapping slightly with the other at about 24 TW, and ...together spanning 14–46 TW. Approximately 20% of this radioactive power (3–8 TW) escapes to space in the form of geoneutrinos. The remaining 11–38 TW heats the planet with significant geodynamical consequences, appearing as the radiogenic component of the 43–49 TW surface heat flow. The nonradiogenic component of the surface heat flow (5–38 TW) is presumably primordial, a legacy of the formation and early evolution of the planet. A constraining measurement of radiogenic heating provides insights to the thermal history of the Earth and potentially discriminates chemical and physical Earth models. Radiogenic heating in the planet primarily springs from unstable nuclides of uranium, thorium, and potassium. The paths to their stable daughter nuclides include nuclear beta decays, producing geoneutrinos. Large subsurface detectors efficiently record the energy but not the direction of the infrequent interactions of the highest‐energy geoneutrinos, originating only from uranium and thorium. The measured energy spectrum of the interactions estimates the relative amounts of these heat‐producing elements, while the intensity estimates planetary radiogenic power. Recent geoneutrino observations in Japan and Italy find consistent values of radiogenic heating. The combined result mildly excludes the lowest model values of radiogenic heating and, assuming whole mantle convection, identifies primordial heat loss. Future observations have the potential to measure radiogenic heating with better precision, further constraining geological models and the thermal evolution of the Earth. This review presents the science and status of geoneutrino observations and the prospects for measuring the radioactive power of the planet.
Key Points
Radiogenic heating and geoneutrino flux are closely related
Earth models predict measurably distinct radioactive heat and geoneutrino flux
Oceanic geoneutrino measurement allows resolution of Earth radioactive power
We present a new pixelized method for the inversion of gravitationally lensed extended source images which we term adaptive semi-linear inversion (SLI). At the heart of the method is an h-means ...clustering algorithm which is used to derive a source plane pixelization that adapts to the lens model magnification. The distinguishing feature of adaptive SLI is that every pixelization is derived from a random initialization, ensuring that data discretization is performed in a completely different and unique way for every lens model parameter set. We compare standard SLI on a fixed source pixel grid with the new method and demonstrate the shortcomings of the former when modelling singular power-law ellipsoid (SPLE) lens profiles. In particular, we demonstrate the superior reliability and efficiency of adaptive SLI which, by design, fixes the number of degrees of freedom (NDOF) of the optimization and thereby removes biases present with other methods that allow the NDOF to vary. In addition, we highlight the importance of data discretization in pixel-based inversion methods, showing that adaptive SLI averages over significant systematics that are present when a fixed source pixel grid is used. In the case of the SPLE lens profile, we show how the method successfully samples its highly degenerate posterior probability distribution function with a single non-linear search. The robustness of adaptive SLI provides a firm foundation for the development of a strong lens modelling pipeline, which will become necessary in the short-term future to cope with the increasing rate of discovery of new strong lens systems.
Synchronous acceleration and thinning of southeast (SE) Greenland glaciers during the early 2000s was the main contributor that resulted in the doubling of annual discharge from the ice sheet. We ...show that this acceleration was followed by a synchronized and widespread slowdown of the same glaciers, in many cases associated with a decrease in thinning rates, and we propose that ice sheet–ocean interactions are the first‐order regional control on these recent mass changes. Sea surface temperature and mooring data show that the preceding dynamic thinning coincides with a brief decline in the cold East Greenland Coastal Current (EGCC) and East Greenland Current. We suggest this decline was partly induced by a reduction in ice sheet runoff, which allowed warm water from the Irminger Current to reach the SE Greenland coast. A restrengthening of the cold waters coincides with the glaciers' subsequent slowdown. We argue that this warming and subsequent cooling of the coastal waters was the cause of the glaciers' dynamic changes. We further suggest that the restrengthening of the EGCC resulted in part from cold water input by increased glacier calving during the speedup and increased ice sheet runoff. We hypothesize that the main mechanism for ice sheet mass loss in SE Greenland is highly sensitive to ocean conditions and is likely subject to negative feedback mechanisms.
SUMMARY
Many recent studies have tried to determine the influence of geometry of faults in earthquake mechanics. However, it still remains largely unknown, and it is explored mainly with numerical ...models. In this paper, we will try to understand how exactly does the geometry come into play in the mechanics of an earthquake from analytical perspective. We suggest a new interpretation of the effect of geometry on the stress on a fault, where the curvatures of the fault that multiply the slip play a major role. Starting from the representation theorem, which links the displacement in a medium to the slip distribution on its boundary, and assuming a 3-D, homogeneous, infinite medium, a regularized boundary-element equation can be obtained. Using this equation, it is possible to separate the influence of geometry, as expressed by the curvatures and torsions of the field lines of slip on the fault surface, which multiply the slip, from the effect of the gradient of slip. This allows us to shed new light on the mechanical effects of geometrical complexities on the fault surface, with the key parameters being the curvatures and torsions of the slip field lines. We have used this new approach to show that, in 2-D static in-plane problems, the shear traction along the fault is mainly controlled by the gradient of slip along the fault, while the normal traction is mainly controlled by the slip that multiplies the curvature along the fault. Finally, we applied this new approach to re-interpret the effect of roughness (why there is a need for a minimum lengthscale in linear elasticity, how to study mechanically the difference of roughness measurements with the direction of slip, scaling of slip distribution versus geometry), bends and kinks (what is the difference between the two, are they sometimes equivalent), as well as to explain further the false paradox between smooth-and-abrupt-bends. This unified framework allows us to improve greatly our understanding of the effect of fault geometry on the mechanics of earthquakes.
Abstract
This paper defines the UK Infra-Red Telescope (UKIRT) Hemisphere Survey (UHS) and release of the remaining ∼12 700 deg2 of J-band survey data products. The UHS will provide continuous J- and ...K-band coverage in the Northern hemisphere from a declination of 0° to 60° by combining the existing Large Area Survey, Galactic Plane Survey and Galactic Clusters Survey conducted under the UKIRT Infra-red Deep Sky Survey (UKIDSS) programme with this new additional area not covered by UKIDSS. The released data include J-band imaging and source catalogues over the new area, which, together with UKIDSS, completes the J-band UHS coverage over the full ∼17 900 deg2 area. 98 per cent of the data in this release have passed quality control criteria. The remaining 2 per cent have been scheduled for re-observation. The median 5σ point source sensitivity of the released data is 19.6 mag (Vega). The median full width at half-maximum of the point spread function across the data set
is 0.75 arcsec. In this paper, we outline the survey management, data acquisition, processing and calibration, quality control and archiving as well as summarizing the characteristics of the released data products. The data are initially available to a limited consortium with a world-wide release scheduled for 2018 August.
Abstract
We investigate the mean star formation rates (SFRs) in the host galaxies of ∼3000 optically selected quasi-stellar objects (QSOs) from the Sloan Digital Sky Survey within the Herschel-ATLAS ...fields, and a radio-luminous subsample covering the redshift range of z = 0.2–2.5. Using Wide-field Infrared Survey Explorer (WISE) and Herschel photometry (12–500 μm) we construct composite spectral energy distributions (SEDs) in bins of redshift and active galactic nucleus (AGN) luminosity. We perform SED fitting to measure the mean infrared luminosity due to star formation, removing the contamination from AGN emission. We find that the mean SFRs show a weak positive trend with increasing AGN luminosity. However, we demonstrate that the observed trend could be due to an increase in black hole (BH) mass (and a consequent increase of inferred stellar mass) with increasing AGN luminosity. We compare to a sample of X-ray selected AGN and find that the two populations have consistent mean SFRs when matched in AGN luminosity and redshift. On the basis of the available virial BH masses, and the evolving BH mass to stellar mass relationship, we find that the mean SFRs of our QSO sample are consistent with those of main sequence star-forming galaxies. Similarly the radio-luminous QSOs have mean SFRs that are consistent with both the overall QSO sample and with star-forming galaxies on the main sequence. In conclusion, on average QSOs reside on the main sequence of star-forming galaxies, and the observed positive trend between the mean SFRs and AGN luminosity can be attributed to BH mass and redshift dependencies.
Abstract
Using high-resolution (sub-kiloparsec scale) data obtained by ALMA, we analyse the star formation rate (SFR), gas content, and kinematics in SDP 81, a gravitationally lensed starburst galaxy ...at redshift 3. We estimate the SFR surface density (ΣSFR) in the brightest clump of this galaxy to be $357^{+135}_{-85}\, \mathrm{M}_{\odot }\, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$, over an area of 0.07 ± 0.02 kpc2. Using the intensity-weighted velocity of CO (5–4), we measure the turbulent velocity dispersion in the plane of the sky and find σv, turb = 37 ± 5 km s−1 for the clump, in good agreement with previous estimates along the line of sight. Our measurements of the gas surface density, freefall time, and turbulent Mach number allow us to compare the theoretical SFR from various star formation models with that observed, revealing that the role of turbulence is crucial to explaining the observed SFR in this clump. While the Kennicutt–Schmidt (KS) relation predicts an SFR surface density of ΣSFR, KS = 52 ± 17 M⊙ yr−1 kpc−2, the single-freefall model by Krumholz, Dekel, and McKee (KDM) predicts ΣSFR, KDM = 106 ± 37 M⊙ yr−1 kpc−2. In contrast, the multifreefall (turbulence) model by Salim, Federrath, and Kewley (SFK) gives $\Sigma _{\mathrm{SFR,SFK}} = 491^{+139 }_{-194}\, \mathrm{M}_{\odot }\, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$. Although the SFK relation overestimates the SFR in this clump (possibly due to the negligence of magnetic fields), it provides the best prediction among the available models. Finally, we compare the star formation and gas properties of this galaxy to local star-forming regions and find that the SFK relation provides the best estimates of SFR in both local and high-redshift galaxies.
We apply a chemical evolution model to investigate the sources and evolution of dust in a sample of 26 high-redshift (z > 1) submillimetre galaxies (SMGs) from the literature, with complete ...photometry from ultraviolet to the submillimetre. We show that dust produced only by low–intermediate-mass stars falls a factor 240 short of the observed dust masses of SMGs, the well-known ‘dust-budget crisis’. Adding an extra source of dust from supernovae can account for the dust mass in 19 per cent of the SMG sample. Even after accounting for dust produced by supernovae the remaining deficit in the dust mass budget provides support for higher supernova yields, substantial grain growth in the interstellar medium or a top-heavy IMF. Including efficient destruction of dust by supernova shocks increases the tension between our model and observed SMG dust masses. The models which best reproduce the physical properties of SMGs have a rapid build-up of dust from both stellar and interstellar sources and minimal dust destruction. Alternatively, invoking a top-heavy IMF or significant changes in the dust grain properties can solve the dust budget crisis only if dust is produced by both low-mass stars and supernovae and is not efficiently destroyed by supernova shocks.
SUMMARY
In a dislocation problem, a paradoxical discordance is known to occur between an original smooth curve and an infinitesimally discretized curve. To solve this paradox, we have investigated a ...non-hypersingular expression for the integral kernel (called the stress Green’s function) which describes the stress field caused by the displacement discontinuity. We first develop a compact alternative expression of the non-hypersingular stress Green’s function for general 2-D and 3-D infinite homogeneous elastic media. We next compute the stress Green’s functions on a curved fault and revisit the paradox. We find that previously obtained non-hypersingular stress Green’s functions are incorrect for curved faults, and that smooth and infinitesimally segmented faults are equivalent. Their compatibility bridges the gap between analytical methods featuring curved faults and numerical methods using subdivided flat patches.