Search NASA⌕ Search

SEARCH · Search NASA

Results for “cloud computing; data analysis”

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 271 records · Page 15

A Quantum Annealing Computer Team Addresses Climate Change Predictability

The near confluence of the successful launch of the Orbiting Carbon Observatory2 on July 2, 2014 and the acceptance on August 20, 2015 by Google, NASA Ames Research Center and USRA of a 1152 qubit D-Wave 2X Quantum Annealing Computer (QAC), offered an exceptional opportunity to explore the potential of this technology to address the scientific prediction of global annual carbon uptake by land surface processes. At UMBC,we have collected and processed 20 months of global Level 2 light CO2 data as well as fluorescence data. In addition we have collected ARM data at 2sites in the US and Ameriflux data at more than 20 stations. J. Dorband has developed and implemented a multi-hidden layer Boltzmann Machine (BM) algorithm on the QAC. Employing the BM, we are calculating CO2 fluxes by training collocated OCO-2 level 2 CO2 data with ARM ground station tower data to infer to infer measured CO2 flux data. We generate CO2 fluxes with a regression analysis using these BM derived weights on the level 2 CO2 data for three Ameriflux sites distinct from the ARM stations. P. Gentine has negotiated for the access of K34 Ameriflux data in the Amazon and is applying a neural net to infer the CO2 fluxes. N. Talik validated the accuracy of the BM performance on the QAC against a restricted BM implementation on the IBM Softlayer Cloud with the Nvidia co-processors utilizing the same data sets. G. Nearing and K. Harrison have extended the GSFC LIS model with the NCAR Noah photosynthetic parameterization and have run a 10 year global prediction of the net ecosystem exchange. C. Pellisier is preparing a BM implementation of the Kalman filter data assimilation of CO2 fluxes. At UMBC, R. Prouty is conducting OSSE experiments with the LISNoah model on the IBM iDataPlex to simulate the impact of CO2 fluxes to improve the prediction of global annual carbon uptake. J. LeMoigne and D. Simpson have developed a neural net image registration system that will be used for MODIS ENVI and will be converted to a BM algorithm implementation on the QAC. The first integer adder has been implemented on the D-Wave 2X by A. Shehab that will perform HAAR wavelets for image compression of MODIS scenes. Finally, based on the next generations of QACs, we are preparing a 5-year performance road map on the scalability of the current QAC algorithms.

Science Data Processing↗

Forward and Inverse Models for Satellite Remote Sensors using Principal Component Analysis

Satellite remote sensors such as AIRS on Aqua, CrIS on S-NPP, NOAA20 and JPSS-2, IASI on Metop A, B, and C make millions of observations each day with thousands of spectral channels for each observation; this poses challenges for efficiently inversion of the inherently large dataset as needed to retrieve atmospheric and surface properties. This presentation will illustrate the use of Principal Component Analysis (PCA) to speed up radiative transfer forward model calculations and to stabilize the inversion algorithms. A Principal Component-based radiative transfer model (PCRTM) developed at NASA Langley Research Center can simulate top of atmosphere (TOA) radiance or reflectance spectra from 50 cm-1 to 50000 cm-1 (200 m to 0.20 m quickly and accurately. PCRTM demonstrated very high accuracy relative to reference line-by-line radiative transfer models and it saves orders of magnitude computational time. Examples of the PCRTM model developed for hyperspectral sensors such as AIRS, CrIS, IASI, NAST-I, SHIS, CPF, TEMPO, SBG, OMI, and SCIAMACHY will be presented. In addition to using the PCRTM as forward model, the NASA Langley developed inversion algorithm also uses PCA to compress the state vector into a compressed dimension to speed up and stabilize the inversion process. Examples of retrieved atmospheric temperature, water vapor, CO2, CO, CH4, N2O, and O3 profiles, cloud properties (optical depth, size, phase, and height), and surface properties (surface emissivity spectra and skin temperatures) will be presented. This algorithm is being transitioned to the NASA Sounder SIPS and NASA's Goddard Earth Sciences Data and Information Services Center (GES DISC).

forward model↗

A Redetermination of the Dissociation Energy of MgO(+)

In 1986, we reported a dissociation energy (D(sub 0) of 2.31 eV for the X(sup 2)Pi ground state of MgO(+). This value was determined by computing the dissociation energy to the Mg(2+) + O(-) limit and adjusting the value to the Mg(+) + O limit using the experimental Ionization Potential (IP) of Mg(+) and the Electron Affinity (EA) of O. The success of this method relies on the assumption that there is little covalent contribution to the bonding. The very small (0.04 eV) correlation contribution to the binding energy was taken as corroboration for the validity of this approach. Our earlier theoretical value was estimated to be accurate to at least 0.2 eV. It is in excellent agreement with the subsequent value of 2.30 +/- 0.13 eV determined by Freiser and co-workers from photodissociation experiments. It is also consistent with the upper (less than 3.1 eV) and lower (greater than 1.1 eV) bounds determined by Rowe obtained by studying the reactions of Mg(+) with 03 and NO2. However, it is inconsistent with an upper bound of 1.7 eV reported by Kappes and Staley based on their failure to observe MgO(+) in the reaction of Mg(+) with N2O. The picture became somewhat clouded, however, by the recent guided-ion beam mass spectrometric studies of Dalleska and Armentrout. Their initial analysis of the reaction data for Mg(+) + O2 lead to a bond dissociation energy of 2.92 +/- 0.25 eV, which is considerably larger than the value of 2.47 +/- 0.06 eV deduced from their studies of the Mg(+)+NO2 reaction.

Bauschlicher, Charles W., Jr.↗

The alpha Centauri Line of Sight: D/H Ratio, Physical Properties of Local Interstellar Gas, and Measurement of Heated Hydrogen (The 'Hydrogen Wall') Near the Heliopause

We analyze high-resolution spectra of the nearby (1.34 pc) stars alpha Cen A (G2 V) and alpha Cen B (K1 V), which were obtained with the Goddard High Resolution Spectrograph on the Hubble Space Telescope. The observations consist of echelle spectra of the Mg II 2800 A and Fe II 2599 A resonance lines and the Lyman-alpha lines of hydrogen and deuterium. The interstellar gas has a velocity (v = - 18.0 +/- 0.2 km/s) consistent with the local flow vector proposed for this line of sight by Lailement & Berlin (1992). The temperature and nonthermal velocity inferred from the Fe II, Mg II, and D I line profiles are T = 5400 +/- 500 K and xi = 1.20 +/- 0.25 km/s, respectively. However, single-component fits to the H I Lyman-alpha lines yield a Doppler parameter (b(sub HI) = 11.80 km/s) that implies a significantly warmer temperature of 8350 K, and the velocity of the H I absorption (v = - 15.8 +/- 0.2 km/s) is redshifted by about 2.2 km/s with respect to the Fe II, Mg II, and D I lines. The one-component model of the interstellar gas suggests natural logarithm N base HI = 18.03 +/- 0.01 and D/H = (5.7 +/- 0.2) x 10(exp -6) . These parameters lead to a good fit to the observed spectra, but this model does not explain the higher temperature and redshift of H I relative to the other interstellar lines. The most sensible way to resolve the discrepancy between H(I) and the other lines is to add a second absorption component to the H(I) lines. This second component is hotter (T approx. equals 30,000 K), is redshifted relative to the primary component by 2-4 km/s, and has a column density too low to be detected in the Fe(II), Mg(II), and D(I) lines. We propose that the gas responsible for this component is located near the heliopause, consisting of the heated H I gas from the interstellar medium that is compressed by the solar wind. This so-called 'hydrogen wall' is predicted by recent multifluid gasdynamical models of the interstellar gas and solar wind interaction. Our data provide the first measurements of the temperature and column density of H(I) in the hydrogen wall. After considering the effects that a corresponding hydrogen wall around alpha Cen would have on our analysis, our best estimates for the parameters of the solar hydrogen wall are natural log N(sup (2))(H(I)) = 14.74 +/- 0.24, b(sup (2))(H(I)) = 21.9 +/- 1.7 km/s (corresponding to T = 29,000 +/- 5000 K), and v(sup (2))(H(I)) greater than -16km/s. Unfortunately, the existence of this heated H(I) reduces our ability to compute the H(I) column density of the interstellar medium accurately because, with slight alterations to our assumed stellar Lyman-alpha profiles, we discovered that acceptable two-component fits also exist with natural log N(H(I))approx. 17.6. We, therefore, quote large error bars for the H I column density along the alpha Cen line of sight, natural log N(H(I)) = 17.80 +/- 0.30. For this range in N(H(I)), n(H(I)) = 0.15 /cu.cm (+/- a factor of 2) and D/H = (0.5-1.9) x 10(exp -5). This is the first direct measurement of the H(I) density in a local cloud and allows us to predict the distance from the Sun to the edge of the local cloud along various lines of sight. This range in D/H is consistent with the value D/H = 1.6 x 10(exp -5) previously derived for the Capella and Procyon lines of sight. We cannot tell whether D/H ratio varies or is constant in the local interstellar medium, but we do find that the D(I)/Mg(II) ratio for the alpha Cen line of sight is about 4 times smaller than for the Capella and Procyon lines of sight. Therefore, either D/H or the Mg depletion varies significantly over distance scales of only a few parsecs.

Linsky, Jeffrey L.↗

Expected Performance of Ozone Climate Data Records from Ozone Mapping and Profiler Suite Limb Profiler

The Ozone Mapping and Profiler Suite Limb Profiler (OMPS/LP) was launched on board of the Soumi NPP space platform in late October 2011. It provides ozone-profiling capability with high-vertical resolution from 60 Ian to cloud top. In this study, an end-to-end Observing System Simulation Experiment (OSSE) of OMPS/LP ozone is discussed. The OSSE was developed at NASA's Global Modeling and Assimilation Office (GMAO) using the Goddard Earth Observing System (GEOS-5) data assimilation system. The "truth" for this OSSE is built by assimilating MLS profiles and OMI ozone columns, which is known to produce realistic three-dimensional ozone fields in the stratosphere and upper troposphere. OMPS/LP radiances were computed at tangent points computed by an appropriate orbital model. The OMPS/LP forward RT model, Instrument Models (IMs) and EDR retrieval model were introduced and pseudo-observations derived. The resultant synthetic OMPS/LP observations were evaluated against the "truth" and subsequently these observations were assimilated into GEOS-5. Comparison of this assimilated dataset with the "truth" enables comparisons of the likely uncertainties in 3-D analyses of OMPS/LP data. This study demonstrated the assimilation capabilities of OMPS/LP ozone in GEOS-5, with the monthly, zonal mean (O-A) smaller than 0.02ppmv at all levels, the nns(O-A) close to O.lppmv from 100hPa to 0.2hPa; and the mean(O-B) around the 0.02ppmv for all levels. The monthly zonal mean analysis generally agrees to within 2% of the truth, with larger differences of 2-4% (0.1-0.2ppmv) around 10hPa close to North Pole and in the tropical tropopause region, where the difference is above 20% due to the very low ozone concentrations. These OSSEs demonstrated that, within a single data assimilation system and the assumption that assimilated MLS observations provide a true rendition of the stratosphere, the OMPS/LP ozone data are likely to produce accurate analyses through much of the stratosphere. Additional studies are needed to evaluate the importance of the extremely high vertical resolution of the OMPS/LP ozone data

Xu, P. Q.↗

Forecasting Lightning at Kennedy Space Center/Cape Canaveral Air Force Station, Florida

The Applied Meteorology Unit (AMU) developed a set of statistical forecast equations that provide a probability of lightning occurrence on Kennedy Space Center (KSC) I Cape Canaveral Air Force Station (CCAFS) for the day during the warm season (May September). The 45th Weather Squadron (45 WS) forecasters at CCAFS in Florida include a probability of lightning occurrence in their daily 24-hour and weekly planning forecasts, which are briefed at 1100 UTC (0700 EDT). This information is used for general scheduling of operations at CCAFS and KSC. Forecasters at the Spaceflight Meteorology Group also make thunderstorm forecasts for the KSC/CCAFS area during Shuttle flight operations. Much of the current lightning probability forecast at both groups is based on a subjective analysis of model and observational data. The objective tool currently available is the Neumann-Pfeffer Thunderstorm Index (NPTI, Neumann 1971), developed specifically for the KSCICCAFS area over 30 years ago. However, recent studies have shown that 1-day persistence provides a better forecast than the NPTI, indicating that the NPTI needed to be upgraded or replaced. Because they require a tool that provides a reliable estimate of the daily thunderstorm probability forecast, the 45 WS forecasters requested that the AMU develop a new lightning probability forecast tool using recent data and more sophisticated techniques now possible through more computing power than that available over 30 years ago. The equation development incorporated results from two research projects that investigated causes of lightning occurrence near KSCICCAFS and over the Florida peninsula. One proved that logistic regression outperformed the linear regression method used in NPTI, even when the same predictors were used. The other study found relationships between large scale flow regimes and spatial lightning distributions over Florida. Lightning, probabilities based on these flow regimes were used as candidate predictors in the equation development. Fifteen years (1 989-2003) of warm season data were used to develop the forecast equations. The data sources included a local network of cloud-to-ground lightning sensors called the Cloud-to-Ground Lightning Surveillance System (CGLSS), 1200 UTC Florida synoptic soundings, and the 1000 UTC CCAFS sounding. Data from CGLSS were used to determine lightning occurrence for each day. The 1200 UTC soundings were used to calculate the synoptic-scale flow regimes and the 1000 UTC soundings were used to calculate local stability parameters, which were used as candidate predictors of lightning occurrence. Five logistic regression forecast equations were created through careful selection and elimination of the candidate predictors. The resulting equations contain five to six predictors each. Results from four performance tests indicated that the equations showed an increase in skill over several standard forecasting methods, good reliability, an ability to distinguish between non-lightning and lightning days, and good accuracy measures and skill scores. Given the overall good performance the 45 WS requested that the equations be transitioned to operations and added to the current set of tools used to determine the daily lightning probability of occurrence.

Lambert, Winfred↗

Icing Frequencies Experienced During Climb and Descent by Fighter-Interceptor Aircraft

Data and analyses are presented on the relative frequencies of occurrence and severity of icing cloud layers encountered by jet aircraft in the climb and descent phases of flights to high altitudes. Fighter-interceptor aircraft operated by the Air Defense Command (USAF) at bases in the Duluth and Seattle areas collected the data with icing meters installed for a l-year period. The project was part of an extensive program conducted by the NACA to collect Icing cloud data for evaluating the icing problem relevant to routine operations. The average frequency of occurrence of icing was found to be about 5 percent of the number of climbs and descents during 1 year of operations The icing encounters were predominantly in the low and middle cloud layers, decreasing above 15,000 feet to practically none above 25,000 feet. The greatest thickness of ice that would accumulate on any aircraft component (as indicated by the accretion on a small object) was measured with the icing meters. The ice thicknesses on a small sensing probe averaged less than 1/32 inch and did not exceed 1/2 inch. Such accumulations are relatively small when compared with those that can form during horizontal flight in icing clouds. The light accretions resulted from relatively steep angles of flight through generally thin cloud layers. Because of the limited statistical reliability of the results, an analysis was made using previous statistics on icing clouds below an altitude of 20,000 feet to determine the general icing severity probabilities. The calculations were made using adiabatic lifting as a basis to establish the liquid-water content. Probabilities of over-all ice accretions on a small object as a function of airspeed and rate of climb were computed from the derived water contents. These results were then combined with the probability of occurrence of icing in order to give the icing severity that can be expected for routine aircraft operations.

Perkins, Porter J.↗

Airborne lidar observations of smoke haze during SCAR-B 1995

The Smoke, Clouds, Aerosol, and Radiation Brazil(SCAR-B) field campaign was conducted to study the effects that widespread and persistent biomass burning have upon radiative and chemical processes in the atmosphere. The radiative transfer characteristics of the atmosphere are altered by the introduction of particulate and gaseous materials which are the products of the combustion of vegetative material at ground level. These substances are transported and distributed horizontally and vertically by atmospheric dynamical processes which may be perturbed by the heat energy from the fires. As the pollutants disperse, their physical and chemical properties change substantially. A complete description of the effects of smoke requires that the evolution back to the natural situation be fully examined. A most important component of smoke haze investigation is finding its vertical and horizontal distribution in relation to the driving factors of dynamics and the related horizontal transport. In this presentation, we employ data from the Cloud Lidar System(CLS), carried aboard the NASA ER-2 aircraft, to provide a unique view of the particulate or aerosol loading produced by fires, especially with regard to the geometrical distribution of the aerosols in the vertical plane. The lidar has the ability to measure aerosol optical properties in a continuous fashion at quite fine vertical and horizontal resolution. The results from the lidar provide measurements that are largely independent of influences that corrupt passive instruments and thus it can serve to corroborate their results. The extended horizontal and vertical range of lidar results can also augment ground based and airborne in situ measurements which have limited horizontal and vertical scope. We present the results of our analysis of CLS observations taken during the SCAR-B field campaign. Observations of the the aerosol optical thickness from the Aerosol Robotic Network(AERONET) of solar photometers are employed in conjunction with CLS data to derive extinction to backscatter ratio values which are used to convert the lidar backscatter coefficient into extinction coefficient. The extinction coefficient is integrated vertically to find aerosol optical thickness along ER-2 flight tracks. We use images of the CLS derived extinction coefficient to depict its horizontal and vertical distribution. Multispectral photometer optical thickness is used to compute the Angstrom exponent. With these, we examine the hypothesis that the values of extinction to backscatter ratio can be related to the Anstrom coefficient since both of these would be a function of the refractive index and size distribution of the aerosols.

Hart, William D.↗

HEASARC - The High Energy Astrophysics Science Archive Research Center

The High Energy Astrophysics Science Archive Research Center (HEASARC) is NASA's archive for high-energy astrophysics and cosmic microwave background (CMB) data, supporting the broad science goals of NASA's Physics of the Cosmos theme. It provides vital scientific infrastructure to the community by standardizing science data formats and analysis programs, providing open access to NASA resources, and implementing powerful archive interfaces. Over the next five years the HEASARC will ingest observations from up to 12 operating missions, while serving data from these and over 30 archival missions to the community. The HEASARC archive presently contains over 37 TB of data, and will contain over 60 TB by the end of 2014. The HEASARC continues to secure major cost savings for NASA missions, providing a reusable mission-independent framework for reducing, analyzing, and archiving data. This approach was recognized in the NRC Portals to the Universe report (2007) as one of the HEASARC's great strengths. This poster describes the past and current activities of the HEASARC and our anticipated developments in coming years. These include preparations to support upcoming high energy missions (NuSTAR, Astro-H, GEMS) and ground-based and sub-orbital CMB experiments, as well as continued support of missions currently operating (Chandra, Fermi, RXTE, Suzaku, Swift, XMM-Newton and INTEGRAL). In 2012 the HEASARC (which now includes LAMBDA) will support the final nine-year WMAP data release. The HEASARC is also upgrading its archive querying and retrieval software with the new Xamin system in early release - and building on opportunities afforded by the growth of the Virtual Observatory and recent developments in virtual environments and cloud computing.

Smale, Alan P.↗

Leveraging Data Intensive Computing to Support Automated Event Services

A large portion of Earth Science investigations is phenomenon- or event-based, such as the studies of Rossby waves, mesoscale convective systems, and tropical cyclones. However, except for a few high-impact phenomena, e.g. tropical cyclones, comprehensive records are absent for the occurrences or events of these phenomena. Phenomenon-based studies therefore often focus on a few prominent cases while the lesser ones are overlooked. Without an automated means to gather the events, comprehensive investigation of a phenomenon is at least time-consuming if not impossible. An Earth Science event (ES event) is defined here as an episode of an Earth Science phenomenon. A cumulus cloud, a thunderstorm shower, a rogue wave, a tornado, an earthquake, a tsunami, a hurricane, or an EI Nino, is each an episode of a named ES phenomenon," and, from the small and insignificant to the large and potent, all are examples of ES events. An ES event has a finite duration and an associated geolocation as a function of time; its therefore an entity in four-dimensional . (4D) spatiotemporal space. The interests of Earth scientists typically rivet on Earth Science phenomena with potential to cause massive economic disruption or loss of life, but broader scientific curiosity also drives the study of phenomena that pose no immediate danger. We generally gain understanding of a given phenomenon by observing and studying individual events - usually beginning by identifying the occurrences of these events. Once representative events are identified or found, we must locate associated observed or simulated data prior to commencing analysis and concerted studies of the phenomenon. Knowledge concerning the phenomenon can accumulate only after analysis has started. However, except for a few high-impact phenomena. such as tropical cyclones and tornadoes, finding events and locating associated data currently may take a prohibitive amount of time and effort on the part of an individual investigator. And even for these high-impact phenomena, the availability of comprehensive records is still only a recent development. A major reason for the lack of comprehensive ,records for the majority of the ES phenomena is the perception that they do not pose immediate and/or severe threat to life and property and are thus not consistently tracked. monitored, and catalogued. Many phenomena even lack commonly accepted criteria for definitions. However. the lack of comprehensive records is also due to the increasingly prohibitive volume of observations and model data that must be examined. NASA Earth Observing System Data Information System (EOSDIS) alone archives several petabytes (PB) of satellite remote sensing data and steadily increases. All of these factors contribute to the difficulty of methodically identifying events corresponding to a given phenomenon and significantly impede systematic investigations. In the following we present a couple motivating scenarios, demonstrating the issues faced by Earth scientists studying ES phenomena.

Clune, Thomas L.↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). V. First Results on the Magnetic Field Orientation of Galaxies

We present the analysis of the magnetic field ( B -field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA of 14 nearby (<20 Mpc) galaxies with resolved (5"-18" ; 90 pc–1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 μm . We compute the magnetic pitch-angle (𝜓 B ) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter ( ζ ) to estimate the disordered-to-ordered ratio of spiral B -fields. We find FIR and radio wavelengths to not generally trace the same B -field morphology in galaxies. The 𝜓 B profiles tend to be more ordered across all galactocentric radii in radio ( ζ 6cm = 0.93 + 0.03) than in FIR ( ζ 154μm = 0.84 + 0.14). For spiral galaxies, FIR B -fields are 2%–75% more turbulent than the radio B -fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B -fields along the galactic outflows than radio polarization. Our results suggest that the B -fields associated with dense, dusty, turbulent star-forming regions (those traced at FIR) are less ordered than warmer, less dense regions (those traced at radio) of the interstellar medium. The FIR B -fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of a few hundred parsecs in the disk of spiral galaxies than the radio B -fields.

magnetic field (B-field)↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Model

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology from the small MVD conditions. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Extragalactic Magnetism with SOFIA (SALSA Legacy Program) - V. First Results on the Magnetic Field Orientation of Galaxies

We present the analysis of the magnetic field (B-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby (< 20 Mpc) galaxies with resolved (5′′ − 18′′; 90 pc–1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 μm. We compute the magnetic pitch angle (Ψ B ) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter (ζ) to estimate the disordered-to-ordered ratio of spiral B -fields. We find FIR and radio wavelengths to not generally trace the same B -field morphology in galaxies. The Ψ B profiles tend to be more ordered across all galactocentric radii in radio (ζ 6cm = 0.93 ± 0.03) than in FIR (ζ 154μm = 0.84 ± 0.14). For spiral galaxies, FIR B-fields are 2−75% more turbulent than the radio B -fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B -fields along the galactic outflows than radio polarization. Our results suggest that the B -fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR B -fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio B -fields.

SALSA↗

Inference of Ice Cloud Properties from High-spectral Resolution Infrared Observations

The theoretical basis is explored for inferring the microphysical properties of ice crystal from high-spectral resolution infrared observations. A radiative transfer model is employed to simulate spectral radiances to address relevant issues. The extinction and absorption efficiencies of individual ice crystals, assumed as hexagonal columns for large particles and droxtals for small particles, are computed from a combination of the finite- difference time-domain (FDTD) technique and a composite method. The corresponding phase functions are computed from a combination of FDTD and an improved geometric optics method (IGOM). Bulk scattering properties are derived by averaging the single- scattering properties of individual particles for 30 particle size distributions developed from in situ measurements and for additional four analytical Gamma size distributions for small particles. The non-sphericity of ice crystals is shown to have a significant impact on the radiative signatures in the infrared (IR) spectrum; the spherical particle approximation for inferring ice cloud properties may result in an overest&ation of the optical thickness and an inaccurate retrieval of effective particle size. Furthermore, we show that the error associated with the use of the Henyey-Greenstein phase function can be as larger as 1 K in terms of brightness temperature for larger particle effective size at some strong scattering wavenumbers. For small particles, the difference between the two phase functions is much less, with brightness temperatures generally differing by less than 0.4 K. The simulations undertaken in this study show that the slope of the IR brightness temperature spectrum between 790-960/cm is sensitive to the effective particle size. Furthermore, a strong sensitivity of IR brightness temperature to cloud optical thickness is noted within the l050-1250/cm region. Based on this spectral feature, a technique is presented for the simultaneous retrieval of the visible optical thickness and effective particle size from high spectral resolution infrared data under ice cloudy con&tion. The error analysis shows that the uncertainty of the retrieved optical thickness and effective particle size has a small range of variation. The error for retrieving particle size in conjunction with an uncertainty of 5 K in cloud'temperature, or a surface temperature uncertainty of 2.5 K, is less than 15%. The corresponding e m r in the uncertainty of optical thickness is within 5-2096, depending on the value of cloud optical thickness. The applicability of the technique is demonstrated using the aircraft-based High- resolution Interferometer Sounder (HIS) data from the Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) in 1996 and the First ISCCP Regional Experiment - Arctic Clouds Experiment (FIRE-ACE) in 1998.

Huang, Hung-Lung↗

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗

Radiative Susceptibility of Cloudy Atmospheres to Droplet Number Perturbations: 1. Theoretical Analysis and Examples from MODIS

Theoretical and satellite-based assessments of the sensitivity of broadband shortwave radiative fluxes in cloudy atmospheres to small perturbations in the cloud droplet number concentration (N) of liquid water clouds under constant water conditions are performed. Two approaches to study this sensitivity are adopted: absolute increases in N, for which the radiative response is referred to as absolute cloud susceptibility, and relative increases in N or relative cloud susceptibility. Estimating the former is more challenging as it requires an assumed value for either cloud liquid water content or geometrical thickness; both susceptibilities require an assumed relationship between the droplet volume and effective radius. Expanding upon previous susceptibility studies, present radiative calculations include the effect of AN perturbations on droplet asymmetry parameter and single-scattering albedo, in addition to extinction. Absolute cloud susceptibility has a strong nonlinear dependence on the droplet effective radius as expected, while relative cloud susceptibility is primarily dependent on optical thickness. Molecular absorption and reflecting surfaces both reduce the relative contribution of the cloud to the top-of-atmosphere (TOA) flux and therefore also reduce the TOA albedo susceptibility. Transmittance susceptibilities are negative with absolute values similar to albedo susceptibility, while atmospheric absorptance susceptibilities are about an order of magnitude smaller than albedo susceptibilities and can be either positive or negative. Observation-based susceptibility calculations are derived from MODIS pixel-level retrievals of liquid water cloud optical thickness, effective radius, and cloud top temperature; two data granule examples are shown. Susceptibility quantifies the aerosol indirect effect sensitivity in a way that can be easily computed from model fields. As such, susceptibilities derived from MODIS observations provide a higher-order test of model cloud properties used for indirect effect studies. MODIS-derived global distributions of cloud susceptibility and radiative forcing calculations are presented in a companion paper.

Platnick, Steven↗

Projected Applications of a "Climate in a Box" Computing System at the NASA Short-Term Prediction Research and Transition (SPoRT) Center

The NASA Short-term Prediction Research and Transition (SPoRT) Center focuses on the transition of unique observations and research capabilities to the operational weather community, with a goal of improving short-term forecasts on a regional scale. Advances in research computing have lead to "Climate in a Box" systems, with hardware configurations capable of producing high resolution, near real-time weather forecasts, but with footprints, power, and cooling requirements that are comparable to desktop systems. The SPoRT Center has developed several capabilities for incorporating unique NASA research capabilities and observations with real-time weather forecasts. Planned utilization includes the development of a fully-cycled data assimilation system used to drive 36-48 hour forecasts produced by the NASA Unified version of the Weather Research and Forecasting (WRF) model (NU-WRF). The horsepower provided by the "Climate in a Box" system is expected to facilitate the assimilation of vertical profiles of temperature and moisture provided by the Atmospheric Infrared Sounder (AIRS) aboard the NASA Aqua satellite. In addition, the Moderate Resolution Imaging Spectroradiometer (MODIS) instruments aboard NASA s Aqua and Terra satellites provide high-resolution sea surface temperatures and vegetation characteristics. The development of MODIS normalized difference vegetation index (NVDI) composites for use within the NASA Land Information System (LIS) will assist in the characterization of vegetation, and subsequently the surface albedo and processes related to soil moisture. Through application of satellite simulators, NASA satellite instruments can be used to examine forecast model errors in cloud cover and other characteristics. Through the aforementioned application of the "Climate in a Box" system and NU-WRF capabilities, an end goal is the establishment of a real-time forecast system that fully integrates modeling and analysis capabilities developed within the NASA SPoRT Center, with benefits provided to the operational forecasting community.

Jedlovec, Gary J.↗