Search NASA⌕ Search

SEARCH · Search NASA

Results for “Centroid”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗
DOE OSTI2023

Photometry, Centroid and Point-spread Function Measurements in the LSST Camera Focal Plane Using Artificial Stars

Abstract The Vera C. Rubin Observatory’s LSST Camera (LSSTCam) pixel response has been characterized using laboratory measurements with a grid of artificial stars. We quantify the contributions to photometry, centroid, point-spread function size, and shape measurement errors due to small anomalies in the LSSTCam CCDs. The main sources of those anomalies are quantum efficiency variations and pixel area variations induced by the amplifier segmentation boundaries and “tree-rings”—circular variations in silicon doping concentration. This laboratory study using artificial stars projected on the sensors shows overall small effects. The residual effects on point-spread function (PSF) size and shape are below 0.1%, meeting the ten-year LSST survey science requirements. However, the CCD mid-line presents distortions that can have a moderate impact on PSF measurements. This feature can be avoided by masking the affected regions. Effects of tree-rings are observed on centroids and PSFs of the artificial stars and the nature of the effect is confirmed by a study of the flat-field response. Nevertheless, further studies of the full-focal plane with stellar data should more completely probe variations and might reveal new features, e.g., wavelength-dependent effects. The results of this study can be used as a guide for the on-sky operation of LSSTCam.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Implementation of a High-Speed Multichannel Data Acquisition System for Magnetic Diagnostics and Plasma Centroid Position Control in ISTTOK

In tokamak and other fusion devices, magnetic control is the main tool that allows to regulate the plasma current, position and shape; it is in charge of actuating the desired plasma current waveform, steering the plasma position to a given set point and maintain the plasma shape close to a prescribed plasma equilibrium. This work describes the application of several physics concepts and computational tools in order to obtain a novel optimal controller for the plasma centroid position, which has been implemented and tested in the real-time plasma control system at the ISTTOK tokamak. A key point for the development of the new control system was the installation of a recently upgraded hardware, that numerically integrates in real-time the magnetic probes signals.

poloidal field coils↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗
DOE OSTI2023

The Zonal Patterns in Late Quaternary Tropical South American Precipitation

Speleothem oxygen isotope records (δ 18 O) of tropical South American rainfall in the late Quaternary show a zonal “South American Precipitation Dipole” (SAPD). The dipole is characterized by opposing east-west precipitation anomalies compared to the present—wetter in the east and drier in the west at the mid-Holocene (~7 ka), and drier in the east and wetter in the west at the Last Glacial Maximum (~21 ka). However, the SAPD remains enigmatic because it is expressed differently in western versus eastern δ18O records and isotope-enabled climate model simulations usually misrepresent the magnitude and/or spatial pattern of δ 18 O change. Here, we address the SAPD enigma in two parts. First, we re-interpret the δ 18 O data to account for upwind rainout effects that are known to be pervasive in tropical South America, but are not always considered in Quaternary paleoclimate studies. Our revised interpretation reconciles the δ 18 O data with cave infiltration and other proxy records, and indicates that the centroid of tropical South American rainfall has migrated zonally over time. Second, using an energy balance model of tropical atmospheric circulation, we hypothesize that zonal migration of the precipitation centroid can be explained by regional energy budget shifts, such as changing Saharan albedo associated with the African Humid Period, that have not been modeled in previous SAPD studies. This hypothesis of a migrating precipitation centroid presents a new framework for interpreting δ18O records from tropical South America and may help explain the zonal rainfall anomalies that predate the late Quaternary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Gated truncated readout system

A gated truncated readout system for position sensitive or imaging detectors that improves resolution over traditional readout systems. The readout system includes two or more amplifiers that receive a multichannel output analog data from the detector. Analog gates control circuitry, included in the readout circuit, receives the signals from the amplifiers, determines a fractional value of the sum-integral of the signals, and enables analog gates operation around an area of interest, disabling all other channels where noise dominates the signal value and thereby improving interpolation accuracy of the signals centroid position and the detector resolution. Filtered signals are transmitted to a centroid interpolation signal processing device for computation of the centroid position. As a result disabling all channels where noise dominates the signal value, the gated truncated readout system provides better accuracy improved detector resolution.

Popov, Vladimir↗
DOE OSTI2024

Fe Kα XANES, Fe Kβ HERFD XANES and EPMA flank method determinations of the oxidation state of Fe in garnet

The ferric to total iron ratios (Fe 3+ /ΣFe) of garnets can be paired with thermodynamic mineral activity models to quantify the oxygen fugacity of garnet-bearing rocks. However, techniques with a high analytical and spatial resolution are necessary to distinguish differences in garnet Fe 3+ /ΣFe ratios at the percent level and to accurately measure garnets that are zoned or contain inclusions. We acquired conventional Fe Kα and high-resolution energy fluorescence detection (HERFD) Fe Kβ X-ray absorption near edge structure (XANES) spectra and electron microprobe flank method analyses on a suite of 27 peridotitic and eclogitic garnets with Fe 3+ /ΣFe ratios previously determined by Mössbauer spectroscopy to evaluate the precision of each technique. We examined variations in the energy and intensity of three XANES spectral features as a function of Fe 3+ /ΣFe ratios: 1) the intensity ratio of two-post edge features (I-ratio; Fe Kα only); 2) the energy of the Fe edge at 90% normalized intensity (E 0.9 ; Fe Kα only) and 3) the pre-edge centroid energy (Fe Kα and HERFD Fe Kβ). In accordance with previous work, we find the energies of garnet pre-edge centroids are relatively insensitive to Fe 3+ /ΣFe ratios. The I-ratios of peridotitic and eclogitic garnets are offset from each other at low Fe 3+ /ΣFe ratios (≤0.13); I-ratio garnet XANES calibrations are composition-specific. The E 0.9 feature is independent of garnet major element composition in spectra that have been corrected for the effects of self-absorption. We produce two Fe Kα garnet XANES calibrations based on variations in the E0.9 feature; one calibration with all garnet reference materials included (Fe 3+ /ΣFe up to 1.0; “all garnet calibration”) and another calibration specific to garnets with low Fe 3+ /ΣFe ratios (“low ferric calibration”). Fe 3+ /ΣFe ratios calculated from the mean of up to 25 flank method measurements on eight garnet reference materials fall within 4% absolute of a one-to-one correlation with Fe 3+ /ΣFe ratios measured by Mössbauer. The standard error of the mean Fe 3+ /ΣFe ratio calculated from flank method approaches the Mössbauer-determined Fe 3+ /ΣFe ratio within estimated error (3%) after three analyses. Flank method precision is enhanced at higher beam current; however, the precision of the flank method does not approach the precision of XANES under any microprobe analytical condition tested here. Garnet reference materials detailed here are available by request to the Smithsonian Institution.

58 GEOSCIENCES↗
DOE OSTI2026

Predicting U 3 O 8 powder processing conditions: An AI/ML approach analyzing deep learning embeddings of SEM micrographs

High-resolution SEM images of uranium-oxide powders encode micro- and nanoscale clues to their synthesis route and calcination temperature. We trained a ResNet-50 model on 11 commercial-scale U₃O₈ classes, ammonium diuranate (ADU) or uranyl peroxide (H₂O₂) precursors calcined at temperatures ranging from 400 to 750 °C and added a 256-D projection head before the classifier to analyze the learned representation. The best of eight seeds reached 92.4 % accuracy on reserved testing data, but our focus is the structure of the embedding space rather than the accuracy and labels. We quantify class relatedness in the original 256-D space using centroid similarity and distributional distances, and we use Uniform Manifold Approximation Projection (UMAP) for visualization. ‘Unknown’ images from different preparation methods, SEM operators, and from the literature localized near the expected classes under a nearest-centroid analysis without retraining, as well as clustered in similar UMAP space. In conclusion, this embedding-centered workflow complements black-box classification by providing quantitative, similarity-based comparisons of U₃O₈ morphologies and reduces storage space by up to 98 % for image data used in millisecond vector search comparisons.

36 MATERIALS SCIENCE↗
DOE OSTI2023

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Locating the missing large-scale emission in the jet of M87* with short EHT baselines

In very-long baseline interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of a large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show how these trivial closure phases become nonzero with a brightness distribution on smaller scales than their short baseline would suggest. When applied to sources that are made up of a bright compact and large-scale diffuse component, the trivial closure phases directly measure the centroid relative to the compact source and higher-order image moments. We present a technique to measure these image moments with minimal model assumptions and validate it on synthetic Event Horizon Telescope (EHT) data. We then apply this technique to 2017 and 2018 EHT observations of M87* and find a weak preference for extended emission in the direction of the large-scale jet. We also apply it to 2021 EHT data and measure the source centroid about 1 mas northwest of the compact ring, which is consistent with the jet observed at lower frequencies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Electron beam characterization via quantum coherent electric field imaging

We demonstrate an all-optical, minimally invasive electron beam (𝑒−beam) characterization using Rydberg electrometry. An 𝑒-beam passes through a localized detector volume containing thermal Rb atoms that are optically excited in a quantum superposition of ground and Rydberg states. The 𝑒-beam's electric field perturbs atomic coherence and modifies resonance fluorescence in a narrow spectral region near the two-photon optical transition. Imaging Rb fluorescence provides a map of the electric field from the 𝑒-beam via Rydberg state Stark shifts and allows us to reconstruct the 𝑒-beam width, centroid position, and current. For an 𝑒-beam of 20 keV and electron currents of 20–40 µ⁢A, we measured the centroid position with an accuracy of ≤ 8 µ⁢m and width to within 100 µ⁢m, limited by the experimental sensitivity of electric field reconstruction of ∼0.02 V/cm. This method is suitable for a wide range of electron energies and currents, and can be adopted for nuclear and high-energy experiments for minimally invasive real-time diagnostics and profiling of charged particle beams.

Behary, Rob [College of William and Mary, Williams↗
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗