SEARCH · Search NASA

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Exploring Stratocumulus Cloud-Top Entrainment Processes and Parameterizations by Using Doppler Cloud Radar Observations

Observations from an upward-pointing Doppler cloud radar are used to examine cloud-top entrainment processes and parameterizations in a non-precipitating continental stratocumulus cloud deck maintained by time varying surface buoyancy fluxes and cloud-top radiative cooling. Radar and ancillary observations were made at the Atmospheric Radiation Measurement (ARM)’s Southern Great Plains (SGP) site located near Lamont, Oklahoma of unbroken, non-precipitating stratocumulus clouds observed for a 14-hour period starting 0900 Central Standard Time on 25 March 2005. The vertical velocity variance and energy dissipation rate (EDR) terms in a parameterized turbulence kinetic energy (TKE) budget of the entrainment zone are estimated using the radar vertical velocity and the radar spectrum width observations from the upward-pointing millimeter cloud radar (MMCR) operating at the SGP site. Hourly averages of the vertical velocity variance term in the TKE entrainment formulation correlates strongly (r=0.72) to the dissipation rate term in the entrainment zone. However, the ratio of the variance term to the dissipation decreases at night due to decoupling of the boundary layer. When the night -time decoupling is accounted for, the correlation between the variance and the EDR term increases (r=0.92). To obtain bulk coefficients for the entrainment parameterizations derived from the TKE budget, independent estimate of entrainment were obtained from an inversion height budget using ARM SGP observations of the local time derivative and the horizontal advection of the cloud-top height. The large-scale vertical velocity at the inversion needed for this budget from EMWF reanalysis. This budget gives a mean entrainment rate for the observing period of 0.76±0.15 cm/s. This mean value is applied to the TKE budget parameterizations to obtain the bulk coefficients needed in these parameterizations. These bulk coefficients are compared with those from previous and are used to in the parameterizations to give hourly estimates of the entrainment rates using the radar derived vertical velocity variance and dissipation rates. Hourly entrainment rates were estimated from a convective velocity w* parameterization depends on the local surface buoyancy fluxes and the calculated radiative flux divergence, parameterization using a bulk coefficient obtained from the mean inversion height budget. The hourly rates from the cloud turbulence estimates and the w* parameterization, which is independent of the radar observations, are compared with the hourly we values from the budget. All show rough agreement with each other and capture the entrainment variability associated with substantial changes in the surface flux and radiative divergence at cloud top. Major uncertainties in the hourly estimates from the height budget and w* are discussed. The results indicate a strong potential for making entrainment rate estimates directly from the radar vertical velocity variance and the EDR measurements—a technique that has distinct advantages over other methods for estimating entrainment rates. Calculations based on the EDR alone can provide high temporal resolution (for averaging intervals as small as 10 minutes) of the entrainment processes and do not require an estimate of the boundary layer depth, which can be difficult to define when the boundary layer is decoupled.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

On the influence of over-parameterization in manifold based surrogates and deep neural operators

Constructing accurate and generalizable approximators (surrogate models) for complex physico-chemical processes exhibiting highly non-smooth dynamics is challenging. The main question is what type of surrogate models we should construct and should these models be under-parameterized or over-parameterized. In this work, we propose new developments and perform comparisons for two promising approaches: manifold-based polynomial chaos expansion (m-PCE) and the deep neural operator (DeepONet), and we examine the effect of over-parameterization on generalization. While m-PCE enables the construction of a mapping by first identifying low-dimensional embeddings of the input functions, parameters, and quantities of interest (QoIs), a neural operator learns the nonlinear mapping via the use of deep neural networks. Here, we demonstrate the performance of these methods in terms of generalization accuracy by solving the 2D time-dependent Brusselator reaction-diffusion system with uncertainty sources, modeling an autocatalytic chemical reaction between two species. We first propose an extension of the m-PCE by constructing a mapping between latent spaces formed by two separate embeddings of the input functions and the output QoIs. To further enhance the accuracy of the DeepONet, we introduce weight self-adaptivity in the loss function. We demonstrate that the performance of m-PCE and DeepONet is comparable for cases of relatively smooth input-output mappings. However, when highly non-smooth dynamics is considered, DeepONet shows higher approximation accuracy. We also find that for m-PCE, modest over-parameterization leads to better generalization, both within and outside of distribution, whereas aggressive over-parameterization leads to over-fitting. In contrast, an even highly over-parameterized DeepONet leads to better generalization for both smooth and non-smooth dynamics. Furthermore, we compare the performance of the above models with another recently proposed operator learning model, the Fourier Neural Operator, and show that its over-parameterization also leads to better generalization. Taken together, our studies show that m-PCE can provide very good accuracy at very low training cost, whereas a highly over-parameterized DeepONet can provide better accuracy and robustness to noise but at higher training cost. In both methods, the inference cost is negligible.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Interparameter tradeoff quantification for isotropic-elastic full-waveform inversion with various model parameterizations

Simultaneous determination of multiple physical parameters using full-waveform inversion (FWI) suffers from interparameter tradeoff diffculty. Analyzing the interparameter tradeoffs in different model parameterizations of isotropic-elastic FWI, and thus determining the appropriate model parameterization, are critical for efficient inversion and obtaining reliable inverted models. Five different model parameterizations are considered and compared including velocity-density, modulus-density, impedance-density and two velocity-impedance parameterizations. The scattering radiation patterns are first used for interparameter tradeoff analysis. Furthermore, a new framework is developed to evaluate the interparameter tradeoff based upon multiparameter Hessian-vector products: multiparameter point spread functions (MPSFs) and interparameter contamination sensitivity kernels (ICSKs), which provide quantitative, second-order measurements of the interparameter contaminations. In the numerical experiments, the interparameter tradeoffs in various model parameterizations are evaluated using the MPSFs and ICSKs. Inversion experiments are then carried out with simple Gaussian-anomaly models and complex Marmousi model. Overall, the parameterization of P-wave velocity, S-wave velocity and density and the parameterization of P-wave velocity, S-wave velocity and S-wave impedance perform best for reconstructing all of the physical parameters. Isotropic-elastic FWI of the Hussar low frequency dataset with various model parameterizations verifies the conclusions.

58 GEOSCIENCES↗

DOE OSTI2023

Memory-based parameterization with differentiable solver: Application to Lorenz ’96

Physical parameterizations (or closures) are used as representations of unresolved subgrid processes within weather and global climate models or coarse-scale turbulent models, whose resolutions are too coarse to resolve small-scale processes. These parameterizations are typically grounded on physically based, yet empirical, representations of the underlying small-scale processes. Machine learning-based parameterizations have recently been proposed as an alternative solution and have shown great promise to reduce uncertainties associated with the parameterization of small-scale processes. Yet, those approaches still show some important mismatches that are often attributed to the stochasticity of the considered process. This stochasticity can be due to coarse temporal resolution, unresolved variables, or simply to the inherent chaotic nature of the process. To address these issues, we propose a new type of parameterization (closure), which is built using memory-based neural networks, to account for the non-instantaneous response of the closure and to enhance its stability and prediction accuracy. We apply the proposed memory-based parameterization, with differentiable solver, to the Lorenz ’96 model in the presence of a coarse temporal resolution and show its capacity to predict skillful forecasts over a long time horizon of the resolved variables compared to instantaneous parameterizations. This approach paves the way for the use of memory-based parameterizations for closure problems.

Mathematics↗

DOE OSTI2024

Cluster-dynamics-based parameterization for sulfuric acid–dimethylamine nucleation: comparison and selection through box and three-dimensional modeling

Clustering of gaseous sulfuric acid (SA) enhanced by dimethylamine (DMA) is a major mechanism for new particle formation (NPF) in polluted atmospheres. However, uncertainty remains regarding the SA–DMA nucleation parameterization that reasonably represents cluster dynamics and is applicable across various atmospheric conditions. This uncertainty hinders accurate three-dimensional (3-D) modeling of NPF and the subsequent assessment of its environmental and climatic impacts. Here we extensively compare different cluster-dynamics-based parameterizations for SA–DMA nucleation and identify the most reliable one through a combination of box model simulations, 3-D modeling, and in situ observations. Results show that the parameterization derived from Atmospheric Cluster Dynamic Code (ACDC) simulations, incorporating the latest theoretical insights (DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97X-D/6-311++G(3df,3pd) level of theory) and adequate representation of cluster dynamics, exhibits dependable performance in 3-D NPF simulation for both winter and summer conditions in Beijing and shows promise for application in diverse atmospheric conditions. Another ACDC-derived parameterization, replacing the level of theory with RI-CC2/aug-cc-pV(T+d)Z//M06-2X/6–311++G(3df,3pd), also performs well in NPF modeling at relatively low temperatures around 280 K but exhibits limitations at higher temperatures due to inappropriate representation of SA–DMA cluster thermodynamics. Additionally, a previously reported parameterization incorporating simplifications is applicable for simulating NPF in polluted atmospheres but tends to overestimate particle formation rates under conditions of elevated temperature (>∼300 K) and low-condensation sink (<∼3×10 -3 s −1 ). Our findings highlight the applicability of the new ACDC-derived parameterization, which couples the latest SA–DMA nucleation theory and holistic cluster dynamics, in 3-D NPF modeling. The ACDC-derived parameterization framework provides a valuable reference for developing parameterizations for other nucleation systems.

Shen, Jiewen↗

DOE OSTI2026

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗

NASA NTRS1992

Methods of testing parameterizations: Vertical ocean mixing

The ocean's velocity field is characterized by an exceptional variety of scales. While the small-scale oceanic turbulence responsible for the vertical mixing in the ocean is of scales a few centimeters and smaller, the oceanic general circulation is characterized by horizontal scales of thousands of kilometers. In oceanic general circulation models that are typically run today, the vertical structure of the ocean is represented by a few tens of discrete grid points. Such models cannot explicitly model the small-scale mixing processes, and must, therefore, find ways to parameterize them in terms of the larger-scale fields. Finding a parameterization that is both reliable and plausible to use in ocean models is not a simple task. Vertical mixing in the ocean is the combined result of many complex processes, and, in fact, mixing is one of the less known and less understood aspects of the oceanic circulation. In present models of the oceanic circulation, the many complex processes responsible for vertical mixing are often parameterized in an oversimplified manner. Yet, finding an adequate parameterization of vertical ocean mixing is crucial to the successful application of ocean models to climate studies. The results of general circulation models for quantities that are of particular interest to climate studies, such as the meridional heat flux carried by the ocean, are quite sensitive to the strength of the vertical mixing. We try to examine the difficulties in choosing an appropriate vertical mixing parameterization, and the methods that are available for validating different parameterizations by comparing model results to oceanographic data. First, some of the physical processes responsible for vertically mixing the ocean are briefly mentioned, and some possible approaches to the parameterization of these processes in oceanographic general circulation models are described in the following section. We then discuss the role of the vertical mixing in the physics of the large-scale ocean circulation, and examine methods of validating mixing parameterizations using large-scale ocean models.

Tziperman, Eli↗

NASA NTRS1994

Evaluation of a GCM cirrus parameterization using satellite observations

This study applies a simple yet effective methodology to validate a general circulation model parameterization of cirrus ice water path. The methodology combines large-scale dynamic and thermodynamic fields from operational analyses with prescribed occurrence of cirrus clouds from satellite observations to simulate a global distribution of ice water path. The predicted cloud properties are then compared with the corresponding satellite measurements of visible optical depth and infrared cloud emissivity to evaluate the reliability of the parameterization. This methodology enables the validation to focus strictly on the water loading side of the parameterization by eliminating uncertainties involved in predicting the occurrence of cirrus internally within the parameterization. Overall the parameterization performs remarkably well in capturing the observed spatial patterns of cirrus optical properties. Spatial correlations between the observed and the predicted optical depths are typically greater than 0.7 for the tropics and northern hemisphere midlatitudes. The good spatial agreement largely stems from the strong dependence of the ice water path upon the temperature of the environment in which the clouds form. Poorer correlations (r approximately 0.3) are noted over the southern hemisphere midlatitudes, suggesting that additional processes not accounted for by the parameterization may be important there. Quantitative evaluation of the parameterization is hindered by the present uncertainty in the size distribution of cirrus ice particles. Consequently, it is difficult to determine if discrepancies between the observed and the predicted optical properties are attributable to errors in the parameterized ice water path or to geographic variations in effective radii.

Soden, B. J.↗

DOE OSTI2015

Effects of pre-existing ice crystals on cirrus clouds and comparison between different ice nucleation parameterizations with the Community Atmosphere Model (CAM5)

In order to improve the treatment of ice nucleation in a more realistic manner in the Community Atmosphere Model version 5.3 (CAM5.3), the effects of pre-existing ice crystals on ice nucleation in cirrus clouds are considered. In addition, by considering the in-cloud variability in ice saturation ratio, homogeneous nucleation takes place spatially only in a portion of the cirrus cloud rather than in the whole area of the cirrus cloud. Compared to observations, the ice number concentrations and the probability distributions of ice number concentration are both improved with the updated treatment. The pre-existing ice crystals significantly reduce ice number concentrations in cirrus clouds, especially at mid- to high latitudes in the upper troposphere (by a factor of ~10). Furthermore, the contribution of heterogeneous ice nucleation to cirrus ice crystal number increases considerably. Besides the default ice nucleation parameterization of Liu and Penner (2005, hereafter LP) in CAM5.3, two other ice nucleation parameterizations of Barahona and Nenes (2009, hereafter BN) and Kärcher et al. (2006, hereafter KL) are implemented in CAM5.3 for the comparison. In-cloud ice crystal number concentration, percentage contribution from heterogeneous ice nucleation to total ice crystal number, and pre-existing ice effects simulated by the three ice nucleation parameterizations have similar patterns in the simulations with present-day aerosol emissions. However, the change (present-day minus pre-industrial times) in global annual mean column ice number concentration from the KL parameterization (3.24 × 10 6 m -2 ) is less than that from the LP (8.46 × 10 6 m -2 ) and BN (5.62 × 10 6 m -2 ) parameterizations. As a result, the experiment using the KL parameterization predicts a much smaller anthropogenic aerosol long-wave indirect forcing (0.24 W m -2 ) than that using the LP (0.46 W m −2 ) and BN (0.39 W m -2 ) parameterizations.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2014

A Flexible Parameterization for Shortwave Optical Properties of Ice Crystals

A parameterization is presented that provides extinction cross section sigma (sub e), single-scattering albedo omega, and asymmetry parameter (g) of ice crystals for any combination of volume, projected area, aspect ratio, and crystal distortion at any wavelength in the shortwave. Similar to previous parameterizations, the scheme makes use of geometric optics approximations and the observation that optical properties of complex, aggregated ice crystals can be well approximated by those of single hexagonal crystals with varying size, aspect ratio, and distortion levels. In the standard geometric optics implementation used here, sigma (sub e) is always twice the particle projected area. It is shown that omega is largely determined by the newly defined absorption size parameter and the particle aspect ratio. These dependences are parameterized using a combination of exponential, lognormal, and polynomial functions. The variation of (g) with aspect ratio and crystal distortion is parameterized for one reference wavelength using a combination of several polynomials. The dependences of g on refractive index and omega are investigated and factors are determined to scale the parameterized (g) to provide values appropriate for other wavelengths. The parameterization scheme consists of only 88 coefficients. The scheme is tested for a large variety of hexagonal crystals in several wavelength bands from 0.2 to 4 micron, revealing absolute differences with reference calculations of omega and (g) that are both generally below 0.015. Over a large variety of cloud conditions, the resulting root-mean-squared differences with reference calculations of cloud reflectance, transmittance, and absorptance are 1.4%, 1.1%, and 3.4%, respectively. Some practical applications of the parameterization in atmospheric models are highlighted.

parameterization↗

DOE OSTI2021

A Parameterization of Turbulent Dissipation and Pressure Damping Time Scales in Stably Stratified Inversions, and its Effects on Low Clouds in Global Simulations

Abstract It is difficult for coarse‐resolution global models of the atmosphere to accurately simulate the observed distribution of low clouds. In particular, it is difficult for moist turbulence closure models to simulate sufficiently bright near‐coastal stratocumulus (Sc) without simulating overly bright marine shallow cumuli (Cu). To parameterize bright Sc, a turbulence parameterization must damp the turbulent fluxes of heat and moisture above cloud top in order to prevent excessive entrainment of dry air into cloud top. To parameterize dim shallow Cu, the subgrid variances of temperature and moisture must remain large, in order to permit partial cloudiness. However, damping the fluxes but not the variances just above cloud top is difficult if a parameterization uses a single “master” time scale to damp both. In nature, the above‐cloud fluxes are damped by pressure fluctuations, whereas scalar variances are damped by a different process, namely, turbulent dissipation. In a stably stratified inversion above cloud, pressure damping is large but turbulent dissipation is small. To avoid this problem, a multitime scale parameterization for damping has been developed. The damping parameterization has been implemented in a global model and evaluated. The parameterization is capable of dimming shallow Cu while producing adequately bright Sc.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Elastic full-waveform inversion and parameterization analysis applied to walk-away vertical seismic profile data for unconventional (heavy oil) reservoir characterization

We report seismic full-waveform inversion (FWI) methods hold strong potential to recover multiple subsurface elastic properties for hydrocarbon reservoir characterization. Simultaneously updating multiple physical parameters introduces the problem of interparameter tradeoff, arising from the covariance between different physical parameters, which increases nonlinearity and uncertainty of multiparameter FWI. The coupling effects of different physical parameters are significantly influenced by model parameterization and acquisition arrangement. An appropriate choice of model parameterization is critical to successful field data applications of multiparameter FWI. The objective of this paper is to examine the performance of various model parameterizations in isotropic-elastic FWI with walk-away vertical seismic profile (W-VSP) dataset for unconventional heavy oil reservoir characterization. Six model parameterizations are considered: velocity-density (α, β and ρ'), modulus-density (κ, μ and ρ), Lamé-density (λ, μ' and ρ'''), impedance-density (IP, IS and ρ''), velocity-impedance-I (α', β' and I' P ), and velocity-impedance-II (α'', β'' and I'S). We begin analyzing the interparameter tradeoff by making use of scattering radiation patterns, which is a common strategy for qualitative parameter resolution analysis. In this paper, we discuss the advantages and limitations of the scattering radiation patterns and recommend that interparameter tradeoffs be evaluated using interparameter contamination kernels, which provide quantitative, second-order measurements of the interparameter contaminations and can be constructed efficiently with an adjoint-state approach. Synthetic W-VSP isotropic-elastic FWI experiments in the time domain verify our conclusions about interparameter tradeoffs for various model parameterizations. Density profiles are most strongly influenced by the interparameter contaminations; depending on model parameterization, the inverted density profile can be over-estimated, under-estimated or spatially distorted. Among the six cases, only the velocity-density parameterization provides stable and informative density features not included in the starting model. Field data applications of multicomponent W-VSP isotropic-elastic FWI in the time domain were also carried out. Finally, the heavy oil reservoir target zone, characterized by low α-to-β ratios and low Poisson’s ratios, can be identified clearly with the inverted isotropic-elastic parameters.

58 GEOSCIENCES↗

DOE OSTI2019

A simplified parameterization of isoprene-epoxydiol-derived secondary organic aerosol (IEPOX-SOA) for global chemistry and climate models: a case study with GEOS-Chem v11-02-rc

Secondary organic aerosol derived from isoprene epoxydiols (IEPOX-SOA) is thought to contribute the dominant fraction of total isoprene SOA, but the current volatility-based lumped SOA parameterizations are not appropriate to represent the reactive uptake of IEPOX onto acidified aerosols. A full explicit modeling of this chemistry is however computationally expensive owing to the many species and reactions tracked, which makes it difficult to include it in chemistry–climate models for long-term studies. Here we present three simplified parameterizations (version 1.0) for IEPOX-SOA simulation, based on an approximate analytical/fitting solution of the IEPOX-SOA yield and formation timescale. The yield and timescale can then be directly calculated using the global model fields of oxidants, NO, aerosol pH and other key properties, and dry deposition rates. The advantage of the proposed parameterizations is that they do not require the simulation of the intermediates while retaining the key physicochemical dependencies. We have implemented the new parameterizations into the GEOS-Chem v11-02-rc chemical transport model, which has two empirical treatments for isoprene SOA (the volatility-basis-set, VBS, approach and a fixed 3% yield parameterization), and compared all of them to the case with detailed fully explicit chemistry. The best parameterization (PAR3) captures the global tropospheric burden of IEPOX-SOA and its spatiotemporal distribution ( R 2 =0.94 ) vs. those simulated by the full chemistry, while being more computationally efficient ( ~5 times faster), and accurately captures the response to changes in NO x and SO 2 emissions. On the other hand, the constant 3 % yield that is now the default in GEOS-Chem deviates strongly ( R 2 =0.66 ), as does the VBS ( R 2 =0.47 , 49% underestimation), with neither parameterization capturing the response to emission changes. With the advent of new mass spectrometry instrumentation, many detailed SOA mechanisms are being developed, which will challenge global and especially climate models with their computational cost. The methods developed in this study can be applied to other SOA pathways, which can allow including accurate SOA simulations in climate and global modeling studies in the future.

58 GEOSCIENCES↗

DOE OSTI2020

A new parameterization of ice heterogeneous nucleation coupled to aerosol chemistry in WRF-Chem model version 3.5.1: evaluation through ISDAC measurements

In the Arctic, during polar night and early spring, ice clouds are separated into two leading types of ice clouds (TICs): (1) TIC1 clouds characterized by a large concentration of very small crystals and TIC2 clouds characterized by a low concentration of large ice crystals. Using a suitable parameterization of heterogeneous ice nucleation is essential for properly representing ice clouds in meteorological and climate models and subsequently understanding their interactions with aerosols and radiation. Here, we describe a new parameterization for ice crystal formation by heterogeneous nucleation in water-subsaturated conditions coupled to aerosol chemistry in the Weather Research and Forecasting model coupled with chemistry (WRF-Chem). The parameterization is implemented in the Milbrandt and Yau (2005a, b) two-moment cloud microphysics scheme, and we assess how the WRF-Chem model responds to the run-time interaction between chemistry and the new parameterization. Well-documented reference cases provided us with in situ data from the spring 2008 Indirect and Semi-Direct Aerosol Campaign (ISDAC) over Alaska. Our analysis reveals that the new parameterization clearly improves the representation of the ice water content (IWC) in polluted or unpolluted air masses and shows the poor performance of the reference parameterization in representing ice clouds with low IWC. The new parameterization is able to represent TIC1 and TIC2 microphysical characteristics at the top of the clouds, where heterogenous ice nucleation is most likely occurring, even with the known bias of simulated aerosols by WRF-Chem over the Arctic.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model

Most global climate models parameterize separate cloud types using separate parameterizations. This approach has several disadvantages, including obscure interactions between parameterizations and inaccurate triggering of cumulus parameterizations. Alternatively, a unified cloud parameterization uses one equation set to represent all cloud types. Such cloud types include stratiform liquid and ice cloud, shallow convective cloud, and deep convective cloud. Vital to the success of a unified parameterization is a general interface between clouds and microphysics. One such interface involves drawing Monte Carlo samples of subgrid variability of temperature, water vapor, cloud liquid, and cloud ice, and feeding the sample points into a microphysics scheme. This study evaluates a unified cloud parameterization and a Monte Carlo microphysics interface that has been implemented in the Community Atmosphere Model (CAM) version 5.3. Model computational expense is estimated, and sensitivity to the number of subcolumns is investigated. Results describing the mean climate and tropical variability from global simulations are presented. The new model shows a degradation in precipitation skill but improvements in shortwave cloud forcing, liquid water path, long-wave cloud forcing, precipitable water, and tropical wave simulation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model

Most global climate models parameterize separate cloud types using separate parameterizations. This approach has several disadvantages, including obscure interactions between parameterizations and inaccurate triggering of cumulus parameterizations. Alternatively, a unified cloud parameterization uses one equation set to represent all cloud types. Such cloud types include stratiform liquid and ice cloud, shallow convective cloud, and deep convective cloud. Vital to the success of a unified parameterization is a general interface between clouds and microphysics. One such interface involves drawing Monte Carlo samples of subgrid variability of temperature, water vapor, cloud liquid, and cloud ice, and feeding the sample points into a microphysics scheme.This study evaluates a unified cloud parameterization and a Monte Carlo microphysics interface that has been implemented in the Community Atmosphere Model (CAM) version 5.3. Results describing the mean climate and tropical variability from global simulations are presented. The new model shows a degradation in precipitation skill but improvements in short-wave cloud forcing, liquid water path, long-wave cloud forcing, precipitable water, and tropical wave simulation. Also presented are estimations of computational expense and investigation of sensitivity to number of subcolumns.

58 GEOSCIENCES↗

NASA NTRS1985

Preliminary GCM Results with a New Radiation Parameterization

A new parameterization of solar and terrestrial radiation was developed and tested. The solar radiation parameterization is based on that of Lacis and Hansen (1974), but with zenith-angle-dependent surface albedoes, and revised treatments of cloudiness. The terrestrial radiation parameterization is based on the work of Chou (1984) for water vapor, Chou and Peng (1983) for carbon dioxide, and Rogers (1968) for ozone, with a new parameterization of the effects of clouds. Results obtained were compared with the new parameterizations to those obtained with the earlier parameterization described by Schlessinger (1976). Several dramatic improvements came to light. For the most part these are related to the fact that the new terrestrial radiation paramerization includes the effects of the water vapor continuum, while the earlier parameterization does not. In the moist tropical planetary boundary layer (PBL) continuum emission leads to much stronger cooling of the PBL over the oceans. Over land, however, the cooling of the PBL is significantly reduced. The latter, somewhat paradoxical result is due to the strong diurnal cycle of the continental PBL. At night the shallow continental PBL is overlain by a moist layer created by mixing during the previous afternoon. The moist upper layer acts as a radiative blanket, reducing the time-averaged radiative cooling of the continental PBL.

Randall, D.↗

NASA NTRS1994

The response of the SSM/I to the marine environment. Part 2: A parameterization of the effect of the sea surface slope distribution on emission and reflection

Based on a geometric optics model and the assumption of an isotropic Gaussian surface slope distribution, the component of ocean surface microwave emissivity variation due to large-scale surface roughness is parameterized for the frequencies and approximate viewing angle of the Special Sensor Microwave/Imager. Independent geophysical variables in the parameterization are the effective (microwave frequency dependent) slope variance and the sea surface temperature. Using the same physical model, the change in the effective zenith angle of reflected sky radiation arising from large-scale roughness is also parameterized. Independent geophysical variables in this parameterization are the effective slope variance and the atmospheric optical depth at the frequency in question. Both of the above model-based parameterizations are intended for use in conjunction with empirical parameterizations relating effective slope variance and foam coverage to near-surface wind speed. These empirical parameterizations are the subject of a separate paper.

Petty, Grant W.↗