Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “corrections”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Extending 3D Near-Cloud Corrections from Shorter to Longer Wavelengths

Satellite observations have shown a positive correlation between cloud amount and aerosol optical thickness (AOT) that can be explained by the humidification of aerosols near clouds, and/or by cloud contamination by sub-pixel size clouds and the cloud adjacency effect. The last effect may substantially increase reflected radiation in cloud-free columns, leading to overestimates in the retrieved AOT. For clear-sky areas near boundary layer clouds the main contribution to the enhancement of clear sky reflectance at shorter wavelengths comes from the radiation scattered into clear areas by clouds and then scattered to the sensor by air molecules. Because of the wavelength dependence of air molecule scattering, this process leads to a larger reflectance increase at shorter wavelengths, and can be corrected using a simple two-layer model. However, correcting only for molecular scattering skews spectral properties of the retrieved AOT. Kassianov and Ovtchinnikov proposed a technique that uses spectral reflectance ratios to retrieve AOT in the vicinity of clouds; they assumed that the cloud adjacency effect influences the spectral ratio between reflectances at two wavelengths less than it influences the reflectances themselves. This paper combines the two approaches: It assumes that the 3D correction for the shortest wavelength is known with some uncertainties, and then it estimates the 3D correction for longer wavelengths using a modified ratio method. The new approach is tested with 3D radiances simulated for 26 cumulus fields from Large-Eddy Simulations, supplemented with 40 aerosol profiles. The results showed that (i) for a variety of cumulus cloud scenes and aerosol profiles over ocean the 3D correction due to cloud adjacency effect can be extended from shorter to longer wavelengths and (ii) the 3D corrections for longer wavelengths are not very sensitive to unbiased random uncertainties in the 3D corrections at shorter wavelengths.

aerosol optical depth↗
NASA NTRSDate not supplied

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets such as sex or age of the model organism used. In the present study, NASA GeneLab-hosted RNAseq datasets from rodent liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC, to determine statistical differences between datasets before and after correction, Principal Component Analysis, to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the standard approach. Thus, the most robust standard correction will be implemented in the GeneLab Visualization 2.0 platform when datasets are combined.

GeneLab, RNA-seq, Batch Correction↗
NASA NTRSDate not supplied

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the respective standard approach. Of the methods tested, standard ComBat and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

GeneLab↗
DOE OSTI2021

Drell-Yan q T resummation of fiducial power corrections at N 3 LL

We consider Drell-Yan production pp → V*X → LX at small q T << Q, where q T and Q are the total transverse momentum and invariant mass of the leptonic final state L. Experimental measurements require fiducial cuts on L, which in general introduce enhanced, linear power corrections in q T /Q. We show that they can be unambiguously predicted from factorization, and resummed to the same order as the leading-power contribution. For the fiducial q T spectrum, they constitute the complete linear power corrections. We thus obtain predictions for the fiducial q T spectrum to N 3 LL and next-to-leading-power in q T /Q. Matching to full NNLO ($α$$^{2}_{s}$), we find that the linear power corrections are indeed the dominant ones, and once included by factorization, the remaining fixed-order corrections become almost negligible below q T ≲ 40 GeV. We also discuss the implications for more complicated observables, and provide predictions for the fiducial Φ* spectrum at N 3 LL+NNLO. We find excellent agreement with ATLAS and CMS measurements of q T and Φ*. We also consider the $p$$^{ℓ}_{T}$ spectrum. We show that it develops leptonic power corrections in q T /(Q – 2$p$$^{ℓ}_{T}$), which diverge near the Jacobian peak $p$$^{ℓ}_{T}$ ~ Q/2 and must be kept to all powers to obtain a meaningful result there. Doing so, we obtain for the first time an analytically resummed result for the $p$$^{ℓ}_{T}$ spectrum around the Jacobian peak at N 3 LL+NNLO. Our method is based on performing a complete tensor decomposition for hadronic and leptonic tensors. We show that in practice this is equivalent to often-used recoil prescriptions, for which our results now provide rigorous, formal justification. Our tensor decomposition yields nine Lorentz-scalar hadronic structure functions, which for Z/γ* → ℓℓ or W → ℓν directly map onto the commonly used angular coefficients, but also holds for arbitrary leptonic final states. In particular, for suitably defined Born-projected leptons it still yields a LO-like angular decomposition even when including QED final-state radiation. Finally, we also discuss the application to q T subtractions. Including the unambiguously predicted fiducial power corrections significantly improves their performance, and in particular makes them applicable near kinematic edges where they otherwise break down due to large leptonic power corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Complete lowest order radiative corrections in semi-inclusive scattering of polarized particles

For this work, the lowest order radiative corrections to the cross section and asymmetries measured in experiments on semi-inclusive deep inelastic scattering of polarized particles were calculated. Both exact and leading log expressions were presented and discussed for the total correction that include the contributions from the processes of i) real photon emission with semi-inclusive processes, ii) loop diagrams, and iii) real photon emission with exclusive processes. Radiative corrections to the Sivers and Collins asymmetries in $\pi^+$ electroproduction were studied numerically within the kinematical conditions of modern experimental environments at Jefferson Laboratory (JLab). The Wandzura-Wilczek approximation for the semi-inclusive structure functions and MAID2007 parameterization for the six amplitudes of exclusive processes were used in numeric analyses. The results show that i) radiative effects can generate a correction comparable to the size of Sivers and Collins asymmetries at the Born level, ii) there is good agreement between the exact and leading-order corrections, iii) external functions (that is, other than the Sivers and Collins functions in the respective asymmetries) can generate a contribution to the radiative correction up to 20%, and iv) there exists a strong dependence of the radiative correction on the models for semi-inclusive and exclusive structure functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Afterpulse correction for micro-pulse lidar to improve middle and upper tropospheric aerosol measurements

Micro-pulse lidar (MPL) measurements have been widely used in atmospheric research over the past few decades. However, the MPL afterpulse noise has a large impact on the MPL aerosol measurement in the middle and upper troposphere, and an effective correction method is still lacking. Here, a new afterpulse correction approach is presented by using measurements with low-level optically thick clouds to act as the lids blocking atmospheric signals beyond the clouds completely. Examples are provided to illustrate the effectiveness of this correction method. Using one-year 2014 MPL measurements at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) North Slope of Alaska (NSA) site, the impact of the correction on the aerosol measurements is quantified. The relative error (RE, %) of attenuated backscattering ratio (ABR) without the afterpulse correction is 30% and 190% at the height of 1.00 km and 9.00 km (AGL), and the RE is larger for weaker aerosol signals (ABR<2). The RE of linear depolarization ratio (LDR), which is significantly higher than that of ABR for a given aerosol layer, is highly non-linear above 3.00 km, with a value of 300%-900% for weak aerosols (ABR<2) above 3.00 km. Therefore, the afterpulse correction is critical for the middle and upper tropospheric aerosol observations. We demonstrated that our afterpulse correction can provide reasonable particle depolarization ratio (PDR) to properly identify dust aerosols. The newly developed method can be applied to long-term MPL measurements to support aerosol and mixed-phase/ice cloud interaction study at the NSA site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Parallax-corrected VISST-derived pixel-level products from satellite GOES-16

The NASA Langley group led by William Smith produced GOES-16 satellite cloud retrievals over an approximate 10 by 10 degree region over the CACTI field campaign location. These retrievals are described here: https://www.arm.gov/capabilities/vaps/visst and are available for download here . They use algorithms historically called VISST that are now referred to as SatCORPS. More information can be found in Trepte et al. (2019), Minnis et al. (2021), and Yost et al. (2021). If using this dataset, please cite these references, the CACTI VISST dataset DOI found at the download link above, and this dataset’s DOI. The CACTI VISST pixel-level retrievals are on a 2 km spatial grid and available every 15 minutes (every 10 minutes late in the campaign), producing 21,765 files for the entire field campaign between October 2018 and April 2019. They are not corrected for parallax error, which is an offset in the actual geographical location of a cloud above the surface due to the satellite viewing the cloud partly from the side off nadir. This dataset applies a correction for parallax using the location relative to the satellite and the retrieved cloud top height above the surface, which allows the dataset to be geo-located with surface-based observations. The parallax correction for each location depends on the longitude, latitude and cloud top height above ground level (AGL) for that longitude and latitude in the original VISST files. The cloud top height AGL requires first computing the surface elevation at each VISST grid point. Data from the Advanced Spaceborne Thermal Emission and Reflection (ASTER) Global Digital Elevation Map Version 3 at 30-m resolution is projected onto the VISST grid using conservative coarsening (conserving surface elevation) in the xESMF Python package. The surface elevation is then subtracted from the VISST-retrieved cloud top height above mean sea level. These cloud top heights AGL are then combined with longitude and latitude to estimate the latitude and longitude corrections. Due to variability in cloud top height, the parallax shifts produce an irregular grid of values since higher cloud tops are shifted further than lower cloud tops. A ball tree-based neighbor search with Haversine distance is performed using the Python-based scikit-learn library to find the nearest VISST grid point to each parallax correction-shifted point. The data value of the shifted point is then assigned to that VISST grid point. In this manner, the irregular geographical shifts to correct for parallax are projected back to the rectilinear VISST grid. Because relatively higher clouds should obscure lower clouds, the variable values for the highest cloud top are preferentially chosen if two or more values are assigned to a grid point. The parallax correction should be viewed as an improved but still imperfect estimation of the cloud top locations, largely because the cloud top height is an imperfect retrieval. Please see the attached README document for further information. Users are encouraged to contact the authors with any additional questions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Absolute oscillator strengths for the iron group - A correction to the data of Warner for the second spectra and a comment on some of the semiempirical results of Kurucz and Peytremann

It has already been shown that the Fe II absolute oscillator strengths (f-values) measured by Warner (1967) contain a systematic error, the magnitude of which depends on the energy of the transition's upper level. The present paper demonstrates that this error affects all of Warner's oscillator strengths for the second spectra of the iron group and develops a correction formula for Warner's results. A comparison of the corrected f-values with other measured f-values for Ti II, V II, and Fe II indicates that the correction improves the usefulness of Warner's data and that the corrected f-values have no large systematic dependence on wavelength, level energy, or line strength. Solar abundances for iron-group elements computed using the original and corrected f-values are compared with recently suggested abundances, and it is found that the corrected f-values yield more consistent results. Semiempirical oscillator strengths calculated by Kurucz and Peytremann (1975) for Fe I, Ti II, and V II are also compared with measured data to evaluate the accuracy of the calculations. It is found that the absolute scales of these results are roughly correct, but some of them are subject to very large errors.

Smith, P. L.↗
NASA NTRS1982

LANDSAT-D accelerated payload correction subsystem output computer compatible tape format

The NASA GSFC LANDSAT-D Ground Segment (GS) is developing an Accelerated Payload Correction Subsystem (APCS) to provide Thematic Mapper (TM) image correction data to be used outside the GS. This correction data is computed from a subset of the TM Payload Correction Data (PCD), which is downlinked from the spacecraft in a 32 Kbps data stream, and mirror scan correction data (MSCD), which is extracted from the wideband video data. This correction data is generated in the GS Thematic Mapper Mission Management Facility (MMF-T), and is recorded on a 9-track 1600 bit per inch computer compatible tape (CCT). This CCT is known as a APCS Output CCT (AOT). The AOT follows standardized corrections with respect to data formats, record construction and record identification. Applicable documents are delineated; common conventions which are used in further defining the structure, format and content of the AOT are defined; and the structure and content of the AOT are described.

Source record↗
NASA NTRS1987

A model-free method for mass spectrometer response correction

A new method for correction of mass spectrometer output signals is described. Response-time distortion is reduced independently of any model of mass spectrometer behavior. The delay of the system is found first from the cross-correlation function of a step change and its response. A two-sided time-domain digital correction filter (deconvolution filter) is generated next from the same step response data using a regression procedure. Other data are corrected using the filter and delay. The mean squared error between a step response and a step is reduced considerably more after the use of a deconvolution filter than after the application of a second-order model correction. O2 consumption and CO2 production values calculated from data corrupted by a simulated dynamic process return to near the uncorrupted values after correction. Although a clean step response or the ensemble average of several responses contaminated with noise is needed for the generation of the filter, random noise of magnitude not above 0.5 percent added to the response to be corrected does not impair the correction severely.

Shykoff, Barbara E.↗
NASA NTRS1992

Present status and future prospects for ionospheric propagation corrections for precise time transfer using GPS

The ionosphere can be the greatest variable source of error in precise time transfer using Global Positioning System (GPS) satellites. For single frequency GPS users, the ionospheric correction algorithm can provide an approximate 50 percent r.m.s. correction to the time delay, but users who desire a more complete correction must make actual measurements of ionospheric time delay along the path to the GPS satellite. Fortunately, at least three commercial GPS receivers, specifically designed to measure and correct for ionospheric time delay, are now, or soon will be, available. Initial operation with two different types of GPS ionospheric receivers demonstrated a high degree of accuracy in measuring the ionospheric group delay. Results of these measurements are presented. For those who use a model to correct for ionospheric time delay, it is tempting to use daily values of solar 10.7 cm radio flux to correct a monthly average ionospheric time delay model for each day's operation. The results of correlation of daily maximum ionospheric time delay against solar radio flux values show a poor correlation will be obtained by this procedure. Prospects for improving ionospheric corrections during the declining phase of the present solar cycle are discussed.

Klobuchar, John A.↗
NASA NTRS1997

Correction of Rayleigh Scattering Effects in Cloud Optical Thickness Retrievals

We present results that demonstrate the effects of Rayleigh scattering on the 9 retrieval of cloud optical thickness at a visible wavelength (0.66 Am). The sensor-measured radiance at a visible wavelength (0.66 Am) is usually used to infer remotely the cloud optical thickness from aircraft or satellite instruments. For example, we find that without removing Rayleigh scattering effects, errors in the retrieved cloud optical thickness for a thin water cloud layer (T = 2.0) range from 15 to 60%, depending on solar zenith angle and viewing geometry. For an optically thick cloud (T = 10), on the other hand, errors can range from 10 to 60% for large solar zenith angles (0-60 deg) because of enhanced Rayleigh scattering. It is therefore particularly important to correct for Rayleigh scattering contributions to the reflected signal from a cloud layer both (1) for the case of thin clouds and (2) for large solar zenith angles and all clouds. On the basis of the single scattering approximation, we propose an iterative method for effectively removing Rayleigh scattering contributions from the measured radiance signal in cloud optical thickness retrievals. The proposed correction algorithm works very well and can easily be incorporated into any cloud retrieval algorithm. The Rayleigh correction method is applicable to cloud at any pressure, providing that the cloud top pressure is known to within +/- 100 bPa. With the Rayleigh correction the errors in retrieved cloud optical thickness are usually reduced to within 3%. In cases of both thin cloud layers and thick ,clouds with large solar zenith angles, the errors are usually reduced by a factor of about 2 to over 10. The Rayleigh correction algorithm has been tested with simulations for realistic cloud optical and microphysical properties with different solar and viewing geometries. We apply the Rayleigh correction algorithm to the cloud optical thickness retrievals from experimental data obtained during the Atlantic Stratocumulus Transition Experiment (ASTEX) conducted near the Azores in June 1992 and compare these results to corresponding retrievals obtained using 0.88 Am. These results provide an example of the Rayleigh scattering effects on thin clouds and further test the Rayleigh correction scheme. Using a nonabsorbing near-infrared wavelength lambda (0.88 Am) in retrieving cloud optical thickness is only applicable over oceans, however, since most land surfaces are highly reflective at 0.88 Am. Hence successful global retrievals of cloud optical thickness should remove Rayleigh scattering effects when using reflectance measurements at 0.66 Am.

Wang, Meng-Hua↗
NASA NTRS2003

Assessment, Validation, and Refinement of the Atmospheric Correction Algorithm for the Ocean Color Sensors

The primary focus of this proposed research is for the atmospheric correction algorithm evaluation and development and satellite sensor calibration and characterization. It is well known that the atmospheric correction, which removes more than 90% of sensor-measured signals contributed from atmosphere in the visible, is the key procedure in the ocean color remote sensing (Gordon and Wang, 1994). The accuracy and effectiveness of the atmospheric correction directly affect the remotely retrieved ocean bio-optical products. On the other hand, for ocean color remote sensing, in order to obtain the required accuracy in the derived water-leaving signals from satellite measurements, an on-orbit vicarious calibration of the whole system, i.e., sensor and algorithms, is necessary. In addition, it is important to address issues of (i) cross-calibration of two or more sensors and (ii) in-orbit vicarious calibration of the sensor-atmosphere system. The goal of these researches is to develop methods for meaningful comparison and possible merging of data products from multiple ocean color missions. In the past year, much efforts have been on (a) understanding and correcting the artifacts appeared in the SeaWiFS-derived ocean and atmospheric produces; (b) developing an efficient method in generating the SeaWiFS aerosol lookup tables, (c) evaluating the effects of calibration error in the near-infrared (NIR) band to the atmospheric correction of the ocean color remote sensors, (d) comparing the aerosol correction algorithm using the singlescattering epsilon (the current SeaWiFS algorithm) vs. the multiple-scattering epsilon method, and (e) continuing on activities for the International Ocean-Color Coordinating Group (IOCCG) atmospheric correction working group. In this report, I will briefly present and discuss these and some other research activities.

Wang, Menghua↗
NASA NTRS2016

Atmospheric Correction for Satellite Ocean Color Radiometry

This tutorial is an introduction to atmospheric correction in general and also documentation of the atmospheric correction algorithms currently implemented by the NASA Ocean Biology Processing Group (OBPG) for processing ocean color data from satellite-borne sensors such as MODIS and VIIRS. The intended audience is graduate students or others who are encountering this topic for the first time. The tutorial is in two parts. Part I discusses the generic atmospheric correction problem. The magnitude and nature of the problem are first illustrated with numerical results generated by a coupled ocean-atmosphere radiative transfer model. That code allow the various contributions (Rayleigh and aerosol path radiance, surface reflectance, water-leaving radiance, etc.) to the topof- the-atmosphere (TOA) radiance to be separated out. Particular attention is then paid to the definition, calculation, and interpretation of the so-called "exact normalized water-leaving radiance" and its equivalent reflectance. Part I ends with chapters on the calculation of direct and diffuse atmospheric transmittances, and on how vicarious calibration is performed. Part II then describes one by one the particular algorithms currently used by the OBPG to effect the various steps of the atmospheric correction process, viz. the corrections for absorption and scattering by gases and aerosols, Sun and sky reflectance by the sea surface and whitecaps, and finally corrections for sensor out-of-band response and polarization effects. One goal of the tutorial-guided by teaching needs- is to distill the results of dozens of papers published over several decades of research in atmospheric correction for ocean color remote sensing.

MODIS↗
NASA NTRS2017

GEO-LEO Reflective Band Inter-Comparison with BRDF and Atmospheric Scattering Corrections

The inter-comparison of the reflective solar bands (RSB) between the instruments onboard a geostationary orbit satellite and a low Earth orbit satellite is very helpful in assessing their calibration consistency. Himawari-8 was launched 7 October 2014 and GOES-R was launched on 19 November 2016. Unlike previous GOES instruments, the Advanced Himawari Imager (AHI) on Himawari-8 and the Advanced Baseline Imager (ABI) on GOES-R have onboard calibrators for the RSB. Independent assessment of calibration is nonetheless important to enhance their product quality. MODIS (Moderate Resolution Imaging Spectroradiometer) and VIIRS (Visible Infrared Imaging Radiometer Suite) can provide good references for sensor calibration. In this work, the inter-comparison between AHI and VIIRS is performed over a pseudo-invariant target. The use of stable and uniform calibration sites provides comparison with accurate adjustment for band spectral difference, reduction of impact from pixel mismatching, and consistency of BRDF (Bidirectional Reflectance Distribution Function) and atmospheric correction. The site used is the Strzelecki Desert in Australia. Due to the difference in solar and view angles, two corrections must be applied in order to compare the measurements. The first is the atmospheric scattering correction applied to the top of atmosphere reflectance measurements. The second correction is applied to correct the BRDF effect. The atmospheric correction is performed using a vector version of the Second Simulation of a Satellite Signal in the Solar Spectrum (6SV) model and the BRDF correction is performed using a semi-empirical model. Our results show that AHI band 1 (0.47 microns) has a good agreement with VIIRS band M3 within 0.15 percent. AHI band 5 (1.61 microns) shows the largest difference (5.09 percent) with VIIRS band M10, while AHI band 5 shows the least difference (1.87 percent) in comparison with VIIRS band I3. The methods developed in this work can also be directly applied to assess GOES-16/ABI (Geostationary Operational Environment Satellite16 / Advanced Baseline Imager) calibration consistency, a topic we will address in the future.

The inter-comparison of the reflective solar bands↗
NASA NTRSDate not supplied

Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods

Background: Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. Methods: In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, the median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

Finsam Samson↗
NASA NTRSDate not supplied

Towards Energy Scale Calibration and Drift Correction of TES Detectors for Athena X-IFU

The Athena X-Ray Integral Field Unit (X-IFU) comprises a 2376-pixel array of transition edge sensors (TES) read out with time-division multiplexing (TDM). X-IFU will provide spatially resolved, high-resolution spectroscopy (2.5 eV full-width-half-maximum up to 7 keV) over the energy range 0.2 to 12 keV, with an absolute energy scale accuracy of 0.4 eV. The energy scale function maps the optimally filtered pulse height, in arbitrary engineering units, to real calibrated energy. Uncertainties in the calibration can result from imperfect fitting of the energy scale between the known calibration points. Furthermore, temporal changes in the TES operating environment, such as heat-sink temperature, magnetic field and bias voltage, can cause significant variations in the detector gain function over time. If not properly corrected, this can result in degradation of the energy resolution, and systematic errors in the absolute energy scale. The non-linear nature of TES detectors, coupled with the possibility of multiple simultaneously occurring sources of drift, can make effective corrections over the full bandpass of the instrument extremely challenging. Athena X-IFU will employ an on-board calibration source that provides known reference x-ray lines. This provides real-time monitoring of the gain stability of the detector system and information that can be used to correct for gain drifts. For X-IFU the baseline approach is to measure a series of calibration curves under different environmental conditions, which bound the expected drifts the instrument is predicted to see over the course of the mission. Using the information from the in-flight calibration source, these energy scale functions can be interpolated to generate a new corrected energy scale as a function of time. In this paper we discuss progress towards demonstrating that the X-IFU energy scale requirements can be met. We present measurements on ~ 200 pixels in a prototype X-IFU array read out with 8-column x 32-row TDM. We use a rotating target source containing 12 fluorescent targets to generate x-ray lines covering the energy range 4 keV (Sc-Kα) to 12 keV (Br-Kα). We present measurements of the non-linear energy scale function and show how variations in heat-sink temperature, TES bias voltage and magnetic field affect the shape of TES energy scale differently and introduce different residual gain errors over the bandpass. We explore different drift correction algorithms that use either a single or multiple referential lines to track and correct the gain from these various sources of drift. In addition to the pulse-height, the DC ‘baseline’ level of the TES can contain information about its bias conditions. Thus, we test a multi-parameter gain correction algorithm that attempts to incorporate both the pulse height and the additional baseline information into the algorithm.

Stephen J Smith↗
NASA NTRS2022

MODIS TEB Electronic Crosstalk Correction Update and Impact on L1B Product Uncertainty

The MODIS instruments onboard the Terra and Aqua satellites have been in operation for over22 and 20 years, respectively. The instruments’ calibration accuracy has been maintained, even with instrument degradation. Electronic crosstalk in the thermal emissive bands (TEB) is a known issue with an increasing impact on the calibration and product. The Terra MODIS photovoltaic (PV) longwave infrared (LWIR) bands crosstalk corrections have been applied in Collection 6.1 (C6.1). However, the electronic crosstalk contamination for some detectors in the mid-wave infrared (MWIR) bands and the Aqua PVLWIR bands affect the Level-1B (L1B) product’s measurement accuracy and image quality. In Collection 7 (C7), crosstalk corrections for select detectors in the Terra and Aqua MWIR and Aqua PV LWIR bands are applied. The entire mission crosstalk coefficients for the select detectors and bands are derived from scheduled lunar observations and populated in the form of look-up tables (LUTs). The Aqua PV-LWIR bands exhibit similar downward crosstalk trends as the Terra PV-LWIR bands, especially in recent years. The crosstalk coefficients and their trends provide a guideline for the correction application. Earth measurement analyses before and after the correction provide contamination and correction assessments. It has been shown that the product quality is enhanced with the crosstalk correction applied in C7. For C7, the crosstalk coefficient uncertainty is derived from the fit residuals between the measured values and a linear fit over a three-year sliding window. The uncertainty propagation is modeled and applied in the total uncertainty calculation intheL1B product. The TEB electronic crosstalk LUTs have been processed over the entire Terra and Aqua MODIS missions. This paper presents the C7 crosstalk correction, as well as its assessment and uncertainty propagation algorithm to the TEB uncertainty.

Tiejun Chang↗