Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “aggregation approach”

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 289 records · Page 16

Dissemination of Earth Remote Sensing Data for Use in the NOAA/NWS Damage Assessment Toolkit

The National Weather Service has developed the Damage Assessment Toolkit (DAT), an application for smartphones and tablets that allows for the collection, geolocation, and aggregation of various damage indicators that are collected during storm surveys. The DAT supports the often labor‐intensive process where meteorologists venture into the storm‐affected area, allowing them to acquire geotagged photos of the observed damage while also assigning estimated EF‐scale categories based upon their observations. Once the data are collected, the DAT infrastructure aggregates the observations into a server that allows other meteorologists to perform quality control and other analysis steps before completing their survey and making the resulting data available to the public. In addition to in‐person observations, Earth remote sensing from operational, polar-orbiting satellites can support the damage assessment process by identifying portions of damage tracks that may be missed due to road limitations, access to private property, or time constraints. Products resulting from change detection techniques can identify damage to vegetation and the land surface, aiding in the survey process. In addition, higher resolution commercial imagery can corroborate ground‐based surveys by examining higher‐resolution commercial imagery. As part of an ongoing collaboration, NASA and NOAA are working to integrate near real‐time Earth remote sensing observations into the NOAA/NWS Damage Assessment Toolkit. This presentation will highlight recent developments in a streamlined approach for disseminating Earth remote sensing data via web mapping services and a new menu interface that has been integrated within the DAT. A review of current and future products will be provided, including products derived from MODIS and VIIRS for preliminary track identification, along with conduits for higher‐resolution Landsat, ASTER, and commercial imagery as they become available. In addition to tornado damage assessments, the team is also investigating the use of near real‐time imagery for identifying hail damage to vegetation, which also results in large swaths of damage, particularly in the central United States during the peak growing season months of June, July, and August. This presentation will present an overview of recent activities, challenges and successes, best practices, and opportunities for future work and collaboration

Molthan, Andrew↗

Resolving Ice Cloud Optical Thickness Biases Between CALIOP and MODIS Using Infrared Retrievals

Despite its importance as one of the key radiative properties that determines the impact of upper tropospheric clouds on the radiation balance, ice cloud optical thickness (IOT) has proven to be one of the more challenging properties to retrieve from space-based remote sensing measurements. In particular, optically thin upper tropospheric ice clouds (cirrus) have been especially challenging due to their tenuous nature, extensive spatial scales, and complex particle shapes and light scattering characteristics. The lack of independent validation motivates the investigation presented in this paper, wherein systematic biases between MODIS Collection 5 (C5) and CALIOP Version 3 (V3) unconstrained retrievals of tenuous IOT (< 3) are examined using a month of collocated A-Train observations. An initial comparison revealed a factor of two bias between the MODIS and CALIOP IOT retrievals. This bias is investigated using an infrared (IR) radiative closure approach that compares both products with MODIS IR cirrus retrievals developed for this assessment. The analysis finds that both the MODIS C5 and the unconstrained CALIOP V3 retrievals are biased (high and low, respectively) relative to the IR IOT retrievals. Based on this finding, the MODIS and CALIOP algorithms are investigated with the goal of explaining and minimizing the biases relative to the IR. For MODIS we find that the assumed ice single scattering properties used for the C5 retrievals are not consistent with the mean IR COT distribution. The C5 ice scattering database results in the asymmetry parameter (g) varying as a function of effective radius with mean values that are too large. The MODIS retrievals have been brought into agreement with the IR by adopting a new ice scattering model for Collection 6 (C6) consisting of a modified gamma distribution comprised of a single habit (severely roughened aggregated columns); the C6 ice cloud optical property models have a constant g approx. = 0.75 in the mid-visible spectrum, 5-15% smaller than C5. For CALIOP, the assumed lidar ratio for unconstrained retrievals is fixed at 25 sr for the V3 data products.This value is found to be inconsistent with the constrained (predominantly nighttime) CALIOP retrievals. An experimental data set was produced using a modified lidar ratio of 32 sr for the unconstrained retrievals (an increase of 28%), selected to provide consistency with the constrained V3 results. These modifications greatly improve the agreement with the IR and provide consistency between the MODIS and CALIOP products. Based on these results the recently released MODIS C6 optical products use the single habit distribution given above, while the upcoming CALIOP V4 unconstrained algorithm will use higher lidar ratios for unconstrained retrievals.

Cloud↗

A Regional CO2 Observing System Simulation Experiment Using ASCENDS Observations and WRF-STILT Footprints

Knowledge of the spatiotemporal variations in emissions and uptake of CO2 is hampered by sparse measurements. The recent advent of satellite measurements of CO2 concentrations is increasing the density of measurements, and the future mission ASCENDS (Active Sensing of CO2 Emissions over Nights, Days and Seasons) will provide even greater coverage and precision. Lagrangian atmospheric transport models run backward in time can quantify surface influences ("footprints") of diverse measurement platforms and are particularly well suited for inverse estimation of regional surface CO2 fluxes at high resolution based on satellite observations. We utilize the STILT Lagrangian particle dispersion model, driven by WRF meteorological fields at 40-km resolution, in a Bayesian synthesis inversion approach to quantify the ability of ASCENDS column CO2 observations to constrain fluxes at high resolution. This study focuses on land-based biospheric fluxes, whose uncertainties are especially large, in a domain encompassing North America. We present results based on realistic input fields for 2007. Pseudo-observation random errors are estimated from backscatter and optical depth measured by the CALIPSO satellite. We estimate a priori flux uncertainties based on output from the CASA-GFED (v.3) biosphere model and make simple assumptions about spatial and temporal error correlations. WRF-STILT footprints are convolved with candidate vertical weighting functions for ASCENDS. We find that at a horizontal flux resolution of 1 degree x 1 degree, ASCENDS observations are potentially able to reduce average weekly flux uncertainties by 0-8% in July, and 0-0.5% in January (assuming an error of 0.5 ppm at the Railroad Valley reference site). Aggregated to coarser resolutions, e.g. 5 degrees x 5 degrees, the uncertainty reductions are larger and more similar to those estimated in previous satellite data observing system simulation experiments.

Wang, James S.↗

Uncertainty of Passive Imager Cloud Optical Property Retrievals to Instrument Radiometry and Model Assumptions: Examples from MODIS

The optical and microphysical structure of clouds is of fundamental importance for understanding a variety of cloud radiation and precipitation processes. With the advent of MODIS on the NASA EOS Terra and Aqua platforms, simultaneous global-daily 1 km retrievals of cloud optical thickness (COT) and effective particle radius (CER) are provided, as well as the derived water path (WP). The cloud product (MOD06/MYD06 for MODIS Terra and Aqua, respectively) provides separate retrieval datasets for various two-channel retrievals, typically a VISNIR channel paired with a 1.6, 2.1, and 3.7 m spectral channel. The MOD06 forward model is derived from on a homogeneous plane-parallel cloud. In Collection 5 processing (completed in 2007 with a modified Collection 5.1 completed in 2010), pixel-level retrieval uncertainties were calculated for the following non-3-D error sources: radiometry, surface spectral albedo, and atmospheric corrections associated with model analysis uncertainties (water vapor only). The latter error source includes error correlation across the retrieval spectral channels. Estimates of uncertainty in 1 aggregated (Level-3) means were also provided assuming unity correlation between error sources for all pixels in a grid for a single day, and zero correlation of error sources from one day to the next. I n Collection 6 (expected to begin in late summer 2013) we expanded the uncertainty analysis to include: (a) scene-dependent calibration uncertainty that depends on new band and detector-specific Level 1B uncertainties, (b) new model error sources derived from the look-up tables which includes sensitivities associated with wind direction over the ocean and uncertainties in liquid water and ice effective variance, (c) thermal emission uncertainties in the 3.7 m band associated with cloud and surface temperatures that are needed to extract reflected solar radiation from the total radiance signal, (d) uncertainty in the solar spectral irradiance at 3.7 m, and (e) addition of stratospheric ozone uncertainty in visible atmospheric corrections. A summary of the approach and example Collection 6 results will be shown.

Cloud Properties Observations↗

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗

Lunar Dust Mitigation: A Guide and Reference: First Edition (2021)

On the surface of the Moon, lunar dust specifically presents unique challenges to operations long-term due to regolith particles’ ubiquitous presence in the lunar environment, potential to electrostatically charge and possible chemical reactivity. Whether to avoid exposure, try to remove or simply to tolerate lunar dust infiltrating a system becomes a complex question involving length of required service life, dust effects and critical risks, mass and complexity trades for the entire system. This publication provides a snapshot of advice from topical experts for specific areas of concern to systems targeted for deployment on the lunar surface. Following introductory overview commentary, dust mitigation approaches appropriate to the lunar surface are first addressed for typically static structures such as optical surfaces, radiators, and other thermal control surfaces, followed by regolith exposure concerns for communications equipment and non-optical sensors. The broad topic of mechanisms and mechanical assemblies is broken down to address relevant component level concerns, such as for bearings and for seals. “Soft goods” components of space suits (fabric) specific concerns are discussed, followed finally by brief coverage of human health issues and concerns, though the emphasis of this publication remains with components directly exposed to the harsh lunar surface environment. A description of some lunar surface hazard details, in particular characteristics of lunar surface dust, follows in Appendix A, and more detailed explanations of quantification issues for particulates are in Appendix B. The simple inertial removal of particles from a surface is discussed in Appendix C, followed by a summary of terrestrial best practices for dust mitigation recommended within select industries in Appendix D. The aggregated bibliography of references, while extensive and very useful, should not be construed to be exhaustive and can serve as a constructive start.

dust↗

Object-Based Mapping of the Circumpolar Taiga-Tundra Ecotone with MODIS Tree Cover

The circumpolar taiga tundra ecotone was delineated using an image-segmentation-based mapping approach with multi-annual MODIS Vegetation Continuous Fields (VCF) tree cover data. Circumpolar tree canopy cover (TCC) throughout the ecotone was derived by averaging MODIS VCF data from 2000 to 2005 and adjusting the averaged values using linear equations relating MODIS TCC to Quickbird-derived tree cover estimates. The adjustment helped mitigate VCF's overestimation of tree cover in lightly forested regions. An image segmentation procedure was used to group pixels representing similar tree cover into polygonal features (segmentation objects) that form the map of the transition zone. Each polygon represents an area much larger than the 500 m MODIS pixel and characterizes the patterns of sparse forest patches on a regional scale. Those polygons near the boreal/tundra interface with either (1) mean adjusted TCC values from5 to 20%, or (2) mean adjusted TCC values greater than 5% but with a standard deviation less than 5% were used to identify the ecotone. Comparisons of the adjusted average tree cover data were made with (1) two existing tree line definitions aggregated for each 1 degree longitudinal interval in North America and Eurasia, (2) Landsat-derived Canadian proportion of forest cover for Canada, and (3) with canopy cover estimates extracted from airborne profiling lidar data that transected 1238 of the TCC polygons. The adjusted TCC from MODIS VCF shows, on average, less than 12% TCC for all but one regional zone at the intersection with independently delineated tree lines. Adjusted values track closely with Canadian proportion of forest cover data in areas of low tree cover. A comparison of the 1238 TCC polygons with profiling lidar measurements yielded an overall accuracy of 67.7%.

Ranson, K. J.↗

Ion loss processes in the stratosphere

Small ions consisting of aggregates of a few molecules determine th stratospheric electrical parameters such as conductivity, mobility, etc. The small ion density is controlled by the ionizing mechanisms for the production of ions and electrons and the loss processes for these charged particles. Ion production in stratosphere is chiefly due to galactic cosmic rays, and the loss processes due to recombination and attachment. Free electrons do not exist at stratospheric heights. The primary positive ion O2+ and the electrons are converted into complex clusters of positive and negative ions. The equilibrium ion density is governed by the equation of continuity for the production and loss of these ions. In the generalized equation of continuity, the gain or loss of ions due to transport is neglected. When aerosols or particulates are not present in the atmosphere, ions are loss due to ion-ion recombination. In the presence of aerosols, small ions can also be lost by attachment to the aerosols, and thus aerosols are likely to cause perturbations in the stratosphere electrification. A simplified model approach is adopted to study the effect of aerosols on the equilibrium ion density. The results of the analysis for the equatorial station Thumba (8.5 deg N) are presented.

Prasad, B. S. N.↗

Thermal and Hydrologic Signatures of Soil Controls on Evaporation: A Combined Energy and Water Balance Approach with Implications for Remote Sensing of Evaporation

The overall goal of this research is to examine the feasibility of applying a newly developed diagnostic model of soil water evaporation to large land areas using remotely sensed input parameters. The model estimates the rate of soil evaporation during periods when it is limited by the net transport resulting from competing effects of capillary rise and drainage. The critical soil hydraulic properties are implicitly estimated via the intensity and duration of the first stage (energy limited) evaporation, removing a major obstacle in the remote estimation of evaporation over large areas. This duration, or "time to drying" (t(sub d)), is revealed through three signatures detectable in time series of remote sensing variables. The first is a break in soil albedo that occurs as a small vapor transmission zone develops near the surface. The second is a break in either surface to air temperature differences or in the diurnal surface temperature range, both of which indicate increased sensible heat flux (and/or storage) required to balance the decrease in latent heat flux. The third is a break in the temporal pattern of near surface soil moisture. Soil moisture tends to decrease rapidly during stage 1 drying (as water is removed from storage), and then become more or less constant during soil limited, or "stage 2" drying (as water is merely transmitted from deeper soil storage). The research tasks address: (1) improvements in model structure, including extensions to transpiration and aggregation over spatially variable soil and topographic landscape attributes; and (2) applications of the model using remotely sensed input parameters.

Salvucci, Guido D.↗

Thermal and Hydrologic Signatures of Soil Controls on Evaporation: A Combined Energy and Water Balance Approach with Implications for Remote Sensing of Evaporation

The overall goal of this research is to examine the feasibility of applying a newly developed diagnostic model of soil water evaporation to large land areas using remotely sensed input parameters. The model estimates the rate of soil evaporation during periods when it is limited by the net transport resulting from competing effects of capillary rise and drainage. The critical soil hydraulic properties are implicitly estimated via the intensity and duration of the first stage (energy limited) evaporation, removing a major obstacle in the remote estimation of evaporation over large areas. This duration, or 'time to drying' (t(sub d)) is revealed through three signatures detectable in time series of remote sensing variables. The first is a break in soil albedo that occurs as a small vapor transmission zone develops near the surface. The second is a break in either surface to air temperature differences or in the diurnal surface temperature range, both of which indicate increased sensible heat flux (and/or storage) required to balance the decrease in latent heat flux. The third is a break in the temporal pattern of near surface soil moisture. Soil moisture tends to decrease rapidly during stage I drying (as water is removed from storage), and then become more or less constant during soil limited, or 'stage II' drying (as water is merely transmitted from deeper soil storage). The research tasks address: (1) improvements in model structure, including extensions to transpiration and aggregation over spatially variable soil and topographic landscape attributes; and (2) applications of the model using remotely sensed input parameters.

Salvucci, Guido D.↗

Metallic Vs Inflatable Primary Structure: A Comprehensive Trade for a Mars Transit Habitat

Large diameter inflatable habitat solutions are hypothesized to have an advantage over traditional metallic habitats of equal volume and functional requirements in terms of minimizing mass. A fair and comprehensive trade of metallic primary structure versus inflatable softgoods primary structure for a Mars Transit Habitat was conducted by The Boeing Company under the NASA Next Space Technologies for Exploration Partnerships-2 (NextSTEP-2)Appendix A contract in order to determine the minimum mass option from a system level. A Mars Transit Habitat is a major system mass driver because small increases in the mass can create large mass changes in the propulsive stage, adding complexity and cost to the mission. Three different vehicle configurations were traded using a fixed volume constraint and the same Ground Rules and Assumptions(GR&A)provided by NASA. The configurations were 1) 5.5 m diameter metallic habitat, 2) 7.0 m diameter metallic habitat, and 3) 7.63 m diameter inflatable habitat. In order to complete the trade, each configuration was analyzed and sized, but a system level approach had to be taken to ensure other aspects and subsystems besides the structure were captured. The primary functions (i.e. restraint layer, bladder, multilayer insulation [MLI], micro-meteoroid orbital debris [MMOD]) and strength advantages of softgoods are well understood, but there were other system level impacts that were discovered or addressed that are unique to softgoods. For the metallic structure, a modern and mass optimal material and manufacturing process was adopted to ensure the lightest weight options were evaluated. Two different aluminium lithium structures were evaluated, a longer length 5.5 m diameter and a shorter length 7 m diameter structure. The trade demonstrates that each configuration is feasible and close to the target dry mass for the Mars Transit Habitat Trans-Martian Injection burn. System level mass was the primary consideration, but since the three concepts were close in total mass, a clear winner did not arise. Secondary factors were evaluated such as architecture habitability (human factors), development considerations, on-orbit activation and outfitting, extensibility to other applications, and configuration flexibility. The aggregation of the lowest mass, most benefits, and fewest challenges resulted in the 5.5m diameter metallic configuration being the most favorable.

Matthew Ziglar↗

Quantification of human responses

Human perception is a complex phenomenon which is difficult to quantify with instruments. For this reason, large panels of people are often used to elicit and aggregate subjective judgments. Print quality, taste, smell, sound quality of a stereo system, softness, and grading Olympic divers and skaters are some examples of situations where subjective measurements or judgments are paramount. We usually express what is in our mind through language as a medium but languages are limited in available choices of vocabularies, and as a result, our verbalizations are only approximate expressions of what we really have in mind. For lack of better methods to quantify subjective judgments, it is customary to set up a numerical scale such as 1, 2, 3, 4, 5 or 1, 2, 3, ..., 9, 10 for characterizing human responses and subjective judgments with no valid justification except that these scales are easy to understand and convenient to use. But these numerical scales are arbitrary simplifications of the complex human mind; the human mind is not restricted to such simple numerical variations. In fact, human responses and subjective judgments are psychophysical phenomena that are fuzzy entities and therefore difficult to handle by conventional mathematics and probability theory. The fuzzy mathematical approach provides a more realistic insight into understanding and quantifying human responses. This paper presents a method for quantifying human responses and subjective judgments without assuming a pattern of linear or numerical variation for human responses. In particular, quantification and evaluation of linguistic judgments was investigated.

Steinlage, R. C.↗

Spatial and Temporal Dust Source Variability in Northern China Identified Using Advanced Remote Sensing Analysis

The aim of this research is to provide a detailed characterization of spatial patterns and temporal trends in the regional and local dust source areas within the desert of the Alashan Prefecture (Inner Mongolia, China). This problem was approached through multi-scale remote sensing analysis of vegetation changes. The primary requirements for this regional analysis are high spatial and spectral resolution data, accurate spectral calibration and good temporal resolution with a suitable temporal baseline. Landsat analysis and field validation along with the low spatial resolution classifications from MODIS and AVHRR are combined to provide a reliable characterization of the different potential dust-producing sources. The representation of intra-annual and inter-annual Normalized Difference Vegetation Index (NDVI) trend to assess land cover discrimination for mapping potential dust source using MODIS and AVHRR at larger scale is enhanced by Landsat Spectral Mixing Analysis (SMA). The combined methodology is to determine the extent to which Landsat can distinguish important soils types in order to better understand how soil reflectance behaves at seasonal and inter-annual timescales. As a final result mapping soil surface properties using SMA is representative of responses of different land and soil cover previously identified by NDVI trend. The results could be used in dust emission models even if they are not reflecting aggregate formation, soil stability or particle coatings showing to be critical for accurately represent dust source over different regional and local emitting areas.

spectral mixing analysis↗

The Dependence of Cirrus Gamma Size Distributions Expressed as Volumes in N(sub 0)-Lambda-Mu Phase Space and Bulk Cloud Properties on Environmental Conditions: Results from the Small Ice Particles in Cirrus Experiment (SPARTICUS)

The variability of cirrus ice microphysical properties is investigated using observations obtained during the Small Particles in Cirrus (SPARTICUS) campaign. An existing approach that represents a size distribution (SD) as a single gamma function using an ellipsoid of equally realizable solutions in (N(sub 0), lambda, mu) phase space is modified to automatically identify multiple modes in SDs and characterize each mode by such an ellipsoid. The modified approach is applied to ice crystals with maximum dimension D greater than15 micrometers collected by the 2-D stereo and 2-D precipitation probes on the Stratton Park Engineering Company Learjet. The dependencies of N(sub 0), mu, and lambda from each mode, total number concentration, bulk extinction, ice water content (IWC), and mass median maximum dimension D(sub mm) as a function of temperature T and cirrus type are then analyzed. The changes in the observed codependencies between N(sub 0), mu, and lambda, bulk extinction, IWC, and D(sub mm) with environmental conditions indicate that particles were larger at higher T during SPARTICUS. At most two modes were observed in any SD during SPARTICUS, with the average boundary between them at 115 micrometers, similar to past studies not using probes with shatter mitigating tips and artifact removal algorithms. The bimodality of the SDs increased with T. This and the differences in N(sub 0), mu, and lambda between the modes suggest that particles with smaller D nucleated more recently than particles with larger D, which grew via vapor deposition and aggregation. Because smaller crystals, whose concentrations are uncertain, make marginal contributions to higher order moments, the use of higher moments for evaluating model fields is suggested.

Ellipsoids↗

Projections of Temperature-Attributable Premature Deaths in 209 U.S. Cities Using a Cluster-Based Poisson Approach

Background: A warming climate will affect future temperature-attributable premature deaths. This analysis is the first to project these deaths at a near national scale for the United States using city and month-specific temperature-mortality relationships. Methods: We used Poisson regressions to model temperature-attributable premature mortality as a function of daily average temperature in 209 U.S. cities by month. We used climate data to group cities into clusters and applied an Empirical Bayes adjustment to improve model stability and calculate cluster-based month-specific temperature-mortality functions. Using data from two climate models, we calculated future daily average temperatures in each city under Representative Concentration Pathway 6.0. Holding population constant at 2010 levels, we combined the temperature data and cluster-based temperature-mortality functions to project city-specific temperature-attributable premature deaths for multiple future years which correspond to a single reporting year. Results within the reporting periods are then averaged to account for potential climate variability and reported as a change from a 1990 baseline in the future reporting years of 2030, 2050 and 2100. Results: We found temperature-mortality relationships that vary by location and time of year. In general, the largest mortality response during hotter months (April - September) was in July in cities with cooler average conditions. The largest mortality response during colder months (October-March) was at the beginning (October) and end (March) of the period. Using data from two global climate models, we projected a net increase in premature deaths, aggregated across all 209 cities, in all future periods compared to 1990. However, the magnitude and sign of the change varied by cluster and city. Conclusions: We found increasing future premature deaths across the 209 modeled U.S. cities using two climate model projections, based on constant temperature-mortality relationships from 1997 to 2006 without any future adaptation. However, results varied by location, with some locations showing net reductions in premature temperature-attributable deaths with climate change.

Schwartz, Joel D.↗

An Optimization Approach to Support Science Decision Making for Lunar Surface Exploration

Introduction: Scientific exploration is one of the three pillars of NASA’s Moon2Mars architecture, with crew surface extra vehicular activities (EVA) serving a critical enabling function. Development of surface EVA operational planning and execution, specifically integrating science and flight control teams (FCT), is currently being explored through analog scenarios. This integration, exercised, for example, through the Joint EVA and Hu-man Surface Mobility Test Team (JETT), allows for science input on EVA activities in near real-time through a Science Evaluation Room (SER), or Arte-mis science backroom, which integrates with the broader FCT through the Science Officer. The SER works within the FCT to support dynamic EVA planning in response to changes in operational constraints as well as science opportunities and re-prioritization, increasing the mission science return and accelerating the accomplishment of the Moon2Mars science objectives. The SER works within the FCT to provide recommendations to traverse execution in near real-time. One challenge is the requirement to deliver SER inputs to the FCT on operationally relevant timelines. Failure to do so may result in suboptimal execution of science exploration EVAs or even loss of key science objectives. To close this gap, we present a network optimization tool to allow the SER to provide rapid input to the FCT in response to changes in operational constraints or science opportunities. Inputs are predicated on approved science objectives, and clear rationale must be provided to the FCT for any requested change. Accordingly, this tool incorporates the Science Traceability Matrix (STM), SER prioritization scheme, and station characterization and action planning with operational constraints such as duration, traverse speed, and distance to maximize science objectives based on SER priorities, consistent with FCT operational requirements. Method: As a proof of concept, we used an existing linear programing software package used to simulate optimal routes through cellular metabolism. We built a Demonstrative Model with three STM objectives and four stations on a region of the Moon. The objectives were given an arbitrary prioritization and mapped to the stations through four possible crew actions. (Figs. 1 and 2). This station to STM mapping is consistent with the method used by the JETT5 Science Team to develop analog surface EVA science planning. We used a grid system with the landing site at the origin and the four stations placed across the positive x,y quadrant. Actions were assigned to each station and the accomplishment of those actions resulted in a numerical “reward” based on the ability of that action to achieve science objectives. The aggregate reward from each individual STM objective contributes to a global score (Science Yield), weighted by its priority. Operational constraints included a requirement to start and end at the landing site, 5 minutes each for initial station characterization and “clean up,” and variable total EVA time, traverse rate (fixed to 0.5 meters per second in our example), and time to perform each action (10, 5, 7, and 15 min for actions 1, 2, 3, and 4, respectively). Additional constraints and variables will be added in the future (e.g., sample mass, number of stations, traverse route constraints, illumination). Optimization. We converted the connections (arcs) between these stations (nodes) into a mixed integer linear programming optimization problem (arcs = constraints, nodes = variables) with the objective to maximize Science Yield. For any action, the Science Yield is equal to the relevance of that action to an STM objective [3, 2, and 1 point(s) for High, Med., and Low relevance, respectively], multiplied by the STM Objective Priority [3, 2, and 1 point(s) for High, Med., and Low priority, respectively]. This resulted in a model that computes the optimal station and action combination to maximize the Science Yield. These weightings can be adjusted by the SER as desired. Results: We explored three test cases for the Demonstrative Model. First, we set the maximum EVA duration to 120 minutes and computed the optimal route (Fig. 3A). The model suggested per-forming Actions 1 and 2 at Station P01, followed by Actions 1 and 2 at Station P02, and finally Actions 1 and 3 at Station P04 before returning to the Landing Site. Second, we adjusted the STM Objective Priori-ty order and computed the new optimal route (Fig. 3B). Under this situation, the model suggested per-forming all Actions at Station P02 followed by all Actions at Station P03. The previous test cases were relevant to SER planning activities. Next, we explored providing mid-EVA replanning input to the FCT. Scenario: While executing the Route in Fig. 3A the crew finishes at Station P01 and FCT decides that the EVA needs to finish in 45 minutes back at the Landing Site. FCT asks SER to recommend changes to the plan to accommodate this operation-al change. Using the model and incorporating these new constraints (start at Station P01, max. time of 45 min), the model suggested performing Actions 2 and 4 at Station P03 (Fig. 4), requiring 41 minutes to complete and return to the Landing Site. Interestingly, Station 3 was not part of the original route. Using the model, we determined the EVA would need 66 minutes, instead of 45, in order for the original Station P04 to yield a larger Science Yield than Station P03. The parametrization and simulation was per-formed in less than a minute, demonstrating the operational relevance of the approach. Future Efforts: The results from the Demonstrative Model suggest this tool can accelerate SER decision making on operationally relevant timelines. Use in analog activities, such as JETT5 or follow-ons, which have over a dozen stations for a crew to explore and over a dozen actions per station, will provide needed validation of the utility of this tool for planning EVAs, replanning mid-EVA, or planning follow-on EVAs based on previous results. Further integration with FCT execution monitoring tools may provide additional efficiency gains, al-lowing rapid and iterative exploration of operation-al and science decision space by the FCT and SER.

Science Operations↗

Above-Aircraft Cirrus Cloud and Aerosol Optical Depth From Hyperspectral Irradiances Measured By a Total-Diffuse Radiometer

This study develops the use of spectral total and diffuse irradiance measurements, made from a prototype hyperspectral total-diffuse sunshine pyranometer (SPN-S), to retrieve layer fine-mode aerosol (τ aer ) and total optical depths from airborne platforms. Additionally, we use spectral analysis in an attempt to partition the total optical depth into its τ aer and cirrus cloud optical depth (τ cld ) components in the absence of coarse-mode aerosols. Two retrieval methods are developed: one leveraging information in the diffuse irradiance and the other using spectral characteristics of the transmitted direct beam, with each approach best suited for specific cloud and aerosol conditions. The SPN-S has advantages over traditional sun photometer systems, including no moving parts and a low cost. However, a significant drawback of the instrument is that it is unable to measure the direct-beam irradiance as accurately as sun photometers. To compensate for the greater measurement uncertainty in the radiometric irradiances, these retrieval techniques employ ratioed inputs or spectral information to reduce output uncertainty. This analysis uses irradiance measurements from the SPN-S and the solar spectral flux radiometer (SSFR) aboard the National Aeronautics and Space Administration’s (NASA) P-3 aircraft during the 2018 deployment of the ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) campaign and the 2019 Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP 2 Ex) mission to quantify above-aircraft cirrus τ cld and derive vertical profiles of layer τ aer . Validation of the τ aer retrieval is accomplished by comparison with co-located measurements of direct solar irradiance made by the Sky-Scanning Sun-Tracking Atmospheric Research (4STAR) and in situ measurements of aerosol optical depth. For the aggregated 2018 ORACLES results, regression between the SPN-S-based method and sun photometer τ aer values yields a slope of 0.96 with an R 2 of 0.96, while the root mean square error (RMSE) is 3.0 x 10 -2 . When comparing the retrieved τ aer to profiles of integrated in situ measurements of optical extinction, the slope, R 2 , and RMSE values for ORACLES are 0.90, 0.96, and 3.4 x 10 -2 , and for CAMP 2 Ex they are 0.94, 0.97, and 3.4 x 10 -2 , respectively. This paper is a demonstration of methods for deriving cloud and aerosol optical properties in environments where both atmospheric constituents may be present. With improvements to the low-cost SPN-S radiometer instrument, it may be possible to extend these methods to a broader set of sampling applications, such as ground-based settings.

spectral instrument (SPN-S)↗

Cloud Microphysics in Hurricane Outflows: Observations in 'Bonnie' (1998) at 12 km Altitude

The water balance of a hurricane is controlled by boundary layer inflow, near vertical motion in the eyewall causing coalescence precipitation at above and residual ice precipitation at below freezing temperatures, and cirrus outflow at below -40 C aloft. In this paper we address the question of efficiency of water removal by this cirrus outflow which is important for the release of latent heat at high altitudes and its role in the dynamic flow at that level. During NASA's 1998 Convection and Moisture Experiment campaign we acquired microphysical outflow data in order to (1) determine the release and redistribution of latent heat near the top of hurricanes, (2) aid in TRMM algorithm development for remote sensing of precipitation, and (3) determine the optical/radiative characteristics of hurricane outflow. The data were acquired with Particle Measuring Systems two dimensional imaging spectrometers. On 23 August and again during the hurricane's landfall on 26 August, 1998, the NASA DC-8 aircraft penetrated hurricane 'Bonnie' four times each near 200 hPa pressure altitude. The eye crossing times were determined by (1) zero counts of cloud particles, (2) approximately 5 C increases in static and potential temperatures, and (3) minima in speeds and changes of direction of horizontal winds. The vertical winds showed shear between -6 m per second and +4 m per second and tangential winds approached 30 m per second in the eyewall. The particle volumes in the eyewall (determined by the pixels the particles shadowed in the direction of flight [x-direction] and normally to it by the number of diodes that they shadowed [y-direction]) ranged between 0.5 and 5.0 cubic centimeters per cubic meter. With a particle density near 0.2 g per cubic centimeter (determined from in situ melting and evaporation on a surface collector), the 1.0 g per meter corresponding mass of cloud ice ranged between 0.27 and 2.7 g per kilograms yielding horizontal fluxes between 8.1 and 81 g per square meters per second. The outflow ice was concentrated in crystals of a modal size of 190 micrometers. The particle size distributions were heavily skewed toward sizes with 98% of all cirrus particles smaller than the modal size comprising, however, only 20% of the mass. Thus the smaller than modal size particles dominantly affected the optical/radiative characteristics of the cloud, whereas the larger than modal size crystals determined the ice mass, hence dominated latent heating. Questions to be addressed relate to the origin of individual ice particles as the hurricane evolved and the likelihood of pristine and aggregate particle formation under the complicated conditions of rotation and outflow in the eyewall.

Pueschel, Rudolf F.↗