Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calibration and fitting methods”

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 253 records · Page 14

Retrievals of Profiles of Fine And Coarse Aerosols Using Lidar And Radiometric Space Measurements

In couple of years we expect the launch of the CALIPSO lidar spaceborne mission designed to observe aerosols and clouds. CALIPSO will collect profiles of the lidar attenuated backscattering coefficients in two spectral wavelengths (0.53 and 1.06 microns). Observations are provided along the track of the satellite around the globe from pole to pole. The attenuated backscattering coefficients are sensitive to the vertical distribution of aerosol particles, their shape and size. However the information is insufficient to be mapped into unique aerosol physical properties and vertical distribution. Infinite number of physical solutions can reconstruct the same two wavelength backscattered profile measured from space. CALIPSO will fly in formation with the Aqua satellite and the MODIS spectro-radiometer on board. Spectral radiances measured by MODIS in six channels between 0.55 and 2.13 microns simultaneously with the CALIPSO observations can constrain the solutions and resolve this ambiguity, albeit under some assumptions. In this paper we describe the inversion method and apply it to aircraft lidar and MODIS data collected over a dust storm off the coast of West Africa during the SHADE experiment. It is shown that the product of the single scattering albedo, omega, and the phase function, P, for backscattering can be retrieved from the synergism between measurements avoiding a priori hypotheses required for inverting lidar measurements alone. The resultant value of (omega)P(180 deg.) = 0.016/sr are significantly different from what is expected using Mie theory, but are in good agreement with recent results obtained from lidar observations of dust episodes. The inversion is robust in the presence of noise of 10% and 20% in the lidar signal in the 0.53 and 1.06 pm channels respectively. Calibration errors of the lidar of 5 to 10% can cause an error in optical thickness of 20 to 40% respectively in the tested cases. The lidar calibration errors cause degradation in the ability to fit the MODIS data. Therefore the MODIS measurements can be used to identify the calibration problem and correct for it. The CALIPSO-MODIS measurements of the profiles of fine and coarse aerosols, together with CALIPSO measurements of clouds vertical distribution, is expected to be critically important in understanding aerosol transport across continents and political boundaries, and to study aerosol-cloud interaction and its effect on precipitation and global forcing of climate.

Kaufman, Yoram↗

Long-wavelength Radar Studies of the Lunar Maria

Radar measurements at 70 cm and 7.5 m wavelengths provide insight into the structure and chemical properties of the upper 5-100 m of the lunar regolith and crust. Past work has identified a number of anomalous regions and changes in echo strength, some attributed to differences in titanium content. There has been little opportunity, however, to compare calibrated long-wavelength backscatter among different units or to theoretical model results. We combine recent high-resolution (3-5 km) 70-cm radar data for the nearside with earlier calibrated full-disk observations to provide a reasonable estimate of the true lunar backscatter coefficient. These data are tested against models for quasi-specular scattering from the surface, echoes from a buried substrate, and Mie scattering from surface and buried rocks. We find that 70 cm echoes likely arise from Mie scattering by distributed rocks within the soil, consistent with earlier hypotheses. Returns from a buried substrate would provide a plausible fit to the observations only if the regolith depth were ~3 m or less and varied little across the maria. Depolarized echoes are due to some combination of single and multiple scattering events, but it appears that single scattering alone could account for the observed echo power, based on comparisons with terrestrial rocky surfaces. Backscatter strength from the regolith is most strongly affected by the loss tangent, whose variation with mineral content is still poorly defined. We compared the backscatter values for the mare deposits to the oxide contents inferred from spectral ratio methods, and found that in general the unit boundaries evident in radar images closely follow those seen in color difference images. The 70-cm data are not well correlated with TiO2 values found using the Charette relationship nor with Fe abundances derived from Clementine observations. The lack of a relationship between radar echo and Fe content is reasonable given the distribution of iron among various mineral phases, but ilmenite content (FeTiO3) has typically been cited as the dominant cause of changes in loss tangent (and thus the radar absorption). The lack of correlation between the radar data and TiO2 estimates may arise from uncertainties in the Charette technique, subtle differences in the upper surface and bulk properties of the regolith, mineralogic effects on the radar not linked to titanium content, or to some combination of these factors. Dark crater haloes in the mare and highlands, and low radar returns from apparent cryptomare regions, are used to illustrate the role radar data can play in identifying changes in regolith composition; low-return haloes around craters such as Petavius may indicate 5-25% contamination of the highlands soil by excavated mare material or a layer of rock-poor ejecta at least several meters deep. The 7.5-m data were shown to correlate to a reasonable degree with estimates of Fe abundance, suggesting that this component of the mare basalts is primarily responsible for attenuation losses at very long wavelengths. The different sensitivities of the two radar wavelengths and multispectral data offers the potential for future deep mapping of the mare lava flows and regolith.

Lunar MariaRadarLunar Studies↗

Long-wavelength Radar Studies of the Lunar Maria

Radar measurements at 70 cm and 7.5 m wavelengths provide insight into the structure and chemical properties of the upper 5-100 m of the lunar regolith and crust. Past work has identified a number of anomalous regions and changes in echo strength, some attributed to differences in titanium content. There has been little opportunity, however, to compare calibrated long-wavelength backscatter among different units or to theoretical model results. We combine recent high-resolution (3-5 km) 70-cm radar data for the nearside with earlier calibrated full-disk observations to provide a reasonable estimate of the true lunar backscatter coefficient. These data are tested against models for quasi-specular scattering from the surface, echoes from a buried substrate, and Mie scattering from surface and buried rocks. We find that 70 cm echoes likely arise from Mie scattering by distributed rocks within the soil, consistent with earlier hypotheses. Returns from a buried substrate would provide a plausible fit to the observations only if the regolith depth were approx.3 m or less and varied little across the maria. Depolarized echoes are due to some combination of single and multiple scattering events, but it appears that single scattering alone could account for the observed echo power, based on comparisons with terrestrial rocky surfaces. Backscatter strength from the regolith is most strongly affected by the loss tangent, whose variation with mineral content is still poorly defined. We compared the backscatter values for the mare deposits to the oxide contents inferred from spectral ratio methods, and found that in general the unit boundaries evident in radar images closely follow those seen in color difference images. The 70-cm data are not well correlated with TiO2 values found using the Charette relationship nor with Fe abundances derived from Clementine observations. The lack of a relationship between radar echo and Fe content is reasonable given the distribution of iron among various mineral phases, but ilmenite content (FeTiO3) has typically been cited as the dominant cause of changes in loss tangent (and thus the radar absorption). The lack of correlation between the radar data and TiO2 estimates may arise from uncertainties in the Charette technique, subtle differences in the upper surface and bulk properties of the regolith, mineralogic effects on the radar not linked to titanium content, or to some combination of these factors. Dark crater haloes in the mare and highlands, and low radar returns from apparent cryptomare regions, are used to illustrate the role radar data can play in identifying changes in regolith composition; low-return haloes around craters such as Petavius may indicate 5-25% contamination of the highlands soil by excavated mare material or a layer of rock-poor ejecta at least several meters deep. The 7.5-m data were shown to correlate to a reasonable degree with estimates of Fe abundance, suggesting that this component of the mare basalts is primarily responsible for attenuation losses at very long wavelengths. The different sensitivities of the two radar wavelengths and multispectral data offers the potential for future deep mapping of the mare lava flows and regolith.

Lunar Studies↗

Quasi-Static and Dynamic Analysis of Composite Panels: Challenges in Engineering-Level Commercial Software

In engineering applications, where analysis time and efficiency are essential, analytical models and commercially available finite element (FE) material models are the two mainstream methods used to analyze composite crush scenarios. A reliable numerical model addressing the crushing response of composite materials requires a nonlinear progressive damage analysis which includes appropriate failure mechanisms such as matrix cracking, fragmentation, delamination, and fiber breakage. One might be able to address all these interacting failure responses through implementation/validation of a user defined material model or by utilization of a commercially available model, the latter being the most efficient engineering approach. In this paper, the quasi-static and dynamic crushing response of composite specimens, consisting of a C-channel and a corrugated geometry manufactured with carbon fiber/epoxy matrix fabric, are modeled using the progressive composite damage model within ABAQUS. The novel part of this work, which has not been previously addressed in the literature, is the utilization of a detailed parametric study to optimize composite damage model parameters for a given quasi-static crush simulation in ABAQUS. C-channel experiments were used to calibrate the model and then refine/optimize a set of material parameters that best fit the experimental results. These identical parameters were then used to simulate the crush response of the corrugated channel specimens, thereby making it a pure prediction and not a double recalibration. Next, the dynamic response of both panels were predicted, further confirming the predictive capability of the developed numerical model. The effect of parameter variation on the crush response is fully documented for future use in engineering applications utilizing this specific ABAQUS damage model.

Paria Naghipour↗

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗

Stellar Mass Calibrations for Local Low-mass Galaxies

The stellar masses of galaxies are measured from integrated light via several methods—however, few of these methods were designed for low-mass (M ⋆ ≲ 10 8 M ⊙ ) “dwarf” galaxies, whose properties (e.g., stochastic star formation, low metallicity) pose unique challenges for estimating stellar masses. In this work, we quantify the precision and accuracy at which stellar masses of low-mass galaxies can be recovered using UV/optical/IR photometry. We use mock observations of 469 low-mass galaxies from a variety of models, including both semi-empirical models (GRUMPY and UniverseMachine-SAGA) and cosmological baryonic zoom-in simulations (MARVELous Dwarfs and FIRE-2), to test literature color–M ⋆ /L relations and multiwavelength spectral energy distribution (SED) mass estimators. We identify a list of “best practices” for measuring stellar masses of low-mass galaxies from integrated photometry. We find that literature color–M ⋆ /L relations are often unable to capture the bursty star formation histories (SFHs) of low-mass galaxies, and we develop an updated prescription for stellar mass based on g − r color that is better able to recover stellar masses for the bursty low-mass galaxies in our sample (with ∼0.1 dex precision). SED fitting can also precisely recover stellar masses of low-mass galaxies, but this requires thoughtful choices about the form of the assumed SFH: Parametric SFHs can underestimate stellar mass by as much as ∼0.4 dex, while nonparametric SFHs recover true stellar masses with insignificant offset (−0.03 ± 0.11 dex). Finally, we also caution that noninformative (wide) dust attenuation priors may introduce M ⋆ uncertainties of up to ∼0.6 dex.

de los Reyes, Mithi A. C. [Amherst College, MA (Un↗

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗

Design of low loss helix circuits for interference fitted and brazed circuits

The RF loss properties and thermal capability of brazed helix circuits and interference fitted circuits were evaluated. The objective was to produce design circuits with minimum RF loss and maximum heat transfer. These circuits were to be designed to operate at 10 kV and at 20 GHz using a gamma a approximately equal to 1.0. This represents a circuit diameter of only 0.75 millimeters. The fabrication of this size circuit and the 0.48 millimeter high support rods required considerable refinements in the assembly techniques and fixtures used on lower frequency circuits. The transition from the helices to the waveguide was designed and the circuits were matched from 20 to 40 GHz since the helix design is a broad band circuit and at a gamma a of 1.0 will operate over this band. The loss measurement was a transmission measurement and therefore had two such transitions. This resulting double-ended match required tuning elements to achieve the broad band match and external E-H tuners at each end to optimize the match for each frequency where the loss measurement was made. The test method used was a substitution method where the test fixture was replaced by a calibrated attenuator.

Jacquez, A.↗

Binning is Sinning (Supernova Version): The Impact of Self-calibration in Cosmological Analyses with Type Ia Supernovae

Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia supernovae (SNe Ia) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic uncertainty vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of ~ 1.5 × with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach indeed yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is solely degenerate with the cosmological model in redshift space or whether it can be described by additional correlations between supernova properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the data set, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically classified samples, we must first address the issue that binning is required in currently used photometric analysis methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Long-Term Trends in the Concentrations of SF6, CHClF2, and COF2 in the Lower Stratosphere from Analysis of High-Resolution Infrared Solar Occultation Spectra

Long-term trends in the concentrations of SF6, CHClF2 (CFC-22), and COF2 in the lower stratosphere have been derived from analysis of ca. 1980 and more recent infrared solar occultation spectra recorded near 32 deg N latitude at approx. 0.02/ cm resolution. Consistent sets of line parameters and spectral calibration methods have been used in the retrievals to minimize systematic error effects. Quoted error limits are 1 sigma estimated precisions. The SF6 and CHClF2 results are based on spectra recorded by balloon-borne interferometers in March 1981 and June 1988 and a comparison of these results with the Atmospheric Trace Molecule Spectroscopy (ATMOS) Experiment/Spacelab 3 measurements obtained in May 1985 near 30 deg N latitude. In the 13-18 km altitude range the mean measured SF6 mixing ratio in parts per trillion by volume (pptv) increased from 1.17 +/- 0.21 in March 1981 to 2.02 +/- 0.20 pptv in June 1988, and the CHClF2 mixing ratio below 15 km altitude increased from 51 +/- 8 pptv in March 1981 to 102 +/- 10 pptv in June 1988. The CHClF2 retrievals used new empirical CHClF2 line parameters derived from 0.03/cm resolution laboratory spectra recorded at six temperatures between 203 and 293 K; the derived mixing ratios are approx. 30% higher than obtained with earlier sets of line parameters, thereby removing a large discrepancy noted previously between IR and in situ measurements of CHClF2. Assuming an exponential growth model for fitting the trends, SF6 and CHClF2 mean increase rates of 7.4% +/- 1.9% and 9.4% +/- 1.3% /year, are obtained, respectively, which correspond to cumulative increases by factors of approx. 1.7 and -2.0 in the concentrations of these gases over the 7.2-year measurement period. Analysis of spectra recorded in October 1979 and April 1989 yields COF2 volume mixing ratios that are respectively 0.44 +/- 0.17 and 1.21 +/- 0.24 times the ATMOS/Spacelab 3 values, from which an average COF2 increase rate of 10.3 +/- 1.8%/ year over this time period has been estimated. The present results are compared with previously reported observations and trends and with one-dimensional model calculations. The model calculated trends are in reasonably good agreement with the observations.

Rinsland, C. P.↗

Effects of Dose Error and Sample Size on Sonic Boom Dose-Response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

Encoding nonlinear and unsteady aerodynamics of limit cycle oscillations using nonlinear sparse Bayesian learning

This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.

97 MATHEMATICS AND COMPUTING↗

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗

X-Ray Spectro-Polarimetry with Photoelectric Polarimeters

We derive a generalization of forward fitting for X-ray spectroscopy to include linear polarization of X-ray sources, appropriate for the anticipated next generation of space-based photoelectric polarimeters. We show that the inclusion of polarization sensitivity requires joint fitting to three observed spectra, one for each of the Stokes parameters, I(E), U(E), and Q(E). The equations for StokesI (E) (the total intensity spectrum) are identical to the familiar case with no polarization sensitivity, and for which the model-predicted spectrum is obtained by a convolution of the source spectrum, F (E), with the familiar energy response function,(E) R(E,E), where (E) and R(E,E) are the effective area and energy redistribution matrix, respectively. In addition to the energy spectrum, the two new relations for U(E) and Q(E) include the source polarization fraction and position angle versus energy, a(E), and 0(E), respectively, and the model-predicted spectra for these relations are obtained by a convolution with the modulated energy response function, (E)(E) R(E,E), where(E) is the energy-dependent modulation fraction that quantifies a polarimeters angular response to 100 polarized radiation. We present results of simulations with response parameters appropriate for the proposed PRAXyS Small Explorer observatory to illustrate the procedures and methods, and we discuss some aspects of photoelectric polarimeters with relevance to understanding their calibration and operation.

polarization↗

Second ARM Aerosol Chemical Speciation Monitor Users’ Meeting Report

The aerosol chemical speciation monitor (ACSM) was developed to adapt the technology of the aerosol mass spectrometer (AMS) to routine, long-term, standalone monitoring. The calculation of particulate mass concentration from ACSM data requires the measurement of the response of the instrument to aerosol of specific size and composition as well as assumptions about the instrument response based on laboratory measurements and field experience acquired over more than two decades of operation of AMS and a decade of operation of the ACSM. Three parameters in the concentration calculations that are particularly important are the NO3 response factor (RFNO3), relative ionization efficiency (RIE), and the collection efficiency (CE). The values of RFNO3 and RIEs are determined from calibration, however the jump calibration method previously used in calibration can result in errors in the RIE for sulfate. This has been corrected by implementing a continuous calibration method. The default collection efficiency is 0.5. This has been shown to result in mass loadings that do not agree with mass determinations from other instruments because of effects of composition on the vaporization of the particles. The previous work of investigators addressing this issue is discussed. After preliminary work on ACSM and scanning mobility particle sizer (SMPS) data from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) observatory collected in late 2016 and 2017 produced a parameterization of composition dependent collection efficiency very different from the results of previous studies, SMPS data were examined and we determined that there was significant mass that the instrument did not capture because the particles with diameters larger than 465 nm are not counted by this instrument. Data for the ultra-high-sensitivity aerosol spectrometer (UHSAS), SMPS, and ACSM are available for nearly all of 2019. The data from UHSAS and SMPS collected in 2019 were compared. We found that the UHSAS data has particle counts and total volumes significantly less than measured by the condensation particle counter (CPC) and SMPS. The SMPS data were extended by fitting the average volume distribution with a log normal curve and using this relationship to estimate a mass value for the SMPS over extended diameter range. The extended SMPS mass values result in a CDCE parameterization that is in better agreement with the results of other investigators, but is still different from other formulations. The working group recommends that the ACSM data be processed with a collective efficiency (CE)=1, that this be documented clearly in the metadata, and the use of the default CE of 0.5 or a formulation of composition-dependent collection efficiencies (CDCE) chosen by the user should be implemented based on the ammonium nitrate mass fraction. This is clearly necessary for the wintertime SGP ACSM data because of the high nitrate concentrations.

47 OTHER INSTRUMENTATION↗

Raman Spectroscopy as a Method for Mineral Identification on Lunar Robotic Exploration Missions

The sharp, nonoverlapping Raman bands for plagioclase, pyroxene, and olivine would be advantageous for on-surface, active mineralogical analysis of lunar materials. A robust, light-weight, low-power, rover-based Raman spectrometer with a laser exciting source, entirely transmission-mode holographic optics, and a charge-coupled device (CCD) detector could fit within a less than 20 cm cube. A sensor head on the end of an optical fiber bundle that carried the laser beam and returned the scattered radiation could be placed against surfaces at any desired angle by a deployment mechanism; otherwise, the instrument would need no moving parts. A modem micro-Raman spectrometer with its beam broadened (to expand the spot to 50-micrometer diameter) and set for low resolution (7/cm in the 100-1400/cm region relative to 514.5-nm excitation), was used to simulate the spectra anticipated from a rover instrument. We present spectra for lunar mineral grains, less than 1 mm soil fines, breccia fragments, and glasses. From frequencies of olivine peaks, we derived sufficiently precise forsterite contents to correlate the analyzed grains to known rock types and we obtained appropriate forsterite contents from weak signals above background in soil fines and breccias. Peak positions of pyroxenes were sufficiently well determined to distinguish among orthorhombic, monoclinic, and triclinic (pyroxenoid) structures; additional information can be obtained from pyroxene spectra, but requires further laboratory calibration. Plagioclase provided sharp peaks in soil fines and most breccias even when the glass content was high.

Wang, Alian↗

Drop Axis Ratio Distributions in Stratiform and Convective Rain

A fully calibrated low profile 2D video disdrometer (2DVD) has been recording many different rainfall events in Northern Alabama (USA) since June 2007. An earlier publication reported drop shapes and axis ratio distributions determined for some of the events. For one of the cases examined, a noticeable shift in the 3.5 - 3.75 mm drop axis ratio distribution was noted. In this paper, we extend the earlier work by separating the 2DVD measurements into stratiform and convective rain. The separation is made possible by using the minute-by-minute drop size distribution (DSD) measured by the 2DVD. The 1-minute DSDs are fitted to a gamma distribution, and using a simple indexing technique which involves two of the fitted parameters, periods of convective and stratiform rain are separated for a given event. The output of the DSD indexing technique is qualitatively confirmed by comparing with simultaneous time series observations from a co-located UHF profiler which continuously records height profiles of reflectivity, Doppler mean and spectral width, all of which enable the identification of bright-band periods and, furthermore, periods of moderate and deep convection. Excellent consistency is found between the output of the DSD-based separation method and the profiler observations. Next, we utilize the output of DSD index-based separation method to flag the periods of severe convection for a given event. Drop axis ratios during the flagged periods are derived and compared with those during stratiform rain periods. Five cases have been considered. Axis ratio distributions do not show appreciable differences between stratiform and convective periods for four of the cases. The fifth case (the same case as reported earlier) shows a shift in the 3.5 - 3.75 mm drop axis ratios during a prolonged period of convection. The contoured shapes for these drops determined from the 2DVD camera data indicate the possibility of non-axisymmetric oscillations, compared with the contoured images for other events which fit well to our reference drop shapes. For all of above cases, observations from a C-band polarimetric radar - situated 15 km away are examined. The variations between the co-polar radar reflectivity and the differential reflectivity as well as the specific differential phase are compared with the 2DVD data based scattering calculations for the 5 events. The implications will be discussed.

Thurai, M.↗

SeaWiFS Technical Report Series. Volume 29: SeaWiFS CZCS-type pigment algorithm

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) mission will provide operational ocean color that will be superior to the previous Coastal Zone Color Sensor (CZCS) proof-of-concept mission. an algorithm is needed that exploits the full functionality of SeaWiFS whilst remaining compatible in concept with algorithms used for the CZCS. This document describes the theoretical rationale of radiance band-radio methods for determining chlorophyll alpha and other important biogeochemical parameters, and their implementation for the SeaWiFS mission. Pigment interrelationships are examined to explain the success of the CZCS algorithms. In the context where chlorophyll alpha absorbs only weakly at 520 nm, the success of the 520 nm to 550 nm CZCS band ratio needs to be explained. This is explained by showing that in pigment data from a range of oceanic provinces chlorophyll alpha (absorbing at less than 490 nm), carotenoids (absorbing at greater than 460 nm), and total pigment are highly correlated. Correlations within pigment groups particularly photoprotectant and photosynthetic carotenoids are less robust. The sources of variability in optical data re examined using the NIMBUS Experiment Team (NET) bio-optical data set and bio-optical model. In both the model and NET data, the majority of the variance in the optical data is attributed to variability in pigment (chlorophyll alpha, and total particulates, with less than 5% of the variability resulting from pigment assemblage. The relationships between band ratios and chlorophyll is examined analytically, and a new formulation based on a dual hyperbolic model is suggested which gives a better calibration curve than the conventional log-log linear regression fit. The new calibration curve shows that 490:555 ratio is the best single-band ratio and is the recommended CZCS-type pigment algorithm. Using both the model and NET data, a number of multiband algorithms are developed; the best of which is an algorithm based on the 443:555 and 490:555 ratios. From model data, the form of potential algorithms for other products, such as total particulates and dissolved organic matter (DOM), are suggested.

Hooker, Stanford B.↗