The Surface Water and Ocean Topography (SWOT) satellite mission planned for launch in 2020 will map river elevations and inundated area globally for rivers >100 m wide. In advance of this launch, we ...here evaluated the possibility of estimating discharge in ungauged rivers using synthetic, daily “remote sensing” measurements derived from hydraulic models corrupted with minimal observational errors. Five discharge algorithms were evaluated, as well as the median of the five, for 19 rivers spanning a range of hydraulic and geomorphic conditions. Reliance upon a priori information, and thus applicability to truly ungauged reaches, varied among algorithms: one algorithm employed only global limits on velocity and depth, while the other algorithms relied on globally available prior estimates of discharge. We found at least one algorithm able to estimate instantaneous discharge to within 35% relative root‐mean‐squared error (RRMSE) on 14/16 nonbraided rivers despite out‐of‐bank flows, multichannel planforms, and backwater effects. Moreover, we found RRMSE was often dominated by bias; the median standard deviation of relative residuals across the 16 nonbraided rivers was only 12.5%. SWOT discharge algorithm progress is therefore encouraging, yet future efforts should consider incorporating ancillary data or multialgorithm synergy to improve results.
Key Points:
SWOT discharge algorithms were tested on synthetic observations for 19 rivers
Algorithms accurately characterized temporal dynamics of river discharge
At least one algorithm estimated discharge to <35% relative RMSE on 14/16 of nonbraided rivers
Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, ...it is possible to detect faint components around bright stars in a few hours of observations. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit using all interferometric observables, i.e., the squared visibilities, closure phases and bi-spectrum amplitudes. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting the Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.
Context.
The nearby and young
β
Pictoris system hosts a well resolved disk, a directly imaged massive giant planet orbiting at ≃9 au, as well as an inner planet orbiting at ≃2.7 au, which was ...recently detected through radial velocity (RV). As such, it offers several unique opportunities for detailed studies of planetary system formation and early evolution.
Aims.
We aim to further constrain the orbital and physical properties of
β
Pictoris b and c using a combination of high contrast imaging, long base-line interferometry, and RV data. We also predict the closest approaches or the transit times of both planets, and we constrain the presence of additional planets in the system.
Methods.
We obtained six additional epochs of SPHERE data, six additional epochs of GRAVITY data, and five additional epochs of RV data. We combined these various types of data in a single Markov-chain Monte Carlo analysis to constrain the orbital parameters and masses of the two planets simultaneously. The analysis takes into account the gravitational influence of both planets on the star and hence their relative astrometry. Secondly, we used the RV and high contrast imaging data to derive the probabilities of presence of additional planets throughout the disk, and we tested the impact of absolute astrometry.
Results.
The orbital properties of both planets are constrained with a semi-major axis of 9.8 ± 0.4 au and 2.7 ± 0.02 au for b and c, respectively, and eccentricities of 0.09 ± 0.1 and 0.27 ± 0.07, assuming the H
IPPARCOS
distance. We note that despite these low fitting error bars, the eccentricity of
β
Pictoris c might still be over-estimated. If no prior is provided on the mass of
β
Pictoris b, we obtain a very low value that is inconsistent with what is derived from brightness-mass models. When we set an evolutionary model motivated prior to the mass of
β
Pictoris b, we find a solution in the 10–11
M
Jup
range. Conversely,
β
Pictoris c’s mass is well constrained, at 7.8 ± 0.4
M
Jup
, assuming both planets are on coplanar orbits. These values depend on the assumptions on the distance of the
β
Pictoris system. The absolute astrometry H
IPPARCOS
-
Gaia
data are consistent with the solutions presented here at the 2
σ
level, but these solutions are fully driven by the relative astrometry plus RV data. Finally, we derive unprecedented limits on the presence of additional planets in the disk. We can now exclude the presence of planets that are more massive than about 2.5
M
Jup
closer than 3 au, and more massive than 3.5
M
Jup
between 3 and 7.5 au. Beyond 7.5 au, we exclude the presence of planets that are more massive than 1–2
M
Jup
.
Conclusions.
Combining relative astrometry and RVs allows one to precisely constrain the orbital parameters of both planets and to give lower limits to potential additional planets throughout the disk. The mass of
β
Pictoris c is also well constrained, while additional RV data with appropriate observing strategies are required to properly constrain the mass of
β
Pictoris b.
Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we ...present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best-fit ExoREM model with solar metallicity and C/O = 0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available on GitHub: https://github.com/sblunt/hip65426.
During their formation phase, stars gain most of their mass in violent episodic accretion events, such as observed in FU Orionis (FUor) and EXor stars. V346 Normae is a well-studied FUor that ...underwent a strong outburst beginning around 1980. Here, we report on photometric and spectroscopic observations, which show that the visual/near-infrared brightness has decreased dramatically between the 1990s and 2010 (ΔR ≈ 10.9 mag, ΔJ ≈ 7.8 mag and ΔK ≈ 5.8 mag). The spectral properties of this fading event cannot be explained by variable extinction alone, but indicate a drop in accretion rate by two to three orders of magnitude. This is the first time that a member of the FUor class has been observed to switch to a very low accretion phase. Remarkably, in the last few years (2011–2015) V346 Nor has brightened again at all near-infrared wavelengths, indicating the onset of a new outburst event. The observed behaviour might be consistent with the clustered luminosity bursts that have been predicted by recent gravitational instability and fragmentation models for the early stages of protostellar evolution. Given V346 Nor's unique characteristics (concerning outburst duration, repetition frequency and spectroscopic diagnostics), our results also highlight the need to revisit the FUor/EXor classification scheme.
Abstract
Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses, can provide a way forward. HD 72946 B is a ...recently imaged brown dwarf companion to a nearby, solar-type star. We achieve ∼100
μ
as relative astrometry of HD 72946 B in the
K
band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using
orbitize!
and derive a strong dynamical mass constraint
M
B
= 69.5 ± 0.5
M
Jup
assuming a strong prior on the host star mass
M
A
= 0.97 ± 0.01
M
⊙
from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent
BT-Settl-CIFIST
model atmospheres, and perform atmospheric retrievals using
petitRADTRANS
. A dynamical mass prior only marginally influences the sampled distribution of effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure–temperature or cloud parameters within our current retrieval setup. Independently of the cloud prescription and prior assumptions, we find agreement within ±2
σ
between the C/O of the host (0.52 ± 0.05) and brown dwarf (0.43–0.63), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high (0.6–0.8) in light of the excellent agreement of the data with the solar-abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis.
We present the results of long-baseline optical interferometry observations using the Precision Astronomical Visual Observations (PAVO) beam combiner at the Center for High Angular Resolution ...Astronomy (CHARA) Array to measure the angular sizes of three bright Kepler stars: θ Cygni, and both components of the binary system 16 Cygni. Supporting infrared observations were made with the Michigan Infrared Combiner (MIRC) and Classic beam combiner, also at the CHARA Array. We find limb-darkened angular diameters of 0.753 ± 0.009 mas for θ Cyg, 0.539 ± 0.007 mas for 16 Cyg A and 0.490 ± 0.006 mas for 16 Cyg B. The Kepler Mission has observed these stars with outstanding photometric precision, revealing the presence of solar-like oscillations. Due to the brightness of these stars the oscillations have exceptional signal-to-noise, allowing for detailed study through asteroseismology, and are well constrained by other observations. We have combined our interferometric diameters with Hipparcos parallaxes, spectrophotometric bolometric fluxes and the asteroseismic large frequency separation to measure linear radii (θ Cyg: 1.48 ± 0.02 R, 16 Cyg A: 1.22 ± 0.02 R, 16 Cyg B: 1.12 ± 0.02 R), effective temperatures (θ Cyg: 6749 ± 44 K, 16 Cyg A: 5839 ± 42 K, 16 Cyg B: 5809 ± 39 K) and masses (θ Cyg: 1.37 ± 0.04 M, 16 Cyg A: 1.07 ± 0.05 M, 16 Cyg B: 1.05 ± 0.04 M) for each star with very little model dependence. The measurements presented here will provide strong constraints for future stellar modelling efforts.
Imaging the direct light signal from a faint exoplanet against the overwhelming glare of its host star presents one of the fundamental challenges to modern astronomical instrumentation. Achieving ...sufficient signal-to-noise ratio for detection by direct imaging is limited by three basic physical processes: aberration of the wavefronts (both instrumental and atmospheric), photon noise and detector noise. In this paper, we advance a novel optical setup which synthesizes the advantages of two different techniques: nulling interferometry to mitigate photon noise, and closure phase to combat optical aberrations. Our design, which employs technology from integrated optics and photonics, is intended to combine the advantageous aspects of both a coronagraph and a non-redundant interferometer inside a single optical device. We show that such an instrument would have a dynamic range limited by either (i) the readout noise (if perfect co-phasing), or (ii) the photon noise due to stellar flux leakage (in the case of imperfect nulling). This concept is optimal when the readout noise is not the main limitation, i.e. for space interferometry or for ground-based observations of bright stellar hosts (apparent magnitude brighter than 10).
Abstract
We present 18 yr of OGLE photometry together with spectra obtained over 12 yr revealing that the early Oe star AzV 493 shows strong photometric (Δ
I
< 1.2 mag) and spectroscopic variability ...with a dominant, 14.6 yr pattern and ∼40 day oscillations. We estimate the stellar parameters
T
eff
= 42,000 K,
log
L
/
L
⊙
=
5.83
±
0.15
,
M
/
M
⊙
= 50 ± 9, and
v
sin
i
= 370 ± 40 km s
−1
. Direct spectroscopic evidence shows episodes of both gas ejection and infall. There is no X-ray detection, and it is likely a runaway star. The star AzV 493 may have an unseen companion on a highly eccentric (
e
> 0.93) orbit. We propose that close interaction at periastron excites ejection of the decretion disk, whose variable emission-line spectrum suggests separate inner and outer components, with an optically thick outer component obscuring both the stellar photosphere and the emission-line spectrum of the inner disk at early phases in the photometric cycle. It is plausible that AzV 493’s mass and rotation have been enhanced by binary interaction followed by the core-collapse supernova explosion of the companion, which now could be either a black hole or a neutron star. This system in the Small Magellanic Cloud can potentially shed light on OBe decretion disk formation and evolution, massive binary evolution, and compact binary progenitors.
Context.
FUors outbursts are a crucial stage of accretion in young stars. However, a complete mechanism at the origin of the outburst still remains missing.
Aims.
We aim to constrain the instability ...mechanism in the star FU Orionis itself by directly probing the size and evolution in time of the outburst region with near-infrared (NIR) interferometry, and to confront it with physical models of this region.
Methods.
As the prototype of the FUors class of objects, FU Orionis has been a regular target of NIR interferometry. In this paper, we analyze more than 20 years of NIR interferometric observations to perform a temporal monitoring of the region of the outburst, and compare it to the spatial structure deduced from 1D magneto-hydrodynamic (MHD) simulations.
Results.
We measure from the interferometric observations that the size variation of the outburst region is compatible with a constant or slightly decreasing size over time: -0.56
-0.36
+0.14
AU/100 yr and -0.30
-0.19
+0.19
AU/100 yr in the
H
and
K
bands, respectively. The temporal variation and the mean size probed by NIR interferometry are consistently reproduced by our 1D MHD simulations. We find that the most compatible scenario is a model of an outburst occurring in a magnetically layered disk, where a magneto-rotational instability (MRI) is triggered by a gravitational instability (GI) at the outer edge of a dead zone. The scenario of a pure thermal instability (TI) fails to reproduce our interferometric sizes because it can only be sustained in a very compact zone of the disk <0.1 AU. The comparison between the data and the MRI-GI models favors MHD parameters of
α
MRI
= 10
−2
,
T
MRI
= 800 K, and Σ
crit
= 10 g cm
−2
, with more work needed in terms of observations and modeling in order to improve the precision of these values. Locally, in the very inner part of the disk, TI can be triggered in addition to MRI-GI, which qualitatively better matches our observation but is not strongly constrained by the currently available data. The scenario of MRI-GI could be compatible with an external perturbation that enhances the GI, such as tidal interaction with a stellar companion, or a planet at the outer edge of the dead zone.
Conclusions.
We favor a layered-disk model driven by MRI turbulence in order to explain the spatial structure and temporal evolution of the outburst region on FU Orionis. Understanding this phase will provide a crucial link between the early phase of disk evolution and the process of planet formation in the first inner astronomical units.