Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-Spherical”

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 127 records · Page 7

Dust in a Type Ia Supernova Progenitor: Spitzer Spectroscopy of Kepler's Supernova Remnant

Characterization of the relatively poorly-understood progenitor systems of Type Ia supernovae is of great importance in astrophysics, particularly given the important cosmological role that these supernovae play. Kepler's Supernova Remnant, the result of a Type Ia supernova, shows evidence for an interaction with a dense circumstellar medium (CSM), suggesting a single-degenerate progenitor system. We present 7.5-38 micron IR spectra of the remnant, obtained with the Spitzer Space Telescope, dominated by emission from warm dust. Broad spectral features at 10 and 18 micron, consistent with various silicate particles, are seen throughout. These silicates were likely formed in the stellar outflow from the progenitor system during the AGB stage of evolution, and imply an oxygen-rich chemistry. In addition to silicate dust, a second component, possibly carbonaceous dust, is necessary to account for the short-wavelength IRS and IRAC data. This could imply a mixed chemistry in the atmosphere of the progenitor system. However, non-spherical metallic iron inclusions within silicate grains provide an alternative solution. Models of collisionally-heated dust emission from fast shocks (> 1000 km/s) propagating into the CSM can reproduce the majority of the emission associated with non-radiative filaments, where dust temperatures are approx 80-100 K, but fail to account for the highest temperatures detected, in excess of 150 K. We find that slower shocks (a few hundred km/s) into moderate density material (n(sub o) approx 50-100 / cubic cm) are the only viable source of heating for this hottest dust. We confirm the finding of an overall density gradient, with densities in the north being an order of magnitude greater than those in the south.

Williams, Brian J.↗

Detection Thresholds of Falling Snow from Satellite-Borne Active and Passive Sensors

Precipitation, including rain and snow, is a critical part of the Earth's energy and hydrology cycles. Precipitation impacts latent heating profiles locally while global circulation patterns distribute precipitation and energy from the equator to the poles. For the hydrological cycle, falling snow is a primary contributor in northern latitudes during the winter seasons. Falling snow is the source of snow pack accumulations that provide fresh water resources for many communities in the world. Furthermore, falling snow impacts society by causing transportation disruptions during severe snow events. In order to collect information on the complete global precipitation cycle, both liquid and frozen precipitation must be collected. The challenges of estimating falling snow from space still exist though progress is being made. These challenges include weak falling snow signatures with respect to background (surface, water vapor) signatures for passive sensors over land surfaces, unknowns about the spherical and non-spherical shapes of the snowflakes, their particle size distributions (PSDs) and how the assumptions about the unknowns impact observed brightness temperatures or radar reflectivities, differences in near surface snowfall and total column snow amounts, and limited ground truth to validate against. While these challenges remain, knowledge of their impact on expected retrieval results is an important key for understanding falling snow retrieval estimations. Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types. For example, can a lake effect snow system with low (approx 2.5 km) cloud tops having an ice water content (IWC) at the surface of 0.25 g / cubic m and dendrite snowflakes be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. Here, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases since simulations provide all the information to determine the measurements from space and the ground truth. Sensitivity analyses were performed to better ascertain the relationships between multifrequency microwave and millimeter-wave sensor observations and the falling snow/underlying field of view. In addition, thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types were studied. Results will be presented for active radar at Ku, Ka, and W-band and for passive radiometer channels from 10 to 183 GHz.

Gail Skofronick-Jackson↗

Steps Toward an EOS-Era Aerosol Type Climatology

We still have a way to go to develop a global climatology of aerosol type from the EOS-era satellite data record that currently spans more than 12 years of observations. We have demonstrated the ability to retrieve aerosol type regionally, providing a classification based on the combined constraints on particle size, shape, and single-scattering albedo (SSA) from the MISR instrument. Under good but not necessarily ideal conditions, the MISR data can distinguish three-to-five size bins, two-to-four bins in SSA, and spherical vs. non-spherical particles. However, retrieval sensitivity varies enormously with scene conditions. So, for example, there is less information about aerosol type when the mid-visible aerosol optical depth (AOD) is less that about 0.15 or 0.2, or when the range of scattering angles observed is reduced by solar geometry, even though the quality of the AOD retrieval itself is much less sensitive to these factors. This presentation will review a series of studies aimed at assessing the capabilities, as well as the limitations, of MISR aerosol type retrievals involving wildfire smoke, desert dust, volcanic ash, and urban pollution, in specific cases where suborbital validation data are available. A synthesis of results, planned upgrades to the MISR Standard aerosol algorithm to improve aerosol type retrievals, and steps toward the development of an aerosol type quality flag for the Standard product, will also be covered.

Kahn, Ralph A.↗

Eyjafjallajokull Volcano Plume Particle-Type Characterization from Space-Based Multi-angle Imaging

The Multi-angle Imaging SpectroRadiometer (MISR) Research Aerosol algorithm makes it possible to study individual aerosol plumes in considerable detail. From the MISR data for two optically thick, near-source plumes from the spring 2010 eruption of the Eyjafjallaj kull volcano, we map aerosol optical depth (AOD) gradients and changing aerosol particle types with this algorithm; several days downwind, we identify the occurrence of volcanic ash particles and retrieve AOD, demonstrating the extent and the limits of ash detection and mapping capability with the multi-angle, multi-spectral imaging data. Retrieved volcanic plume AOD and particle microphysical properties are distinct from background values near-source, as well as for overwater cases several days downwind. The results also provide some indication that as they evolve, plume particles brighten, and average particle size decreases. Such detailed mapping offers context for suborbital plume observations having much more limited sampling. The MISR Standard aerosol product identified similar trends in plume properties as the Research algorithm, though with much smaller differences compared to background, and it does not resolve plume structure. Better optical analogs of non-spherical volcanic ash, and coincident suborbital data to validate the satellite retrieval results, are the factors most important for further advancing the remote sensing of volcanic ash plumes from space.

Kahn, Ralph A.↗

Tropical Atlantic Dust and Smoke Aerosol Variabilities Related to the Madden-Julian Oscillation in MODIS and MISR Observations

In this study, MODIS fine mode fraction and MISR non-spherical fraction are 2used to derive dust and smoke AOT components (tau(sub dust) and tau(sub smoke)) over the tropical Atlantic, and their variabilities related to the Madden-Julian Oscillation (MJO) are then investigated. Both MODIS and MISR show a very similar dust and smoke winter climatology. tau(sub dust) is found to be the dominant aerosol component over the tropical Atlantic while tau(sub smoke) is significantly smaller than tau(sub dust). The daily MODIS and MISR tau(sub dust) are overall highly correlated, with the correlation coefficients typically about 0.7 over the North Atlantic. The consistency between the MODIS and MISR dust and smoke aerosol climatology and daily variations give us confidence to use these two data sets to investigate their relative contributions to the total AOT variation associated with the MJO. However, unlike the MISR dust discrimination, which is based on particle shape retrievals, the smoke discrimination is less certain, based on assumed partitioning of maritime aerosol for both MISR and MODIS. The temporal evolution and spatial patterns of the tau(sub dust) anomalies associated with the MJO are consistent between MODIS and MISR. The tau(sub dust) anomalies are very similar to those of tau anomalies, and are of comparable magnitude. In contrast, the MJO-related tau(sub smoke) anomalies are rather small, and the tau(sub mar) anomalies are negligible. The consistency between the MODIS and MISR results suggests that dust aerosol is the dominant component on the intra-seasonal time scale over the tropical Atlantic Ocean.

Guo, Yanjuan↗

Steps Toward an EOS-Era Aerosol Air Mass Type Climatology

We still have a way to go to develop a global climatology of aerosol type from the EOS-era satellite data record that currently spans more than 12 years of observations. We have demonstrated the ability to retrieve aerosol type regionally, providing a classification based on the combined constraints on particle size, shape, and single-scattering albedo (SSA) from the MISR instrument. Under good but not necessarily ideal conditions, the MISR data can distinguish three-to-five size bins, two-to-four bins in SSA, and spherical vs. non-spherical particles. However, retrieval sensitivity varies enormously with scene conditions. So, for example, there is less information about aerosol type when the mid-visible aerosol optical depth (AOD) is less that about 0.15 or 0.2.

Kahn, Ralph A.↗

Evaluation of Model Microphysics Within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is ~0.25 meters per second too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were ~0.25 meters per second too slow, while the SBU-YLIN was 0.25 to 0.5 meters per second too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗

Evaluation of Model Microphysics within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is approx 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were approx 0.25 m/s too slow, while the SBU-YLIN was 0.25 to 0.5 m/s too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗

Absorption Efficiencies of Forsterite. I: DDA Explorations in Grain Shape and Size

We compute the absorption efficiency (Q(sub abs)) of forsterite using the discrete dipole approximation (DDA) in order to identify and describe what characteristics of crystal grain shape and size are important to the shape, peak location, and relative strength of spectral features in the 8 - 40 micron wavelength range. Using the DDSCAT code, we compute Q(sub abs) for non-spherical polyhedral grain shapes with a(sub eff) = 0.1 micron. The shape characteristics identified are: 1) elongation/reduction along one of three crystallographic axes; 2) asymmetry, such that all three crystallographic axes are of different lengths; and 3) the presence of crystalline faces that are not parallel to a specific crystallographic axis, e.g., non-rectangular prisms and (di)pyramids. Elongation/reduction dominates the locations and shapes of spectral features near 10, 11, 16, 23.5, 27, and 33.5 micron, while asymmetry and tips are secondary shape effects. Increasing grain sizes (0.1 - 1.0 micron) shifts the 10, 11 micron features systematically towards longer wavelengths and relative to the 11 micron feature increases the strengths and slightly broadens the longer wavelength features. Seven spectral shape classes are established for crystallographic a-, b-, and c-axes and include columnar and platelet shapes plus non-elongated or equant grain shapes. The spectral shape classes and the effects of grain size have practical application in identifying or excluding columnar, platelet or equant forsterite grain shapes in astrophysical environs. Identification of the shape characteristics of forsterite from 8 - 40 micron spectra provides a potential means to probe the temperatures at which forsterite formed.

circumstellar matter↗

Impact of Radiatively Interactive Dust Aerosols in the NASA GEOS-5 Climate Model: Sensitivity to Dust Particle Shape and Refractive Index

We investigate the radiative effects of dust aerosols in the NASA GEOS-5 atmospheric general circulation model. GEOS-5 is improved with the inclusion of a sectional aerosol and cloud microphysics module, the Community Aerosol and Radiation Model for Atmospheres (CARMA). Into CARMA we introduce treatment of the dust and sea salt aerosol lifecycle, including sources, transport evolution, and sinks. The aerosols are radiatively coupled to GEOS-5, and we perform a series of multi-decade AMIP-style simulations in which dust optical properties (spectral refractive index and particle shape distribution) are varied. Optical properties assuming spherical dust particles are from Mie theory, while those for non-spherical shape distributions are drawn from a recently available database for tri-axial ellipsoids. The climatologies of the various simulations generally compare well to data from the MODIS, MISR, and CALIOP space-based sensors, the ground-based AERONET, and surface measurements of dust deposition and concentration. Focusing on the summertime Saharan dust cycle we show significant variability in our simulations resulting from different choices of dust optical properties. Atmospheric heating due to dust enhances surface winds over important Saharan dust sources, and we find a positive feedback where increased dust absorption leads to increased dust emissions. We further find that increased dust absorption leads to a strengthening of the summertime Hadley cell circulation, increasing dust lofting to higher altitudes and strengthening the African Easterly Jet. This leads to a longer atmospheric residence time, higher altitude, and generally more northward transport of dust in simulations with the most absorbing dust optical properties. We find that particle shape, although important for radiance simulations, is a minor effect compared to choices of refractive index, although total atmospheric forcing is enhanced by greater than 10 percent for simulations incorporating a spheroidal shape distribution versus ellipsoidal or spherical shapes.

radiative forcing↗

On Spectral Invariance of Single Scattering Albedo for Water Droplets and Ice Crystals at Weakly Absorbing Wavelengths

The single scattering albedo omega(sub O lambda) in atmospheric radiative transfer is the ratio of the scattering coefficient to the extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, thus the single scattering albedo, are functions of wavelength lambda and droplet size r. This note shows that for water droplets at weakly absorbing wavelengths, the ratio omega(sub O lambda)(r)/omega(sub O lambda)(r (sub O)) of two single scattering albedo spectra is a linear function of omega(sub O lambda)(r). The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo spectrum omega(sub O lambda)(r) via one known spectrum omega(sub O lambda)(r (sub O)). We provide a simple physical explanation of the discovered relationship. Similar linear relationships were found for the single scattering albedo spectra of non-spherical ice crystals.

single scattering albedo↗

Aerosol Airmass Type Mapping Over the Urban Mexico City Region From Space-based Multi-angle Imaging

Using Multi-angle Imaging SpectroRadiometer (MISR) and sub-orbital measurements from the 2006 INTEX-B/MILAGRO field campaign, in this study we explore MISR's ability to map different aerosol air mass types over the Mexico City metropolitan area. The aerosol air mass distinctions are based on shape, size and single scattering albedo retrievals from the MISR Research Aerosol Retrieval algorithm. In this region, the research algorithm identifies dust-dominated aerosol mixtures based on non-spherical particle shape, whereas spherical biomass burning and urban pollution particles are distinguished by particle size. Two distinct aerosol air mass types based on retrieved particle microphysical properties, and four spatially distributed aerosol air masses, are identified in the MISR data on 6 March 2006. The aerosol air mass type identification results are supported by coincident, airborne high-spectral-resolution lidar (HSRL) measurements. Aerosol optical depth (AOD) gradients are also consistent between the MISR and sub-orbital measurements, but particles having single-scattering albedo of approx. 0.7 at 558 nm must be included in the retrieval algorithm to produce good absolute AOD comparisons over pollution-dominated aerosol air masses. The MISR standard V22 AOD product, at 17.6 km resolution, captures the observed AOD gradients qualitatively, but retrievals at this coarse spatial scale and with limited spherical absorbing particle options underestimate AOD and do not retrieve particle properties adequately over this complex urban region. However, we demonstrate how AOD and aerosol type mapping can be accomplished with MISR data over complex urban regions, provided the retrieval is performed at sufficiently high spatial resolution, and with a rich enough set of aerosol components and mixtures.

aerosol optical depth↗

A Numerical Testbed for Remote Sensing of Aerosols, and its Demonstration for Evaluating Retrieval Synergy from a Geostationary Satellite Constellation of GEO-CAPE and GOES-R

We present a numerical testbed for remote sensing of aerosols, together with a demonstration for evaluating retrieval synergy from a geostationary satellite constellation. The testbed combines inverse (optimal-estimation) software with a forward model containing linearized code for computing particle scattering (for both spherical and non-spherical particles), a kernel-based (land and ocean) surface bi-directional reflectance facility, and a linearized radiative transfer model for polarized radiance. Calculation of gas absorption spectra uses the HITRAN (HIgh-resolution TRANsmission molecular absorption) database of spectroscopic line parameters and other trace species cross-sections. The outputs of the testbed include not only the Stokes 4-vector elements and their sensitivities (Jacobians) with respect to the aerosol single scattering and physical parameters (such as size and shape parameters, refractive index, and plume height), but also DFS (Degree of Freedom for Signal) values for retrieval of these parameters. This testbed can be used as a tool to provide an objective assessment of aerosol information content that can be retrieved for any constellation of (planned or real) satellite sensors and for any combination of algorithm design factors (in terms of wavelengths, viewing angles, radiance and/or polarization to be measured or used). We summarize the components of the testbed, including the derivation and validation of analytical formulae for Jacobian calculations. Benchmark calculations from the forward model are documented. In the context of NASA's Decadal Survey Mission GEOCAPE (GEOstationary Coastal and Air Pollution Events), we demonstrate the use of the testbed to conduct a feasibility study of using polarization measurements in and around the O2 A band for the retrieval of aerosol height information from space, as well as an to assess potential improvement in the retrieval of aerosol fine and coarse mode aerosol optical depth (AOD) through the synergic use of two future geostationary satellites, GOES-R (Geostationary Operational Environmental Satellite R-series) and TEMPO (Tropospheric Emissions: Monitoring of Pollution). Strong synergy between GEOS-R and TEMPO are found especially in their characterization of surface bi-directional reflectance, and thereby, can potentially improve the AOD retrieval to the accuracy required by GEO-CAPE.

Aerosols↗

Evaluation of Mixed-Phase Microphysics Within Winter Storms Using Field Data and In Situ Observations

Snow prediction within models is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33-km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Cold-season Precipitation Experiment (GCPEx) experiment over southern Ontario, as well as a few years (12 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. The GCPex presentation will focus on verification using aircraft spirals through warm frontal snow band event on 18 February 2012. All the BMPs realistically simulated the structure of the band and the vertical distribution of snow/ice aloft, except the SBU-YLIN overpredicted slightly and Thompson (THOM) underpredicted somewhat. The Morrison (MORR) scheme produced the best slope size distribution for snow, while the Stony Brook (SBU) underpredicted and the THOM slightly overpredicted. Those schemes that have the slope intercept a function of temperature (SBU and WSM6) tended to perform better for that parameter than others, especially the fixed intercept in Goddard. Overall, the spread among BMPs was smaller than in other studies, likely because there was limited riming with the band. For the 15 cases at SBNY, which include moderate and heavy riming events, the non-spherical snow assumption (THOM and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM is approx. 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM, and MORR schemes were approx. 0.25 m/s too slow, while the SBU was 0.25 to 0.5 m/s too fast because of some excessive cloud water issues.

Colle, Brian A.↗

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗

9+ Years of CALIPSO PSC Observations: An Evolving Climatology

Polar stratospheric clouds (PSCs) play a crucial role in the springtime chemical depletion of ozone at high latitudes. PSC particles (primarily supercooled ternary solution, or STS droplets) provide sites for heterogeneous chemical reactions that transform stable chlorine and bromine reservoir species into highly reactive ozone-destructive forms. Furthermore, large nitric acid trihydrate (NAT) PSC particles can irreversibly redistribute odd nitrogen through gravitational sedimentation (a process commonly known as denitrification), which prolongs the ozone depletion process by slowing the reformation of the stable chlorine reservoirs. Spaceborne observations from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) lidar on the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) satellite are providing a rich new dataset for studying PSCs. CALIPSO is an excellent platform for studying polar processes with CALIOP acquiring, on average, over 300,000 backscatter profiles daily at latitudes between 55o and 82o in both hemispheres. PSCs are detected in the CALIOP backscatter profiles using a successive horizontal averaging scheme that enables detection of strongly scattering PSCs (e.g., ice) at the finest possible spatial resolution (5 km), while enhancing the detection of very tenuous PSCs (e.g., low number density NAT) at larger spatial scales (up to 135 km). CALIOP PSCs are separated into composition classes (STS; liquid/NAT mixtures; and ice) based on the ensemble 532-nm scattering ratio (the ratio of total-to-molecular backscatter) and 532-nm particulate depolarization ratio (which is sensitive to the presence of non-spherical, i.e. NAT and ice particles). In this paper, we will provide an overview of the CALIOP PSC detection and composition classification algorithm and then examine the vertical and spatial distribution of PSCs in the Arctic and Antarctic on vortex-wide scales for entire PSC seasons over the more than nine-year data record from 2006- 2015.

Pitts, Michael C.↗

CALIPSO Polar Stratospheric Cloud Observations from 2006-2015

Polar stratospheric clouds (PSCs) play a crucial role in the springtime chemical depletion of ozone at high latitudes. PSC particles (primarily supercooled ternary solution, or STS droplets) provide sites for heterogeneous chemical reactions that transform stable chlorine and bromine reservoir species into highly reactive ozone-destructive forms. Furthermore, large nitric acid trihydrate (NAT) PSC particles can irreversibly redistribute odd nitrogen through gravitational sedimentation (a process commonly known as denitrification), which prolongs the ozone depletion process by slowing the reformation of the stable chlorine reservoirs. Spaceborne observations from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) lidar on the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) satellite are providing a rich new dataset for studying PSCs. CALIPSO is an excellent platform for studying polar processes with CALIOP acquiring, on average, over 300,000 backscatter profiles daily at latitudes between 55o and 82o in both hemispheres. PSCs are detected in the CALIOP backscatter profiles using a successive horizontal averaging scheme that enables detection of strongly scattering PSCs (e.g., ice) at the finest possible spatial resolution (5 km), while enhancing the detection of very tenuous PSCs (e.g., low number density NAT) at larger spatial scales (up to 135 km). CALIOP PSCs are separated into composition classes (STS; liquid/NAT mixtures; and ice) based on the ensemble 532-nm scattering ratio (the ratio of total-to-molecular backscatter) and 532-nm particulate depolarization ratio (which is sensitive to the presence of non-spherical, i.e. NAT and ice particles). In this paper, we will provide an overview of the CALIOP PSC detection and composition classification algorithm and then examine the vertical and spatial distribution of PSCs in the Arctic and Antarctic on vortex-wide scales for entire PSC seasons over the more than nine-year data record from 2006- 2015.

Pitts, Michael C.↗

Orbit Stability of OSIRIS-REx in the Vicinity of Bennu Using a High-Fidelity Solar Radiation Model

Solar radiation pressure is one of the largest perturbing forces on the OSIRISRex trajectory as it orbits the asteroid Bennu. In this work, we investigate how forces due to solar radiation perturb the OSIRIS-REx trajectory in a high-fidelity model. The model accounts for Bennu's non-spherical gravity field, third-body gravity forces from the Sun and Jupiter, as well as solar radiation forces acting on a simplified spacecraft model. Such high-fidelity simulations indicate significant solar radiation pressure perturbations from the nominal orbit. Modifications to the initial design of the nominal orbit are found using a variation of parameters approach that reduce the perturbation in eccentricity by a factor of one-half.

Williams, Trevor W.↗