Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability analysis”

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 685 records · Page 38

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗

Nuclear Thermal Energy Storage Configurations for Industrial Combined Heat and Power Supply: Conceptual Study and Engineering Designs

The industries examined in this report primarily rely on moderate-temperature heat provided by gas- or coal-fired boilers and combined heat and power (CHP) plants, delivered through standard process steam systems. High-temperature energy demands are often industry-specific and typically exceed the capabilities of high-temperature gas-cooled reactors (HTGRs). While it is technically feasible to replace process steam from fossil-based heat sources with nuclear energy, certain industries, such as methanol production and pulp and paper, face technoeconomic challenges in integrating nuclear energy without major changes or a technological shift. This is mainly due to the limited external energy demand remaining after the use of internal byproducts, waste heat recovery, and simple efficiency improvements. Achieving full decarbonization of these processes with nuclear energy would require significant technological advancements, involving experimental technology and substantial investments, making widespread adoption in existing industrial plants unlikely in the near term. This study reviews TES options in the context of enabling a flexible CHP supply while maintaining a steady nuclear heat input. Heat storage systems that interface between the reactor primary fluid and the CHP system offer superior performance and flexibility. Specifically, steam extraction downstream of the reheater with a two-tank molten-salt TES appears as the best solution regarding thermodynamic system benefits and system drawbacks. Using selected system configurations, a conceptual design of an industrial energy park was developed for industries with varying energy demands, such as steel production plants utilizing electric arc furnaces (EAFs) and chemical plants, as well as for those with constant energy demands, like petroleum refineries. This design highlights the capabilities of TES and explores its potential business cases. The study also conceptually develops the potential for integrating additional energy sources with nuclear systems through the implementation of TES. The potential of the HTGR-TES-CHP system was also evaluated considering key uncertainties such as industrial demand profiles, external grid access availability, and eligible tax credit levels, using the Holistic Energy Resource Optimization Network. Sensitivity of net present value to these uncertainties was analyzed to determine the optimal number of nuclear reactors (and CHP systems) and the suitable TES capacity. The results were interpreted from a decision-maker’s perspective, focusing on three key areas: deployment strategy (oversized units vs. undersized units with TES support), industrial process characteristics (thermal-intensive single profiles vs. electricity-intensive combined profiles), and operational goals (maximizing profits vs. minimizing natural gas (NG) consumption or external grid dependence). The optimization results indicate that the HTGR-TES-CHP system significantly reduces reliance on NG boilers for individual industrial processes by 9-60% (in NG capacity factor), with an average reduction of 38%, compared to standalone NG boiler operation case (Business As Usual [BAU]). For combined industrial processes, the reduction ranges from 37-77%, with an average of 60%. Additionally, the system greatly reduces dependence on external grids. In meeting industrial electrical demands, a 33-100% self-sufficient internal electricity supply is achieved for single industrial process, with an average of 74%, compared to the BAU scenario, where 100% of electricity is imported. For combined processes, 35-100% of internal electricity demands are met by the reactor, with an average of 73%. At last, the relative NG price levels at which the proposed HTGR-TES-CHP system can cost-effectively enter the market currently dominated by existing NG boilers were estimated. For a moderate HTGR CAPEX level ($\$$2500/kWth, $\$$6329/kWe), the analysis suggests that NG prices must be 2.5 to 7 times higher than HTGR variable operating and maintenance costs for single industrial process, and 5.5 to 9.5 times higher for a combined process scenario. Tax credit modeling shows that the Investment Tax Credit significantly reduces the price threshold needed to break even, making the system competitive with NG boilers in certain cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Unsupervised Process Anomaly Detection and Identification Using the Leave-One-Variable-Out Approach

Automated anomaly detection and identification can signal equipment issues and pinpoint causes in large-scale industrial systems. For systems with limited failure history, unsupervised machine learning methods can be utilized as they do not require past failures. This study introduces the leave-one-variable-out (LOVO) model, which masks one variable at a time to predict the others, learning underlying process correlations. Detection performance was assessed with synthetic and experimental data, while identification performance used only synthetic data due to its ability to generate labeled anomaly types. For detection using synthetic data, the LOVO model generally outperformed comparative models; while using experimental data, the comparative methods outperformed the LOVO model. However, the comparative methods required selecting a latent size, and these conclusions pertain to using the optimal size. In practice, it would not be feasible to always select the optimal value, and incorrect selections impacted performance. In contrast, the LOVO model does not require a latent space. For identification using synthetic data, the LOVO model was slightly outperformed in interpretability and repeatability but still demonstrated impressive results. These outcomes suggest that the LOVO model is an effective model and may be more easily implemented without the challenging tuning process of selecting a latent size.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simplification of the Grid Model and its Impact on the Analysis of Electrical Power Systems

Here, this paper analyzes the impact of the use of Kron reduction on the state variables of a three-phase electrical system, even when it does not meet the necessary conditions for its application. Reduction is applied to a power line model to eliminate the equation corresponding to the neutral conductor of the line. The ATP program is used to model and simulate the behavior of an electrical system considering different degrees of disequilibrium as a reference for the comparison of results. The results show that under certain conditions of disequilibrium the Kron reduction can lead to significant errors in the state variables of the system.

Electric Power Systems↗

Plant Bioengineering Atlas: A Knowledge Graph of Genes, DNA Constructs, and Plant Traits.

Plant bioengineering has generated tens of thousands of genotype-to-phenotype relationships, but this knowledge remains fragmented across narrative literature and difficult to use computationally. Inconsistent descriptions of DNA constructs, host species, and traits, including variable species names, omitted regulatory elements, and inconsistent gene symbols, impede data reuse, comparative analysis, and design-build-test-learn cycles. Here, we present the Plant Bioengineering Atlas, a literature-mined, ontology-grounded knowledge base assembled using an artificial intelligence (AI)-aided extraction pipeline. A large language model parsed open-access primary research articles to generate structured, provenance-anchored records of engineered genes, modification types, promoter-gene-terminator constructs, host species, target traits, and reported phenotypes, with every record traceable to its source. The current release contains 14,358 curated records encompassing 6,998 distinct genes across 436 plant species from 6,452 papers published between 2000 and 2026. Corpus analysis reveals that experiments are concentrated in a small group of model and crop species, disease and pathogen resistance is the most frequently engineered trait class, and constitutive regulatory parts (particularly the CaMV 35S promoter and NOS terminator) remain pervasive. Two in five records omit one or both flanking regulatory elements (i.e., promoter and terminator), while only 23.4% describe cassettes in which both elements resolve to named part classes, exposing a systematic reproducibility gap. We organize these data into a knowledge graph linking genes, constructs, species, and traits; provide access through an interactive web portal; and propose an AI-compatible documentation standard for AI-ready reporting. The Plant Bioengineering Atlas provides a foundation for data-driven hypothesis generation and AI-aided plant biodesign.

, Genes, DNA Constructs↗

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗

Equivalence Ratio Scans in a Rotating Detonation Engine

Published data for run times exceeding a few seconds in rotating detonation engines (RDEs) is rare, but long duration tests are necessary to understand both transient and persistent phenomena in RDE operation. Previous work has noted different detonation wave modes at different equivalence ratios. In the current study, multiple long-duration runs with constant air flow rate and variable equivalence ratio are conducted in a hydrogen-fueled, water-cooled RDE. Analysis of the capillary tube averaged pressure (CTAP) data reveals a strong correlation with equivalence ratio as well as a small time dependency due to transient thermal processes. Ion probe measurements are used to analyze detonation wave frequencies and modes, both of which are strongly dependent on the equivalence ratio. Mode transitions are observed to occur consistently as the equivalence ratio is traversed between 0.55 and 0.60. These transitions are also found to be a function of the direction of changing equivalence ratio (ascending or descending) and time.

Boyette, Wesley R.↗

Equivalence Ratio Scans in a Rotating Detonation Engine

Published data for run times exceeding a few seconds in rotating detonation engines (RDEs) is rare, but long duration tests are necessary to understand both transient and persistent phenomena in RDE operation. Previous work has noted different detonation wave modes at different equivalence ratios. In the current study, multiple long-duration runs with constant air flow rate and variable equivalence ratio are conducted in a hydrogen-fueled, water-cooled RDE. Analysis of the capillary tube averaged pressure (CTAP) data reveals a strong correlation with equivalence ratio as well as a small time dependency due to transient thermal processes. Ion probe measurements are used to analyze detonation wave frequencies and modes, both of which are strongly dependent on the equivalence ratio. Mode transitions are observed to occur consistently as the equivalence ratio is traversed between 0.55 and 0.60. These transitions are also found to be a function of the direction of changing equivalence ratio (ascending or descending) and time.

Boyette, Wesley R.↗

Energy Storage and Decarbonization Analysis for Energy Regulators: Technical Analysis for the Illinois Commerce Commission

Jurisdictions around the world are enacting and enforcing an increasing number of policies to fight climate change, leading to higher penetration of variable renewable energy (VRE) and energy storage systems (ESSs) in the power grid. One of the biggest challenges associated with this process is the evaluation of the appropriate amount of ESS required to mitigate the variability of the VREs and achieve decarbonization goals of a particular jurisdiction. This report presents methodologies developed and results obtained for determining the minimum amount of ESS required to adequately serve load in a system where fossil fueled generators are being replaced by VREs over the next two decades. This technical analysis is performed by Sandia National Laboratories for the DOE Office of Electricity Energy Storage Program in collaboration with the Illinois Commerce Commission (ICC). The Illinois MISO Zone 4 is used as a case study. Several boundary conditions are investigated in this analysis including capacity adequacy and energy adequacy to determine the quantity of ESS required for MISO Zone 4. Multiple scenarios are designed and evaluated to incorporate the impact of varying capacity values of VREs and on the resource adequacy of the system. Several retirement scenarios involving fossil-fueled assets are also considered. Based on the current plans of new additions and retirements of generating assets, the results of the technical analysis indicate that Illinois MISO Zone 4 will require a significant quantity of ESS to satisfy their electricity demand over the next two decades.

25 ENERGY STORAGE↗

Data and code for: "Soil respiration response to simulated precipitation change depends on ecosystem type and study duration"

Data, code, and figures from our meta-analysis of the effect of precipitation change on soil respiration, with an emphasis on ecosystem variability and the effect of study duration. We combined the results of 81 separate studies to determine the effect of altered rainfall on the amount of carbon that leaves soil in the form of carbon dioxide, a major component of global carbon cycling. In addition we looked at how long the changes lasted, and how different soil properties and the intensity of precipitation changes affected the study results. We found that more precipitation resulted in greater amounts of carbon dioxide leaving the soil, and less precipitation resulted in less. However, the changes weakened over time in ecosystems that typically receive plenty of rainfall (e.g., forests), in contrast to ecosystems that typically receive little rainfall (e.g., deserts) where changes strengthened over time. Changes in the amount of carbon dioxide that left the soil were also affected by the amount of carbon in the soil, which impacts how much water soil can hold.

54 ENVIRONMENTAL SCIENCES↗

Understanding and Modeling Pooled Rideshare Acceptance: Influential Factors, Preferred User Experiences, and Implications

This dissertation explores factors influencing pooled rideshare (PR) adoption to provide actionable insights for transportation network companies (TNCs) and policymakers. PR allows travelers to share rides with unknown passengers, offering benefits such as cost reduction and congestion relief. However, adoption remains limited due to safety concerns, privacy issues, and trust in rideshare platforms. A national U.S. survey with 5,385 respondents examined transportation preferences and barriers to PR adoption. Exploratory and confirmatory factor analyses identified five key factors influencing PR consideration—safety, service experience, privacy, traffic/environment, and time/cost. Second factor analyses examined ways to optimize PR experiences, revealing four factors—comfort/ease of use, convenience, vehicle technology/accessibility, and passenger safety. Privacy concerns, for instance, using regression analysis, were found to reduce the likelihood of PR adoption by 77%, and convenience had the potential to increase it by 156%. The Pooled Rideshare Acceptance Model (PRAM), based on the Technology Acceptance Model, assessed the impact of these factors using the Structural Equation Model (SEM). Privacy, safety, trust, and convenience had a large effect (Cohen's f2 > 0.35) on PR acceptance, while multigroup analyses (PRAMMA) explored 16 demographic variables such as gender, generation, and income, emphasizing the need for tailored strategies. Based on all the statistical analysis and workshops using descriptive statistics, 95 actionable recommendations were made from the riders' perspective. Findings highlight the importance of customized services, user experience improvements, and policy interventions to enhance PR adoption. This dissertation provides a roadmap for future research and policy development, ensuring evidence-based, practical strategies to improve PR services in the U.S. and beyond.

Gangadharaiah, Rakesh↗

HFTS-2 Completions Design and State-of-the-Art Diagnostics Results

Various completion schemes were tested as part of the Hydraulic Fracture Test Site-2 (HFTS-2) project. The wells were instrumented with permanent fiber optic (FO) cable, and unique near-wellbore and far-field region surveys were acquired to evaluate spatial and temporal hydraulic fracture characteristics. The state-of-the-art diagnostics included substantial monitoring of both the well stimulation and longer-term production performance. The objective was to evaluate and improve stimulation distribution effectiveness (SDE) for different completion designs based on the diagnostics results. A hydraulic fracturing (HF) profile (i.e., proppant/fluid allocation among clusters) was obtained using near-wellbore Distributed Acoustic Sensing (DAS) data acquired during stimulation. Statistical analysis provided a quantitative assessment of SDE for the different completion designs. A novel workflow was developed to assess time-lapse fracture property changes at the cluster and stage levels based on Distributed Strain Sensing (DSS) data acquired during the production phase. These time-lapse fracture property changes were compared with the HF profile. Based on the statistical analysis, Rate Per Cluster (RPC) during proppant slurry placement appears to be a primary completion design variable, rather than the typical Rate Per Perf (RPP). If RPC was not maintained at a certain level, the likelihood of uneven distribution increased dramatically. The stage configurations that created low RPC while maintaining high RPP did not show improved SDE. The strain change results based on DSS data during the production phase indicated near-wellbore fracture property changes (e.g., aperture changes during well shut-in and flow periods) and aligned well with near-wellbore DAS signals acquired during stimulation. More importantly, fracture property changes based on DSS data acquired during two tests over a 7-month period show high consistency and good alignment with the near-wellbore HF profile. This workflow quantitatively evaluates SDE for different completion designs using diagnostic results and could help optimize completion strategy for future development projects. Here, the analysis and interpretation of results demonstrate the technical feasibility of using permanent fiber optic cable for stimulation and production monitoring and indic

58 GEOSCIENCES↗

Application of Atmospheric Gases and Particulate Matter to the Assessment of Urban Heat Island

Background: Urban heat island (UHI), where built areas are warmer compared to non-urban regions, increases human related diseases and mortality. A key challenge in UHI analysis is the designation of sites as urban or suburban/rural; however, the growing complexity of green spaces in urban areas and the predominance of the transportation sector in nonurban areas creates a dilemma for distinct delineation. Objectives: This study aims to utilize the variability of atmospheric components such as particulate matter (PM), inorganic gases, and volatile organic compounds (VOCs) as direct tracers of the degree of urbanization for ground-based measurements to fully comprehend UHI in convoluted regions with indistinct delineation of urban and nonurban environments. Methods: Atmospheric gases and aerosols were used as direct tracers of urbanization for UHI analysis. Inorganic gases and particulate matter were monitored in two sites in a southeastern US city with varying degrees of urbanization. VOCs were analyzed using a proton transfer reaction time-of-flight mass spectrometer. Results: The more-urbanized site exhibited warmer night conditions and elevated total oxidant levels, leading to the formation of nanometer-sized particles. Machine learning analysis revealed similar atmospheric pollutant profiles for both sites, suggesting comparable sources and variability. Biogenic VOCs were enhanced at the less-urbanized site; however, levels of anthropogenic aromatic VOCs were comparable for both sites. A comprehensive mass spectra analysis revealed distinct molecular backbones per site that further affirmed the applicability of VOCs as indicators of urbanization. Conclusion: This study concludes that VOCs provide more direct and accurate information than typical inorganic gases and PM parameters for characterizing the degree of urbanization. Further exploration of VOCs can enhance our understanding of UHI dynamics and its interaction with vegetation in urban green spaces.

Air quality sensor↗

Uncertainty in El Niño-like warming and California precipitation changes linked by the Interdecadal Pacific Oscillation

Abstract Marked uncertainty in California (CA) precipitation projections challenges their use in adaptation planning in the region already experiencing severe water stress. Under global warming, a westerly jet extension in the North Pacific analogous to the El Niño-like teleconnection has been suggested as a key mechanism for CA winter precipitation changes. However, this teleconnection has not been reconciled with the well-known El Niño-like warming response or the controversial role of internal variability in the precipitation uncertainty. Here we find that internal variability contributes > 70% and > 50% of uncertainty in the CA precipitation changes and the El Niño-like warming, respectively, based on analysis of 318 climate simulations from several multi-model and large ensembles. The Interdecadal Pacific Oscillation plays a key role in each contribution and in connecting the two via the westerly jet extension. This unifying understanding of the role of internal variability in CA precipitation provides critical guidance for reducing and communicating uncertainty to inform adaptation planning.

54 ENVIRONMENTAL SCIENCES↗

In-situ imaging of flexure-induced fracture in tape-laminate composites using high-resolution X-ray computed tomography

Here, this study presents development of an in-situ flexure test for imaging progressive inter- and intralaminar fracture in tape-laminate composites using X-ray computed tomography (CT). The intent of this test is to provide detailed experimental observations of ply-level damage that can be used to validate existing, and develop new, progressive damage analysis (PDA) tools. The test consists of a vertically mounted specimen which is flexed using two eccentric compressive loads using an in-situ uniaxial load stage. The flexure specimen contains a starter notch mid-span which promotes initiation of composite failure within the X-ray field of view. Specimens with two different laminate stacking sequences were tested, imaged, analyzed. For a quasi-isotropic laminate with large angle changes between adjacent plies, there was near simultaneous growth of transverse cracks and delaminations below the midplane of the laminate. For a laminate with small angles between adjacent plies, there was extensive formation of transverse crack networks which penetrated the laminate thickness without delamination growth. In addition to imaging fracture, the X-ray CT data from both specimen types were used to quantify the variability in thickness and analyze the local orientation of individual plies. Overall, the proposed test and the image-data analysis methodology provided an important insight into the fracture processes in tape laminates and highlights the inherent ply-level geometrical variabilities that should be accounted for in PDA simulations.

42 ENGINEERING↗

Efficient and robust stability analysis in the internal energy, volume, and moles (UVN) space

Hydrodynamic simulation of fluids in vessels is often based on the internal energy (U), volume (V), and the number of moles of all components ( re N = N 1 , N 2 , ..., N n ) as the working variables to determine the thermodynamic state by finding the global maximum of entropy. The procedure of finding the equilibrium state may be divided into two steps: stability analysis and phase-split computations. Stability analysis is performed first, and the phase-split computation proceeds if the system is unstable. Most of the past work on UVN space has been based on phase-split computations, while the stability analysis is barely investigated. Furthermore, previous studies for the stability analysis use only Newton's method to solve the non-linear algebraic equations. We present an efficient and robust approach for stability analysis where UVN is specified. The successive substitution iteration (SSI) is used to provide good initial guesses for Newton's method. The proposed approach results in a reduced number of unknowns and does not require a large number of iterations to achieve convergence in Newton's method. Overall, our proposed formulation is compatible with different equations of state and is applicable to both pure component and multi-components. The robustness and efficiency of the algorithm for stability analysis are demonstrated in various examples.

02 PETROLEUM↗

Applications of the Dulmage–Mendelsohn decomposition for debugging nonlinear optimization problems

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

42 ENGINEERING↗

Quantifying the Role of Ocean Dynamics in SST Variability across GCMs and Observations

Abstract Midlatitude SSTs forced by mesoscale oceanic processes can affect the large-scale atmosphere, pointing to the ocean’s crucial role outside the tropics. Previous studies have shown oceanic mesoscale processes’ effect on global and regional climate variability. This study quantifies the local contribution of ocean dynamics to mixed-layer temperature across the globe by directly estimating the ocean heat flux divergence resolved by state-of-the-art ocean reanalysis, eddy-resolving, and eddy-parameterized versions of two U.S. national climate models and indirectly from air–sea flux satellite-based estimates. Our results show that the eddy-resolving climate simulations resolve mixed-layer temperature variances that are larger and closer to those inferred from observations than both their eddy-parameterized counterparts and ECCO over much of the extratropics. The observations and the eddy-resolving models indicate a more significant role of ocean dynamics in the mixed-layer temperature variability than the surface fluxes over most extratropics compared to their eddy-parameterized versions. A frequency domain analysis shows that the better-resolved ocean mesoscale and thermal gradients enhance the variance over a time scale from 2 months to 30 years. Results show agreement in the ocean’s contribution among satellite-based estimates, ocean reanalysis products, and ocean eddy-resolving simulations. At the same time, differences emerge for ECCO and the eddy-parameterized models, suggesting that surface fluxes account for a larger fraction of the mixed-layer temperature variability in most of the extratropics.

Siqueira, L.↗