Search NASA⌕ Search

SEARCH · Search NASA

Results for “a priori estimates”

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 433 records · Page 24

Representation of Probability Density Functions from Orbit Determination using the Particle Filter

Statistical orbit determination enables us to obtain estimates of the state and the statistical information of its region of uncertainty. In order to obtain an accurate representation of the probability density function (PDF) that incorporates higher order statistical information, we propose the use of nonlinear estimation methods such as the Particle Filter. The Particle Filter (PF) is capable of providing a PDF representation of the state estimates whose accuracy is dependent on the number of particles or samples used. For this method to be applicable to real case scenarios, we need a way of accurately representing the PDF in a compressed manner with little information loss. Hence we propose using the Independent Component Analysis (ICA) as a non-Gaussian dimensional reduction method that is capable of maintaining higher order statistical information obtained using the PF. Methods such as the Principal Component Analysis (PCA) are based on utilizing up to second order statistics, hence will not suffice in maintaining maximum information content. Both the PCA and the ICA are applied to two scenarios that involve a highly eccentric orbit with a lower apriori uncertainty covariance and a less eccentric orbit with a higher a priori uncertainty covariance, to illustrate the capability of the ICA in relation to the PCA.

Mashiku, Alinda K.↗

Local Lunar Gravity Field Analysis over the South Pole-aitken Basin from SELENE Farside Tracking Data

We present a method with which we determined the local lunar gravity field model over the South Pole-Aitken (SPA) basin on the farside of the Moon by estimating adjustments to a global lunar gravity field model using SELENE tracking data. Our adjustments are expressed in localized functions concentrated over the SPA region in a spherical cap with a radius of 45deg centered at (191.1 deg E, 53.2 deg S), and the resolution is equivalent to a 150th degree and order spherical harmonics expansion. The new solution over SPA was used in several applications of geophysical analysis. It shows an increased correlation with high-resolution lunar topography in the frequency band l = 40-70, and admittance values are slightly different and more leveled when compared to other, global gravity field models using the same data. The adjustments expressed in free-air anomalies and differences in Bouguer anomalies between the local solution and the a priori global solution correlate with topographic surface features. The Moho structure beneath the SPA basin is slightly modified in our solution, most notably at the southern rim of the Apollo basin and around the Zeeman crater

Lunar Gravity↗

A first-order estimate of shock heating and vaporization in oceanic impacts

The vaporization of water in oceanic impacts of asteroids or comets of multikilometer dimensions is estimated by a semianalytical modeling of impact heating and shock isobar geometry that is based on computer code calculations. The mass of water vaporized in an infinitely deep ocean by the impact of a 10 km diameter asteroid at 25 km/sec (these values have been proposed for the Cretaceous/Tertiary extinction bolide) is approximately equal to the total mass of water vapor present in the earth's atmosphere, and 3-4 orders of magnitude larger than the mass of water vapor in the stratosphere. For projectiles of this size, however, the finite depth of the ocean becomes significant and may considerably reduce the amount of water vapor initially generated by the impact. Climatological models and extinction scenarios invoking the effects of impact-generated water vapor may critically depend on the a priori ambiguous details of the hypothesized impact.

Croft, S. K.↗

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

LiDAR-Inertial Based Navigation and Mapping for Precision Landing

Future lander missions will travel to ambitious, scientifically interesting locations near rough and dangerous terrain. They will need to operate with limited prior information about the terrain, and under varying lighting conditions. Landing safely and precisely in the face of these challenges is difficult for existing vision-based landing systems, which require detailed orbital reconnaissance, a priori hazard maps, and impose time-of-day restrictions on landing to ensure similar lighting conditions in orbital and descent imagery. Advanced 3D imaging LiDAR systems currently under development, and originally intended for single-scan hazard detection, have the potential to be operated continuously from altitudes of up to 5 km. Used together with existing inertial measurement units (IMUs), these sensors open a path-to-flight for a full navigation and mapping system, which could replace or augment a traditional landing sensor suite. A landing system based around these sensors can perform accurate altimetry, map-relative localization (MRL), LiDAR-inertial odometry, and map refinement in an illumination-insensitive manner, over unknown or partially known terrain. This paper outlines preliminary work on a LiDAR-inertial landing system that: estimates the spacecraft trajectory during entry, descent, and landing (EDL); and maps the topography of the terrain below, for future use in hazard detection and avoidance. An incremental, factor graph based, smoothing approach is used to solve for the maximum a posteriori trajectory of spacecraft states. Integrated IMU measurements and features tracked in adjacent range and intensity images are used to estimate motion (LiDAR-inertial odometry). LiDAR scans are binned into motion-corrected digital elevation models (DEMs), which are matched to an existing orbital topographic map to provide absolute position information (MRL). The estimated trajectory is then used to project the LiDAR scans into the map frame, creating a variable-resolution quadtree topographic map suitable for hazard detection and avoidance. Existing topographic maps from throughout the solar system (i.e., Earth, the Moon, Mars, Ceres, Vesta, Europa, Enceladus, and Eros) are upsampled for use in EDL simulations. The Mars 2020 Lander Vision System Simulator (LVSS) is extended to simulate LiDAR-inertial data for realistic EDL trajectories. Results of the algorithm operating on the simulated data are presented. Estimated spacecraft trajectory and refined map are compared to ground truth to assess estimation accuracy.

Katake, Anup↗

Variability in Surface BRDF at Different Spatial Scales (30 m-500 m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75 off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertain ties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure.

Roman, Miguel O.↗

Variability in Surface BRDF at Different Spatial Scales (30m-500m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle 1 remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75deg off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular 18 characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertainties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure.

Roman, Miguel O.↗

Extreme Lake-Effect Snow from a GPM Microwave Imager Perspective: Observational Analysis and Precipitation Retrieval Evaluation

This study focuses on the ability of the Global Precipitation Measurement (GPM) passive microwave sensors to detect and provide quantitative precipitation estimates (QPE) for extreme lake-effect snowfall events over the United States lower Great Lakes region. GPM Microwave Imager (GMI) high frequency channels can clearly detect intense shallow convective snowfall events. However, GMI Goddard PROfiling (GPROF) QPE retrievals produce inconsistent results when compared against the Multi-Radar/Multi-Sensor (MRMS) ground-based radar reference dataset. While GPROF retrievals adequately capture intense snowfall rates and spatial patterns of one event, GPROF systematically underestimates intense snowfall rates in another event. Furthermore, GPROF produces abundant light snowfall rates that do not conform with MRMS observations. Ad-hoc precipitation rate thresholds are suggested to partially mitigate GPROF’s overproduction of light snowfall rates. The sensitivity and retrieval efficiency of GPROF to key parameters (2-meter temperature, total precipitable water, and background surface type) used to constrain the GPROF a-priori retrieval database are investigated. Results demonstrate that typical lake-effect snow environmental and surface conditions, especially coastal surfaces, are underpopulated in the database and adversely affect GPROF retrievals. For the two presented case studies, using snow cover a-priori database in the locations of originally deemed as coastline improves retrieval. This study suggests that it is particularly important to have more accurate GPROF surface classifications and better representativeness of the a-priori databases to improve intense lake-effect snow detection and retrieval performance.

Lisa Milani↗

Lidar-Based Safe Site Relative Navigation

There has been a renewed focus in exploration of the lunar surface and maximizing scientific potential of such missions is made possible in part by minimizing the time required to set up operations; that is, reducing transit time on the surface by increasing the landing precision with respect to the intended target. Established SPLICE project precision landing requirements necessitate a navigation filter architecture and underlying models developed specifically with these needs in mind. To date, test flights to characterize SPLICE GNC system performance have not provided a means to divert from the a priori selected landing site due to hazardous conditions. With the inclusion of a new sensor, the HDL, coupled with safe landing site selection algorithms, GNC can divert from the originally planned trajectory and navigate relative to the new targeted landing site. This work presents a novel hazard relative measurement model and covariance transformation methodology that enable the navigation system to inform a safe-site relative guidance profile to meet project precision landing goals.

Navigation↗

GPM Mission Overview and U.S. Science Status

The Global Precipitation Measurement (GPM) Mission is an international satellite mission to unify and advance precipitation measurements from a constellation of research and operational sensors to provide "next-generation" precipitation products [1-2]. Water is fundamental to life on Earth. Knowing where and how much rain and snow falls globally is vital to understanding how weather and climate impact both our environment and Earth's water and energy cycles, including effects on agriculture, fresh water availability, and responses to natural disasters. Since rainfall and snowfall vary greatly from place to place and over time, satellites can provide more uniform observations of rain and snow around the globe than ground instruments, especially in areas where surface measurements are difficult. Relative to current global rainfall products, GPM data products will be characterized by: (l) more accurate instantaneous precipitation measurements (especially for light rain and cold-season solid precipitation), (2) more frequent sampling by an expanded constellation of domestic and international microwave radiometers including operational humidity sounders, (3) intercalibrated microwave brightness temperatures from constellation radiometers within a unified framework, and (4) physical-based precipitation retrievals from constellation radiometers using a common a priori cloud/hydrometeor database derived from GPM Core sensor measurements. The cornerstone of the GPM mission is the deployment of a Core Observatory in a unique 65 non-Sun-synchronous orbit to serve as a physics observatory and a reference standard to unify precipitation measurements by a constellation of dedicated and operational passive microwave sensors. The design of the GPM Core Observatory is an advancement of the Tropical Rainfall Measuring Mission (TRMM)'s highly successful rain-sensing package. The Core Observatory will carry a Ku/Ka-band Dual-frequency Precipitation Radar (DPR) and a multichannel (l0-183 GHz) GPM Microwave Radiometer (GMI). Since light rain and falling snow account for a significant fraction of precipitation occurrence in middle and high latitudes, the GPM instruments extend the capabilities of the TRMM sensors to detect falling snow, measure light rain, and provide, for the first time, quantitative estimates of microphysical properties of precipitation particles. The combined use of DPR and GMI measurements will place greater constraints on possible solutions to radiometer retrievals to improve the accuracy and consistency of precipitation retrievals from all constellation radiometers. The GMI uses 13 different microwave channels to observe energy from the different types of precipitation through clouds for estimating everything from heavy to light rain and for detecting falling snow. As the satellite passes over Earth, the GMI constantly scans a region 885 kilometers across. The Ball Aerospace and Technology Corporation built the GMI under contract with NASA Goddard Space Flight Center. The DPR provides three-dimensional information about precipitation particles derived from reflected energy by these particles at different heights within the cloud system. The two frequencies of the DPR also allow the radar to infer the sizes of precipitation particles and offer insights into a storm's physical characteristics. The Ka-band frequen~y scans across a region of 125 kilometers and is nested within the wider scan of the Ku-band frequency of 245 kilometers. The Japan Aerospace and Exploration Agency (JAXA) and Japan's National Institute of Information and Communications Technology (NICT) built the DPR. The Core Observatory satellite will fly at an altitude of 253 miles (407 kilometers) in a non-Sun-synchronous orbit that covers the Earth from 65 S to 65 N - from about the Antarctic Circle to the Arctic Circle. The GPM Core Observatory is being developed and tested at NASA Goddard Space Flight Center. Once complete, a Japanese H-lIA rocket will carry thPM Core Observatory into orbit from Tanegashima Island, Japan in 2014. The GPM constellation is envisioned to comprise 8 or more microwave sensors provided by partners, including both conical imagers and cross-track sounders. GPM is currently a partnership between NASA and the Japan Aerospace Exploration Agency (JAXA). Additional partnerships are under development to include microwave radiometers on the French-Indian Megha-Tropiques satellite and U.S. Defense Meteorological Satellite Program (DMSP) satellites, as well as humidity sounders or precipitation sensors on operational satellites such as the National Polar-orbiting Operational Environmental Satellite System (NPOESS) Preparatory Project (NPP), NOAA-NASA Joint Polar Satellite System (JPSS) satellites, European MetOp satellites, and DMSP follow-on sensors. In addition, data from Chinese and Russian microwave radiometers may be available through international cooperation under the auspices of the Committee on Earth Observation Satellites (CEOS) and Group on Earth Observations (GEO). GPM's next-generation global precipitation data will lead to scientific advances and societal benefits in the following areas: (1) Improved knowledge of the Earth's water cycle and its link to climate change (2) New insights into precipitation microphysics, storm structures and large-scale atmospheric processes (3) Better understanding of climate sensitivity and feedback processes (4) Extended capabilities in monitoring and predicting hurricanes and other extreme weather events (5) Improved forecasting capabilities for natural hazards, including floods, droughts and landslides. (6) Enhanced numerical prediction skills for weather and climate (7) Better agricultural crop forecasting and monitoring of freshwater resources. An overview of the GPM mission concept and science activities in the United States, together with an update on international collaborations in radiometer intercalibration and ground validation, will be presented.

Hou, Arthur Y.↗

Simplified Approach for Scoping Assessment of Non-LWR Source Terms

This report describes a structure to aid in evaluation of release mitigation strategies across a range of reactor technologies. The assessment performed for example reactor concepts utilizes previous studies of postulated accident sequences for each reactor concept. This simplified approach classifies release mitigation strategies based on a range of barriers, physical attenuation processes, and system performance. It is not, however, intended to develop quantitative estimates of radiological release magnitudes and compositions to the environment. Rather, this approach is intended to identify the characteristics of a reactor design concepts release mitigation strategies that are most important to different classes of accident scenarios. It uses a scoping methodology to provide an approximate, order-of-magnitude, estimate of the radiological release to the environment and associated off-site consequences. This scoping method is applied to different reactor concepts, considering the performance of barriers to fission product release for these concepts under sample accident scenarios. The accident scenarios and sensitivity evaluations are selected in this report to evaluate the role of different fission product barriers in ameliorating the source term to the environment and associated off-site consequences. This report applies this structure to characterize how release mitigation measures are integrated to define overall release mitigation strategies for High Temperature Gas Reactors (HTGRs), Sodium Fast Reactors (SFRs), and liquid fueled Molten Salt Reactors (MSRs). To support this evaluation framework, factors defining a chain of release attenuation stages, and thus an overall mitigation strategy, must be established through mechanistic source term calculations. This has typically required the application of an integral plant analysis code such as MELCOR. At present, there is insufficient evidence to support a priori evaluation of the effectiveness of a release mitigation strategy for advanced reactor concepts across the spectrum of events that could challenge the radiological containment function. While it is clear that these designs have significant margin to radiological release to the environment for the scenarios comprising the design basis, detailed studies have not yet been performed to assess the risk profile for these plants. Such studies would require extensive evaluation across a reasonably complete spectrum of accident scenarios that could lead to radiological release to the environment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A new continuum model for suspensions of gyrotactic micro-organisms

A new continuum model is formulated for dilute suspensions of swimming micro-organisms with asymmetric mass distributions. Account is taken of randomness in a cell's swimming direction, p, by postulating that the probability density function for p satisfies a Fokker-Planck equation analogous to that obtained for colloid suspensions in the presence of rotational Brownian motion. The deterministic torques on a cell, viscous and gravitational, are balanced by diffusion, represented by an isotropic rotary diffusivity Dr, which is unknown a priori, but presumably reflects stochastic influences on the cell's internal workings. When the Fokker-Planck equation is solved, macroscopic quantities such as the average cell velocity Vc, the particle diffusivity tensor D and the effective stress tensor sigma can be computed; Vc and D are required in the cell conservation equation, and sigma in the momentum equation. The Fokker-Planck equation contains two dimensionless parameters, lambda and epsilon; lambda is the ratio of the rotary diffusion time Dr-1 to the torque relaxation time B (balancing gravitational and viscous torques), while epsilon is a scale for the local vorticity or strain rate made dimensionless with B. In this paper we solve the Fokker-Planck equation exactly for epsilon = 0 (lambda arbitrary) and also obtain the first-order solution for small epsilon. Using experimental data on Vc and D obtained with the swimming alga, Chlamydomonas nivalis, in the absence of bulk flow, the epsilon = 0 results can be used to estimate the value of lambda for that species (lambda approximately 2.2; Dr approximately 0.13 s-1). The continuum model for small epsilon is then used to reanalyse the instability of a uniform suspension, previously investigated by Pedley, Hill & Kessler (1988). The only qualitatively different result is that there no longer seem to be circumstances in which disturbances with a non-zero vertical wavenumber are more unstable than purely horizontal disturbances. On the way, it is demonstrated that the only significant contribution to sigma, other than the basic Newtonian stress, is that derived from the stresslets associated with the cells' intrinsic swimming motions.

NASA Discipline Number 28-20↗

Hybrid phenology matching model for robust crop phenological retrieval

Crop phenology regulates seasonal agroecosystem carbon, water, and energy exchanges, and is a key component in empirical and process-based crop models for simulating biogeochemical cycles of farmlands, assessing gross and net primary production, and forecasting the crop yield. The advances in phenology matching models provide a feasible means to monitor crop phenological progress using remote sensing observations, with a priori information of reference shapes and reference phenological transition dates. Yet the underlying geometrical scaling assumption of models, together with the challenge in defining phenological references, hinders the applicability of phenology matching in crop phenological studies. The objective of this study is to develop a novel hybrid phenology matching model to robustly retrieve a diverse spectrum of crop phenological stages using satellite time series. The devised hybrid model leverages the complementary strengths of phenometric extraction methods and phenology matching models. It relaxes the geometrical scaling assumption and can characterize key phenological stages of crop cycles, ranging from farming practice-relevant stages (e.g., planted and harvested) to crop development stages (e.g., emerged and mature). To systematically evaluate the influence of phenological references on phenology matching, four representative phenological reference scenarios under varying levels of phenological calibrations in terms of time and space are further designed with publicly accessible phenological information. The results indicate that the hybrid phenology matching model can achieve high accuracies for estimating corn and soybean phenological growth stages in Illinois, particularly with the year- and region-adjusted phenological reference (R-squared higher than 0.9 and RMSE less than 5 days for most phenological stages). The inter-annual and regional phenological patterns characterized by the hybrid model correspond well with those in the crop progress reports (CPRs) from the USDA National Agricultural Statistics Service (NASS). Compared to the benchmark phenology matching model, the hybrid model is more robust to the decreasing levels of phenological reference calibrations, and is particularly advantageous in retrieving crop early phenological stages (e.g., planted and emerged stages) when the phenological reference information is limited. This innovative hybrid phenology matching model, together with CPR-enabled phenological reference calibrations, holds 3 considerable promise in revealing spatio-temporal patterns of crop phenology over extended geographical regions.

Phenology↗

A new method for determining the CO to H2 conversion factor for translucent clouds

In this paper we describe a new technique for obtaining the conversion factor between the molecular hydrogen column density and the CO(J = 1-0) integrated antenna temperature. This factor, typically known as X(sub CO) is often to be of order a few times 10(exp 20)/sq cm/K km/s) for the molecular clouds in the Galaxy and is one of the primary means of determining the molecular cloud mass from CO observations. However, for the low-extinction interstellar clouds known as the translucent molecular clouds, estimates of X(sub CO) vary by up to a factor of 60 depending on the object and techniques employed to calibrate X(sub CO). Since the cloud mass is directly proportional to X(sub CO) uncertainties in mass estimates of translucent clouds can be more than an order of magnitude. We calibrate the H2 content in translucent clouds by using the linear relationship between the CH and H2 column densities. The CH column density is readily determined from observations of the CH ground-state hyperfine main-line transition at 3335 MHz. Using CH as a surrogate tracer for H2 and CO(J = 1-0) observations of a sample of translucent and dark molecular clouds, we find a wide variation in values for X(sub CO). For translucent clouds, X(sub CO) ranges from 0.3 to 6.8 x 10(exp 20) and for dark clouds the values range from 0.8 to 8.6. Although the average values for both types of cloud are similar to the canonical value determined for the Galactic molecular cloud ensemble (2-4 x 10(exp 20)), the scatter in individual X(sub CO) values may indicate that X(sub CO) for a given translucent cloud cannot be determined a priori and must be obtained for each cloud so that a reliable mass determination may be made.

Magnani, L.↗

Evaluation of Precipitation Detection over Various Surfaces from Passive Microwave Imagers and Sounders

During the middle part of this decade a wide variety of passive microwave imagers and sounders will be unified in the Global Precipitation Measurement (GPM) mission to provide a common basis for frequent (3 hr), global precipitation monitoring. The ability of these sensors to detect precipitation by discerning it from non-precipitating background depends upon the channels available and characteristics of the surface and atmosphere. This study quantifies the minimum detectable precipitation rate and fraction of precipitation detected for four representative instruments (TMI, GMI, AMSU-A, and AMSU-B) that will be part of the GPM constellation. Observations for these instruments were constructed from equivalent channels on the SSMIS instrument on DMSP satellites F16 and F17 and matched to precipitation data from NOAA's National Mosaic and QPE (NMQ) during 2009 over the continuous United States. A variational optimal estimation retrieval of non-precipitation surface and atmosphere parameters was used to determine the consistency between the observed brightness temperatures and these parameters, with high cost function values shown to be related to precipitation. The minimum detectable precipitation rate, defined as the lowest rate for which probability of detection exceeds 50%, and the detected fraction of precipitation, are reported for each sensor, surface type (ocean, coast, bare land, snow cover) and precipitation type (rain, mix, snow). The best sensors over ocean and bare land were GMI (0.22 mm/hr minimum threshold and 90% of precipitation detected) and AMSU (0.26 mm/hr minimum threshold and 81% of precipitation detected), respectively. Over coasts (0.74 mm/hr threshold and 12% detected) and snow-covered surfaces (0.44 mm/hr threshold and 23% detected), AMSU again performed best but with much lower detection skill, whereas TMI had no skill over these surfaces. The sounders (particularly over water) benefited from the use of re-analysis data (vs. climatology) to set the a-priori atmospheric state and all instruments benefit from the use of a conditional snow cover emissivity database over land. It is recommended that real-time sources of these data be used in the operational GPM precipitation algorithms.

GPM↗

A Priori Analysis of Flamelet-Based Modeling for a Dual-Mode Scramjet Combustor

An a priori investigation of the applicability of flamelet-based combustion models to dual-mode scramjet combustion was performed utilizing Reynolds-averaged simulations (RAS). For this purpose, the HIFiRE Direct Connect Rig (HDCR) flowpath, fueled with a JP-7 fuel surrogate and operating in dual- and scram-mode was considered. The chemistry of the JP-7 fuel surrogate was modeled using a 22 species, 18-step chemical reaction mechanism. Simulation results were compared to experimentally-obtained, time-averaged, wall pressure measurements to validate the RAS solutions. The analysis of the dual-mode operation of this flowpath showed regions of predominately non-premixed, high-Damkohler number, combustion. Regions of premixed combustion were also present but associated with only a small fraction of the total heat-release in the flow. This is in contrast to the scram-mode operation, where a comparable amount of heat is released from non-premixed and premixed combustion modes. Representative flamelet boundary conditions were estimated by analyzing probability density functions for temperature and pressure for pure fuel and oxidizer conditions. The results of the present study reveal the potential for a flamelet model to accurately model the combustion processes in the HDCR and likely other high-speed flowpaths of engineering interest.

Quinlan, Jesse R.↗

Evaluating Long-Term Changes in Atmospheric Ozone

Ozonesondes have made inexpensive, accurate measurements of ozone from the ground to 30km for more than 50 years. The data are used extensively for trend analyses and for evaluation of satellite and model data products, and are also part of climatologies that are used as a priori data for satellite retrievals. They are essential as a transfer standard when merging shorter satellite-derived time series, and are the most important source of trend-quality long-term records below about 18 km. The importance of long-term ozonesonde records as a stable reference has led to increased attention to quantifying uncertainties and changes in ozonesonde data. Based on past intercomparison data, ECC sondes show a modest (~1-5%) high bias with respect to UV-absorption measurements (including MOZAIC-IAGOS) in the troposphere, with an uncertainty of 5%, but no evidence of a change with time. Other sonde types show an increase of 5-20% in sensitivity to tropospheric ozone from 1970-1995. Agreement in the stratosphere is much better. In the past 20 years ozonesonde precision has improved by a factor of 2, primarily through the adoption of strict standard operating procedures. In addition, many of the most important ozonesonde records have been re-evaluated and corrected, and detailed estimates made of their biases and uncertainties. Corrections to historical data for known issues may reduce biases but simultaneously introduce additional uncertainties. The uncertainty budget for the ozone partial pressure reading has contributions from stoichiometry, cell background current, pump efficiency and temperature, sensing solution type and volume. The much-discussed “background current” has recently been shown to be primarily related to reaction stoichiometry. Several quality assurance issues remain, but are tractable problems that can be addressed with further research. This will be required if the present goal of better than 5% overall uncertainty throughout the global ozonesonde network is to be achieved. Finally, the importance of regular sonde intercomparisons, employing UV photometers traceable to the modern UV-absorption standard, and of regular comparison of operational station records with multiple satellite sensors, is emphasized.

Ozonesondes↗

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗