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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The X-Ray Core of the Low-Luminosity Radio Galaxy 3C346 and ASCA Spectroscopy to Test BL LAC/Radio Galaxy Unification

Radio galaxies are relatively faint sources for Advanced Spacecraft for Cosmology Astrophysics (ASCA), and so in order to get the best possible results from the observations two things have been necessary, both of which delayed the fast preparation of papers. Firstly, the best possible data screening and background subtraction were necessary to improve the signal-to-noise, and all our several initial analysis trials were discarded in favor of using FTOOLS versions 4.1 and above. Secondly, we found that the ASCA spectra were statistically too poor to discriminate well between non-thermal and thermal models, never mind the mixture of the two which we expected on the basis of our ROSAT spatial separation of components in radio galaxies. This means that in each case we have needed to combine the ASCA spectroscopy with analysis of data from other X-ray or radio observations in order to exploit the ASCA data to the full. Our analysis for 3C 346 has yielded the cleanest final result. This powerful radio galaxy at a redshift of 0.161, lies in a poor cluster, which we have separated well from the dominant X-ray component of unresolved emission using a spatial analysis of archival ROSAT data. We were then able to fix the thermal component in our ASCA spectral analysis, and have found evidence that the unresolved emission varied by 32 +/- 13% over the 18 months between the ROSAT and ASCA observations. The unresolved X-ray emission does not suffer from intrinsic absorption, and we have related it to radio structures on both milliarcsecond scales and the arcsecond scales which Chandra can resolve. The source is a target of a Chandra AO2 proposal which we have recently submitted to follow up on our ASCA (and ROSAT) work. 3C 346's orientation to the line of sight is uncertain. However, the absence of X-ray absorption, and the radio/optical/X-ray colors, when combined with with previous radio evidence that the source is a foreshortened radio galaxy of the FRII class, suggest that the radio jets are seen at an angle to the line of sight of about 30 deg, intermediate between the radio-galaxy and quasar classes. The relatively hard ASCA response has allowed us to place an upper limit of 5.6 x 10(exp 43) ergs/ s on the 2-10 keV luminosity of any central X-ray component absorbed bN, gas which might be obscuring the broad-line emission region. Attached to this report is an almost final draft of a paper which we have prepared for submission to the Astrophysical Journal. Our combined ASCA and ROSAT results for NGC 6251 rule out our previously preferred flat-spectrum model and inverse-Compton interpretation for the source based on ROSAT data alone, but a softer X-ray spectrum and moderate absorption bring all the available data (including our early VLA HI measurements) into consistency, and we are reasonably confident that we understand the processes responsible for the X-ray emission. We have made some more sensitive HI absorption measurements which are currently being analyzed, and our plans are to publish our ASCA analysis in conjunction with the new HI results. The ASCA data for NGC 4261 have been difficult to interpret. A re-analysis of our ROSAT data with a wider range of physical parameters brings the ROSAT and ASCA results into reasonable agreement only if the emission from hot gas dominates more than suggested by our earlier work, which is itself unexpected since the radio core is bright and a large jet-related X-ray component would bring the source into agreement with results for others of its type. However, we have recently received our Chandra A01 data for this source, with the spatial resolution which allows us to separate thermal and non-thermal emission components. Our ASCA results will be re-interpreted once the analysis of our Chandra data is complete. The interpretation of the ASCA data for BL Lac object 3C 371 is ongoing, in conjunction with analysis of archival multifrequency data. Radio galaxies are complex in their X-ray properties, and hindsight has shown that the spatial resolution of ASCA is too poor for a reliable interpretation of the data without drawing on other observations. However, the ASCA spectra have made a useful contribution to the interpretation of these sources, and the groundwork is now there for more sensitive work using Chandra and XMM-Newton.

Worrall, Diana↗
NASA NTRS2000

Satellite Data Sets in the Polar Regions

We have generated about two decades of consistently derived geophysical parameters in the polar regions. The key parameters are sea ice concentration, surface temperature, albedo, and cloud cover statistics. Sea ice concentrations were derived from the Scanning Multichannel Microwave Radiometer (SMMR) data and the Special Scanning Cl Microwave Imager (SSM/I) data from several platforms using the enhanced Bootstrap Algorithm for the period 1978 through 1999. The new algorithm reduces the errors associated with spatial and temporal variations in the emissivity and surface temperatures of sea ice. Also, bad data at ocean/land interfaces are identified and deleted in an unsupervised manner. Surface ice temperature, albedo and cloud cover statistics are derived simultaneously from the Advanced Very High Resolution Radiometer (AVHRR) data from 1981 through 1999 and mapped at a higher resolution but the same format as the ice concentration data. The technique makes use these co-registered ice concentration maps to enable cloud masking to be done separately for open ocean, sea ice and land areas. The effect of inversion is minimized by taking into consideration the expected changes in the effect of inversion with altitude, especially in the Antarctic. A technique for ice type regional classification has also been developed using multichannel cluster analysis and a neural network. This provide a means to identify large areas of thin ice, first year ice, and older ice types. The data sets have been shown to be coherent with each other and provide a powerful tool for in depth studies of the currently changing Arctic and Antarctic environment.

Comiso, Josefino C.↗
NASA NTRS1994

Statistical scaling relationships and size effects in the strength and creep rupture of fibrous composites

In this presentation we discuss a new theoretical model and supporting experimental results for the strength and lifetime in creep rupture of unidirectional, carbon fiber/epoxy matrix composites at ambient conditions. First we review the 'standard' Weibull/power-law methodology that has been standard practice. Then we discuss features of a recent model which build on the statistical aspects of fiber strength, micromechanical aspects of stress transfer around fiber breaks, and time-dependent creep of the matrix. The model is applied to 'microcomposites' consisting of seven fibers in a matrix for which strength and creep-rupture data are available. The model yields Weibull distributions in an envelope format for both strength and lifetime. The respective shape, scale and power-law parameters depend on such parameters as the Weibull shape parameter for fiber strength, the exponent for matrix creep, the effective load transfer length (which grows in time due to matrix creep) and the critical cluster size for failed fibers. The experimental results are consistent with the theory, though time-dependent debonding appears to be part of the failure process.

Phoenix, S. Leigh↗
NASA NTRS2015

Real-Time Kennedy Space Center and Cape Canaveral Air Force Station High-Resolution Model Implementation and Verification

Customer: NASA's Launch Services Program (LSP), Ground Systems Development and Operations (GSDO), and Space Launch System (SLS) programs. NASA's LSP, GSDO, SLS and other programs at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) use the daily and weekly weather forecasts issued by the 45th Weather Squadron (45 WS) as decision tools for their day-to-day and launch operations on the Eastern Range (ER). For example, to determine if they need to limit activities such as vehicle transport to the launch pad, protect people, structures or exposed launch vehicles given a threat of severe weather, or reschedule other critical operations. The 45 WS uses numerical weather prediction models as a guide for these weather forecasts, particularly the Air Force Weather Agency (AFWA) 1.67 kilometer Weather Research and Forecasting (WRF) model. Considering the 45 WS forecasters' and Launch Weather Officers' (LWO) extensive use of the AFWA model, the 45 WS proposed a task at the September 2013 Applied Meteorology Unit (AMU) Tasking Meeting requesting the AMU verify this model. Due to the lack of archived model data available from AFWA, verification is not yet possible. Instead, the AMU proposed to implement and verify the performance of an ER version of the AMU high-resolution WRF Environmental Modeling System (EMS) model (Watson 2013) in real-time. The tasking group agreed to this proposal; therefore the AMU implemented the WRF-EMS model on the second of two NASA AMU modeling clusters. The model was set up with a triple-nested grid configuration over KSC/CCAFS based on previous AMU work (Watson 2013). The outer domain (D01) has 12-kilometer grid spacing, the middle domain (D02) has 4-kilometer grid spacing, and the inner domain (D03) has 1.33-kilometer grid spacing. The model runs a 12-hour forecast every hour, D01 and D02 domain outputs are available once an hour and D03 is every 15 minutes during the forecast period. The AMU assessed the WRF-EMS 1.33-kilometer domain model performance for the 2014 warm season (May-September). Verification statistics were computed using the Model Evaluation Tools, which compared the model forecasts to observations. The mean error values were close to 0 and the root mean square error values were less than 1.8 for mean sea-level pressure (millibars), temperature (degrees Kelvin), dewpoint temperature (degrees Kelvin), and wind speed (per millisecond), all very small differences between the forecast and observations considering the normal magnitudes of the parameters. The precipitation forecast verification results showed consistent under-forecasting of the precipitation object size. This could be an artifact of calculating the statistics for each hour rather than for the entire 12-hour period. The AMU will continue to generate verification statistics for the 1.33-kilometer WRF-EMS domain as data become available in future cool and warm seasons. More data will produce more robust statistics and reveal a more accurate assessment of model performance. Once the formal task was complete, the AMU conducted additional work to better understand the wind direction results. The results were stratified diurnally and by wind speed to determine what effects the stratifications would have on the model wind direction verification statistics. The results are summarized in the addendum at the end of this report. In addition to verifying the model's performance, the AMU also made the output available in the Advanced Weather Interactive Processing System II (AWIPS II). This allows the 45 WS and AMU staff to customize the model output display on the AMU and Range Weather Operations AWIPS II client computers and conduct real-time subjective analyses. In the future, the AMU will implement an updated version of the WRF-EMS model that incorporates local data assimilation. This model will also run in real-time and be made available in AWIPS II.

Numerical Weather Prediction↗
NASA NTRS1993

Quons, an interpolation between Bose and Fermi oscillators

After a brief mention of Bose and Fermi oscillators and of particles which obey other types of statistics, including intermediate statistics, parastatistics, paronic statistics, anyon statistics, and infinite statistics, I discuss the statistics of 'quons' (pronounced to rhyme with muons), particles whose annihilation and creation operators obey the q-deformed commutation relation (the quon algebra or q-mutator) which interpolates between fermions and bosons. I emphasize that the operator for interaction with an external source must be an effective Bose operator in all cases. To accomplish this for parabose, parafermi and quon operators, I introduce parabose, parafermi, and quon Grassmann numbers, respectively. I also discuss interactions of non-relativistic quons, quantization of quon fields with antiparticles, calculation of vacuum matrix elements of relativistic quon fields, demonstration of the TCP theorem, cluster decomposition, and Wick's theorem for relativistic quon fields, and the failure of local commutativity of observables for relativistic quon fields. I conclude with the bound on the parameter q for electrons due to the Ramberg-Snow experiment.

Greenberg, O. W.↗
NASA NTRS2011

The Transitional Protoplanetary Disk Frequency as a Function of Age: Disk Evolution in the Coronet Cluster, Taurus, and Other 1--8 Myr-old Regions

We present Spitzer 3.6-24 micron photometry and spectroscopy for stars in the 1-3 Myr-old Coronet Cluster, expanding upon the survey of Sicilia-Aguilar et al. (2008). Using sophisticated radiative transfer models, we analyze these new data and those from Sicilia-Aguilar et al. (2008) to identify disks with evidence for substantial dust evolution consistent with disk clearing: transitional disks. We then analyze data in Taurus and others young clusters - IC 348, NGC 2362, and eta Cha -- to constrain the transitional disk frequency as a function of time. Our analysis confirms previous results finding evidence for two types of transitional disks -- those with inner holes and those that are homologously depleted. The percentage of disks in the transitional phase increases from approx.15-20% at 1-2 Myr to > 50% at 5-8 Myr; the mean transitional disk lifetime is closer to approx. 1 Myr than 0.1-0.5 Myr, consistent with previous studies by Currie et al. (2009) and Sicilia-Aguilar et al. (2009). In the Coronet Cluster and IC 348, transitional disks are more numerous for very low-mass M3--M6 stars than for more massive K5-M2 stars, while Taurus lacks a strong spectral type-dependent frequency. Assuming standard values for the gas-to-dust ratio and other disk properties, the lower limit for the masses of optically-thick primordial disks is Mdisk approx. 0.001-0.003 M*. We find that single color-color diagrams do not by themselves uniquely identify transitional disks or primordial disks. Full SED modeling is required to accurately assess disk evolution for individual sources and inform statistical estimates of the transitional disk population in large samples using mid-IR colors.

Currie, Thayne↗
NASA NTRS2003

Transverse Dimensions of Chorus in the Source Region

We report measurement of whistler-mode chorus by the four Cluster spacecraft at close separations. We focus our analysis on the generation region close to the magnetic equatorial plane at a radial distance of 4.4 Earth's radii. We use both linear and rank correlation analysis to define perpendicular dimensions of the sources of chorus elements below one half of the electron cyclotron frequency. Correlation is significant throughout the range of separation distances of 60-260 km parallel to the field line and 7-100 km in the perpendicular plane. At these scales, the correlation coefficient is independent for parallel separations, and decreases with perpendicular separation. The observations are consistent with a statistical model of the source region assuming individual sources as gaussian peaks of radiated power with a common half-width of 35 km perpendicular to the magnetic field. This characteristic scale is comparable to the wavelength of observed waves.

Santolik, O.↗
NASA NTRS1991

Geography of cretaceous extinctions: Data base development

Data bases built from the source literature are plagued by problems of data quality. Unless the data acquisition is done by experts, working slowly, the data base may contain so much garbage that true signals and patterns cannot be detected. On the other hand, high quality data bases develop so slowly that satisfactory statistical analysis may never be possible due to the small sample sizes. Results of a test are presented of the opposite strategy: rapid data acquisition by non-experts with minimal control on data quality. A published list of 186 species and genera of fossil invertibrates of the latest Cretaceous Age (Maestrichtian) were located through a random search of the paleobiological and geological literature. The geographic location for each faunal list was then transformed electronically to Maestrichtian latitude and longitude and the lists were further digested to identify the genera occurring in each ten-degree, latitude-longitude block. The geographical lists were clustered using the Otsuka similarity coefficient and a standard unweight-pair-group method. The resulting clusters are remarkably consistent geographically, indicating that a strong biogeographic signal is visible despite low-quality data. A further test evaluated the geographic pattern of end-Cretaceaous extinctions. All genera in the data base were compared with Sepkoski's compendium of time ranges of genera to determine which of the reported genera survived the Cretaceous mass extinction. In turn, extinction rates for the ten-degree, latitude-longitude blocks were mapped. The resulting distribution is readily interpretable as a robust pattern of the geography of the mass extinction. The study demonstrates that a low-quality data base, built rapidly, can provide a basis for meaningful analysis of past biotic events.

Raup, D. M.↗
NASA NTRS2004

Chemical models for simulating single-walled nanotube production in arc vaporization and laser ablation processes

Chemical kinetic models for the nucleation and growth of clusters and single-walled carbon nanotube (SWNT) growth are developed for numerical simulations of the production of SWNTs. Two models that involve evaporation and condensation of carbon and metal catalysts, a full model involving all carbon clusters up to C80, and a reduced model are discussed. The full model is based on a fullerene model, but nickel and carbon/nickel cluster reactions are added to form SWNTs from soot and fullerenes. The full model has a large number of species--so large that to incorporate them into a flow field computation for simulating laser ablation and arc processes requires that they be simplified. The model is reduced by defining large clusters that represent many various sized clusters. Comparisons are given between these models for cases that may be applicable to arc and laser ablation production. Solutions to the system of chemical rate equations of these models for a ramped temperature profile show that production of various species, including SWNTs, agree to within about 50% for a fast ramp, and within 10% for a slower temperature decay time.

Nanotubes, Carbon/chemistry↗
NASA NTRS2016

A Neural Network Approach for Identifying Particle Pitch Angle Distributions in Van Allen Probes Data

Analysis of particle pitch angle distributions (PADs) has been used as a means to comprehend a multitude of different physical mechanisms that lead to flux variations in the Van Allen belts and also to particle precipitation into the upper atmosphere. In this work we developed a neural network-based data clustering methodology that automatically identifies distinct PAD types in an unsupervised way using particle flux data. One can promptly identify and locate three well-known PAD types in both time and radial distance, namely, 90deg peaked, butterfly, and flattop distributions. In order to illustrate the applicability of our methodology, we used relativistic electron flux data from the whole month of November 2014, acquired from the Relativistic Electron-Proton Telescope instrument on board the Van Allen Probes, but it is emphasized that our approach can also be used with multiplatform spacecraft data. Our PAD classification results are in reasonably good agreement with those obtained by standard statistical fitting algorithms. The proposed methodology has a potential use for Van Allen belt's monitoring.

Souza, V. M.↗
NASA NTRS1976

Shuttle onboard IMU alignment methods

The current approach to the shuttle IMU alignment is based solely on the Apollo Deterministic Method. This method is simple, fast, reliable and provides an accurate estimate for the present cluster to mean of 1,950 transformation matrix. If four or more star sightings are available, the application of least squares analysis can be utilized. The least squares method offers the next level of sophistication to the IMU alignment solution. The least squares method studied shows that a more accurate estimate for the misalignment angles is computed, and the IMU drift rates are a free by-product of the analysis. Core storage requirements are considerably more; estimated 20 to 30 times the core required for the Apollo Deterministic Method. The least squares method offers an intermediate solution utilizing as much data that is available without a complete statistical analysis as in Kalman filtering.

Henderson, D. M.↗
NASA NTRS1986

Statistical methods for astronomical data with upper limits. II - Correlation and regression

Statistical methods for calculating correlations and regressions in bivariate censored data where the dependent variable can have upper or lower limits are presented. Cox's regression and the generalization of Kendall's rank correlation coefficient provide significant levels of correlations, and the EM algorithm, under the assumption of normally distributed errors, and its nonparametric analog using the Kaplan-Meier estimator, give estimates for the slope of a regression line. Monte Carlo simulations demonstrate that survival analysis is reliable in determining correlations between luminosities at different bands. Survival analysis is applied to CO emission in infrared galaxies, X-ray emission in radio galaxies, H-alpha emission in cooling cluster cores, and radio emission in Seyfert galaxies.

Isobe, T.↗
NASA NTRS1988

IRAS observations of galaxies in the Virgo cluster area

IRAS data on 196 galaxies in the area of the Virgo Cluster are presented. The data derive from combining all available survey observations for each object and therefore achieve greater sensitivity than the IRAS Point Source Catalog (PSC). The enhanced sensitivity allows 78 galaxies to be detected at 12 microns, 82 at 25 microns, 139 at 60 microns, and 135 at 100 microns, compared to 16, 23, 88, and 95 detections listed in the PSC. From the blue compact dwarf galaxy sample, 23 and 19 objects are detected at 60 and 100 microns, compared to three and two detections listed in the PSC. The emission in three close pairs of galaxies which are reported as single sources in the PSC are separated here. These statistics demonstrate the importance and potential of a detailed examination of IRAS data, especially for possibly resolved sources and, in particular, for galaxies out to redshifts of 0.008 or galaxies with D(25) of 3 arcmin or greater.

Helou, George↗
NASA NTRS2023

Multiyear Dry Periods in Southern Africa

Characteristics and physical features related to low precipitation across many years in Southern Africa that lead to societal disruptions are diagnosed using observed analyses and an ensemble of historical coupled climate model simulations during 1921 to 2014. Four regions are evaluated, as identified through a hierarchical clustering algorithm applied to the Standardized Precipitation Index (SPI) during the October–April precipitation season. Although dryness spanning many October–April occurs periodically in each region, they seldom occur simultaneously, consistent with largely insignificant SPI cross-correlations between them. However, characteristics relevant to low precipitation across many years are generalizable between the four regions, including the serial persistence of October–April precipitation, the likelihood of consecutive dry October–April, and the likelihood of dry October–April in temporal extents of up to 10 consecutive such 7-month seasons. Systematic precipitation persistence is not a feature in any of the four Southern Africa regions, as serial correlations of October–April SPI are not statistically significant at any time lags. It follows that there is an exponential-folding decay in the likelihood of consecutive October–April for various SPI thresholds and that there is a large spread in the likelihood of low October–April SPI across many years. In terms of physical features, low October–April SPI in each Southern Africa region is closely related to local atmospheric circulations; however, they are not as closely related to sea surface temperatures (SSTs). These results suggest that dryness spanning many years is determined primarily by persistent local circulations related to atmospheric variability and to a lesser extent variability related to SST anomalies, including the El Niño–Southern Oscillation.

subtropical Indian Ocean dipole↗
NASA NTRS2006

Automated Classification of Thermal Infrared Spectra Using Self-organizing Maps

Existing and planned space missions to a variety of planetary and satellite surfaces produce an ever increasing volume of spectral data. Understanding the scientific informational content in this large data volume is a daunting task. Fortunately various statistical approaches are available to assess such data sets. Here we discuss an automated classification scheme based on Kohonen Self-organizing maps (SOM) we have developed. The SUM process produces an output layer were spectra having similar properties lie in close proximity to each other. One major effort is partitioning this output layer into appropriate regions. This is prefonned by defining dosed regions based upon the strength of the boundaries between adjacent cells in the SOM output layer. We use the Davies-Bouldin index as a measure of the inter-class similarities and intra-class dissimilarities that determines the optimum partition of the output layer, and hence number of SOM clusters. This allows us to identify the natural number of clusters formed from the spectral data. Mineral spectral libraries prepared at Arizona State University (ASU) and John Hopkins University (JHU) are used to test and evaluate the classification scheme. We label the library sample spectra in a hierarchical scheme with class, subclass, and mineral group names. We use a portion of the spectra to train the SOM, i.e. produce the output layer, while the remaining spectra are used to test the SOM. The test spectra are presented to the SOM output layer and assigned membership to the appropriate cluster. We then evaluate these assignments to assess the scientific meaning and accuracy of the derived SOM classes as they relate to the labels. We demonstrate that unsupervised classification by SOMs can be a useful component in autonomous systems designed to identify mineral species from reflectance and emissivity spectra in the therrnal IR.

Roush, Ted L.↗
NASA NTRS1990

The nova rate in the elliptical component of NGC 5128

The results are presented of a five year nova survey of NGC 5128 performed with an H-alpha filter at the prime focus of the CTIO 4 m telescope. In all, 16 novae were detected, of which 12 are members of a statistically complete and homogeneous sample unaffected by the absorption lane of the galaxy. When normalized to B luminosity, the nova rate in the elliptical component of NGC 5128 is 2.7 + or - 1.0 times less than that observed in the bulge of M31. However, when normalized in the infrared K bandpass, the two nova rates are identical, and comparable to that measured in the Large Magellanic Cloud. These data suggest that (1) the percentage of binary stars is approximately constant over a large range of star-forming environments, (2) the nova rate in our own galaxy is between 11 and 46/yr, and (3) nova rates are not strongly affected by tidal captures in bright globular clusters.

Ciardullo, Robin↗
NASA NTRS1994

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully nonlinear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate that that from ZA to TZA, (b) the error in the phase angle of Fourier components is worse that that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗
NASA NTRS2017

An Orbital Meteoroid Stream Survey Using the Southern Argentina Agile Meteor Radar (SAAMER) Based on a Wavelet Approach

Over a million individually measured meteoroid orbits were collected with the Southern Argentina Agile MEteor Radar (SAAMER) between 2012-2015. This provides a robust statistical database to perform an initial orbital survey of meteor showers in the Southern Hemisphere via the application of a 3D wavelet transform. The method results in a composite year from all 4 years of data, enabling us to obtain an undisturbed year of meteor activity with more than one thousand meteors per day. Our automated meteor shower search methodology identified 58 showers. Of these showers, 24 were associated with previously reported showers from the IAU catalogue while 34 showers are new and not listed in the catalogue. Our searching method combined with our large data sample provides unprecedented accuracy in measuring meteor shower activity and description of shower characteristics in the Southern Hemisphere. Using simple modeling and clustering methods we also propose potential parent bodies for the newly discovered showers.

Meteor showers↗