Engineering Papers⌕ Search

Engineering topics

Sakaguchi, Koichi

Publications and source records attributed to Sakaguchi, Koichi.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Root Dynamics Mitigate Warm and Dry Biases over the Central United States

The central United States frequently exhibits warm and dry biases in simulations of summertime conditions, a persistent feature that remains unresolved. While previous studies linked these biases to misrepresented surface energy exchanges, the role of belowground processes remains poorly understood. Here, we demonstrate that inadequate representation of root water uptake in land surface models contributes to this bias. Using both offline Noah-MP and coupled WRF-Noah-MP simulations with static and dynamic root water uptake schemes, we show that the inclusion of dynamic root processes reduces 2-m air temperature biases and enhances precipitation, primarily by increasing the rain rate of convective systems. Offline and coupled simulations further reveal that the cooling effects and precipitation increases are amplified through positive land-atmosphere feedback, active only in the coupled model. These findings highlight an important role of root in modulating land-atmosphere interactions and underscore the need to refine root-zone processes to improve regional atmospheric simulations.

Yang, Zhao (ORCID:0000000288027130)↗
DOE OSTI2024

The first ensemble of kilometer-scale simulations of a hydrological year over the third pole

An accurate understanding of the current and future water cycle over the Third Pole is of great societal importance, given the role this region plays as a water tower for densely populated areas downstream. An emerging and promising approach for skillful climate assessments over regions of complex terrain is kilometer-scale climate modeling. As a foundational step towards such simulations over the Third Pole, we present a multi-model and multi-physics ensemble of kilometer-scale regional simulations for the hydrological year of October 2019 to September 2020. The ensemble consists of 13 simulations performed by an international consortium of 10 research groups, configured with a horizontal grid spacing ranging from 2.2 to 4 km covering all of the Third Pole region. These simulations are driven by ERA5 and are part of a Coordinated Regional Climate Downscaling EXperiment Flagship Pilot Study on Convection-Permitting Third Pole. The simulations are compared against available gridded and in-situ observations and remote-sensing data, to assess the performance and spread of the model ensemble compared to the driving reanalysis during the cold and warm seasons. Although ensemble evaluation is hindered by large differences between the gridded precipitation datasets used as a reference over this region, we show that the ensemble improves on many warm-season precipitation metrics compared with ERA5, including most wet-day and hour statistics, and also adds value in the representation of wet spells in both seasons. As such, the ensemble will provide an invaluable resource for future improvements in the process understanding of the hydroclimate of this remote but important region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

pinacles

Predicting Interactions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) is a high-performance, massively parallel, anelastic atmospheric model written in Python

Pressel, Kyle↗
DOE OSTI2022

A Machine‐Learning‐Assisted Stochastic Cloud Population Model as a Parameterization of Cumulus Convection

Abstract A machine‐learning‐assisted stochastic cloud population model is coupled with the Advanced Research Weather Research and Forecasting (WRF) model to represent fluctuations in the cloud‐base mass flux associated with the life cycles and interactions among cumulus convection cells. In this cloud population model, the size distribution and the associated cloud‐base mass flux of the convective cells are related to their previous state and to the change in the total convective area via a transition function. The convective area tendency in turn is assumed to depend on the cloud‐base mass flux that is resolved by the host WRF model. The transition function is represented by a single hidden‐layer neural network trained by the evolution of convective cell size distributions in a 1‐km grid‐spacing WRF simulation run over the Australian Monsoon region. At every grid point of the host model, the cloud population model predicts the cell size and cloud‐base mass flux distributions from which a random sample of cells is fed to an entraining parcel model that calculates precipitation as well as the associated liquid water potential temperature and total moisture tendencies. These tendencies are averaged over the cells and provided to the host model. Several regional simulations are performed over tropical and midlatitude domains to test this as a potential approach to scale‐aware parameterization. It is shown that such an approach could be a new promising path to simulating realistic precipitation statistics and propagation of precipitation associated with the Madden‐Julian Oscillation while maintaining realistic depictions of the diurnal cycle over both land and ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effect of the Grell-Freitas Deep Convection Scheme in Quasi-Uniform and Variable-Resolution Aquaplanet CAM Simulations

The Grell-Freitas deep convection scheme (GF) has been implemented in the Community Atmosphere Model (CAM) with the nonhydrostatic Model for Prediction Across Scales (MPAS) dynamical core. Aquaplanet simulations have been performed with a global quasi-uniform resolution mesh with ~120 km grid spacing and a variable resolution mesh employing a circular refined region with ~30 km grid spacing centering at the equator and ~120 km elsewhere. GF produces overall patterns of the general circulations comparable to those of other aquaplanet experiments, including those produced using the default Zhang-McFarlane convection scheme (ZM) in CAM Version 5.4. GF alleviates the overestimated occurrence of weak intensity precipitation characteristic of ZM due to GF’s different closure methods and deep-convection trigger conditions. Compared to ZM, the GF characteristics of the equatorial Kelvin waves are closer to observations. GF also better simulates precipitation at variable resolution, with no unrealistic wave distortion by the refinement. Overall, GF shows less sensitivity of convective precipitation and heating/moistening tendencies to changes in horizontal resolution compared to ZM. However, vertically-overlapped total cloud fraction simulated with GF in the tropics is more sensitive to the resolution than ZM due to the larger resolution-sensitivity of high-level cloud fraction in the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Conservation of Dry Air, Water, and Energy in CAM and Its Potential Impact on Tropical Rainfall

Abstract For the Community Atmosphere Model version 6 (CAM6), an adjustment is needed to conserve dry air mass. This adjustment exposes an inconsistency in how CAM6’s energy budget incorporates water—in CAM6 water in the vapor phase has energy, but condensed phases of water do not. When water vapor condenses, only its latent energy is retained in the model, while its remaining internal, potential, and kinetic energy are lost. A global fixer is used in the default CAM6 model to maintain global energy conservation, but locally the energy tendency associated with water changing phase violates the divergence theorem. This error in energy tendency is intrinsically tied to the water vapor tendency, and reaches its highest values in regions of heavy rainfall, where the error can be as high as 40 W m −2 annually averaged. Several possible changes are outlined within this manuscript that would allow CAM6 to satisfy the divergence theorem locally. These fall into one of two categories: 1) modifying the surface flux to balance the local atmospheric energy tendency and 2) modifying the local atmospheric tendency to balance the surface plus top-of-atmosphere energy fluxes. To gauge which aspects of the simulated climate are most sensitive to this error, the simplest possible change—where condensed water still does not carry energy and a local energy fixer is used in place of the global one—is implemented within CAM6. Comparing this experiment with the default configuration of CAM6 reveals precipitation, particularly its variability, to be highly sensitive to the energy budget formulation. Significance Statement This study examines and explains spurious regional sources and sinks of energy in a widely used climate model. These energy errors result from not tracking energy associated with water after it transitions from the vapor phase to either liquid or ice. Instead, the model used a global fixer to offset the energy tendency related to the energy sources and sinks associated with condensed water species. We replace this global fixer with a local one to examine the model sensitivity to the regional energy error and find a large sensitivity in the simulated hydrologic cycle. This work suggests that the underlying thermodynamic assumptions in the model should be revisited to build confidence in the model-simulated regional-scale water and energy cycles.

Harrop, Bryce E.↗
DOE OSTI2022

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multiscale Simulation of Precipitation Over East Asia by Variable Resolution CAM-MPAS

This study evaluates the precipitation over East Asia simulated by the variable resolution CAM-MPAS featuring the Model for Prediction Across Scales (MPAS) atmospheric dynamical core coupled with the physics parameterizations of the Community Atmosphere Model (CAM) version 5.4. Two CAM-MPAS experiments, one with a global quasi-uniform resolution mesh of 120 km (MPAS-UR) and the other with a variable resolution mesh of 30 to 120 km refined over East Asia (MPAS-VR), are conducted from 1989 to 2005 following the Atmospheric Model Intercomparison Project protocol. Two regional climate model simulations at ~25 km resolution from the Coordinated Regional Downscaling Experiment East Asia second plan are also analyzed for comparison. Results show that CAM-MPAS performs better than the selected regional models in simulating the precipitation climatology over East Asia. Compared with MPAS-UR, MPAS-VR with refinement over East Asia better simulates the precipitation over and around the Tibetan Plateau (TP), the frequency distribution of light, moderate, and heavy precipitation in various regions of China, and the interannual precipitation pattern associated with the East Asian summer monsoon. The improvements in MPAS-VR relative to MPAS-UR largely come from the decreased grid spacing rather than the use of finer-scale terrain information. Particularly, decreased grid spacing is more essential for simulating the precipitation magnitude over TP and the location of peak precipitation south of TP, likely due to the better resolved physical and dynamical processes associated with orographic precipitation. The 30-120 km CAM-MPAS model is shown to be a promising tool for precipitation simulation over East Asia.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2021

Crucial roles of eastward propagating environments in the summer MCS initiation over the U.S. Great Plains

Mesoscale convective systems (MCSs) frequently occur over the U.S. Great Plains during summer. An analysis using self-organizing map is conducted to identify four types of summer MCS initiation environments during 2004-2017. The first two types feature favorable large-scale environments at both upper and low levels, while Type-3 has favorable lower-level and surface conditions but unfavorable upper-level circulation, and Type-4 features the most unfavorable large-scale environments for MCS initiation. Despite the unfavorable large-scale environment, the convection-centered environments in Type-4 are favorable for MCS initiation and similar to the first two types, suggesting a role of sub-synoptic disturbances as an MCS precursor. All four types of the MCS initiation delineate a clear eastward propagating feature in many fields, such as upper-level potential vorticity/geopotential height, surface pressure and surface equivalent potential temperature, upstream up to 25°-longitude west of and ~36 hours before the MCS initiation. The propagating environments and local, non-propagating low-level moisture are found to be important in MCS initiation at the foothill of the Rocky Mountains, but over the central Great Plains, it is the coupling of dynamical and moisture anomalies associated with propagating waves that results in the MCS initiation. By tracking MCSs and mid-tropospheric perturbations (MPs), a type of sub-synoptic disturbances with Rocky Mountains origin, ~30% of MPs is associated with MCS initiation, mostly in Type-4. Although MPs are related to a small fraction of MCS initiation, MCSs that are associated with MPs tend to produce more rainfall in a larger area with a stronger convective intensity, suggesting MPs to be a source of intense MCSs in summer.

Song, Fengfei↗
DOE OSTI2021

Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study

As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Climate change impacts on wind power generation

Wind energy is a virtually carbon-free and pollution-free electricity source, with global wind resources greatly exceeding electricity demand. Accordingly, the installed capacity of wind turbines grew at an annualized rate of >20% from 2000 to 2019 and is projected to increase by a further 50% by the end of 2023. Here, in this Review, we describe the factors that dictate the wind resource magnitude and variability and illustrate the tools and techniques that are being used to make projections of wind resources and wind turbine operating conditions. Natural variability due to the action of internal climate modes appears to dominate over global-warming-induced non-stationarity over most areas with large wind energy installations or potential. However, there is evidence for increased wind energy resources by the end of the current century in northern Europe and the US Southern Great Plains. New technology trends are changing the sensitivity of wind energy to global climate non-stationarity and, thus, present new challenges and opportunities for innovative research. The evolution of climate modelling to increasingly address mesoscale processes is providing improved projections of both wind resources and wind turbine operating conditions, and will contribute to continued reductions in the levelized cost of energy from wind power generation.

17 WIND ENERGY↗
DOE OSTI2020

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

FACETS Dynamical Downscaling Simulations over North America by the CAM-MPAS Variable-Resolution Model

Comprehensive assessment of climate datasets created by statistical or dynamical models is important for effectively communicating model projection and associated uncertainty to stakeholders and decision-makers. The Department of Energy FACETS project aims to foster such communication through development of metrics and their demonstration on a hierarchy of downscaled climate datasets to quantify aspects of climate change projections that are credible, particularly for supporting decisions related to the energy-water-land nexus. As a part of this effort, we have produced a regional climate dataset using the Model for Prediction Across Scales coupled to the Community Atmosphere Model (CAM-MPAS). This global modeling framework is configured with variable-resolution meshes featuring higher resolutions over North America, as well as quasi-uniform resolution meshes across the globe. The variable-resolution configurations allow fine-scale features to be better resolved inside the refinement and interact with global-scale circulations. The dataset includes multiple uniform- (240km and 120km) and variable-resolution (200-50km, 100-25km, and 46-12km) simulations that are designed to be compatible with other regional climate simulations that contribute to the hierarchy of downscaled climate datasets of the project. Furthermore, the dataset consists of simulations for both the present-day (1989-2010) and future (2079-2100) climate and post-processing of the model output has been coordinated across the project for consistency to facilitate common analysis across the hierarchy of datasets. Altogether, this CAM-MPAS model dataset provides a unique opportunity to assess the influence of resolutions and modeling framework on model credibility and climate change projection.

Sakaguchi, Koichi↗
DOE OSTI2020

Understanding Irrigation Impacts on Low-Level Jets over the Great Plains

Low-level jets (LLJs) over the Great Plains are well-known for their importance in transporting moisture from the Gulf of Mexico and forming precipitation in the Central US. However, the impact of irrigation practice on LLJs is rarely studied. To understand the irrigation impact on the Great Plains LLJs, the Weather Research and Forecasting (WRF) model is employed with and without an irrigation scheme. The simulated wind field is evaluated against the observed wind from Midlatitude Continental Convective Clouds Experiment (MC3E). The results show that adding irrigation improves the simulated wind, especially at lower altitudes below 800 hPa. Irrigation reduces the LLJ frequency over the Great Plains, mainly by reducing the intense LLJ frequency. The reduced LLJ frequency is caused by the irrigation-induced weakening of meridional wind. The favorable conditions for LLJs have been identified by examining the mechanism associated with LLJs: a ridge-like pattern with weak and calm synoptic-scale conditions in the mid-levels and a strong west-east temperature gradient at lower levels that promote meridional winds. Irrigation-induced cooling leads to a baroclinic pattern with high pressure at lower altitudes over the irrigation region and low pressure in the mid- and upper levels in the downstream region, resulting in a northerly component that hinders the formation of LLJs. Low-level surface cooling also decreases the west-east temperature gradient and hence weakens the thermal wind forcing. The irrigation-induced mid-level low and weaker west-east temperature gradient are responsible for the reduced meridional winds and the reduced LLJ frequencies.

Yang, Zhao↗
DOE OSTI2020

The ongoing need for high-resolution regional climate models: Process understanding and stakeholder information

Regional climate modeling addresses our need to understand and simulate climatic processes and phenomena unresolved in global models. High resolution models are generally more skillful in simulating extremes, such as heavy precipitation, strong winds, and severe storms. In addition, research has shown that fine-scale features such as mountains, coastlines, lakes, irrigation, land use, and urban heat islands can substantially influence a region’s climate and its response to changing forcings. Regional climate simulations explicitly simulating convection are now being performed, providing an opportunity to illuminate new physical behavior that previously was represented by parameterizations with large uncertainties. Regional and global models are both advancing toward higher resolution, as computational capacity increases. However, the resolution and ensemble size necessary to produce a sufficient statistical sample of these processes in global models has proven too costly for contemporary supercomputing systems. Regional climate models are thus indispensable tools that complement global models for understanding regional climate variability and change, and are critical for supporting societal responses to changing climate.

Gutowski, William↗
DOE OSTI2019

The Impact of Variable Land-Atmosphere Coupling on Convective Cloud Populations Observed During the 2016 HI-SCALE Field Campaign

We use a Large Eddy Simulation (LES) model with a nested domain configuration (297 and 120 km wide) and an interactive land surface parameterization to simulate the complex population of shallow clouds observed on August 30, 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) campaign conducted in north-central Oklahoma. Shallow convective clouds first formed over southeast Oklahoma and then spread towards the northwest into southern Kansas. By the early afternoon, the relatively uniform shallow cloud field became more complex in which some regions became nearly cloud free and in other regions larger shallow clouds developed with some transitioning to deeper, precipitating convection. We show that the model reproduces the observed heterogeneity in the cloud populations only when realistic variations in soil moisture are used to initialize the model. While more variable soil moisture and to a lesser extend cool lake temperatures drive the initial spatial heterogeneity in the cloud populations, precipitation-driven cold pools become an important factor after 1300 CST. When smoother soil moisture variations are used in the model, more uniform shallow cloud populations are predicted with far fewer clouds that transition to deeper, precipitating convection that produce cold pools. An algorithm that tracks thousands of individual cumulus show that the more realistic soil moisture distributions produces clouds that are larger and have a longer lifetime. The results suggest that shallow and deep convection parameterizations used by mesoscale models need to account for the effects of variable land-atmosphere interactions and cold pools.

Fast, Jerome D.↗