SEARCH · Engineering Papers

Results for “Filter error”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,531 records · Page 85

All sources NASA NTRS DOE OSTI

NASA NTRS1999

A Diagnostic Analysis of the Kennedy Space Center LDAR Network: Data Characteristics - 1

An analytical framework is developed in which to analyze climatological VHF (66 MHz) radiation measurements taken by the Kennedy Space Center Lightning Detection and Ranging (LDAR) network. A 19-month noise-filtered sample of LDAR observations is examined using this framework. It is found that the climatological VHF source density as observed by LDAR falls off approximately 10 dB every 71 km of ground range away from the network centroid (a 31 km e-folding scale). From this framework it is inferred that the underlying source distribution is likely inverse exponential in amplitude, with a log-slope of 126-147 mW(sup -1/2). The underlying vertical distribution of VHF sources is approximately normally distributed with a mean altitude of 9 km and a standard deviation of 2.7 km; this implies that the loss of below-horizon sources has a negligible effect on column-integrated source densities within 200 km ground range. At medium to far ranges, location errors are primarily radial and have a slightly asymmetric distribution whose first moment increases as range squared. Error moments estimated from observed lightning are significantly higher than those from aircraft-based signal generator or analytic solution estimates, suggesting that timing errors arising from poor signal identification and discrimination may dominate over timing errors arising from nominal sensor resolution. The VHF source properties of individual LDAR-observed flashes are computed and an analytic expression for flash detection efficiency vs. range is derived. This reveals nearly constant flash detection efficiency to 80-94 km range from the network centroid.

Boccippio, D. J.↗

DOE OSTI2021

Early season prediction of within-field crop yield variability by assimilating CubeSat data into a crop model

Accurate early season predictions of crop yield at the within-field scale can be used to address a range of crop production, management, and precision agricultural challenges. While the remote sensing of within-field insights has been a research goal for many years, it is only recently that observations with the required spatio-temporal resolutions, together with efficient assimilation methods to integrate these into modeling frameworks, have become available to advance yield prediction efforts. Here we explore a yield prediction approach that combines daily high-resolution CubeSat imagery with the APSIM crop model. The approach employs APSIM to train a linear regression that relates simulated yield to simulated leaf area index (LAI). That relationship is then used to identify the optimal regression date at which the LAI provides the best prediction of yield: in this case, approximately 14 weeks prior to harvest. Instead of applying the regression on satellite imagery that is coincident, or closest to, the regression date, our method implements a particle filter that integrates CubeSat-based LAI into APSIM to provide end-of-season high-resolution (3 m) yield maps weeks before the optimal regression date. The approach is demonstrated on a rainfed maize field located in Nebraska, USA, where suitable collections of both imagery and in-situ data were available for assessment. The procedure does not require in-field data to calibrate the regression model, with results showing that even with a single assimilation step, it is possible to provide yield estimates with good accuracy up to 21 days before the optimal regression date. Yield spatial variability was reproduced reasonably well, with a strong correlation to independently collected measurements (R 2 = 0.73 and rRMSE = 12%). When the field averaged yield was compared, our approach reduced yield prediction error from 1 Mg/ha (control case based on a calibrated APSIM model), to 0.5 Mg/ha (using satellite imagery alone), and then to 0.2 Mg/ha (results with assimilation up to three weeks prior to the optimal regression date). Such a capacity to provide spatially explicit yield predictions early in the season has considerable potential to enhance digital agricultural goals and improve end-of-season yield predictions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Advanced Multilevel Monte Carlo Methods

This article reviews the application of some advanced Monte Carlo techniques in the context of multilevel Monte Carlo (MLMC). MLMC is a strategy employed to compute expectations, which can be biassed in some sense, for instance, by using the discretization of an associated probability law. The MLMC approach works with a hierarchy of biassed approximations, which become progressively more accurate and more expensive. Additionally, using a telescoping representation of the most accurate approximation, the method is able to reduce the computational cost for a given level of error versus i.i.d. sampling from this latter approximation. All of these ideas originated for cases where exact sampling from couples in the hierarchy is possible. This article considers the case where such exact sampling is not currently possible. We consider some Markov chain Monte Carlo and sequential Monte Carlo methods, which have been introduced in the literature, and we describe different strategies that facilitate the application of MLMC within these methods.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2004

Implicit and Explicit Spacecraft Gyro Calibration

This paper presents a comparison between two approaches to sensor calibration. According to one approach, called explicit, an estimator compares the sensor readings to reference readings, and uses the difference between the two to estimate the calibration parameters. According to the other approach, called implicit, the sensor error is integrated to form a different entity, which is then compared with a reference quantity of this entity, and the calibration parameters are inferred from the difference. In particular this paper presents the comparison between these approaches when applied to in-flight spacecraft gyro calibration. Reference spacecraft rate is needed for gyro calibration when using the explicit approach; however, such reference rates are not readily available for in-flight calibration. Therefore the calibration parameter-estimator is expanded to include the estimation of that reference rate, which is based on attitude measurements in the form of attitude-quaternion. A comparison between the two approaches is made using simulated data. It is concluded that the performances of the two approaches are basically comparable. Sensitivity tests indicate that the explicit filter results are essentially insensitive to variations in given spacecraft dynamics model parameters.

Bar-Itzhack, Itzhack Y.↗

NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗

DOE OSTI2023

The effect of relative humidity on eddy covariance latent heat flux measurements and its implication for partitioning into transpiration and evaporation

While the eddy covariance (EC) technique is a well-established method for measuring water fluxes (i.e., evaporation or 'evapotranspiration’, ET), the measurement is susceptible to many uncertainties. One such issue is the potential underestimation of ET when relative humidity (RH) is high (> 70%), due to low-pass filtering with some EC systems. Yet, this underestimation for different types of EC systems (e.g. open-path or closed-path sensors) has not been characterized for synthesis datasets such as the widely used FLUXNET2015 dataset. Here, we assess the RH-associated underestimation of latent heat fluxes (LE, or ET) from different EC systems for 163 sites in the FLUXNET2015 dataset. We found that the LE underestimation is most apparent during hours when RH is higher than 70%, predominantly observed at sites using closed-path EC systems, but the extent of the LE underestimation is highly site-specific. We then propose a machine learning based method to correct for this underestimation, and compare it to two energy balance closure based LE correction approaches (Bowen ratio correction, BRC, and attributing all errors to LE). Our correction increases LE by 189% for closed-path sites at high RH (> 90%), while BRC increases LE by around 30% for all RH conditions. Additionally, we assess the influence of these corrections on ET-based transpiration (T) estimates using two different ET partitioning methods. Results show opposite responses (increasing vs. slightly decreasing T-to-ET ratios, T/ET) between the two methods when comparing T based on corrected and uncorrected LE. Overall, our results demonstrate the existence of a high RH bias in water fluxes in the FLUXNET2015 dataset and suggest that this bias is a pronounced source of uncertainty in ET measurements to be considered when estimating ecosystem T/ET and WUE.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1973

Application of remote sensing in the study of vegetation and soils in Idaho

The author has identified the following significant results. Comparison of ERTS-1 imagery and USGS 1:250,000 scale maps of study areas with known ground points revealed significant map errors. These errors were sufficient to render impractical the projection of ERTS-1 imagery directly onto maps of the area. Marked differences were found in the delineation of ground features by different MSS bands. Generally, Band 4 was least useful, while Band 5 proved valuable for indicating patterns of native vegetation, cultivated areas - both dry and irrigated, lava fields, drainage basins, and deep bodies of water. Band 6 was better for landforms and drainages and for shallow bodies of water than Band 5 but inferior for indicating patterns in native vegetation and most types of cultivated land. Band 7 was best of all for indicating lava flows, water bodies, and landform features. Use of a additive color viewer-projector aided greatly in separation of images. A combination of Bands 5 and 7 with appropriate color filters proved best for separating most types of native vegetation and cultivated crops. Landform features and water bodies also showed well with this combination. The addition of Band 4 imagery to these further enhanced the identification of semi-dormant vegetation.

Tisdale, E. W.↗

NASA NTRS2009

Metrology System for a Large, Somewhat Flexible Telescope

A proposed metrology system would be incorporated into a proposed telescope that would include focusing optics on a rigid bench connected via a deployable mast to another rigid bench holding a focal-plane array of photon counting photodetectors. Deformations of the deployable mast would give rise to optical misalignments that would alter the directions (and, hence, locations) of incidence of photons on the focal plane. The metrology system would measure the relative displacement of the focusing- optics bench and the focal-plane array bench. The measurement data would be used in post-processing of the digitized photodetector outputs to compensate for the mast-deformation-induced changes in the locations of incidence of photons on the focal plane, thereby making it possible to determine the original directions of incidence of photons with greater accuracy. The proposed metrology system is designed specifically for the Nuclear Spectroscopic Telescope Array (NuSTAR) a proposed spaceborne x-ray telescope. The basic principles of design and operation are also applicable to other large, somewhat flexible telescopes, both terrestrial and spaceborne. In the NuSTAR, the structural member connecting the optical bench and the photodetector array would be a 10-m-long deployable mast, and there is a requirement to keep errors in measured directions of incidence of photons below 10 arc seconds (3 sigma). The proposed system would include three diode lasers that would be mounted on the focusing-optics bench. For clarity, only one laser is shown in the figure, which is a greatly simplified schematic diagram of the system. Each laser would be aimed at a position-sensitive photodiode that would be mounted on the detector bench alongside the aforementioned telescope photodetector array. The diode lasers would operate at a wavelength of 830 nm, each at a power of 200 mW. Each laser beam would be focused to a spot of .1-mm diameter on the corresponding position-sensitive photodiode. To reduce the effect of sunlight on the measurements, a one-stage light baffle and an 830-nm transmission filter of 10-nm bandwidth would be placed in front of the position- sensitive photodiode. For each metrology reading, the output of the position-sensitive detector would be sampled and digitized twice: once with the lasers turned on, then once with the lasers turned off. The data from these two sets of samples would be subtracted from each other to further reduce the effects of sun glints or other background light sources.

Liebe, Carl Christian↗

DOE OSTI2020

GridOPTICS/GridPACK

GridPACK is a software framework consisting of a set of modules designed to simplify the development of programs that model the power grid and run on parallel, high performance computing platforms. It also contains several fully developed applications, including powerflow, dynamic simulation, state estimation, Kalman filter analysis (dynamic state estimation), contingency analysis and real time path rating. These applications can be used either standalone or as components in more complicated workflows that combine several different types of application together. The framework modules are available as a combination of libraries and software templates and consist of components for setting up and distributing power grid networks, support for modeling the behavior of individual buses and branches in the network, converting the network models to the corresponding algebraic equations, and parallel routines for manipulating and solving large algebraic systems. The framework also contains a module for distributing tasks evenly amongst computing resources, even if individual tasks vary widely in their execution times. Additional modules support input and output, basic statistical analysis of contingency based calculations, distributed data structures, as well as basic profiling and error management.

Palmer, Bruce↗

NASA NTRS1984

Specular reflectance of optical-black coatings in the far infrared

Far-infrared specular reflectance spectra of seven optically black coatings near normal incidence are presented. Seven photometric spectra were obtained using eleven bandpass transmission filters in the wavelength range between 12 and 500 microns, and three interferometric spectra were obtained for corroboration. Data on the construction, thickness, and rms surface roughness of the coatings are also presented. The chemical composition of three coatings can be distinguished from that of the others by a strong absorption feature between 20 and 40 microns, which can be largely attributed to amorphous silicate material. At 100 microns, the most and least reflective coatings differ by nearly 3 orders of magnitude. Inverse relationships observed between the spectra and the roughness and thickness of the coatings led to development of a reflecting-layer model for the measured reflectance. The model successfully describes the spectra at wavelengths outside the silicate absorption, and optical constants are deduced from a nonlinear least squares fit to the data. Parametric errors are estimated by chi-square analysis, and sensitivity tests are performed to determine which parameters control reflectance in different spectral regions.

Smith, S. M.↗

NASA NTRS2010

Solar Array Disturbances to Spacecraft Pointing During the Lunar Reconnaissance Orbiter (LRO) Mission

The Lunar Reconnaissance Orbiter (LRO), the first spacecraft to support NASA s return to the Moon, launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle. It was initially inserted into a direct trans-lunar trajectory to the Moon. After a five day transit to the Moon, LRO was inserted into the Lunar orbit and successfully lowered to a low altitude elliptical polar orbit for spacecraft commissioning. Successful commissioning was completed in October 2009 when LRO was placed in its near circular mission orbit with an approximate altitude of 50km. LRO will spend at least one year orbiting the Moon, collecting lunar environment science and mapping data, utilizing a suite of seven instruments to enable future human exploration. The objective is to provide key science data necessary to facilitate human return to the Moon as well as identification of opportunities for future science missions. LRO's instrument suite will provide the high resolution imaging data with sub-meter accuracy, highly accurate lunar cartographic maps, mineralogy mapping, amongst other science data of interest. LRO employs a 3-axis stabilized attitude control system (ACS) whose primary control mode, the "Observing Mode", provides Lunar nadir, off-nadir, and inertial fine pointing for the science data collection and instrument calibration. This controller combines the capability of fine pointing with on-demand large angle full-sky attitude reorientation. It provides simplicity of spacecraft operation as well as additional flexibility for science data collection. A conventional suite of ACS components is employed in the Observing Mode to meet the pointing and control objectives. Actuation is provided by a set of four reaction wheels developed in-house at NASA Goddard Space Flight Center (GSFC). Attitude feedback is provided by a six state Kalman filter which utilizes two SELEX Galileo Star Trackers for attitude updates, and a single Honeywell Miniature Inertial Measurement Unit (MIMU) to provide body rates for attitude propagation. Rate is computed by differentiating accumulated angle provided by the MIMU. The Observing Mode controller is required to maintain fine pointing while a large fully-articulated solar array (SA) maintains its panel normal to the solar incidence. This paper describes the disturbances to the attitude control resulting from the SA articulation. Observing Mode performance in the presence of this disturbance was assessed while the spacecraft was in an initial elliptical low altitude orbit during the commissioning phase, which started about two weeks after launch and lasted for 90 days. LRO demonstrated excellent pointing performance during Observing Mode nadir and inertial attitude target operations during this phase. Transient LRO attitude errors observed during commissioning resulted primarily from three sources, Diviner instrument calibrations, RW zero crossings, and SA articulation. Even during times of considerable disturbance from SA articulation, the attitude errors were maintained below the statistical attitude error requirement level of 15 arc-sec (3 sigma).

Calhoun, Philip↗

DOE OSTI2020

The good, the bad, and the ugly: Data-driven load profile discord identification in a large building portfolio

Reducing the overall energy consumption and associated greenhouse gas emissions in the building sector is essential for meeting our future sustainability goals. Recently, smart energy metering facilities have been deployed to enable monitoring of energy consumption data with hourly or subhourly temporal resolution. This unprecedented data collection has created various opportunities for advanced data analytics involving load profiles (e.g., building energy benchmarking programs, building-to-grid integration, and calibration of urban-scale energy models). These applications often need preprocessing steps to detect daily load profile discords, such as: 1) outliers due to system malfunctions (the bad) and 2) irregular energy consumption patterns, such as those resulting from holidays (the ugly) compared to normal consumption patterns (the good). However, current preprocessing methods predominantly focus on filtering using statistical threshold values, which fail to capture the contextual discords of daily profiles. In addition, discord detection algorithms in building research are often aimed at finding individual building-level discords, which are not suitable at a large scale. Thus, here, we develop a method for automated load profile discord identification (ALDI) in a large portfolio of buildings (more than 100 buildings). Specifically, ALDI 1) uses the matrix profile (MP) method to quantify the similarities of daily subsequences in time series meter data, 2) compares daily MP values with typical-day MP distributions using the Kolmogorov-Smirnov test, and 3) identifies daily load profile discords in a large building portfolio. We evaluate ALDI using the metering data of both an academic campus and a residential neighborhood. Our results demonstrate that ALDI efficiently discovers measurement errors by system malfunctions and low energy consumption days in the academic campus portfolio, and it detects unique load shape patterns likely driven by occupant behavior and extreme weather conditions in the residential neighborhood.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1977

The 4th order GISS model of the global atmosphere

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Kalnay-Rivas, E.↗

NASA NTRS1999

HRDI Observations of Inertia-Gravity Waves in the Mesosphere and Lower Thermosphere

Vertical profiles of High-resolution Doppler imager (HRDI) mesospheric winds have small-scale structure (vertical wavelengths between 10 and 20 km) that is virtually always present. Fourier analysis of HRDI zonal and meridional wind profiles have been carried out, and the spectral characteristics are sorted by latitude, month and local time. Power spectral density (PSD) exhibits a universal exp(-km) structure in the 10-20km wavelength regime, with K lying between 2 and 3. The observed PSD for wavelengths between 10 and 20 km is a factor of 3 higher than a null spectrum constructed from HRDI reported error bars multiplied by randomly varying numbers between -1 and +1. Stokes parameters were consolidated by month into Northern and Southern hemisphere middle and high latitudes belts (40-72 degrees), tidal belts (32-16 degrees) and a tropical belt (8S-8N). Vertical waves between 10 and 15 km in wavelength are about 10-15% polarized everywhere. The inferred propagation direction in the middle and high latitude Southern hemisphere is predominantly meridional during solstice, and significantly more zonal during equinoxes. In the tropical belt, the wave orientations are nearly North-South during solstices, with a slightly higher east-west component during equinox. In the tidal belts where the background wind includes a strong meridional tidal wind, the preferred wave orientation has a significant zonal component during equinox. These findings are consistent with the interpretation of wave filtering by the background wind.

Lieberman, Ruth S.↗

NASA NTRS1982

Results from the Orbiting Astronomical Observatory. XLI - Photometry of M31 and M33

Filter photometry observations from the Orbiting Astronomical Observatory, in six wavelength bands between 1550 and 4250 A, are reported for several positions along the major and minor axes in M31 and at some positions in M33. The spatial resolution is 10 arcmin. The colors along the major axis of M31 within about 20 arcmin of the nucleus are similar to the nucleus; at 45 arcmin the colors are much bluer, indicating a spiral arm population. This behavior is expected from other observations of H I, CO, and early-type stars, which are all concentrated in a ring of material whose inner edge is about 7 kpc (= 35 arcmin) from the nucleus. The center of M33 is bluer than the ring of M31, and the arms of M33 even bluer. There are no colors too blue to be explained (within observational error) by unreddened B stars.

Davis, J.↗

NASA NTRS2021

Estimation of Snow Mass Information via Assimilation of C-Band Synthetic Aperture Radar Backscatter Observations Into an Advanced and Surface Model

This study assimilated Sentinel-1 C-band backscatter observations over snow-covered terrain into the Noah-Multiparameterization land surface model using support vector machine (SVM) regression and an ensemble Kalman filter to improve the modeled terrestrial snow mass estimates. The data assimilation (DA) experiment was conducted across Western Colorado from September 2016 to August 2017. As part of the DA experiments, the impact of a rule-based update was evaluated by comparing snow water equivalent (SWE) estimates via DA (with [ DAv1 ] and without [ DAv2 ] the rule-based update) against SNOTEL SWE measurements. Results confirmed that rule-based update helped minimize SVM controllability issues, and in turn, improved the accuracy of SWE estimates relative to both open loop (OL) and DAv2 . Comparison of SWE estimates from Sentinel-1 DAv1 against SNOTEL SWE revealed that 75% of stations showed improvements in bias and correlation coefficient relative to the OL. Assimilated SWE estimates also showed statistical improvements during both the snow accumulation and snow ablation periods. However, unbiased root mean square error showed a slight increase during the snow ablation period due to the large variability in the electromagnetic response of C-band backscatter over deep and/or wet snow. Improvement of the SWE estimates also resulted in improving river discharge estimates compared to in situ measurements. River discharge using Sentinel-1 DAv1 improved the Nash–Sutcliffe efficiency at all available stations. These results suggest that physically constrained SVM can serve as an efficient observation operator for snow mass DA through explicit consideration of the first-order C-band scattering mechanisms over different terrestrial snow conditions.

Jongmin Park↗

NASA NTRS1999

The Distribution of Titanium in Lunar Soils on the Basis of Sensor and In Situ Data Fusion

A variety of remote-sensing measurements have been used to map the distribution of elements on the Moon as a means of providing constraints on the processes from which its crust and major terranes originated. Discussed here is Ti, which is incorporated into refractory minerals such as ilmenite during the latter stages of differentiation, and is thus a most useful element for understanding mare basalt petrogenesis. One of the earliest Ti maps showed Ti variations in nearside maria on the basis of groundbased spectral reflectance measurements. A map of Ti derived from gamma-ray measurements on Apollo 15 and 16 was produced at about the same time, and was improved upon considerably by Davis and coworkers, who effectively removed sources of spurious variation from Fe and Al or REE (e.g., Th) interference, and calibrated Ti on the bases of landing-site soil averages. In recent years, spectral reflectance measurements from Clementine have been used by Lucey and coworkers to produce global Ti distribution maps as well. As we indicated previously, the Lucey and Davis maps agree to first order. Meanwhile, we are using the concept of sensor data fusion to combine measurements from the AGR (Apollo gamma-ray) and CSR (Clementine Spectral Reflectance) techniques with ground truth from lunar soils to utilize the differences between the two maps to understand the distribution of Ti within lunar soil components, as we have done with Fe. This technique should be verified and applied on Lunar Prospector gamma-ray measurements of Ti, as the calibrated data become available within the next couple of years. Lunar Ti is found principally in the mineral ilmenite, and is associated with certain components of lunar soil: crystalline Ilmenite mineral fragments and high Ti-bearing glass. All data indicate that Ti is associated with maria and mafic minerals. In AGR and CSR datasets, Ti is highest on the nearside and in the maria, particularly in southern Serenitatis/northern Tranquillitatis and northwestern Procellarum. Unlike CSR-derived values, AGR Ti values show modest increase (up to 0.7%) on the northern farside. Both techniques show a primarily unimodal distribution with a shoulder in the high Ti direction and a primary mode at approximately 0.2% (0.15% for CSR and 0.25% for AGR data), but the AGR Ti data have more structure in the shoulder, including an apparent minor mode at 1.2% Ti, representing non-highland areas. Recalibrating the two datasets on the basis of matching the peaks and range on the histograms would not account for this additional structure in the AGR data. Low Ti areas, which occur predominantly on the farside, are not represented well in the lunar sample collection, although several meteorites thought to be of lunar highland origin show Ti abundances averaging about 0.2%. (Error bars on AGR Ti values are about 0.5%.) Gamma-ray measurements reflect intrinsic Ti surface composition, regardless of the physical or chemical state of Ti. Lucey and coworkers have attempted to produce an equivalent bulk Ti map by normalizing the spectral feature at 415 nm to remove the effects of physical variations from soil to soil. The normalization is based on laboratory measurements of available samples, primarily nearside maria, and is undoubtedly optimized for the proportions in which Ti- components are found in these soils. In particular, mare soils tend to have a much higher proportion of opaque mineral grains, which appears to be well correlated with Ti abundance. The relationship between Ti abundance and proportion of Ti-bearing glass (weighted for Ti abundance) is not as direct. Also shows the relationship between landing site soil average, AGR, and CSR Ti abundances (with CSR filtered to match resolution and field of view of AGR data). AGR Ti values show the best agreement with Ti soil averages, CSR Ti values being somewhat higher in high-Ti soils, and lower in low-Ti soils. These results would tend to support the argument that the technique for deriving Ti from the CSR measurements is optimized for one of the primary Ti-bearing components, most likely opaque mineral grains, and is thus over-compensating and reporting higher Ti values in mare areas, and underreporting Ti in highland areas with much lower proportions of opaque mineral grains. More information is contained in original.

Clark, P. E.↗

DOE OSTI2023

AutoSourceID-Classifier: Star-galaxy classification using a convolutional neural network with spatial information

Aims.Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification’s reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. Methods.The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results. Results.We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C’s direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.

Astronomy & Astrophysics↗