SEARCH · Search NASA
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.
Explore the source record for details and available documents.
A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.
Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.
The middle atmosphere is an important component of the climate system, primarily because of the radiative forcing of ozone. Middle atmospheric ozone can change, over long times, because of changes in the abundance of anthropogenic pollutants which catalytically destroy it, and because of the temperature sensitivity of kinetic reaction rates. There is thus a complex interaction between ozone, involving chemical and climatic mechanisms. One question of interest is how ozone will change over the next decades , as the "greenhouse-gas cooling" of the middle atmosphere increases but the concentrations of chlorine species decreases (because of policy changes). concerns the climate biases in current middle atmosphere-climate models, especially their ability to simulate the correct seasonal cycle at high latitudes, and the existence of temperature biases in the global mean. A major obstacle when addressing this question This paper will present a summary of recent results from the "GCM-Reality Intercomparison Project for SPARC" (GRIPS) initiative. A set of middle atmosphere-climate models has been compared, identifying common biases. Mechanisms for these biases are being studied in some detail, including off-line assessments of the radiation transfer codes and coordinated studies of the impacts of gravity wave drag due to sub-grid-scale processes. ensemble of models will be presented, along with numerical experiments undertaken with one or more models, designed to investigate the mechanisms at work in the atmosphere. The discussion will focus on dynamical and radiative mechanisms in the current climate, but implications for coupled ozone chemistry and the future climate will be assessed.
This project advanced our understanding of the processes governing cumulus cloud formation and provides an improved observational basis for validating earth system models. To be specific, the project used laser- and radar-remote sensors to examine the physical properties of updrafts that rise from earth’s surface and initiate clouds deeper in the atmosphere. These updrafts comprise “thermals” and “plumes” and occur at small spatial and temporal scales (e.g., 10s of minutes, 100s of meters). These small scales preclude explicit representation in most earth system and climate models, and thus necessitate “sub-grid-scale” parameterization of updraft processes. The innovation of this project was to directly measure the size, shape, strength and water vapor content of these updrafts with Doppler and Raman lidars, respectively, and to link these updraft properties to cloud processes using vertically pointed weather radars. The resulting data sets comprise 100s of thousands of updrafts and thousands of clouds, which far exceeds previous efforts, and thereby provides a robust statistical and physical representation of these processes. From these large datasets the project produced a sequence of scientific analyses that: (1) Elucidate how variations in the turbulent structure of the convective boundary layer control shallow cumulus convection, (2) Quantify the upward transport of water vapor to cloud base via thermals and plumes, (3) Validate large-eddy simulations of updrafts and shallow convective clouds, (4) Demonstrate a size-to-strength relationship between updraft width and updraft speed, and (5) Demonstrate how updrafts interact with the stability at the top of the convective boundary layer to modulate the depth and vigor of convective clouds. These results have been disseminated via several published journal articles, academic theses, and conference presentations. Collectively these results contribute to the Atmospheric System Research (ASR) program’s goal to “improve understanding of the key cloud, aerosol, precipitation, and radiation processes that affect the Earth’s radiative balance and hydrological cycle, particularly processes that limit the predictive ability of regional and global earth system models”.
Several experiments are described in which the sub-grid-scale vertical eddy viscosity in the GISS global general circulation model was varied. The results show that large viscosities suppress large-scale eddies in middle and high latitudes, but enhance the circulation in the tropical Hadley cell and increase the extent of the tropical easterlies. Comparison with observations shows that the GISS model requires eddy viscosities about 1 sq m per sec or less to give realistic results for middle and high latitudes, and eddy viscosities about 100 sq m per sec to give realistic results for low latitudes. A plausible mechanism for the implied increase in small-scale mixing in low latitudes is cumulus convection.
NASA's GMAO produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time 5-day forecasts of atmospheric chemical composition at unprecedented horizontal resolution of 0.25 degrees (~25 km), freely available to the public. This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). GEOS-CF also assimilated satellite observations into the system for improved representation of weather and smoke. The assimilation system is currently being expanded to include chemically reactive trace gases. While the main focus of this new product is on tropospheric air quality information, the GEOS-Chem chemistry model used in this system includes the unified tropospheric stratospheric chemistry mechanism for improved forecasts of total column ozone during anomalous dynamical and chemical events. I will discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions for GEOS-CF and reanalysis products. In addition, I will show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.
NASA's Global Modeling and Assimilation Office (GMAO) produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time5-day forecasts of atmospheric composition at unprecedented horizontal resolution of 0.25 degrees (~25 km). This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module (version 12) to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). Satellite observations are assimilated into the system for improved representation of weather and smoke. The assimilation system is being expanded to include chemically reactive trace gases. We discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions, notably incorporating new observations (TROPOMI, geostationary satellites) and machine learning techniques. We show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.
Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.
We report turbulence intensity (TI) is often used to quantify the strength of turbulence in wind energy applications and serves as the basis of standards in wind turbine design. Thus, accurately characterizing the spatiotemporal variability in TI should lead to improved predictions of power production. Nevertheless, turbulence measurements over the ocean are far less prevalent than over land due to challenges in instrumental deployment, maintenance, and operation. Atmospheric models such as mesoscale (weather prediction) and large-eddy simulation (LES) models are commonly used in the wind energy industry to assess the spatial variability of a given site. However, the TI derivation from atmospheric models has not been well examined. An algorithm is proposed in this study to realize online calculation of TI in the Weather Research and Forecasting (WRF) model. Simulated TI is divided into two components depending on scale, including sub-grid (parameterized based on turbulence kinetic energy (TKE)) and grid resolved. The sensitivity of sea surface temperature (SST) on simulated TI is also tested. An assessment is performed by using observations collected during a field campaign conducted from February to June 2020 near the Woods Hole Oceanographic Institution Martha's Vineyard Coastal Observatory. Results show that while simulated TKE is generally smaller than the lidar-observed value, wind speed bias is usually small. Overall, this leads to a slight underestimation in sub-grid-scale estimated TI. Improved SST representation subsequently reduces model biases in atmospheric stability as well as wind speed and sub-grid TI near the hub height. Large TI events in conjunction with mesoscale weather systems observed during the studied period pose a challenge to accurately estimating TI from models. Due to notable uncertainty in accurately simulating those events, this suggests summing up sub-grid and resolved TI may not be an ideal solution. Efforts in further improving skills in simulating mesoscale flow and cloud systems are necessary as the next steps.
Earth system models (ESMs) are a common tool for estimating local and global greenhouse gas emissions under current and projected future conditions. Efforts are underway to expand the representation of wetlands in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) by resolving the simultaneous contributions to greenhouse gas fluxes from multiple, different, sub-grid-scale patch-types, representing different eco-hydrological patches within a wetland. However, for this effort to be effective, it should be coupled with the detection and mapping of within-wetland eco-hydrological patches in real-world wetlands, providing models with corresponding information about vegetation cover. In this short communication, we describe the application of a recently developed NDVI-based method for within-wetland vegetation classification on a coastal wetland in Louisiana and the use of the resulting yearly vegetation cover as input for ELM simulations. Processed Harmonized Landsat and Sentinel-2 (HLS) datasets were used to drive the sub-grid composition of simulated wetland vegetation each year, thus tracking the spatial heterogeneity of wetlands at sufficient spatial and temporal resolutions and providing necessary input for improving the estimation of methane emissions from wetlands. Our results show that including NDVI-based classification in an ELM reduced the uncertainty in predicted methane flux by decreasing the model’s RMSE when compared to Eddy Covariance measurements, while a minimal bias was introduced due to the resampling technique involved in processing HLS data. Our study shows promising results in integrating the remote sensing-based classification of within-wetland vegetation cover into earth system models, while improving their performances toward more accurate predictions of important greenhouse gas emissions.
The fluvial system is a major concern in modeling landform evolution in response to tectonic deformation. Three stream bed types (bedrock, coarse-bed alluvial, and fine-bed alluvial) differ in factors controlling their occurrence and evolution and in appropriate modeling approaches. Spatial and temporal transitions among bed types occur in response to changes in sediment characteristics and tectonic deformation. Erosion in bedrock channels depends upon the ability to scour or pluck bed material; this detachment capacity is often a power function of drainage area and gradient. Exposure of bedrock in channel beds, due to rapid downcutting or resistant rock, slows the response of headwater catchments to downstream baselevel changes. Sediment routing through alluvial channels must account for supply from slope erosion, transport rates, abrasion, and sorting. In regional landform modeling, implicit rate laws must be developed for sediment production from erosion of sub-grid-scale slopes and small channels.
In GCM-physics evaluations with the currently available ARM-CART SCM datasets, McRAS produced very similar character of near surface errors of simulated temperature and humidity containing typically warm and moist biases near the surface and cold and dry biases aloft. We argued it must have a common cause presumably rooted in the model physics. Lack of vertical adjustment of horizontal transport was thought to be a plausible source. Clearly, debarring such a freedom would force the incoming air to diffuse into the grid-cell which would naturally bias the surface air to become warm and moist while the upper air becomes cold and dry, a characteristic feature of McRAS biases. Since, the errors were significantly larger in the two winter cases that contain potentially more intense episodes of cold and warm advective transports, it further reaffirmed our argument and provided additional motivation to introduce the corrections. When the horizontal advective transports were suitably modified to allow rising and/or sinking following isentropic pathways of subgrid scale motions, the outcome was to cool and dry (or warm and moisten) the lower (or upper) levels. Ever, crude approximations invoking such a correction reduced the temperature and humidity biases considerably. The tests were performed on all the available ARM-CART SCM cases with consistent outcome. With the isentropic corrections implemented through two different numerical approximations, virtually similar benefits were derived further confirming the robustness of our inferences. These results suggest the need for insentropic advective transport adjustment in a GCM due to subgrid scale motions.
The Land and Land Ice Theme in the Coupled Model Intercomparison Project Phase 7 (CMIP7) represents the current understanding of physical processes in land surface ecosystems, hydrology, cryosphere, and their physical interactions with other Earth system components. Simulations from Earth system models (ESMs) could provide crucial information for assessing planetary safety, such as critical tipping elements, and be used to inform climate risks for improving climate impact assessments and policy decisions. This paper presents a collaborative effort to identify scientific opportunities in the Land and Land Ice Theme of the CMIP7 Data Request. The proposed opportunities build upon advances in ESMs, including new freshwater system and land ice processes being included in CMIP7, as well as the scientific community's demand for high-frequency and sub-grid-scale land surface outputs. In total, 25 variable groups that contain 716 variables have been identified to be potentially available to the broad scientific audience for performing analysis in land–atmosphere coupling, hydrological processes and freshwater systems, glacier and ice sheet mass balance and their influence on the sea levels, land use, and plant phenology. Key reflections from this data request effort include advocacy for closer engagement between the user community and modeling groups, reduction in the technical barriers to tracking existing parameters and defining new variables, and more streamlined variable management. These will be essential to enhance the usability and reliability of CMIP7 outputs for climate and Earth system research and applications to a broad audience that relies on the CMIP7 endeavor.
Sudden Stratospheric Warmings (SSWs) and Polar Stratospheric Clouds (PSCs) exhibit striking inter-hemispheric asymmetries: SSWs are frequent in the Arctic but rare in the Antarctic, while PSCs are more persistent in the Antarctic. Although land–sea thermal contrast and orography (LSCO) have been traditionally invoked to explain these asymmetries, here we show using semi-idealized model experiments that even though LSCO strongly impacts the mean state of SSWs and PSCs, it alone cannot fully account for the observed differences. Using model experiments, we reveal that midlatitude oceanic sea surface temperature (SST) fronts represent a crucial supporting additional driver of the hemispheric stratospheric differences. Like LSCO, SST fronts enhance stratospheric convergence of resolved waves, strengthening the Brewer–Dobson Circulation and inducing high-latitude adiabatic warming. This warming significantly enhances Arctic SSW frequency and strongly suppresses PSC formation. Sub-grid-scale non–resolved wave forcing modulates the stratospheric resolved waves effect. The oceanic impact is dominated by North Pacific SST fronts. Our results highlight the indispensable role of SST fronts in shaping Arctic–Antarctic asymmetries in stratospheric dynamics and associated extremes.
Tropical tropospheric ozone fields from two different experiments performed with an off-line ozone assimilation system developed in NASA's Data Assimilation Office (DAO) are examined. Assimilated ozone fields from the two experiments are compared with the collocated ozone profiles from the Southern Hemispheric Additional Ozonesondes (SHADOZ) network. Results are presented for 1998. The ozone assimilation system includes a chemistry-transport model, which uses analyzed winds from the Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). The two experiments use wind fields from different versions of GEOS DAS: an operational version of the GEOS-2 system and a prototype of the GEOS-4 system. While both versions of the DAS utilize the Physical-space Statistical Analysis System and use comparable observations, they use entirely different general circulation models and data insertion techniques. The shape of the annual-mean vertical profile of the assimilated ozone fields is sensitive to the meteorological analyses, with the GEOS-4-based ozone being closest to the observations. This indicates that the resolved transport in GEOS-4 is more realistic than in GEOS-2. Remaining uncertainties include quantification of the representation of sub-grid-scale processes in the transport calculations, which plays an important role in the locations and seasons where convection dominates the transport.
Use of the Smagorinsky eddy-viscosity formulation and related schemes for subgrid-scale parameterization of large eddy simulation models requires specification of a single length scale, earlier related by Lilly to the scale of filtering and/or numerical resolution. An anisotropic integration of the Kolmogoroff enstrophy spectrum allows generalization of that relationship to anisotropic resolution. It is found that the Deardorff assumption is reasonably accurate for small anisotropies and can be simply improved for larger values.
Sub-grid-scale processes occurring at or near the surface of an ice sheet have a potentially large impact on local and integrated net accumulation of snow via redistribution and sublimation. Given observational complexity, they are either ignored or parameterized over large-length scales. Here, we train random forest (RF) models to predict variability in net accumulation over the Antarctic Ice Sheet using atmospheric variables and topographic characteristics as predictors at 1 km resolution. Observations of net snow accumulation from both in situ and airborne radar data provide the input observable targets needed to train the RF models. We find that local net accumulation deviates by as much as 172% of the atmospheric model mean. The correlation in space between the predicted net accumulation variability and satellite-derived surface-height change indicates that surface processes operate differently through time, driven largely by the seasonal anomalies in snow accumulation.