Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Forecasting Remote Atmospheric Responses to Decadal Kuroshio Stability Transitions

Based on observational estimates and global ocean eddy-resolving coupled retrospective initialized predictions, we show that Kuroshio Extension variability affects rainfall variability along the west coast of North America. We show that the teleconnection between the current undulations and downstream rainfall can lead to improved subseasonal to seasonal predictions of precipitation over California, and we demonstrate that capturing these teleconnections requires coupled systems with sufficient ocean resolution (i.e., eddy-resolving), especially over time scales longer than one season. The improved forecast skill is diagnosed in terms of 35 years of retrospective initialized ensemble forecasts with an ocean eddy-resolving and an ocean eddy-parameterized coupled model. Not only does the ocean eddy-resolving model show sensitivity to Kuroshio Extension variability in terms of western North America precipitation, but the ocean eddy-resolving forecasts also show improved forecast skill compared to the ocean eddy-parameterized model. Additionally, the ocean eddy-parameterized coupled model shows no sensitivity to Kuroshio Extension variability. We also find near-decadal variability associated with a progression of a lower-tropospheric height dipole around the North Pacific and how these height anomalies lead to wind-driven Rossby waves that affect the eddy activity in the Kuroshio Extension with a time lag on the order of four years. This decadal-scale variability (~10 years) opens the possibility of multiyear predictability of western North American rainfall.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗
DOE OSTI2015

Prognostic residual mean flow in an ocean general circulation model and its relation to prognostic Eulerian mean flow

Recent work has shown that taking the thickness-weighted average (TWA) of the Boussinesq equations in buoyancy coordinates results in exact equations governing the prognostic residual mean flow where eddy–mean flow interactions appear in the horizontal momentum equations as the divergence of the Eliassen–Palm flux tensor (EPFT). It has been proposed that, given the mathematical tractability of the TWA equations, the physical interpretation of the EPFT, and its relation to potential vorticity fluxes, the TWA is an appropriate framework for modeling ocean circulation with parameterized eddies. The authors test the feasibility of this proposition and investigate the connections between the TWA framework and the conventional framework used in models, where Eulerian mean flow prognostic variables are solved for. Using the TWA framework as a starting point, this study explores the well-known connections between vertical transfer of horizontal momentum by eddy form drag and eddy overturning by the bolus velocity, used by Greatbatch and Lamb and Gent and McWilliams to parameterize eddies. After implementing the TWA framework in an ocean general circulation model, we verify our analysis by comparing the flows in an idealized Southern Ocean configuration simulated using the TWA and conventional frameworks with the same mesoscale eddy parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Controls on Turbulent Mixing in a Strongly Stratified and Sheared Tidal River Plume

Considerable effort has been made to parameterize turbulent kinetic energy (TKE) dissipation rate ..epsilon.. and mixing in buoyant plumes and stratified shear flows. Here, a parameterization based on Kunze et al. is examined, which estimates ..epsilon.. as the amount of energy contained in an unstable shear layer (Ri < Ric) that must be dissipated to increase the Richardson number Ri = N2/S2 to a critical value Ric within a turbulent decay time scale. Observations from the tidal Columbia River plume are used to quantitatively assess the relevant parameters controlling ..epsilon.. over a range of tidal and river discharge forcings. Observed ..epsilon.. is found to be characterized by Kunze et al.'s form within a factor of 2, while exhibiting slightly decreased skill near Ri = Ric. Observed dissipation rates are compared to estimates from a constant interfacial drag formulation that neglects the direct effects of stratification. This is found to be appropriate in energetic regimes when the bulk-averaged Richardson number Rib is less than Ric/4. However, when Rib > Ric/4, the effects of stratification must be included. Similarly, ..epsilon.. scaled by the bulk velocity and density differences over the plume displays a clear dependence on Rib, decreasing as Rib approaches Ric. The Kunze et al. ..epsilon.. parameterization is modified to form an expression for the nondimensional dissipation rate that is solely a function of Rib, displaying good agreement with the observations. It is suggested that this formulation is broadly applicable for unstable to marginally unstable stratified shear flows.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

Improving Wind Predictions in the Marine Atmospheric Boundary Layer through Parameter Estimation in a Single-Column Model

A current barrier to greater deployment of offshore wind turbines is the poor quality of numerical weather prediction model wind and turbulence forecasts over open ocean. The bulk of development for atmospheric boundary layer (ABL) parameterization schemes has focused on land, partly due to a scarcity of observations over ocean. The 100-m FINO1 tower in the North Sea is one of the few sources worldwide of atmospheric profile observations from the sea surface to turbine hub height. These observations are crucial to developing a better understanding and modeling of physical processes in the marine ABL. In this study, we use the WRF single column model (SCM), coupled with an ensemble Kalman filter from the Data Assimilation Research Testbed (DART), to create 100-member ensembles at the FINO1 location. The goal of this study is to determine the extent to which model parameter estimation can improve offshore wind forecasts.

17 WIND ENERGY↗
DOE OSTI2025

Data Files for Runoff Evaluation in an Earth System Land Model for Permafrost Regions

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce.This data set includes all files that were produced and applied in the paper Runoff Evaluation in an Earth System Land Model for Permafrost Regions [Xiang et al. in review]. The paper is in review as of July 1 2025 in Geoscientific Model Development (GMD). In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-rich integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM’s parameterized representation of total runoff. This dataset contains 2 figure image files (*.png, *jpg) that describe the study site and methods, as well as folders (Figure*.zip) that contain the associated data files (*.csv, *.dat) and python code notebooks (*.ipynb) for figures 3-7 in the paper. Jupyter notebook (*.ipynb) files that produce the figure files using the associated data files will run within a python environment configured with Jupyter Lab or Notebook packages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗
DOE OSTI2018

Characterization of Grain‐Size Distribution, Thermal Conductivity, and Gas Diffusivity in Variably Saturated Binary Sand Mixtures

Core Ideas Parameterized models are relevant to applications where different porous media mixtures are used. Models include a grain‐size distribution function to describe bimodal behavior for binary mixtures. Improved model describes observed thermal conductivity–saturation relations for binary mixtures. Combined Buckingham–Penman model used to describe observed gas diffusivity–air content relations. This work is relevant for proper simulation of mixed porous media. Characterization of differently textured porous materials, as well as different volumetric porous media mixtures, in relation to mass and heat transport is vital for many engineering and research applications. Functional relations describing physical (e.g., grain‐size distribution, total porosity), thermal, and gas diffusion properties of porous media and mixtures are necessary to optimize the design of porous systems that involve heat and gas transport processes. However, only a limited number of studies provide characterization of soil physical, thermal, and gas diffusion properties and the functional relationships of these properties under varying soil water contents, especially for soil mixtures, complicating optimization efforts. To better understand how mixing controls the physical, thermal, and gas diffusion properties of porous media, a set of laboratory experiments was performed using five volumetric mixtures of coarse‐ and fine‐grained sand particles. For each mixture, the grain‐size distribution (GSD), thermal conductivity, and gas diffusivity were obtained and parameterized using existing and suggested parametric models. Results show that the extended, two‐region Rosin–Rammler particle‐size distribution model proposed in this study could successfully describe the bimodal behavior of the GSD of binary mixtures. Further, the modified Côté and Konrad thermal conductivity model adequately described the thermal conductivity–water saturation relations observed in different mixtures. The proposed simple soil‐gas diffusivity descriptive model parameterized the upper limit, average, and lower limit behavior in gas diffusivity–air content relations in apparently texture‐invariant gas diffusivity data. Results further show a close analogy between gas diffusivity and thermal conductivity and their variation with saturation across different binary mixtures. Overall, the results of the study provide useful numerical insight into the physical, thermal, and gas transport characteristics of binary mixtures, with wide implications for future engineering and research applications that involve multicomponent porous systems.

Deepagoda, T. K. K. Chamindu↗
DOE OSTI2015

Final Report for Collaborative Project: Sensitivity of Atmospheric Parametric Formulations to Regional Mesh Refinement in Global Climate Simulations Using CESM-HOMME

In this project we analyze climate simulations using the Community Earth System Model (CESM) in order to determine the modeled response and sensitivity to horizontal resolution. Simple aqua-planet configurations were used to provide a clean comparison of the response to resolution in CESM. This enables us to easily examine all aspects of the model sensitivity to resolution including mean quantities, variability and physical parameterization tendencies: the chief reflection of resolution sensitivity. An extension to the global resolution sensitivity study is the examination of regional grid refinement where resolution changes are prescribed in a single global simulation. We examine the relevance of the global resolution sensitivity results as applied to these regional refinement simulations. In particular we examine how variations in the grid resolution, centered on different parts of the globe, lead to differences in the parameterized response and the potential to generate residual circulations as a result. Given the potential to generate this resolution sensitivity we examine simple modifications to the parameterized physics that are able to moderate any residual circulations. Finally, we transfer the framework to the standard AMIP configuration to examine the resolution sensitivity in the presence of compounding effects such as land-sea distributions, orography and seasonal variation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Characterization of 3D Cirrus Cloud and Radiation Fields Using ARS/AIRS/MODIS data and its Application to Climate Model

During the report period, we have made the following research accomplishments. First, we performed analysis for a number of MODIS scenes comprising of heavy dust events and ice clouds, covering regions of frequent dust outbreaks in East Asia, Middle East, and West Africa, as well as areas associated with long-range dust transports over the Equatorial Tropical Atlantic Ocean. These scenes contain both dust/aerosols and clouds. We collected suitable aerosol/ice-cloud data, correlated ice cloud and aerosol parameters by means of statistical analysis, and interpreted resulting correlation trends based on the physical principles governing cloud microphysics. Aerosol and cloud optical depths and cloud effective particle size inferred from MODIS for selected domains were analyzed from which the parameters including dust aerosol number concentration, ice cloud water path, and ice particle number concentration were subsequently derived. We illustrated that the Twomey (solar albedo) effect can be statistically quantified based on the slope of best-fit straight lines in the correlation study. Analysis of aerosol and cloud retrieval products revealed that for all cases, the region with a larger dust aerosol optical depth is always characterized by a smaller cloud particle size, consistent with the Twomey hypothesis for aerosol-cloud interactions. Second, we developed mean correlation curves with uncertainties associated with small ice-crystal concentration observations for the mean effective ice crystal size (De) and ice water content (IWC) by dividing the atmosphere into three characteristic regions: Tropics cirrus, Midlatitude cirrus, including a temperature classification to improve correlation, and Arctic ice clouds. We illustrated that De has a high correlation with IWC based on theoretical consideration and analysis of thousands of observed ice crystal data obtained from a number of ARM-DOE field campaigns and other experiments. The correlation has the form: ln(De) = a + b ln(IWC) + c ((ln(IWC))2, where a, b, and c are fitting coefficients and are functions of three regions. We demonstrated that this correlation can be effectively incorporated in GCMs and climate models that predict IWC - a significant advance in ice microphysics parameterization for interactive cloud-radiation analysis and feedback. Substantial July mean differences are shown in the OLR (W/m2) and precipitation (mm/day) patterns between UCLA GCM simulations based on Des determined from the De-IWC correlations and the control run using a fixed ice crystal size. Third, in order to improve the computation of spectral radiative transfer processes in the WRF model, we developed a consistent and efficient radiation scheme that can better resolve the spectral bands, determine the cloud optical properties, and provide more reliable and accurate radiative heating fields. In the newly developed radiation module, we have implemented in WRF a modified and improved version referred to as the Fu-Liou-Gu scheme, which includes a combination of delta-four-stream and delta-two-stream approximations for solar and IR flux calculations, respectively. This combination has been proven to be computationally efficient and at the same time to produce a high degree of accuracy. The incorporation of nongray gaseous absorption in multiple scattering atmospheres was based on the correlated k-distribution method. The solar and IR spectra are divided into 6 and 12 bands, respectively, according to the location of absorption bands of H2O, CO2, O3, CH4, N2O, and CFCs. We further included absorption by the water vapor continuum and a number of minor absorbers in the solar spectrum leading to an additional absorption of solar flux in a clear atmosphere on the order of 1-3 W/m2. Additionally, we incorporated the ice microphysics parameterization that includes an interactive mean effective ice crystal size in association with radiation parameterizations. The Fu-Liou-Gu scheme is an ideal tool for the simulation of radiative transfer and ice microphysics within the domain of WRF. It is particularly useful for studying direct and indirect aerosol radiative effects associated with ice cloud formation. The newly implemented radiation module has been demonstrated to work well in WRF and can be effectively used for studies related to cirrus cloud formation and evolution as well as aerosol-cloud-radiation interactions. With the newly implemented radiation scheme, the simulations of cloud cover and ice water path have been improved for cirrus clouds, with a more consistent comparison with the corresponding MODIS observations, especially for optically thin cirrus with an improvement of about 20% in the simulated mean ice water path.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Site Scientist for the North Slope of Alaska Site

Under this grant our team contributed scientific support to the Department of Energy Atmospheric Radiation Program’s (DOE-ARM) Infrastructure team to maintain high quality research data at the DOE-ARM North Slope of Alaska with special emphasis on the radars. Under our guidance two major field campaigns focusing on mixed-phase Arctic clouds were conducted that greatly increased the community’s understanding of the many processes working together to control the evolution of single-layer cloud mixed-phase clouds. A series of modeling and observational studies revealed that the longevity of the radiatively important liquid phase is strongly dependent on how the ice phase develops in mixed-phase clouds. A new ice microphysics parameterization was developed to capture better the natural evolution of ice particle growth in evolving environments. An ice particle scattering database was developed for all ARM radar frequencies. This database was used in a radar simulator (Doppler spectrum and polarimetric variables) to aid in the interpretation of the advanced ARM radars. At the conclusion of this project our team was poised to develop a complete radar simulator consistent with the new microphysical parameterization, taking advantage of parameterization’s advanced characterization of the ice shape and ice density.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Evaluation of Simulated Marine Aerosol Production Using the WaveWatchIII Prognostic Wave Model Coupled to the Community Atmosphere Model within the Community Earth System Model

Primary marine aerosol (PMA) is emitted into the atmosphere via breaking wind waves on the ocean surface. Most parameterizations of PMA emissions use 10-meter wind speed as a proxy for wave action. This investigation coupled the 3 rd generation prognostic WAVEWATCH-III wind-wave model within a coupled Earth system model (ESM) to drive PMA production using wave energy dissipation rate – analogous to whitecapping – in place of 10-meter wind speed. The wind speed parameterization did not capture basin-scale variability in relations between wind and wave fields. Overall, the wave parameterization did not improve comparison between simulated versus measured AOD or Na + , thus highlighting large remaining uncertainties in model physics. Results confirm the efficacy of prognostic wind-wave models for air-sea exchange studies coupled with laboratory- and field-based characterizations of the primary physical drivers of PMA production. No discernible correlations were evident between simulated PMA fields and observed chlorophyll or sea surface temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Multiscale Methods for Accurate, Efficient, and Scale-Aware Models

The goal of UWM’s portion of the Multiscale project was to develop a unified cloud parameterization that could simulate all cloud types --- including stratocumulus, shallow cumulus, and deep cumulus --- using the single equation set implemented in CLUBB. An advantage of a unified parameterization methodology is that it avoids the difficult task of interfacing different cloud parameterizations for different cloud types. To interface CLUBB’s clouds to the microphysics, a Monte Carlo interface, SILHS, was further developed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Major improvements on the longwave radiative interactions between surface and clouds in the Polar Regions in atmospheric global circulation model (GCM). Final Project Report

1. Summary of Work and Accomplishments in Year 1 (01/15/2015-01/14/2016). We developed new ice cloud parameterizations based on state-of-the-art MODIS Collection 6 (MC6) ice particle model. Features of the model include the consideration of surface roughness of ice crystals and spectral consistency in cloud property retrievals. Using more than 14,000 particle size distributions (PSDs) from in-situ measurements and the latest ice particle single-scattering property library, our parameterizations provide band-averaged bulk-scattering properties. According to statistics of fittings, the parameterization schemes represent the band-averaged bulk-scattering properties very well. Compared to the ice particle model developed by Fu (1998), flux simulations with the MC6 particle model result in larger upward flux at the top of the atmosphere (TOA) and smaller at the surface for typical cirrus clouds. In the assessment of the effect of longwave light scattering, our results also support the conclusions of previous studies. The error in flux calculations reaches about 10 W/m 2 when a cirrus cloud is thin and its particle size is small. 2. Summary of Work and Accomplishments in Year 2 (01/15/2016-01/14/2017) Clouds are a major modulator of the global radiation budgets. However, representations of clouds in general circulation models (GCMs) neglect many important processes, particularly those related to ice clouds. One major uncertainty source is the radiative transfer scheme in GCMs. To reduce the computational burden in most GCMs, absorption is the only radiative transfer process considered in longwave spectral bands. Recently, a state-of-the-art MC6 ice model has been available, which provides spectral consistency in ice cloud retrievals. With an advanced light scattering library for ice particles, we performed sensitivity studies of the potential impacts on the climate including light scattering for MC6 model and other 11 ice cloud models. Specifically, we modified RRTMG_LW (the GCM version of Longwave Rapid Radiative Transfer Model), a radiative transfer module widely used in GCMs and numerical prediction models, to identify flux simulation differences with and without considering light scattering processes. The simulations with light scattering come from a rigorous radiative solver, DISORT (Discrete Ordinates Radiative Transfer Program for a Multi-Layered Plane-Parallel Medium), are used as benchmarks. The results show that the weighted annual mean biases of RRTMG_LW for the upward flux at the top of the atmosphere, the downward flux at the surface, and the net flux into the atmosphere are about 0.8 ± 0.3, -0.1 ± 0.05, and -0.7 ± 0.3 W m-2, respectively. According to the different the ice crystal shapes, the weighted annual mean bias for heating rate is about -0.006 ± 0.02 to 0.04 K/day. 3. Summary of Work and Accomplishments in Year 3 (01/15/2017-01/14/2018) Since absorption dominates optical properties of clouds in the longwave (LW) spectrum, most of general circulation models (GCMs) only take absorption properties of clouds into account in the radiative transfer modules in order to reduce computational costs. To quantify the biases of excluding LW scattering when clouds exist, we simulated fluxes and heating rates by using satellite data from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO), CloudSat, Clouds and the Earth’s Radiant Energy System (CERES) and Moderate Resolution Imaging Spectrometer (MODIS) merged products (CCCM) in 2010. In annual global average, neglecting LW scattering causes overestimation by 2.6 W/m 2 for the top-of-atmosphere (TOA) upward fluxes and underestimation by 1.2 W/m 2 for surface downward fluxes. In particular, regional extreme biases are approximately 12 W/m 2 for TOA upward fluxes and are approximately -3.6 W/m 2 for surface downward fluxes, or are approximately 5% of the global averaged outgoing longwave radiation (OLR) and 1% of the global averaged surface downward fluxes, respectively. In terms of heating rates, on average, LW scattering heats the whole atmosphere column by 0.0045 K/day, heats the cloud layers by 0.0420 K/day, heats the surface by 0.028 K/day, and cools the tropopause by 0.018 K/day.

58 GEOSCIENCES↗
DOE OSTI2018

2014-2018 Summary of NCAR activities in the Earth Energy Exascale System Model (E3SM) project

NCAR participation in E3SM activities as part of this funded collaborative project has been successful in leveraging the existing atmosphere model development experience from CESM scientists and applying it to initial v1 atmosphere model (EAM) development and analysis. NCAR partners collectively had many decades of experience in parameterization development, running and tuning coupled models and building the required software infrastructure to support ongoing model developments. Furthermore, collaboration between E3SM and CESM developers was of mutual benefit to both modeling endeavors. In the earlier periods of the project, participation in E3SM both complemented core critical path DOE activities and advanced individual parameterization development; including orographic drag, convection parameterization, convective gustiness and cloud and microphysical processes. Significant progress has been made implementing model capabilities for separating out the physics and dynamics grids, and providing the infrastructure for including orographic elevation class communication between the land and atmosphere models. Our activities working with the E3SM workflow and software engineering teams have established efficient protocols and infrastructure for the sharing of CESM and E3SM model codes of interest to both modeling groups, a large component of which is the mutual communication of common code bugs. As well as these specific project tasks all NCAR project members have been regular participants in atmosphere and coupled group telecons, as well as PI meetings. We have been fully engaged with and contributed to decisions regarding v1 atmosphere model developments, coupled model simulations, model tunings, and software engineering infrastructure. NCAR’s contribution to more recent E3SM development activities has centered on analysis of v1 release and pre-release simulations in order to address outstanding systematic biases and communicate climate improvements. Of particularly note is the work in understanding clouds and cloud feedbacks, shortcomings in ENSO and sub-seasonal variability; in particular the Madden Julian Oscillation (MJO). Contributions have been made to a number of papers that are currently published or in preparation detailing model performance in CMIP6-protocal experiments. They outline the simulations of mean climate, clouds and convection processes, convective gustiness, atmospheric blocking coupled climates. The contributions of NCAR’s funded activities to the E3SM atmosphere in the initial years of the project have had a positive impact on many aspects of the model simulations described in these papers.

58 GEOSCIENCES↗
DOE OSTI2019

Expanding the computational frontier of multi-scale atmospheric simulation to advance understanding of low cloud / climate feedbacks: Final Technical Report

Marine boundary layer clouds are a challenge for climate models. They cover much of the oceans, but they are driven by small-scale turbulent eddies only a few hundred meters across, which in turn respond to the cloud formation. Global climate models have a grid spacing that is far too large to simulate such fluid motions, so they parameterize these cloud formation processes. Different climate modeling groups have designed different parameterizations in which the clouds turn out to be differently sensitive to a warming climate. This drives uncertainties in our best guess at the sensitivity of global warming to greenhouse gas increases. If these cloud-forming eddies could be directly simulated using the well-known equations of fluid motion, they would no longer need to be parameterized, removing a major source of climate modeling uncertainty. In this project, we overcame software engineering challenges to successfully implemented ‘ultraparameterization’ (UP), the first global model that does this, and we tested how well it works. We simulated five-year periods with present-day temperatures and with a warmer climate, and we investigated how the UP-simulated clouds responded to climate – the ‘cloud feedback’ problem. We found little response of clouds at all latitudes to the imposed climate change, which is within the range of predictions of conventional global climate models. UP is a variation on superparameterization, in which small cloud-resolving models (CRMs) are embedded in each column of the global model. In UP, the CRM grid is fine enough (250 m horizontal × 20 m vertical) to explicitly capture boundary-layer turbulent eddies and associated clouds. Because only one small columnar patch is simulated within each climate model grid cell, this is a million-fold more efficient than simulating the entire globe on this same CRM grid., but achieves much of the same effect for the cloud properties. It doesn’t work perfectly. For instance, like conventional climate models, UP simulates too little subtropical stratocumulus cloud, a bias that we are continuing to work to reduce. However, because it directly simulates the turbulent cloud-forming processes, UP is inherently more plausible for simulating how clouds will change in a perturbed climate. In future, we hope to apply UP to another key climate modeling issue: cloud-aerosol interaction and the effect of human-produced aerosols on the climate change we have already experienced and that which is likely to come.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Evaluation of NCAR CAM5 simulated marine boundary layer cloud properties using a combination of satellite and surface observations

Marine boundary layer (MBL) clouds are a strong modulator of Earth’s radiative energy budget. The feedback of MBL clouds in a changing climate remains one of the largest uncertainties in our understanding of climate sensitivity, and recent studies have also suggested that radiative feedbacks related to MBL clouds may play an important role in low-frequency climate variabilities such as the Atlantic multidecadal oscillation. Despite their importance in the climate system, simulating MBL clouds in Earth system models (ESMs) has proved to be extremely challenging. The overarching objective of this proposal is to evaluate the MBL cloud simulation skills of the CAM5 with different cloud parameterization schemes through systematic comparison with both ground- and satellite based cloud observations. In particular, we seek to understand the advantages and limitations of the newly developed subgrid cloud parameterization scheme CLUBB (Cloud Layers Unified By Binormals) for simulation MBL clouds through comparisons of CAM5-CLUBB simulations with ARM (CAP-MBL campaign data) and satellite observations, as well as the “conventional” CAM5 and a CAM5 with the so-called super-parameterization (i.e., CAM5-MMF). In the rest of this report, we will briefly highlight our major research achievements.

58 GEOSCIENCES↗
DOE OSTI2017

Gradient-based Optimization for Regression in the Functional Tensor-Train Format

We consider the task of low-multilinear-rank functional regression, i.e., learning a low-rank parametric representation of functions from scattered real-valued data. Our first contribution is the development and analysis of gradient-based optimization procedures, including stochastic gradient descent and quasi-Newton methods, for learning the parameters of a functional tensor-train (FT) that yields improved accuracy over standard alternating least squares methods. The functional tensor-train uses the tensor-train (TT) representation of low-rank arrays as an ansatz for a class of low-multilinear-rank functions. The FT is represented by a set of matrix-valued functions that contain a set of univariate functions, and the regression task is to learn the parameters of these univariate functions. Our second contribution demonstrates that using nonlinearly parameterized univariate functions, e.g., symmetric kernels with moving centers, within each core can outperform the standard approach of using a linear expansion of basis functions. Our final contributions are new rank adaptation and group-sparsity regularization procedures to minimize overfitting. We use several benchmark problems to demonstrate at least an order of magnitude lower accuracy with gradient-based optimization methods than ALS in the low-sample number regime. We also demonstrate an order of magnitude reduction in accuracy on a test problem resulting from nonlinear parameterizations over linear parameterizations. Finally we compare regression performance with 22 other nonparametric and parametric regression methods on 10 real-world data sets. We achieve top-five accuracy for seven of the data sets and best accuracy for two of the data sets. These rankings are the best amongst parametric models and competetive with the best non-parametric methods.

97 MATHEMATICS AND COMPUTING↗