Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Advancements in the Aerosol Robotic Network (AERONET) Version 3 database – automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements

The Aerosol Robotic Network (AERONET) has provided highly accurate, ground-truth measurements of the aerosol optical depth (AOD) using Cimel Electronique Sun–sky radiometers for more than 25 years. In Version 2 (V2) of the AERONET database, the near-real-time AOD was semiautomatically quality controlled utilizing mainly cloud-screening methodology, while additional AOD data contaminated by clouds or affected by instrument anomalies were removed manually before attaining quality-assured status (Level 2.0). The large growth in the number of AERONET sites over the past 25 years resulted in significant burden to the manual quality control of millions of measurements in a consistent manner. The AERONET Version 3 (V3) algorithm provides fully automatic cloud screening and instrument anomaly quality controls. All of these new algorithm updates apply to near-real-time data as well as post-field-deployment processed data, and AERONET reprocessed the database in 2018. A full algorithm redevelopment provided the opportunity to improve data inputs and corrections such as unique filter-specific temperature characterizations for all visible and near-infrared wavelengths, updated gaseous and water vapor absorption coefficients, and ancillary data sets. The Level 2.0 AOD quality-assured data set is now available within a month after post-field calibration, reducing the lag time from up to several months. Near-real-time estimated uncertainty is determined using data qualified as V3 Level 2.0 AOD and considering the difference between the AOD computed with the pre-field calibration and AOD computed with pre-field and post-field calibration. This assessment provides a near-real-time uncertainty estimate for which average differences of AOD suggest a +0.02 bias and one sigma uncertainty of 0.02, spectrally, but the bias and uncertainty can be significantly larger for specific instrument deployments. Long-term monthly averages analyzed for the entire V3 and V2 databases produced average differences (V3–V2) of +0.002 with a ±0.02 SD (standard deviation), yet monthly averages calculated using time-matched observations in both databases were analyzed to compute an average difference of −0.002 with a ±0.004 SD. The high statistical agreement in multiyear monthly averaged AOD validates the advanced automatic data quality control algorithms and suggests that migrating research to the V3 database will corroborate most V2 research conclusions and likely lead to more accurate results in some cases.

David M. Giles↗
NASA NTRS2023

Atmospheric Correction Inter-comparison eXercise, ACIX-II Land: An Assessment of Amospheric Correction Processors for Landsat 8 and Sentinel-2 Over Land

The correction of the atmospheric effects on optical satellite images is essential for quantitative and multi-temporal remote sensing applications. In order to study the performance of the state-of-the-art methods in an integrated way, a voluntary and open-access benchmark Atmospheric Correction Inter-comparison eXercise (ACIX) was initiated in 2016 in the frame of Committee on Earth Observation Satellites (CEOS) Working Group on Calibration & Validation (WGCV). The first exercise was extended in a second edition wherein twelve atmospheric correction (AC) processors, a substantially larger testing dataset and additional validation metrics were involved. The sites for the inter-comparison analysis were defined by investigating the full catalogue of the Aerosol Robotic Network (AERONET) sites for coincident measurements with satellites' overpass. Although there were more than one hundred sites for Copernicus Sentinel-2 and Landsat 8 acquisitions, the analysis presented in this paper concerns only the common matchups amongst all processors, reducing the number to 79 and 62 sites respectively. Aerosol Optical Depth (AOD) and Water Vapour (WV) retrievals were consequently validated based on the available AERONET observations. The processors mostly succeeded in retrieving AOD for relatively light to medium aerosol loading (AOD < 0.2) with uncertainties <0.08, while the overall uncertainty values were typically 0.23 ± 0.15. Better performances were observed for WV retrievals with >90% of the results falling within the suggested empirical specifications and with the Root Mean Square Error (RMSE) being mostly <0.25 g/cm2. Regarding Surface Reflectance (SR) validation two main approaches were followed. For the first one, a simulated SR reference dataset was computed over all of the test sites by using the 6SV (Second Simulation of the Satellite Signal in the Solar Spectrum vector code) full radiative transfer modelling (RTM) and AERONET measurements for the required aerosol variables and water vapour content. The performance assessment demonstrated that the retrievals were not biased for most of the bands. The uncertainties ranged from approximately 0.003 to 0.01 (excluding B01) for the best performing processors in both sensors' analyses. For the second one, measurements from the radiometric calibration network RadCalNet over La Crau (France) and Gobabeb (Namibia) were involved in the validation. The performance of the processors was in general consistent across all bands for both sensors and with low standard deviations (<0.04) between on-site and estimated surface reflectance. Overall, our study provides a good insight of AC algorithms' performance to developers and users, pointing out similarities and differences for AOD, WV and SR retrievals. Such validation though still lacks of ground-based measurements of known uncertainty to better assess and characterize the uncertainties in SR retrievals.

Atmospheric correction↗
DOE OSTI2024

Evaluation of Autoconversion Representation in E3SMv2 Using an Ensemble of Large-Eddy Simulations of Low-Level Warm Clouds

In numerical atmospheric models that treat cloud and rain droplet populations as separate condensate categories, precipitation initiation in warm clouds is often represented by an autoconversion rate (Au), which is the rate of formation of new rain droplets through the collisions of cloud droplets. Being a function of the cloud droplet size distribution (DSD), the local Au is commonly parameterized as a function of DSD moments: cloud droplet number (n c ) and mass (q c ) concentrations. When applied in a large-scale model, the grid-mean Au must also include a correction, or enhancement factor, to account for the horizontal variability of the cloud properties across the model grid. In this study, we evaluate the Au representation in the Energy Exascale Earth System Model version 2 (E3SMv2) climate model using large-eddy simulations (LES), which explicitly resolve cloud droplet spectra, and therefore the local Au, as well as its spatial variability. The analysis of an ensemble of warm low-level cloud cases shows that the E3SMv2 formulation represents the Au reasonably well compared to the horizontally averaged explicitly computed rate from LES. The agreement, however, comes from a combination of an underestimated E3SM-tuned local Au rate and an overestimated subgrid cloud variability enhancement factor. The latter bias is traced to neglecting the horizontal variability of n c and its co-variability with q c in parameterizing the grid-mean Au.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

The GEWEX Process Evaluation Study: GEWEX-PROES

Progress in realistically simulating the energy and water cycles in weather and climate models has been slower than desirable. Climate models submitted to the Coupled Model Intercomparison Project Phase-5 (CMIP-5), while more comprehensive than their predecessors, have shown little to no improvement in their biases in simulating key features of the climate system. As a result, uncertainties in global climate and hydrological sensitivities, and in the simulation of regional climate change, have not been reduced significantly. It is time to make use of these opportunities to significantly advance our understanding of key energy and water cycle processes at a wide range of space and time scales, and to provide a more insightful evaluation of the representation of these processes within models. This will require new ways of both analyzing the observations and diagnosing model behavior. The keys to success will be in skillfully combining different data sets and exploring relationships between them, as well as in the ability of models to reproduce those relationships correctly. In response to this challenge, a new GEWEX-wide activity, the GEWEX Process Evaluation Study, will take advantage of the opportunities that the combination of many of the existing data sets provide.

climate models↗
DOE OSTI2023

Core-collapse supernovae in the Dark Energy Survey: luminosity functions and host galaxy demographics

ABSTRACT We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically confirmed supernovae over a redshift range 0.045 < z < 0.25. We fit the observed DES griz CCSN light curves and K-correct to produce rest-frame R-band light curves. We compare the sample with lower redshift CCSN samples from Zwicky Transient Facility (ZTF) and Lick Observatory Supernova Search (LOSS). Comparing luminosity functions, the DES and ZTF samples of SNe II are brighter than that of LOSS with significances of 3.0σ and 2.5σ, respectively. While this difference could be caused by redshift evolution in the luminosity function, simpler explanations such as differing levels of host extinction remain a possibility. We find that the host galaxies of SNe II in DES are on average bluer than in ZTF, despite having consistent stellar mass distributions. We consider a number of possibilities to explain this – including galaxy evolution with redshift, selection biases in either the DES or ZTF samples, and systematic differences due to the different photometric bands available – but find that none can easily reconcile the differences in host colour between the two samples and thus its cause remains uncertain.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Laboratory instrumentation and techniques for characterizing multi-junction solar cells for space applications

Multi-junction solar cells are attractive for space applications because they can be designed to convert a larger fraction of AMO into electrical power at a lower cost than single-junction cells. The performance of multi-junction cells is much more sensitive to the spectral irradiance of the illuminating source than single-junction cells. The design of high efficiency multi-junction cells for space applications requires matching the optoelectronic properties of the junctions to AMO spectral irradiance. Unlike single-junction cells, it is not possible to carry out quantum efficiency measurements using only a monochromatic probe beam and determining the cell short-circuit current assuming linearity of the quantum efficiency. Additionally, current-voltage characteristics can not be calculated from measurements under non-AMO light sources using spectral-correction methods. There are reports in the literature on characterizing the performance of multi junction cells by measuring and convoluting the quantum efficiency of each junction with the spectral irradiance; the technique is of limited value for the characterization of cell performance under AMO power-generating conditions. We report the results of research to develop instrumentation and techniques for characterizing multi junction solar cells for space . An integrated system is described which consists of a standard lamp, spectral radiometer, dual-source solar simulator, and personal computer based current-voltage and quantum efficiency equipment. The spectral radiometer is calibrated regularly using the tungsten-halogen standard lamp which has a calibration based on NIST scales. The solar simulator produces the light bias beam for current-voltage and cell quantum efficiency measurements. The calibrated spectral radiometer is used to 'fit' the spectral irradiance of the dual-source solar simulator to WRL AMO data. The quantum efficiency apparatus includes a monochromatic probe beam for measuring the absolute cell quantum efficiency at various voltage biases, including the voltage bias corresponding to the maximum-power point under AMO light bias. The details of the procedures to 'fit' the spectral irradiance to AMO will be discussed. An assessment of the role of the accuracy of the 'fit' of the spectral irradiance and probe beam intensity on measured cell characteristics will be presented. quantum efficiencies were measured with both spectral light bias and AMO light bias; the measurements show striking differences. Spectral irradiances were convoluted with cell quantum efficiencies to calculate cell currents as function of voltage. The calculated currents compare with measured currents at the 1% level. Measurements on a variety of multi-junction cells will be presented. The dependence of defects in junctions on cell quantum efficiencies measured under light and voltage bias conditions will be presented. Comments will be made on issues related to standards for calibration, and limitations of the instrumentation and techniques. Expeditious development of multi-junction solar cell technology for space presents challenges for cell characterization in the laboratory.

Woodyard, James R.↗
NASA NTRS2021

Cloud Detection over Sunglint Regions with Observations from the Earth Polychromatic Imaging Camera (EPIC)

With the ability to observe the entire sunlit side of the Earth, EPIC data have become an important resource for studying cloud daily variability. Inaccurate cloud masking is a great source of uncertainty. One main region that is prone to error in cloud masking is the sunglint area over ocean surfaces. Cloud detection over these regions is challenging for the EPIC instrument because of its limited spectral channels. Clear sky ocean surface reflectance from visible channels over sunglint is much larger than that over the non-glint areas and can exceed reflectance from thin clouds. This paper presents an improved EPIC ocean cloud masking algorithm (Version 3). Over sunglint regions (glint angle ≤ 25°), the algorithm utilizes EPIC’s oxygen (O2) A-band ratio (764nm to 780 nm) in addition to the 780 nm reflectance observations in masking tests. Outside the sunglint regions, a dynamic reflectance threshold for the Rayleigh corrected 780 nm reflectance is applied. The thresholds are derived as a function of glint angle. When compared with co-located data from the geosynchronous Earth orbit (GEO) and the low Earth orbit (LEO) observations, the consistency of the new ocean cloud mask algorithm has increased by 4~10% and 4~6% in the glint center and granule edges respectively. The false positive rate is reduced by 10~17%. Overall global ocean cloud detection consistency increases by 2%. This algorithm, along with other improvements to the EPIC cloud masks, has been implemented in the EPIC cloud products Version 3. This algorithm will improve the cloud daily variability analysis by removing the artificial peak at local noon time in the glint center latitudes and reducing biases in the early morning and late afternoon cloud fraction over ocean surfaces.

EPIC↗
DOE OSTI2023

Near-term tropical cyclone risk and coupled Earth system model biases

Most current climate models predict that the equatorial Pacific will evolve under greenhouse gas–induced warming to a more El Niño-like state over the next several decades, with a reduced zonal sea surface temperature gradient and weakened atmospheric Walker circulation. Yet, observations over the last 50 y show the opposite trend, toward a more La Niña-like state. Recent research provides evidence that the discrepancy cannot be dismissed as due to internal variability but rather that the models are incorrectly simulating the equatorial Pacific response to greenhouse gas warming. This implies that projections of regional tropical cyclone activity may be incorrect as well, perhaps even in the direction of change, in ways that can be understood by analogy to historical El Niño and La Niña events: North Pacific tropical cyclone projections will be too active, North Atlantic ones not active enough, for example. Other perils, including severe convective storms and droughts, will also be projected erroneously. While it can be argued that these errors are transient, such that the models’ responses to greenhouse gases may be correct in equilibrium, the transient response is relevant for climate adaptation in the next several decades. Given the urgency of understanding regional patterns of climate risk in the near term, it would be desirable to develop projections that represent a broader range of possible future tropical Pacific warming scenarios—including some in which recent historical trends continue—even if such projections cannot currently be produced using existing coupled earth system models.

Science & Technology - Other Topics↗
NASA NTRS1999

Atmospheric Multiple Scattering Effects on GLAS Altimetry: Calculations of Single Pulse Bias - Part 1

Estimates of the effect of pulse stretching on satellite laser altimetry in particular the Geoscience Laser Altimeter System (GLAS), by cloud multiple scattering were made from an analytical method and from Monte Carlo calculations. The path delay of the return pulse was found to be largest for low-level clouds with particle radii (3-20 microns). The magnitude of the path delay was affected by several factors including cloud height, cloud optical depth, cloud particle size, particle shape, and receiver field of view. Polar aerosol and Rayleigh scattering usually added less than 1 cm to the overall path delay. Path delay estimates for all cloud conditions would be less if a simple Gaussian fit of the return pulse peak were used to measure the pulse's centroid. However, care must be taken in determining the centroid as factors such as pulse width, surface slope and the fitting method used will affect the estimate. A planned application for laser altimetry is high precision monitoring of the height change of polar ice sheets. In the absence of a correction algorithm, the required GLAS altimetry accuracies will not be achieved unless atmospheric channel information is used to remove profiles that are likely to be heavily contaminated by multiple scattering.

Duda, David P.↗
DOE OSTI2020

A Chip-Scale Oscillation-Mode Optomechanical Inertial Sensor Near the Thermodynamical Limits

Modern navigation systems integrate the global positioning system (GPS) with an inertial navigation system (INS), which complement each other for correct attitude and velocity determination. Here, the core of the INS integrates accelerometers and gyroscopes used to measure forces and angular rate in the vehicular inertial reference frame. With the help of gyroscopes and by integrating the acceleration to compute velocity and distance, precision and compact accelerometers with sufficient accuracy can provide small–error location determination. Solid–state implementations, through coherent readout, can provide a platform for high performance acceleration detection. In contrast to prior accelerometers using piezoelectric or capacitive readout techniques, optical readout provides narrow–linewidth high–sensitivity laser detection along with low–noise resonant optomechanical transduction near the thermodynamical limits. Here an optomechanical inertial sensor with an 8.2 µg Hz –1/2 velocity random walk (VRW) at an acquisition rate of 100 Hz and 50.9 µg bias instability is demonstrated, suitable for applications, such as, inertial navigation, inclination sensing, platform stabilization, and/or wearable device motion detection. Driven into optomechanical sustained–oscillation, the slot photonic crystal cavity provides radio–frequency readout of the optically–driven transduction with an enhanced 625 µg Hz –1 sensitivity. Measuring the optomechanically–stiffened oscillation shift, instead of the optical transmission shift, provides a 220× VRW enhancement over pre–oscillation mode detection.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Polarized and unpolarized gluon PDFs: Generative machine learning applications for lattice QCD matrix elements at short distance and large momentum

Lattice quantum chromodynamics (QCD) calculations share a defining challenge by requiring a small finite range of spatial separation z between quark/gluon bilinears for controllable power corrections in the perturbative QCD factorization, and a large hadron boost p z for a successful determination of collinear parton distribution functions (PDFs). However, these two requirements make the determination of PDFs from lattice data very challenging. We present the application of generative machine learning algorithms to estimate the polarized and unpolarized gluon correlation functions utilizing short-distance data and extending the correlation up to z p z ≲ 14 , surpassing the current capabilities of lattice QCD calculations. We train physics-informed machine learning algorithms to learn from the short-distance correlation at z ≲ 0.36 fm and take the limit, p z → ∞ , thereby minimizing possible contamination from the higher-twist effects for a successful reconstruction of the polarized gluon PDF. We also expose the bias and problems with underestimating uncertainties associated with the use of model-dependent and overly constrained functional forms, such as x α ( 1 − x ) β and its variants to extract PDFs from the lattice data. We propose the use of generative machine learning algorithms to mitigate these issues and present our determination of the polarized and unpolarized gluon PDFs in the nucleon. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Detector alignment for X-ray crystallography using Millepede-II

I describe a method for accurately refining the geometrical parameters of segmented X-ray area detectors on the basis of serial crystallography data, using 'Millepede' – an algorithm created for a very similar problem in high-energy physics. The Millepede method for serial crystallography builds on the approach of Brewster et al. [Acta Cryst. (2018), D74, 877–894], in which the detector parameters are refined simultaneously with the parameters for each individual crystal. This accounts for the mutual dependency between the parameters and thereby avoids the bias and slow convergence problems that have afflicted older approaches in which the deviations between observed and calculated Bragg peak positions were taken directly as the updates for the detector panel positions. The Millepede method uses the special structure of the least-squares normal equations to reduce them to a much smaller form that can be solved very quickly, even compared with the sparse matrix methods used previously. This makes it practical to refine the detector geometry frequently and thereby maintain accurate calibration without specialized alignment campaigns. Tilts of detector panels out of the plane can be reliably refined, as can the overall distance of the detector in the beam direction. With a simulated test case, the new method produced panel shifts within 7% of the correct values with only one iteration, and produced almost exactly correct shifts after a second iteration. A simulated out-of-plane panel rotation was correctly determined to within 0.001°. Applied to experimental data from an X-ray free-electron laser, the method increased the indexable fraction of frames from 30% to 91% in a single iteration, and to 96% after two further iterations. Computing the geometry updates on the basis of 2060 crystals took only 0.819 s on desktop computing hardware, including the time taken to read the required data from disk. The scaling was found to be very close to linear for up to 100 980 sets of crystal parameters, which took only 78.2 s to process under the same conditions. The method has been applied as part of a real-time feedback system at a synchrotron radiation beamline, in which an out-of-plane detector tilt of 0.04° was detected and corrected. Possible further applications are also described here.

Millepede-II↗
NASA NTRS2015

The Global Structure of UTLS Ozone in GEOS-5: A Multi-Year Assimilation of EOS Aura Data

Eight years of ozone measurements retrieved from the Ozone Monitoring Instrument (OMI) and the Microwave Limb Sounder, both on the EOS Aura satellite, have been assimilated into the Goddard Earth Observing System version 5 (GEOS-5) data assimilation system. This study thoroughly evaluates this assimilated product, highlighting its potential for science. The impact of observations on the GEOS-5 system is explored by examining the spatial distribution of the observation-minus-forecast statistics. Independent data are used for product validation. The correlation coefficient of the lower-stratospheric ozone column with ozonesondes is 0.99 and the bias is 0.5%, indicating the success of the assimilation in reproducing the ozone variability in that layer. The upper-tropospheric assimilated ozone column is about 10% lower than the ozonesonde column but the correlation is still high (0.87). The assimilation is shown to realistically capture the sharp cross-tropopause gradient in ozone mixing ratio. Occurrence of transport-driven low ozone laminae in the assimilation system is similar to that obtained from the High Resolution Dynamics Limb Sounder (HIRDLS) above the 400 K potential temperature surface but the assimilation produces fewer laminae than seen by HIRDLS below that surface. Although the assimilation produces 5 - 8 fewer occurrences per day (up to approximately 20%) during the three years of HIRDLS data, the interannual variability is captured correctly. This data-driven assimilated product is complementary to ozone fields generated from chemistry and transport models. Applications include study of the radiative forcing by ozone and tracer transport near the tropopause.

GEOS-5↗
NASA NTRS2016

Dynamic Characterization of Galfenol (Fe81.6Ga18.4)

Galfenol has the potential to transform the smart materials industry by allowing for the development of multifunctional, load-bearing devices. One of the primary technical challenges faced by this development is the very limited experimental data on Galfenol's frequency-dependent response to dynamic stress, which is critically important for the design of such devices. This report details a novel and precise characterization of the constitutive behavior of polycrystalline Galfenol (Fe81.6Ga18.4) under quasi-static (1 Hz) and dynamic (4 to 1000 Hz) stress loadings. Mechanical loads are applied using a high-frequency load frame. Quasi-static minor and major hysteresis loop measurements of magnetic flux density and strain are presented for constant electromagnet currents (0 to 1.1 A) and constant magnetic fields 0 to 14 kA/m (0 to 180 Oe). The dynamic stress amplitude for minor and major loops is 2.88 and 31.4 MPa (418 and 4550 psi), respectively. Quasi-static material properties closely match published values for similar Galfenol materials. Quasi-static actuation responses are also measured and compared to quasi-static sensing responses; the high degree of reversibility seen in the comparison is consistent with published measurements and modeling results. Dynamic major and minor loops are measured for dynamic stresses up to 31 MPa (4496 psi) and 1 kHz, and the bias condition resulting in maximum, quasi-static sensitivity. Eddy current effects are quantified by considering solid and laminated Galfenol rods. Three key sources of error in the dynamic measurements are accounted for: (1) electromagnetic noise in strain signals due to Galfenol's magnetic response, (2) error in load signals due to the inertial force of fixturing, and (3) phase misalignment between signals due to conditioning electronics. For dynamic characterization, strain error is kept below 1.2 percent of full scale by wiring two collocated gauges in series (noise cancellation) and through leadwire weaving. Inertial force error is kept below 0.41 percent by measuring the dynamic force in the specimen using a nearly collocated piezoelectric load washer. The phase response of all conditioning electronics is explicitly measured and corrected for. In general, as frequency is increased, the sensing response becomes more linear because of an increase in eddy currents. As frequency increases above approximately 100 Hz, the elbow in the strain-versus-stress response disappears as the active (soft) regime stiffens toward the passive (hard) regime. Under constant-field conditions, the loss factors of the solid rod peak between 200 and 600 Hz, rather than exhibiting a monotonic increase. Compared to the solid rod, the laminated rod exhibits much slower increases in hysteresis with frequency, and its quasi-static behavior extends to higher frequencies. The elastic modulus of the laminated rod decreases between 100 and 300 Hz; this trend is currently unexplained.

Dynamic sensing↗
NASA NTRS2020

Synergistic Use of Hyperspectral UV-Visible OMI and Broadband Meteorological Imager MODIS Data for a Merged Aerosol Product

The retrieval of optimal aerosol datasets by the synergistic use of hyperspectral ultraviolet(UV)–visible and broadband meteorological imager (MI) techniques was investigated. The Aura Ozone Monitoring Instrument (OMI) Level 1B (L1B) was used as a proxy for hyperspectral UV–visible instrument data to which the Geostationary Environment Monitoring Spectrometer (GEMS) aerosol algorithm was applied. Moderate-Resolution Imaging Spectroradiometer (MODIS) L1B and dark target aerosol Level 2 (L2) data were used with a broadband MI to take advantage of the consistent time gap between the MODIS and the OMI. First, the use of cloud mask information from the MI infrared (IR) channel was tested for synergy. High-spatial-resolution and IR channels of the MI helped mask cirrus and sub-pixel cloud contamination of GEMS aerosol, as clearly seen in aerosol optical depth (AOD) validation with Aerosol Robotic Network (AERONET) data. Second, dust aerosols were distinguished in the GEMS aerosol-type classification algorithm by calculating the total dust confidence index (TDCI) from MODIS L1B IR channels. Statistical analysis indicates that the Probability of Correct Detection (POCD) between the forward and inversion aerosol dust models (DS) was increased from 72% to 94% by use of the TDCI for GEMS aerosol-type classification, and updated aerosol types were then applied to the GEMS algorithm. Use of the TDCI for DS type classification in the GEMS retrieval procedure gave improved single-scattering albedo (SSA) values for absorbing fine pollution particles (BC) and DS aerosols. Aerosol layer height (ALH) retrieved from GEMS was compared with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data, which provides high-resolution vertical aerosol profile information. The CALIOP ALH was calculated from total attenuated backscatter data at 1064 nm, which is identical to the definition of GEMS ALH. Application of the TDCI value reduced the median bias of GEMS ALH data slightly. The GEMS ALH bias approximates zero, especially for GEMS AOD values of>~0.4 and GEMS SSA values of<~0.95.Finally, the AOD products from the GEMS algorithm and MI were used in aerosol merging with the maximum-likelihood estimation method, based on a weighting factor derived from the standard deviation of the original AOD products. With the advantage of the UV–visible channel in retrieving aerosol properties over bright surfaces, the combined AOD products demonstrated better spatial data availability than the original AOD products, with comparable accuracy. Furthermore, pixel-level error analysis of GEMS AOD data indicates improvement through MI synergy.

aerosol↗
NASA NTRS2022

Low-Level Jets and the Convergence of Mars Data Assimilation Algorithms

Data assimilation is an increasingly popular technique in Mars atmospheric science, but its effect on the mean states of the underlying atmosphere models has not been thoroughly examined. The robustness of results to the choice of model and assimilation algorithm also warrants further study. We investigate these issues using two Mars general circulation models (MGCMs), with particular emphasis on zonal wind and temperature fields. When temperature retrievals from the Mars Global Surveyor Thermal Emission Spectrometer (TES) are assimilated into the U.K.-Laboratoire de M ́et ́eorologie Dynamique (UK-LMD) MGCM to create the Mars Analysis Correction Data Assimilation (MACDA) reanalysis, low-level zonal jets in the winter northern hemisphere shift equatorward and weaken relative to a free-running control simulation from the same MGCM. The Ensemble Mars Atmosphere Reanalysis System (EMARS) reanalysis, which is also based on TES temperature retrievals, also shows jet weakening (but less if any shifting) relative to a control simulation performed with the underlying Geophysical Fluid Dynamics Laboratory (GFDL) MGCM. Examining higher levels of the atmosphere, monthly mean three-dimensional temperature and zonal wind fields are in generally better agreement between the two reanalyses than between the two control simulations. In conjunction with information about the MGCMs’ physical parametrizations, intercomparisons between the various reanalyses and control simulations suggest that overall the EMARS control run is plausibly less biased (relative to the true state of the Martian atmosphere) than the MACDA control run. Implications for future observational studies are discussed. Plain Language Summary An increasingly popular way to study Martian weather and climate is to combine atmospheric temperature observations with a computer model (specifically, a global climate model). The process of combining model and observations is called “data assimilation”, and the resulting merged data set is called a “reanalysis”. One advantage of reanalyses is that they include variables (such as wind) that are not directly observed. For scientific and practical applications we want these variables to be reasonably accurate— however, it is not clear how well data assimilation algorithms compute them. Our study investigates this issue using two Mars reanalyses and two model simulations that do not assimilate temperature data. We focus on slowly-varying atmospheric phenomena (timescales from 10 Mars days to a season). Assimilating temperature data into two different global climate models changes the strength and/or spatial pattern of east-west winds at low altitudes. Furthermore, monthly mean three-dimensional temperature and east-west wind fields agree better between reanalyses than between non-assimilating model simulations. This suggests that the data assimilation process is basically successful. One non-assimilating model simulation has less realistic representations of atmospheric physical processes than the other—we argue that this plausibly gives it larger biases relative to the true state of the atmosphere.

Todd A Mooring↗
DOE OSTI2025

Uranium particle age dating, aggregation, and model age best estimators

We present important aspects of uranium particle age dating by Large-Geometry Secondary Ion Mass Spectrometry (LG-SIMS) that can introduce bias and increase model age uncertainties, especially for small, young, and/or low-enriched particles. This metrology is important for applications related to International Nuclear Safeguards. We explore influential factors related to model age estimation, including the effects of evolving surface chemistry on inter-element measurements of particles (e.g., Th and U), detector background, and aggregation methods using simulated and actual particle samples. We introduce a new model age estimator, called “mid68”, that supplements 95% confidence intervals, providing a “best estimate” and uncertainty about the most likely age. The mid68 estimator can be calculated using the Feldman and Cousins method or Bayesian methods and provides a value with a symmetric uncertainty that can be used for calculations and approximate aggregation of processed model age values when the raw data and correction factors are not available. For particles yielding low 230 Th counts amidst nonzero detector background, their underlying model age probability distributions are asymmetric, so the mid68 estimator provides additional robust information regarding the underlying model age likelihood. This study provides a comprehensive and timely examination of critical aspects of uranium particle age dating as more laboratories establish particle chronometry capabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

The Pilatus Unmanned Aircraft System for Lower Atmospheric Research

This paper presents details of the University of Colorado (CU) "Pilatus" unmanned research aircraft, assembled to provide measurements of aerosols, radiation and thermodynamics in the lower troposphere. This aircraft has a wingspan of 3.2 m and a maximum take-off weight of 25 kg, and it is powered by an electric motor to reduce engine exhaust and concerns about carburetor icing. It carries instrumentation to make measurements of broadband up- and downwelling shortwave and longwave radiation, aerosol particle size distribution, atmospheric temperature, relative humidity and pressure and to collect video of flights for subsequent analysis of atmospheric conditions during flight. In order to make the shortwave radiation measurements, care was taken to carefully position a high-quality compact inertial measurement unit (IMU) and characterize the attitude of the aircraft and its orientation to the upward-looking radiation sensor. Using measurements from both of these sensors, a correction is applied to the raw radiometer measurements to correct for aircraft attitude and sensor tilt relative to the sun. The data acquisition system was designed from scratch based on a set of key driving requirements to accommodate the variety of sensors deployed. Initial test flights completed in Colorado provide promising results with measurements from the radiation sensors agreeing with those from a nearby surface site. Additionally, estimates of surface albedo from onboard sensors were consistent with local surface conditions, including melting snow and bright runway surface. Aerosol size distributions collected are internally consistent and have previously been shown to agree well with larger, surface-based instrumentation. Finally the atmospheric state measurements evolve as expected, with the near-surface atmosphere warming over time as the day goes on, and the atmospheric relative humidity decreasing with increased temperature. No directional bias on measured temperature, as might be expected due to uneven heating of the sensor housing over the course of a racetrack pattern, was detected. The results from these flights indicate that the CU Pilatus platform is capable of performing research-grade lower tropospheric measurement missions.

de Boer, Gijs↗