Search NASA⌕ Search

SEARCH · Search NASA

Results for “Astronomy image processing”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

TESS Hunt for Young and Maturing Exoplanets (THYME). II. A 17 Myr Old Transiting Hot Jupiter in the Sco-Cen Association

We present the discovery of a transiting hot Jupiter orbiting HIP 67522 (T {sub eff} ∼ 5650 K; M {sub *} ∼ 1.2M {sub ⊙}) in the 10–20 Myr old Sco-Cen OB association. We identified the transits in the TESS data using our custom notch filter planet search pipeline and characterize the system with additional photometry from Spitzer; spectroscopy from SOAR/Goodman, SALT/HRS, LCOGT/NRES, and SMARTS/CHIRON; and speckle imaging from SOAR/HRCam. We model the photometry as a periodic Gaussian process with transits to account for stellar variability and find an orbital period of 6.9596{sub −0.000015}{sup +0.000016} days and radius of 10.02{sub −0.53}{sup +0.54} R {sub ⊕}. We also identify a single transit of an additional candidate planet with radius 8.01{sub −0.71}{sup +0.75} R {sub ⊕} that has an orbital period of ≳23 days. The validated planet HIP 67522b is currently the youngest transiting hot Jupiter discovered and is an ideal candidate for transmission spectroscopy and radial velocity follow-up studies, while also demonstrating that some young giant planets either form in situ at small orbital radii or else migrate promptly from formation sites farther out in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

DeepCMB: Lensing reconstruction of the cosmic microwave background with deep neural networks

Next-generation cosmic microwave background (CMB) experiments will have lower noise and therefore increased sensitivity, enabling improved constraints on fundamental physics parameters such as the sum of neutrino masses and the tensor-to-scalar ratio r. Achieving competitive constraints on these parameters requires high signal-to-noise extraction of the projected gravitational potential from the CMB maps. Standard methods for reconstructing the lensing potential employ the quadratic estimator (QE). However, the QE performs suboptimally at the low noise levels expected in upcoming experiments. Other methods, like maximum likelihood estimators (MLE), are under active development. In this work, we demonstrate reconstruction of the CMB lensing potential with deep convolutional neural networks (CNN) - ie, a ResUNet. The network is trained and tested on simulated data, and otherwise has no physical parametrization related to the physical processes of the CMB and gravitational lensing. We show that, over a wide range of angular scales, ResUNets recover the input gravitational potential with a higher signal-to-noise ratio than the QE method, reaching levels comparable to analytic approximations of MLE methods. We demonstrate that the network outputs quantifiably different lensing maps when given input CMB maps generated with different cosmologies. We also show we can use the reconstructed lensing map for cosmological parameter estimation. This application of CNN provides a few innovations at the intersection of cosmology and machine learning. First, while training and regressing on images, we predict a continuous-variable field rather than discrete classes. Second, we are able to establish uncertainty measures for the network output that are analogous to standard methods. We expect this approach to excel in capturing hard-to-model non-Gaussian astrophysical foreground and noise contributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Investigating the Magnetic Imprints of Major Solar Eruptions with SDO /HMI High-cadence Vector Magnetograms

The solar active region photospheric magnetic field evolves rapidly during major eruptive events, suggesting appreciable feedback from the corona. Previous studies of these “magnetic imprints” are mostly based on line of sight only or lower-cadence vector observations; a temporally resolved depiction of the vector field evolution is hitherto lacking. Here, we introduce the high-cadence (90 s or 135 s) vector magnetogram data set from the Helioseismic and Magnetic Imager, which is well suited for investigating the phenomenon. These observations allow quantitative characterization of the permanent, step-like changes that are most pronounced in the horizontal field component (B {sub h}). A highly structured pattern emerges from analysis of an archetypical event, SOL2011-02-15T01:56, where B {sub h} near the main polarity inversion line increases significantly during the earlier phase of the associated flare with a timescale of several minutes, while B {sub h} in the periphery decreases at later times with smaller magnitudes and a slightly longer timescale. The data set also allows effective identification of the “magnetic transient” artifact, where enhanced flare emission alters the Stokes profiles and the inferred magnetic field becomes unreliable. Our results provide insights on the momentum processes in solar eruptions. The data set may also be useful to the study of sunquakes and data-driven modeling of the corona.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Long-lived Dust Asymmetries at Dead Zone Edges in Protoplanetary Disks

A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after 10{sup 4} orbits, or 2.5 Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

DeepCMB: Lensing Reconstruction of the CMB with Deep Neural Networks

Next-generation cosmic microwave background (CMB) experiments will have lower noise and therefore increased sensitivity, enabling improved constraints on fundamental physics parameters such as the sum of neutrino masses and the tensor-to-scalar ratio r. Achieving competitive constraints on these parameters requires high signal-to-noise extraction of the projected gravitational potential from the CMB maps. Standard methods for reconstructing the lensing potential employ the quadratic estimator (QE). However, the QE performs suboptimally at the low noise levels expected in upcoming experiments. Other methods, like maximum likelihood estimators (MLE), are under active development. In this work, we demonstrate reconstruction of the CMB lensing potential with deep convolutional neural networks (CNN) - ie, a ResUNet. The network is trained and tested on simulated data, and otherwise has no physical parametrization related to the physical processes of the CMB and gravitati onal lensing. We show that, over a wide range of angular scales, ResUNets recover the input gravitational potential with a higher signal-to-noise ratio than the QE method, reaching levels comparable to analytic approximations of MLE methods. We demonstrate that the network outputs quantifiably different lensing maps when given input CMB maps generated with different cosmologies. We also show we can use the reconstructed lensing map for cosmological parameter estimation. This application of CNN provides a few innovations at the intersection of cosmology and machine learning. First, while training and regressing on images, we predict a continuous-variable field rather than discrete classes. Second, we are able to establish uncertainty measures for the network output that are analogous to standard methods. We expect this approach to excel in capturing hard-to-model non-Gaussian astrophysical foreground and noise contributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

AN INVESTIGATION OF THE SOURCES OF EARTH-DIRECTED SOLAR WIND DURING CARRINGTON ROTATION 2053

In this work we analyze multiple sources of solar wind through a full Carrington Rotation (CR 2053) by analyzing the solar data through spectroscopic observations of the plasma upflow regions and the in situ data of the wind itself. Following earlier authors, we link solar and in situ observations by a combination of ballistic backmapping and potential-field source-surface modeling. We find three sources of fast solar wind that are low-latitude coronal holes. The coronal holes do not produce a steady fast wind, but rather a wind with rapid fluctuations. The coronal spectroscopic data from Hinode ’s Extreme Ultraviolet Imaging Spectrometer show a mixture of upflow and downflow regions highlighting the complexity of the coronal hole, with the upflows being dominant. There is a mix of open and multi-scale closed magnetic fields in this region whose (interchange) reconnections are consistent with the up- and downflows they generate being viewed through an optically thin corona, and with the strahl directions and freeze-in temperatures found in in situ data. At the boundary of slow and fast wind streams there are three short periods of enhanced-velocity solar wind, which we term intermediate based on their in situ characteristics. These are related to active regions that are located beside coronal holes. The active regions have different magnetic configurations, from bipolar through tripolar to quadrupolar, and we discuss the mechanisms to produce this intermediate wind, and the important role that the open field of coronal holes adjacent to closed-field active regions plays in the process.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Partially Transparent Petaled Mask/Occulter for Visible-Range Spectrum

The presence of the Poisson Spot, also known as the spot of Arago, has been known since the 18th century. This spot is the consequence of constructive interference of light diffracted by the edge of the obstacle where the central position can be determined by symmetry of the object. More recently, many NASA missions require the suppression of this spot in the visible range. For instance, the exoplanetary missions involving space telescopes require telescopes to image the planetary bodies orbiting central stars. For this purpose, the starlight needs to be suppressed by several orders of magnitude in order to image the reflected light from the orbiting planet. For the Earth-like planets, this suppression needs to be at least ten orders of magnitude. One of the common methods of suppression involves sharp binary petaled occulters envisioned to be placed many thousands of miles away from the telescope blocking the starlight. The suppression of the Poisson Spot by binary sharp petal tips can be problematic when the thickness of the tips becomes smaller than the wavelength of the incident beam. First they are difficult to manufacture and also it invalidates the laws of physical optics. The proposed partially transparent petaled masks/occulters compensate for this sharpness with transparency along the surface of the petals. Depending on the geometry of the problem, this transparency can be customized such that only a small region of the petal is transparent and the remaining of the surface is opaque. This feature allows easy fabrication of this type of occultation device either as a mask or occulter. A partially transparent petaled mask/ occulter has been designed for the visible spectrum range. The mask/occulter can suppress the intensity along the optical axis up to ten orders of magnitude. The design process can tailor the mask shape, number of petals, and transparency level to the near-field and farfield diffraction region. The mask/occulter can be used in space astronomy, ground-based telescope, and high-energy laser systems, and optical lithography to eliminate the Poisson Spot.

Shiri, Ron Shahram↗
DOE OSTI2017

Long-term Study of the Solar Filaments from the Synoptic Maps as Derived from H{sub α} Spectroheliograms of the Kodaikanal Observatory

The century long (1914–2007) H{sub α} (656.28 nm) spectroheliograms from the Kodaikanal Solar Observatory (KSO) have been recently digitized. Using these newly calibrated, processed images we study the evolution of dark elongated on-disk structures called filaments, which are potential representatives of magnetic activities on the Sun. To our knowledge, this is the oldest uniform digitized data set with daily images available today in H{sub α}. We generate Carrington maps for the entire time duration and try to find the correlations with maps of the same Carrington rotation from the Ca ii K KSO data. Filaments are segmented from the Carrington maps using a semi-automated technique and are studied individually to extract their centroids and tilts. We plot the time-latitude distribution of the filament centroids, producing a butterfly diagram which clearly shows the presence of poleward migration. We separate polar filaments for each cycle and try to estimate the delay between the polar filament number cycle and the sunspot number cycle peaks. We correlate this delay with the delay between polar reversal and sunspot number maxima. This provides new insight on the role of polar filaments on polar reversal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Properties of tree rings in LSST sensors

Images of uniformly illuminated sensors for the Large Synoptic Survey Telescope have circular periodic patterns with an appearance similar to tree rings. Furthermore, these patterns are caused by circularly symmetric variations of the dopant concentration in the monocrystal silicon boule induced by the manufacturing process. Non-uniform charge density results in the parasitic electric field inside the silicon sensor, which may distort shapes of astronomical sources. Here, we analyzed data from fifteen LSST sensors produced by ITL to determine the main parameters of the tree rings: amplitude and period, and also variability across the sensors tested at Brookhaven National Laboratory. Tree ring pattern has a weak dependence on the wavelength. But the ring amplitude gets smaller as wavelength gets longer, since longer wavelengths penetrate deeper into the silicon. Tree ring amplitude gets larger as it gets closer to the outer part of the wafer, from 0.1 to 1.0%, indicating that the resistivity variation is larger for larger radii.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Lensing substructure quantification in RXJ1131-1231: a 2 keV lower bound on dark matter thermal relic mass

We study the substructure content of the strong gravitational lens RXJ1131-1231 through a forward modelling approach that relies on generating an extensive suite of realistic simulations. We use a semi-analytic merger tree prescription that allows us to stochastically generate substructure populations whose properties depend on the dark matter particle mass. These synthetic halos are then used as lenses to produce realistic mock images that have the same features, e.g. luminous arcs, quasar positions, instrumental noise and PSF, as the data. We then analyse the data and the simulations in the same way with summary statistics that are sensitive to the signal being targeted and are able to constrain models of dark matter statistically using Approximate Bayesian Computing (ABC) techniques. (In this work, we focus on the thermal relic mass estimate and fix the semi-analytic descriptions of the substructure evolution based on recent literature.) We are able, based on the HST data for RXJ1131-1231, to rule out a warm dark matter thermal relic mass below 2 keV at the 2σ confidence level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Redshift-space distortions in simulations of the 21-cm signal from the cosmic dawn

ABSTRACT The 21-cm signal from the Cosmic Dawn (CD) is likely to contain large fluctuations, with the most extreme astrophysical models on the verge of being ruled out by observations from radio interferometers. It is therefore vital that we understand not only the astrophysical processes governing this signal, but also other inherent processes impacting the signal itself, and in particular line-of-sight effects. Using our suite of fully numerical radiative transfer simulations, we investigate the impact on the redshifted 21-cm from the CD from one of these processes, namely the redshift-space distortions (RSDs). When RSDs are added, the resulting boost to the power spectra makes the signal more or equally detectable for our models for all redshifts, further strengthening hopes that a power spectra measurement of the CD will be possible. RSDs lead to anisotropy in the signal at the beginning and end of the CD, but not while X-ray heating is underway. The inclusion of RSDs, however, decreases detectability of the non-Gaussianity of fluctuations from inhomogeneous X-ray heating as measured by the skewness and kurtosis. On the other hand, mock observations created from all our simulations that include telescope noise corresponding to 1000 h of observation with the Square Kilometre Array telescope show that we may be able to image the CD for all heating models considered and suggest RSDs dramatically boost fluctuations coming from the inhomogeneous Ly α background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ODYSSEY: A PUBLIC GPU-BASED CODE FOR GENERAL RELATIVISTIC RADIATIVE TRANSFER IN KERR SPACETIME

General relativistic radiative transfer calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra, and light curves from matter in the vicinity of black holes. Such studies are especially important for ongoing and upcoming millimeter/submillimeter very long baseline interferometry observations of the supermassive black holes at the centers of Sgr A* and M87. To this end we introduce Odyssey, a graphics processing unit (GPU) based code for ray tracing and radiative transfer in the Kerr spacetime. On a single GPU, the performance of Odyssey can exceed 1 ns per photon, per Runge–Kutta integration step. Odyssey is publicly available, fast, accurate, and flexible enough to be modified to suit the specific needs of new users. Along with a Graphical User Interface powered by a video-accelerated display architecture, we also present an educational software tool, Odyssey-Edu, for showing in real time how null geodesics around a Kerr black hole vary as a function of black hole spin and angle of incidence onto the black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

The Statistical Consulting Center for Astronomy (SCCA)

The process by which raw astronomical data acquisition is transformed into scientifically meaningful results and interpretation typically involves many statistical steps. Traditional astronomy limits itself to a narrow range of old and familiar statistical methods: means and standard deviations; least-squares methods like chi(sup 2) minimization; and simple nonparametric procedures such as the Kolmogorov-Smirnov tests. These tools are often inadequate for the complex problems and datasets under investigations, and recent years have witnessed an increased usage of maximum-likelihood, survival analysis, multivariate analysis, wavelet and advanced time-series methods. The Statistical Consulting Center for Astronomy (SCCA) assisted astronomers with the use of sophisticated tools, and to match these tools with specific problems. The SCCA operated with two professors of statistics and a professor of astronomy working together. Questions were received by e-mail, and were discussed in detail with the questioner. Summaries of those questions and answers leading to new approaches were posted on the Web (www.state.psu.edu/ mga/SCCA). In addition to serving individual astronomers, the SCCA established a Web site for general use that provides hypertext links to selected on-line public-domain statistical software and services. The StatCodes site (www.astro.psu.edu/statcodes) provides over 200 links in the areas of: Bayesian statistics; censored and truncated data; correlation and regression, density estimation and smoothing, general statistics packages and information; image analysis; interactive Web tools; multivariate analysis; multivariate clustering and classification; nonparametric analysis; software written by astronomers; spatial statistics; statistical distributions; time series analysis; and visualization tools. StatCodes has received a remarkable high and constant hit rate of 250 hits/week (over 10,000/year) since its inception in mid-1997. It is of interest to scientists both within and outside of astronomy. The most popular sections are multivariate techniques, image analysis, and time series analysis. Hundreds of copies of the ASURV, SLOPES and CENS-TAU codes developed by SCCA scientists were also downloaded from the StatCodes site. In addition to formal SCCA duties, SCCA scientists continued a variety of related activities in astrostatistics, including refereeing of statistically oriented papers submitted to the Astrophysical Journal, talks in meetings including Feigelson's talk to science journalists entitled "The reemergence of astrostatistics" at the American Association for the Advancement of Science meeting, and published papers of astrostatistical content.

Akritas, Michael↗
DOE OSTI2016

GLOBAL ENERGETICS OF SOLAR FLARES. IV. CORONAL MASS EJECTION ENERGETICS

This study entails the fourth part of a global flare energetics project, in which the mass m {sub cme}, kinetic energy E {sub kin}, and the gravitational potential energy E {sub grav} of coronal mass ejections (CMEs) is measured in 399 M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission, using a new method based on the EUV dimming effect. EUV dimming is modeled in terms of a radial adiabatic expansion process, which is fitted to the observed evolution of the total emission measure of the CME source region. The model derives the evolution of the mean electron density, the emission measure, the bulk plasma expansion velocity, the mass, and the energy in the CME source region. The EUV dimming method is truly complementary to the Thomson scattering method in white light, which probes the CME evolution in the heliosphere at r ≳ 2 R {sub ⊙}, while the EUV dimming method tracks the CME launch in the corona. We compare the CME parameters obtained in white light with the LASCO/C2 coronagraph with those obtained from EUV dimming with the Atmospheric Imaging Assembly onboard the SDO for all identical events in both data sets. We investigate correlations between CME parameters, the relative timing with flare parameters, frequency occurrence distributions, and the energy partition between magnetic, thermal, nonthermal, and CME energies. CME energies are found to be systematically lower than the dissipated magnetic energies, which is consistent with a magnetic origin of CMEs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Evolution of Occator Crater on (1) Ceres

The dwarf planet Ceres (diameter 939 km) is the largest object in the main asteroid belt. Recent investigations suggest that Ceres is a thermally evolved, volatile-rich body with potential geological activity, a body which was never completely molten but possibly differentiated into a rocky core, an ice-rich mantle, and which may contain remnant internal liquid water. Thermal alteration and exogenic material infall contribute to producing a (dark) carbonaceous chondritic-like surface containing ammoniated phyllosilicates. Here we report imaging and spectroscopic analyses of Occator crater derived from the Framing Camera and the Visible and Infrared Spectrometer onboard Dawn. We found that the central bright spot (Cerealia Facula) of Occator is ∼30 Myr younger than the crater itself. The central spot is located in a central pit which contains a dome that is spectrally homogenous, exhibiting absorption features that are consistent with carbonates. Multiple radial fractures across the dome indicate an extrusive formation process. Our results lead us to conclude that the floor region was subject to past endogenic activity. Dome and bright material in its vicinity formed likely due to a long-lasting, periodic, or episodic ascent of bright material from a subsurface reservoir rich in carbonates. Originally triggered by an impact event, gases, possibly dissolved from a subsurface water/brine layer, enabled material rich in carbonates to ascend through fractures and be deposited onto the surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Formation and Maintenance of the Dominant Southern Polar Crown Cavity of Cycle 24

In this article, we report a study of the longest-lived polar crown cavity of Solar Cycle 24, using an observation from 2013, and propose a physical mechanism to explain its sustained existence. We used high temporal and spatial resolution observations from the Atmospheric Imaging Assembly (AIA) and the Helioseismic Magnetic Imager (HMI) instruments on board the Solar Dynamics Observatory ( SDO ) to explore the structure and evolution of the cavity. Although it existed for more than a year, we examined the circumpolar cavity in great detail from 2013 March 21 to 2013 October 31. Our study reinforces the existing theory of formation of polar crown filaments that involves two basic processes to form any polar crown cavity as well as the long-lived cavity that we studied here. First, the underlying polarity inversion line (PIL) of the circumpolar cavity is formed between (1) the trailing part of dozens of decayed active regions distributed in different longitudes and (2) the unipolar magnetic field in the polar coronal hole. Second, the long life of the cavity is sustained by the continuing flux cancellation along the PIL. The flux is persistently transported toward the polar region through surface meridional flow and diffusion. The continuing flux cancellation leads to the shrinking of the polar coronal hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Solar Energetic Particle Event of 2010 August 14: Connectivity with the Solar Source Inferred from Multiple Spacecraft Observations and Modeling

We analyze one of the first solar energetic particle (SEP) events of solar cycle 24 observed at widely separated spacecraft in order to assess the reliability of models currently used to determine the connectivity between the sources of SEPs at the Sun and spacecraft in the inner heliosphere. This SEP event was observed on 2010 August 14 by near-Earth spacecraft, STEREO-A (∼80° west of Earth) and STEREO-B (∼72° east of Earth). In contrast to near-Earth spacecraft, the footpoints of the nominal magnetic field lines connecting STEREO-A and STEREO-B with the Sun were separated from the region where the parent fast halo coronal mass ejection (CME) originated by ∼88° and ∼47° in longitude, respectively. We discuss the properties of the phenomena associated with this solar eruption. Extreme ultraviolet and white-light images are used to specify the extent of the associated CME-driven coronal shock. We then assess whether the SEPs observed at the three heliospheric locations were accelerated by this shock or whether transport mechanisms in the corona and/or interplanetary space provide an alternative explanation for the arrival of particles at the poorly connected spacecraft. A possible scenario consistent with the observations indicates that the observation of SEPs at STEREO-B and near Earth resulted from particle injection by the CME shock onto the field lines connecting to these spacecraft, whereas SEPs reached STEREO-A mostly via cross-field diffusive transport processes. The successes, limitations, and uncertainties of the methods used to resolve the connection between the acceleration sites of SEPs and the spacecraft are evaluated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

Abstract The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm , ipole , and igrmonty , paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

79 ASTRONOMY AND ASTROPHYSICS↗