SEARCH · Engineering Papers

Results for “pair distribution function analysis”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Studies of BONuS12 Radial GEM Detector and TCS Beam Spin Asymmetry in CLAS12

The Barely Offshell Nucleon Structure (BONuS12) experiment adopted the concept of spectator tagging technique to study the nearly-free neutron structure function Fn2 in the CLAS12 of Jefferson Lab. A novel Radial Time Projection Chamber (RTPC) detector was built, tested and integrated into the CLAS12 system to detect a back-moving low momentum tagged proton in d(e,ep)X deep-inelastic scattering. It was a 40 cm long gaseous detector consisting of 3 layers of cylindrical GEM foils for the charge amplification, with the data readout directly from the surrounding padboard. The RTPC detected the recoiling spectator proton, in coincidence with the scattered electron in the CLAS12. Nucleon structure functions are directly related to the partonic functions, quarks momentum distribution in one dimension. A Generalized Parton Distribution (GPD) came to the lime-light as it encodes the information of both longitudinal momentum and transverse position of partons inside the nucleons. Factorization of hard process such as DVCS allows to access GPDs. Timelike Compton Scattering (TCS), ?p ? ??p, is another process that allows to access the GPDs. TCS is studied experimentally in the CLAS12 of Jefferson lab using the quasi-real photoproduction of time-like photon which eventually decays to lepton pair. This dissertation presents the concept of spectator tagging in BONuS12, and the research and development efforts during the BONuS12 preparation leading up to the successful data-taking during spring and summer 2020. In addition, analysis framework to extract the beam spin asymmetry of TCS events through the CLAS12 Run group A data is presented.

Poudel, Jiwan↗

DOE OSTI2024

Observation of $t\overline{t}$ production in the lepton+jets and dilepton channels in p+Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 8.16 TeV with the ATLAS detector

This paper reports the observation of top-quark pair production in proton-lead collisions in the ATLAS experiment at the Large Hadron Collider. The measurement is performed using 165 nb -1 of p+Pb data collected at $\sqrt{^sNN}$ = 8.16 TeV in 2016. Events are categorised in two analysis channels, consisting of either events with exactly one lepton (electron or muon) and at least four jets, or events with two opposite-charge leptons and at least two jets. In both channels at least one b-tagged jet is also required. Top-quark pair production is observed with a significance over five standard deviations in each channel. The top-quark pair production cross-section is measured to be σ $t\overline{t}$ =58.1 ± 2.0(stat.)$^{+4.8}_{-4.4}$(syst.) nb, with a total uncertainty of 9%. In addition, the nuclear modification factor is measured to be R pA = 1.090 ± 0.039(stat.)$^{+0.094}_{-0.087}$(syst.). The measurements are found to be in good agreement with theory predictions involving nuclear parton distribution functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Bacterial hemophilin homologs and their specific type eleven secretor proteins have conserved roles in heme capture and are diversifying as a family

Cellular life relies on enzymes that require metals, which must be acquired from extracellular sources. Bacteria utilize surface and secreted proteins to acquire such valuable nutrients from their environment. These include the cargo proteins of the type eleven secretion system (T11SS), which have been connected to host specificity, metal homeostasis, and nutritional immunity evasion. This Sec-dependent, Gram-negative secretion system is encoded by organisms throughout the phylum Proteobacteria, including human pathogens Neisseria meningitidis, Proteus mirabilis, Acinetobacter baumannii, and Haemophilus influenzae. Experimentally verified T11SS-dependent cargo include transferrin-binding protein B (TbpB), the hemophilin homologs heme receptor protein C (HrpC), hemophilin A (HphA), the immune evasion protein factor-H binding protein (fHbp), and the host symbiosis factor nematode intestinal localization protein C (NilC). Here, we examined the specificity of T11SS systems for their cognate cargo proteins using taxonomically distributed homolog pairs of T11SS and hemophilin cargo and explored the ligand binding ability of those hemophilin cargo homologs. In vivo expression in Escherichia coli of hemophilin homologs revealed that each is secreted in a specific manner by its cognate T11SS protein. Sequence analysis and structural modeling suggest that all hemophilin homologs share an N-terminal ligand-binding domain with the same topology as the ligand-binding domains of the Haemophilus haemolyticus heme binding protein (Hpl) and HphA. We term this signature feature of this group of proteins the hemophilin ligand-binding domain. Network analysis of hemophilin homologs revealed five subclusters and representatives from four of these showed variable heme-binding activities, which, combined with sequence-structure variation, suggests that hemophilins are diversifying in function.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Measurement of the top quark pole mass using $ \textrm{t}\overline{\textrm{t}} $+jet events in the dilepton final state in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A measurement of the top quark pole mass $\mathcal{m}^\text{pole}_\text{t}$ in events where a top quark-antiquark pair ($\text{t}\bar{\text{t}}$) is produced in association with at least one additional jet ($\text{t}\bar{\text{t}}$ +jet) is presented. This analysis is performed using proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the CMS experiment at the CERN LHC, corresponding to a total integrated luminosity of 36.3 fb -1 . Events with two opposite-sign leptons in the final state (e + e – , μ + μ – , e ± μ ∓ ) are analyzed. The reconstruction of the main observable and the event classification are optimized using multivariate analysis techniques based on machine learning. The production cross section is measured as a function of the inverse of the invariant mass of the $\text{t}\bar{\text{t}}$ +jet system at the parton level using a maximum likelihood unfolding. Given a reference parton distribution function (PDF), the top quark pole mass is extracted using the theoretical predictions at next-to-leading order. For the ABMP16NLO PDF, this results in $\mathcal{m}^\text{pole}_\text{t}$ = 172.93 ± 1.36 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Search for resonant production of dark quarks in the dijet final state with the ATLAS detector

This paper presents a search for a new Z' resonance decaying into a pair of dark quarks which hadronise into dark hadrons before promptly decaying back as Standard Model particles. This analysis is based on proton-proton collision data recorded at $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider between 2015 and 2018, corresponding to an integrated luminosity of 139 fb –1 . After selecting events containing large-radius jets with high track multiplicity, the invariant mass distribution of the two highest-transverse-momentum jets is scanned to look for an excess above a data-driven estimate of the Standard Model multijet background. No significant excess of events is observed and the results are thus used to set 95% confidence-level upper limits on the production cross-section times branching ratio of the Z' to dark quarks as a function of the Z' mass for various dark-quark scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Local structural mechanism for enhanced energy storage properties in heterovalent doped NaNbO 3 ceramics

In recent years, there is a growing interest for new lead-free oxides with reversible antiferroelectric (AFE)-ferroelectric (FE) phase transition for high-power energy-storage applications. NaNbO 3 -based ceramics are particularly attractive due to their easy synthesis and cost-effectiveness. In order to stabilize reversible AFE-FE phase transition, NaNbO 3 is doped with a combination of heterovalent substitutions, although the underlying structural mechanism for the same is poorly understood. Here, we investigated local and average structures of Ca/Zr doped NaNbO 3 using neutron total scattering. Here, we show that Ca/Zr doping increases the average AFE phase (Pbma) fraction, however, the material remains as a composite of both FE (P2 1 ma) and AFE regions. Analysis of local structure suggests that increase in the long-range AFE phase results from more extensive twinning of local FE regions, due to introduced charge disorder. We propose that enhanced energy-storage properties of Ca/Zr-doped NaNbO 3 arises from localized twin boundary motion between the defect-induced pinning centers.

25 ENERGY STORAGE↗

DOE OSTI2021

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Magnetic anisotropy and two-dimensional short-range chemical ordering in Ba 1−𝑥 ⁢Na 𝑥 ⁢Fe 2 ⁢As 2

A true understanding of the properties of pnictide superconductors requires the development of high-quality materials and performing measurements designed to unravel their intrinsic properties and short-range nematic correlations which are often obscured by extrinsic effects such as poor crystallinity, inhomogeneity, domain formation, and twinning. Here, in this paper, we report the systematic growth of high-quality Na-substituted BaFe 2 As 2 single crystals and their characterization using pulsed magnetic fields x-ray diffraction and x-ray diffuse scattering. Analysis of the properties and compositions of the highest-quality crystals shows that their actual Na stoichiometry is about 50-60% of the nominal content and that the targeted production of crystals with specific compositions is accessible. We derived a reliable equation to estimate the Na stoichiometry based on the measured superconducting T c of these materials. Attempting to force spin reorientation and induce tetragonality, orthorhombic Ba 1-x Na x Fe 2 As 2 single crystals subjected to out-of-plane magnetic fields up to 31.4T are found to exhibit strong in-plane magnetic anisotropy demonstrated by the insufficiency of such high fields in manipulating the relative population of their twinned domains or in suppressing the orthorhombic order. Broad x-ray diffuse-intensity rods observed at temperatures between 30 and 300 K uncover short-range structural correlations. Local structure modeling together with 3D-Δ pair-distribution function mapping of real-space interatomic vectors show that the diffuse scattering arises from in-plane short-range chemical correlations of the Ba and Na atoms coupled with short-range atomic displacements within the same plane due to an effective size difference between the two atomic species.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Local structure of zinc–indium–tin oxide films via grazing-incidence x-ray pair-distribution functions and theoretical methods

A detailed experimental and theoretical study on the local (r ≤ 4.5 Å) atomic structure of amorphous and crystalline zinc–indium–tin oxide (ZITO) thin films using grazing-incidence x-ray Pair-Distribution Functions (PDFs), ab initio Molecular Dynamics (MD), and Empirical Potential Structure Refinement (EPSR) Monte Carlo simulations is presented. High-energy synchrotron x rays, a two-dimensional detector, and different incident angles were used to probe the depth uniformity of five (ZnO) 0.15 (In 2 O 3 ) 0.70 (SnO 2 ) 0.15 films that were deposited via pulsed-laser deposition at growth temperatures (T G ) ranging from 25 to 300 °C. Films deposited at T G ≤ 150 °C were amorphous. The partially crystalline (T G = 200 °C) and fully crystalline (T G = 300 °C) films were highly textured. Both crystalline and amorphous structures were investigated using ab initio MD and EPSR Monte Carlo simulations. The density of the amorphous films determined from the experimental data agreed with MD calculations. Coordination numbers, bond lengths, and distortion for metal–oxygen and for both the edge- and corner-shared In–metal shells up to 4.5 Å obtained from PDF analysis closely agreed with MD and EPSR simulations. There is a pronounced decrease in the edge- and corner-shared In–Zn distances arising from the shorter Zn–O bond length, Zn–O tetrahedral coordination, and In–O–Zn angle in amorphous ZITO compared to its crystalline counterpart. A maximum in electrical mobility was observed for the amorphous film just before crystallization occurred. While the peak is broad, consistent with nearly unchanged overall cation–oxygen coordination in the amorphous films, ESPR results indicate that the tetrahedral coordination follows the conductivity trend.

Grazing Incidence X-ray↗

DOE OSTI2022

Information content of the differences in the charge radii of mirror nuclei

Differences in the charge radii of mirror nuclei have been recently suggested to contain information on the slope of the symmetry energy L. To test this hypothesis, we perform statistical correlation analysis using quantified energy density functionals that are consistent with our previous knowledge on global nuclear observables such as binding energies and charge radii. We conclude that the difference in charge radii between a mirror pair, $ΔR$$^{mir}_{ch}$, is an inferior isovector indicator compared to other observables, such at the neutron skin or electric dipole polarizability α D . In particular, this quantity correlates poorly with both the neutron skin and L. Here, we demonstrate that $ΔR$$^{mir}_{ch}$ is influenced by pairing correlations in the presence of low-lying proton continuum in the proton-rich mirror-partner nucleus. Considering the large theoretical uncertainties on $ΔR$$^{mir}_{ch}$, we conclude that the precise data on mirror charge radii cannot provide a stringent constraint on L.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

eDNAjoint: An R package for interpreting paired or semi‐paired environmental DNA and traditional survey data in a Bayesian framework

Abstract Environmental DNA (eDNA) sampling is increasingly used in surveys of species distribution as a potentially sensitive and efficient monitoring method. Yet access to modelling tools designed specifically for interpreting this new data type lags behind its ubiquity. While occupancy modelling software has dominated the analytical landscape for eDNA data analysis of single species, this type of model may not always be the most appropriate. The rate of eDNA detection often corresponds to species density, rather than just occupancy, and researchers often have access to observations from non‐genetic sampling methods at the same sites. To provide users access to a modelling framework designed to maximize the use of all available data, we developed an R package, eDNAjoint . The package provides an easy‐to‐use interface for fitting a ‘joint’ model that integrates data from paired or semi‐paired eDNA and traditional surveys in a Bayesian framework. The model can be used to estimate parameters like the probability of a false positive eDNA detection and mean catch rate at a site, and the package allows access to multiple model variations and Bayesian prior customization. Additional functionality can be used for model selection, summarising posteriors and comparing the relative sensitivities of the two survey methods. We demonstrate the use of eDNAjoint by fitting a variation of the model with site‐level covariates that scale the sensitivity of eDNA sampling relative to traditional sampling. The example workflow uses binary eDNA and seine count data for the endangered tidewater goby ( Eucyclogobius newberryi ) from a study by Schmelzle and Kinziger (2016). This use case includes a prior sensitivity analysis and an evaluation of the relationship between detection rates and environmental variables. eDNAjoint has the potential to greatly increase the range of users who will be able to rigorously analyse eDNA and traditional survey data in a Bayesian framework, understand if and how eDNA can improve monitoring practices, and gain confidence in the interpretability of eDNA data.

Keller, Abigail G. [Department of Environment Scie↗

DOE OSTI2020

Elevation distributed micro-climatology data in a coastal glaciated watershed

Air temperature, ground temperature and relative humidity data were collected in a longitudinal transect of the Nooksack watershed at varying elevations from 500 to 1800 m above sea level. Data were collected by anchoring sensors from trees above winter snow levels and shaded from direct solar radiation. Paired sensors were also buried 3 cm under ground near each air temperature sensor to determine snow absence or presence. Select sites included relative humidity sensors to indicate whether precipitation was occurring. Data were collected every 3-4 h from December 2015 to Sept 2018 (with ongoing collection). Code for analysis of daily mean, minimum, maximum, and temperature change with elevation (lapse rates) are available on Github (https://doi.org/10.5281/zenodo.3239539). The sensor download and intermediate data products are available on HydroShare at (http://www.hydroshare.org/resource/222e832d3df24dea9bae9bbeb6f4219d) with publicly accessible visualization available from the Nooksack Observatory at data.cuahsi.org. Hydrologic models are generally structured with a single annual average lapse rate parameter which assumes a linear temperature gradient with elevation. The daily data (2016-2018) is used as part of ongoing studies on the non-linear dynamics and temporal variability of temperature with elevation to improve assessments of watershed function and salmon habitat.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Gene expression of functionally-related genes coevolves across fungal species: detecting coevolution of gene expression using phylogenetic comparative methods

Researchers often measure changes in gene expression across conditions to better understand the shared functional roles and regulatory mechanisms of different genes. Analogous to this is comparing gene expression across species, which can improve our understanding of the evolutionary processes shaping the evolution of both individual genes and functional pathways. One area of interest is determining genes showing signals of coevolution, which can also indicate potential functional similarity, analogous to co-expression analysis often performed across conditions for a single species. However, as with any trait, comparing gene expression across species can be confounded by the non-independence of species due to shared ancestry, making standard hypothesis testing inappropriate. We compared RNA-Seq data across 18 fungal species using a multivariate Brownian Motion phylogenetic comparative method (PCM), which allowed us to quantify coevolution between protein pairs while directly accounting for the shared ancestry of the species. Our work indicates proteins which physically-interact show stronger signals of coevolution than randomly-generated pairs. Interactions with stronger empirical and computational evidence also showing stronger signals of coevolution. We examined the effects of number of protein interactions and gene expression levels on coevolution, finding both factors are overall poor predictors of the strength of coevolution between a protein pair. Simulations further demonstrate the potential issues of analyzing gene expression coevolution without accounting for shared ancestry in a standard hypothesis testing framework. Furthermore, our simulations indicate the use of a randomly-generated null distribution as a means of determining statistical significance for detecting coevolving genes with phylogenetically-uncorrected correlations, as has previously been done, is less accurate than PCMs, although is a significant improvement over standard hypothesis testing. These methods are further improved by using a phylogenetically-corrected correlation metric. Our work highlights potential benefits of using PCMs to detect gene expression coevolution from high-throughput omics scale data. This framework can be built upon to investigate other evolutionary hypotheses, such as changes in transcription regulatory mechanisms across species.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Multimodal Analysis of Reaction Pathways of Cathode Materials for Lithium Ion Batteries

Conversion mechanism in lithium ion batteries provides higher capacity than intercalation mechanism since multiple numbers of electrons and lithium ions are associated. However, poor cycling stability, large voltage hysteresis, and low energy efficiency have been great challenges of conversion reaction. To address those issues, conversion-type electrode materials have been reformed via doping or substituting other elements. For example, iron fluorides (FeF 2 , FeF 3 ) have modified as iron oxyfluorides (FeF 1-x O x ), showing enhanced long-term stability. Furthermore, co-substituted (both anion and cation substituted) Fe 0.9 Co 0.1 OF (FeCoOF) was demonstrated excellent cycling stability. (capacity of 350 mAh g -1 at a current of 500 mA g -1 for 1000 cycles). Substituting anion and cation in iron fluoride has been suggested as an effective method to achieve better reversibility but understanding of lithiation reactions in co-substituted FeCoOF is not clear. This work takes advantage of ex-situ/ in-situ synchrotron X-ray based techniques and transmission electron microscopy to elucidate structural changes with lithium ion, which may provide fundamental insights into modifying conversion-type materials. Figure 1 presents discharge-charge curves and pair distribution function patterns acquired at each potential. As lithium ions were inserted, structural changes were noticed both at short-range and long-range. However, long-range ordering was nearly maintained even at 1 V, indicating absence of conversion reaction. Figure 2 shows lithiation induced structural evolution of a single FeCoOF nanorod observed in real time. As lithiation proceeds, the width of nanorods shows a stepwise increase, particularly A in figure 2d, which may indicate multiple steps of lithiation occur. Considering that conversion reaction takes place around 2 V in FeF 3 , co-substitution Co and O into iron fluoride may change thermodynamic features of lithiation reactions by lowering the initiation potential for conversion reaction. Instead, phase transformations occur at long-range order, which may help maintaining structural integrity during operation, eventually, achieving cycling stability.

25 ENERGY STORAGE↗

NASA NTRS1998

A Compendium of Wind Statistics and Models for the NASA Space Shuttle and Other Aerospace Vehicle Programs

The wind profile with all of its variations with respect to altitude has been, is now, and will continue to be important for aerospace vehicle design and operations. Wind profile databases and models are used for the vehicle ascent flight design for structural wind loading, flight control systems, performance analysis, and launch operations. This report presents the evolution of wind statistics and wind models from the empirical scalar wind profile model established for the Saturn Program through the development of the vector wind profile model used for the Space Shuttle design to the variations of this wind modeling concept for the X-33 program. Because wind is a vector quantity, the vector wind models use the rigorous mathematical probability properties of the multivariate normal probability distribution. When the vehicle ascent steering commands (ascent guidance) are wind biased to the wind profile measured on the day-of-launch, ascent structural wind loads are reduced and launch probability is increased. This wind load alleviation technique is recommended in the initial phase of vehicle development. The vehicle must fly through the largest load allowable versus altitude to achieve its mission. The Gumbel extreme value probability distribution is used to obtain the probability of exceeding (or not exceeding) the load allowable. The time conditional probability function is derived from the Gumbel bivariate extreme value distribution. This time conditional function is used for calculation of wind loads persistence increments using 3.5-hour Jimsphere wind pairs. These increments are used to protect the commit-to-launch decision. Other topics presented include the Shuttle Shuttle load-response to smoothed wind profiles, a new gust model, and advancements in wind profile measuring systems. From the lessons learned and knowledge gained from past vehicle programs, the development of future launch vehicles can be accelerated. However, new vehicle programs by their very nature will require specialized support for new databases and analyses for wind, atmospheric parameters (pressure, temperature, and density versus altitude), and weather. It is for this reason that project managers are encouraged to collaborate with natural environment specialists early in the conceptual design phase. Such action will give the lead time necessary to meet the natural environment design and operational requirements, and thus, reduce development costs.

Smith, O. E.↗

NASA NTRS2019

Surface erythemal UV irradiance in the continental United States derived from ground-based and OMI observations: quality assessment, trend analysis and sampling issues

Surface full-sky erythemal dose rate (EDR) from the Ozone Monitoring Instrument (OMI) at both satellite overpass time and local noon time is evaluated against ground measurements at 31 sites from the US Department of Agriculture's (USDA) UV-B Monitoring and Research Program (UVMRP) over the period of 2005–2017. We find that both OMI overpass and solar noon time EDR are highly correlated with the measured counterparts (with a linear correlation coefficient of 0.90 and 0.88, respectively). Although the comparison statistics are improved with a longer time window (0.5–1.0 h) for pairing surface and OMI measurements, both OMI overpass and local noon time EDRs have 7 % overestimation that is larger than 6 % uncertainty in the ground measurements and show different levels of dependence on solar zenith angle (SZA) and to lesser extent on cloud optical depth. The ratio of EDR between local noon and OMI overpass time is often (95 % in frequency) larger than 1 with a mean of 1.18 in the OMI product; in contrast, the same ratio from surface observation is normally distributed with 22 % of the times less than 1 and a mean of 1.38. This contrast in part reflects the deficiency in the OMI surface UV algorithm that assumes constant atmospheric conditions between overpass and noon time. The probability density functions (PDFs) for both OMI and ground measurements of noontime EDR are in statistically significant agreement, showing dual peaks at ∼20 and ∼200 mW/sq. m, respectively; the latter is lower than 220 mW/ sq. m, the value at which the PDF of daily EDR from ground measurements peaks, and this difference indicates that the largest EDR value for a given day may not often occur at local noon. Lastly, statistically significant positive trends of EDR are found in the northeastern US in OMI data, but opposite trends are found within ground-based data (regardless of sampling for either noontime or daily averages). While positive trends are consistently found between OMI and surface data for EDR over the southern Great Plains (Texas and Oklahoma), their values are within the uncertainty of ground measurements. Overall, no scientifically sound trends can be found among OMI data for aerosol total and absorbing optical depth, cloud optical depth and total ozone to explain coherently the surface UV trends revealed either by OMI or ground-based estimates; these data also cannot reconcile trend differences between the two estimates (of EDR from OMI and surface observations). Future geostationary satellites with better spatiotemporal resolution data should help overcome spatiotemporal sampling issues inherent in OMI data products and therefore improve the estimates of surface UV flux and EDR from space.

Ozone Monitoring Instrument (OMI)↗

DOE OSTI2021

Multimodal Bayesian registration of noisy functions using Hamiltonian Monte Carlo

Functional data registration is a necessary processing step for many applications. The observed data can be inherently noisy, often due to measurement error or natural process uncertainty; which most functional alignment methods cannot handle. A pair of functions can also have multiple optimal alignment solutions, which is not addressed in current literature. In this paper, a flexible Bayesian approach to functional alignment is presented, which appropriately accounts for noise in the data without any pre-smoothing required. Additionally, by running parallel MCMC chains, the method can account for multiple optimal alignments via the multi-modal posterior distribution of the warping functions. To most efficiently sample the warping functions, the approach relies on a modification of the standard Hamiltonian Monte Carlo to be well-defined on the infinite-dimensional Hilbert space. In this work, this flexible Bayesian alignment method is applied to both simulated data and real data sets to show its efficiency in handling noisy functions and successfully accounting for multiple optimal alignments in the posterior; characterizing the uncertainty surrounding the warping functions.

97 MATHEMATICS AND COMPUTING↗