Search NASA⌕ Search

SEARCH · Search NASA

Results for “subgrid 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Biophysical Impact of Land-Use and Land-Cover Change on Subgrid Temperature in CMIP6 Models

In this study, we investigate the air temperature response to land-use and land-cover change (LULCC; cropland expansion and deforestation) using subgrid land model output generated by a set of CMIP6 model simulations. Our study is motivated by the fact that ongoing land-use activities are occurring at local scales, typically significantly smaller than the resolvable scale of a grid cell in Earth system models. It aims to explore the potential for a multimodel approach to better characterize LULCC local climatic effects. On an annual scale, the CMIP6 models are in general agreement that croplands are warmer than primary and secondary land (psl; mainly forests, grasslands, and bare ground) in the tropics and cooler in the mid–high latitudes, except for one model. The transition from warming to cooling occurs at approximately 40°N. Although the surface heating potential, which combines albedo and latent heat flux effects, can explain reasonably well the zonal mean latitudinal subgrid temperature variations between crop and psl tiles in the historical simulations, it does not provide a good prediction on subgrid temperature for other land tile configurations (crop vs forest; grass vs forest) under Shared Socioeconomic Pathway 5–8.5 (SSP5–8.5) forcing scenarios. A subset of simulations with the CESM2 model reveals that latitudinal subgrid temperature variation is positively related to variation in net surface shortwave radiation and negatively related to variation in the surface energy redistribution factor, with a dominant role from the latter south of 30°N. We suggest that this emergent relationship can be used to benchmark the performance of land surface parameterizations and for prediction of local temperature response to LULCC.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

Sensitivity of Middle Atmospheric Temperature and Circulation in the UIUC Mesosphere-Stratosphere-Troposphere GCM to the Treatment of Subgrid-Scale Gravity-Wave Breaking

The sensitivity of the middle atmospheric temperature and circulation to the treatment of mean- flow forcing due to breaking gravity waves was investigated using the University of Illinois at Urbana-Champaign 40-layer Mesosphere-Stratosphere-Troposphere General Circulation Model (MST-GCM). Three GCM experiments were performed. The gravity-wave forcing was represented first by Rayleigh friction, and then by the Alexander and Dunkerton (AD) parameterization with weak and strong breaking effects of gravity waves. In all experiments, the Palmer et al. parameterization was included to treat the breaking of topographic gravity waves in the troposphere and lower stratosphere. Overall, the experiment with the strong breaking effect simulates best the middle atmospheric temperature and circulation. With Rayleigh friction and the weak breaking effect, a large warm bias of up to 60 C was found in the summer upper mesosphere and lower thermosphere. This warm bias was linked to the inability of the GCM to simulate the reversal of the zonal winds from easterly to westerly crossing the mesopause in the summer hemisphere. With the strong breaking effect, the GCM was able to simulate this reversal, and essentially eliminated the warm bias. This improvement was the result of a much stronger meridional transport circulation that possesses a strong vertical ascending branch in the summer upper mesosphere, and hence large adiabatic cooling. Budget analysis indicates that 'in the middle atmosphere the forces that act to maintain a steady zonal-mean zonal wind are primarily those associated with the meridional transport circulation and breaking gravity waves. Contributions from the interaction of the model-resolved eddies with the mean flow are small. To obtain a transport circulation in the mesosphere of the UIUC MST-GCM that is strong enough to produce the observed cold summer mesopause, gravity-wave forcing larger than 100 m/s/day in magnitude is required near the summer mesopause. In the tropics, only with the AD parameterization can the model produce realistic semiannual oscillations.

Yang, Fanglin↗
NASA NTRS1994

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗
NASA NTRS2004

Temporal Large-Eddy Simulation

In 1999, Stolz and Adams unveiled a subgrid-scale model for LES based upon approximately inverting (defiltering) the spatial grid-filter operator and termed .the approximate deconvolution model (ADM). Subsequently, the utility and accuracy of the ADM were demonstrated in a posteriori analyses of flows as diverse as incompressible plane-channel flow and supersonic compression-ramp flow. In a prelude to the current paper, a parameterized temporal ADM (TADM) was developed and demonstrated in both a priori and a posteriori analyses for forced, viscous Burger's flow. The development of a time-filtered variant of the ADM was motivated-primarily by the desire for a unifying theoretical and computational context to encompass direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaged Navier-Stokes simulation (RANS). The resultant methodology was termed temporal LES (TLES). To permit exploration of the parameter space, however, previous analyses of the TADM were restricted to Burger's flow, and it has remained to demonstrate the TADM and TLES methodology for three-dimensional flow. For several reasons, plane-channel flow presents an ideal test case for the TADM. Among these reasons, channel flow is anisotropic, yet it lends itself to highly efficient and accurate spectral numerical methods. Moreover, channel-flow has been investigated extensively by DNS, and a highly accurate data base of Moser et.al. exists. In the present paper, we develop a fully anisotropic TADM model and demonstrate its utility in simulating incompressible plane-channel flow at nominal values of Re(sub tau) = 180 and Re(sub tau) = 590 by the TLES method. The TADM model is shown to perform nearly as well as the ADM at equivalent resolution, thereby establishing TLES as a viable alternative to LES. Moreover, as the current model is suboptimal is some respects, there is considerable room to improve TLES.

Pruett, C. D.↗
DOE OSTI2019

The Second Wind Forecast Improvement Project (WFIP 2): General Overview

In 2015 the U.S. Department of Energy initiated a four-year study, the second Wind Forecast Improvement Project (WFIP 2), to improve the representation of boundary-layer physics and related processes in mesoscale models for better wind and wind power forecasts. This goal challenges numerical weather prediction models in complex terrain due to inherent assumptions underlying their boundary-layer parameterizations. The WFIP 2 effort involved the wind industry, universities, NOAA, and DOE's national laboratories in an integrated observational and modeling study. Observations spanned 18-months to assure a full annual cycle of continuously recorded observations from remote-sensing and in situ measurement systems. The study area comprised the Columbia Basin of eastern Washington and Oregon, containing more than 4 GW of installed wind capacity; nests of observational systems captured important atmospheric scales from mesoscale to NWP sub-grid scale. Model improvements targeted NOAA's High-Resolution Rapid Refresh (HRRR) model to facilitate transfer of improvements to NWS operational forecast models. This paper describes the general WFIP 2 scope and objectives, the particular scientific challenges of improving wind forecasts in complex terrain, early successes of the project, and an integrated approach to archiving observations and model output. It provides an introduction for a set of more detailed BAMS papers addressing WFIP 2 observational science, modeling challenges and solutions, practical transfer of boundary layer advances to industry use, and advances in coupling improved mesoscale models to microscale models that can represent interactions between wind plants and the atmosphere.

WFIP2, boundary layer, complex terrain, forecastin↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Exposing and Reducing Biases of Simulating Mixed-Phase Clouds in the Convection-Permitting E3SM Atmosphere Model: Lessons From an Arctic Cold-Air Outbreak

Mixed-phase clouds modulate the water and energy cycles of high-latitude regions, yet their liquid-ice phase partitioning has long been poorly simulated in climate models. Here, simulations of Arctic mixed-phase clouds by the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) are assessed against large-eddy simulations, satellite data, and ground-based observations during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment field campaign. SCREAM simulates nearly completely frozen clouds, which is attributed largely to the unreasonably strong Wegener–Bergeron–Findeisen (WBF) process that converts liquid to ice excessively and partly to the early over-abundant ice production at cold temperatures from a temperature-deterministic deposition ice nucleation scheme. Assuming no subgrid variation for the WBF process in the original formulation particularly conflicts with the instantaneous saturation adjustment assumption in the condensation scheme that assumes subgrid variability, leading to exaggerated WBF process rates. A proposed simple physically-based improvement on the treatment of subgrid cloud overlap substantially increases supercooled liquid water content and notably improves cloud-top phase partitioning, aligning better with observations. Improvement of supercooled liquid water content also converges with increasing horizontal resolution. The deposition ice nucleation scheme is found responsible for a falsely-produced ice cloud aloft that is not observed, biasing the simulated cloud radiative effects and top-of-atmosphere radiative fluxes. This study identifies key deficiencies in cloud parameterizations that continue to challenge convection-permitting models.

Geosciences↗
DOE OSTI2026

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗
NASA NTRS1995

The importance of mesoscale circulations generated by subgrid-scale landscape heterogeneities in general circulation models

A mesoscale atmospheric model was used to evaluate the impact of subgrid-scale landscape discontinuities on the vertical profiles of resolved temperature, moisture, and moist static energy in the planetary boundary layer (PBL) of general circulation models (GCMs). These profiles were produced with a 3D version of the model (using a horizontal grid resolution of 7.5 km and 13 vertical layers in the PBL) by averaging horizontally the various atmospheric variables over a 180 x 180 sq km domain-about the size of the horizontal domain represented by a single grid element in a GCM. They were compared to corresponding vertical profiles produced with a 1D version of the model, which simulates the PBL, as in a GCM, over a single horizontal grid element. Differences obtained between the horizontally averaged atmospheric variables produced with the 3D simulations and the 1D simulations emphasize the impact of subgrid-scale landscape discontinuities on GCM-resolved variables. Various types of landscape discontinuities, characterized by horizontal contrasts of surface wetness and size of land patches, were simulated under various background-wind conditions. Differences of temperature, specific humidity, and moist static energy as large as 4 K, 6 g/kg, and 10 kJ/kg were obtained in some cases. These differences were not affected significantly by moderate winds but were sensitive to the spatial distribution of surface wetness. These results emphasize the need to parameterize mesoscale processes induced by landscape discontinuities in GCMss.

Lynn, Barry H.↗
DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗
DOE OSTI2022

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Crossing Multiple Gray Zones in the Transition from Mesoscale to Microscale Simulation over Complex Terrain

This review paper explores the field of mesoscale to microscale modeling over complex terrain as it traverses multiple so-called gray zones. In an attempt to bridge the gap between previous large-scale and small-scale modeling efforts, atmospheric simulations are being run at an unprecedented range of resolutions. The gray zone is the range of grid resolutions where particular features are neither subgrid nor fully resolved, but rather are partially resolved. The definition of a gray zone depends strongly on the feature being represented and its relationship to the model resolution. This paper explores three gray zones relevant to simulations over complex terrain: turbulence, convection, and topography. Taken together, these may be referred to as the gray continuum. The focus is on horizontal grid resolutions from ∼10 km to ∼10 m. In each case, the challenges are presented together with recent progress in the literature. A common theme is to address cross-scale interaction and scale-awareness in parameterization schemes. How numerical models are designed to cross these gray zones is critical to complex terrain applications in numerical weather prediction, wind resource forecasting, and regional climate modeling, among others.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Parameterization of Mixed Layer and Deep-Ocean Mesoscales Including Nonlinearity

In 2011, Chelton et al. carried out a comprehensive census of mesoscales using altimetry data and reached the following conclusions: "essentially all of the observed mesoscale features are nonlinear" and "mesoscales do not move with the mean velocity but with their own drift velocity," which is "the most germane of all the nonlinear metrics." Accounting for these results in a mesoscale parameterization presents conceptual and practical challenges since linear analysis is no longer usable and one needs a model of nonlinearity. A mesoscale parameterization is presented that has the following features: 1) it is based on the solutions of the nonlinear mesoscale dynamical equations, 2) it describes arbitrary tracers, 3) it includes adiabatic (A) and diabatic (D) regimes, 4) the eddy-induced velocity is the sum of a Gent and McWilliams (GM) term plus a new term representing the difference between drift and mean velocities, 5) the new term lowers the transfer of mean potential energy to mesoscales, 6) the isopycnal slopes are not as flat as in the GM case, 7) deep-ocean stratification is enhanced compared to previous parameterizations where being more weakly stratified allowed a large heat uptake that is not observed, 8) the strength of the Deacon cell is reduced. The numerical results are from a stand-alone ocean code with Coordinated Ocean-Ice Reference Experiment I (CORE-I) normal-year forcing.

Baroclinic flows; Mesoscale processes; Mixing; Sub↗
NASA NTRS2000

A Catchment-Based Approach to Modeling Land Surface Processes in a GCM: Model Structure - Part 1

A new strategy for modeling the land surface component of the climate system is described. The strategy is motivated by an arguable deficiency in most state-of-the-art land surface models (LSMs), namely the disproportionately higher emphasis given to the formulation of one-dimensional, vertical physics relative to the treatment of horizontal heterogeneity in surface properties -- particularly subgrid soil moisture variability and its effects on runoff generation. The new strategy calls for the partitioning of the continental surface into a mosaic of hydrologic catchments, delineated through analysis of high-resolution surface elevation data. The effective "grid" used for the land surface is therefore not specified by the overlying atmospheric grid. Within each catchment, the variability of soil moisture is related to characteristics of the topography and to three bulk soil moisture variables through a well-established model of catchment processes. This modeled variability allows the partitioning of the catchment into several areas representing distinct hydrological regimes, wherein distinct (regime-specific) evaporation and runoff parameterizations are applied. Care is taken to ensure that the deficiencies of the catchment model in regions of little to moderate topography are minimized.

Koster, Randal D.↗
DOE OSTI2020

Impacts of Stochastic Mixing in Idealized Convection-Permitting Simulations of Squall Lines

This study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δ h ) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δ h simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δ h . It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δ h compared to 0.5-km Δ h . Further, the stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

On the Subgrid-Scale Modeling of Compressible Turbulence

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Squires, Kyle↗
DOE OSTI2021

Tropical Cirrus in Global Storm-Resolving Models. Part II: Cirrus Life Cycle and Top-of-Atmosphere Radiative Fluxes

Cirrus clouds of various thicknesses and radiative characteristics extend over much of the tropics, especially around deep convection. They are difficult to observe due to their high altitude and sometimes small optical depths. They are also difficult to simulate in conventional global climate models, which have coarse grid spacings and simplified parameterizations of deep convection and cirrus formation. We investigate the representation of tropical cirrus in global storm-resolving models (GSRMs), which have higher spatial resolution and explicit convection and could more accurately represent cirrus cloud processes. This study uses GSRMs from the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project. The aggregate life cycle of tropical cirrus is analyzed using joint albedo and outgoing longwave radiation (OLR) histograms to assess the fidelity of models in capturing the observed cirrus cloud populations over representative tropical ocean and land regions. The proportions of optically-thick deep convection, anvils, and cirrus vary across models and are portrayed in the vertical distribution of cloud cover and top-of-atmosphere radiative fluxes. Model differences in cirrus populations, likely driven by subgrid processes such as ice microphysics, dominate over regional differences between convectively-active tropical land and ocean locations.

54 ENVIRONMENTAL SCIENCES↗