Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model uncertainty”

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 937 records · Page 52

Energy balance climate models

An introductory survey of the global energy balance climate models is presented with an emphasis on analytical results. A sequence of increasingly complicated models involving ice cap and radiative feedback processes are solved, and the solutions and parameter sensitivities are studied. The model parameterizations are examined critically in light of many current uncertainties. A simple seasonal model is used to study the effects of changes in orbital elements on the temperature field. A linear stability theorem and a complete nonlinear stability analysis for the models are developed. Analytical solutions are also obtained for the linearized models driven by stochastic forcing elements. In this context the relation between natural fluctuation statistics and climate sensitivity is stressed.

North, G. R.↗

The role of global cloud climatologies in validating numerical models

Global maps of the monthly mean net upward longwave radiation flux at the ocean surface were obtained for April, July, October 1985 and January 1986. These maps were produced by blending information obtained from a combination of general circulation model cloud radiative forcing fields, the top of the atmosphere cloud radiative forcing from ERBE and TOVS profiles and sea surface temperature on ISCCP C1 tapes. The fields are compatible with known meteorological regimes of atmospheric water vapor content and cloudiness. There is a vast area of high net upward longwave radiation flux (greater than 80/sq Wm) in the eastern Pacific Ocean throughout most of the year. Areas of low net upward longwave radiation flux ((less than 40/sq Wm) are the tropical convective regions and extra tropical regions that tend to have persistent low cloud cover.The technique used relies on General Circulation Model simulations and so is subject to some of the uncertainties associated with the model. However, all input information regarding temperature, moisture, and cloud cover is from satellite data having near global coverage. This feature of the procedure alone warrants its consideration for further use in compiling global maps of longwave radiation.

HARSHVARDHAN↗

Object-Oriented MDAO Tool with Aeroservoelastic Model Tuning Capability

An object-oriented multi-disciplinary analysis and optimization (MDAO) tool has been developed at the NASA Dryden Flight Research Center to automate the design and analysis process and leverage existing commercial as well as in-house codes to enable true multidisciplinary optimization in the preliminary design stage of subsonic, transonic, supersonic and hypersonic aircraft. Once the structural analysis discipline is finalized and integrated completely into the MDAO process, other disciplines such as aerodynamics and flight controls will be integrated as well. Simple and efficient model tuning capabilities based on optimization problem are successfully integrated with the MDAO tool. More synchronized all phases of experimental testing (ground and flight), analytical model updating, high-fidelity simulations for model validation, and integrated design may result in reduction of uncertainties in the aeroservoelastic model and increase the flight safety.

Pak, Chan-gi↗

Assessment of NASA GISS CMIP5 and Post-CMIP5 Simulated Clouds and TOA Radiation Budgets Using Satellite Observations. Part I: Cloud Fraction and Properties

Although many improvements have been made in phase 5 of the Coupled Model Intercomparison Project (CMIP5), clouds remain a significant source of uncertainty in general circulation models (GCMs) because their structural and optical properties are strongly dependent upon interactions between aerosol/cloud microphysics and dynamics that are unresolved in such models. Recent changes to the planetary boundary layer (PBL) turbulence and moist convection parameterizations in the NASA GISS Model E2 atmospheric GCM(post-CMIP5, hereafter P5) have improved cloud simulations significantly compared to its CMIP5 (hereafter C5) predecessor. A study has been performed to evaluate these changes between the P5 and C5 versions of the GCM, both of which used prescribed sea surface temperatures. P5 and C5 simulated cloud fraction (CF), liquid water path (LWP), ice water path (IWP), cloud water path (CWP), precipitable water vapor (PWV), and relative humidity (RH) have been compared to multiple satellite observations including the Clouds and the Earth's Radiant Energy System-Moderate Resolution Imaging Spectroradiometer (CERES-MODIS, hereafter CM), CloudSat- Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO; hereafter CC), Atmospheric Infrared Sounder (AIRS), and Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E). Although some improvements are observed in the P5 simulation on a global scale, large improvements have been found over the southern midlatitudes (SMLs), where correlations increased and both bias and root-mean-square error (RMSE) significantly decreased, in relation to the previous C5 simulation, when compared to observations. Changes to the PBL scheme have resulted in improved total column CFs, particularly over the SMLs where marine boundary layer (MBL) CFs have increased by nearly 20% relative to the previous C5 simulation. Globally, the P5 simulated CWPs are 25 gm22 lower than the previous C5 results. The P5 version of the GCM simulates PWV and RH higher than its C5 counterpart and agrees well with the AMSR-E and AIRS observations. The moister atmospheric conditions simulated by P5 are consistent with the CF comparison and provide a strong support for the increase in MBL clouds over the SMLs. Over the tropics, the P5 version of the GCM simulated total column CFs and CWPs are slightly lower than the previous C5 results, primarily as a result of the shallower tropical boundary layer in P5 relative to C5 in regions outside the marine stratocumulus decks.

General circulation models↗

Sensitivity to Factors Underlying the Hiatus

Recent trends in global mean surface air temperature fall outside the 90 range predicted by models using the CMIP5 forcings and scenarios; this recent period of muted warming is dubbed the hiatus. The hiatus has attracted broad attention in both the popular press and the scientific literature, primarily because of its perceived implications for understanding long-term trends. Many hypotheses have been offered to explain the warming slowdown during the hiatus, and comprehensive studies of this period across multiple variables and spatial scales will likely improve our understanding of the physical mechanisms driving global temperature change and variability.We argue, however, that decadal temperature trends by themselves are unlikely to constrain future trajectories of global mean temperature and that the hiatus does not significantly revise our understanding of overall climate sensitivity. Instead, we demonstrate that, because of the poorly constrained nature of the hiatus, model-observation disagreements over this period may be resolvable via uncertainties in the observations, modeled internal variability, forcing estimates, or (more likely) some combination of all three factors. We define the hiatus interval as 1998-2012, endpoints judiciously chosen to minimize observed warming by including the large 1998 El Nio event and excluding 2014, an exceptionally warm year. Such choices are fundamentally subjective and cannot be considered random, so any probabilistic statements regarding the likelihood of this occurring need to be made carefully. Using this definition, the observed global temperature trend estimates from four datasets fall outside the 5-95 interval predicted by the CMIP5 models. Here we explore some of the plausible explanations for this discrepancy, and show that no unique explanation is likely to fully account for the hiatus.

trends↗

Overview of the New Version 3 NASA Micro-Pulse Lidar Network (MPLNET) Automatic Precipitation Detection Algorithm

Precipitation modifies atmospheric column thermodynamics through the process of evaporation and serves as a proxy for latent heat modulation. For this reason, a correct precipitation parameterization (especially for low-intensity precipitation) within global scale models is crucial. In addition to improving our modeling of the hydrological cycle, this will reduce the associated uncertainty of global climate models in correctly forecasting future scenarios, and will enable the application of mitigation strategies. In this manuscript we present a proof of concept algorithm to automatically detect precipitation from lidar measurements obtained from the National Aeronautics and Space Administration Micropulse lidar network (MPLNET). The algorithm, once tested and validated against other remote sensing instruments, will be operationally implemented into the network to deliver a near real time (latency <1.5 h) rain masking variable that will be publicly available on MPLNET website as part of the new Version 3 data products. The methodology, based on an image processing technique, detects only light precipitation events (defined by intensity and duration) such as light rain, drizzle, and virga. During heavy rain events, the lidar signal is completely extinguished after a few meters in the precipitation or it is unusable because of water accumulated on the receiver optics. Results from the algorithm, in addition to filling a gap in light rain, drizzle, and virga detection by radars, are of particular interest for the scientific community as they help to fully characterize the aerosol cycle, from emission to deposition, as precipitation is a crucial meteorological phenomenon accelerating atmospheric aerosol removal through the scavenging effect. Algorithm results will also help the understanding of long term aerosol–cloud interactions, exploiting the multi-year database from several MPLNET permanent observational sites across the globe. The algorithm is also applicable to other lidar and/or ceilometer network infrastructures in the framework of the Global Aerosol Watch (GAW) aerosol lidar observation network (GALION).

Simone Lolli↗

Surface and aerosol models for use in radiative transfer codes

Absolute reflectance-based radiometric calibrations of Landsat-5 Thematic Mapper (TM) are improved with the inclusion of a method to invert optical-depth measurements to obtain aerosol-particle size distributions, and a non-Lambertian surface reflectance model. The inverted size distributions can predict radiances varying from the previously assumed jungian distributions by as much as 5 percent, though the reduction in the estimated error is less than one percent. Comparison with measured diffuse-to-global ratios show that neither distribution consistently predicts the ratio accurately, and this is shown to be a large contributor to calibration uncertainties. An empirical model for the surface reflectance of White Sands, using a two-degree polynomial fit as a function of scattering angle, was employed. The model reduced estimated errors in radiance predictions by up to one percent. Satellite calibrations dating from October, 1984 were reprocessed using the improved methods and linear estimations of satellite counts per unit radiance versus time since launch were determined which showed a decrease over time for the first four bands.

Hart, Quinn J.↗

A stochastic damage model for the rupture prediction of a multi-phase solid. I - Parametric studies. II - Statistical approach

A stochastic damage model for predicting the rupture of a brittle multiphase material is developed, based on the microcrack-macrocrack interaction. The model, which incorporates uncertainties in locations, orientations, and numbers of microcracks, characterizes damage by microcracking and fracture by macrocracking. A parametric study is carried out to investigate the change of the stress intensity at the macrocrack tip by the configuration of microcracks. The inherent statistical distribution of the fracture toughness arising from the intrinsic random nature of microcracks is explored using a statistical approach. For this purpose, a computer simulation model is introduced, which incorporates a statistical characterization of geometrical parameters of a random microcrack array.

Lua, Yuan J.↗

Assimilation of Gridded Terrestrial Water Storage Observations from GRACE into a Land Surface Model

Observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission have a coarse resolution in time (monthly) and space (roughly 150,000 km(sup 2) at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This work proposes a variant of existing ensemble-based GRACE-TWS data assimilation schemes. The new algorithm differs in how the analysis increments are computed and applied. Existing schemes correlate the uncertainty in the modeled monthly TWS estimates with errors in the soil moisture profile state variables at a single instant in the month and then apply the increment either at the end of the month or gradually throughout the month. The proposed new scheme first computes increments for each day of the month and then applies the average of those increments at the beginning of the month. The new scheme therefore better reflects submonthly variations in TWS errors. The new and existing schemes are investigated here using gridded GRACE-TWS observations. The assimilation results are validated at the monthly time scale, using in situ measurements of groundwater depth and soil moisture across the U.S. The new assimilation scheme yields improved (although not in a statistically significant sense) skill metrics for groundwater compared to the open-loop (no assimilation) simulations and compared to the existing assimilation schemes. A smaller impact is seen for surface and root-zone soil moisture, which have a shorter memory and receive smaller increments from TWS assimilation than groundwater. These results motivate future efforts to combine GRACE-TWS observations with observations that are more sensitive to surface soil moisture, such as L-band brightness temperature observations from Soil Moisture Ocean Salinity (SMOS) or Soil Moisture Active Passive (SMAP). Finally, we demonstrate that the scaling parameters that are applied to the GRACE observations prior to assimilation should be consistent with the land surface model that is used within the assimilation system.

GRACE↗

NASA GEOS-3/TRMM Re-analysis: Capturing Observed Tropical Rainfall Variability in Global Analysis for Climate Research

Understanding climate variability over a wide range of space-time scales requires a comprehensive description of the earth system. Global analyses produced by a fixed assimilation system (i.e., re-analyses) - as their quality continues to improve - have the potential of providing a vital tool for meeting this challenge. But at the present time, the usefulness of re-analyses is limited by uncertainties in such basic fields as clouds, precipitation, and evaporation - especially in the tropics, where observations are relatively sparse. Analyses of the tropics have long been shown to be sensitive to. the treatment of cloud precipitation processes, which remains a major source of uncertainty in current models. Yet, for many climate studies it is crucial that analyses can accurately reproduce the observed rainfall intensity and variability since a small error of 1 mm/d in surface rain translates into an error of approx. 30 W/sq m in energy (latent heat) flux. Currently, discrepancies between the observed and analyzed monthly-mean rain rates averaged to 100 km x 100 km resolution can exceed 4 mm/d (or 120 W/sq m ), compared to uncertainties in surface radiative fluxes of approx. 10-20 W/sq m . Improving precipitation in analyses would reduce a major source of uncertainty in the global energy budget. Uncertainties in tropical precipitation have also been a major impediment in understanding how the tropics interact with other regions, including the remote response to El Nino/Southern Oscillation (ENSO) variability on interannual time scales, the influence of Madden-Julian Oscillation (MJO) and monsoons on intraseasonal time scales. A global analysis that can replicate the observed precipitation variability together with physically consistent estimates of other atmospheric variables provides the key to breaking this roadblock. NASA Goddard Space Flight Center has been exploring the use of satellite-based microwave rainfall measurements in improving global analyses and has recently produced a multi-year, 1 x 1 TRMM re-analysis , which assimilates 6-hourly TMI and SSM/I surface rain rates over tropical oceans using a ID variational continuous assimilation (VCA) procedure in the GEOS-3 global data assimilation system. The analysis period extends from 1 November 1997 through 3 1 December 2002. The goal is to produce a multi-year global analysis that is dynamically consistent with available tropical precipitation observations for the community to assess its utility in climate applications and identify areas for further improvements. A distinct feature of the GEOS-3RRMh4 re-analysis is that its precipitation analysis is not derived from a short-term forecast (as done in most operational systems) but is given by a time- continuous model integration constrained by precipitation observations within a 6-h analysis window, while the wind, temperature, and pressure fields are allowed to directly respond to the improved precipitation and associated latent heating structures within the same analysis window. In this talk, I will assess the impact VCA precipitation assimilation on analyses of climate signals ranging from a few weeks to interannual time scales and compare results against other operational and reanalysis products.

Hou, Arthur Y.↗

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

Assessments of Hollow Cathode Wear in the Xenon Ion Propulsion System (XIPs(c)) by Numerical Analyses and Wear Tests

The standard approach presently followed by NASA to qualify electric propulsion for the required mission throughput has been based largely on life tests, which can be costly and time consuming. Revised electric propulsion lifequalification approaches are being formulated that combine analytical and/or computational methods with (shorter-duration) wear tests. As a model case, a wear test is being performed at JPL to assess the lifetime of the discharge hollow cathode in the Xenon Ion Propulsion System (XIPS(c)), a 25-cm ion engine developed by L-3 Communications Electron Technologies, Inc. for commercial applications. Wear and plasma data accumulated throughout this life-assessment program are being used to validate the existing 2-D hollow cathode code OrCa2D. We find that the OrCa2D steady-state solution predicts very well the time-averaged plasma data and the keeper voltage after 5500 hrs of operation in high-power mode. When the wave motion that occurs naturally in these devices is accounted for, based on an estimate of the maximum wave amplitude, the molybdenum-keeper erosion profile observed in the XIPS(c) discharge cathode is also reproduced within a factor of two of the observation. When the same model is applied to predict the erosion of a tantalum keeper we find that erosion is reduced by more than two orders of magnitude compared to the molybdenum keeper due the significantly lower sputtering yield of tantalum. A tantalum keeper would therefore allow keeper lifetimes that greatly exceed the present requirements for deep-space robotic missions considered by NASA. Moreover, such large reduction of the erosion renders the largest uncertainties in the models, which are associated with the wave amplitude estimates and the electron transport model, negligible.

plasma simulations↗

Developing Atmospheric Retrieval Methods for Direct Imaging Spectroscopy of Gas Giants in Reflected Light I: Methane Abundances and Basic Cloud Properties

Reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets are expected to be performed in the next decade by space-based telescopes equipped with optical wavelength coronagraphs and integral field spectrographs, such as the Wide-Field Infrared Survey Telescope (WFIRST). We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs an albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler (emcee) and a multimodal nested sampling algorithm (MultiNest) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model, and highlights possible discrepancies in the likelihood maps. Here we apply this methodology to simulated spectra of cool giant planets. As a proof-of-concept, our current atmospheric model contains 1 or 2 cloud layers, methane as a major absorber, and a H2-He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise, in the presence of spectral noise correlations. After internal validation, the method is applied to realistic reflected-light spectra of Jupiter, Saturn, and HD 99492 c, a likely observing target. We find that the presence or absence of clouds and methane can be determined with high accuracy, while parameters uncertainties are model-dependent.

spectroscopy↗

Assessment of the Quality of the Version 1.07 Temperature-Versus-Pressure Profiles of the Middle Atmosphere from TIMED/SABER

The quality of the retrieved temperature-versus-pressure (or T(p)) profiles is described for the middle atmosphere for the publicly available Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) Version 1.07 (V1.07) data set. The primary sources of systematic error for the SABER results below about 70 km are (1) errors in the measured radiances, (2) biases in the forward model, and (3) uncertainties in the corrections for ozone and in the determination of the reference pressure for the retrieved profiles. Comparisons with other correlative data sets indicate that SABER T(p) is too high by 1-3 K in the lower stratosphere but then too low by 1 K near the stratopause and by 2 K in the middle mesosphere. There is little difference between the local thermodynamic equilibrium (LTE) algorithm results below about 70 km from V1.07 and V1.06, but there are substantial improvements/differences for the non-LTE results of V1.07 for the upper mesosphere and lower thermosphere (UMLT) region. In particular, the V1.07 algorithm uses monthly, diurnally averaged CO2 profiles versus latitude from the Whole Atmosphere Community Climate Model. This change has improved the consistency of the character of the tides in its kinetic temperature (T(sub k)). The T(sub k) profiles agree with UMLT values obtained from ground-based measurements of column-averaged OH and O2 emissions and of the Na lidar returns, at least within their mutual uncertainties. SABER T(sub k) values obtained near the mesopause with its daytime algorithm also agree well with the falling sphere climatology at high northern latitudes in summer. It is concluded that the SABER data set can be the basis for improved, diurnal-to-interannual-scale temperatures for the middle atmosphere and especially for its UMLT region.

Remsberg, E. E.↗

Stochastic Analysis of Orbital Lifetimes of Spacecraft

A document discusses (1) a Monte-Carlo-based methodology for probabilistic prediction and analysis of orbital lifetimes of spacecraft and (2) Orbital Lifetime Monte Carlo (OLMC)--a Fortran computer program, consisting of a previously developed long-term orbit-propagator integrated with a Monte Carlo engine. OLMC enables modeling of variances of key physical parameters that affect orbital lifetimes through the use of probability distributions. These parameters include altitude, speed, and flight-path angle at insertion into orbit; solar flux; and launch delays. The products of OLMC are predicted lifetimes (durations above specified minimum altitudes) for the number of user-specified cases. Histograms generated from such predictions can be used to determine the probabilities that spacecraft will satisfy lifetime requirements. The document discusses uncertainties that affect modeling of orbital lifetimes. Issues of repeatability, smoothness of distributions, and code run time are considered for the purpose of establishing values of code-specific parameters and number of Monte Carlo runs. Results from test cases are interpreted as demonstrating that solar-flux predictions are primary sources of variations in predicted lifetimes. Therefore, it is concluded, multiple sets of predictions should be utilized to fully characterize the lifetime range of a spacecraft.

Sasamoto, Washito↗

Attribution of Flux Partitioning Variations Between Land Surface Models over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

uncertainty↗

Attribution of Flux Partitioning Variations Between Land Surface Models Over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

evapotranspiration partitioning; transpiration; so↗

Optimized Jupiter, Saturn, and Uranus interior models

A method for inverting gravitational data is the basis for the present models of Saturn, Jupiter, and Uranus, which have been computed to the third-order in rotational disturbance to the total potential and yield exact matches for recent, highly-accurate determinations of gravitational harmonics for these planets. These results indicate that a major He depletion has either not occurred in the outer layers of Jupiter or of Saturn, or that such a depletion's effect on the density profile has been masked by the presence of denser components. Considerable uncertainty characterizes Uranian models; the one presently preferred has substantial heavy-element enhancement at depth and little ice-rock separation.

Hubbard, W. B.↗