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 973 records · Page 54

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.↗

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↗

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.↗

Joint Estimation of Topology and Injection Statistics in Distribution Grids with Missing Nodes

Optimal operation of distribution grid resources relies on accurate estimation of its state and topology. Practical estimation of such quantities is complicated by the limited presence of real-time meters. This article discusses a theoretical framework to jointly estimate the operational topology and statistics of injections in radial distribution grids under limited availability of nodal voltage measurements. In particular, we show that our proposed algorithms are able to provably learn the exact grid topology and injection statistics at all unobserved nodes as long as they are not adjacent. The algorithm design is based on novel ordered trends in voltage magnitude fluctuations at node groups, that are independently of interest for radial physical flow networks. The complexity of the designed algorithms is theoretically analyzed and their performance is validated using both linearized and nonlinear ac power flow samples in test distribution grids.

97 MATHEMATICS AND COMPUTING↗

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↗

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully non-linear, 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 than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than 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.↗

Hyaloscypha finlandica Metabolome Repository

This repository provides the curated data tables, manuscript figure and table exports, dependency records, and workflow scripts supporting an integrated comparative genomics and untargeted LC-MS/MS metabolomics analysis of Hyaloscypha finlandica strain PMI 746, a root-associated dark septate endophyte of poplar. The repository includes genome-mining summaries from antiSMASH, FunBGCeX, BGC-Prophet, and BiG-SCAPE; processed metabolomics inputs; metabolite annotation evidence; statistical outputs; and publication-facing figures and tables. Raw LC-MS/MS spectra, full genome/protein downloads, and large generated tool outputs are referenced through public archive/accession records and are not stored in Git.

59 BASIC BIOLOGICAL SCIENCES↗

Astrometric Detection of Extrasolar Planets: Results of a Feasibility Study with the Palomar 5 Meter Telescope

The detection of extrasolar planets around stars like the Sun remains an important goal of astronomy. We present results from Palomar 5 m observations of the open cluster NGC 2420 in which we measure some of the sources of noise that will be present in an astrometric search for extrasolar planets. This is the first time that such a large aperture has been used for high-precision astrometry. We find that the atmospheric noise is 150 micro-arcsec hr(exp 1/2) across a 90 sec field of view and that differential chromatic refraction (DCR) can be calibrated to 128 micro-arcsec for observations within 1 hr of the meridian and 45 deg of zenith. These results confirm that a model for astrometric measurements can be extrapolated to large apertures. We demonstrate, based upon these results, that a large telescope achieves the sensitivity required to perform a statistically significant search for extra solar planets. We describe an astrometric technique to detect planets, the astrometric signals expected, the role of reference stars, and the sources of measurement noise: photometric noise, atmospheric motion between stars, sky background, instrumental noise, and DCR. For the latter, we discuss a method to reduce the noise further to 66 micro-arcsecond for observations within 1 hr of the meridian and 45 deg of zenith. We discuss optimal lists of target stars taken from the latest Gliese & Jahreiss catalog of nearby stars with the largest potential astrometric signals, declination limits for both telescope accessibility and reduced DCR, and galactic latitude limits for a sufficiant number of reference stars. Two samples are described from which one can perform statistically significant searches for gas giant planets around nearby stars. One sample contains 100 "solar class" stars with an average stellar mass of 0.82 solar mass; the other maximizes the number of stars, 574, by searching mainly low-mass M stars. We perform Monte Carlo simulations of the statistical significance of the expected results by using measured and estimated noise quantities. We show the semimajor axis parameter spaces that are searched for each star and how an increase in the length of the observing program expands these spaces. The search over semimajor axis parameter space relates to the theory of gas giant planet formation.

Pravdo, Steven H.↗

Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators

We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high k modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss.

79 ASTRONOMY AND ASTROPHYSICS↗

Lensing without borders – I. A blind comparison of the amplitude of galaxy–galaxy lensing between independent imaging surveys

ABSTRACT Lensing without borders is a cross-survey collaboration created to assess the consistency of galaxy–galaxy lensing signals (ΔΣ) across different data sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of ΔΣ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported systematic errors which agree to better than 2.3σ in four lens bins and three radial ranges. For lenses with zL > 0.43 and considering statistical errors, we detect a 3–4σ correlation between lensing amplitude and survey depth. This correlation could arise from the increasing impact at higher redshift of unrecognized galaxy blends on shear calibration and imperfections in photometric redshift calibration. At zL > 0.54, amplitudes may additionally correlate with foreground stellar density. The amplitude of these trends is within survey-defined systematic error budgets that are designed to include known shear and redshift calibration uncertainty. Using a fully empirical and conservative method, we do not find evidence for large unknown systematics. Systematic errors greater than 15 per cent (25 per cent) ruled out in three lens bins at 68 per cent (95 per cent) confidence at z < 0.54. Differences with respect to predictions based on clustering are observed to be at the 20–30 per cent level. Our results therefore suggest that lensing systematics alone are unlikely to fully explain the ‘lensing is low’ effect at z < 0.54. This analysis demonstrates the power of cross-survey comparisons and provides a promising path for identifying and reducing systematics in future lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Smooth Pursuit Saccade Amplitude Modulation During Exposure to Microgravity

Russian investigators have reported changes in pursuit tracking of a vertically moving point stimulus during space flight. Early in microgravity, changes were manifested by decreased eye movement amplitude (undershooting) and the appearance of correction saccades. As the flight progressed, pursuit of the moving point stimulus deteriorated while associated saccadic movements were unchanged. Immediately postflight there was an improved execution of active head movements indicating that the deficiencies in pursuit function noted in microgravity may be of central origin. In contrast, tests of two cosmonauts showed that horizontal and vertical smooth pursuit were unchanged inflight. However, results of corresponding saccadic tasks showed a tendency toward the overshooting of a horizontal target early inflight with high accuracy developing later inflight, accompanied by an increased saccade velocity and a trend toward decreased saccade latency. Based on these equivocal results, we have further investigated the effects of space flight on the smooth pursuit mechanism during and after short duration flight, and postflight on returning MIR crewmembers. Sinusoidal target movement was presented horizontally at frequencies of 0.33 and 1.0 Hz. Subjects were asked to perform two trials for each stimulus combination: (1) moving eyes-only (EO) and (2) moving eyes and head (EH) with the target motion. Peak amplitude was 30 deg for 0.33 Hz trials and 15 deg for the 1.0 Hz trials. The relationship between saccade amplitude and peak velocity were plotted as a main sequence for each phase of flight, and linear regression analysis allowed us to determine the slope of each main sequence plot. The linear slopes were then combined for each flight phase for each individual subject. The main sequence for both EO and EH trials at both the 0.33 and 1.0 Hz frequencies during flight for the short duration flyers showed a reduction in saccade velocity and amplitude when compared to the preflight main sequence . This difference in the regression slopes between flight phase, head/eye condition (EO or EH), and pursuit target frequency was observed across all subjects (statistically significant at the p<0.02, df= 2). It is interesting to note that postflight for the short duration flyers there was an immediate recovery to the preflight main sequence across all trials. There were no significant differences observed between the preflight slopes for either head movement condition (EO vs. EH). When the immediate postflight (R+O) regression slopes were compared with the preflight slopes, there was a tendency (not significant) for both saccade amplitude and peak velocity to increase during the postflight testing. This tendency had vanished by R+ 1. Of particular interest was the redistribution of saccades during the latter stages of the flight and immediately postflight in the EO condition. At the 1.0 Hz frequency the saccades tended to be clustered near the lowest target velocity. It was also interesting to note that gaze performance (eye in skull + head in space) was consistently better during the EH condition; a finding also observed by our Russian colleagues. As the results of the long duration flight become available we expect that they will not only show that postflight effects will be similar to those observed during the short duration flights, but will also last for a greater period of time following flight. It is not clear what mechanism is responsible for the decreased peak saccadic velocity during flight unless the change is related to the control of retinal slip. For example, it is possible that saccades will tend to initially undershoot their targets by a small percentage and these saccades are then followed, if vision is available, by a small augmenting corrective saccade. It has been postulated that the functional significance of this undershooting tendency is to maintain the spatial representation of the target on the same side of the fovea (as opposedo racing across the fovea) and hence in the same cerebral hemisphere that initiated the primary saccade thus minimizing delays caused by an intra-hemispheric transfer of information . One could also speculate that with saccade velocities greater than normal, additional corrective saccades would be required to bring the target back on the fovea. A less plausible explanation of our findings could be fatigue. Yet it seems unlikely that our subjects would show lower velocities on all inflight test days while showing increased saccade velocities immediately following space flight where fatigue is usually the greatest. Finally, the redistribution effect noted late inflight is likely caused by adaptive changes. Overall, corrective saccades appeared to be used in maintaining gaze on target; reducing retinal slip and assisting space travelers in maintaining clear vision throughout the different phases of the space flight.

Reschke, M. F.↗

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey quasar sample: testing observational systematics on the Baryon Acoustic Oscillation measurement

ABSTRACT Baryon Acoustic Oscillations (BAO) are considered to be a very robust standard ruler against various systematics. This premise has been tested against observational systematics, but not to the level required for the next generation of galaxy surveys such as the Dark Energy Spectroscopic Instrument (DESI) and Euclid. In this paper, we investigate the effect of observational systematics on the BAO measurement of the final sample of quasars from the extended Baryon Oscillation Spectroscopic Survey Data Release 16 in order to prepare and hone a similar analysis for upcoming surveys. We employ catalogues with various treatments of imaging systematic effects using linear and neural network-based non-linear approaches and consider how the BAO measurement changes. We also test how the variations to the BAO fitting model respond to the observational systematics. As expected, we confirm that the BAO measurements obtained from the DR16 quasar sample are robust against imaging systematics well within the statistical error, while reporting slightly modified constraints that shift the line-of-sight BAO signal by less than 1.1 per cent. We use realistic simulations with similar redshift and angular distributions as the DR16 sample to conduct statistical tests for validating the pipeline, quantifying the significance of differences, and estimating the expected bias on the BAO scale in future high-precision data sets. Although we find a marginal impact for the eBOSS QSO data, the work presented here is of vital importance for constraining the nature of dark energy with the BAO feature in the new era of big data cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Crash Risks Evaluation of Urban Expressways: A Case Study in Shanghai

We report that proactive traffic safety management systems can reduce crashes by identifying crash precursors, evaluating real-time crash risks, and implementing suitable interventions. The basic prerequisite for developing such a system is to propose a reliable crash risk evaluation model that takes real-time traffic flow data as input. Previous studies have primarily focused on real-time crash prediction using some statistical or machine-learning methods. However, further quantitative evaluation and classification of crash risks have been ignored. In this study, we conduct a systematic crash risk evaluation workflow, including crash risk prediction, crash risk quantification, and crash risk classification. Specifically, the crash risk prediction using an extended logit model is proposed, from which CAS, CSD, UAS, DAS, DTV are identified to be contributing factors of crash risks. Then a crash risk quantification model based on the parameter evaluation of the extended logit model is developed. The crash risks of urban expressways and their spatial-temporal evolution trends are quantified. Finally, the crash risks are classified into high crash risk level, moderate crash risk level, and low crash risk level by the k-means cluster algorithm. Then the threshold boundaries of different crash risk levels are determined. The research results provide a proactive guidance for traffic safety management of urban expressways.

33 ADVANCED PROPULSION SYSTEMS↗

Condition-Based Maintenance of a Circulating Water System of a Canadian Nuclear Power Plant using Machine Learning and Statistical Tools

Canada Deuterium Uranium pressurized-heavy-water reactors (PHWR) are a type of nuclear power plant that generate clean and reliable energy. The scope of this work is to automate data analysis methodologies to inform a condition-based maintenance strategy of a circulating water system (CWS) of a PHWR. The multiunit CWS provides a continuous supply of water to cool steam condensers, even during transient scenarios, thereby improving the thermal efficiency. This work aims to develop a machine learning (ML) based approach to detect anomalies in heterogeneous data of a CWS in a PHWR to help inform a predictive maintenance strategy. The heterogeneous data include textual and numeric time series data for a PHWR. Natural-language-processing (NLP)-based models are used to analyze textual data contained in work orders and operator logs and an event-timeseries correlation detection method is applied to assist anomalies diagnoses for CWS. An ML model Robust Linear Model (RLM) is also used to remove the seasonal variations in the system variable distributions based on distributions of environmental variables. A machine learning model, Density-Based Spatial Clustering of Applications with Noise (DBSCAN), trained on both original data and data without any seasonal variations will then be used to detect if an anomaly exists. Thus, by moving to an automated methodology to detect, classify, and forecast anomalies, the maintenance strategy would be based on component condition instead of a time-based schedule.

97 - MATHEMATICS AND COMPUTING↗

The formation of cosmic structure in a texture-seeded cold dark matter cosmogony

The growth of density fluctuations induced by global texture in an Omega = 1 cold dark matter (CDM) cosmogony is calculated. The resulting power spectra are in good agreement with each other, with more power on large scales than in the standard inflation plus CDM model. Calculation of related statistics (two-point correlation functions, mass variances, cosmic Mach number) indicates that the texture plus CDM model compares more favorably than standard CDM with observations of large-scale structure. Texture produces coherent velocity fields on large scales, as observed. Excessive small-scale velocity dispersions, and voids less empty than those observed may be remedied by including baryonic physics. The topology of the cosmic structure agrees well with observation. The non-Gaussian texture induced density fluctuations lead to earlier nonlinear object formation than in Gaussian models and may also be more compatible with recent evidence that the galaxy density field is non-Gaussian on large scales. On smaller scales the density field is strongly non-Gaussian, but this appears to be primarily due to nonlinear gravitational clustering. The velocity field on smaller scales is surprisingly Gaussian.

Gooding, Andrew K.↗

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the minimal bias model

Here, we present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg 2 , and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, w p ⁡(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δ⁢Σ⁡(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ξ ± ⁡($\vartheta$) measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat Λ⁢CDM cosmology. We use conservative scale cuts, R >12 and 8⁢ℎ -1 Mpc for Δ⁢Σ and w p , respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S 8 = $0.77⁢5^{+0.043}_{-0.038}$ (68% C.I.) for the Λ⁢CDM model, after marginalizing over uncertainties in other parameters. Our value of S 8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δ⁢Σ⁡(R) and w p ⁡(R) down to R >3 and 2⁢ℎ -1 Mpc, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational Inference of Vibratory System with Incomplete Modal Information Using Parallel, Interactive and Adaptive Markov Chains

Inverse analysis of vibratory system is an important subject in fault identification, model updating, and robust design and control. It is challenging subject because 1) the problem is oftentimes underdetermined while the measurements are limited and/or incomplete; 2) many combinations of parameters may yield results that are similar with respect to actual response measurements; and 3) uncertainties inevitably exist. The aim of this research is to leverage upon computational intelligence through statistical inference to facilitate an enhanced, probabilistic framework using incomplete modal response measurement. This new framework is built upon efficient inverse identification through optimization, whereas Bayesian inference is employed to account for the effect of uncertainties. To overcome the computational cost barrier, we adopt Markov chain Monte Carlo (MCMC) to characterize the target function/distribution. Instead of using single Markov chain in conventional Bayesian approach, we develop a new sampling theory with multiple parallel, interactive and adaptive Markov chains and incorporate it into Bayesian inference. This can harness the collective power of these Markov chains to realize the concurrent search of multiple local optima. The number of required Markov chains and their respective initial model parameters are automatically determined via Monte Carlo simulation-based sample pre-screening followed by K-means clustering analysis. These enhancements can effectively address the aforementioned challenges in finite element inverse analysis. The validity of this framework is systematically demonstrated through case studies.

K Zhou↗

4.32 Billion Year Old Impact Melts at Apollo 14: Dating the Procellarum Basin?

Introduction: Recent work has improved our understanding of lunar crustal structure and basin evolution on the Moon, but its early impact history and implications for solar system dynamics remain poorly established [1,2]. In particular, there is a derth of absolute ages, especially for the period >4.0 Ga, that could be related to basin-forming events. Here we present U-Pb isotopic data for Zr-rich minerals found in impact-melt fragments from Apollo 14 soil sample 14163 that, when combined with previously published data, dates an impact event at 4324±15 Ma. The compositions and ages of these impact-melt fragments provide unique information about the timing of early impact events and the composition of the lunar crust. Results: Fourteen rocklets ranging in size from 1 to 3 mm were extracted from soil 14163. Most are impactmelt rocks with 10-30% clasts of mostly pyroxene and plagioclase (100 to 500 μm) in a crystalline matrix. Textures of the matrix vary from subophitic, formed by intergrowths of 10-20 μm plagioclase and pyroxene crystals, to poikilitic, with plagioclase and pyroxene reaching 20-50 μm size. Some fragments that contain no visible clasts have similar textures are also interpreted as impact-melt rocks. All fragments contain notably large proportions of euhedral to subhedral ilmenite, zircon, apatite/merrillite and less abundant zirconolite and baddeleyite in the melt matrices. These grains often form intricate intergrowths with each other and rock-forming minerals, indicating their crystallization from the melt (Fig. 1). However, some slightly larger (~50 μm) zircon and phosphate grains can be interpreted as relict clasts based on their granular textures and relationships with the surrounding phases (Fig. 1). This implies that zircon and phosphate minerals were present in the target rocks. Some of these grains were profoundly remelted during the impact, which resulted in an oversaturation of the melt in Zr and P and crystallization of new grains of Zrrich minerals and phosphates during melt solidification. The presence of zircon and phosphates indicates that the melts were significantly enriched in KREEP components and ties their provenance to the Procellarum- KREEP Terrane (PKT) [3]. Combined U-Pb zircon data obtained for all fragments form two main clusters on a concordia diagram (Fig. 2), between about 4.3 and 3.9 Ga. Phosphate data mostly concentrate on the younger end of this age range but some analyses are almost as old as the older zircon grains (Fig. 2). Figure 1: Zircon (Zr) and phosphate (Apt; Mer) grains in impact melt fragments from Apollo 14 soil 14161. A-zircon grains crystallized from impact melt; B-phosphate grain inherited from the target; C- granular zircon grains Our interpretation of these data is that all U-rich minerals experienced variable resetting of the U-Pb system, with phosphate, where closure temperature is significantly lower than that in zircon, affected more profoundly by Pb loss than zircon. Our best estimate of the time of formation of zircon and phosphate is based on statistically valid analysis of grains interpreted texturally as grown from the impact melt (Fig. 2). Ten of these analyses define an age of 4324±15 Ma (MSWD=3.0, probability of fit P=0.002). Our best estimate for the time of resetting is obtained by combining data from phosphates that are statistically indistinguishable from 3.9 Ga within the analytical uncertainties. This group is represented by 18 analyses of 15 phosphate grains from different fragments (Fig. 2) and defines an age of 3922±6 Ma (MSWD=1.2, P=0.23). Conversely, a minimum age of the target lithologies, remelted in the impact that produced the rocklets studied here (Fig. 2), can be determined from the five oldest analyses of zircon clasts at 4338±13 Ma (MSWD=1.5, P=0.2), which is indistinguishable from the age of the impact melt within the uncertainties. The obtained ages, combined with textural evidence, imply that the impact melt was formed at 4324±15 Ma, and that it occurred in a zircon-rich target with a minimum age of 4338±13 Ma. Further reworking occurred during a second impact event at 3922±6 Ma. Discussion: The younger age of 3922±6 Ma can be interpreted as the time of the Imbrium impact. The older age of 4324±15 Ma would then be the time of formation of the impact melt, which was then caught in the Imbrium ejecta either at the Apollo 14 landing site or within the target rocks of the Imbrium impact. The presence of abundant accessory phases such as zircons and phosphates is consistent with a substantial KREEP component in the analyzed particles. The current distribution of KREEP on the lunar surface appears to be strongly influenced by Imbrum ejecta [4] so interpretation of the 4.32 Ga age depends in part on assumed structure of the pre-impact crust and distribution of KREEP within the crust. If KREEP was present only in the deep crust at 4.32 Ga, then a basin-scale impact possibly analogous to Imbrium or larger may be necessary to excavate a KREEPy impact melt at this time. However, if KREEPy materials were present closer to the surface perhaps due to redistribution related to Mg-suite magmatism, then smaller impacts might be able to rework KREEPy compositions at shallower depths. The coherence of the data on the particles analysed here suggests a large volume of melt that has been preserved since 4.32 Ga, consistent with a large impact event. Relict zircons and mineral clasts suggests that the igneous crust in the vicinity of this impact was well developed by at least 4.34 Ga, similar to the model age of KREEP and older than the isochron ages of many Mg-suite cumulates [5]. A problem that confronts all lunar sample studies using the current collection is that the pre-Imbrium geology of the PKT (the source of these 4.32 Ga impact melts) is not well constrained. The South Pole-Aitken basin contains regions that are moderately enriched in Th, but its ejecta is Th-poor [6]; therefore these A14 fragments are probably not SPA ejecta. Alternatively, these fragments may represent formation of a hypothesized Procellarum basin [7] although the lack of a clearly defined basin ejecta signature is a potential problem with that interpretation. In any case, they provide a younger limit on the age of lunar differentiation and formation of KREEP within the lunar crust. Figure 2: U-Pb data for zircon and phosphate grains from impact melt fragments. A-all data; B-data used for age calculations. References: [1] Orgel C., Michael G., Fassett C. I., van der Bogert C. H., Riedel C., Kneissl T., and Hiesinger H. (2018) J. Geophys. Res. Planets 123, 748- 762. [2] Evans A. J., Andrews-Hanna J. C., Head J. W., Soderblom J. M., Solomon S. C., and Zuber M. T. (2018) J. Geophys. Planets. 123, 1596-1617. [3] Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L. and Wieczorek, M.A. (2000) J. Geophys. Res: Planets 105, 4197-4216. [4] Haskin L. A. (1998) J. Geophys. Res. Planets 103, 1679-1689. [5] Borg L.E., Gaffney A.M., and Shearer C.K. (2015) MAPS 50, 715-732. [6] Moriarty, D.P., Watkins, R.N., Valencia, S.N., Kendall, J.D., Evans, A.J., Dygert, N. and Petro, N.E. (2021) J. Geophys. Res. Planets 126. [7] Zhu, M.H., Wünnemann, K., Potter, R.W., Kleine, T. and Morbidelli, A. (2019) J. Geophys. Res: Planets 124, 2117-2140

M D Norman↗