Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variables”

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 127 records · Page 7

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↗

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↗

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↗

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↗

Climatology and natural variability of the global hydrologic cycle in the GLA atmospheric general circulation model

Time average climatology and low-frequency variabilities of the global hydrologic cycle (GHC) in the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) were investigated in the present work. A 730-day experiment was conducted with the GLA GCM forced by insolation, sea surface temperature, and ice-snow undergoing climatological annual cycles. Ifluences of interactive soil moisture on time average climatology and natural variability of the GHC were also investigated by conducting 365-day experiments with and without interactive soil moisture. Insolation, sea surface temperature, and ice-snow were fixed at their July levels in the latter two experiments. Results show that the model's time average hydrologic cycle variables for July in all three experiments agree reasonably well with observations. Except in the case of precipitable water, the zonal average climates of the annual cycle experiment and the two perpetual July experiments are alike, i.e., their differences are within limits of the natural variability of the model's climate. Statistics of various components of the GHC, i.e., water vapor, evaporation, and precipitation, are significantly affected by the presence of interactive soil moisture. A long-term trend is found in the principal empirical modes of variability of ground wetness, evaporation, and sensible heat. Dominant modes of variability of these quantities over land are physically consistent with one another and with land surface energy balance requirements. The dominant mode of precipitation variability is found to be closely related to organized convection over the tropical western Pacific Ocean. The precipitation variability has timescales in the range of 2 to 3 months and can be identified with the stationary component of the Madden-Julian Oscillation. The precipitation mode is not sensitive to the presence of interactive soil moisture but is closely linked to both the rotational and divergent components of atmospheric moisture transport. The present results indicate that globally coherent natural variability of the GHC in the GLA GCM has two basic timescales in the absence of annual cycles of external forcings: a long-term trend associated with atmosphere-soil moisture interaction which affects the model atmosphere mostly over midlatitude continental regions and a large-scale 2- to 3-month variability associated with atmospheric moist processes over the western Pacific Ocean.

Lau, K.-M.↗

Decadal variability of the Tropical Atlantic Ocean Surface Temperature in shipboard measurements and in a Global Ocean-Atmosphere model

Sea surface temperature (SST) variability was investigated in a 200-yr integration of a global model of the coupled oceanic and atmospheric general circulations developed at the Geophysical Fluid Dynamics Laboratory (GFDL). The second 100 yr of SST in the coupled model's tropical Atlantic region were analyzed with a variety of techniques. Analyses of SST time series, averaged over approximately the same subregions as the Global Ocean Surface Temperature Atlas (GOSTA) time series, showed that the GFDL SST anomalies also undergo pronounced quasi-oscillatory decadal and multidecadal variability but at somewhat shorter timescales than the GOSTA SST anomalies. Further analyses of the horizontal structures of the decadal timescale variability in the GFDL coupled model showed the existence of two types of variability in general agreement with results of the GOSTA SST time series analyses. One type, characterized by timescales between 8 and 11 yr, has high spatial coherence within each hemisphere but not between the two hemispheres of the tropical Atlantic. A second type, characterized by timescales between 12 and 20 yr, has high spatial coherence between the two hemispheres. The second type of variability is considerably weaker than the first. As in the GOSTA time series, the multidecadal variability in the GFDL SST time series has approximately opposite phases between the tropical North and South Atlantic Oceans. Empirical orthogonal function analyses of the tropical Atlantic SST anomalies revealed a north-south bipolar pattern as the dominant pattern of decadal variability. It is suggested that the bipolar pattern can be interpreted as decadal variability of the interhemispheric gradient of SST anomalies. The decadal and multidecadal timescale variability of the tropical Atlantic SST, both in the actual and in the GFDL model, stands out significantly above the background 'red noise' and is coherent within each of the time series, suggesting that specific sets of processes may be responsible for the choice of the decadal and multidecadal timescales. Finally, it must be emphasized that the GFDL coupled ocean-atmosphere model generates the decadal and multidecadal timescale variability without any externally applied force, solar or lunar, at those timescales.

Mehta, Vikram M.↗

An Expanded RXTE Survey of Long-Term X-ray Variability in Seyfert 1 Galaxies

The first seven years of RXTE monitoring of Seyfert 1 active galactic nuclei have been systematically analyzed to yield five homogenous samples of 2-12 keV light curves, probing hard X-ray variability on successively longer durations from approx. 1 day to approx. 3.5 years. 2-10 keV variability on time scales of approx. 1 day, as probed by ASCA, are included. All sources exhibit stronger X-ray variability towards longer time scales, with variability amplitudes saturating at the longest time scales, but the increase is greater for relatively higher luminosity sources. The well-documented anticorrelation between variability amplitude and luminosity is confirmed on all time scales. However, anticorrelations between variability amplitude and black hole mass estimate are evident on only the shortest time scales probed. The data are consistent with the models of power spectral density (PSD) movement described in Markowitz et al. (2003) and McHardy et al. (2004), whereby Seyfert 1 galaxies variability can be described by a single, universal PSD shape whose cutoff frequency scales with black hole mass. The best-fitting scaling relations between variability time scale, black hole mass and X-ray luminosity support an average accretion rate of 2% of the Eddington limit for the sample. Nearly all sources exhibit stronger variability in the relatively soft 2-4 keV band compared to the 7-12 keV band on all time scales. Color-flux diagrams support also Seyfert 1s' softening as they brighten. There are indications that relatively less luminous or less massive sources exhibit a greater degree of spectral variability for a given increase in overall flux.

Markowitz, A.↗

Variability Selected Low-Luminosity Active Galactic Nuclei in the 4 Ms Chandra Deep Field-South

The 4 Ms Chandra Deep Field-South (CDF-S) and other deep X-ray surveys have been highly effective at selecting active galactic nuclei (AGN). However, cosmologically distant low-luminosity AGN (LLAGN) have remained a challenge to identify due to significant contribution from the host galaxy. We identify long-term X ray variability (approx. month years, observed frame) in 20 of 92 CDF-S galaxies spanning redshifts approx equals 00.8 - 1.02 that do not meet other AGN selection criteria. We show that the observed variability cannot be explained by X-ray binary populations or ultraluminous X-ray sources, so the variability is most likely caused by accretion onto a supermassive black hole. The variable galaxies are not heavily obscured in general, with a stacked effective power-law photon index of Gamma(sub Stack) approx equals 1.93 +/- 0.13, and arc therefore likely LLAGN. The LLAGN tend to lie it factor of approx equal 6-89 below the extrapolated linear variability-luminosity relation measured for luminous AGN. This may he explained by their lower accretion rates. Variability-independent black-hole mass and accretion-rate estimates for variable galaxies show that they sample a significantly different black hole mass-accretion-rate space, with masses a factor of 2.4 lower and accretion rates a factor of 22.5 lower than variable luminous AGNs at the same redshift. We find that an empirical model based on a universal broken power-law power spectral density function, where the break frequency depends on SMBH mass and accretion rate, roughly reproduces the shape, but not the normalization, of the variability-luminosity trends measured for variable galaxies and more luminous AGNs.

Young, M.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Implementation of a feature selection algorithm in FARM to identify important state variables and time-invariant matrices

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM aids the HERON software module in the evaluation of the optimal dispatch for the different IES components. Set-point trajectories are required to meet limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To evaluate the feasibility of HERON generated set-points and to do so in an acceptable time, FARM employs reduced order models to represent the dynamic behavior of the systems to be dispatched. These surrogate models take the form of a linear dynamic system with sets of Linear Parameter Varying (LPV) matrices that are mapped to the system operating space. These matrices are derived from the trajectories of system state variables and system output variables during transients. The accuracy of LPV matrices depends on the selection of state variables. In previous reports, state variables were selected by adopting a complicated workflow requiring multiple software licenses and an advanced level of user expertise. In this report, a new workflow that automates the state variable selection process is presented. It significantly reduces the frequency of user interventions and does not require multiple software licenses. Each module in the new workflow is described in detail, and the input / output examples in each step of the workflow are provided. It was demonstrated that this workflow can greatly reduce the complexity of the state variable selection process, and that the updated FARM-Gamma and FARM-Delta validators can benefit from this workflow when solving the power dispatch problem of a representative IES test case. Finally, some code improvements that can further enhance the efficiency are suggested.

42 ENGINEERING↗

Cooperative Agreement To Analyze variabiLity, change and predictabilitY in the earth SysTem (CATALYST)

CATALYST proposes to perform foundational coordinated research in a team-oriented collaborative effort aimed at advancing a robust understanding of modes of Earth system variability and change using models, observations and process studies. The proposed research will address the DOE/BER mission by exploring the limits to predictability, identifying fundamental underlying mechanisms, quantifying interactions among modes of variability, and discovering tipping points in the Earth system to understand the current and future impacts of these phenomena on regional and global climate. Four fundamental gaps are identified in our knowledge of the Earth system: 1) What are the limits to predictability on various timescales? 2) What are the interactions among modes of Earth system variability? 3) How may modes of Earth system variability change in response to changes in external forcing, and what are the tipping points involved with those changes? 4) How are high impact events connected to modes of Earth system variability and how may they change in the future? Related to those gaps in our knowledge, we formulate four research objectives to address those gaps using a combination of Earth system models (ESMs) and machine learning (ML) methods. Research Objective 1 (RO1) addresses the first gap above and proposes to understand modes of variability and their limits of predictability on subseasonal to decadal timescales using ESMs and ML. Research Objective 2 (RO2) addresses the second gap and proposes to use a hierarchy of models to understand relevant processes and feedbacks related to how modes of variability interact with each other. Research Objective 3 (RO3) is designed to study the third gap and proposes to examine the role of external forcings in changes of modes of Earth system variability and their interactions, and the likelihood and predictability of tipping points and irreversible changes. Research Objective 4 (RO4) will address the fourth gap and proposes to use high resolution ESMs, regionally refined models (RRMs), and ML methods to investigate the relationships between high impact events (e.g. flash droughts and precipitation extremes, atmospheric rivers (ARs), tropical cyclones (TCs), storm surge/sea level rise), the synoptic systems that produce them, and their changes related to modes of Earth system variability. The research will involve the use of the Community Earth System Model (CESM), Energy Exascale Earth System Model (E3SM), CMIP multi-model data sets, a hierarchy of simpler models, and numerous observational data sets. In the course of the proposed research, CATALYST will contribute to metrics and diagnostics that will be integrated in Coordinated Model Evaluation Capabilities (CMEC), particularly with regards to the Quasi-biennial Oscillation (QBO) and its interactions with the Madden-Julian Oscillation (MJO), high atmospheric pressure blocking, and new precipitation metrics.

54 ENVIRONMENTAL SCIENCES↗

Baseline Climate Variables for Earth System Modelling

The Baseline Climate Variables for Earth System Modelling (ESM-BCVs) are defined as a list of 135 variables which have high utility for the evaluation and exploitation of climate simulations. The list reflects the most frequently used elements of the Coupled Model Intercomparison Project Phase 6 (CMIP6) archive. Successive phases of CMIP have supported strong results in science and substantially influence international climate policy formulation. This paper responds to both interest in exploiting CMIP data standards in a broader range of climate modelling activities and a need to achieve greater clarity about the significance and intention of variables in the CMIP Data Request. As Earth system modelling archives grow in scale and complexity, there are emerging problems associated with weak standardisation at the variable collection level. That is, there are good standards covering how specific variables should be archived, but this paper fills a gap in the standardisation of which variables should be archived. The ESM-BCV list is intended as a resource for ESM intercomparison projects (MIPs) developing requests to enable greater consistency among MIPs and as a reference for modelling centres to enhance consistency within MIPs. Provisional planning for the CMIP7 Data Request exploits the ESM-BCVs as a core element. The baseline variable list includes 98 variables which have modest or minor data volume footprints and could be generated systematically when simulations are produced and archived for exploitation by the World Climate Research Programme (WCRP) community. A further 35 variables are classed as “high volume” and are only suitable for production when the resource implications are justified.

Juckes, Martin [University of Oxford (United Kingd↗