skip to main content


Title: Photometry on Structured Backgrounds: Local Pixel-wise Infilling by Regression
Abstract

Photometric pipelines struggle to estimate both the flux and flux uncertainty for stars in the presence of structured backgrounds such as filaments or clouds. However, it is exactly stars in these complex regions that are critical to understanding star formation and the structure of the interstellar medium. We develop a method, similar to Gaussian process regression, which we term local pixel-wise infilling (LPI). Using a local covariance estimate, we predict the background behind each star and the uncertainty of that prediction in order to improve estimates of flux and flux uncertainty. We show the validity of our model on synthetic data and real dust fields. We further demonstrate that the method is stable even in the crowded field limit. While we focus on optical-IR photometry, this method is not restricted to those wavelengths. We apply this technique to the 34 billion detections in the second data release of the Dark Energy Camera Plane Survey. In addition to removing many >3σoutliers and improving uncertainty estimates by a factor of ∼2–3 on nebulous fields, we also show that our method is well behaved on uncrowded fields. The entirely post-processing nature of our implementation of LPI photometry allows it to easily improve the flux and flux uncertainty estimates of past as well as future surveys.

 
more » « less
Award ID(s):
2019786
NSF-PAR ID:
10369014
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
933
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 155
Size(s):
["Article No. 155"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The bright starλSer hosts a hot Neptune with a minimum mass of 13.6Mand a 15.5 day orbit. It also appears to be a solar analog, with a mean rotation period of 25.8 days and surface differential rotation very similar to the Sun. We aim to characterize the fundamental properties of this system and constrain the evolutionary pathway that led to its present configuration. We detect solar-like oscillations in time series photometry from the Transiting Exoplanet Survey Satellite, and we derive precise asteroseismic properties from detailed modeling. We obtain new spectropolarimetric data, and we use them to reconstruct the large-scale magnetic field morphology. We reanalyze the complete time series of chromospheric activity measurements from the Mount Wilson Observatory, and we present new X-ray and ultraviolet observations from the Chandra and Hubble space telescopes. Finally, we use the updated observational constraints to assess the rotational history of the star and estimate the wind braking torque. We conclude that the remaining uncertainty on the stellar age currently prevents an unambiguous interpretation of the properties ofλSer, and that the rate of angular momentum loss appears to be higher than for other stars with a similar Rossby number. Future asteroseismic observations may help to improve the precision of the stellar age.

     
    more » « less
  2. Abstract

    Calcium-rich (Ca-rich) transients are a new class of supernovae (SNe) that are known for their comparatively rapid evolution, modest peak luminosities, and strong nebular calcium emission lines. Currently, the progenitor systems of Ca-rich transients remain unknown. Although they exhibit spectroscopic properties not unlike core-collapse Type Ib/c SNe, nearly half are found in the outskirts of their host galaxies, which are predominantly elliptical, suggesting a closer connection to the older stellar populations of SNe Ia. In this paper, we present a compilation of publicly available multiwavelength observations of all known and/or suspected host galaxies of Ca-rich transients ranging from far-UV to IR, and use these data to characterize their stellar populations withprospector. We estimate several galaxy parameters including integrated star formation rate, stellar mass, metallicity, and age. For nine host galaxies, the observations are sensitive enough to obtain nonparametric star formation histories, from which we recover SN rates and estimate probabilities that the Ca-rich transients in each of these host galaxies originated from a core-collapse versus Type Ia-like explosion. Our work supports the notion that the population of Ca-rich transients do not come exclusively from core-collapse explosions, and must either be only from white dwarf stars or a mixed population of white dwarf stars with other channels, potentially including massive star explosions. Additional photometry and explosion site spectroscopy of larger samples of Ca-rich host galaxies will improve these estimates and better constrain the ratio of white dwarf versus massive star progenitors of Ca-rich transients.

     
    more » « less
  3. Abstract

    Estimating phenotypic distributions of populations and communities is central to many questions in ecology and evolution. These distributions can be characterized by their moments (mean, variance, skewness and kurtosis) or diversity metrics (e.g. functional richness). Typically, such moments and metrics are calculated using community‐weighted approaches (e.g. abundance‐weighted mean). We propose an alternative bootstrapping approach that allows flexibility in trait sampling and explicit incorporation of intraspecific variation, and show that this approach significantly improves estimation while allowing us to quantify uncertainty.

    We assess the performance of different approaches for estimating the moments of trait distributions across various sampling scenarios, taxa and datasets by comparing estimates derived from simulated samples with the true values calculated from full datasets. Simulations differ in sampling intensity (individuals per species), sampling biases (abundance, size), trait data source (local vs. global) and estimation method (two types of community‐weighting, two types of bootstrapping).

    We introduce thetraitstrapR package, which contains a modular and extensible set of bootstrapping and weighted‐averaging functions that use community composition and trait data to estimate the moments of community trait distributions with their uncertainty. Importantly, the first function in the workflow,trait_fill, allows the user to specify hierarchical structures (e.g. plot within site, experiment vs. control, species within genus) to assign trait values to each taxon in each community sample.

    Across all taxa, simulations and metrics, bootstrapping approaches were more accurate and less biased than community‐weighted approaches. With bootstrapping, a sample size of 9 or more measurements per species per trait generally included the true mean within the 95% CI. It reduced average percent errors by 26%–74% relative to community‐weighting. Random sampling across all species outperformed both size‐ and abundance‐biased sampling.

    Our results suggest randomly sampling ~9 individuals per sampling unit and species, covering all species in the community and analysing the data using nonparametric bootstrapping generally enable reliable inference on trait distributions, including the central moments, of communities. By providing better estimates of community trait distributions, bootstrapping approaches can improve our ability to link traits to both the processes that generate them and their effects on ecosystems.

     
    more » « less
  4. Abstract

    We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the Hubble Space Telescope–based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300–9400 Å as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all presubmission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures ≳5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12 mag range. Overall, the mean of our system is calibrated to the mean of CALSPEC at the level of ∼3 mmag. With our large number of observations, careful cross-checks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys.

     
    more » « less
  5. ABSTRACT

    We present SCUBA-2 850 $\mathrm{ \mu}$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850 $\mathrm{ \mu}$m sources in 9 of 13 of these candidate protoclusters show significant overdensities compared to the field, with the probability <10−2 assuming the sources are randomly distributed in the sky. Using the 250, 350, 500, and 850 $\mathrm{ \mu}$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution templates with an MCMC method. The photometric redshift distribution, peaking at 2 < z < 3, is consistent with that of known z > 2 protoclusters and the peak of the cosmic star formation rate density (SFRD). We find that the 850 $\mathrm{ \mu}$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim 10^{12}\, \mathrm{L}_{\odot }$ and star formation rates of SFR  = (500–1500) M⊙ yr−1. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorized into four groups: six of them being high-redshift starbursting protoclusters, one being a lower redshift cluster or protocluster, three being protoclusters that contain lensed dusty star-forming galaxies or are rich in 850 $\mathrm{ \mu}$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850 $\mathrm{ \mu}$m sources.

     
    more » « less