Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 343 records · Page 19

Carbon Monoxide Accumulation in the Extravehicular Mobility Unit

Introduction: Life support technology in large closed systems like submarines and space stations catalyzes carbon monoxide (CO) to carbon dioxide, which is easily removed. However, in a small system like the Extravehicular Mobility Unit (EMU), spacesuit, CO from exogenous (contaminated oxygen (O (sub 2) supply) and endogenous (human metabolism) sources will accumulate in the free suit volume. The free volume becomes a sink for CO that is rebreathed by the astronaut. The accumulation through time depends on many variables: the amount absorbed by the astronaut, the amount produced by the astronaut (between 0.28 and 0.34 moles per hour per kilogram)[1], the amount that enters the suit from contaminated O (sub 2), the amount removed through suit leak, the free volume of the suit, and the O (sub 2) partial pressure[2], just to list a few. Contamination of the EMU O (sub 2) supply with no greater than 1 part per million CO was the motivation for empirical measurements from CO pulse oximetry (SpCO) as well as mathematical modeling of the EMU as a rebreather for CO. Methods: We developed a first-order differential mixing equation as well as an iterative method to compute CO accumulation in the EMU. Pre-post measurements of SpCO (Rad-57, Masimo Corporation) from EMU ground training and on-orbit extravehicular activities (EVAs) were collected. Results: Initial modeling without consideration of the astronaut as a sink but only the source of CO showed that after 8 hours breathing 100 percent O (sub 2) with a 10 milliliter per minute (760 millimeters Hg at 21 degrees Centigrade standard) suit leak, an endogenous production rate of 0.23 moles per hour per kilogram for a 70 kilogram person with 42 liters (1.5 cubic feet) free suit volume resulted in a peak CO partial pressure (pCO) of 0.047 millimeters Hg at 4.3 pounds per square inch absolute (222 millimeters Hg). Preliminary results based on a 2008 model[3] with consideration of the astronaut as a sink and source of CO suggests that most of the rebreathed CO stays bound to hemoglobin and myoglobin and; therefore, pCO only increased to 0.002 millimeters Hg in the EMU. Hemoglobin saturation after 8 hours was an insignificant 0.4 percent compared to about 4 percent for cigarette smokers in the general population. This preliminary modeling result supplements 11 pre-post index finger SpCO measurements from EMU ground training (mean 0.5 percent versus 1 percent, probability equal to 0.41 from paired t-test) and 10 on-orbit pre-post EVAs (mean 1.5 percent versus 1.1 percent, probability equal to 0.17 from paired t-test) that showed no consistent increase, at least no increase outside the accuracy of the oximeter (1 percent display resolution with plus or minus 3 percent Standard Deviation). Discussion: Simulations continue, but a preliminary conclusion is that rebreathing endogenous CO accumulating in the EMU is not a serious medical issue during EVA. The absence of CO poisoning signs or symptoms following hundreds of EVAs is also good empirical evidence that corroborates the limited SpCO measurements and preliminary modeling results.

Conkin, J.↗

Trace Gas/Aerosol Interactions and GMI Modeling Support

Current global aerosol models use different physical and chemical schemes and parameters, different meteorological fields, and often different emission sources. Since the physical and chemical parameterization schemes are often tuned to obtain results that are consistent with observations, it is difficult to assess the true uncertainty due to meteorology alone. Under the framework of the NASA global modeling initiative (GMI), the differences and uncertainties in aerosol simulations (for sulfate, organic carbon, black carbon, dust and sea salt) solely due to different meteorological fields are analyzed and quantified. Three meteorological datasets available from the NASA DAO GCM, the GISS-II' GCM, and the NASA finite volume GCM (FVGCM) are used to drive the same aerosol model. The global sulfate and mineral dust burdens with FVGCM fields are 40% and 20% less than those with DAO and GISS fields, respectively due to its heavier rainfall. Meanwhile, the sea salt burden predicted with FVGCM fields is 56% and 43% higher than those with DAO and GISS, respectively, due to its stronger convection especially over the Southern Hemispheric Ocean. Sulfate concentrations at the surface in the Northern Hemisphere extratropics and in the middle to upper troposphere differ by more than a factor of 3 between the three meteorological datasets. The agreement between model calculated and observed aerosol concentrations in the industrial regions (e.g., North America and Europe) is quite similar for all three meteorological datasets. Away from the source regions, however, the comparisons with observations differ greatly for DAO, FVGCM and GISS, and the performance of the model using different datasets varies largely depending on sites and species. Global annual average aerosol optical depth at 550 nm is 0.120-0.131 for the three meteorological datasets.

Penner, Joyce E.↗

On the Role of Trend and Variability of Hydroxyl Radical (OH) in the Global Methane Budget

Decadal trends and interannual variations in the hydroxyl radical (OH), while poorly constrained at present, are critical for understanding the observed evolution of atmospheric methane (CH4). Through analyzing the OH fields simulated by the model ensemble of the Chemistry-Climate Model Initiative (CCMI), we find (1) the negative OH anomalies during the El Niño years mainly corresponding to the enhanced carbon monoxide (CO) emissions from biomass burning and (2) a positive OH trend during 1980-2010 dominated by the elevated primary production and the reduced loss of OH due to decreasing CO after 2000. Both two-box model inversions and variational 4D inversions suggest that ignoring the negative anomaly of OH during the El Niño years leads to a large overestimation of the increase in global CH4 emissions by up to10±3Tg yr-1to match the observed CH4 increase over these years. Not accounting for the increasing OH trends given by the CCMI models leads to an underestimation of the CH4 emission increase by 23±9Tg yr-1from 1986 to 2010. The variational inversion estimated CH4 emissions show that the tropical regions contribute most to the uncertainties related to OH. This study highlights the significant impact of climate and chemical feedbacks related to OH on the top-down estimates of the global CH4 budget.

Hydroxyl radical (OH)↗

Global High Resolution Crustal Magnetic Field at the Surface of the Moon from Low-Altitude Lunar Prospector Magnetic Gradient Data

We derived new vector gradients based models of crustal magnetic field at the lunar surface with data from the Lunar Prospector (LP) satellite using two model parameterization approaches: a global set of 35820 1° spaced (~30 km) equal area monopoles at 20 km below the surface (O’Brien and Parker, 1994; Olsen et al., 2017) and combined results of subsets of 100000 0.66° spaced monopoles at the same depth. We use the scheme of iteratively reweighted least-squares inversion to compute the initial model. Then the amplitudes of these monopoles are determined by minimizing the misfit to the components together with the global average of |Br| at the ellipsoid surface (i.e. applying a L1 model regularization of Br). In previous approaches using vector fields for modeling, we found that external field contamination leads to spurious anomalies in the downward continued field models even with stringent data selection criteria and ad-hoc noise removal techniques (e.g., satellite’s position in the Moon’s wake w.r.t. the solar wind and in the Earth’s magnetotail, internal/external dipoles fields removal, low-order polynomial removal, joint equivalent source cross-validation technique and visually removing remaining anomalous segments). On the other hand, with the use of gradients-only data (along-track first finite differences), we were able to completely bypass the ad-hoc techniques. Similar processing of Kaguya magnetic data, which have only higher altitude coverage over the most of the Moon except in the region of the South Pole–Aitken (SPA) basin, completely misses some of the anomalies seen in the Lunar Prospector data. The combined Lunar Prospector and Kaguya gradient-based models also severely degrade the derivation of the anomaly fields in many regions. Euler analysis of isolated anomaly features from the Reiner Gamma swirl suggests top depths of about 0.3 to 1.5 km and center depths of 10-14 km; in the region between Stein and Vallier craters north of the SPA basin our analysis suggests top depths of around 1.5 km and center depths of 13-15 km. With the spectral depth determination techniques, the SPA basin region yields depths to the base of magnetization ranging between 15 and 40 km. Three-dimensional modeling and the bulk magnetization determinations of the sources constrained by the Euler and spectral methods is underway.

Gradient↗

Improving Seasonal Forecast Using Probabilistic Deep Learning

The path toward realizing the potential of seasonal forecasting and its socioeconomic benefits relies on improving general circulation model (GCM) based dynamical forecast systems. To improve dynamical seasonal forecasts, it is crucial to set up forecast benchmarks, and clarify forecast limitations posed by model initialization errors, formulation deficiencies, and internal climate variability. With huge costs in generating large forecast ensembles, and limited observations for forecast verification, the seasonal forecast benchmarking and diagnosing task proves challenging. Here, we develop a probabilistic deep learning-based statistical forecast methodology, drawing on a wealth of climate simulations to enhance seasonal forecast capability and forecast diagnosis. By explicitly modeling the internal climate variability and GCM formulation differences, the proposed Conditional Generative Forecasting (CGF) methodology enables bypassing crucial barriers in dynamical forecast, and offers a top-down viewpoint to examine how complicated GCMs encode the seasonal predictability information. We apply the CGF methodology for global seasonal forecast of precipitation and 2 m air temperature, based on a unique data set consisting 52,201 years of climate simulation. Results show that the CGF methodology can faithfully represent the seasonal predictability information encoded in GCMs. We successfully apply this learned relationship in real-world seasonal forecast, achieving competitive performance compared to dynamical forecasts. Using this CGF as benchmark, we reveal the impact of insufficient forecast spread sampling that limits the skill of the considered dynamical forecast system. Finally, we introduce different strategies for composing ensembles using the CGF methodology, highlighting the potential for leveraging the strengths of multiple GCMs to achieve advantgeous seasonal forecast.

54 ENVIRONMENTAL SCIENCES↗

An Objective Verification of the North American Mesoscale Model for Kennedy Space Center and Cape Canaveral Air Force Station

The 45th Weather Squadron (45 WS) Launch Weather Officers (LWO's) use the 12-km resolution North American Mesoscale (NAM) model (MesoNAM) text and graphical product forecasts extensively to support launch weather operations. However, the actual performance of the model at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) has not been measured objectively. In order to have tangible evidence of model performance, the 45 WS tasked the Applied Meteorology Unit (AMU; Bauman et ai, 2004) to conduct a detailed statistical analysis of model output compared to observed values. The model products are provided to the 45 WS by ACTA, Inc. and include hourly forecasts from 0 to 84 hours based on model initialization times of 00, 06, 12 and 18 UTC. The objective analysis compared the MesoNAM forecast winds, temperature (T) and dew pOint (T d), as well as the changes in these parameters over time, to the observed values from the sensors in the KSC/CCAFS wind tower network shown in Table 1. These objective statistics give the forecasters knowledge of the model's strengths and weaknesses, which will result in improved forecasts for operations.

Bauman, William H., III↗

Local-scale Arctic tundra heterogeneity affects regional-scale carbon dynamics

In northern Alaska nearly 65% of the terrestrial surface is composed of polygonal ground, where geomorphic tundra landforms disproportionately influence carbon and nutrient cycling over fine spatial scales. Process-based biogeochemical models used for local to Pan-Arctic projections of ecological responses to climate change typically operate at coarse-scales (1km 2 –0.5°) at which fine-scale (<1km 2 ) tundra heterogeneity is often aggregated to the dominant land cover unit. Here, we evaluate the importance of tundra heterogeneity for representing soil carbon dynamics at fine to coarse spatial scales. We leveraged the legacy of data collected near Utqiagvik, Alaska between 1973 and 2016 for model initiation, parameterization, and validation. Simulation uncertainty increased with a reduced representation of tundra heterogeneity and coarsening of spatial scale. Hierarchical cluster analysis of an ensemble of 21 st -century simulations reveals that a minimum of two tundra landforms (dry and wet) and a maximum of 4km 2 spatial scale is necessary for minimizing uncertainties (<10%) in regional to Pan-Arctic modeling applications.

54 ENVIRONMENTAL SCIENCES↗

IDAES-PSE 2.5.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.5.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New diagnostics check for near-parallel variables and constraints. New diagnostics tools for identifying causes of infeasibility in models. New example for creating a custom model of a liquid-liquid extractor unit operation. Bug Fixes Fixed bug in Gibbs reactor that caused it to appear to have additional spurious degrees of freedom. Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Testing and Robustness Deployed the IDAES Diagnostics Toolbox to confirm that there are no structural or numerical issues in the core model libraries. Additional robustness tests for core model, and some associated improvements in the converge tester class. Fixed a number of issues that were causing unexpected warnings to be emitted during testing. Deprecations and Removals Removed examples for RIPE tool which has not been supported for a number of releases.

AS↗

Machine Learning for Rapid Biomarker Discovery via Image-Omic Fusion

In this project we proposed to accelerate biomarker discovery and improve disease prediction by combining image and omic data using a single, unified modelling approach. We hypothesized that much richer information could be derived from both the image and omic data when they are analyzed together. Using imaging mass spectrometry (IMS) data collected from mice we fused image and omic data to increase the resolution of the IMS data. Using human data provided to us by Oregon Health and Science University, we developed an initial modeling framework that can increase the accuracy of both imaging and MS datasets.

59 BASIC BIOLOGICAL SCIENCES↗

Synergism of Saturn, Enceladus and Titan and Formation of HCNO Exobiological Molecules

Saturn as a system has two very exotic moons Titan and Enceladus. Titan with energy input from Saturn's magnetosphere, solar UV irradiation, and cosmic rays can make HCN based molecules as discussed in earlier paper by [1]. Space radiation effects at both moons, and as coupled by the Saturn magnetosphere could cause an unexpected series of events leading to the evolution of biological models at Titan composed of HCNO with oxygen as the new ingredient. The "Old Faithful" model by [2] suggests that Enceladus, highly irradiated by Saturn magnetospheric electrons, has episodic ejections of water vapor driven by radiolytic oxidation gas products into Saturn's magnetosphere. At Titan Cassini discovered 1) that keV oxygen ions, evidently from Enceladus, are bombarding Titan's upper atmosphere [3] and 2) the discovery of heavy positive and negative ions within Titan's upper atmosphere [4]. Initial models of heavy ion formation in Titan's upper atmosphere invoked polymerization of aromatics such as Benzenes and their radicals to make PAHs [5], while a more recent model by [6] has raised the possibility of carbon chains forming from the polymerization of acetylene and its radicals to eventually make fullerenes. Laboratory measurements indicate that fullerenes, which are hollow carbon shells, can trap the keV oxygen and with the clustering of fullerenes and possible mixture with PAHs, some with nitrogen molecules, can make the larger aerosols with oxygen within them. Then with further ionizing irradiation from cosmic rays deep in the atmosphere "tholin" molecules are produced with all the molecular components present from which organic molecules can form. Among the molecular components are amino acids, the fundamental building blocks of life as we know it. This process maybe a common chemical pathway, both at the system level and at the molecular level, to form prebiotic and perhaps even biotic molecules. Such processes can be occurring throughout our universe, such as molecular clouds in the ISM.

Sittler, Edward C., Jr.↗

Quantum-Assisted Variational Segmentation for Image-to-Image Wildfire Detection Using Satellite Data

The quantum computing community has been searching for suitable applications to demonstrate the potential of near-term quantum devices. Quantum machine learning is a potential candidate, particularly using models that cannot be efficiently simulated with classical computers [1, 2]. This work focuses on a transition phase of quantum computers where the quantum machine learning model is still simulable classically but projected not to be simulable as the size of the model grows. Ultimately quantum computers may have advantages for high-dimensional real-world problems. Due to the limited number of qubits in current noisy intermediate-scale quantum (NISQ) devices, the direct application of quantum computers in high dimensional data is not feasible. To remedy this problem, an encoder-decoder architecture can be utilized. The encoder model would transform the high-dimensional data into a compact representation, to a level that small quantum computers can be used today (or in the near future), and the decoder would take the quantum processed outputs back to the high-dimensional space. Addressing the two challenges of quantum machine learning, this work investigates a hybrid supervised generative model with a quantum Ising Born machine embedded as the latent distribution. The model contains four main parts (Figure 1.a.): (1) a U-NET architecture responsible for learning segmentation flow, (2) a Prior network responsible for learning an encoded latent distribution of the input data, (3) a Born machine which represents the latent distribution, and (4) a Posterior network in charge of learning the joint encoded latent distribution of inputs and target data. The initial model, proposed by [3], is optimized by (1) maximizing the overlap of the prior and posterior latent distributions, and (2) minimizing the segmentation loss. The proposed model is designed to be investigated in a simulation environment applied to the real-world application of wildfire segmentation. Specifically, the model is designed to solve the patchy wildfire segmentations of Moderate Resolution Imaging Spectroradiometer (MODIS) by taking the MODIS observations and using Visible Infrared Imaging Radiometer Suite’s (VIIRS) consistent wildfire product as the target. The model solves patchy wildfire segmentations and provides insight into the epistemic errors sourced from model variation. The model utilizes the Born machine as a QUBO solver to represent the latent space as a Bernoulli distribution. The proposed configuration allows the variational segmentation model to leverage the true quantum probabilistic nature and derive a more expressive latent configuration, increasing the model performance in describing wildfire segmentations. The quantum probabilistic information of the Born machine is directly incorporated in the Kullback-Leibler divergence loss in the prior and posterior distributions, forcing the Bernoulli latent distribution to maximize the overlap of input and joint input-target distributions. The proposed model is then trained and compared with a baseline only consisting of direct Bernoulli latent distribution with no Born machine representing the latent space. The models are evaluated based on the segmentation metrics, such as precision, recall, intersect of union, with uncertainty boundaries accounting for the stochastic nature of the model. Our findings show that even in low latent-dimensional space (due to the limit in computational power of the classical quantum simulator), we are able to effectively capture the latent representation and hence the model performs better than the baseline. The findings are a projection for scaling the model into higher dimensional latent space with the Born machine surpassing the baseline performance. Figure 1. Sub-figure (a) demonstrates the architecture for the training phase. The model consists of a Prior and Posterior network that encode inputs and joint input-target data into compact representations, respectively. The Born machine represents the latent distribution, and the U-NET branch learns the segmentation patterns of the data. The stochasticity is introduced to the U-NET through its last layer to create meaningful but stochastic segmentations. Sub-figure (b) represents the inference phase where the model takes the stochastic behavior from the prior network and injects that into the U-NET. Each attempt of inference will generate different but similar segmentations from the same distribution of the wildfire event. REFERENCES [1] Coyle, B., Mills, D., Danos, V., & Kashefi, E. (2020). The Born supremacy: quantum advantage and training of an Ising Born machine. npj Quantum Information, 6(1), 1-11. [2] Liu, J. G., & Wang, L. (2018). Differentiable learning of quantum circuit born machines. Physical Review A, 98(6), 062324. [3] Kohl, S., Romera-Paredes, B., Meyer, C., De Fauw, J., Ledsam, J. R., Maier-Hein, K., ... & Ronneberger, O. (2018). A probabilistic u-net for segmentation of ambiguous images. Advances in neural information processing systems, 31.

quantum machine learning↗

Wrinkle ridges, reverse faulting, and the depth penetration of lithospheric stress in lunae planum, Mars

Tectonic features on a planetary surface are commonly used as constraints on models to determine the state of stress at the time the features formed. Quantitative global stress models applied to understand the formation of the Tharsis province on Mars constrained by observed tectonics have calculated stresses at the surface of a thin elastic shell and have neglected the role of vertical structure in influencing the predicted pattern of surface deformation. Wrinkle ridges in the Lunae Planum region of Mars form a conentric pattern of regularly spaced features in the eastern and southeastern part of Tharsis; they are formed due to compressional stresses related to the response of the Martian lithosphere to the Tharsis bulge. As observed in the exposures of valley walls in areas such as the Kasei Valles, the surface plains unit is underlain by an unconsolidated impact-generated megaregolith that grades with depth into structurally competent lithospheric basement. The ridges have alternatively been hypothesized to reflect deformation restricted to the surface plains unit ('thin skinned deformation') and deformation that includes the surface unit, megaregolith and basement lithosphere ('thick skinned deformation'). We have adopted a finite element approach to quantify the nature of deformation associated with the development of wrinkle ridges in a vertically stratified elastic lithosphere. We used the program TECTON, which contains a slippery node capability that allowed us to explicitly take into account the presence of reverse faults believed to be associated with the ridges. In this study we focused on the strain field in the vicinity of a single ridge when slip occurs along the fault. We considered two initial model geometries. In the first, the reverse fault was assumed to be in the surface plains unit, and in the second the initial fault was located in lithospheric basement, immediately beneath the weak megaregolith. We are interested in the conditions underwhich strain in the surface layer and basement either penetrates or fails to penetrate through the megaregolith. We thus address the conditions required for an initial basement fault to propagate through the megaregolith to the surface, as well as the effect of the megareolith on the strain tensor in the vicinity of a fault that nucleates in the surface plains unit.

Zuber, M. T.↗

A review of the thermonuclear runaway model of a nova outburst

Kraft's (1963) model for a nova outburst caused by a thermonuclear runaway in the hydrogen-rich matter accreted onto the white dwarf in a binary system is reviewed. Hydrostatic and hydrodynamic studies of this phenomenon are summarized, analyses of shock ejection based on hydrodynamic computer codes are discussed, and one specific hydrodynamic code is outlined. Results obtained with this code are presented and evaluated for an initial model containing a white dwarf with a hydrogen-rich envelope of 0.0001 to 0.001 solar mass in hydrostatic and thermal equilibrium. It is shown that an implicit hydrodynamic computer code is required in order to study the thermonuclear-runaway phenomenon. The early evolution of three models with different intrinsic luminosities is illustrated, and enhancement of CNO nuclei at the bottom of the hydrogen-rich envelope is investigated. It is suggested that the concentration of C-13, N-15, and O-17 should be greatly enhanced in nova ejecta.

Sparks, W. M.↗

A new approach to turbulence modeling

A new approach to Reynolds averaged turbulence modeling is proposed which has a computational cost comparable to two equation models but a predictive capability approaching that of Reynolds stress transport models. This approach isolates the crucial information contained within the Reynolds stress tensor, and solves transport equations only for a set of 'reduced' variables. In this work, Direct Numerical Simulation (DNS) data is used to analyze the nature of these newly proposed turbulence quantities and the source terms which appear in their respective transport equations. The physical relevance of these quantities is discussed and some initial modeling results for turbulent channel flow are presented.

Perot, B.↗

Tropical Rainfall Measuring Mission (TRMM) and the Future of Rainfall Estimation from Space

Tropical rainfall is important in the hydrological cycle and to the lives and welfare of humans. Three-fourths of the energy that drives the atmospheric wind circulation comes from the latent heat released by tropical precipitation. Recognizing the importance of rain in the tropics, NASA for the U.S.A. and NASDA for Japan have partnered in the design, construction and flight of a satellite mission to measure tropical rainfall and calculate the associated latent heat release. The Tropical Rainfall Measuring Mission (TRMM) satellite was launched on November 27, 1997, and data from all the instruments first became available approximately 30 days after launch. Since then, much progress has been made in the calibration of the sensors, the improvement of the rainfall algorithms and applications of these results to areas such as Data Assimilation and model initialization. TRMM has reduced the uncertainty of climatological rainfall in tropics by over a factor of two, therefore establishing a standard for comparison with previous data sets and climatologies. It has documented the diurnal variation of precipitation over the oceans, showing a distinct early morning peak and this satellite mission has shown the utility of precipitation information for the improvement of numerical weather forecasts and climate modeling. This paper discusses some promising applications using TRMM data and introduces a measurement concept being discussed by NASA/NASDA and ESA for the future of rainfall estimation from space.

Kakar, Ramesh↗

A study of GLAS and UCLA model forecasts of explosive maritime cyclogenesis

The explosive development of a maritime cyclone was examined with a numeric model of a cyclone which occurred during the FGGE Special Observing Period One in January 1979. The storm exhibited a Petterssen type-system frontal wave development. The model used was a nine-level, fourth-order coarse-mesh GLAS model initialized with FGGE level II-B data. Two versions of 72 hr synoptic-scale predictions were performed for the storm: with and without the input of meteorological satellite data. Attention was focused on the heat, mass and vorticity budgets of the storm, starting 24 hr before storm onset. The inclusion of satellite data, due to the poor quality of remotely sensed wind fields, degraded the performance of the model.

Peak, J. E.↗

GRO: Black hole models for gamma-ray bursts

The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRB's) is isotropic but is bound radially. This finding suggests that the bursts are either cosmological or they originate from an extended Galactic halo. The implied luminosities and the observed variability of the GRB's on time scales as short as one millisecond suggest that they originate from compact objects. We are presently studying black hole models for GRB's. Any such model must produce a non-thermal photon spectrum to agree with the observed properties. For a wide range of burst parameters the assumed bursting source consists of a non-thermal electron-positron-photon plasma of very high density. It seems possible to produce such a plasma in accretion onto black holes. In our on-going work, we are developing the kinetic theory for a non-equilibrium pair plasma. The main new features of our work are as follows: (1) We do not assume the presence of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an arbitrary concentration and energy distribution. (3) There is no soft photon source in our model; initially all the photons in the plasma are either energetic X-rays or gamma-rays. (4) The initial energy distribution of the pairs as well as photons is arbitrary. (5) We collect the analytical expressions for the kinetic kernels for all relevant processes. And (6) we present a different approach to finding the time-evolution of pair and photon spectra, which is a combination of the kinetic-theory and the non-linear Monte-Carlo schemes. We have developed many Monte-Carlo programs to model various process, to take into account the time evolution, and to incorporate various physical effects which are unique to non-thermal plasmas. The hydrodynamics of fireballs in GRB's was studied before. Applying results from kinetic theory will improve our understanding of these systems.

Ruderman, Malvin↗

Evaluation of MERRA-2-Based Ozone Profile Simulations with the Global Ozonesonde Network

Chemical transport model (CTM) hindcasts of ozone (O3) are useful for filling in observational gaps and providing context for observed O3 variability and trends. We use global networks of ozonesonde stations to evaluate the O3 profiles in two simulations running versions of the NASA Global Modeling Initiative (GMI) chemical mechanism. Both simulations are tied to the NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological reanalysis: 1) The GMI CTM, and 2) The MERRA-2 GMI Replay (M2 GMI). Both simulations start in 1980, and are compared against >50,000 ozonesonde profiles from 37 global stations from the tropics to the poles. The comparisons allow us to evaluate how the Replay technique affects modeled O3 distribution, how an updated chemical mechanism in the GMI CTM affects simulated tropospheric O3 amounts, and how observed O3 distributions compare to the full set of model output. In general, M2 GMI O3 is ~10% higher than in the GMI CTM, and shows global near-surface and tropical upper troposphere/lower stratosphere (UT/LS) high biases. The updated chemical mechanism in the GMI CTM reduces these high biases. Both simulations show similar negative biases in tropical free-tropospheric O3, especially during typical biomass burning seasons. The simulations are highly-correlated with ozonesonde measurements, particularly in the UT/LS (r > 0.8), showing the ability of MERRA-2 to capture tropopause height variations. Both simulations show improved correlations with ozonesonde data and smaller O3 biases in recent years. We expect to use the sonde/model comparisons to diagnose causes of disagreement and to gauge the feasibility of calculating multidecadal O3 trends from the model output.

Stauffer, Ryan↗