Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temperature Gradient”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

A Qualitative Investigation of Selected Infrared Flow Visualization Image Processing Techniques

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel testenvironment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is discussed as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity.This method makes use of dynamic mode decomposition and discrete-time filtering to separate the background reflections that dominate low emissivity surfaces from the aerodynamic driven surface temperature gradients.This process will be shown to resolve the onset of boundary layer transition on a bare metal wing as well as identify and resolve hidden features in the image data. While the implementation of this technique is very preliminary it demonstrates the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Image Processeing↗
DOE OSTI2022

First-Principles Calculations of Thermoelectric Transport Properties of Quaternary and Ternary Bulk Chalcogenide Crystals

Chalcogenide crystals have a wide range of applications, especially as thermoelectric materials for energy conversion. Thermoelectric materials can be used to generate an electric current from a temperature gradient based on the Seebeck effect and based on the Peltier effect, and they can be used in cooling applications. Using first-principles calculations and semiclassical Boltzmann theory, we have computed the Seebeck coefficient, electrical conductivity, electronic thermal conductivity, power factor, and figure of merit of 30 chalcogenide crystals. A Quantum Espresso package is used to calculate the electronic properties and locate the Fermi level. The transport properties are then calculated using the BoltzTraP code. The 30 crystals are divided into two groups. The first group has four crystals with quaternary composition (A 2 BCQ 4 ) (A = Tl; B = Cd, Hg; C = Si, Ge, Sn; Q = S, Se, Te). The second group contains 26 crystals with the ternary composition (A’B’Q 2 ) (A’ = Ag, Cu, Au, Na; B’ = B, Al, Ga, In; Q = S, Se, Te). Among these 30 chalcogenide crystals, the results for 11 crystals: Tl 2 CdGeSe 4 , Tl 2 CdSnSe 4 , Tl 2 HgSiSe 4 , Tl 2 HgSnS 4 , AuBSe 2 , AuBTe 2 , AuAlTe 2 , AuGaTe 2 , AuInTe 2 , AgAlSe 2 , and AgAlTe 2 are revealed for the first time. In addition, temperature-dependent transport properties of pure and doped AgSbSe 2 and AgSbTe 2 crystals with dopant compositions of AgSb 0.94 Cd 0.06 Te 2 and AgSbTe 1.85 Se 0.15 were explored. These results provide an excellent database for bulk chalcogenides crucial for a wide range of potential applications in renewable energy fields.

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

Probing the Response of Tropical Deep Convection to Aerosol Perturbations Using Idealized Cloud-Resolving Simulations with Parameterized Large-Scale Dynamics

Abstract A framework is introduced to investigate the indirect effect of aerosol loading on tropical deep convection using three-dimensional limited-domain idealized cloud-system-resolving model simulations coupled with large-scale dynamics over fixed sea surface temperature. The large-scale circulation is parameterized using the spectral weak temperature gradient (WTG) approximation that utilizes the dominant balance between adiabatic cooling and diabatic heating in the tropics. The aerosol loading effect is examined by varying the number of cloud condensation nuclei (CCN) available to form cloud droplets in the two-moment bulk microphysics scheme over a wide range of environments from 30 to 5000 cm−3. The radiative heating is held at a constant prescribed rate in order to isolate the microphysical effects. Analyses are performed over the period after equilibrium is achieved between convection and the large-scale environment. Mean precipitation is found to decrease modestly and monotonically when the aerosol number concentration increases as convection gets weaker, despite the increase in cloud liquid water in the warm-rain region and ice crystals aloft. This reduction is traced down to the reduction in surface enthalpy fluxes as an energy source to the atmospheric column induced by the coupling of the large-scale motion, though the gross moist stability remains constant. Increasing CCN concentration leads to 1) a cooler free troposphere because of a reduction in the diabatic heating and 2) a warmer boundary layer because of suppressed evaporative cooling. This dipole temperature structure is associated with anomalously descending large-scale vertical motion above the boundary layer and ascending motion at lower levels. Sensitivity tests suggest that changes in convection and mean precipitation are unlikely to be caused by the impact of aerosols on cloud droplets and microphysical properties but rather by accounting for the feedback from convective adjustment with the large-scale dynamics. Furthermore, a simple scaling argument is derived based on the vertically integrated moist static energy budget, which enables estimation of changes in precipitation given known changes in surfaces enthalpy fluxes and the constant gross moist stability. The impact on cloud hydrometeors and microphysical properties is also examined, and it is consistent with the macrophysical picture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Setting the H-mode pedestal structure: variations of particle source location using gas puff and pellet fueling

Experiments in DIII-D show that a particle source location inside the top of the H-mode pedestal (pellets) maintains a higher pedestal pressure than an edge source (gas fueling) through a widening of the electron temperature pedestal with reduction of the temperature gradient. The effect of these two fueling schemes on the H-mode pedestal structure was examined in DIII-D by comparing controlled pellet-fueled and gas-fueled discharges across a fueling scan up to 40 torr 1 s –1 . High resolution electron profiles reveal that gas fueling lowers the pedestal pressure as the density profile shifts radially outwards and the separatrix density increases, while pellet fueling maintains a constant pedestal pressure. The neutral source locations from pellets and gas are determined with the PELLET and UEDGE codes, respectively, and quantify the particle source localization. Pellets provide significant ionization inside the pedestal top while gas puffing localizes ionization in the scrape-off-layer and pedestal foot, broadly consistent with the density profile structure influenced by the source. ELMs are observed to increase in frequency and reduce impurity content as fueling is increased. Stability analysis with ELITE shows that both conditions are near the type-I ELM corner of the peeling–ballooning stability diagram, which is altered significantly by the introduction of pellets. Since transport mechanisms are not observed to change substantially with particle source location, wider pedestals allow the pellet-fueled discharges to retain higher pedestal temperatures at similar pedestal densities. EPED1 is tested to capture the pedestal pressure, under-predicting the height with pellets and over-predicting the height with gas by ~15%. Overall, these results have important implications for future reactors where pellet fueling will be the primary particle source due to an opaque scrape-off-layer by showing that the neutral source location plays a role in setting the structure of the H-mode pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

The Response of the Quasi-Biennial Oscillation to Increased CO and its Modulation by Composition Feedbacks

The Quasi-Biennial Oscillation (QBO) impacts the large-scale circulation by altering lower stratospheric meridional temperature gradients which alter the propagation of upward propagating planetary waves. Proposed teleconnections of the QBO include impacts on the stratospheric polar vortices, extratropical surface winter climate, and the Madden-Julian Oscillation. However, long-term projections of the QBO remain highly uncertain. While recent multi-model investigations (Richter et al. 2019, Butchart et al. 2020) show that the amplitude of the QBO weakens robustly among models, changes in QBO period and stability remain highly uncertain. Here we examine the QBO response to increased greenhouse gases using the NASA Goddard Institute for Space Studies Middle Atmosphere Model E2.2 (Rind et al., 2020; Orbe et al. 2020). Compared to lower vertical resolution versions of ModelE, E2.2 also has a higher model top (0.002 hPa) and employs additional interactive non-orographic gravity wave drag sources from convection and shear, which produce a sufficiently realistic QBO, thus rendering it suitable for use in climate change studies. Overall, we find that both the QBO period and amplitude decrease in response to increased CO , the former related to increased lower stratospheric momentum fluxes (associated with convection) and the latter associated partly with a stronger residual mean circulation. Experiments integrated using fully interactive chemistry also reveal that ozone feedbacks significantly impact the magnitude of the QBO amplitude response. Finally, integrations using fixed (pre-industrial) sea surface temperatures (SST) show that the QBO amplitude responds differently to rapid adjustments versus SST feedbacks, as compared to the QBO period.

Quasi-Biennial Oscillation↗
DOE OSTI2023

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Fiber Optic Sensing Technologies for Battery Management Systems and Energy Storage Applications

Applications of fiber optic sensors to battery monitoring have been increasing due to the growing need of enhanced battery management systems with accurate state estimations. The goal of this review is to discuss the advancements enabling the practical implementation of battery internal parameter measurements including local temperature, strain, pressure, and refractive index for general operation, as well as the external measurements such as temperature gradients and vent gas sensing for thermal runaway imminent detection. A reasonable matching is discussed between fiber optic sensors of different range capabilities with battery systems of three levels of scales, namely electric vehicle and heavy-duty electric truck battery packs, and grid-scale battery systems. The advantages of fiber optic sensors over electrical sensors are discussed, while electrochemical stability issues of fiber-implanted batteries are critically assessed. This review also includes the estimated sensing system costs for typical fiber optic sensors and identifies the high interrogation cost as one of the limitations in their practical deployment into batteries. Finally, future perspectives are considered in the implementation of fiber optics into high-value battery applications such as grid-scale energy storage fault detection and prediction systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗
DOE OSTI2024

Importance of gas heating in capacitively coupled radiofrequency plasma-assisted synthesis of carbon nanomaterials

In pursuit of diamond nanoparticles, a capacitively-coupled radio frequency flow-through plasma reactor was operated with methane-argon gas mixtures. Signatures of the final product obtained microscopically and spectroscopically indicated that the product was an amorphous form of graphite. This result was consistent irrespective of combinations of the macroscopic reactor settings. To explain the observed synthesis output, measurements of C 2 and gas properties were carried out by laser-induced fluorescence and optical emission spectroscopy. Strikingly, the results indicated a strong gas temperature gradient of 100 K per mm from the center of the reactor to the wall. Based on additional plasma imaging, a model of hot constricted region (filamentation region) was then formulated. It illustrated that, while the hot constricted region was present, the bulk of the gas was not hot enough to facilitate diamond sp 3 formation: characterized by much lower reaction rates, when compared to sp 2 , sp 3 formation kinetics are expected to become exponentially slow. This result was further confirmed by experiments under identical conditions but with a H 2 /CH 4 mixture, where no output material was detected: if graphitic sp 2 formation was expected as the main output material from the methane feedstock, atomic hydrogen would then be expected to etch it away in situ, such that the net production of that sp 2 -hybridized solid material is nearly a zero. Finally, the crucial importance of gas heating was corroborated by replacing RF with microwave source whereby facile sp 3 production was attained with H 2 /CH 4 gas mixture.

36 MATERIALS SCIENCE↗
NASA NTRS2022

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗
DOE OSTI2015

Black carbon aerosol-induced Northern Hemisphere tropical expansion

Global climate models (GCMs) underestimate the observed trend in tropical expansion. Recent studies partly attribute it to black carbon (BC) aerosols, which are poorly represented in GCMs. In this paper, we conduct a suite of idealized experiments with the Community Atmosphere Model version 4 coupled to a slab ocean model forced with increasing BC concentrations covering a large swath of the estimated range of current BC radiative forcing while maintaining their spatial distribution. The Northern Hemisphere (NH) tropics expand poleward nearly linearly as BC radiative forcing increases (0.7° W -1 m 2 ), indicating that a realistic representation of BC could reduce GCM biases. We find support for the mechanism where BC-induced midlatitude tropospheric heating shifts the maximum meridional tropospheric temperature gradient poleward resulting in tropical expansion. Finally, we also find that the NH poleward tropical edge is nearly linearly correlated with the location of the Intertropical Convergence Zone, which shifts northward in response to increasing BC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Numerical modeling of cold magmatic CO 2 flux measurements for the exploration of hidden geothermal systems

The most accepted conceptual model to explain surface degassing of cold magmatic CO 2 in volcanic-geothermal systems involves the presence of a gas reservoir. Here, numerical simulations using the TOUGH2-ECO2N V2.0 package are performed to get quantitative insights into how cold CO 2 soil flux measurements are related to reservoir and fluid properties. While the modeling is based on flux data measured at a specific geothermal site, the Acoculco caldera (Mexico), some general insights have been gained. Both the CO 2 fluxes at the surface and the depth at which CO 2 exsolves are highly sensitive to the dissolved CO 2 content of the deep fluid. If CO 2 mainly exsolves above the reservoir within a fracture zone, the surface CO 2 fluxes are not sensitive to the reservoir size but depend on the CO 2 dissolved content and the rock permeability. For gas exsolution below the top of the reservoir, surface CO 2 fluxes also depend on the gas saturation of the deep fluid as well as the reservoir size. The absence of thermal anomalies at the surface is mainly a consequence of the low enthalpy of CO 2 . The heat carried by CO 2 is efficiently cooled down by heat conduction and to a certain extent by isoenthalpic volume expansion depending on the temperature gradient. Thermal anomalies occur at higher CO 2 fluxes (>37,000 g m -2 d -1 ) when the heat flux of the rising CO 2 is not balanced anymore. Lastly, specific results are obtained for the Acoculco area (reservoir depth, CO 2 dissolved content, and gas saturation state).

58 GEOSCIENCES↗
DOE OSTI2016

Uncertainty in future projections of the North Pacific subtropical high and its implication for California winter precipitation change: FUTURE CALIFORNIA PRECIPITATION

This study examines future projections of sea level pressure change in the North Pacific and its impact on winter precipitation changes in California. The multi-model analysis, based on the Coupled Model Inter-comparison Project phase 5 (CMIP5) models under the Representative Concentration Pathway 8.5 (RCP8.5) scenario, shows a robust sea-level pressure change in the late 21st century over the western North Pacific in which both the Aleutian Low and North Pacific Subtropical High (NPSH) shift poleward in concert with a widening of the Hadley Cell. This change is partly explained by a systematic increase of static stability in the subtropics. However, over the eastern North Pacific, the projected NPSH changes exhibit a substantial inter-model spread, resulting in uncertain projections of precipitation changes in California. This inter-model spread in the eastern North Pacific is associated with a Pacific Decadal Oscillation-like surface temperature change in the western North Pacific and the resulting meridional temperature gradient change. This study points to a major source of uncertainty for the response of winter precipitation to global warming over the West Coast of North America: atmosphere-ocean coupling in the North Pacific.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2015

A microbial functional group-based module for simulating methane production and consumption: Application to an incubated permafrost soil

In this study, accurately estimating methane (CH 4 ) flux is critically important for investigating and predicting the biogeochemistry-climate feedback. Better simulating CH 4 flux requires explicit representations of microbial processes on CH 4 dynamics because all processes for CH 4 production and consumption are actually carried out by microbes. A microbial functional group based module was developed and tested against an incubation experiment. The module considers four key mechanisms for CH 4 production and consumption: methanogenesis from acetate or single-carbon compounds and CH 4 oxidation using molecular oxygen or other inorganic electron acceptors. These four processes were carried out by four microbial functional groups: acetoclastic methanogens, hydrogenotrophic methanogens, aerobic methanotrophs, and anaerobic methanotrophs. This module was then linked with the decomposition subroutine of the Community Land Model , and was further used to simulate dynamics of carbon dioxide (CO 2 ) and CH 4 concentrations from an incubation experiment with permafrost soils. The results show that the model could capture the dynamics of CO 2 and CH 4 concentrations in microcosms with top soils, mineral layer soils and permafrost soils under natural and saturated moisture conditions and a temperature gradient of -2°C, 3°C, and 5°C. Sensitivity analysis confirmed the importance of acetic acid's direct contribution as substrate and indirect effects through pH feedback on CO 2 and CH 4 production and consumption. This study suggests that representing the microbial mechanisms is critical for modeling CH 4 production and consumption; it is urgent to incorporate microbial mechanisms into Earth system models for better predicting the behavior of the climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Intercomparison of methods of coupling between convection and large-scale circulation. 1. Comparison over uniform surface conditions

Here, as part of an international intercomparison project, a set of single-column models (SCMs) and cloud-resolving models (CRMs) are run under the weak-temperature gradient (WTG) method and the damped gravity wave (DGW) method. For each model, the implementation of the WTG or DGW method involves a simulated column which is coupled to a reference state defined with profiles obtained from the same model in radiative-convective equilibrium. The simulated column has the same surface conditions as the reference state and is initialized with profiles from the reference state. We performed systematic comparison of the behavior of different models under a consistent implementation of the WTG method and the DGW method and systematic comparison of the WTG and DGW methods in models with different physics and numerics. CRMs and SCMs produce a variety of behaviors under both WTG and DGW methods. Some of the models reproduce the reference state while others sustain a large-scale circulation which results in either substantially lower or higher precipitation compared to the value of the reference state. CRMs show a fairly linear relationship between precipitation and circulation strength. SCMs display a wider range of behaviors than CRMs. Some SCMs under the WTG method produce zero precipitation. Within an individual SCM, a DGW simulation and a corresponding WTG simulation can produce different signed circulation. When initialized with a dry troposphere, DGW simulations always result in a precipitating equilibrium state. The greatest sensitivities to the initial moisture conditions occur for multiple stable equilibria in some WTG simulations, corresponding to either a dry equilibrium state when initialized as dry or a precipitating equilibrium state when initialized as moist. Multiple equilibria are seen in more WTG simulations for higher SST. In some models, the existence of multiple equilibria is sensitive to some parameters in the WTG calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

The effect of large-scale model time step and multiscale coupling frequency on cloud climatology, vertical structure, and rainfall extremes in a superparameterized GCM

The effect of global climate model (GCM) time step—which also controls how frequently global and embedded cloud resolving scales are coupled—is examined in the Superparameterized Community Atmosphere Model ver 3.0. Systematic bias reductions of time-mean shortwave cloud forcing (~10 W/m 2 ) and longwave cloud forcing (~5 W/m 2 ) occur as scale coupling frequency increases, but with systematically increasing rainfall variance and extremes throughout the tropics. An overarching change in the vertical structure of deep tropical convection, favoring more bottom-heavy deep convection as a global model time step is reduced may help orchestrate these responses. The weak temperature gradient approximation is more faithfully satisfied when a high scale coupling frequency (a short global model time step) is used. These findings are distinct from the global model time step sensitivities of conventionally parameterized GCMs and have implications for understanding emergent behaviors of multiscale deep convective organization in superparameterized GCMs. Lastly, the results may also be useful for helping to tune them.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Intercomparison of methods of coupling between convection and large-scale circulation: 2. Comparison over nonuniform surface conditions

As part of an international intercomparison project, the weak temperature gradient (WTG) and damped gravity wave (DGW) methods are used to parameterize large-scale dynamics in a set of cloud-resolving models (CRMs) and single column models (SCMs). The WTG or DGW method is implemented using a configuration that couples a model to a reference state defined with profiles obtained from the same model in radiative-convective equilibrium. We investigated the sensitivity of each model to changes in SST, given a fixed reference state. We performed a systematic comparison of the WTG and DGW methods in different models, and a systematic comparison of the behavior of those models using the WTG method and the DGW method. The sensitivity to the SST depends on both the large-scale parameterization method and the choice of the cloud model. In general, SCMs display a wider range of behaviors than CRMs. All CRMs using either the WTG or DGW method show an increase of precipitation with SST, while SCMs show sensitivities which are not always monotonic. CRMs using either the WTG or DGW method show a similar relationship between mean precipitation rate and column-relative humidity, while SCMs exhibit a much wider range of behaviors. DGW simulations produce large-scale velocity profiles which are smoother and less top-heavy compared to those produced by the WTG simulations. Lastly, these large-scale parameterization methods provide a useful tool to identify the impact of parameterization differences on model behavior in the presence of two-way feedback between convection and the large-scale circulation.

54 ENVIRONMENTAL SCIENCES↗