Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixed-effects models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Linear Mixed-Effects Models for Human-in-the-Loop Tracking Experiment Data

Linear mixed-effects models provide several benefits over more traditional statistical inference tests that are particularly useful for most human-in-the-loop tracking experiment data. However, surprisingly, mixed models are virtually not used for the analysis of tracking experiment data. This paper uses of linear mixed-effects models to analyze combined tracking data from two previous human-in-the-loop roll tracking experiments that compared control behavior metrics collected in both a research aircraft and a motion-base simulator. In the experiments, pilots' behavior under 10 different motion configurations with varying motion filter gains and break frequencies was evaluated and compared to that in the real aircraft. The linear mixed-effects model analysis on the combined dataset confirmed the main statistical outcomes of the individual experiments. The main benefits of mixed models for this type of data were demonstrated by successfully combining data from two experiments that used different experimental conditions and of which one had an additional apparatus and the other a missing participant. Finally, the mixed-model analysis was able to explicitly test for scientifically relevant statistical differences in the dependent measures between the aircraft and simulator, as well as both experiments.

manual control↗
DOE OSTI2025

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm’s behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl .

Gorstein, Evan↗
DOE OSTI2025

The temporal onset of associations of cortical proteins with cognitive resilience vary during late life

Background: Cortical proteins associated with cognitive resilience have been identified but their temporal onset in older adults is unknown. We present a multistage approach to first identify cortical proteins associated with cognitive resilience and then examine their associated temporal onset. Methods: We used data from a subset of 1088 decedents from two cohort-studies who had selected reaction monitoring proteomics from the dorsolateral prefrontal cortex, and at least 3 cognitive assessments. Cognition was assessed using a composite derived from 19 tests. We first used linear mixed-effects models to identify cortical proteins associated with cognitive resilience. We then used functional mixed-effects models to examine non-linear associations between proteins and cognitive resilience to identify their temporal onset. Results: Mean age at death was 90 years (SD = 6.4); 69 % were female. On average, cognition started to decline at around 15 years before death, with accelerated decline in the last 7 years. We identified 40 proteins associated with cognitive resilience, of which 17 proteins also showed non-linear associations. Non-linear associations indicated that higher levels of 10 proteins were associated with slower cognitive decline between 23 and 4 years before death. In contrast, higher levels of 7 proteins were associated with faster decline only within the last 7 years before death. Conclusions: Cognitive resilience proteins are differentially related to late-life cognitive aging; the onset of proteins that maintain cognition may begin many years before the onset of proteins that hasten cognitive decline. The temporal onset of cognitive resilience proteins may be crucial for timing efficacious interventions.

Zammit, Andrea↗
NASA NTRS2017

Testing Extensions of Our Quantitative Daily of San Joaquin Wintertime Aerosols Using MAIAC and Meteorology Without Transport/Transformation Assumptions

The Western US and many regions globally present daunting difficulties in understanding and mapping PM2.5 episodes. We evaluate extensions of a method independent of source-description and transport/transformation. These regions suffer frequent few-day episodes due to shallow mixing; low satellite AOT and bright surfaces complicate the description. Nevertheless, we expect residual errors in our maps of less than 8 ug/m^3 in episodes reaching 60-100 ug/m^3; maps which detail pollution from Interstate 5. Our current success is due to use of physically meaningful functions of MODIS-MAIAC-derived AOD, afternoon mixed-layer height, and relative humidity for a basin in which the latter are correlated. A mixed-effects model then describes a daily AOT-to-PM2.5 relationship. (Note: in other published mixed-effects models, AOT contributes minimally. We seek to extend on these to develop useful estimation methods for similar situations. We evaluate existing but more spotty information on size distribution (AERONET, MISR, MAIA, CALIPSO, other remote sensing). We also describe the usefulness of an equivalent mixing depth for water vapor vs meteorological boundary layer height. Each has virtues and limitations. Finally, we begin to evaluate methods for removing the complications due to detached but polluted layers (which don't mix to the surface) using geographical, meteorological, and remotely sensed data.

Chatfield, Robert B.↗
NASA NTRS2016

Maps Suggest Transport and Source Processes of PM2.5 at 1 km x 1 km for the Whole San Joaquin Valley, Winter 2011 (Generalizations from DISCOVER-AQ)

We present interpreted data analysis using MAIAC (Multiangle implementation of Atmospheric Correction) retrievals and appropriate RAPid Update Cycle (RAP) meteorology to map respirable aerosol (PM2.5) for the period January and February, 2011. The San Joaquin Valley is one of the unhealthiest regions in the USA for PM2.5 and related morbidity. The methodology evaluated can be used for the entire moderate-resolution imaging spectrometer (MODIS, VIIRS) data record. Other difficult areas of the West: Riverside, CA, Salt Lake City, UT, and Doa Ana County, NM share similar difficulties and solutions. The maps of boundary layer depth for 1116 hr local time from RAP allows us to interpret aerosol optical thickness as a concentration of particles in a nearly well-mixed box capped by clean air. That mixing is demonstrated by DISCOVER-AQ data and afternoon samples from the airborne measurements, P3B (on-board) and B200 (HSRL2 lidar). This data and the PM2.5 gathered at the deployment sites allowed us to estimate and then evaluate consistency and daily variation of the AOT to PM2.5 relationship. Mixed-effects modeling allowed a refinement of that relation from day to day; RAP mixed layers explained the success of previous mixed-effects modeling. Compositional, size-distribution, and MODIS angle-of-regard effects seem to describe the need for residual daily correction beyond ML depth.We report on an extension method to the entire San Joaquin Valley for all days with MODIS imagery using the permanent PM2.5 stations, evaluated for representativeness. Resulting map movies show distinct sources, particularly Interstate-5 (at approx. 1km x 1km resolution) and the broader Bakersfield area. Accompanying winds suggest transport effects and variable pathways of pollution cleanout. Such estimates should allow morbiditymortality studies. They should be also useful for actual model assimilations, where composition and sources are uncertain. We conclude with a description of new work to extend these insights to similar regions, e.g. interior valleys of California, the Po Valley, the Mediterranean litoral, and the Ganges Plain.This work show generalizable use of remote sensing, a major goal of DISCOVER-AQ, Deriving Information on Surface Conditions from COlumn and VERtically Resolved Observations Relevant to Air Quality.

Chatfield, R.↗
NASA NTRS2020

Bayesian Statistical Models for Community Annoyance Survey Data

This paper demonstrates the use of two Bayesian statistical models to analyze single-event sonic boom exposure and human annoyance data from community response surveys. Each model is fit to data from a NASA pilot study.Unlike many community noise surveys, this study used a panel sample to collect multiple observations per participant instead of a single observation. Thus, a multilevel (also known as hierarchical or mixed-effects) model is used to account for the within-subject correlation in the panel sample data. This paper describes a multilevel logistic regression model and a multilevel ordinal regression model. The paper also proposes a method for calculating a summary dose-response curve from the multilevel models that represents the population. The two models’ summary dose-response curves are visually similar. However, their estimates differ when calculating the noise dose at a fixed percent highly annoyed.

Musical instruments↗
NASA NTRS2015

Spatiotemporal Prediction of Fine Particulate Matter Using High-Resolution Satellite Images in the Southeastern US 2003-2011

Numerous studies have demonstrated that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with adverse health outcomes. The use of ground monitoring stations of PM(sub 2.5) to assess personal exposure, however, induces measurement error. Land-use regression provides spatially resolved predictions but land-use terms do not vary temporally. Meanwhile, the advent of satellite-retrieved aerosol optical depth (AOD) products have made possible to predict the spatial and temporal patterns of PM(sub 2.5) exposures. In this paper, we used AOD data with other PM(sub 2.5) variables, such as meteorological variables, land-use regression, and spatial smoothing to predict daily concentrations of PM(sub 2.5) at a 1 sq km resolution of the Southeastern United States including the seven states of Georgia, North Carolina, South Carolina, Alabama, Tennessee, Mississippi, and Florida for the years from 2003 to 2011. We divided the study area into three regions and applied separate mixed-effect models to calibrate AOD using ground PM(sub 2.5) measurements and other spatiotemporal predictors. Using 10-fold cross-validation, we obtained out of sample R2 values of 0.77, 0.81, and 0.70 with the square root of the mean squared prediction errors of 2.89, 2.51, and 2.82 cu micrograms for regions 1, 2, and 3, respectively. The slopes of the relationships between predicted PM2.5 and held out measurements were approximately 1 indicating no bias between the observed and modeled PM(sub 2.5) concentrations. Predictions can be used in epidemiological studies investigating the effects of both acute and chronic exposures to PM(sub 2.5). Our model results will also extend the existing studies on PM(sub 2.5) which have mostly focused on urban areas because of the paucity of monitors in rural areas.

PM2.5↗
NASA NTRSDate not supplied

Light Improves Alertness and Mood during the Sleep Inertial Period following Slow Wave Sleep

Introduction: Waking from sleep, especially slow wave sleep (SWS), is associated with reduced alertness known as sleep inertia. Light improves alertness during sleep deprivation and circadian misalignment. In this study, we assessed the efficacy of light to improve alertness and mood immediately after waking from SWS. Methods: Twelve participants kept a sleep schedule of 8.5 h for 5 nights and 5 h for one night prior to the overnight laboratory visit (confirmed by actigraphy). Participants went to bed at their scheduled habitual bedtime in the laboratory and were monitored by standard polysomnography. After at least 5 min of SWS, participants were awoken and exposed to either red ambient light (control) or blue-enriched bright light (light) for 1 h. During this time, participants completed a subjective scale of alertness (Karolinska Sleepiness Scale, KSS) and visual analogue scales (VAS) of mood at 2 min, 17 min, 32 min, and 47 min after waking. Following this sleep inertia measurement period, all lights were turned off and participants were allowed to return to sleep. They were then awoken again from their subsequent SWS period and exposed to the opposite condition (control or light). A linear mixed-effects model with fixed effects of condition, time, and condition*time and a random effect of participant was used to determine the impact of light across the testing period. An average of baseline responses (pre-sleep) was included as a covariate. Results: Compared to the control condition, participants exposed to bright blue-enriched light reported feeling more alert (KSS: F1,77=4.955, p=.029; VASalert: F1,77=8.226, p=.005), more cheerful (VAScheerful: F1,77=8.615, p=.004), less depressed (VASdepressed: F1,77=4.649, p=.034), and less lethargic (VASlethargic: F1,77=5.652, p=.020). Discussion: Exposure to blue-enriched bright light immediately after waking from SWS may help to improve subjective alertness and mood. Future analyses will explore whether these findings extend to effects on cognitive performance.

alertness↗
NASA NTRS2011

Reliability of Upright and Supine Power Measurements Using an Inertial Load Cycle Ergometer

Practical, reliable, and time efficient methods of measuring muscular power are desirable for both research and applied testing situations. The inertial-load cycling method (ILC; Power/Cycle, Austin, TX) requires subjects to pedal as fast as possible against the inertial load of a flywheel for only 3-5 seconds, which could help reduce the time and effort required for maximal power testing. PURPOSE: 1) To test the intramachine reliability of ILC over 3 separate sessions, 2) to compare postural stance (upright vs. supine) during testing, and 3) to compare the maximal power (Pmax) output measured using ILC to that obtained from traditional isokinetic and leg press testing. METHODS: Subjects (n = 12) were tested on 4 non-consecutive days. The following tests were done on the first day of testing: isometric knee extension, isokinetic knee extension at several speeds, isokinetic power/endurance at 180/sec (Biodex System 4), leg press maximal isometric force, and leg press power/endurance. The other 3 days consisted exclusively of ILC testing. Subjects performed 6 ILC tests in an upright position and 6 ILC tests in a supine position on each day. The starting position was counterbalanced. Mixed-effects linear modeling was used to determine if any differences existed between testing days and between upright and supine for Pmax and revolutions per minute at Pmax (RPMpk). Mixed-modeling was also used to calculate intraclass correlation coefficients (ICC) to determine the reliability of the ILC on each testing day for Pmax and RPMpk (ICCs were calculated separately for upright and supine). gKendall fs Tau a h was used to determine the association between ILC Pmax and isokinetic and leg press data. RESULTS: For Pmax, significant differences were found between days 1 and 2 (upright: p = 0.018; supine: p = 0.014) and between days 1 and 3 (upright: p = 0.001; supine: p = 0.002), but not between days 2 and 3 (upright: p = 0.422; supine: p = 0.501). Pmax ICC values were greater than or equal to 0.97 for all days in both positions. Also, no significant differences between upright and supine postures were found for Pmax. No significant differences between days were found for RPMpk; however, there was a significant posture effect (upright greater than supine). Moderate correlations were observed between ILC Pmax and isokinetic and leg press tests (upright: 0.64-0.79, supine: 0.52-0.82). CONCLUSIONS: Overall, ILC is a very reliable test. Since a significant difference was found between day 1 and the other ILC testing days, it is suggested that day 1 of ILC testing should be used as a familiarization session to allow for subject learning. No significant difference in Pmax was seen from test 3 to test 6. However, an increase of 1.3% was observed from test 4 to test 6. Therefore, although 4 tests may be sufficient for most subjects to produce Pmax, in some cases 6 tests may be required. PRACTICAL APPLICATIONS: No differences were seen in Pmax between upright and supine positions despite differing RPMpk. This suggests that ILC testing can be used to provide reliable testing both in an upright position (appropriate for athletes) and in research (e.g., bed rest) or rehabilitation settings where supine testing is necessary. Future research should evaluate whether peak power measurements obtained with the ILC are sensitive to changes such as that observed with training and de-training.

Wickwire, P. J.↗
NASA NTRS2016

Comparison of Structural and Functional Ocular Outcomes Between 14- and 70 Day Bed Rest

Purpose: To compare structural and functional ocular outcomes in healthy human subjects undergoing 14- and/or 70-day head-down-tilt bed rest (HDTBR). We hypothesized the amount of HDTBR-induced ocular changes be affected by the HDTBR duration. Methods: The studies were conducted at the NASA Flight Analogs Research Unit, The University of Texas Medical Branch at Galveston, Galveston, TX. Participants were selected using NASA standard screening procedures. Standardized NASA screening procedures and bed rest conditions (e.g., strict sleep-wake cycle, standardized diet, continuous video monitoring) were implemented in both studies. Participants maintained a 6deg HDTBR position for 14 and/or 70 consecutive days and did not engage in exercise. Weekly ophthalmological examinations were conducted in the sitting (pre/post-bed rest only) and HDT positions. Ocular outcomes of interest included: near best-corrected visual acuity (BCVA); spherical equivalent, as determined by cycloplegic autorefraction; Goldmann applanation tonometry and iCare (Icare Finland Oy, Espoo, Finland) intraocular pressure (IOP) measurement; color vision; red dot test; modified Amsler grid test; confrontational visual field; stereoscopic color fundus photography; Spectralis OCT (Heidelberg Engineering, GmbH, Heidelberg, Germany) retinal nerve fiber layer thickness (RNFLT), peripapillary and macular retinal thicknesses. Mixed-effects linear models were used to compare pre- and post-HDTBR observations between 14- and 70-day HDTBR for our continuously scaled outcomes.

Cromwell, R. L.↗
NASA NTRS2018

Quantitative Daily Maps of PM 2.5 Episodes for California and Other Regions: Satellite Column Water and Optical Depth as Allied Tracers of Dilution

The Western US and many regions globally present daunting difficulties in understanding PM 2.5 episodes. We evaluate extensions of a method independent of modeled source-description and transport/transformation and using several satellite remote sensing products from imaging spectrometers. The San Joaquin Valley (SJV) especially suffers few-day episodes due to shallow mixing; PM 2.5 retrieval suffers low satellite AOT (Aerosol Optical Thickness) and bright surfaces.Nevertheless, we find residual errors in our maps of of typically 5-8 micrograms per cubic meter. Episodes in the Valley reaching 60-100 micrograms per cubic meter. These maps detail pollution from Interstate 5 at the scale of a few kilometers. The maps are based on NASA's MODerate resolution Imaging Spectrometer (MODIS) data at circa 1 kilometer as processed with the Multi-Angle Implementation of Atmospheric Correction. The Bay Area Air Quality Management District has requested that we test our methods in their challenging environment characterized by multiple sub-basins defined by complex topography. Our tests suggest that nearly similar precision may be expected for wintertime conditions with high PM 2.5 . We note difficulties when measured PM 2.5 is less than 8-10 micrograms per cubic meter, but good relative precision when PM 2.5 rises above 20; i.e. in episodes of concern for morbidity and mortality. Our method stresses physically meaningful functions of MODIS-MAIAC (Multi-Angle Implementation of Atmospheric Correction)-derived AOD (Aerosol Optical Depth) and total water vapor column. A mixed-effects statistical model exploiting existing station data works powerfully to allow us daily AOT-to-PM 2.5 relationships that allow a calibration of the map. In those cases where water vapor and particles have generally similar surface sources, using the ratio of AOT / Column_water can improve the daily calibrations so as to reach our quoted precision. We briefly present some cartoon idealizations that explain this success and also the likely reasons that our mixed effects model (or "daily calibration") works; also when it should not work. The combined satellite/mixed-effects model works best for wintertime San Joaquin Valley episodes, where the meteorology of particle and H2O(v) dilution is quite appropriate. We extended and tested the methodology (a) for the Bay Area wintertime situations and (b) for smoke plume events (e.g. the October 2017 fire events of the Sonoma area). Our SJV work was evaluated using NASA's DISCOVER-AQ (Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality) airborne measurements, and by season- long measurements in Fresno. If the composition and size distribution of the aerosols can be assessed for the regions we describe, retrievals should have improved accuracy.

Chatfield, Robert B.↗
NASA NTRS2018

Evaluation of the Danish Aerospace Corporation Portable Pulmonary Function System

A research project designed to investigate changes in maximal oxygen consumption (VO2max) during and following long duration flight on the International Space Station (ISS) has recently been completed. The device used to measure VO2 on board ISS, the Portable Pulmonary Function System (PPFS) manufactured by the Danish Aerospace Corporation (DAC), is based on previous-generation devices manufactured by DAC, but the PPFS has not been validated for analyzing metabolic gases or measuring cardiac output (Qc). The purpose of the present evaluation is to compare PPFS metabolic gas analysis measurements to measurements obtained using a clinically-validated system (ParvoMedics TrueOne(c) 2400 system; Parvo). In addition, Qc data collected with the PPFS were compared to Qc measurements from echocardiography. METHODS: Ten subjects completed three cycle exercise tests to maximal exertion. The first test was conducted to determine each subject's VO2max and set the work rates for the second and third (comparison) tests. The protocol for the two comparison tests consisted of three 5-minute stages designed to elicit 25%, 50%, and 75% VO2max (based upon results from the initial test), followed by 1-minute stages of increasing work rate (25 watt/minute) until the subject reached maximal effort. During one of the two comparison tests, metabolic gases and Qc were assessed with the PPFS; metabolic gases and Qc were assessed with the Parvo and by echocardiography, respectively, during the other test. The order of the comparison tests was counterbalanced. VO2max and maximal work rate during the comparison tests were compared using t tests. Mixed-effects regression modeling was used to analyze submaximal data. RESULTS: All of the data were within normal physiological ranges. The PPFS-measured values for VO2max were 6% lower than values obtained with the Parvo (PPFS: 3.11 +/- 0.75 L/min; Parvo: 3.32 +/- 0.87 L/min; mean +/- standard deviation; P = 0.02); this difference is probably due to flow restriction imposed by the PPFS Qc accessories. Submaximal VO2 values were slightly lower when measured with the PPFS, although differences were not physiologically relevant. The PPFS-measured values of submaximal carbon dioxide production (VCO2) were lower than the data obtained from Parvo, which could be attributed to lower fractions of expired carbon dioxide measured by the PPFS. The PPFS Qc values tended to be lower than echocardiography-derived values. CONCLUSIONS: The results of the present study indicate a need to further examine the PPFS and to better quantify its reproducibility; however, none of the findings of the current evaluation indicate that the PPFS needs to be replaced or modified.

Portable Pulmonary Function System (PPFS)↗
NASA NTRS2015

NASA's Functional Task Test: High Intensity Exercise Improves the Heart Rate Response to a Stand Test Following 70 Days of Bedrest

Cardiovascular adaptations due to spaceflight are modeled with 6deg head-down tilt bed rest (BR) and result in decreased orthostatic tolerance. We investigated if high-intensity resistive and aerobic exercise with and without testosterone supplementation would improve the heart rate (HR) response to a 3.5-min stand test and how quickly these changes recovered following BR. During 70 days of BR male subjects performed no exercise (Control, n=10), high intensity supine resistive and aerobic exercise (Exercise, n=9), or supine exercise plus supplemental testosterone (Exercise+T, n=8; 100 mg i.m., weekly in 2-week on/off cycles). We measured HR for 2 min while subjects were prone and for 3 min after standing twice before and 0, 1, 6, and 11 days after BR. Mixed-effects linear regression models were used to evaluate group, time, and interaction effects. Compared to pre-bed rest, prone HR was elevated on BR+0 and BR+1 in Control, but not Exercise or Exercise+T groups, and standing HR was greater in all 3 groups. The increase in prone and standing HR in Control subjects was greater than either Exercise or Exercise+T groups and all groups recovered by BR+6. The change in HR from prone to standing more than doubled on BR+0 in all groups, but was significantly less in the Exericse+T group compared to the Control, but not Exercise group. Exercise reduces, but does not prevent the increase in HR observed in response to standing. The significantly lower HR response in the Exercise+T group requires further investigation to determine physiologic significance.

Laurie, Steven S.↗
NASA NTRS2019

Satellite Mapping of PM2.5 Episodes in the Wintertime San Joaquin Valley: A "Static" Model Using Column Water Vapor

The use of satellite Aerosol Optical Thickness (AOT) from imaging spectrometers has been successful in quantifying and mapping high PM2.5 (particulate matter mass <2.5µm diameter) episodes for pollution abatement and health studies. However, some regions have high PM2.5 but poor estimation success. The challenges in using Aerosol Optical Thickness (AOT) from imaging spectrometers to characterize PM2.5 worldwide was especially evident in the wintertime San Joaquin Valley (SJV). The SJV's attendant difficulties of high-albedo surfaces and very shallow, variable vertical mixing also occur in other significantly polluted regions around the world. We report on more accurate PM2.5 maps for the whole-winter period in the SJV, Nov 14, 2012?Dec 11, 2013. Intensive measurements by including NASA aircraft were made for several weeks in that winter, the DISCOVER-AQ California mission.We found success with a relatively simple method based on calibration and checking with surface monitors and a characterization of vertical mixing, and incorporating specific understandings of the region's climatology. We estimate PM2.5 to within ~7µgm?3 RMSE and with R values of ~0.9, based on remotely sensed MAIAC (Multi-Angle Implementation of Atmospheric Correction) observations, and that certain further work will improve that accuracy. Mapping is at 1km resolution. This allows a time sequence of mapped aerosols at 1km for cloud-free days. We describe our technique as a "static estimation". Estimation procedures like this one, not dependent on well-mapped source strengths or on transport error, should help full source-driven simulations by deconstructing processes. They also provide a rapid method to create a long-term climatology.Essential features of the technique are (a) daily calibration of the AOT to PM2.5 using available surface monitors, and (b) characterization of mixed-layer dilution using column water vapor (CWV, otherwise "precipitable water"). We noted that on multi-day timescales both water vapor and particles share near-surface sources and both fall to very low values with altitude; indeed, both are largely removed by precipitation. The existence of layers of H2O or aerosol not within the mixed layer adds complexity, but mixed-effects statistical regression captures essential proportionality of PM2.5 and the ratio variable (AOT/CWV). Accuracy is much higher than previous statistical models, and can be extended to the whole Aqua-satellite data record. The maps and time-series we show suggest a repeated pattern for large valleys like the SJV ? progressive stabilization of the mixing height after frontal passages: PM2.5 is somewhat more determined by day-by-day changes in mixing than it is by the progressive accumulation of pollutants (revealed as increasing AOT).

Chatfield, Robert B.↗
NASA NTRS2015

Techniques for Estimating Emissions Factors from Forest Burning: ARCTAS and SEAC4RS Airborne Measurements Indicate Which Fires Produce Ozone

Previous studies of emission factors from biomass burning are prone to large errors since they ignore the interplay of mixing and varying pre-fire background CO2 levels. Such complications severely affected our studies of 446 forest fire plume samples measured in the Western US by the science teams of NASA's SEAC4RS and ARCTAS airborne missions. Consequently we propose a Mixed Effects Regression Emission Technique (MERET) to check techniques like the Normalized Emission Ratio Method (NERM), where use of sequential observations cannot disentangle emissions and mixing. We also evaluate a simpler "consensus" technique. All techniques relate emissions to fuel burned using C(sub burn) = delta C(sub tot) added to the fire plume, where C(sub tot) approximately equals (CO2 + CO). Mixed-effects regression can estimate pre-fire background values of Ctot (indexed by observation j) simultaneously with emissions factors indexed by individual species i, delta epsilon lambda tau alpha−x(sub i)/(C(sub burn))i,j., MERET and "consensus" require more than two emissions indicators. Our studies excluded samples where exogenous CO or CH4 might have been fed into a fire plume, mimicking emission. We sought to let the data on 13 gases and particulate properties suggest clusters of variables and plume types, using non-negative matrix factorization (NMF). While samples were mixtures, the NMF unmixing suggested purer burn types. Particulate properties (bscat, babs, SSA, AAE) and gas-phase emissions were interrelated. Finally, we sought a simple categorization useful for modeling ozone production in plumes. Two kinds of fires produced high ozone: those with large fuel nitrogen as evidenced by remnant CH3CN in the plumes, and also those from very intense large burns. Fire types with optimal ratios of delta-NOy/delta- HCHO associate with the highest additional ozone per unit Cburn, Perhaps these plumes exhibit limited NOx binding to reactive organics. Perhaps these plumes exhibit limited NOx binding to reactive organics.

fires↗
NASA NTRS2016

Techniques for Estimating Emissions Factors from Forest Burning: ARCTAS and SEAC4RS Airborne Measurements Indicate which Fires Produce Ozone

Previous studies of emission factors from biomass burning are prone to large errors since they ignore the interplay of mixing and varying pre-fire background CO2 levels. Such complications severely affected our studies of 446 forest fire plume samples measured in the Western US by the science teams of NASA's SEAC4RS and ARCTAS airborne missions. Consequently we propose a Mixed Effects Regression Emission Technique (MERET) to check techniques like the Normalized Emission Ratio Method (NERM), where use of sequential observations cannot disentangle emissions and mixing. We also evaluate a simpler "consensus" technique. All techniques relate emissions to fuel burned using C(burn) = delta C(tot) added to the fire plume, where C(tot) approximately equals (CO2 = CO). Mixed-effects regression can estimate pre-fire background values of C(tot) (indexed by observation j) simultaneously with emissions factors indexed by individual species i, delta, epsilon lambda tau alpha-x(sub I)/C(sub burn))I,j. MERET and "consensus" require more than emissions indicators. Our studies excluded samples where exogenous CO or CH4 might have been fed into a fire plume, mimicking emission. We sought to let the data on 13 gases and particulate properties suggest clusters of variables and plume types, using non-negative matrix factorization (NMF). While samples were mixtures, the NMF unmixing suggested purer burn types. Particulate properties (b scant, b abs, SSA, AAE) and gas-phase emissions were interrelated. Finally, we sought a simple categorization useful for modeling ozone production in plumes. Two kinds of fires produced high ozone: those with large fuel nitrogen as evidenced by remnant CH3CN in the plumes, and also those from very intense large burns. Fire types with optimal ratios of delta-NOy/delta- HCHO associate with the highest additional ozone per unit Cburn, Perhaps these plumes exhibit limited NOx binding to reactive organics. Perhaps these plumes exhibit limited NOx binding to reactive organics

ozone↗
NASA NTRS2014

Revised (Mixed-Effects) Estimation for Forest Burning Emissions of Gases and Smoke, Fire/Emission Factor Typology, and Potential Remote Sensing Classification of Types for Ozone and Black-Carbon Simulation

We summarize recent progress (a) in correcting biomass burning emissions factors deduced from airborne sampling of forest fire plumes, (b) in understanding the variability in reactivity of the fresh plumes sampled in ARCTAS (2008), DC3 (2012), and SEAC4RS (2013) airborne missions, and (c) in a consequent search for remotely sensed quantities that help classify forest-fire plumes. Particle properties, chemical speciation, and smoke radiative properties are related and mutually informative, as pictures below suggest (slopes of lines of same color are similar). (a) Mixed-effects (random-effects) statistical modeling provides estimates of both emission factors and a reasonable description of carbon-burned simultaneously. Different fire plumes will have very different contributions to volatile organic carbon reactivity; this may help explain differences of free NOx(both gas- and particle-phase), and also of ozone production, that have been noted for forest-fire plumes in California. Our evaluations check or correct emission factors based on sequential measurements (e.g., the Normalized Ratio Enhancement and similar methods). We stress the dangers of methods relying on emission-ratios to CO. (b) This work confirms and extends many reports of great situational variability in emissions factors. VOCs vary in OH reactivity and NOx-binding. Reasons for variability are not only fuel composition, fuel condition, etc., but are confused somewhat by rapid transformation and mixing of emissions. We use "unmixing" (distinct from mixed-effects) statistics and compare briefly to approaches like neural nets. We focus on one particularly intense fire the notorious Yosemite Rim Fire of 2013. In some samples, NOx activity was not so suppressed by binding into nitrates as in other fires. While our fire-typing is evolving and subject to debate, the carbon-burned delta(CO2+CO) estimates that arise from mixed effects models, free of confusion by background-CO2 variation, should provide a solid base for discussion. (c) We report progress using promising links we find between emissions-related "fire types" and promising features deducible from remote observations of plumes, e.g., single scatter albedo, Angstrom exponent of scattering, Angstrom exponent of absorption, (CO column density)/(aerosol optical depth).

Remote sensing↗
DOE OSTI2026

Data driven investigation to understand the influence of total solids on biological biogas upgrading

In situ biogas upgrading achieves CO 2 conversion to CH 4 via hydrogenotrophic methanogenesis; however, gas-liquid mass transfer constraints limit the upgrading performance. Recognizing that optimization studies often underrepresent the effects of total solids (TS) and organic loading rate (OLR), this study undertook a holistic, statistics driven assessment of operating conditions for in situ H 2 assisted biogas upgrading, centering the analysis on TS and OLR. A dataset of 31 studies was compiled and comprised 99 observations. A rigorous analytical framework was employed, combining data standardization, fixed- and random-effects (REML) weighted regressions with cluster-robust errors, stratified analyses, and machine learning. Mixed-effects meta regression indicated that TS was the main factor explaining differences of methane fraction (CH 4 %) when considering the between studies heterogeneity. Focusing on a near-stoichiometric subset (H 2 /CO 2 ≈ 4:1), TS remained significant. Stratified results showed a stronger negative relationship between TS and CH 4 % in UASB reactors than in CSTRs, with a negative effect under mesophilic conditions and no significant effect under thermophilic conditions. A Random Forest model corroborated the statistical findings, consistently ranking H 2 /CO 2 ratio, OLR, TS, and hydrogen injection rate (HIR) as the most influential predictors. These findings delineate trends across increasing TS levels, particularly between 1% and 10%, and provide preliminary insights for TS above 15% in in situ biogas upgrading. They further provide insights for the influence of TS by reactor type and temperature, thereby advancing the evidence base for implementing biological CO 2 conversion to CH 4 in practice.

In situ biogas upgrading↗