Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical cloud generation”

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 37 records · Page 2

On the cosmic ray diffusion in a violent interstellar medium

A variety of the available observational data on the cosmic ray (CR) spectrum, anisotropy and composition are in good agreement with a suggestion on the diffusion propagation of CR with energy below 10(15) eV in the interstellar medium. The magnitude of the CR diffusion coefficient and its energy dependence are determined by interstellar medium (ISM) magnetic field spectra. Direct observational data on magnetic field spectra are still absent. A theoretical model to the turbulence generation in the multiphase ISM is resented. The model is based on the multiple generation of secondary shocks and concomitant large-scale rarefactions due to supernova shock interactions with interstellar clouds. The distribution function for ISM shocks are derived to include supernova statistics, diffuse cloud distribution, and various shock wave propagation regimes. This permits calculation of the ISM magnetic field fluctuation spectrum and CR diffusion coefficient for the hot phase of ISM.

Bykov, A. M.↗

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗

A statistical data analysis and plotting program for cloud microphysics experiments

The analysis software developed for atmospheric cloud microphysics experiments conducted in the laboratory as well as aboard a KC-135 aircraft is described. A group of four programs was developed and implemented on a Hewlett Packard 1000 series F minicomputer running under HP's RTE-IVB operating system. The programs control and read data from a MEMODYNE Model 3765-8BV cassette recorder, format the data on the Hewlett Packard disk subsystem, and generate statistical data (mean, variance, standard deviation) and voltage and engineering unit plots on a user selected plotting device. The programs are written in HP FORTRAN IV and HP ASSEMBLY Language with the graphics software using the HP 1000 Graphics. The supported plotting devices are the HP 2647A graphics terminal, the HP 9872B four color pen plotter, and the HP 2608A matrix line printer.

Jordan, A. J.↗

Intercomparison of Sea-surface Temperature Anomaly Fields

Sea-surface temperature anomaly fields produced from analysis of HIRS2/MSU data at GLAS for the months November 1979, December 1981, March 1982, and July 1982 were compared to those produced from data from ships, AVHRR, and SMMR in a series of NASA sponsored sea-surface temperature workshops. The methods used by GLAS to analyze the HIRS2/MSU data were made to the algorithm used to retrieve sea and land surface temperatures and in the cloud filtering algorithm. One of the improvements involved generating weights for a given sounding, which decreased with increasing cloudiness and increasing humidity. These weights should be used in generating monthly mean fields. The workshop generated statistics for the weighted HIRS2/MSU field only for the month of July 1982.

Susskind, J.↗

Electrical Breakdown of Anodized Structures in a Low Earth Orbital Environmental

A comprehensive set of investigations involving arcing on a negatively biased anodized aluminum plate immersed in a low density argon plasma at low pressures (P(sub O), 7.5 x 10(exp -5) Torr) have been performed. These arcing experiments were designed to simulate electrical breakdown of anodized coatings in a Low Earth Orbital (LEO) environment. When electrical breakdown of an anodized layer occurs, an arc strikes, and there is a sudden flux of electrons accelerated into the ambient plasma. This event is directly followed by ejection of a quasi-neutral plasma cloud consisting of ejected material blown out of the anodized layer. Statistical analysis of plasma cloud expansion velocities have yielded a mean propagation velocity, v = (19.4 +/- 3.5) km/s. As the plasma cloud expands into the ambient plasma, energy in the form of electrical noise is generated. The radiated electromagnetic noise is detected by means of an insulated antenna immersed in the ambient plasma. The purpose of the investigations is (1) to observe and record the electromagnetic radiation spectrum resulting from the arcing process. (2) Make estimates of the travel time of the quasi-neutral plasma cloud based on fluctuations to several Langmuir probes mounted in the ambient plasma. (3) To study induced arcing between two anodized aluminum structures in close proximity.

Galofaro, J. T.↗

Development of Gridded Innovations and Observations Supplement to MERRA-2

Atmospheric reanalysis have become an important source of data for weather and climate research, owing to the continuity of the data, but especially because of the multitude of observational data included (radiosondes, commercial aircraft, retrieved data products and radiances). However, the presence of assimilated observations can vary based on numerous factors, and so it is difficult or impossible for a researcher to say with any degree of certainty how many and what type of observations contributed to the reanalysis data they are using at any give point in time or space, or their contribution to the eventual analyzed fields. For example, quality control, transmission interruptions, and station outages can occasionally affect data availability. While orbital paths are known, drift in certain instruments, cloud clearing, and the large number of available instruments makes it challenging to know which satellite is observing any region at any point in the diurnal cycle. Furthermore, there is information from the statistics generated by the data assimilation that can help understand the model and the quality of the reanalysis. Typically, the assimilated observations and their innovations are saved in observation-space data formats and are not easily made available to reanalysis users.Here, we present an early version of a data set has been developed to make the MERRA-2 assimilated observations available for rapid and general use, by simplifying the data format. The observations are binned to a grid similar as MERRA-2 and saved as netCDF. This data collection includes the mean and number of observations in the bin as well as its variance. The data will also include the innovations from the data assimilation, the forecast departure and the analysis increment, as well as bias correction (for satellite radiances). In this paper, we present the data format (called MERRA-2 Gridded Innovations and Observations or GIO) and its strengths and limitations with some initial testing and validation of the methodology.

Bosilovich, Michael G.↗

Monitoring Areal Snow Cover Using NASA Satellite Imagery

The objective of this project is to develop products and tools to assist in the hydrologic modeling process, including tools to help prepare inputs for hydrologic models and improved methods for the visualization of streamflow forecasts. In addition, this project will facilitate the use of NASA satellite imagery (primarily snow cover imagery) by other federal and state agencies with operational streamflow forecasting responsibilities. A GIS software toolkit for monitoring areal snow cover extent and producing streamflow forecasts is being developed. This toolkit will be packaged as multiple extensions for ArcGIS 9.x and an opensource GIS software package. The toolkit will provide users with a means for ingesting NASA EOS satellite imagery (snow cover analysis), preparing hydrologic model inputs, and visualizing streamflow forecasts. Primary products include a software tool for predicting the presence of snow under clouds in satellite images; a software tool for producing gridded temperature and precipitation forecasts; and a suite of tools for visualizing hydrologic model forecasting results. The toolkit will be an expert system designed for operational users that need to generate accurate streamflow forecasts in a timely manner. The Remote Sensing of Snow Cover Toolbar will ingest snow cover imagery from multiple sources, including the MODIS Operational Snowcover Data and convert them to gridded datasets that can be readily used. Statistical techniques will then be applied to the gridded snow cover data to predict the presence of snow under cloud cover. The toolbar has the ability to ingest both binary and fractional snow cover data. Binary mapping techniques use a set of thresholds to determine whether a pixel contains snow or no snow. Fractional mapping techniques provide information regarding the percentage of each pixel that is covered with snow. After the imagery has been ingested, physiographic data is attached to each cell in the snow cover image. This data can be obtained from a digital elevation model (DEM) for the area of interest.

Harshburger, Brian J.↗

IRIS Product Recommendations

This report presents the Applied Meteorology Unit's (AMU) evaluation of SIGMET Inc.'s Integrated Radar Information System (IRIS) Product Generator and recommendations for products emphasizing lightning and microburst tools. The IRIS Product Generator processes radar reflectivity data from the Weather Surveillance Radar, model 74C (WSR-74C), located on Patrick Air Force Base. The IRIS System was upgraded from version 6.12 to version 7.05 in late December 1999. A statistical analysis of atmospheric temperature variability over the Cape Canaveral Air Force Station (CCAFS) Weather Station provided guidance for the configuration of radar products that provide information on the mixed-phase (liquid and ice) region of clouds, between 0 C and -20 C. Mixed-phase processes at these temperatures are physically linked to electrification and the genesis of severe weather within convectively generated clouds. Day-to-day variations in the atmospheric temperature profile are of sufficient magnitude to warrant periodic reconfiguration of radar products intended for the interpretation of lightning and microburst potential of convectively generated clouds. The AMU also examined the radar volume-scan strategy to determine the scales of vertical gaps within the altitude range of the 0 C to -20 C isotherms over the Kennedy Space Center (KSC)/CCAFS area. This report present's two objective strategies for designing volume scans and proposes a modified scan strategy that reduces the average vertical gap by 37% as a means for improving radar observations of cloud characteristics in the critical 0 C to -20 C layer. The AMU recommends a total of 18 products, including 11 products that require use of the IRIS programming language and the IRIS User Product Insert feature. Included is a cell trends product and display, modeled after the WSR-88D cell trends display in use by the National Weather Service.

Short, David A.↗

The Earth Model Column Collaboratory (EMC2) v1.1: An Open-Source Ground-Based Lidar and Radar Instrument Simulator and Subcolumn Generator for Large-Scale Models

Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

Earth Model Column Collaboratory↗

Deriving Severe Hail Likelihood from Satellite Observations and Model Reanalysis Parameters using a Deep Neural Network

Geostationary satellite imagers, such as those of the Geostationary Operational Environmental Satellite (GOES) series, have been observing severe convection at 15–60-minute intervals for over 40 years. When properly assessed, such a data record can be valuable in efforts of estimating severe storm risk throughout the diurnal cycle based on automated detection of patterns consistently found atop severe storms. Furthermore, environmental conditions favorable for severe weather are well-known and are thought to be represented well by modern reanalysis products. Promoting resilience against such hazards on local and global scales is a chief goal the NASA Disasters program, which seeks to encourage use of satellite observations to mitigate risk. For instance, hail is the costliest severe weather hazard across the globe in terms of insured loss, but reporting inconsistencies for hail events globally make it difficult to develop models that can quantify the risk. Satellite observation and model reanalysis taken together have the potential to, with reasonable skill and specificity, characterize environmental conditions that are favorable for hazardous weather, and thereby enable creation of hazard climatologie. Such climatologies are particularly useful over regions without extensive radar networks or storm reporting. By mapping the multivariate combination of observed cloud features and reanalysis environmental parameters/indices to United States Next Generation Weather Radar (NEXRAD) radar-estimated Maximum Expected Size of Hail (MESH) by way of a deep neural network (DNN), estimates of likelihood for potentially severe hail can be produced. Such estimates are of greater complexity and efficiency than could be performed with previous multivariate or logistic regression analyses for observed points within convective systems. Statistical distributions of convective parameters from satellite and reanalysis are shown to highlight non-severe/severe class separation for well-known hailstorm predictors, e.g., overshooting cloud top characteristics, deep-layer wind shear, mid-level stability, helicity, and convective inhibition. These complex, multivariate predictor relationships are exploited within a DNN, which can efficiently produce a quantitative hail risk metric with better than 70% detection rate and under 30% false alarms. These hail classifications can then be aggregated across the satellite record to yield a hazard climatology for hail frequency and severity – knowledge of which is of particular interest to those who manage risk (e.g., insurers) and are seeking opportunities to identify hail-prone regions, particularly in developing nations. This NASA study uses satellite observations and model parameters in a DNN to perform climatological hailstorm analysis in support of catastrophe model development, with the hope of promoting risk resilience particularly in regions without adequate weather radar coverage.

Passive Remote Sensing↗

Toward The Development of Hailstorm Climatologies Derived From Reanalyses and Infared/Passive Microwave Satellite Imagers

Geostationary satellite imagers, such as those of the Geostationary Operational Environmental Satellite (GOES) and Meteosat series, provide both historical and near-real-time observations of cloud top patterns that are commonly associated with severe convection. Environmental conditions favorable for severe weather are thought to be represented well by reanalyses. Predicting exactly where convection and costly storm hazards like hail will occur using models or satellite imagery alone, however, is extremely challenging. The multivariate combination of satellite-observed cloud patterns with reanalysis environmental parameters, linked to United States Next Generation Weather Radar- (NEXRAD-) estimated Maximum Expected Size of Hail (MESH) using a deep neural network (DNN), enables estimation of potentially severe hail likelihood for any observed storm cell. These estimates are specifically designed to make hail likelihood distinctions based on satellite-indicated points of deep convection within environments favorable for storm development. We seek an approach that can be used to estimate climatological hailstorm frequency and risk throughout the historical satellite data record. This presentation demonstrates that statistical distributions of convective parameters from satellite and reanalysis show separation between non-severe/severe hailstorm classes for predictors including overshooting cloud top temperature and area characteristics, convective available potential energy, vertical wind shear, 500 hPa temperature, mid-level lapse rate, precipitable water, and convective inhibition. These complex, multivariate predictor relationships are exploited within a DNN to produce a hail likelihood metric with a critical success index of 0.504 and Heidke skill score of 0.403, which is exceptional among recent analogous hail studies. Furthermore, applications of the DNN to select case studies demonstrate good qualitative agreement between hail likelihood and MESH. These hail classifications are aggregated across an 11-year GOES-12/13 image database to derive a hail frequency and severity climatology, which denotes the Central Plains, the Midwest, and northwestern Mexico as being the most hail-prone regions within the domain studied. Opportunities for training and applying DNN-based hailstorm predictions to recently developed GOES-8/10/12/13/16 and Meteosat Second Generation convective storm detection and characterization climatologies over South America and South Africa, respectively, will also be presented.

Kristopher Bedka↗

Dynamics of the Impacts

The SL9 impacts are best known by their plumes, several of which were imaged towering 3000 km above Jupiter's limb. The heat released when they fen back on the atmosphere produced the famous infrared main events. The reentry shock was observed directly in CO emission and indirectly through thermally glowing dust. The relatively low color temperature of the latter implies that most of the energy released on reentry was radiated. If so, the R impact released some 3 x 10(exp 26) ergs, and the larger L, K, and G impacts as much as 3 x 10(exp 27) ergs. These estimates agree well with impact energies deduced from tidal disruption computations. The duration of the infrared events measures time aloft and hence provides a second measure of plume height. These indicate that plume height was only weakly dependent on impact energy. Evidently all the plumes were launched at roughly 10-13 km/s. Using a semi-analytic model for the deceleration, disintegration, and destruction of intruding bodies by an ever-vigilant atmosphere, we find that similar plume heights is a direct consequence of smaller impactors exploding at higher altitudes, in such a way that the different explosions were geometrically similar. We then compare the predictions of our model to the Venerian cratering record, which provides a good statistical test of impact deceleration by a thick atmosphere. Chemistry should have provided an independent measure of explosion altitude: abundant shock-generated CS, CS2 and HCN indicates a source in dry jovian air, above the putative water clouds. However, the Galileo Probe results seem to imply that we should expect no more. Observed water and S2 are consistent with a somewhat oxidized gas (presumably the comet itself), but the absence of SO2 and CO2 shows that conditions were neither too oxidizing nor the shocks too hot.

Zahnle, Kevin↗

Using High Frequency Passive Microwave, A-train, and TRMM Data to Evaluate Hydrometer Structure in the NASA GEOS-5 Data Assimilation System

Validating water vapor and prognostic condensate in global models remains a challenging research task. Model parameterizations are still subject to a large number of tunable parameters; furthermore, accurate and representative in situ observations are very sparse, and satellite observations historically have significant quantitative uncertainties. Progress on improving cloud / hydrometeor fields in models stands to benefit greatly from the growing inventory ofA-Train data sets. ill the present study we are using a variety of complementary satellite retrievals of hydrometeors to examine condensate produced by the emerging NASA Modem Era Retrospective Analysis for Research and Applications, MERRA, and its associated atmospheric general circulation model GEOS5. Cloud and precipitation are generated by both grid-scale prognostic equations and by the Relaxed Arakawa-Schubert (RAS) diagnostic convective parameterization. The high frequency channels (89 to 183.3 GHz) from AMSU-B and MRS on NOAA polar orbiting satellites are being used to evaluate the climatology and variability of precipitating ice from tropical convective anvils. Vertical hydrometeor structure from the Tropical Rainfall Measuring Mission (TRMM) and CloudSat radars are used to develop statistics on vertical hydrometeor structure in order to better interpret the extensive high frequency passive microwave climatology. Cloud liquid and ice water path data retrieved from the Moderate Resolution Imaging Spectroradiometer, MODIS, are used to investigate relationships between upper level cloudiness and tropical deep convective anvils. Together these data are used to evaluate cloud / ice water path, gross aspects of vertical hydrometeor structure, and the relationship between cloud extent and surface precipitation that the MERRA reanalysis must capture.

Robertson, Franklin↗

Applications and Improvement of a Coupled, Global and Cloud-Resolving Modeling System

Recently Grabowski (2001) and Khairoutdinov and Randall (2001) have proposed the use of 2D CFWs as a "super parameterization" [or multi-scale modeling framework (MMF)] to represent cloud processes within atmospheric general circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization used in conventional GCMs. A prototype Goddard MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM) is now being developed. The prototype includes the fvGCM run at 2.50 x 20 horizontal resolution with 32 vertical layers from the surface to 1 mb and the 2D (x-z) GCE using 64 horizontal and 32 vertical grid points with 4 km horizontal resolution and a cyclic lateral boundary. The time step for the 2D GCE would be 15 seconds, and the fvGCM-GCE coupling frequency would be 30 minutes (i.e. the fvGCM physical time step). We have successfully developed an fvGCM-GCE coupler for this prototype. Because the vertical coordinate of the fvGCM (a terrain-following floating Lagrangian coordinate) is different from that of the GCE (a z coordinate), vertical interpolations between the two coordinates are needed in the coupler. In interpolating fields from the GCE to fvGCM, we use an existing fvGCM finite- volume piecewise parabolic mapping (PPM) algorithm, which conserves the mass, momentum, and total energy. A new finite-volume PPM algorithm, which conserves the mass, momentum and moist static energy in the z coordinate, is being developed for interpolating fields from the fvGCM to the GCE. In the meeting, we will discuss the major differences between the two MMFs (i.e., the CSU MMF and the Goddard MMF). We will also present performance and critical issues related to the MMFs. In addition, we will present multi-dimensional cloud datasets (i.e., a cloud data library) generated by the Goddard MMF that will be provided to the global modeling community to help improve the representation and performance of moist processes in climate models and to improve our understanding of cloud processes globally (the software tools needed to produce cloud statistics and to identify various types of clouds and cloud systems from both high-resolution satellite and model data will be also presented).

Tao, W.-K.↗

Polarbrdf: A General Purpose Python Package for Visualization Quantitative Analysis of Multi-Angular Remote Sensing Measurements

The Bidirectional Reflectance Distribution Function (BRDF) is a fundamental concept for characterizing the reflectance property of a surface, and helps in the analysis of remote sensing data from satellite, airborne and surface platforms. Multi-angular remote sensing measurements are required for the development and evaluation of BRDF models for improved characterization of surface properties. However, multi-angular data and the associated BRDF models are typically multidimensional involving multi-angular and multi-wavelength information. Effective visualization of such complex multidimensional measurements for different wavelength combinations is presently somewhat lacking in the literature, and could serve as a potentially useful research and teaching tool in aiding both interpretation and analysis of BRDF measurements. This article describes a newly developed software package in Python (PolarBRDF) to help visualize and analyze multi-angular data in polar and False Color Composite (FCC) forms. PolarBRDF also includes functionalities for computing important multi-angular reflectance/albedo parameters including spectral albedo, principal plane reflectance and spectral reflectance slope. Application of PolarBRDF is demonstrated using various case studies obtained from airborne multi-angular remote sensing measurements using NASA's Cloud Absorption Radiometer (CAR). Our visualization program also provides functionalities for untangling complex surface/atmosphere features embedded in pixel-based remote sensing measurements, such as the FCC imagery generation of BRDF measurements of grasslands in the presence of wild fire smoke and clouds. Furthermore, PolarBRDF also provides quantitative information of the angular distribution of scattered surface/atmosphere radiation, in the form of relevant BRDF variables such as sunglint, hotspot and scattering statistics.

CAR↗

Comparison of Orbital Parameters for GEO Debris Predicted by LEGEND and Observed by MODEST: Can Sources of Orbital Debris be Identified?

Since 2002 the National Aeronautics and Space Administration (NASA) has carried out an optical survey of the debris environment in the geosynchronous Earth-orbit (GEO) region with the Michigan Orbital Debris Survey Telescope (MODEST) in Chile. The survey coverage has been similar for 4 of the 5 years allowing us to follow the orbital evolution of Correlated Targets (CTs), both controlled and un-controlled objects, and Un-Correlated Targets (UCTs). Under gravitational perturbations the distributions of uncontrolled objects, both CTs and UCTs, in GEO orbits will evolve in predictable patterns, particularly evident in the inclination and right ascension of the ascending node (RAAN) distributions. There are several clusters (others have used a "cloud" nomenclature) in observed distributions that show evolution from year to year in their inclination and ascending node elements. However, when MODEST is in survey mode (field-of-view approx.1.3deg) it provides only short 5-8 minute orbital arcs which can only be fit under the assumption of a circular orbit approximation (ACO) to determine the orbital parameters. These ACO elements are useful only in a statistical sense as dedicated observing runs would be required to obtain sufficient orbital coverage to determine a set of accurate orbital elements and then to follow their evolution. Identification of the source(s) for these "clusters of UCTs" would be advantageous to the overall definition of the GEO orbital debris environment. This paper will set out to determine if the ACO elements can be used to in a statistical sense to identify the source of the "clustering of UCTs" roughly centered on an inclination of 12deg and a RAAN of 345deg. The breakup of the Titan 3C-4 transtage on February 21, 1992 has been modeled using NASA s LEGEND (LEO-to-GEO Environment Debris) code to generate a GEO debris cloud. Breakup fragments are created based on the NASA Standard Breakup Model (including fragment size, area-to-mass (A/M), and delta-V distributions). Once fragments are created, they are propagated forward in time with a subroutine GEOPROP. Perturbations included in GEOPROP are those due to solar/lunar gravity, radiation pressure, and major geopotential terms. The question to be addressed: are the UCTs detected by MODEST in this inclination/RAAN region related to the Titan 3C-4 breakup? Discussion will include the observational biases in attempting to detect a specific, uncontrolled target during given observing session. These restrictions include: (1) the length of the observing session which is 8 hours or less at any given date or declination; (2) the assumption of ACO elements for detected object when the breakup model predicts debris with non-zero eccentricities; (3) the size and illumination or brightness of the debris predicted by the model and the telescope/sky limiting magnitude.

Barker, Edwin S.↗

Satellite Remote Sensing of Snow/Ice Albedo over the Himalayas

The Himalayan glaciers and snowpacks play an important role in the hydrological cycle over Asia. The seasonal snow melt from the Himalayan glaciers and snowpacks is one of the key elements to the livelihood of the downstream densely populated regions of South Asia. During the pre-monsoon season (April-May-June), South Asia not only experiences the reversal of the regional meridional tropospheric temperature gradient (i.e., the onset of the summer monsoon), but also is being bombarded by dry westerly airmass that transports mineral dust from various Southwest Asian desert and arid regions into the Indo-Gangetic Plains in northern India. Mixed with heavy anthropogenic pollution, mineral dust constitutes the bulk of regional aerosol loading and forms an extensive and vertically extended brown haze lapping against the southern slopes of the Himalayas. Episodic dust plumes are advected over the Himalayas, and are discernible in satellite imagery, resulting in dust-capped snow surface. Motivated by the potential implications of accelerated snowmelt, we examine the changes in radiative energetics induced by aerosol transport over the Himalayan snow cover by utilizing space borne observations. Our objective lies in the investigation of potential impacts of aerosol solar absorption on the Top-of-Atmosphere (TOA) spectral reflectivity and the broadband albedo, and hence the accelerated snowmelt, particularly in the western Himalayas. Lambertian Equivalent Reflectivity (LER) in the visible and near-infrared wavelengths, derived from Moderate Resolution Imaging Spectroradiometer radiances, is used to generate statistics for determining perturbation caused due to dust layer over snow surface in over ten years of continuous observations. Case studies indicate significant reduction of LER ranging from 5 to 8% in the 412-860nm spectra. Broadband flux observations, from the Clouds and the Earth's Radiant Energy System, are also used to investigate changes in shortwave TOA flux over dust-laden and dust-free snow covered regions. Additionally, spatio-temporal and intra-seasonal variations of LER, along with snow cover information, are used to characterize the seasonal melt pattern and thus to distinguish the outstanding aerosol-induced snowmelt signal. Results from this observational work are expected to provide better understanding of the radiative impact of aerosols over snow surface, especially its role in the Himalayan hydro-glacialogical variability.

Hsu, N. Christina↗