SEARCH · Engineering Papers

Results for “Variability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Data for Greenhouse Gas Accounting Procedures in Low Carbon Fuel Policies Overlook the Spatial Variability of Miscanthus-Derived Sustainable Aviation Fuel

Low carbon fuel policies such as the U.S. Renewable Fuel Standard (RFS), Canada Clean Fuel Regulations (CFR), and California Low Carbon Fuel Standard (LCFS) as well as the 45Z tax credit are intended to reduce greenhouse gas (GHG) emissions from transportation. Cellulosic feedstocks, optimized biorefineries, and favorable farming locations can significantly reduce biofuel carbon intensity (CI). Despite advances in field-to-fuel GHG monitoring and flexibility in resource allocation within biorefineries (e.g., governing net electricity production), rigid CI accounting procedures in current policies may limit CI responsiveness across candidate sites and processing facilities. This work examines a hypothetical biomass-to-sustainable aviation fuel (SAF) pathway using miscanthus and alcohol-to-jet (i) to demonstrate how GHG accounting requirements drive estimates of biofuel CIs and (ii) to explore potential CI and financial implications of scenario-specific life cycle assessment (LCA). Results demonstrate that GHG accounting using the CFR/LCFS can reasonably account for distinct levels of net electricity production by a biorefinery, but only the CFR yields similar CI sensitivity to spatially explicit factors (feedstock CI, grid electricity CI) as scenario-specific LCA: most GHG accounting frameworks do not capture CI variation across candidate sites in the United States. Ultimately, this work demonstrates the importance of LCA methodological specifications in low carbon fuel policies and tax credits.

Miscanthus↗

DOE OSTI2020

Consistency-Enhanced Evolution for Variable Selection Can Identify Key Chemical Information from Spectroscopic Data

In the last few decades, spectroscopic techniques such as near-infrared (NIR) spectroscopy have gained wide applications in several industries, such as the pharmaceutical, agricultural, oil, and gas industries. As a result, various soft sensors have been developed to predict sample properties from spectroscopic readings. Because the spectroscopic readings at different wavelengths, especially at the adjacent wavelengths, are highly correlated, it has been shown that variable selection could significantly improve a soft sensor’s prediction performance while reducing the model complexity. To improve the prediction performance, most variable selection methods focus on identifying the variables (i.e., wavelengths or wavelength segments) that are strongly correlated with the dependent variable. Although many successful applications have been reported, these variable selection methods do have their limitations. Specifically, the selected wavelengths sometimes show little connection to the chemical bounds or functional groups presenting in the sample. In addition, the selected variables can be quite sensitive to the choice of the training samples. In this work, we address these limitations from a different perspective: if a variable selection algorithm can identify the truly relevant input variables, it should consistently identify the same subset of variables regardless of the choice of the training samples. Therefore, we propose a variable selection method that aims to improve the consistency of variable selection resulting from different training samples. Furthermore, the new algorithm is termed consistency-enhanced evolution for variable selection (CEEVS). To demonstrate the performance and robustness of CEEVS, we compare the proposed method with three representative variable selection methods using five published NIR data sets. These case studies clearly demonstrate that by improving the variable selection consistency, we can not only achieve improved prediction performance, but also identify key chemical information from spectroscopic data.

42 ENGINEERING↗

DOE OSTI2022

Benchmarking Simulated Precipitation Variability Amplitude across Time Scales

Here, objective performance metrics that measure precipitation variability across time scales from subdaily to interannual are presented and applied to Historical simulations of Coupled Model Intercomparison Project phase 5 and 6 (CMIP5 and CMIP6) models. Three satellite-based precipitation estimates (IMERG, TRMM, and CMORPH) are used as reference data. We apply two independent methods to estimate temporal variability of precipitation and compare the consistency in their results. The first method is derived from power spectra analysis of 3-hourly precipitation, measuring forced variability by solar insolation (diurnal and annual cycles) and internal variability at different time scales (subdaily, synoptic, subseasonal, seasonal, and interannual). The second method is based on time averaging and facilitates estimating the seasonality of subdaily variability. Supporting the robustness of our metric, we find a near equivalence between the results obtained from the two methods when examining simulated-to-observed ratios over large domains (global, tropics, extratropics, land, or ocean). Additionally, we demonstrate that our model evaluation is not very sensitive to the discrepancies between observations. Our results reveal that CMIP5 and CMIP6 models in general overestimate the forced variability while they underestimate the internal variability, especially in the tropical ocean and higher-frequency variability. The underestimation of subdaily variability is consistent across different seasons. The internal variability is overall improved in CMIP6, but remains underestimated, and there is little evidence of improvement in forced variability. Increased horizontal resolution results in some improvement of internal variability at subdaily and synoptic time scales, but not at longer time scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Variability and Diversity Load Model Tool [SWR-20-03]

The motivation for the development of this tool and the underlying algorithms and methods was to enable the development of high-temporal resolution, realistic time-series data for quasi-static time-series (QSTS) analysis of distribution systems. Often, aggregated load profile data for a distribution circuit is available (e.g. feeder loading data collected via SCADA at the utility substation) and, while this data is typically accurate it masks the considerable variability of the 100’s or 1000’s of individual loads connected on the circuit. This tool was developed to model both the increased variability expected for these individual loads (e.g. the load of a single distribution transformer connected to 8-12 houses) and the expected diversity between loads on the circuit. It is important to note that the difference in variability and diversity, in the context of this tool, is that variability modeling only adds representative variability due to disaggregated load characteristics (e.g. the presence in the load profile of loads turning off and on like an air conditioner/oven) while the average energy profile remains the same as the user supplied power profile. Diversity modeling generates multiple individual load profiles which, in aggregate, sum to the user supplied power profile. Diversity is effectively variability in the energy usage over longer periods of time than seen in the variability model. Put another way, variability modeling supplies the expected variability due to the operation of various end-use loads and diversity modeling supplies the usage differences due to human behavior, schedules, etc. This load modeling tool was developed for use in generating data for distribution systems. Modeling is summarized by two major functions: 1) taking low resolution load profiles and adding intra-seconds variability onto the profiles, and 2) taking a user supplied load profile and distribution factors and adding both diversity and variability to the user supplied profile.

Zhu, Xiangqi↗

DOE OSTI2023

Quantifying spatiotemporal variability in occupant exposure to an indoor airborne contaminant with an uncertain source location

Well-mixed zone models are often employed to compute indoor air quality and occupant exposures. While effective, a potential downside to assuming instantaneous, perfect mixing is underpredicting exposures to high intermittent concentrations within a room. When such cases are of concern, more spatially resolved models, like computational-fluid dynamics methods, are used for some or all of the zones. But, these models have higher computational costs and require more input information. A preferred compromise would be to continue with a multi-zone modeling approach for all rooms, but with a better assessment of the spatial variability within a room. Here to do so, we present a quantitative method for estimating a room's spatiotemporal variability, based on influential room parameters. Our proposed method disaggregates variability into the variability in a room's average concentration, and the spatial variability within the room relative to that average. This enables a detailed assessment of how variability in particular room parameters impacts the uncertain occupant exposures. To demonstrate the utility of this method, we simulate contaminant dispersion for a variety of possible source locations. We compute breathing-zone exposure during the releasing (source is active) and decaying (source is removed) periods. Using CFD methods, we found after a 30 minutes release the average standard deviation in the spatial distribution of exposure was approximately 28% of the source average exposure, whereas variability in the different average exposures was lower, only 10% of the total average. We also find that although uncertainty in the source location leads to variability in the average magnitude of transient exposure, it does not have a particularly large influence on the spatial distribution during the decaying period, or on the average contaminant removal rate. By systematically characterizing a room's average concentration, its variability, and the spatial variability within the room important insights can be gained as to how much uncertainty is introduced into occupant exposure predictions by assuming a uniform in-room contaminant concentration. We discuss how the results of these characterizations can improve our understanding of the uncertainty in occupant exposures relative to well-mixed models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Long-term spatial and temporal solar resource variability over America using the NSRDB version 3 (1998–2017)

The study assesses the long-term spatial and temporal solar resource variability in America using the 20-year National Renewable Energy Laboratory's (NREL's) National Solar Radiation Database (NSRDB). Specifically, the coefficient of variation (COV) is used to analyze the spatial and temporal (interannual and seasonal) variability. Further, both spatial and temporal long-term variability are analyzed using the Köppen-Geiger climate classification. The temporal variability is found that, on average, the continental United States (CONUS) COV reaches up to 5% for global horizontal irradiance (GHI) and 10% for direct normal irradiance (DNI), and that the NSRDB domain's COV is roughly twice that of CONUS. For the seasonal variability analysis, the winter months are found to exhibit higher COV than the other seasons. In particular, December exhibits the highest variability, reaching on average 30% for DNI and 20% for GHI over various areas. On the other hand, the summer months demonstrate significantly lower variability, reaching only less than 20% for DNI and 10% for GHI, on average. Similarly, the spatial variability is analyzed by comparing each pixel to its neighbors. The long-term spatial variability is found to increase with the number of neighboring pixels being considered, which is equivalent to an increase in distance (within a 100-km x 100-km square grid). As expected, the DNI spatial variability is higher than that of GHI. Moreover, the annual solar irradiance anomalies are found to reach ±25% for both GHI and DNI (and even exceed those value in some instances) during each year of the 20-year period.

14 SOLAR ENERGY↗

DOE OSTI2022

Descriptor Aided Bayesian Optimization for Many-Level Qualitative Variables With Materials Design Applications

Abstract Engineering design often involves qualitative and quantitative design variables, which requires systematic methods for the exploration of these mixed-variable design spaces. Expensive simulation techniques, such as those required to evaluate optimization objectives in materials design applications, constitute the main portion of the cost of the design process and underline the need for efficient search strategies—Bayesian optimization (BO) being one of the most widely adopted. Although recent developments in mixed-variable Bayesian optimization have shown promise, the effects of dimensionality of qualitative variables have not been well studied. High-dimensional qualitative variables, i.e., with many levels, impose a large design cost as they typically require a larger dataset to quantify the effect of each level on the optimization objective. We address this challenge by leveraging domain knowledge about underlying physical descriptors, which embody the physics of the underlying physical phenomena, to infer the effect of unobserved levels that have not been sampled yet. We show that physical descriptors can be intuitively embedded into the latent variable Gaussian process approach—a mixed-variable GP modeling technique—and used to selectively explore levels of qualitative variables in the Bayesian optimization framework. This physics-informed approach is particularly useful when one or more qualitative variables are high dimensional (many-level) and the modeling dataset is small, containing observations for only a subset of levels. Through a combination of mathematical test functions and materials design applications, our method is shown to be robust to certain types of incomplete domain knowledge and significantly reduces the design cost for problems with high-dimensional qualitative variables.

Engineering↗

DOE OSTI2023

Evaluation of the Near-Surface Variables in the HRRR Weather Model Using Observations from the ARM SGP Site

Abstract The performance of version 4 of the NOAA High-Resolution Rapid Refresh (HRRR) numerical weather prediction model for near-surface variables, including wind, humidity, temperature, surface latent and sensible fluxes, and longwave and shortwave radiative fluxes, is examined over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) region. The study evaluated the model’s bias and bias-corrected mean absolute error relative to the observations on different time scales. Forecasts of near-surface geophysical variables at five SGP sites (HRRR at 3-km scale) were found to agree well with observations, but some consistent observation–forecast differences also occurred. Sensible and latent heat fluxes are the most challenging variables to be reproduced. The diurnal cycle is the main temporal scale affecting observation–forecast differences of the near-surface variables, and almost all of the variables showed different biases throughout the diurnal cycle. Results show that the overestimation of downward shortwave and the underestimation of downward longwave radiative flux are the two major biases found in this study. The timing and magnitude of downward longwave flux, wind speed, and sensible and latent heat fluxes are also different with contributions from model representations, data assimilation limitations, and differences in scales between HRRR and SGP sites. The positive bias in downward shortwave and negative bias in longwave radiation suggests that the model is underestimating cloud fraction in the study domain. The study concludes by showing a brief comparison with version 3 of the HRRR and shows that version 4 has better performance in almost all near-surface variables. Significance Statement A correct representation of the near-surface variables is important for numerical weather prediction models. This study investigates the capability of the latest NOAA High-Resolution Rapid Refresh (HRRRv4) model in simulating the near-surface variables by comparing against the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) in situ observations. Among others, we find that the surface heat fluxes, such as sensible and latent heat fluxes, are the most difficult variables to be reproduced. This study also shows that the diurnal cycle has the dominant impact on the model’s performance, which means the majority of the outputted near-surface variables have the strong diurnal cycle in their bias errors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A Comparison of Properties of Quasars with and without Rapid Broad Absorption Line Variability

We investigate the correlation between rest-frame UV flux variability of broad absorption line (BAL) quasars and their variability in BAL equivalent widths (EWs) in a various timescale from <10 days to a few years in the quasar rest frame. We use the data sets of BAL EWs taken by the Sloan Digital Sky Survey Reverberation Mapping project and photometric data taken by the intermediate Palomar Transient Factory in the g and R bands and the Panoramic Survey Telescope and Rapid Response System in grizy bands. Our results are summarized as follows: (1) the distributions of flux variability versus BAL variability show weak, moderate, or a strong positive correlation; (2) there is no significant difference in flux variability amplitudes between BAL quasars with significant short timescale EW variability (called class S1) and without (class S2); (3) in all timescales considered in this paper, the class S1 quasars show systematically larger BAL variability amplitudes than those of the class S2 quasars; and (4) there are possible correlations between BAL variability and physical parameters of the quasars such as black hole masses (moderate positive), Eddington ratios, and accretion disk temperature (strong negative) in the class S2 quasars. These results indicate that the BAL variability requires changing in the ionizing continuum and an ancillary mechanism such as variability in X-ray shielding gas located at the innermost region of an accretion disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Variability in NGC 3201 Giant Stars and Its Impact on Their Spectroscopic [Fe/H] Determination

We present the analysis of 34 light curves in V and I of 17 giant stars in the globular cluster NGC 3201 to check if such stars are variable and if their variability has some kind of impact on the iron abundance as obtained from spectroscopic measurements. First, we computed the generalized Lomb–Scargle and phase dispersion minimization periodograms of the sample to check if the stars were variables. In this way, seven stars of the sample were found to be nonvariable, two stars are considered as possible variables, and eight stars were found to be variable, with periods ranging from 0.0881 ± 0.0001 to 0.5418 ± 0.0027 days. According to the literature, the variables have distinct values of [Fe I/H]: the three most metal-rich stars are in the red giant branch (RGB) stage, one has [Fe I/H] = −1.37 dex, while the other two have [Fe I/H] = −1.31 dex. The two most metal-poor variables have [Fe I/H] = −1.61 dex and [Fe I/H] = −1.62 dex, and are AGB stars; the remaining variables have [Fe I/H] = −1.44, −1.48, and −1.50 dex, the first two being RGB stars while the last is an AGB star. On the other hand, stars that appear to be nonvariable have −1.56 ≤ [Fe I/H] ≤ −1.40. We conclude that variability somehow affects the spectroscopic determination of the iron content of giant stars in NGC 3201, increasing the iron spread of the cluster. If variability is not taken into account, this spread could be incorrectly interpreted as due to an intrinsic iron spread affecting the stars of the cluster.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Assessing Photovoltaic Capacity Factor Variability Using Long-Term Satellite Derived Solar Resource Data Under Brazilian Climate

Accurate estimation of photovoltaic (PV) energy yield and its variability is essential for reducing financial risk and supporting reliable system planning for rapidly expanding PV markets. In Brazil, high solar adoption and increasing levels of distributed energy resources are beginning to introduce operational challenges such as curtailment and evolving grid requirements. Understanding how natural variability in solar resource propagates into PV system performance is therefore increasingly important for both project design and grid integration. Modern PV yield assessments commonly rely on multi-year meteorological datasets and probabilistic exceedance metrics (e.g., P50/P90) to quantify energy yield uncertainty for project financing. However, the implications of long-term solar resource variability for PV system design choices and high-adoption grid conditions remain less well characterized for rapidly expanding markets such as Brazil. In particular, understanding how weather-driven variability propagates into PV production distributions and capacity factor expectations is important for evaluating curtailment exposure, deployment strategies, and storage requirements in regions experiencing rapid growth of distributed and utility-scale PV. Seasonal and interannual variability in atmospheric conditions can produce substantial fluctuations in monthly PV energy production, which propagate into uncertainty in annual energy yield and capacity factor expectations. Characterizing this variability using long-term meteorological datasets allows probabilistic estimation of PV system performance and provides improved insight into the range of expected PV energy outcomes. This study explores the use of long-term satellite-derived meteorological data from the National Solar Radiation Database (NSRDB) to evaluate the variability of photovoltaic system performance across multiple locations in Brazil. Using a 27-year dataset (1998-2024), PV system simulations are performed to characterize the distribution of annual and seasonal capacity factors and energy yield outcomes, while propagating key sources of meteorological variability and model uncertainty through the PV modeling chain. The analysis also investigates the sensitivity of PV performance outcomes to key system design assumptions within the PV modeling chain, including tracking configuration and system sizing parameters. The resulting probabilistic performance characterization provides insight into how weather-driven variability influences PV production expectations and capacity factor distributions. These results provide a foundation for evaluating how weather-driven variability interacts with high PV adoption and potential storage or curtailment mitigation strategies.

14 SOLAR ENERGY↗

DOE OSTI2020

Quantifying Solar PV Variability at Multiple Timescales for Power Systems Studies: Preprint

Increasing penetration of solar photovoltaics (PV) has concerned power system operators about their ability to balance load and generation over multiple timescales and maintain acceptable frequency profile because of the variability and uncertainty associated with solar PV. This paper focuses on developing metrics to quantify solar PV variability at multiple temporal scales ranging from seconds to an hour to help power system operators and planners better understand and account for solar PV variability under high solar PV penetration. Fourier coefficients of solar irradiation spectrum and average size of solar irradiation ramps over a day are the two metrics that are used to quantity the variability. It is shown that by using the two metrics the days with highest and lowest solar irradiation variability can be identified, which is useful information for operations planning or long-term power system planning studies. The key findings of the paper are (i) single sensor-based variability quantification can significantly exaggerate the variability of a solar PV plant, particularly at seconds to 1-min temporal resolution, (ii) at 1-minute or lower temporal resolutions, single sensor-based variability quantification of an entire plant can be fairly accurate, with the accuracy increasing as the temporal resolution reduces from 1-minute to 1-hour, (iii) at each 1-minute and lower temporal resolutions, maximum solar irradiation ramps appear to be independent of the geographical location of the sensor, and (iv) significant reduction in solar PV plant variability over a region with multiple plants is possible at all temporal scales, but more high quality, high temporal resolution data is needed to quantify the reduction with high confidence.

photovoltaics↗

DOE OSTI2021

An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis

Uncertainty in predicting solar energy resources introduces major challenges in power system management and necessitates the development of reliable probabilistic solar forecasts. As the first part of the development of probabilistic forecasts based on the Weather Research and Forecasting model with solar extensions (WRF-Solar), this study presents a tangent linear approach to identify input variables responsible for the largest uncertainties in predicting surface solar irradiance and clouds. A tangent linear analysis is capable of efficiently investigating sensitivities of output variables with respect to various input variables of WRF-Solar because this approach avoids the computational burden of perturbing the initial conditions of individual input variables. We develop tangent linear models (TLMs) for six WRF-Solar physics packages that control the formation and dissipation of clouds and solar radiation, and we evaluate the validity of TLMs using a linearity test. The tangent linear sensitivity analysis is conducted under various scenarios based on satellite observations and model simulations to consider realistic input conditions. A simple method is used to quantify the impact of the uncertainty of input variables on the output variables from the TLMs. The results demonstrate that uncertainties in the output variables that are the focus of this study—including global horizontal irradiance, direct normal irradiance, cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes—are highly sensitive to uncertainties in 14 input variables. This study indicates that the tangent linear method can identify key variables of physics modules in WRF-Solar that can be stochastically perturbed to generate ensemble-based probabilistic forecasts.

14 SOLAR ENERGY↗

DOE OSTI2025

Key Environmental and Ecological Variables of Wetland CH 4 and CO 2 Fluxes Change With Warming

Wetlands are important ecosystems for the global carbon cycle, impacting regional and global methane (CH 4 ) and carbon dioxide (CO 2 ) budgets. This study examines how environmental and ecological variables impact wetland CH 4 flux and net ecosystem exchange of CO 2 (NEE) across 17 sites globally. We also quantified the importance of variables for each wetland type and site at monthly scale under normal and warm temperatures using dominance analysis. We identified soil and air temperature (TS, TA, respectively) as key variables influencing wetland CH4, and latent heat (LE) and shortwave radiation (SW) for NEE under normal and warm conditions. However, the importance of some variables shifted with warming. For predicting the variability of wetland CH4 flux under warming, gross primary productivity (GPP) and LE, replacing wind direction (WD), were dominant variables for tropical swamps, while NEE was important for high-latitude fens and bogs under warm temperatures. For wetland NEE, the role of TA and TS decreased across all wetland types with warming, while vapor pressure deficit (VPD) became more important for mid and high-latitude wetlands. Our results reveal the complex responses of wetland carbon flux to environmental and ecological variables with warming and provide new insights into improving wetland models by incorporating additional variables and accounting for the changing roles of variables in carbon flux under warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Time variability of ultra fast BAL outflows using SALT: C IV equivalent width analysis

We study the time variability (over ≤7.3 yr) of ultra fast outflows (UFOs) detected in a sample of 64 C IV broad absorption line (BAL) quasars (with 80 distinct BAL components) monitored using the Southern African Large Telescope. By comparing the properties of the quasar in our sample with those of a control sample of non-BAL quasars, we show that the distributions of black hole mass are different and the bolometric luminosities and optical photometric variations of UFO BAL quasars are slightly smaller compared to that of non-BAL quasars. The detection fraction of C IV equivalent width (W) variability (∼95 per cent), the fractional variability amplitude $\left(\frac{\Delta W}{W}\right)$ and the fraction of ‘highly variable’ BAL (i.e. $\big|$$\frac{\Delta W}{W}$$\big|$ > 0.67) components (∼33 per cent) are higher in our sample compared to the general BAL population. The scatter in $\frac{\Delta W}{W}$ and the fraction of ‘highly variable’ BALs increase with the time-scale probed. The $\frac{\Delta W}{W}$ distribution is asymmetric at large time scales. We attribute this to the BAL strengthening time-scales being shorter than the weakening time-scales. The BAL variability amplitude correlates strongly with the BAL properties compared to the quasar properties. BALs with low W, high-velocity, shallow profiles, and low-velocity width tend to show more variability. When multiple BAL components are present, a correlated variability is seen between low- and high-velocity components with the latter showing a larger amplitude variations. We find an anticorrelation between the fractional variations in the continuum flux and W. While this suggests photoionization induced variability, the scatter in continuum flux is much smaller than that of W.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Benchmarking Variables for Checkpointing in HPC Applications

Checkpoint/Restart (C/R) is a widely used fault tolerance mechanism in converged systems of cloud, edge, and HPC. However, users often rely on their experience to determine which variables to checkpoint, as there is currently no benchmark that can provide a reference. This can result in checkpointing redundant or even incorrect variables. To address this issue, we propose a benchmark suite that includes critical variables for checkpointing, which have been manually identified, and a method for identifying those critical variables, with 20 representative HPC applications. Our method involves analyzing data dependency between variables to identify critical variables analytically. We verify the identified variables' correctness with a widely used C/R library FTI by an ablation study. With our benchmark suite and data dependency analysis, HPC practitioners now have a reference for identifying checkpointing variables and better knowledge of what kind of variables to checkpoint.

Fu, Xiang↗

DOE OSTI2024

Impacts of stochastic coalescence variability on warm rain initiation using Lagrangian microphysics in box and large-eddy simulations

Various coalescence methods for Lagrangian microphysics schemes are tested in box and large-eddy simulation (LES) models, including the stochastic all-or-nothing super-droplet method (SDM) and a version of SDM (dSDM) that applies a fractional approach similar to the average impact method. In LES, variabilities driven by microphysics and by flow realizations are separated using the “piggybacking” technique. Rain initiation averaged over many realizations of the box model is delayed and rain variability increases as the number of super-drops per collision volume ( N SD ) is decreased using SDM. In contrast, rain initiation time using SDM in LES is insensitive to N SD for 32 ≤ N SD ≤ 512. This is explained through the interaction between LES grid boxes, each acting as a separate collision volume. Variability across the ensemble of LES collision volumes using SDM results in rain quickly initiating in some of the LES grid cells at low N SD and leading to a similar overall timing of rain initiation from the cloud compared to simulations with high N SD . There is a ~20% decrease in the total rain mass and mean rain flux as N SD is increased from 32 to 256, with little additional change as N SD is increased from 256 to 512. The fractional coalescence approach in dSDM leads to reduced microphysical variability and a 15-18 min delay in rain initiation compared to SDM. Finally, an additional LES ensemble with microphysical variability feeding back to the dynamics shows that flow variability dominates the impact of microphysical variability on rain properties. Thus, flow variability must be constrained to isolate impacts of microphysical variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Atmospheric Variability Drives Anomalies in the Bering Sea Air–Sea Heat Exchange

Abstract High latitudes, including the Bering Sea, are experiencing unprecedented rates of change. Long-term Bering Sea warming trends have been identified, and marine heatwaves (MHWs), event-scale elevated sea surface temperature (SST) extremes, have also increased in frequency and longevity in recent years. Recent work has shown that variability in air–sea coupling plays a dominant role in driving Bering Sea upper-ocean thermal variability and that surface forcing has driven an increase in the occurrence of positive ocean temperature anomalies since 2010. In this work, we characterize the drivers of the anomalous surface air–sea heat fluxes in the Bering Sea over the period 2010–22 using ERA5 fields. We show that the surface turbulent heat flux dominates the net surface heat flux variability from September to April and is primarily a result of near-surface air temperature and specific humidity anomalies. The airmass anomalies that account for the majority of the turbulent heat flux variability are a function of wind direction, with southerly (northerly) wind advecting anomalously warm (cool), moist (dry) air over the Bering Sea, resulting in positive (negative) surface turbulent flux anomalies. During the remaining months of the year, anomalies in the surface radiative fluxes account for the majority of the net surface heat flux variability and are a result of anomalous cloud coverage, anomalous lower-tropospheric virtual temperature, and sea ice coverage variability. Our results indicate that atmospheric variability drives much of the Bering Sea upper-ocean temperature variability through the mediation of the surface heat fluxes during the analysis period. Significance Statement A long-term ocean warming trend and a recent increase in marine heatwaves in the Bering Sea have been identified. Previous work showed that anomalies in the exchange of heat between the ocean and the atmosphere were the primary driver of Bering Sea temperature variability, but the processes responsible for the heat exchange anomalies were unknown. In this work, we show that the atmosphere is the primary driver of anomalies in the Bering Sea air–sea heat exchange and therefore plays an important role in altering the thermal state of the Bering Sea. Our results highlight the importance of understanding more about how the ocean and the atmosphere interact at high latitudes and how this relationship will be affected by future climate change.

Hayden, Emily E. [a College of Earth, Ocean, and A↗