Search NASA⌕ Search

SEARCH · Search NASA

Results for “robust model predictive control”

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 145 records · Page 8

Assessing the Impact of Cybersecurity Attacks on Energy Systems

This paper investigates the cyber resiliency of future power systems with high penetration of distributed energy resources using advanced distributed and (or) hierarchical control architectures. Specifically, we simulate cyberattacks on three prototypical use cases, and we identify attack scenarios that are the most damaging to the overall system performance. We show that these attacks can have a significant impact on grid operation. Results provide additional insights into the robustness of the system to the most common cyberattacks.

buildings↗

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Remote Americium Detection Using an Optical Sensor: A D-Optimal Strategy for Efficient PLS-Based Modeling

A fiber-optic visible–near-infrared absorption spectroscopy system in a glove box was demonstrated for remote quantification of Am(III) (0–500 µM) and HNO 3 (0.1–9 M) using partial least squares regression (PLSR) models. The sensor platform, featuring a simple plug-and-play spectrophotometer, can enable noninvasive, real-time monitoring of actinide process solutions. To establish a flexible PLSR model calibration strategy, a D-optimal design developed using Nd(III) in previous studies was successfully extended to an actinide system with Am(III) to effectively minimize sample set size while maintaining robust prediction performance. The results suggest strong spectral similarities between Nd(III) and Am(III) and validate Nd(III) as an effective optical surrogate for trivalent actinide species. This work also supports the generalizability of a D-optimal training set selection approach for two-factor systems. The PLS1 models for Am(III) and HNO 3 outperformed a PLS2 model and maintained reasonable performance in the presence of interfering U(VI). The resulting sensor system and multivariate approach provides a flexible and scalable solution for process monitoring, control, and safety in diverse nuclear applications.

actinide↗

Real-Time Implementation of Smart Wireless Charging of On-Demand Shuttle Service for Demand Charge Mitigation

This paper presents a smart charge management strategy for an on-demand electric shuttle operating at the National Renewable Energy Laboratory (NREL) campus and supported by an inductive charger at the vehicle's waiting spot. A new control algorithm has been proposed for mitigating the demand charges incurred from the wireless charger. It monitors the shuttle, wireless charger, renewable energy generation, and other loads and regulates charging behavior for demand charge mitigation. Within the control algorithm, an energy prediction is made to estimate the mobility needs of the vehicle and maintain uninterrupted service during operation while still minimizing peak demand. The proposed controller is designed and optimized using a Simulink model for the entire system. It is then implemented and tested in real time at the NREL campus using online cloud services. Two vehicle-use cases' charge-sustaining and charge-depletion operation' are tested under different campus power profiles and drive cycles to assess the controller's performance. In this work, the proposed controller showed a robust performance under different driving scenarios, with high correlation between simulation and experimental data. The results show that proper demand response can be achieved, with an average of 94% reduction of charging loads during peak demand events.

33 ADVANCED PROPULSION SYSTEMS↗

Charge Management for an Inductively Charged On-Demand Battery-Electric Shuttle Service with High Penetration of Renewable Energy

This paper presents a charge management control strategy for an on-demand battery-electric shuttle van operating at the National Renewable Energy Laboratory (NREL) campus and supported by day-time inductive charging at the vehicle's waiting spot. A new control algorithm has been proposed for reducing the demand charge costs incurred from wireless charging of the on-demand shuttle. A custom controller has been developed to monitor the shuttle, wireless charger, renewable energy generation, and various loads at NREL's campus, and regulate charging behavior for demand response. The intermittent renewable generation and sporadic operation of the on-demand shuttle service contribute to a high level of uncertainty in expected campus load profile, which must be carefully managed. The control algorithm predicts energy profile to estimate the mobility needs of the vehicle and maintain uninterrupted service during operation while still minimizing peak demand. The proposed controller has been designed and optimized using a Simulink model for the entire system. Next, it has been implemented and tested in real-time on the NREL campus. Two primary vehicle-use cases, charge sustaining and charge depletion operation, are tested under different load profiles and drive cycles to assess the controller's effectiveness at reducing peak demand and therefore demand charges. The proposed controller showed robust performance under different driving scenarios with high correlation between simulated and experimental data. The results showed that proper demand response can be achieved with an average of 94% reduction of charging loads during peak demand events.

ADVANCED PROPULSION SYSTEMS,ENERGY STORAGE↗

Establishment of a Vertically Integrated Domestic Manufacturing Process for Production of Substrates Needed for Manufacture of Gas Diffusion Layers

In this project, AvCarb, LLC evaluated the baseline performance metrics of commercial carbon veils and their corresponding Gas Diffusion Layers (GDLs) with the goal of establishing an optimized, vertically integrated production system for wet-laid nonwoven substrates used in gas diffusion media for electrochemical energy storage and conversion devices. Mechanical testing and microstructural characterization were conducted and used to develop a multiscale computational model capable of simulating and predicting the performance of GDLs in fuel cells. Although the project successfully generated foundational transport and modeling data, it was terminated prior to identifying the critical GDL design parameters necessary for full optimization. The program aimed to improve carbon veil fabrication through enhanced fiber dispersion, fiber-fiber adhesion control, and improved web formation, enabling the production of high-quality, uniform substrates. Simulations were intended to guide mixing and solution delivery system design and process conditions, followed by production-scale trials to evaluate fiber dispersion, web uniformity, and mechanical robustness. At full deployment, the proposed production line would have been capable of producing approximately 650,000 m² of carbon veil annually. This capability remains strategically important, as the United States currently lacks a domestic source of wet-laid nonwoven carbon substrates that satisfy the stringent quality requirements for fuel-cell GDLs and electrolyzers representing an ongoing supply-chain vulnerability. Beyond supply-chain benefits, the project established a robust benchmarking dataset for existing commercial carbon veils while advancing next-generation material concepts targeting improved performance and manufacturing consistency.

Olson, Cynthia Lemay↗

Intrinsic chemical reactivity of solid-electrolyte interphase components in silicon–lithium alloy anode batteries probed by FTIR spectroscopy

Herein we report the solid reaction products from the chemical reaction of aprotic battery electrolyte and three purported components of the Si-based anode SEI : SiO 2 nanoparticles (NPs), lithium silicate (Li x SiO y ) powders, and Si NPs. We use FTIR and classical molecular dynamics/density functional perturbation theory to assess the solid products remaining with these model materials after exposure to electrolyte. The absence of electrochemical bias provides a view of the chemical speciation resulting from early-stage chemical reactivity during battery assembly as well as under open circuit storage conditions. We believe these species represent the initial stages of SEI growth and predict they likely drive subsequent chemical and electrochemical reactions by controlling molecular interactons at the Si active material interface. We find that nominally equivalent materials react differently even before any electrochemistry is performed (e.g., acidic SiO 2 dissolves whereas alkaline SiO 2 is relatively robust), and derive new understanding of the chemical species that could and could not form stable SEI components in Si-based anodes. These results can be used to inform how to passivate Si anode surfaces and potentially generate an artificially engineered SEI that would be stable and enable next-generation battery anodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Field testing and validation of a low-cost MPC for demand flexibility for grid-interactive K-12 schools

K-12 school buildings account for the highest energy consumption within the public sector. Implementing advanced HVAC controls in grid-interactive K-12 schools could bring substantial economic advantages and grid flexibility. Our previous study demonstrated that a low-cost model predictive control (MPC) solution, which coordinates multiple packaged units, can enable demand flexibility without major hardware upgrades. However, a significant gap remains between academic pilots and market-ready scalable solutions. This paper extends the previous single-site pilot to a multi-site demonstration involving three school campuses (95 total units) through a commercial technology transfer process. Addressing the challenge of verifying performance with sparse field data, we present a new statistical approach using Bayesian methods to estimate the MPC’s effect on peak demand. Unlike traditional methods, this approach robustly quantifies uncertainty in non-normal, limited datasets. The results confirm the solution’s replicability, achieving a 21.6–38.9% reduction in HVAC peak demand (10.8–22.1% at the site-level) with > 98% probability across diverse locations. Finally, we document critical barriers to scaling software-as-a-service (SaaS) solutions–such as API instability and diverse legacy systems–and offer practical strategies to accelerate the commercial adoption of grid-interactive efficient buildings.

Ham, Sang Woo↗

Neural message-passing for objective-based uncertainty quantification and optimal experimental design

Various real-world scientific applications involve the mathematical modeling of complex uncertain systems with numerous unknown parameters. Accurate parameter estimation is often practically infeasible in such systems, as the available training data may be insufficient and the cost of acquiring additional data may be high. In such cases, based on a Bayesian paradigm, we can design robust operators retaining the best overall performance across all possible models and design optimal experiments that can effectively reduce uncertainty to enhance the performance of such operators maximally. While objective-based uncertainty quantification (objective-UQ) based on MOCU (mean objective cost of uncertainty) provides an effective means for quantifying uncertainty in complex systems, the high computational cost of estimating MOCU has been a challenge in applying it to real-world scientific/engineering problems. In this work, we propose a novel scheme to reduce the computational cost for objective-UQ via MOCU based on a data-driven approach. We adopt a neural message-passing model for surrogate modeling, incorporating a novel axiomatic constraint loss that penalizes an increase in the estimated system uncertainty. As an illustrative example, we consider the optimal experimental design (OED) problem for uncertain Kuramoto models, where the goal is to predict the experiments that can most effectively enhance robust synchronization performance through uncertainty reduction. We show that our proposed approach can accelerate MOCU-based OED by four to five orders of magnitude, without any visible performance loss compared to the state-of-the-art. The proposed approach applies to general OED tasks, beyond the Kuramoto model.

97 MATHEMATICS AND COMPUTING↗

From Data to Insights: A Covariate Analysis of the IARPA BRIAR Dataset for Multimodal Biometric Recognition Algorithms at Altitude and Range

This paper examines covariate effects on fused whole body biometrics performance in the IARPA BRIAR dataset, specifically focusing on UAV platforms, elevated positions, and distances up to 1000 meters. The dataset includes outdoor videos compared with indoor images and controlled gait recordings. Normalized raw fusion scores relate directly to predicted false accept rates (FAR), offering an intuitive means for interpreting model results. A linear model is developed to predict biometric algorithm scores, analyzing their performance to identify the most influential covariates on accuracy at altitude and range. Weather factors like temperature, wind speed, solar loading, and turbulence are also investigated in this analysis. The study found that resolution and camera distance best predicted accuracy and findings can guide future research and development efforts in long-range/elevated/UAV biometrics and support the creation of more reliable and robust systems for national security and other critical domains.

Bolme, David↗

An approach for fast and accurate simulation of phase change material based thermal energy storage in buildings

Latent heat thermal energy storage (LHTES) has significant potential for mitigating peak electricity demand and enabling load shifting in buildings. Phase Change Material embedded heat exchangers (PCM-HX) can significantly improve energy demand management due to high storage capacity. However, PCM-HX evaluation typically depends on computationally expensive fully transient simulations, posing significant challenges for scalable system- and building-level energy assessments across different climates and system architectures. This paper presents a generalized, accurate, and computationally efficient methodology for simulating building energy systems integrated with LHTES. The PCM-HX transient performance is represented by performance maps generated using a Generalized Resistance-Capacitance Model (GRCM) that enables accurate predictions of arbitrary PCM-HXs at low computational cost. The feasibility of the proposed approach was verified using a case study considering a dual-mode heat pump-thermal energy storage (HP-TES) system simulated in Modelica with Spawn of EnergyPlus™ for a DOE prototype small office building in two locations: Tampa, FL, and International Falls, MN. The PCM-HX performance maps provided accurate predictions of PCM-HX transient behavior, with mean absolute percentage deviations within 2–4% compared to GRCM while also achieving at least 1800× reduction in computational time. Moreover, the HP-TES system achieved energy savings of up to 17.4% in Tampa, FL, and 62.2% in International Falls, MN, demonstrating the broader applicability of the proposed methodology across different climate zones. This work highlights the importance of robust PCM-HX models in enabling accurate and computationally efficient building-level simulations and enabling future research opportunities for investigating optimized HP-TES designs and advanced control strategies for grid-interactive buildings.

Modelica↗

Thermal Energy Transport in Oxide Nuclear Fuel

To efficiently capture the energy of the nuclear bond, advanced nuclear reactor concepts seek solid fuels that must withstand unprecedented temperature and radiation extremes. In these advanced fuels, thermal energy transport under irradiation is directly related to reactor performance as well as reactor safety. The science of thermal transport in nuclear fuel is a grand challenge as a result of both computational and experimental complexities. Here we provide a comprehensive review of thermal transport research on two actinide oxides: one currently in use in commercial nuclear reactors, uranium dioxide (UO 2 ), and one advanced fuel candidate material, thorium dioxide (ThO 2 ). In both materials, heat is carried by lattice waves or phonons. Crystalline defects caused by fission events effectively scatter phonons and lead to a degradation in fuel performance over time. Bolstered by new computational and experimental tools, researchers are now developing the foundational work necessary to accurately model and ultimately control thermal transport in advanced nuclear fuels. We begin by reviewing research aimed at understanding thermal transport in perfect single crystals. The absence of defects enables studies that focus on the fundamental aspects of phonon transport. Next, we review research that targets defect generation and evolution. Here the focus is on ion irradiation studies used as surrogates for damage caused by fission products. We end this review with a discussion of modeling and experimental efforts directed at predicting and validating mesoscale thermal transport in the presence of irradiation defects. While efforts in these research areas have been robust, challenging work remains in developing holistic tools to capture and predict thermal energy transport across widely varying environmental conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

QC-GN 2 oMS 2 : a Graph Neural Net for High Resolution Mass Spectra Prediction

Predicting the mass spectrum of a molecular ion is often accomplished via three generalized approaches: rules-based methods for bond breaking, deep learning, or quantum chemical (QC) modeling. Rules-based approaches are often limited by the conditions for different chemical subspaces and perform poorly under chemical regimes with few defined rules. QC modeling is theoretically robust but requires significant amounts of computational time to produce a spectrum for a given target. Among deep learning techniques, graph neural networks (GNNs) have performed better than previous work with fingerprint-based neural networks in mass spectra prediction. To explore this technique further, we investigate the effects of including quantum chemically derived information as edge features in the GNN to increase predictive accuracy. The models we investigated include categorical bond order, bond force constants derived from extended tight-binding (xTB) quantum chemistry, and acyclic bond dissociation energies. Throughout this work, we evaluated these models against a control GNN with no edge features in the input graphs. Bond dissociation enthalpies yielded the best improvement with a cosine similarity score of 0.462 relative to the baseline model (0.437). In this work we also apply dynamic graph attention which improves performance on benchmark problems and supports the inclusion of edge features. Between implementations, we investigate the nature of the molecular embedding for spectra prediction and discuss the recognition of fragment topographies in distinct chemistries for further development in tandem mass spectrometry prediction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predictive analytics of selections of russet potatoes

We explore the application of machine learning algorithms specifically to enhance the selection process of Russet potato (Solanum tuberosum L.) clones in breeding trials by predicting their suitability for advancement. This study addresses the challenge of efficiently identifying high-yield, disease-resistant, and climate-resilient potato varieties that meet processing industry standards. Leveraging manually collected data from trials in the state of Oregon, we investigate the potential of a wide variety of state-of-the-art binary classification models. The dataset includes 1086 clones, with data on 38 attributes recorded for each clone, focusing on yield, size, appearance, and frying characteristics, with several control varieties planted consistently across four Oregon regions from 2013 to 2021. We conduct a comprehensive analysis of the dataset that includes preprocessing, feature engineering, and imputation to address missing values. We focus on several key metrics such as accuracy, F1-score, and Matthews correlation coefficient (MCC) for model evaluation. The top-performing models, namely a feedforward neural network classifier (Neural Net), a histogram-based gradient boosting classifier (HGBC), and a support vector machine classifier (SVM), demonstrate consistent and significant results. To further validate our findings, we conducted a simulation study using the aims, data-generating mechanisms, estimands, methods, and performance measures (ADEMP) framework, simulating different data-generating scenarios to assess model robustness and performance through true positive, true negative, false positive, and false negative distributions, area under the receiver operating characteristic curve (AUC-ROC) and MCC. The simulation results highlight that non-linear models like SVM and HGBC consistently show higher AUC-ROC and MCC than logistic regression, thus outperforming the traditional linear model across various distributions, and emphasizing the importance of model selection and tuning in agricultural trials. Variable selection further enhances model performance and identifies influential features in predicting trial outcomes. The findings emphasize the potential of machine learning in streamlining the selection process for potato varieties, offering benefits such as increased efficiency, substantial cost savings, and judicious resource utilization. Our study contributes insights into precision agriculture and showcases the relevance of advanced technologies for informed decision-making in breeding programs.

60 APPLIED LIFE SCIENCES↗

Boundary Layer Controls on the Shallow-to-Deep Cumulus Transition (Final Technical Report)

This project advanced our understanding of the processes governing cumulus cloud formation and provides an improved observational basis for validating earth system models. To be specific, the project used laser- and radar-remote sensors to examine the physical properties of updrafts that rise from earth’s surface and initiate clouds deeper in the atmosphere. These updrafts comprise “thermals” and “plumes” and occur at small spatial and temporal scales (e.g., 10s of minutes, 100s of meters). These small scales preclude explicit representation in most earth system and climate models, and thus necessitate “sub-grid-scale” parameterization of updraft processes. The innovation of this project was to directly measure the size, shape, strength and water vapor content of these updrafts with Doppler and Raman lidars, respectively, and to link these updraft properties to cloud processes using vertically pointed weather radars. The resulting data sets comprise 100s of thousands of updrafts and thousands of clouds, which far exceeds previous efforts, and thereby provides a robust statistical and physical representation of these processes. From these large datasets the project produced a sequence of scientific analyses that: (1) Elucidate how variations in the turbulent structure of the convective boundary layer control shallow cumulus convection, (2) Quantify the upward transport of water vapor to cloud base via thermals and plumes, (3) Validate large-eddy simulations of updrafts and shallow convective clouds, (4) Demonstrate a size-to-strength relationship between updraft width and updraft speed, and (5) Demonstrate how updrafts interact with the stability at the top of the convective boundary layer to modulate the depth and vigor of convective clouds. These results have been disseminated via several published journal articles, academic theses, and conference presentations. Collectively these results contribute to the Atmospheric System Research (ASR) program’s goal to “improve understanding of the key cloud, aerosol, precipitation, and radiation processes that affect the Earth’s radiative balance and hydrological cycle, particularly processes that limit the predictive ability of regional and global earth system models”.

54 ENVIRONMENTAL SCIENCES↗

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass

Abstract Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

59 BASIC BIOLOGICAL SCIENCES↗

CO2 response screen in grass Brachypodium reveals the key role of a MAP kinase in CO2-triggered stomatal closure

Abstract Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized B. distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a “chill1” mutant exhibited cooler leaf temperatures than wild-type Bd21-3 parent control plants after exposure to increased CO2. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3− sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

Lopez, Bryn N. K. (ORCID:0009000937288216)↗

Sea Ice and Cloud Processes Mediating Compensation between Atmospheric and Oceanic Meridional Heat Transports across the CMIP6 Preindustrial Control Experiment

Abstract Bjerknes compensation (BJC) refers to the anticorrelation observed between atmospheric and oceanic heat transport (AHT/OHT) variability, particularly on decadal to longer time scales that may be important to the predictability of the climate system. This study investigates the spread in BJC across fully coupled simulations of phase 6 of the Coupled Model Intercomparison Project (CMIP6) and critical processes (particularly related to sea ice and clouds) that may contribute to that spread. BJC on decadal to longer time scales is confirmed across all the simulations evaluated, and it is strongest in the Northern Hemisphere (NH) between 60° and 70°N. At these latitudes, BJC appears to be primarily driven by the exchange of turbulent fluxes (sensible and latent) in the Greenland, Iceland, and Barents Seas. Metrics to break down how sea ice and clouds uniquely modify the radiative balance of the polar atmosphere during anomalous OHT events are presented. These metrics quantify the impacts of sea ice and clouds on surface and top of atmosphere (latent, sensible, longwave, and shortwave radiative) energy fluxes. Cloud responses tend to counter the clear sky impacts over the Marginal Ice Zone (MIZ). It is further shown that the degree of BJC present in a simulation at high latitudes is heavily influenced by the sensitivity of the sea ice to OHT, which is most influential over the MIZ. These results are qualitatively robust across models and explain the intermodel spread in NH BJC in the preindustrial control experiment.

54 ENVIRONMENTAL SCIENCES↗