Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clustering statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

Measuring cosmological distances using cluster edges as a standard ruler

ABSTRACT The line-of-sight velocity dispersion profile of galaxy clusters exhibits a ‘kink’ corresponding to the spatial extent of orbiting galaxies. Because the spatial extent of a cluster is correlated with the amplitude of the velocity dispersion profile, we can utilize this feature as a gravity-calibrated standard ruler. Specifically, the amplitude of the velocity dispersion data allows us to infer the physical cluster size. Consequently, observations of the angular scale of the ‘kink’ in the profile can be translated into a distance measurement to the cluster. Assuming the relation between cluster radius and cluster velocity dispersion can be calibrated from simulations, we forecast that with existing data from the Sloan Digital Sky Survey we will be able to measure the Hubble constant with 3.0 per cent precision. Implementing our method with data from the Dark Energy Spectroscopic Instrument (DESI) will result in a 1.3 per cent measurement of the Hubble constant. Adding cosmological supernova data improves the uncertainty of the DESI measurement to 0.7 per cent. While these error estimates are statistical only, they provide strong motivation for pursuing the necessary simulation program required to characterize and calibrate the systematic uncertainties impacting our proposed measurement. Whether or not our proposed measurement can in fact result in competitive H0 constraints will depend on what the eventual systematics floor for this method is.

79 ASTRONOMY AND ASTROPHYSICS↗

Fuzzy set methods for object recognition in space applications

Progress on the following tasks is reported: (1) fuzzy set-based decision making methodologies; (2) feature calculation; (3) clustering for curve and surface fitting; and (4) acquisition of images. The general structure for networks based on fuzzy set connectives which are being used for information fusion and decision making in space applications is described. The structure and training techniques for such networks consisting of generalized means and gamma-operators are described. The use of other hybrid operators in multicriteria decision making is currently being examined. Numerous classical features on image regions such as gray level statistics, edge and curve primitives, texture measures from cooccurrance matrix, and size and shape parameters were implemented. Several fractal geometric features which may have a considerable impact on characterizing cluttered background, such as clouds, dense star patterns, or some planetary surfaces, were used. A new approach to a fuzzy C-shell algorithm is addressed. NASA personnel are in the process of acquiring suitable simulation data and hopefully videotaped actual shuttle imagery. Photographs have been digitized to use in the algorithms. Also, a model of the shuttle was assembled and a mechanism to orient this model in 3-D to digitize for experiments on pose estimation is being constructed.

Keller, James M.↗

Constraints on the Correlation of IceCube Neutrinos with a Tracer of Nearby Large-scale Structure

The IceCube Neutrino Observatory has observed extragalactic astrophysical neutrinos with an apparently isotropic distribution. Only a small fraction of the observed astrophysical neutrinos can be explained by known sources. Neutrino production is thought to occur in energetic environments that are ultimately powered by the gravitational collapse of dense regions of the large-scale mass distribution in the universe. Whatever their identity, neutrino sources likely trace this large-scale mass distribution. The clustering of neutrinos with a tracer of the large-scale structure may provide insight into the distribution of neutrino sources with respect to redshift and the identity of neutrino sources. We implement a two-point angular cross correlation of the Northern sky track events with an infrared galaxy catalog derived from the Wide-field Infrared Survey Explorer (WISE) and Two Micron All Sky Survey (2MASS) source catalogs, which trace the nearby large-scale structure. No statistically significant correlation is found between the neutrinos and this infrared galaxy catalog. We find that ≤54% of the diffuse muon neutrino flux can be attributed to sources correlated with the galaxy catalog with 90% confidence. Additionally, when assuming that the neutrino source comoving density evolves following a power law in redshift, dN s /dV ∝ (1 + z) k , we find that sources with negative evolution, in particular k < −1.75, are disfavored at the 90% confidence level.

Abbasi, R. [Loyola University, Chicago, IL (United↗

Automated Storm Tracking and the Lightning Jump Algorithm Using GOES-R Geostationary Lightning Mapper (GLM) Proxy Data

This study develops a fully automated lightning jump system encompassing objective storm tracking, Geostationary Lightning Mapper proxy data, and the lightning jump algorithm (LJA), which are important elements in the transition of the LJA concept from a research to an operational based algorithm. Storm cluster tracking is based on a product created from the combination of a radar parameter (vertically integrated liquid, VIL), and lightning information (flash rate density). Evaluations showed that the spatial scale of tracked features or storm clusters had a large impact on the lightning jump system performance, where increasing spatial scale size resulted in decreased dynamic range of the system's performance. This framework will also serve as a means to refine the LJA itself to enhance its operational applicability. Parameters within the system are isolated and the system's performance is evaluated with adjustments to parameter sensitivity. The system's performance is evaluated using the probability of detection (POD) and false alarm ratio (FAR) statistics. Of the algorithm parameters tested, sigma-level (metric of lightning jump strength) and flash rate threshold influenced the system's performance the most. Finally, verification methodologies are investigated. It is discovered that minor changes in verification methodology can dramatically impact the evaluation of the lightning jump system.

lightning jump↗

The Diurnal Cycle of Precipitation and Organized Convection in a Set of Global Mesoscale Simulations with the NASA GEOS AGCM

The NASA Global Modeling and Assimilation Office (GMAO) has conducted a series of 40-day nonhydrostatic global simulations with horizontal grid spacing ranging from 200 km to 3 km, as part of the DYAMOND model intercomparison project. The Goddard Earth Observing System (GEOS) model was run with the Grell-Freitas scale-aware convection scheme, which smoothly reduces parameterized deep convection with increasing resolution. Here we evaluate the diurnal cycle and other statistics of precipitation and organized convection as a function of resolution. For validation we use the 0.1 degree IMERG precipitation and 4 km Merged IR brightness temperature datasets, focusing on four regions: the continental United States, Amazonia, the equatorial Indian ocean, and the Maritime Continent. Early results indicate good phase agreement but excessive magnitude of the continental diurnal cycle of precipitation at coarser resolutions, with improved magnitude as resolution increases and the role of parameterization is reduced. Convective cloud clusters are identified with a brightness temperature threshold, and we find realistic numbers of the largest clusters (>10^4 km sq) at all resolutions, while the number of smaller clusters increases with resolution, approaching observations when dx=3 km. The observed diurnal cycle in the cluster size distribution is also reproduced, with realistic magnitude in the highest resolution runs. Precipitation characteristics across cluster sizes are also examined. The results show the potential for global mesoscale simulations as a community science resource, and we invite collaboration to explore these runs in greater detail.

Arnold, Nathan↗

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

Here, we investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core–valence separated (FC-CVS) EOM-CCSD approach – a successful variant of EOM-CCSD tailored towards core excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Venus global radar reflectivity and correlations with elevation

A geological analysis of the global radar reflectivity data for Venus is presented. The statistical and map pattern distribution of the radar reflectivity rho is first analyzed in the context of the porosity structure of the surface. Rho is then correlated with elevation in order to examine the degree of topographic control on material properties which could depend on PT-sensitive chemical weathering phenomena. Comparison of the rms slope alpha and rho correlations with elevation demonstrates that there are distinct topographic zones defined on the basis of the rho and alpha trends. Finally, hierarchical clustering analysis is applied to the global rho, alpha, and 100-km baseline regional slope data sets in order to define radar property-related subregions.

Garvin, J. B.↗

A structural-stochastic model for the analysis and synthesis of cloud images

A structural-stochastic image model is developed for the analysis and synthesis of cloud images. The ability of the model to characterize the visual appearance of cloud fields observed by satellite with a limited number of parameters is demonstrated. The model merges structural and stochastic information, the stochastic model acting as a local statistical operator applied to the output of the structural model. The structural or large-scale organization of the scene is retrieved from the two-dimensional Fourier representation of the digital image. The pattern generated by the major Fourier components provides a first guess of the scene. The stochastic aspect is described by a Markov model of texture that assumes a binomial probability distribution for the local grey-level variability. This Markov model provides four parameters that represent the clustering strength in the horizontal, vertical and diagonal directions. These parameters are estimated by a standard maximum-likelihood technique. The image can be reproduced with a fair degree of verisimilitude from these parameters. The data compression factor is of the order of one hundred to several hundreds.

Garand, L.↗

A Shift Selection Strategy for Parallel Shift-invert Spectrum Slicing in Symmetric Self-consistent Eigenvalue Computation

The central importance of large-scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world’s largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions that will be explored in future work.

97 MATHEMATICS AND COMPUTING↗

Near-infrared Census of RR Lyrae Variables in the Messier 3 Globular Cluster and the Period–Luminosity Relations

We present new near-infrared (NIR), JHK {sub s}, time-series observations of RR Lyrae variables in the Messier 3 (NGC 5272) globular cluster using the WIRCam instrument at the 3.6 m Canada–France–Hawaii Telescope. Our observations cover a sky area of ∼21′ × 21′ around the cluster center and provide an average of 20 epochs of homogeneous JHK {sub s}-band photometry. New homogeneous photometry is used to estimate robust mean magnitudes for 175 fundamental-mode (RRab), 47 overtone-mode (RRc), and 11 mixed-mode (RRd) variables. Our sample of 233 RR Lyrae variables is the largest thus far obtained in a single cluster with time-resolved, multiband NIR photometry. NIR-to-optical amplitude ratios for RR Lyrae in Messier 3 exhibit a systematic increase moving from RRc to short-period (P < 0.6 day) and long-period (P ≳ 0.6 day) RRab variables. We derive JHK {sub s}-band period–luminosity relations for RRab, RRc, and the combined sample of variables. Absolute calibrations based on the theoretically predicted period–luminosity–metallicity relations for RR Lyrae stars yield a distance modulus, μ=15.041±0.017 (statistical)±0.036 (systematic) mag, to Messier 3. When anchored to trigonometric parallaxes for nearby RR Lyrae stars from the Hubble Space Telescope and the Gaia mission, our distance estimates are consistent with those resulting from the theoretical calibrations, albeit with relatively larger systematic uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Weak lensing trispectrum and Kurt-spectra

We introduce two kurt-spectra to probe fourth-order statistics of weak lensing convergence maps. Using state-of-the-art numerical simulations, we study the shapes of these kurt-spectra as a function of source redshifts and smoothing angular scales. We employ a pseudo-Cℓ approach to estimate the spectra from realistic convergence maps in the presence of an observational mask and noise for stage-IV large-scale structure surveys. We compare these results against theoretical predictions calculated using the FFTLog formalism, and find that a simple nonlinear clustering model — the hierarchical ansatz — can reproduce the numerical trends for the kurt-spectra in the nonlinear regime. In addition, we provide estimators for beyond fourth-order spectra where no definitive analytical results are available, and present corresponding results from numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating China's fossil-fuel CO2 emissions from a comprehensive dataset of nine inventories

China's fossil-fuel CO2 (FFCO2) emissions accounted for approximately 28 % of the global total FFCO2 in 2016. An accurate estimate of China's FFCO2 emissions is a prerequisite for global and regional carbon budget analyses and the monitoring of carbon emission reduction efforts. However, significant uncertainties and discrepancies exist in estimations of China's FFCO2 emissions due to a lack of detailed traceable emission factors (EFs) and multiple statistical data sources. Here, we evaluated China's FFCO2 emissions from nine published global and regional emission datasets. These datasets show that the total emissions increased from 3.4 (3.0–3.7) in 2000 to 9.8 (9.2–10.4) Gt CO2 per yr in 2016. The variations in these estimates were largely due to the different EF (0.491–0.746 t C per t of coal) and activity data. The large-scale patterns of gridded emissions showed a reasonable agreement, with high emissions being concentrated in major city clusters, and the standard deviation mostly ranged from 10 % to 40 % at the provincial level. However, patterns beyond the provincial scale varied significantly, with the top 5 % of the grid level accounting for 50 %–90 % of total emissions in these datasets. Our findings highlight the significance of using locally measured EF for Chinese coal. To reduce uncertainty, we recommend using physical CO2 measurements and use these values for dataset validation, key input data sharing (e.g., point sources), and finer-resolution validations at various levels.

Pengfei Han↗

Wide-angle effects for peculiar velocities

ABSTRACT The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as ‘wide-angle’ effects in redshift-space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work, we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in autocorrelation and in cross-correlation with the density field. Using both linear theory and the Zel'dovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low-redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration space, Fourier space, and spherical analyses; providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on dark matter from cosmic background anisotropies

The major stages in the linear evolution of the statistical ensemble of adiabatic fluctuations in radiation, baryons, and dark matter are discussed. If it is assumed that the distribution of light emitters (i.e., of galaxies) follows the distribution of mass (i.e., of dark matter), then universes dominated by massive collisionless relics of the Big Bang must have Omega larger than 0.2 h exp -4/3 to avoid exceeding the current observational limits on small-scale anisotropies in the microwave background. However, values of Omega of about 0.2 are indicated by dynamical studies of galaxy clustering. It is concluded that universes dominated by cold dark matter in which light traces mass are probably not viable models.

Bond, J. R.↗

Point mass deflectors in gravitational lenses

The observed properties of the six systems generally accepted as gravitational lenses are compared with the properties predicted for the images of point-mass and distributed mass deflectors. Taken as a statistical distribution, the characteristics of the six known systems strongly suggest the existence of a super-massive black hole as the effective deflecting mass in one or more of these systems. If the deflecting mass in a gravitational lens is a black hole, the distance to the deflecting mass can be determined directly from the difference in light travel time along the separate image paths. (No direct solution for the distance to the deflecting mass exists when the deflecting mass is a galaxy or a cluster of galaxies). Geometric parallaxes of objects at 1000 Mpc distance are of obvious importance in a wide variety of cosmological studies.

Dolan, Joseph F.↗

Anomaly detection in Hyper Suprime-Cam galaxy images with generative adversarial networks

ABSTRACT The problem of anomaly detection in astronomical surveys is becoming increasingly important as data sets grow in size. We present the results of an unsupervised anomaly detection method using a Wasserstein generative adversarial network (WGAN) on nearly one million optical galaxy images in the Hyper Suprime-Cam (HSC) survey. The WGAN learns to generate realistic HSC-like galaxies that follow the distribution of the data set; anomalous images are defined based on a poor reconstruction by the generator and outlying features learned by the discriminator. We find that the discriminator is more attuned to potentially interesting anomalies compared to the generator, and compared to a simpler autoencoder-based anomaly detection approach, so we use the discriminator-selected images to construct a high-anomaly sample of ∼13 000 objects. We propose a new approach to further characterize these anomalous images: we use a convolutional autoencoder to reduce the dimensionality of the residual differences between the real and WGAN-reconstructed images and perform UMAP clustering on these. We report detected anomalies of interest including galaxy mergers, tidal features, and extreme star-forming galaxies. A follow-up spectroscopic analysis of one of these anomalies is detailed in the Appendix; we find that it is an unusual system most likely to be a metal-poor dwarf galaxy with an extremely blue, higher-metallicity H ii region. We have released a catalogue with the WGAN anomaly scores; the code and catalogue are available at https://github.com/kstoreyf/anomalies-GAN-HSC; and our interactive visualization tool for exploring the clustered data is at https://weirdgalaxi.es.

79 ASTRONOMY AND ASTROPHYSICS↗

Using action space clustering to constrain the recent accretion history of Milky Way-like galaxies

ABSTRACT In the currently favoured cosmological paradigm galaxies form hierarchically through the accretion of satellites. Since a satellite is less massive than the host, its stars occupy a smaller volume in action space. Actions are conserved when the potential of the host halo changes adiabatically, so stars from an accreted satellite would remain clustered in action space as the host evolves. In this paper, we identify recently disrupted accreted satellites in three Milky Way-like disc galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellites through simulation snapshots. We try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the z = 0 snapshot. Even with completely error-free mock data we find that only 35 per cent (14/39) satellites are well recovered while the rest (25/39) are poorly recovered (i.e. either contaminated or split up). Most (10/14 ∼70 per cent) of the well-recovered satellites have infall times <7.1 Gyr ago and total mass >4 × 108M⊙ (stellar mass more than 1.2 × 106 M⊙, although our upper mass limit is likely to be resolution dependent). Since cosmological simulations predict that stellar haloes include a population of in situ stars, we test our ability to recover satellites when the data include 10–50 per cent in situ contamination. We find that most previously well-recovered satellites stay well recovered even with 50 per cent contamination. With the wealth of 6D phase space data becoming available we expect that cluster analysis in action space will be useful in identifying the majority of recently accreted and moderately massive satellites in the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Is Terzan 5 the remnant of a building block of the Galactic bulge? Evidence from APOGEE

ABSTRACT It has been proposed that the globular cluster-like system Terzan 5 is the surviving remnant of a primordial building block of the Milky Way bulge, mainly due to the age/metallicity spread and the distribution of its stars in the α–Fe plane. We employ Sloan Digital Sky Survey data from the Apache Point Observatory Galactic Evolution Experiment to test this hypothesis. Adopting a random sampling technique, we contrast the abundances of 10 elements in Terzan 5 stars with those of their bulge field counterparts with comparable atmospheric parameters, finding that they differ at statistically significant levels. Abundances between the two groups differ by more than 1σ in Ca, Mn, C, O, and Al, and more than 2σ in Si and Mg. Terzan 5 stars have lower [α/Fe] and higher [Mn/Fe] than their bulge counterparts. Given those differences, we conclude that Terzan 5 is not the remnant of a major building block of the bulge. We also estimate the stellar mass of the Terzan 5 progenitor based on predictions by the Evolution and Assembly of GaLaxies and their Environments suite of cosmological numerical simulations, concluding that it may have been as low as ∼3 × 108 M⊙ so that it was likely unable to significantly influence the mean chemistry of the bulge/inner disc, which is significantly more massive (∼1010 M⊙). We briefly discuss existing scenarios for the nature of Terzan 5 and propose an observational test that may help elucidate its origin.

79 ASTRONOMY AND ASTROPHYSICS↗