skip to main content


Title: Spatially Variable Advection Correction of Doppler Radial Velocity Data
Abstract Techniques to mitigate analysis errors arising from the nonsimultaneity of data collections typically use advection-correction procedures based on the hypothesis (frozen turbulence) that the analyzed field can be represented as a pattern of unchanging form in horizontal translation. It is more difficult to advection correct the radial velocity than the reflectivity because even if the vector velocity field satisfies this hypothesis, its radial component does not—but that component does satisfy a second-derivative condition. We treat the advection correction of the radial velocity ( υ r ) as a variational problem in which errors in that second-derivative condition are minimized subject to smoothness constraints on spatially variable pattern-translation components ( U , V ). The Euler–Lagrange equations are derived, and an iterative trajectory-based solution is developed in which U , V , and υ r are analyzed together. The analysis code is first verified using analytical data, and then tested using Atmospheric Imaging Radar (AIR) data from a band of heavy rainfall on 4 September 2018 near El Reno, Oklahoma, and a decaying tornado on 27 May 2015 near Canadian, Texas. In both cases, the analyzed υ r field has smaller root-mean-square errors and larger correlation coefficients than in analyses based on persistence, linear time interpolation, or advection correction using constant U and V . As some experimentation is needed to obtain appropriate parameter values, the procedure is more suitable for non-real-time applications than use in an operational setting. In particular, the degree of spatial variability in U and V , and the associated errors in the analyzed υ r field are strongly dependent on a smoothness parameter.  more » « less
Award ID(s):
1623626 1823478
NSF-PAR ID:
10212195
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
Volume:
78
Issue:
1
ISSN:
0022-4928
Page Range / eLocation ID:
167 to 188
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    2MASS J20395358+4222505  is an obscured early B supergiant near the massive OB star association Cygnus OB2. Despite its bright infrared magnitude (Ks = 5.82) it has remained largely ignored because of its dim optical magnitude (B = 16.63, V = 13.68). In a previous paper, we classified it as a highly reddened, potentially extremely luminous, early B-type supergiant. We obtained its spectrum in the U, B and R spectral bands during commissioning observations with the instrument MEGARA at the Gran Telescopio CANARIAS. It displays a particularly strong Hα emission for its spectral type, B1 Ia. The star seems to be in an intermediate phase between supergiant and hypergiant, a group that it will probably join in the near (astronomical) future. We observe a radial velocity difference between individual observations and determine the stellar parameters, obtaining Teff = 24 000 K and log gc = 2.88 ± 0.15. The rotational velocity found is large for a B supergiant, $v$ sin i = 110 ± 25 $\rm km\, s^{-1}$. The abundance pattern is consistent with solar, with a mild C  underabundance (based on a single line). Assuming that J20395358+4222505  is at the distance of Cyg OB2, we derive the radius from infrared photometry, finding R = 41.2 ± 4.0 R⊙, log(L/L⊙) = 5.71 ± 0.04 and a spectroscopic mass of 46.5 ± 15.0 M⊙. The clumped mass-loss rate (clumping factor 10) is very high for the spectral type, $\dot{M}$ = 2.4 × 10−6 M⊙ a−1. The high rotational velocity and mass-loss rate place the star at the hot side of the bi-stability jump. Together with the nearly solar CNO abundance pattern, they may also point to evolution in a binary system, J20395358+4222505  being the initial secondary.

     
    more » « less
  2. . (Ed.)
    This paper proposes a new approach for the adaptive functional estimation of second order infinite dimensional systems with structured perturbations. First, the proposed observer is formulated in the natural second order setting thus ensuring the time derivative of the estimated position is the estimated velocity, and therefore called natural adaptive observer. Assuming that the system does not yield a positive real system when placed in first order form, then the next step in deriving parameter adaptive laws is to assume a form of input-output collocation. Finally, to estimate structured perturbations taking the form of functions of the position and/or velocity outputs, the parameter space is not identified by a finite dimensional Euclidean space but instead is considered in a Reproducing Kernel Hilbert Space. Such a setting allows one not to be restricted by a priori assumptions on the dimension of the parameter spaces. Convergence of the position and velocity errors in their respective norms is established via the use of a parameter-dependent Lyapunov function, specifically formulated for second order infinite dimensional systems that include appropriately defined norms of the functional errors in the reproducing kernel Hilbert spaces. Boundedness of the functional estimates immediately follow and via an appropriate definition of a persistence of excitation condition for functional estimation, a functional convergence follows. When the system is governed by vector second order dynamics, all abstract spaces for the state evolution collapse to a Euclidean space and the natural adaptive observer results simplify. Numerical results of a second order PDE and a multi-degree of freedom finite dimensional mechanical system are presented. 
    more » « less
  3. The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$ . Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$  is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$ . The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs. 
    more » « less
  4. Abstract We present photometric evidence for multiple stellar populations (MPs) in 14 globular clusters (GCs) toward the southern Galactic bulge. The photometric data come as part of the Blanco DECam Bulge Survey, which is a deep, wide-field near-UV-near-IR ( ugriz Y) survey of the southern Galactic bulge. Here, we present the first systematic study of bulge GC multiple populations with deep photometry including the u band, which is a crucial indicator of the abundance of CNO-bearing molecules in stellar atmospheres. We identify cluster members using Gaia EDR3 proper motion measurements, and then isolate red giant branch stars using r versus u − r color–magnitude diagrams. We find evidence suggesting all 14 clusters host at least two populations, and NGC 6441, NGC 6626, and NGC 6656 appear to have at least three populations. Many of these clusters are not part of the Hubble Space Telescope (HST) surveys nor do they have comprehensive spectroscopic analyses so we are presenting the first evidence of MPs in several clusters. Not only do we find a strong anticorrelation between the fraction of first-generation stars and cluster absolute V magnitude, but the correlation coefficient and cluster-to-cluster scatter are similar to the results obtained from HST. Our ground-based data extend to much larger radial distances than similar HST observations, enabling a reliable estimate of the global fraction of first-generation stars in each cluster. This study demonstrates that ground-based u -band photometry as provided by DECam will prove powerful in the study of multiple populations in resolved GCs. 
    more » « less
  5. Abstract In this Letter, we report observations of magnetic switchback (SB) features near 1 au using data from the Wind spacecraft. These features appear to be strikingly similar to the ones observed by the Parker Solar Probe mission closer to the Sun: namely, one-sided spikes (or enhancements) in the solar-wind bulk speed V that correlate/anticorrelate with the spikes seen in the radial-field component B R . In the solar-wind streams that we analyzed, these specific SB features near 1 au are associated with large-amplitude Alfvénic oscillations that propagate outward from the Sun along a local background (prevalent) magnetic field B 0 that is nearly radial. We also show that, when B 0 is nearly perpendicular to the radial direction, the large-amplitude Alfvénic oscillations display variations in V that are two sided (i.e., V alternately increases and decreases depending on the vector Δ B = B − B 0 ). As a consequence, SBs may not always appear as one-sided spikes in V , especially at larger heliocentric distances where the local background field statistically departs from the radial direction. We suggest that SBs can be well described by large-amplitude Alfvénic fluctuations if the field rotation is computed with respect to a well-determined local background field that, in some cases, may deviate from the large-scale Parker field. 
    more » « less