skip to main content


Title: Inferring Latent Velocities from Weather Radar Data using Gaussian Processes
Archived data from the US network of weather radars hold detailed information about bird migration over the last 25 years, including very high-resolution partial measurements of velocity. Historically, most of this spatial resolution is discarded and velocities are summarized at a very small number of locations due to modeling and algorithmic limitations. This paper presents a Gaussian process (GP) model to reconstruct high-resolution full velocity fields across the entire US. The GP faithfully models all aspects of the problem in a single joint framework, including spatially random velocities, partial velocity measurements, station-specific geometries, measurement noise, and an ambiguity known as aliasing. We develop fast inference algorithms based on the FFT; to do so, we employ a creative use of Laplace's method to sidestep the fact that the kernel of the joint process is non-stationary.  more » « less
Award ID(s):
1661259
NSF-PAR ID:
10092572
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Advances in neural information processing systems
Volume:
31
ISSN:
1049-5258
Page Range / eLocation ID:
8984–8993
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. We present Glacier Image Velocimetry (GIV), an open-source and easy-to-use software toolkit for rapidly calculating high-spatial-resolutionglacier velocity fields. Glacier ice velocity fields reveal flow dynamics, ice-flux changes, and (with additional data and modelling) icethickness. Obtaining glacier velocity measurements over wide areas with field techniques is labour intensive and often associated with safetyrisks. The recent increased availability of high-resolution, short-repeat-time optical imagery allows us to obtain ice displacement fields using“feature tracking” based on matching persistent irregularities on the ice surface between images and hence, surface velocity over time. GIV isfully parallelized and automatically detects, filters, and extracts velocities from large datasets of images. Through this coupled toolchain and aneasy-to-use GUI, GIV can rapidly analyse hundreds to thousands of image pairs on a laptop or desktop computer. We present four example applicationsof the GIV toolkit in which we complement a glaciology field campaign (Glaciar Perito Moreno, Argentina) and calculate the velocity fields of smallmid-latitude (Glacier d'Argentière, France) and tropical glaciers (Volcán Chimborazo, Ecuador), as well as very large glaciers (Vavilov Ice Cap,Russia). Fully commented MATLAB code and a stand-alone app for GIV are available from GitHub and Zenodo (see https://doi.org/10.5281/zenodo.4624831, Van Wyk de Vries, 2021a). 
    more » « less
  2. Context. The recent and exquisite astrometric, photometric, and radial velocity measurements of the Gaia mission resulted in a substantial advancement of the determination of the orbits for old star clusters, including the oldest Milky Way globular clusters (MW GCs). Aims. The main goal of the present paper is to use the new Gaia data release 3 (DR3) and the VISTA Variables in the Via Láctea Extended Survey (VVVX) measurements to obtain the orbits for nearly a dozen new MW GC candidates that have been poorly studied or previously unexplored. Methods. We use the Gaia DR3 and VVVX databases to identify bona fide MW GC candidates, namely VVV-CL160, Patchick 122, Patchick 125, Patchick 126, Kronberger 99, Kronberger 119, Kronberger 143, ESO 92-18, ESO 93-08, Gaia 2, and Ferrero 54. The relevant mean cluster physical parameters are derived (distances, Galactic coordinates, proper motions, radial velocities). We also measure accurate mean radial velocities for the GCs VVV-CL160 and Patchick 126 using observations acquired at the Gemini-South telescope with the Immersion GRating INfrared Spectrometer (IGRINS) high-resolution spectrograph. Orbits for each cluster are then computed using the GravPot16 model, assuming typical Galactic bar pattern speeds. Results. We reconstruct the orbits for these 11 star clusters for the first time. These include star clusters with retrograde and prograde orbital motions, both in the Galactic bulge and disk. We obtain orbital properties for this sample, such as the mean time-variations of perigalactic and apogalactic distances, eccentricities, vertical excursions from the Galactic plane, and Z -components of the angular momentum. Conclusions. Our main conclusion is that, based on the orbital parameters, Patchick 125 and Patchick 126 are genuine MW bulge or halo GCs; and Ferrero 54, Gaia 2, and Patchick 122 are MW disk GCs. In contrast, the orbits of Kronberger 99, Kronberger 119, Kronberger 143, ESO 92-18, and ESO 93-08 are more consistent with old MW disk open clusters, in agreement with previous results. VVV-CL160 falls very close to the Galactic centre, but reaches larger distances beyond the Solar orbit, and therefore its origin is still unclear. 
    more » « less
  3. Abstract

    To investigate the controls on continental rifting in the western branch of the East Africa Rift System, we conduct shear velocity imaging of the crust and uppermost mantle beneath the weakly extended Malawi Rift and the Rungwe Volcanic Province (RVP). We use local‐scale measurements of Rayleigh wave phase velocities between 9‐ and 100‐s periods combined with constraints on basin architecture and crustal thickness to invert for shear velocity from the surface to ~135 km. Our resulting 3‐D model reveals a localized low‐velocity anomaly associated with the RVP extending from the crust and through the upper mantle, which can be explained with modestly elevated temperatures. Away from the RVP, velocities within mantle flanking the rift are fast (>4.6 ± 0.1 km/s), suggesting depleted lithospheric mantle to depths of ~100 and >135 km to the west and east of the rift, respectively. The upper mantle beneath the rift axis is characterized by thinned lithosphere with slower velocities than the surrounding plateau, suggestive of thermal and/or chemical modification by the rifting process. Slowest velocities are mildly asymmetric about the rift axis, with the lowest velocities observed beneath the rift and adjacent footwall escarpments. The underlying asthenosphere is only moderately slow (~4.25 ± 0.1 km/s), including beneath the RVP, precluding the presence of significant volumes of partial melt. The positions of localized lithospheric modification and basin‐bounding border faults correlate with the location of Proterozoic mobile belts, suggesting that these sutures provide lithospheric‐scale weakening mechanisms necessary for localizing strain and allowing extension to occur in the Malawi Rift.

     
    more » « less
  4. null (Ed.)
    ABSTRACT It has been proposed that gravothermal collapse due to dark matter self-interactions (i.e. self-interacting dark matter, SIDM) can explain the observed diversity of the Milky Way (MW) satellites’ central dynamical masses. We investigate the process behind this hypothesis using an N-body simulation of a MW-analogue halo with velocity-dependent SIDM (vdSIDM) in which the low-velocity self-scattering cross-section, $\sigma _{\rm T}/m_{\rm x}$, reaches 100 cm2 g−1; we dub this model the vd100 model. We compare the results of this simulation to simulations of the same halo that employ different dark models, including cold dark matter (CDM) and other, less extreme SIDM models. The masses of the vd100 haloes are very similar to their CDM counterparts, but the values of their maximum circular velocities, Vmax, are significantly higher. We determine that these high Vmax subhaloes were objects in the mass range [5 × 106, 1 × 108] M⊙ at z = 1 that undergo gravothermal core collapse. These collapsed haloes have density profiles that are described by single power laws down to the resolution limit of the simulation, and the inner slope of this density profile is approximately −3. Resolving the ever decreasing collapsed region is challenging, and tailored simulations will be required to model the runaway instability accurately at scales <1 kpc. 
    more » « less
  5. Exploiting the sensitivity of the IRAM NOrthern Extended Millimeter Array (NOEMA) and its ability to process large instantaneous bandwidths, we have studied the morphology and other properties of the molecular gas and dust in the star forming galaxy, H-ATLAS J131611.5+281219 (HerBS-89a), at z = 2.95. High angular resolution (0 . ″3) images reveal a partial 1 . ″0 diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12 CO(9−8) and H 2 O(2 02  − 1 11 ). Together with lower angular resolution (0 . ″6) images, we report the detection of a series of molecular lines including the three fundamental transitions of the molecular ion OH + , namely (1 1  − 0 1 ), (1 2  − 0 1 ), and (1 0  − 0 1 ), seen in absorption; the molecular ion CH + (1 − 0) seen in absorption, and tentatively in emission; two transitions of amidogen (NH 2 ), namely (2 02  − 1 11 ) and (2 20  − 2 11 ) seen in emission; and HCN(11 − 10) and/or NH(1 2  − 0 1 ) seen in absorption. The NOEMA data are complemented with Very Large Array data tracing the 12 CO(1 − 0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present Hubble Space Telescope imaging that reveals the foreground lensing galaxy in the near-infrared (1.15  μ m). Together with photometric data from the Gran Telescopio Canarias, we derive a photometric redshift of z phot = 0.9 −0.5 +0.3 for the foreground lensing galaxy. Modeling the lensing of HerBS-89a, we reconstruct the dust continuum (magnified by a factor μ  ≃ 5.0) and molecular emission lines (magnified by μ  ∼ 4 − 5) in the source plane, which probe scales of ∼0 . ″1 (or 800 pc). The 12 CO(9 − 8) and H 2 O(2 02  − 1 11 ) emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter, which reveals two compact rotating components with sizes of ≈1 kpc that are likely merging, more naturally accounts for the broad line widths observed in HerBS-89a. In the core of HerBS-89a, very dense gas with n H 2  ∼ 10 7 − 9 cm −3 is revealed by the NH 2 emission lines and the possible HCN(11 − 10) absorption line. HerBS-89a is a powerful star forming galaxy with a molecular gas mass of M mol  = (2.1 ± 0.4) × 10 11   M ⊙ , an infrared luminosity of L IR  = (4.6 ± 0.4) × 10 12   L ⊙ , and a dust mass of M dust  = (2.6 ± 0.2) × 10 9   M ⊙ , yielding a dust-to-gas ratio δ GDR  ≈ 80. We derive a star formation rate SFR = 614 ± 59  M ⊙ yr −1 and a depletion timescale τ depl  = (3.4 ± 1.0) × 10 8 years. The OH + and CH + absorption lines, which trace low (∼100 cm −3 ) density molecular gas, all have their main velocity component red-shifted by Δ V  ∼ 100 km s −1 relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow toward the center of a galaxy, indicating that HerBS-89a is accreting gas from its surroundings. 
    more » « less