Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

WE-A-17A-04: Development of An Ultra-Fast Monte Carlo Dose Engine for High Dose Rate Brachytherapy

Purpose: Current clinical brachytherapy dose calculations are based on AAPM TG43 guidelines, which approximate the patient geometry as a large water phantom. This ignores heterogeneities and tends to overestimate skin dose. Although Monte Carlo (MC) dose calculations have been recognized as the most accurate method, its associated long computational time is a major bottleneck for routine clinical applications. This work aims to develop a GPU-based ultra-fast MC dose engine (gBMC) for HDR brachytherapy to provide clinical users with accurate sub-minute dose calculations. Methods: Standard photon transport with discrete events including: Compton scattering, Rayleigh scattering and Photoelectric effect, was implemented. Secondary electrons were transported under the continuous slowing down approximation. To reduce the GPU thread divergence, photons and electrons were separately transported. Transport of photons was grouped according to energy. The source model in gBMC can be either a phase-space file generated using Geant4 or a parameterized source model. This dose engine was validated against TG43 in a water phantom and against Geant4 calculations in heterogeneous patient geometries. Results: A phase space file was generated for the Varian VS2000 Ir-192 source. In a water phantom, the calculated radial dose function was within 0.6% of the TG43 calculations for radial distances from 1 cm to 20 cm. The anisotropy functions were within 1% for radial distances from 1 cm to 20 cm except for polar angles larger than 173°. Local point-dose differences were within 2%. In a Mammosite breast cancer case with 22 dwell locations, gBMC and Geant4 isodose lines compared well. The computation time was about 28 seconds using the phase-space file source and 20 seconds using the parameterized source to simulate 1 billion particles, yielding less than 1% statistical uncertainty. Conclusion: The gBMC dose engine makes it possible to use fast and accurate MC dose calculations for clinical work.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Importance of different parameterization changes for the updated dust cycle modeling in the Community Atmosphere Model (version 6.1)

Abstract. The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Recasting a model atomistic glassformer as a system of icosahedra

We consider a binary Lennard-Jones glassformer whose super-Arrhenius dynamics are correlated with the formation of icosahedral structures. Upon cooling, these icosahedra organize into mesoclusters. We recast this glassformer as an effective system of icosahedra which we describe with a population dynamics model. This model we parameterize with data from the temperature regime accessible to molecular dynamics simulations. We then use the model to determine the population of icosahedra in mesoclusters at arbitrary temperature. Using simulation data to incorporate dynamics into the model, we predict relaxation behavior at temperatures inaccessible to conventional approaches. Our model predicts super-Arrhenius dynamics whose relaxation time remains finite for non-zero temperature.

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

Evolutionary reinforcement learning of dynamical large deviations

In this work, we show how to bound and calculate the likelihood of dynamical large deviations using evolutionary reinforcement learning. An agent, a stochastic model, propagates a continuous-time Monte Carlo trajectory and receives a reward conditioned upon the values of certain path-extensive quantities. Evolution produces progressively fitter agents, potentially allowing the calculation of a piece of a large-deviation rate function for a particular model and path-extensive quantity. For models with small state spaces, the evolutionary process acts directly on rates, and for models with large state spaces, the process acts on the weights of a neural network that parameterizes the model's rates. This approach shows how path-extensive physics problems can be considered within a framework widely used in machine learning.

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

Water and vapor transport in algal‐fungal lichen: Modeling constrained by laboratory experiments, an application for Flavoparmelia caperata

ABSTRACT Algal‐fungal symbionts share water, nutrients, and gases via an architecture unique to lichens. Because lichen activity is controlled by moisture dynamics, understanding water transport is prerequisite to understand their fundamental biology. We propose a model of water distributions within foliose lichens governed by laws of fluid motion. Our model differentiates between water stored in symbionts, on extracellular surfaces, and in distinct morphological layers. We parameterize our model with hydraulic properties inverted from laboratory measurements of Flavoparmelia caperata and validate for wetting and drying. We ask: (1) Where is the bottleneck to water transport? (2) How do hydration and dehydration dynamics differ? and (3) What causes these differences? Resistance to vapor flow is concentrated at thallus surfaces and acts as the bottleneck for equilibrium, while internal resistances are small. The model captures hysteresis in hydration and desiccation, which are shown to be controlled by nonlinearities in hydraulic capacitance. Muting existing nonlinearities slowed drying and accelerated wetting, while exaggerating nonlinearities accelerated drying and slowed wetting. The hydraulic nonlinearity of F. caperata is considerable, which may reflect its preference for humid and stable environments. The model establishes the physical foundation for future investigations of transport of water, gas, and sugar between symbionts.

Potkay, Aaron↗
DOE OSTI2013

Model-based Analysis of HER Activation in Cells Co-Expressing EGFR, HER2 and HER3.

The HER/ErbB family of receptor tyrosine kinases drive critical responses in normal physiology and cancer, and the expression levels of the various HER receptors are critical determinants of clinical outcomes. HER activation is driven by the formation of various dimer complexes between members of this receptor family. The HER dimer types can have differential effects on downstream signaling and phenotypic outcomes. We constructed an integrated mathematical model of HER activation and trafficking to quantitatively link receptor expression levels to dimerization and activation. We parameterized the model with a comprehensive set of HER phosphorylation and abundance data collected in a panel of human mammary epithelial cells expressing varying levels of EGFR, HER2 and HER3. Although parameter estimation yielded multiple solutions, predictions for dimer phosphorylation were in agreement with each other. We validated the model using experiments where pertuzumab was used to block HER2 dimerization. We used the model to predict HER dimerization and activation patterns in a panel of epithelial cells lines with known HER expression levels. Simulations over the range of expression levels seen in various cell lines indicate that: i) EGFR phosphorylation is driven by HER1/1 and HER1/2 dimers, and not HER1/3 dimers, ii) HER1/2 and HER2/3 dimers both contribute significantly to HER2 activation with the EGFR expression level determining the relative importance of these species, and iii) the HER2/3 dimer is largely responsible for HER3 activation. The model can be used to predict phosphorylated dimer levels for any given HER expression profile. This information in turn can be used to quantify the potencies of the various HER dimers, and can potentially inform personalized therapeutic approaches.

ErbB↗
DOE OSTI2015

Moderate forest disturbance as a stringent test for gap and big-leaf models

Disturbance-induced tree mortality is a key factor regulating the carbon balance of a forest, but tree mortality and its subsequent effects are poorly represented processes in terrestrial ecosystem models. It is thus unclear whether models can robustly simulate moderate (non-catastrophic) disturbances, which tend to increase biological and structural complexity and are increasingly common in aging US forests. We tested whether three forest ecosystem models – Biome-BGC (BioGeochemical Cycles), a classic big-leaf model, and the ZELIG and ED (Ecosystem Demography) gap-oriented models – could reproduce the resilience to moderate disturbance observed in an experimentally manipulated forest (the Forest Accelerated Succession Experiment in northern Michigan, USA, in which 38% of canopy dominants were stem girdled and compared to control plots). Each model was parameterized, spun up, and disturbed following similar protocols and run for 5 years post-disturbance. The models replicated observed declines in aboveground biomass well. Biome-BGC captured the timing and rebound of observed leaf area index (LAI), while ZELIG and ED correctly estimated the magnitude of LAI decline. None of the models fully captured the observed post-disturbance C fluxes, in particular gross primary production or net primary production (NPP). Biome-BGC NPP was correctly resilient but for the wrong reasons, and could not match the absolute observational values. ZELIG and ED, in contrast, exhibited large, unobserved drops in NPP and net ecosystem production. The biological mechanisms proposed to explain the observed rapid resilience of the C cycle are typically not incorporated by these or other models. It is thus an open question whether most ecosystem models will simulate correctly the gradual and less extensive tree mortality characteristic of moderate disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Large Divergence of Projected High Latitude Vegetation Composition and Productivity Due To Functional Trait Uncertainty

Abstract Vegetation distribution and composition are expected to change in northern high latitudes under rapid warming, which regulates ecosystem functions but remains challenging to predict. Vegetation change arises from the interplay of chronic climate trends such as warming and transient demographic processes of recruitment, growth, competition, and mortality. Most predictive models overlooked the role of demographic dynamics controlled by plant traits. Here, we simulate vegetation dynamics at the Kougarok Hillslope site in Alaska under historical and future climates using the E3SM Land Model coupled to the Functionally Assembled Terrestrial Simulator (ELM‐FATES). To evaluate the roles of plant traits, we parameterize the model with 5,265 trait configurations representing diverse physiological and demographic strategies. Results show current modeled biomass, composition, and productivity are most sensitive to traits controlling photosynthetic capacity, carbon allocation, allometry, and phenology. Among all trait configurations, ∼5% reproduce in situ biomass and plant functional type (PFT) composition measured in 2016, that are indistinguishable from these two observed ecosystem states. Notably, these same trait configurations produce diverging biomass, composition, and productivity under future climate, where the uncertainty attributable to traits is twice the change attributable to climate change. The variation of projected productivity arises from emerging PFT composition under novel climate regimes, primarily explained by traits controlling cold‐induced mortality, recruitment, and allometry. Our findings highlight the importance and uncertainty of demographic dynamics and its interaction with climate change in shaping Arctic vegetation change. Improved model predictions will likely benefit from explicit consideration of vegetation demography and better constraints of critical traits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Physically-based strength model of tantalum incorporating effects of temperature, strain rate and pressure

In this work, we develop a tantalum strength model that incorporates e ects of temperature, strain rate and pressure. Dislocation kink-pair theory is used to incorporate temperature and strain rate e ects while the pressure dependent yield is obtained through the pressure dependent shear modulus. Material constants used in the model are parameterized from tantalum single crystal tests and polycrystalline ramp compression experiments. It is shown that the proposed strength model agrees well with the temperature and strain rate dependent yield obtained from polycrystalline tantalum experiments. Furthermore, the model accurately reproduces the pressure dependent yield stresses up to 250 GPa. The proposed strength model is then used to conduct simulations of a Taylor cylinder impact test and validated with experiments. This approach provides a physically-based multi-scale strength model that is able to predict the plastic deformation of polycrystalline tantalum through a wide range of temperature, strain and pressure regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Global soil consumption of atmospheric carbon monoxide: an analysis using a process-based biogeochemistry model

Carbon monoxide (CO) plays an important role in controlling the oxidizing capacity of the atmosphere by reacting with OH radicals that affect atmospheric methane (CH 4 ) dynamics. We develop a process-based biogeochemistry model to quantify the CO exchange between soils and the atmosphere with a 5 min internal time step at the global scale. The model is parameterized using the CO flux data from the field and laboratory experiments for 11 representative ecosystem types. The model is then extrapolated to global terrestrial ecosystems using monthly climate forcing data. Global soil gross consumption, gross production, and net flux of the atmospheric CO are estimated to be from -197 to -180, 34 to 36, and -163 to -145 Tg CO yr -1 (1 Tg = 10 12 g), respectively, when the model is driven with satellite-based atmospheric CO concentration data during 2000–2013. Tropical evergreen forest, savanna and deciduous forest areas are the largest sinks at 123 Tg CO yr -1 . The soil CO gross consumption is sensitive to air temperature and atmospheric CO concentration, while the gross production is sensitive to soil organic carbon (SOC) stock and air temperature. By assuming that the spatially distributed atmospheric CO concentrations (~ 128 ppbv) are not changing over time, the global mean CO net deposition velocity is estimated to be 0.16–0.19 mm s -1 during the 20th century. Under the future climate scenarios, the CO deposition velocity will increase at a rate of 0.0002–0.0013 mm s -1 yr -1 during 2014–2100, reaching 0.20–0.30 mm s -1 by the end of the 21st century, primarily due to the increasing temperature. Areas near the Equator, the eastern US, Europe and eastern Asia will be the largest sinks due to optimum soil moisture and high temperature. The annual global soil net flux of atmospheric CO is primarily controlled by air temperature, soil temperature, SOC and atmospheric CO concentrations, while its monthly variation is mainly determined by air temperature, precipitation, soil temperature and soil moisture.

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

Modeling the spatiotemporal variability in subsurface thermal regimes across a low-relief polygonal tundra landscape

Vast carbon stocks stored in permafrost soils of Arctic tundra are under risk of release to the atmosphere under warming climate scenarios. Ice-wedge polygons in the low-gradient polygonal tundra create a complex mosaic of microtopographic features. This microtopography plays a critical role in regulating the fine-scale variability in thermal and hydrological regimes in the polygonal tundra landscape underlain by continuous permafrost. Modeling of thermal regimes of this sensitive ecosystem is essential for understanding the landscape behavior under the current as well as changing climate. Here, we present an end-to-end effort for high-resolution numerical modeling of thermal hydrology at real-world field sites, utilizing the best available data to characterize and parameterize the models. We also develop approaches to model the thermal hydrology of polygonal tundra and apply them at four study sites near Barrow, Alaska, spanning across low to transitional to high-centered polygons, representing a broad polygonal tundra landscape. A multiphase subsurface thermal hydrology model (PFLOTRAN) was developed and applied to study the thermal regimes at four sites. Using a high-resolution lidar digital elevation model (DEM), microtopographic features of the landscape were characterized and represented in the high-resolution model mesh. The best available soil data from field observations and literature were utilized to represent the complex heterogeneous subsurface in the numerical model. Simulation results demonstrate the ability of the developed modeling approach to capture – without recourse to model calibration – several aspects of the complex thermal regimes across the sites, and provide insights into the critical role of polygonal tundra microtopography in regulating the thermal dynamics of the carbon-rich permafrost soils. Moreover, areas of significant disagreement between model results and observations highlight the importance of field-based observations of soil thermal and hydraulic properties for modeling-based studies of permafrost thermal dynamics, and provide motivation and guidance for future observations that will help address model and data gaps affecting our current understanding of the system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Modeling the Spatio-Temporal Variability in Subsurface Thermal Regimes Across a Low-Relief Polygonal Tundra Landscape: Modeling Archive

Vast carbon stocks stored in permafrost soils of Arctic tundra are under risk of release to atmosphere under warming climate. Ice--wedge polygons in the low-gradient polygonal tundra create a complex mosaic of microtopographic features. The microtopography plays a critical role in regulating the fine scale variability in thermal and hydrological regimes in the polygonal tundra landscape underlain by continuous permafrost. Modeling of thermal regimes of this sensitive ecosystem is essential for understanding the landscape behaviour under current as well as changing climate. We present here an end-to-end effort for high resolution numerical modeling of thermal hydrology at real-world field sites, utilizing the best available data to characterize and parameterize the models. We develop approaches to model the thermal hydrology of polygonal tundra and apply them at four study sites at Barrow, Alaska spanning across low to transitional to high-centered polygon and representative of broad polygonal tundra landscape. A multi--phase subsurface thermal hydrology model (PFLOTRAN) was developed and applied to study the thermal regimes at four sites. Using high resolution LiDAR DEM, microtopographic features of the landscape were characterized and represented in the high resolution model mesh. Best available soil data from field observations and literature was utilized to represent the complex heterogeneous subsurface in the numerical model. This data collection provides the complete set of input files, forcing data sets and computational meshes for simulations using PFLOTRAN for four sites at Barrow Environmental Observatory. It also documents the complete computational workflow for this modeling study to allow verification, reproducibility and follow up studies. This dataset includes two zipped files and one .pdf file.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Water vapor distribution in protoplanetary disks

Water vapor has been detected in protoplanetary disks. In this work, we model the distribution of water vapor in protoplanetary disks with a thermo-chemical code. For a set of parameterized disk models, we calculate the distribution of dust temperature and radiation field of the disk with a Monte Carlo method, and then solve the gas temperature distribution and chemical composition. The radiative transfer includes detailed treatment of scattering by atomic hydrogen and absorption by water of Lyα photons, since the Lyα line dominates the UV spectrum of accreting young stars. In a fiducial model, we find that warm water vapor with temperature around 300 K is mainly distributed in a small and well-confined region in the inner disk. The inner boundary of the warm water region is where the shielding of UV field due to dust and water itself become significant. The outer boundary is where the dust temperature drops below the water condensation temperature. A more luminous central star leads to a more extended distribution of warm water vapor, while dust growth and settling tends to reduce the amount of warm water vapor. Based on typical assumptions regarding the elemental oxygen abundance and the water chemistry, the column density of warm water vapor can be as high as 10{sup 22} cm{sup –2}. A small amount of hot water vapor with temperature higher than ∼300 K exists in a more extended region in the upper atmosphere of the disk. Cold water vapor with temperature lower than 100 K is distributed over the entire disk, produced by photodesorption of the water ice.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Final Report – Study of Shortwave Spectra in Fully 3D Environment. Synergy Between Scanning Radars and Spectral Radiation Measurements

ARM set out 20 years ago to “close” the radiation problem, that is, to improve radiation models to the point where they could routinely predict the observed spectral radiation fluxes knowing the optical properties of the surface and of gases, clouds and aerosols in the atmosphere. Only then could such radiation models form a proper springboard for global climate model (GCM) parameterizations of spectral radiation. Sustained efforts have more or less achieved that goal with regard to longwave radiation; ASR models now routinely predict ARM spectral longwave radiances to 1–2%. Similar efforts in the shortwave have achieved far less; the successes are mainly for carefully selected 1D stratiform cloud cases. Such cases amount, even with the most optimistic interpretation, to no more than 30% of all cases at SGP. The problem has not been lack of effort but lack of appropriate instruments.The new ARM stimulus-funded instruments, with their new capabilities, will dramatically improve this situation and once again make progress possible on the shortwave problem. The new shortwave spectrometers will provide a reliable, calibrated record including the near infrared – and for other climatic regimes than SGP. The new scanning radars will provide the 3D cloud view, making it possible to tackle fully 3D situations. Thus, our main theme for the project is the understanding and closure of the surface spectral shortwave radiation problem in fully 3D cloud situations by combining the new ARM scanning radars and shortwave spectrometers with the arsenal of radiative transfer tools.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessing climate change impacts on live fuel moisture and wildfire risk using a hydrodynamic vegetation model

Abstract. Live fuel moisture content (LFMC) plays a critical role in wildfire dynamics, but little is known about responses of LFMC to multivariate climate change, e.g., warming temperature, CO2 fertilization, and altered precipitation patterns, leading to a limited prediction ability of future wildfire risks. Here, we use a hydrodynamic demographic vegetation model to estimate LFMC dynamics of chaparral shrubs, a dominant vegetation type in fire-prone southern California. We parameterize the model based on observed shrub allometry and hydraulic traits and evaluate the model's accuracy through comparisons between observed and simulated LFMC of three plant functional types (PFTs) under current climate conditions. Moreover, we estimate the number of days per year of LFMC below 79 % (which is a critical threshold for wildfire danger rating of southern California chaparral shrubs) from 1960 to 2099 for each PFT and compare the number of days below the threshold for medium and high greenhouse gas emission scenarios (RCP4.5 and 8.5). We find that climate change could lead to more days per year (5.2 %–14.8 % increase) with LFMC below 79 % between the historical (1960–1999) and future (2080–2099) periods, implying an increase in wildfire danger for chaparral shrubs in southern California. Under the high greenhouse gas emission scenario during the dry season, we find that the future LFMC reductions mainly result from a warming temperature, which leads to 9.1 %–18.6 % reduction in LFMC. Lower precipitation in the spring leads to a 6.3 %–8.1 % reduction in LFMC. The combined impacts of warming and precipitation change on fire season length are equal to the additive impacts of warming and precipitation change individually. Our results show that the CO2 fertilization will mitigate fire risk by causing a 3.5 %–4.8 % increase in LFMC. Our results suggest that multivariate climate change could cause a significant net reduction in LFMC and thus exacerbate future wildfire danger in chaparral shrub systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗
DOE OSTI2021

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗
DOE OSTI2021

Improved Convective Ice Microphysics Parameterization in the NCAR CAM Model

Partitioning deep convective cloud condensates into components that sediment and detrain, known to be a challenge for global climate models, is important for cloud vertical distribution and anvil cloud formation. In this study, we address this issue by improving the convective microphysics scheme in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 (CAM5.3). The improvements include: (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. Results from the default and modified CAM5.3 models were evaluated against observations from the U.S. Department of Energy Tropical Warm Pool-International Cloud Experiment (TWP-ICE) field campaign. The default model overestimates ice amount, which is largely attributed to the underestimation of convective ice particle sedimentation. By considering cloud ice sedimentation and rimed ice particles and applying a new convective snow terminal velocity parameterization, the vertical distribution of ice amount is much improved in the midtroposphere and upper troposphere when compared to observations. The vertical distribution of ice condensate also agrees well with observational best estimates upon considering snow detrainment. Further, comparison with observed convective updrafts reveals that current bulk model fails to reproduce the observed updraft magnitude and occurrence frequency, suggesting spectral distributions be required to simulate the subgrid updraft heterogeneity.

54 ENVIRONMENTAL SCIENCES↗