Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coefficient of variation”

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 OSTI2025

Evaluation of historical precipitation interannual variability in CMIP6 over the United States

Interannual precipitation variability profoundly influences society via its effects on agriculture, water resources, infrastructure, and disaster risks. In this study, we use daily in situ precipitation observations from the global historical climatology network-daily (GHCN-D) to assess the ability of 21 Coupled Model Intercomparison Project Phase 6 (CMIP6) models, including the 50-member fifth-generation Canadian Earth System Model single model initial-condition large ensemble (CanESM5_SMILE), to realistically simulate historical interannual precipitation variability trends within 17 regions of the contiguous United States (CONUS). We assess how accurately the CMIP6 simulations align with observational data across annual, summer, and winter periods, focusing on four key hydrometeorological metrics, including interannual precipitation variability, relative interannual precipitation variability (coefficient of variation), annual mean precipitation, and annual wet day frequency. Our findings reveal that CMIP6 ensemble members generally reproduce the spatial patterns of observed trends in annual mean precipitation. In most regions, models agree well with the signs of observed changes in annual mean precipitation, though discrepancies in trend magnitude are evident. Further, observed trends in winter mean precipitation broadly exhibit a spatial pattern similar to that of the observed annual mean. However, analysis of the CanESM5_SMILE shows that trends in precipitation variability may primarily be the result of model-simulated internal variability, suggesting caution in interpreting multi-model single-realization ensemble results. Challenges in accurately simulating interannual precipitation variability underscore the need for ongoing model refinement and validation to enhance climate projections, especially in regions vulnerable to extreme precipitation events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings

This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

CHARACTERISTICS OF CYCLE-TO-CYCLE COMBUSTION VARIABILITY AT PARTIAL-BURN LIMITED AND MISFIRE LIMITED SPARK TIMING UNDER HIGHLY DILUTED CONDITIONS

Diluted combustion with exhaust gas recirculation (EGR) has been widely employed to improve the fuel economy of spark ignition engines. The combustion kinetics, however, are affected and the flame propagation speed is decreased. In order to compensate for this adverse effect, the spark timing needs to be recalibrated to achieve maximum brake torque (MBT). At high levels of EGR dilution, the spark timing is constrained by two ignition limits: 1) the partial-bum limit where the spark timing is retarded from MBT and 2) the misfire limit where the spark timing is too advanced. This work uses a probabilistic framework to capture the differences between both ignition limits. In particular, it introduces the concept of a nominal indicated mean effective pressure (IMEP) distribution based on the stochastic properties of the cycle-to-cycle variability (CCV) at nominal stable conditions. By defining a nominal band where fully burned cycles occur with high probability, we introduce a cycle classification method that can be used to 1) determine the level of randomness of misfire and partial-burn events, and 2) measure CCV. The new CCV metric based on the density of the nominal band is compared with the traditional coefficient of variation of IMEP (CoVIMEP). It is shown that the nominal band concept, together with the CoVIMEP, can help to discern between partial-burn limited and misfire limited conditions. Furthermore, the Kullback-Leibler divergence is used to demonstrate that the IMEP distribution is significantly different between nominal and partial-burn/misfire limited conditions. Experiments are carried at various EGR levels and spark timings while recoding in-cylinder pressure at steady state. Although the emphasis of this work is to characterize the differences of both ignition limits from a probabilistic point of view, similarities between partial-burn cycles at either limiting conditions are also discussed.

Maldonado, Bryan↗
DOE OSTI2022

EDD Basic Stats and Graphs Notebook analysis (EDD BSG Notebook) v1.0

This jupyter notebook calculates basic statistics (e.g., mean, standard deviation, coefficient of variation) and simple graphs (e.g., bar graphs, line plots) for data from the Experiment Data Depot (EDD) to provide rapid and reproducible assessment of data quality to aid research efforts across the JBEI and ABF projects. It rapidly and reproducibly calculates basic statistical values for data stored in the EDD which aids researchers and strengthens comparisons across different experiments and projects.

Petzold, ChristopherJ↗
DOE OSTI2020

Cyclic dynamics of misfires and partial burns in a dilute spark-ignition engine

Here, this study investigates the cyclic dynamics of cumulative heat release data past the edge of stability in a dilute spark-ignition engine. Emphasis is placed on analyzing the cyclic dynamics near the dilute limit where partial burns and misfires are frequent. These events are often followed by a higher-energy cycle due to the feed-forward mechanism present in the residual gases. These patterns are deterministic and increase the coefficient of variation to undesirable levels. Symbol sequence analysis was used to investigate the cyclic dynamics of these low–high patterns. The heat release was partitioned on an energy basis to give physical meaning to each partition and each sequence created when analyzing the symbol sequence results. This partitioning method provided insight into the differences in the dynamics when operating in the misfire or partial burn regime. These differences could impact the control method used.

42 ENGINEERING↗
DOE OSTI2024

Genome-wide Association Study Identifies Candidate Loci with Major Contributions to the Genetic Control of Pod Morphological Traits in Snap Bean

Snap beans are cultivars of common bean (Phaseolus vulgaris) that are cultivated for their fleshy immature pods that exhibit a wide diversity of pod shapes and sizes. The genetic basis of the snap bean pod shape is complex and involves the interaction of multiple genes. This study used a snap bean diversity panel composed of heirloom and improved cultivars used in North America and genome-wide association studies (GWAS) to investigate the genetic basis of pod morphological characteristics, including length, width, height, width/height ratio, and coefficients of variation (CVs). The GWAS detected multiple genomic regions associated with each pod trait, with a total of 20 quantitative trait loci (QTLs) for pod length, 9 for pod width, 14 for pod height, and 10 for pod width/height ratio. Regarding theCVof each pod trait, genome-wide association analyses detected six QTL for lengthCVs, five for widthCVs, 15 for heightCVs, and six for width/height ratioCVs. Thirteen regions in seven chromosomes were associated with two or more pod traits. Eighteen QTLs for pod traits in this study colocated with previously reported QTLs for pod and seed traits. The QTL intervals encompass gene models with homologues in other species that are involved in the control of developmental processes. These results capture the complex nature of the genetic control of snap bean pod traits and confirm the significance of genomic regions harboring overlapping QTLs identified in this and other studies. The phenotypic expression of pod traits in snap bean appears to be under the control of a few genomic regions with a strong effect with additional contributions of multiple small-effect regions. Validation of the function of the candidate genes identified in associated regions will contribute to our understanding of legume pod development.

Agriculture↗
DOE OSTI2021

Results of the 2018 Wood Stove Design Challenge

The 2018 Wood Stove Design Challenge was an international design competition which sought to identify top performing residential wood stoves based on automation. To prepare for the event, the Northeast States for Coordinated Air Use Management (NESCAUM) developed a testing protocol that challenged the stoves by testing them in more field-like conditions—capturing emissions from start-up, reloading a stove, and the use of larger piece sizes. Flue gas emissions were measured using a combination of in-stack and novel dilution sampling methods, as the event occurred on the National Mall in Washington D.C. in a non-laboratory setting. Particulate matter (PM), carbon monoxide (CO), carbon dioxide (CO2), and methane (CH4) were measured from three stoves in real-time. A combustion efficiency was also calculated for each stove in real-time. Measured PM emission rates ranged between 1.8 and 8.0 g/hr for all stoves and operating conditions, with the highest emissions being measured during cold start in all cases. The test average PM emission rates were 2.4, 4.0, and 2.4 g/hr for stoves A, B, and C, respectively. Results showed emissions measured during transient operations, which are often excluded in current certification methods, can be significantly higher than steady-state periods—echoing other studies and highlighting the importance of testing in various operational modes for more realistic emission estimates. Additionally, the stoves had overall CO emission rates of 48.8, 284.8, and 100.8 g/hr, for stoves A, B and C, respectively. Calculated PM and CO emission factors overall were low considering the testing protocol sought to follow more challenging yet realistic practices in terms of fuel loading and operating procedures. The overall estimated combustion efficiency for stoves A, B, and C was 85%, 78%, 76%, respectively and in all cases was lowest during the cold start period.This report details the successfully proposed and assembled instrumentation that was relatively portable for field-site testing and the results of the emissions measurements for the 2018 Wood Stove Design Challenge competition. Overall, the repeatability of each stove was favorable with coefficients of variation (COV) ranging from 7 to 15% for measured PM concentrations.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Mach Conference, 2022.

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiC/SiC) composite tubes due to their superior mechanical properties. The SiC/SiC tubes are known to remain stable even after prolonged exposure to radiation. However, there are many variables in the manufacturing of SiC/SiC composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. The intention is to identify the distribution of elastic constants of SiC/SiC tubes which would help model the performance of the novel cladding material. To do so, a sensitivity analysis will be carried out to determine the dominant variables that influence the elastic constants of SiC/SiC tubes. A chosen number of combinations of dominant variables will be sampled through the Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in the circumferential and longitudinal direction, Poisson's ratio, and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used as a surrogate to FE analysis. Once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The uncertainty can be quantified by calculating the coefficient of variation (CV) based on MCS-generated data which will be presented at the conference.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Characterization of B255 R183 C-Cell Cobalt-60 Irradiator Pool

Characterization of the Radiation Calibration and Spectroscopy Laboratory (RCS) Co-60 source radiation pool was successfully completed. The characterization was done by performing several measurements of dose rate in the center of the Co-60 source radiation pool and at multiple heights from the bottom of the pool. The measurements were taken between November 2018 through January 2019. Opti-Chromic FWT-70-40M dosimeters were used for this measurement. A 7.3 cm inner diameter stainless steel cylinder (vessel) was used to lower the dosimeters to the bottom of the pool inside the 16 cm dimeter Co-60 source basket. The dose rate was measured at heights of 5.5, 12, 19, and 23 cm from the bottom of the pool. The report “Tour of the Standards and Calibration Laboratory” includes information about the Co-60 source radiation pool. The importance of performing the Co-60 source characterization is to confirm previously reported Co-60 source dose rate measurements and to provide a documented measuring procedure for researchers and others who want to use the Co-60 source to evaluate the effects of radiation on space-based equipment for other applications. The measured dose rate at the abovementioned heights was 0.65, 0.69, 0.53, 0.23 krad/min respectively with the highest coefficient of variation of less than 10%.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Fast Semi-automated Filtration Method for Non-targeted LC-QTOF Data of Aged Nitroplasticizer Samples

A full dataset of aged nitroplasticizer (NP) is composed of more than 2000 unique mass-to-charges (m/z) when combining the non-targeted data obtained from both positive and negative electrospray ionization modes in time-of-flight mass spectrometry. Therefore, manual processing of these data often takes days, weeks, or even months to scrutinize for mechanistic insights. To effectively extract meaningful signals that represent vital degradation intermediates in the early NP degradation mechanism, a semi-automated postprocessing workflow for data filtering, tailored to the aging experiment of NP, has been developed. The automated portion of this workflow is written in a Python code (using pandas, numpy, and matplotlib libraries), which removes more than 65% of potential false signals within seconds via four threshold-based adjustable filters: signal sensitivity, coefficient of variation, number of measurements, and retention time variability. As for the manual portion, a pattern-based inspection method is employed to reduce another 23% or more false positives, which greatly simplifies data visualization and results in less than 3% of potential candidate m/z needing in-depth data interpretation. As a positive control, known compounds are verified. Using this semi-automated data reduction method, the amount of time required is reduced to a matter of hours for data filtering in the non-targeted datasets of aged NP, which saves more time and effort for compound identification.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Annual Performance Testing of Tracer Gas Detectors for use in Radionuclide NESHAP Compliance Testing

Los Alamos National Laboratory’s (LANL’s) compliance with Radionuclide NESHAP0F 1 regulations is managed by the Radioactive Air Emissions Management (RAEM) team, part of LANL’s Compliance Programs group (EPC-CP). One area of the Radionuclide NESHAP addresses requirements for siting a stack sample system. Prior to commissioning a new stack sampling system, the ANSI Standard for stack sampling requires that the stack sample location must meet several criteria, including uniform mixing of tracer gas (sulfur hexafluoride, SF 6 ) and tracer aerosol (liquid oil droplets) in the air stream. This memorandum will specifically focus on the tracer gas portion of the ANSI Standard. For these mix tests, tracer gas is injected into the stack air stream and the resulting air concentrations are measured across the plane of the stack at the proposed sampling location. The coefficient of variation of the media concentration must be under 20% when evaluated over the central 2/3 area of the stack or duct.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗
DOE OSTI2023

Hanford Waste Treatment Plant LAB Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Waste Treatment Plant laboratory (LAB) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAB stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing, determines the range of conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rates along with the average and maximum qualified stack flow rates for each of the LAB facility stacks. For each stack, the operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack tests at the LAB facility. Flow angle results were primarily less than 10°, except for the LB-C2 Fan A results, which were an average of 13.7°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition averaged between 1.5 and 4.1% COV, which were all within the range of the target percent coefficient of variation values from the scale model tests. Based on these stack verification test results, the three LAB filtered exhaust stack sampling locations meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes single-fan as well as dual-fan operations for LB-C2, each of the dual-fan operating conditions for LB-S1, and each single-fan operating condition for LB-S2. Further changes to the system configuration or operating conditions that are outside the bounds described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Probabilistic Assessment and Uncertainty Analysis of CO2 Storage Capacity of the Morrow B Sandstone—Farnsworth Field Unit

This paper presents probabilistic methods to estimate the quantity of carbon dioxide (CO2) that can be stored in a mature oil reservoir and analyzes the uncertainties associated with the estimation. This work uses data from the Farnsworth Field Unit (FWU), Ochiltree County, Texas, which is currently undergoing a tertiary recovery process. The input parameters are determined from seismic, core, and fluid analyses. The results of the estimation of the CO2 storage capacity of the reservoir are presented with both expectation curve and log probability plot. The expectation curve provides a range of possible outcomes such as the P90, P50, and P10. The deterministic value is calculated as the statistical mean of the storage capacity. The coefficient of variation and the uncertainty index, P10/P90, is used to analyze the overall uncertainty of the estimations. A relative impact plot is developed to analyze the sensitivity of the input parameters towards the total uncertainty and compared with Monte Carlo. In comparison to the Monte Carlo method, the results are practically the same. The probabilistic technique presented in this paper can be applied in different geological settings as well as other engineering applications.

03 NATURAL GAS↗
DOE OSTI2024

A Comparative Experimental Analysis of Natural Gas Dual Fuel Combustion Ignited by Diesel and Poly OxyMethylene Dimethyl Ether

Dual-fuel low-temperature combustion is a possible solution for alleviating the tradeoff between oxides of nitrogen and soot emissions in conventional diesel combustion, albeit with poor combustion stability, high carbon monoxide, and unburned hydrocarbon emissions at low engine loads. The present work compares emissions and combustion (heat release and other metrics) of both diesel and poly-oxy methylene dimethyl ether as high-reactivity fuels to ignite natural gas while leveraging spray-targeted reactivity stratification, which involved multiple injections of the high-reactivity fuels. The experiments included six parametric sweeps of: (1) start of first injection, (2) start of second injection, (3) percentage of energy substitution of natural gas, (4) commanded injection duration ratio, (5) rail pressure, and (6) intake pressure. The experiments were performed on a 1.8 L heavy-duty single-cylinder research engine operating at a medium speed of 1339 rev/min. Not-to-exceed limits for the indicated oxides of nitrogen emissions, maximum pressure rise rate, and the coefficient of variation of the indicated mean effective pressure were set to 1 g/kWh, 10 bar/CAD, and 10%, respectively. The indicated emissions decreased and combustion improved significantly for both fueling combinations when the experimental procedure was applied.

Partridge, Kendyl Ryan↗
DOE OSTI2020

Comparative Analysis of Thrombin Calibration Algorithms and Correction for Thrombin-α2macroglobulin Activity

Background: The thrombin generation (TG) test is useful for characterizing global hemostasis potential, but fluorescence substrate artifacts, such as thrombin-α2macroglobulin (T-α2MG) signal, inner filter effect (IFE), substrate consumption, and calibration algorithms have been suggested as sources of intra- and inter-laboratory variance, which may limit its clinical utility. Methods: Effects of internal vs. external normalization, IFE and T-α2MG on TG curves in normal plasma supplemented with coagulation factors, thrombomodulin, and tissue factor were studied using the Calibrated Automated Thrombinography (CAT; Diagnostica Stago, Parsippany, NJ, USA) and in-house software. Results: The various calibration methods demonstrated no significant difference in producing TG curves, nor increased the robustness of the TG assay. Several TG parameters, including thrombin peak height (TPH), produced from internal linear calibration did not differ significantly from uncalibrated TG parameters. Further, TPH values from internal linear and nonlinear calibration with or without T-α2MG correction correlated well with TPH from external calibration. Higher coefficients of variation (CVs) for TPH values were observed in both platelet-free and platelet-rich plasma with added thrombomodulin. Conclusions: Our work suggests minimal differences between distinct computational approaches toward calibrating and correcting fluorescence signals into TG levels, with most samples returning similar or equivalent TPH results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Fast, Nondestructive and Precise Biomass Measurements Are Possible Using Lidar-Based Convex Hull and Voxelization Algorithms

Light detection and ranging (lidar) scanning tools are available that can make rapid digital estimations of biomass. Voxelization and convex hull are two algorithms used to calculate the volume of the scanned plant canopy, which is correlated with biomass, often the primary trait of interest. Voxelization splits the scans into regular-sized cubes, or voxels, whereas the convex hull algorithm creates a polygon mesh around the outermost points of the point cloud and calculates the volume within that mesh. In this study, digital estimates of biomass were correlated against hand-harvested biomass for field-grown corn, broom corn, and energy sorghum. Voxelization (r = 0.92) and convex hull (r = 0.95) both correlated well with plant dry biomass. Lidar data were also collected in a large breeding trial with nearly 900 genotypes of energy sorghum. In contrast to the manual harvest studies, digital biomass estimations correlated poorly with yield collected from a forage harvester for both voxel count (r = 0.32) and convex hull volume (r = 0.39). However, further analysis showed that the coefficient of variation (CV, a measure of variability) for harvester-based estimates of biomass was greater than the CV of the voxel and convex-hull-based biomass estimates, indicating that poor correlation was due to harvester imprecision, not digital estimations. Overall, results indicate that the lidar-based digital biomass estimates presented here are comparable or more precise than current approaches.

Environmental Sciences & Ecology↗
DOE OSTI2021

Statistical Characterization of Temperature and Pressure Vertical Profiles for the Analysis of Laser Heterodyne Radiometry Data

The statistical analysis of historic pressure and temperature profiles from radiosonde launches for use in the fitting of molecular oxygen line shapes is presented. As the O 2 mixing ratio is nearly constant throughout the lower atmosphere, only variations in pressure and temperature profiles will affect the fit of observed O 2 features in Laser Heterodyne Radiometry (LHR) spectra. Radiosonde temperature and pressure data are extracted from the Integrated Global Radiosonde Archive (IGRA) for a given station, date, and launch time. Data may be extracted for a single launch, for the same date over several years, and/or within a window centered on a target date. The temperature and pressure profiles are further characterized by the statistical variation in coefficients of polynomial fits in altitude. The properties of the probability distributions for each coefficient are used to constrain fits of O 2 line shapes through Nelder–Mead optimization. The refined temperature and pressure profiles are then used in the retrieval of vertically resolved mixing ratios for greenhouse gases (GHGs) measured in the same instrument. In continuous collections, each vertical profile determination may be treated as a Bayesian prior to inform subsequent measurements and provide an estimate of uncertainties.

47 OTHER INSTRUMENTATION↗