Search NASA⌕ Search

SEARCH · Search NASA

Results for “data parameterization”

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 127 records · Page 7

Representing the Sub-Grid Heterogeneity of Surface Precipitation in A General Circulation Model

Precipitation variability on spatial scales smaller than a typical general circulation model (GCM) grid box is often neglected, with the grid-mean precipitation rate being applied uniformly to underlying surface tiles. This reduces the extrema seen by the surface, with corresponding reductions in surface runoff and altered land-atmosphere fluxes. Here we present a novel approach to stochastically distribute precipitation across sub-grid surface tiles within a GCM. Based on 4 km Stage IV precipitation data, the scheme parameterizes the dry area fraction as a function of grid mean precipitation rate, and defines the relative distribution of intensities across non-dry surface tiles. To incorporate memory and mimic the persistence of precipitating storms, the relative intensity assigned to each sub-grid tile is determined by an autoregressive process. Using single column experiments, the scheme is shown to reproduce observed precipitation statistics at the scale of model surface tiles. We also document impacts on surface hydrology and energy partitioning, with notable increases in precipitation runoff, surface temperature variance, and the Bowen ratio.

GCM↗

Parameterization of Vertical Cloud Distribution from C3M and MERRA Data Using ML Method

Clouds play a key role in regulating the hydrological cycle and the Earth's radiative energy budget. However, global climate models (GCMs) with a horizontal grid spacing on the order of 100 km have limitations in representing sub-grid cloud dynamics with spatial scales on the order of 1 km, leading to potential uncertainties in cloud radiative feedback on the global scale. In our research, we will leverage the capabilities of Deep Machine Learning (DML) methods to construct parameterizations of sub-grid volumetric cloud fraction (VCF), which is the frequency of occurrence on a grid volume accumulated in the horizontal and vertical directions. Our investigation delves into the intricate relationship between VCF obtained from the NASA CALIPSO-CloudSat-CERES-MODIS (CCCM) satellite observation data and 3-D MERRA-2 reanalysis meteorological profiling data (e.g., wind, relative humidity, temperature). Through a comprehensive one-year data training utilizing the Sequence to Sequence DML method, we have successfully disentangled the complicated cloud formation dynamics across diverse meteorological conditions through a day-to-day analysis framework. Preliminary findings reveal promising statistical agreements in geographical and vertical distributions and seasonal variations of volumetric cloud fraction between ML prediction and satellite measurements. These results underscore the aptitude of our DML model to discern underlying cloud physical processes and accurately represent sub-grid cloud formation dynamics. Additionally, we have also employed trained neural network to analyze uncertainties arising from errors in meteorological data, further enhancing the robustness of our VCF parameterization.

Shan Zeng↗

A statistical analysis of the dependency of closure assumptions in cumulus parameterization on the horizontal resolution

Simulated data from the UCLA cumulus ensemble model are used to investigate the quasi-universal validity of closure assumptions used in existing cumulus parameterizations. A closure assumption is quasi-universally valid if it is sensitive neither to convective cloud regimes nor to horizontal resolutions of large-scale/mesoscale models. The dependency of three types of closure assumptions, as classified by Arakawa and Chen, on the horizontal resolution is addressed in this study. Type I is the constraint on the coupling of the time tendencies of large-scale temperature and water vapor mixing ratio. Type II is the constraint on the coupling of cumulus heating and cumulus drying. Type III is a direct constraint on the intensity of a cumulus ensemble. The macroscopic behavior of simulated cumulus convection is first compared with the observed behavior in view of Type I and Type II closure assumptions using 'quick-look' and canonical correlation analyses. It is found that they are statistically similar to each other. The three types of closure assumptions are further examined with simulated data averaged over selected subdomain sizes ranging from 64 to 512 km. It is found that the dependency of Type I and Type II closure assumptions on the horizontal resolution is very weak and that Type III closure assumption is somewhat dependent upon the horizontal resolution. The influences of convective and mesoscale processes on the closure assumptions are also addressed by comparing the structures of canonical components with the corresponding vertical profiles in the convective and stratiform regions of cumulus ensembles analyzed directly from simulated data. The implication of these results for cumulus parameterization is discussed.

Xu, Kuan-Man↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Earth system dynamics↗

Timber Resources Inventory and Monitoring Joint Research Project

Primary objectives were to develop remote sensing analysis techniques for extracting forest related information from LANDSAT Multispectral Scanner (MMS) and Thematic Mapper data and to determine the extent to which International Paper Company information needs can be addressed with remote sensing information. The company actively manages 8.4 million acres of forest land. Traditionally, their forest inventories, updated on a three year cycle, are conducted through field surveys and aerial photography. The results reside in a digital forest data base containing 240 descriptive parameteres for individual forest stands. The information in the data base is used to develop seasonal and long range management strategies. Forest stand condition assessements (species composition, age, and density stratification) and identification of silvicultural activities (site preparation, planting, thinning, and harvest) are addressed.

Hill, C. L.↗

Utilization of satellite-derived cloud cover to improve the estimation of latent heat release in a winter extratropical cyclone

A simple scheme for utilizing satellite cloud pattern information to improve the estimation of parameterized stable and convective latent heat release is proposed and tested for an extratropical cyclone occurring over the United States during January 9-11, 1975. Calculations are done over the standard North American rawinsonde data region and over a sparse data region simulated by ignoring one-third of the original rawinsonde stations. Results indicate 15 percent and 30 percent improvements in latent heating estimates derived from the total and sparse data sets, respectively. A larger improvement was found using the sparse data set because the parameterized latent heat release estimates obtained using these data were poorer than those obtained using the total data set.

Lin, S.-J.↗

Surface radiation in the tropical Pacific

Monthly surface radiative fluxes have been calculated for the tropical Pacific between January 1970 and February 1978, using a radiative transfer parameterization. The radiative transfer parameterization included detailed treatments of the molecular and droplet absorptions, and surface and cloud reflections. The input data used in the calculations were obtained from the National Climatic Center (NCC), the National Center for Atmospheric Research (NCAR), and from the University of Hawaii. The results show that the distribution of surface radiation closely follows the distribution of cloudiness, and, to a lesser degree, humidity. The rms net error in the surface radiation estimates was about 15 W per sq m, with the largest contribution from uncertainties in the cloud cover and humidity data. The sensitivity of surface radiation parameterizations to input data errors is discussed, and some accuracy requirements for satellite retrievals of atmospheric and cloud parameters are proposed. The calculations are presented in the annually-averaged maps of surface radiation variations.

Chou, M.-D.↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Jet stream winds - Enhanced aircraft data acquisition and analysis over Southwest Asia

A project is described for providing the accurate initial and verification analyses for the jet stream in regions where general circulation models are known to have large systematic errors, either due to the extreme sparsity of data or to incorrect physical parameterizations. For this purpose, finely spaced aircraft-based meteorological data for the Southwest Asian regions collected for three 10-day periods in the winter of 1988-1989 will be used, together with corresponding data for the North American regions used as a control, to rerun the assimilation cycles and forecast models of the NMC and the NASA Goddard Laboratory for Atmospheres. Data for Southeast Asia will be collected by three carriers with extensive wide-body routes crossing the total region, while data for the North American region will be obtained from the archives of ACARS and GTS.

Tenenbaum, J.↗

Relativistic Three-Dimensional Lippmann-Schwinger Cross Sections for Space Radiation Applications

Radiation transport codes require accurate nuclear cross sections to compute particle fluences inside shielding materials. The Tripathi semi-empirical reaction cross section, which includes over 60 parameters tuned to nucleon-nucleus (NA) and nucleus-nucleus (AA) data, has been used in many of the world’s best-known transport codes. Although this parameterization fits well to reaction cross section data, the predictive capability of any parameterization is questionable when it is used beyond the range of the data to which it was tuned. Using uncertainty analysis, it is shown that a relativistic Three- Dimensional Lippmann-Schwinger (LS3D) equation model based on Multiple Scatter- ing Theory (MST) that uses 5 parameterizations—3 fundamental parameterizations to nucleon-nucleon (NN) data and 2 nuclear charge density parameterizations—predicts NA and AA reaction cross sections as well as the Tripathi cross section parameterization for reactions in which the kinetic energy of the projectile in the laboratory frame (TLab) is greater than 220 MeV/n. The relativistic LS3D model has the additional advantage of being able to predict highly accurate total and elastic cross sections. Consequently, it is recommended that the relativistic LS3D model be used for space radiation applications in which TLab > 220 MeV/n.

C M Werneth↗

Tests of a Semi-Analytical Case 1 and Gelbstoff Case 2 SeaWiFS Algorithm with a Global Data Set

A semi-analytical algorithm was tested with a total of 733 points of either unpackaged or packaged-pigment data, with corresponding algorithm parameters for each data type. The 'unpackaged' type consisted of data sets that were generally consistent with the Case 1 CZCS algorithm and other well calibrated data sets. The 'packaged' type consisted of data sets apparently containing somewhat more packaged pigments, requiring modification of the absorption parameters of the model consistent with the CalCOFI study area. This resulted in two equally divided data sets. A more thorough scrutiny of these and other data sets using a semianalytical model requires improved knowledge of the phytoplankton and gelbstoff of the specific environment studied. Since the semi-analytical algorithm is dependent upon 4 spectral channels including the 412 nm channel, while most other algorithms are not, a means of testing data sets for consistency was sought. A numerical filter was developed to classify data sets into the above classes. The filter uses reflectance ratios, which can be determined from space. The sensitivity of such numerical filters to measurement resulting from atmospheric correction and sensor noise errors requires further study. The semi-analytical algorithm performed superbly on each of the data sets after classification, resulting in RMS1 errors of 0.107 and 0.121, respectively, for the unpackaged and packaged data-set classes, with little bias and slopes near 1.0. In combination, the RMS1 performance was 0.114. While these numbers appear rather sterling, one must bear in mind what mis-classification does to the results. Using an average or compromise parameterization on the modified global data set yielded an RMS1 error of 0.171, while using the unpackaged parameterization on the global evaluation data set yielded an RMS1 error of 0.284. So, without classification, the algorithm performs better globally using the average parameters than it does using the unpackaged parameters. Finally, the effects of even more extreme pigment packaging must be examined in order to improve algorithm performance at high latitudes. Note, however, that the North Sea and Mississippi River plume studies contributed data to the packaged and unpackaged classess, respectively, with little effect on algorithm performance. This suggests that gelbstoff-rich Case 2 waters do not seriously degrade performance of the semi-analytical algorithm.

Carder, Kendall L.↗

A parameterization for the shortwave transmissivity of stratiform water clouds based on empirical data and radiative transfer theory

Surface measurements of solar flux and total integrated liquid-water content, radiosonde data, and infrared satellite images are analyzed in conjunction with radiative transfer calculations to derive an empirical parameterization for the shortwave transmissivity of continental stratiform water clouds. The data were collected near Denver, Colorado, over a period of six years. A multiple reflection radiative transfer model is employed to estimate stratus cloud transmissivity in terms of the measurable liquid-water path (LWP). A nonlinear fit of estimated transmissivities to the corresponding observations of LWP yields close agreement with a previous, more complicated parameterization. The derived expression for cloud transmissivity is used to predict mean daily surface fluxes for 61 days during which periods of stratiform clouds were observed over the Denver area. A comparison between predicted and measured fluxes shows agreement to within + or - 4 percent, with best agreement for clouds of moderate optical thickness. Potential sources of error are identified with sensitivity studies.

Derr, V. E.↗

Characteristics of the Martian atmosphere surface layer

Elements of various terrestrial boundary layer models are extended to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface ('constant flux') layer. The atmospheric surface layer consists of an interfacial sublayer immediately adjacent to the ground and an overlying fully turbulent surface sublayer where wind-shear production of turbulence dominates buoyancy production. Within the interfacial sublayer, sensible and latent heat are transported by non-steady molecular diffusion into small-scale eddies which intermittently burst through this zone. Both the thickness of the interfacial sublayer and the characteristics of the turbulent eddies penetrating through it depend on whether airflow is aerodynamically smooth or aerodynamically rough, as determined by the Roughness Reynold's number. Within the overlying surface sublayer, similarity theory can be used to express the mean vertical windspeed, temperature, and water vapor profiles in terms of a single parameter, the Monin-Obukhov stability parameter. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed using data from the TPRC Data Series and the first-order Chapman-Cowling expressions; the required collision integrals were approximated using the Lenard-Jones potential. Parameterizations for specific heat and binary diffusivity were also determined. The Brutsart model for sensible and latent heat transport within the interfacial sublayer for both aerodynamically smooth and rough airflow was experimentally tested under similar conditions, validating its application to Martian conditions. For the surface sublayer, the definition of the Monin-Obukhov length was modified to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. It was found that under most Martian conditions, the interfacial and surface sublayers offer roughly comparable resistance to sensible heat and water vapor transport and are thus both important in determining the associated fluxes.

Clow, G. D.↗

Enhanced Representation of Soil NO Emissions in the Community Multiscale Air Quality (CMAQ) Model Version 5.0.2

Modeling of soil nitric oxide (NO) emissions is highly uncertain and may misrepresent its spatial and temporal distribution. This study builds upon a recently introduced parameterization to improve the timing and spatial distribution of soil NO emission estimates in the Community Multiscale Air Quality (CMAQ) model. The parameterization considers soil parameters, meteorology, land use, and mineral nitrogen (N) availability to estimate NO emissions. We incorporate daily year-specific fertilizer data from the Environmental Policy Integrated Climate (EPIC) agricultural model to replace the annual generic data of the initial parameterization, and use a 12km resolution soil biome map over the continental USA. CMAQ modeling for July 2011 shows slight differences in model performance in simulating fine particulate matter and ozone from Interagency Monitoring of Protected Visual Environments (IMPROVE) and Clean Air Status and Trends Network (CASTNET) sites and NO2 columns from Ozone Monitoring Instrument (OMI) satellite retrievals. We also simulate how the change in soil NO emissions scheme affects the expected O3 response to projected emissions reductions.

remote sensing↗

Ground-based grasslands data to support remote sensing and ecosystem modeling of terrestrial primary production

Estimating terrestrial net primary production (NPP) using remote-sensing tools and ecosystem models requires adequate ground-based measurements for calibration, parameterization, and validation. These data needs were strongly endorsed at a recent meeting of ecosystem modelers organized by the International Geosphere-Biosphere Program's (IGBP's) Data and Information System (DIS) and its Global Analysis, Interpretation, and Modelling (GAIM) Task Force. To meet these needs, a multinational, multiagency project is being coordinated by the IGBP DIS to compile existing NPP data from field sites and to regionalize NPP point estimates to various-sized grid cells. Progress at Oak Ridge National Laboratory (ORNL) on compiling NPP data for grasslands as part of the IGBP DIS data initiative is described. Site data and associated documentation from diverse field studies are being acquired for selected grasslands and are being reviewed for completeness, consistency, and adequacy of documentation, including a description of sampling methods. Data are being compiled in a database with spatial, temporal, and thematic characteristics relevant to remote sensing and global modeling. NPP data are available from the ORNL Distributed Active Archive Center (DAAC) for biogeochemical dynamics. The ORNL DAAC is part of the Earth Observing System Data and Information System, of the US National Aeronautics and Space Administration.

Olson, R. J.↗

A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials

Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.

Kerisit, Sebastien N.↗

A comparison between satellite-defined and parameterized land-water differences in emitted longwave radiation

An analysis is performed to qualitatively compare the seasonal variation in emitted longwave radiation over land and over water areas as determined from 12 months of Nimbus 6 satellite data with that defined from parameterizations of this radiation budget component. These variations are noted when land and water surface areas are mapped to corresponding areas at the 'top' of the atmosphere. Variations of a surface-temperature-dependent parameterization of emitted longwave radiation originally suggested by Budyko (1969) are considered. The longwave radiation parameterizations indicate small differences between land and water profiles of emitted longwave radiation at the top of an atmospheric column in low latitudes in comparison to large differences in this feature shown to exist in the satellite data. The small differences are noted in linear parameterizations of emitted flux when zonally-averaged satellite data are used to define equation coefficients.

Stephenson-Graves, D.↗

Evaluation of Cloud Physical Properties of ECMWF Analysis and Re-Analysis (ERA-40 and ERA Interim) against CERES Tropical Deep Convective Cloud Object Observations

This study presents an approach that converts the vertical profiles of grid-averaged cloud properties from large-scale models to probability density functions (pdfs) of subgrid-cell cloud physical properties measured at satellite footprints. Cloud physical and radiative properties, rather than just cloud and precipitation occurrences, of assimilated cloud systems by the European Center for Medium-range Weather Forecasts (ECMWF) operational analysis (EOA) and ECMWF Re-Analyses (ERA-40 and ERA Interim) are validated against those obtained from Earth Observing System satellite cloud object data for January-August 1998 and March 2000 periods. These properties include ice water path (IWP), cloud-top height and temperature, cloud optical depth and solar and infrared radiative fluxes. Each cloud object, a contiguous region with similar cloud physical properties, is temporally and spatially matched with EOA and ERA-40 data. Results indicate that most pdfs of EOA and ERA-40 cloud physical and radiative properties agree with those of satellite observations of the tropical deep convective cloud-object type for the January-August 1998 period. There are, however, significant discrepancies in selected ranges of the cloud property pdfs such as the upper range of EOA cloud top height. A major discrepancy is that the dependence of the pdfs on the cloud object size for both EOA and ERA-40 is not as strong as in the observations. Modifications to the cloud parameterization in ECMWF that occurred in October 1999 eliminate the clouds near the tropopause but shift power of the pdf to lower cloud-top heights and greatly reduce the ranges of IWP and cloud optical depth pdfs. These features persist in ERA-40 due to the use of the same cloud parameterizations. The downgrade of data assimilation technique and the lack of snow water content information in ERA-40, not the coarser horizontal grid resolution, are also responsible for the disagreements with observed pdfs of cloud physical properties although the detection rates of cloud object occurrence are improved for small size categories. A possible improvement to the convective parameterization is to introduce a stronger dependence of updraft penetration heights with grid-cell dynamics. These conclusions will be rechecked using the ERA Interim data, due to recent changes in the ECMWF convective parameterization (Bechtold et al. 2004, 2008). Results from the ERA Interim will be presented at the meeting.

Xu, Kuan-Man↗