Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical distributions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Bulk and position-specific isotope geochemistry of natural gases from the Late Cretaceous Eagle Ford Shale, south Texas

Here, a dataset of bulk and position-specific isotope compositions of shale gases from the Late Cretaceous Eagle Ford Shale, south Texas is reported. The chemical and bulk isotopic compositions of the seven samples in this study and those available from the literature show that across the play, produced natural gas from the Eagle Ford Shale preserves a wide range of gas compositions (wetness, ~5 to >90%) and δ 13 C values of C 1 –C 3 . The depth profiles of the isotopic compositions suggest that the western region of the Eagle Ford Shale deposit had experienced as much as 700–800 m of uplift-erosion or gas migration. The timing of and lost-gas fractions by gas expulsion events very likely affected their gas compositions and bulk isotope compositions. Some deep (>3000–3500 m), matured (%Ro > 1.5) gases with heavy δ 13 C(C 1 –C 3 ) values indicate significant loss (>50%) of the early-stage gases. For the position-specific isotope deviations of propane, the ΔC 2-1 values of the five samples show small decrease with well depths, while the ΔH 2-1 values have a general increasing trend. Although some samples fall very close to the equilibrium model trajectories with reasonable calculated temperatures (138–148 °C), the position-specific isotope compositions of propane from the Eagle Ford Shale gases are likely results of thermal cracking of various organic molecules within the source rocks, which have different activation energy for cracking and non-statistical distributions of C/H isotope within them.

58 GEOSCIENCES↗
DOE OSTI2025

Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness

A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

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

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

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

Visualization of Multifractal Superconductivity in a Two-Dimensional Transition Metal Dichalcogenide in the Weak-Disorder Regime

Eigenstate multifractality is a distinctive feature of noninteracting disordered metals close to a metal–insulator transition, whose properties are expected to extend to superconductivity. While multifractality in three dimensions (3D) only develops near the critical point for specific strong-disorder strengths, multifractality in 2D systems is expected to be observable even for weak disorder. Here we provide evidence for multifractal features in the superconducting state of an intrinsic, weakly disordered single-layer NbSe 2 by means of low-temperature scanning tunneling microscopy/spectroscopy. The superconducting gap, characterized by its width, depth, and coherence peaks’ amplitude, shows a characteristic spatial modulation coincident with the periodicity of the quasiparticle interference pattern. The strong spatial inhomogeneity of the superconducting gap width, proportional to the local order parameter in the weak-disorder regime, follows a log-normal statistical distribution as well as a power-law decay of the two-point correlation function, in agreement with our theoretical model. Furthermore, the experimental singularity spectrum f(α) shows anomalous scaling behavior typical from 2D weakly disordered systems.

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

Hyperspectral Detection of the Fluorescence Shift between Chirality-Sorted Empty and Water-Filled Single-Wall Carbon Nanotube Enantiomers

Single-wall carbon nanotubes (SWCNTs) have extraordinary electronic and optical properties that depend strongly on their exact chiral structure and their interaction with their inner and outer environment. The fluorescence (PL) of semiconducting SWCNTs, for instance, will shift depending on the molecules with which the SWCNT’s hollow core is filled. These interaction-induced shifts are challenging to resolve on the ensemble level in samples containing a mixture of different filling contents due to the relatively large inhomogeneous line width of the ensemble SWCNT PL compared to the size of these shifts. To circumvent this inhomogeneous broadening, single-tube spectroscopy and hyperspectral imaging are often applied, which until now required time-consuming statistical studies. Here, we present hyperspectral PL microscopy combined with automated SWCNT segmenting based on either principal component analysis or a convolutional neural network, capable of both spatially and spectrally resolving the PL along the length of many individual SWCNTs at the same time and automatically fitting peak positions and line widths of individual SWCNTs. The methodology is demonstrated by accurately determining the emission shifts and line widths of thousands of left- and right-handed empty and water-filled SWCNTs coated with a chiral surfactant, resulting in four statistical distributions which cannot be resolved in ensemble spectroscopy of unsorted samples. The results demonstrate a robust method to quickly probe ensemble properties with single-enantiomer spectral resolution. Moreover, it promises to be an absolute quantitative method to characterize the relative abundances of SWCNTs with different handedness or filling content in macroscopic samples, simply by counting individual species.

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

“Which Projections Do I Use?” Strategies for Climate Model Ensemble Subset Selection Based on Regional Stakeholder Needs

Climate model (or earth system model) projections are increasingly used for climate adaptation planning and impact assessments. As part of this process, many end‐users evaluate a subset of downscaled climate projections without being aware of the implications of downscaling methodology for statistics or event outcomes. Approaches for determining a subset of global climate models to use often focus on values from the raw models, rather than from their downscaled counterparts, in other words assuming that the statistical distribution of the multi‐model ensemble does not change post downscaling. This study demonstrates that a downscaled ensemble will typically retain the change distribution as a raw ensemble, but individual models can differ dramatically post‐downscaling. We recommend that subset‐selection methods account for this possibility and that decision‐relevant downscaled climate projections provide proper descriptions of fitness‐for‐purpose and essential caveats, so that non‐specialists can interpret the results with an appropriate level of confidence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessing the potential of composite confining systems for secure and long-term CO 2 retention in geosequestration

A potential geologic target for CO 2 storage should ensure secure containment of injected CO 2 . Traditionally, this objective has been achieved by targeting reservoirs with overlying seals-regionally extensive, low permeability units that have been proven capable of retaining buoyant fluid accumulations over geologic time. However, considering that the amount of CO 2 is limited by a decadal injection period, vertical migration of CO 2 can be effectively halted by a composite system of discontinuous shale/silt/mudstone barriers in bedded sedimentary rocks. Here, we studied the impact of depositional architectures in a composite confining system on CO 2 migration and confinement at reservoir scale. We stochastically generated lithologically heterogeneous reservoir models containing discontinuous barriers consistent with statistical distributions of net-sand-to-gross-shale ratio (NTG) and horizontal correlation lengths derived from well log data and observations of producing hydrocarbon fields in Southern Louisiana. We then performed an extensive suite of reservoir simulations of CO 2 injection and post-injection to evaluate the sensitivity of CO 2 plume migration and pressure response of the composite system to a series of geologic and fluid parameters including the lateral continuity of barriers, NTG, permeability anisotropy within the sand body, and capillary pressure contrast between the sand and shale facies. The results indicate that lateral continuity of barriers and NTG are the dominant controls on CO 2 plume geometry and pressure build-up in the reservoir, while the impact of NTG is particularly pronounced. The significance of intraformational barriers becomes apparent as they facilitate the local capillary trapping of CO 2 . Those barriers improve the pore space occupancy by promoting a more dispersed shape of the plume and ultimately retard the buoyancy-driven upward migration of the plume post injection.

58 GEOSCIENCES↗
DOE OSTI2025

X-Ray Microtomography Measurements of Coevolving Particle Size and Shape

This study measures simultaneous changes in particle size and shape during grain crushing and proposes quantitative metrics to track their coevolving statistical distribution. Idealized granular materials with widely different initial particle shapes (i.e., spheres, quasi-ellipsoids, plates, and rods) are compressed oedometrically at vertical stress levels sufficient to induce particle breakage. X-ray microtomography and digital image analysis are then used to track changes in particle-scale characteristics. The observations indicate that more irregular particles tend to exhibit a higher degree of breakage compared to more spherical grains subjected to an equivalent stress level. Despite these differences, the measurements also indicate that all the tested particle sets approach a common attractor (expressed in terms of aspect ratio, flatness, and elongation), regardless of their initial morphology. A new shape evolution index is proposed to quantify such trends. With this new metric, it is shown that in all the tested materials, the particle shape evolves faster than size, especially for the more regular morphologies, thus implying that the shape tends to reach its attractor point before the grading reaches its ultimate particle size distribution. Furthermore, this finding underpins an ultimate stage of crushing-driven compression that occurs with particle size reduction but minimal shape changes (i.e., self-replicating shapes due to fractures).

Oedometric compression↗
DOE OSTI2023

An upper bound for extreme temperatures over midlatitude land

Heatwaves damage societies worldwide and are intensifying with global warming. Several mechanistic drivers of heatwaves, such as atmospheric blocking and soil moisture-atmosphere feedback, are well-known for their ability to raise surface air temperature. However, what limits the maximum surface air temperature in heatwaves remains unclear; this became evident during recent Northern Hemisphere heatwaves which achieved temperatures far beyond the upper tail of the observed statistical distribution. Here, we present evidence for the hypothesis that convective instability limits annual maximum surface air temperatures (TXx) over midlatitude land. We provide a theory for the corresponding upper bound of midlatitude temperatures, which accurately describes the observed relationship between temperatures at the surface and in the midtroposphere. We show that known heatwave drivers shift the position of the atmospheric state in the phase space described by the theory, changing its proximity to the upper bound. This theory suggests that the upper bound for midlatitude TXx should increase 1.9 times as fast as 500-hPa temperatures at the time and location of TXx occurrences. Using empirical 500-hPa warming, we project that the upper bound of TXx over Northern Hemisphere midlatitude land (40°N to 65°N) will increase about twice as fast as global mean surface air temperature, and TXx will increase faster than this bound over regions that dry on the hottest days.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reconstructing Hanford worker external doses from photons for epidemiology

The accurate reconstruction of external photon doses is essential for credible radiation epidemiology. This article presents the methodology used to derive dose estimates for 37 012 Hanford Site workers included in the Million Person Study. The approach employs historical dose records from the Hanford Radiation Exposure database and a previous epidemiology study. Bias correction factors specific to dosimeter type and period of use were applied and missing annual doses were estimated using a hierarchical nearby method to estimate deep dose equivalent for each worker. For early years with limited detection sensitivity, missed doses were quantified based on expected time-period-specific, low-dose statistical distributions. The revised dose estimates resulted in lower median and mean career doses than unadjusted data, while increasing the number of person-years with nonzero dose. Sensitivity analyses assessed the influence of bias in dosimetry measurements, missed doses and gap years on dose estimates. Differences in cumulative dose estimates between unadjusted and revised annual estimates are most prominent in the early operational years due to the highest bias during that time period.

dose reconstruction↗
DOE OSTI2024

Reliability and resilience of environmental flows under uncertainty: reconsidering water year types and inconsistent flow requirements in California

Abstract Environmental water allocation in California is a complex legal process involving various government agencies and stakeholders. E-flow requirements can be based on annual runoff typologies called water year types (WYTs), which dictate water volume, timing, and duration. In this study, we examined hydropower licensing documents of the major water and power projects in the Central Sierra Nevada to catalog e-flow requirements by WYT. In this study case, we identify how WYT classification systems and categories vary across and within different basins. Additionally, we assessed the impacts of climate change on hydrology, the frequency of WYTs identified, and the reliability and resilience of e-flows using future projections (2031–2060) of 10 Global Circulation Models (GCMs). We then propose a potential adaptation strategy using a 30 year moving percentiles approach to recalculate WYTs. We identified eight WYT classifications systems were identified, and their WYT distributions statistically significantly changes across all GCMs, even though most GCMs indicate no statistically significant change in hydrology. Disparities in future impacts are observed among and within hydropower projects, with some river reaches showing negative impacts on reliability and resilience. The adaptation strategy can generally boost resilience and improve reliability, but simply updating existing WYT thresholds without flexible regulatory frameworks reconsidering WYTs and e-flows thresholds, may not yield substantial improvements. Challenges in managing e-flows in California within regulatory and hydroclimatic contexts are intricate due to the lack of standardized approaches, leading to inconsistencies and potential conflicts among stakeholders, that will likely be exacerbated by climate change. Thus, we emphasize that targeted, site-specific, and adaptive management strategies are crucial, besides the need for a harmonized and consistent approach to defining and applying WYT categories and methods and/or e-flow assessments.

Environmental Sciences & Ecology↗
DOE OSTI2024

Uncertainty analysis of Raman spectra for measuring ortho-parahydrogen compositions

Determining ortho-parahydrogen compositions at cryogenic temperatures is an important quality control for liquid hydrogen custody exchange. High orthohydrogen compositions lead to an exothermic reaction resulting in increased boil off and increased venting losses of liquid hydrogen product by either the supplier or consumer. Traditional methods for measuring ortho-parahydrogen compositions such as hot-wire anemometry, nuclear magnetic resonance, and infrared spectroscopy are typically inconvenient at best. As the number of cryogenic hydrogen systems continue to increase, there is an increased need for flexible, standardized techniques for post processing composition measurements. A Raman spectrometer implemented in the Cryo-Catalysis Hydrogen Experimental Facility (CHEF) has demonstrated incredible flexibility for composition measurements of cryogenic hydrogen flows through a wide range of temperatures and pressures with the Raman probe in-situ at cryogenic temperatures. This article analyzes the uncertainty of Raman spectroscopy data of cryogenic hydrogen equilibrated at temperatures in IONEX catalyst. Spectra from hydrogen catalyzed to the equilibrium ortho-parahydrogen composition are processed and compared to expected statistical distributions for method verification.

08 HYDROGEN↗
DOE OSTI2025

Deep-learning-based canopy height model generation from sub-meter resolution panchromatic satellite imagery

Canopy height models (CHMs) with sufficient resolution to distinguish individual trees are useful for a variety of applications. However, standard techniques to acquire such data, such as airborne lidar surveying, are often prohibitively expensive. Deep learning techniques for generating CHMs from high-resolution imagery are an attractive option to reduce costs. To date, success with these methods has been demonstrated using multichannel aerial photography and specialized satellite data products derived from multiple sensors, neither of which is commonly available at temporal resolutions finer than one year. Here we demonstrate a method to generate sub-meter resolution CHMs in three forests in California using a more abundant data source: sub-meter resolution, panchromatic satellite imagery from a single sensor. We show that phenology and species composition play important roles in model transferability; when trained using imagery from a single conifer forest in autumn, the model performs well on autumn imagery from a second conifer forest several hundred kilometers distant with no re-training. With modest additions to the training dataset, the same model generates minimally biased estimates of canopy height in both conifer and deciduous forests during multiple seasons. Because the model operates on satellite data with global coverage and a relatively short return interval, we propose its suitability to extrapolate tree-level canopy height data to remote regions and conduct high-temporal resolution monitoring of forest structure. We furthermore demonstrate the workflow’s applicability to fire modeling by conducting simulations in forests populated by trees measured using both this approach and airborne lidar surveying. We find minimal differences in fire behavior relative to a baseline case in which only statistical distributions of tree height and crown area are known. This result underscores the value of forest structural information derived from our workflow for improving the fidelity of wildland fire simulations, among other ecological applications.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2021

Convolutional neural network identification of galaxy post-mergers in UNIONS using IllustrisTNG

ABSTRACT The Canada–France Imaging Survey (CFIS) will consist of deep, high-resolution r-band imaging over ∼5000 deg2 of the sky, representing a first-rate opportunity to identify recently merged galaxies. Because of the large number of galaxies in CFIS, we investigate the use of a convolutional neural network (CNN) for automated merger classification. Training samples of post-merger and isolated galaxy images are generated from the IllustrisTNG simulation processed with the observational realism code RealSim. The CNN’s overall classification accuracy is 88 per cent, remaining stable over a wide range of intrinsic and environmental parameters. We generate a mock galaxy survey from IllustrisTNG in order to explore the expected purity of post-merger samples identified by the CNN. Despite the CNN’s good performance in training, the intrinsic rarity of post-mergers leads to a sample that is only ∼6 per cent pure when the default decision threshold is used. We investigate trade-offs in purity and completeness with a variable decision threshold and find that we recover the statistical distribution of merger-induced star formation rate enhancements. Finally, the performance of the CNN is compared with both traditional automated methods and human classifiers. The CNN is shown to outperform Gini–M20 and asymmetry methods by an order of magnitude in post-merger sample purity on the mock survey data. Although the CNN outperforms the human classifiers on sample completeness, the purity of the post-merger sample identified by humans is frequently higher, indicating that a hybrid approach to classifications may be an effective solution to merger classifications in large surveys.

Bickley, Robert W.↗
DOE OSTI2022

Stellar multiplicity and stellar rotation: insights from APOGEE

ABSTRACT We measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterize the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds vsin i determined from our own pipeline, and distances, masses, and ages measured by Sanders & Das. We use the statistical distribution of the maximum shift in the radial velocities, ΔRVmax, as a proxy for the close binary fraction to explore the interplay between stellar evolution, rotation, and multiplicity. Assuming that the minimum orbital period allowed is the critical period for Roche Lobe overflow and rotational synchronization, we calculate theoretical upper limits on expected vsin i and ΔRVmax values. These expectations agree with the positive correlation between the maximum ΔRVmax and vsin i values observed in our sample as a function of log(g). We find that the fast rotators in our sample have a high occurrence of short-period [log(P/d) ≲ 4] companions. We also find that old, rapidly rotating main-sequence stars have larger completeness-corrected close binary fractions than their younger peers. Furthermore, rapidly rotating stars with large ΔRVmax consistently show differences of 1–10 Gyr between the predicted gyrochronological and measured isochronal ages. These results point towards a link between rapid rotation and close binarity through tidal interactions. We conclude that stellar rotation is strongly correlated with stellar multiplicity in the field, and caution should be taken in the application of gyrochronology relations to cool stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantifying concordance of correlated cosmological data sets

We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the 1% probability threshold. Furthermore, this shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

79 ASTRONOMY AND ASTROPHYSICS↗