SEARCH · Engineering Papers

Results for “Linear response”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Kerr, Faraday, and magnetoelectric effects in MnBi 2 Te 4 thin films

The topological magnetoelectric effect (TME) is a characteristic property of topological insulators. In this paper, we use a simplified coupled-Dirac-cone electronic structure model to theoretically evaluate the THz and far infrared Kerr and Faraday responses of thin films of MnBi 2 ⁢Te 4 with up to N = 10 septuple layers with the goal of clarifying the relationship between these convenient magnetooptical observables and the TME. Here we find that for even N, the linear Kerr and Faraday responses to an electric field vanish in the low-frequency limit, even though the magnetoelectric response is large and approximately quantized.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

The Role of Unit-Cell Topology in Modulating the Compaction Response of Additively Manufactured Cellular Materials using Simulations and Validation Experiments

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.

36 MATERIALS SCIENCE↗

DOE OSTI2024

The role of unit cell topology in modulating the compaction response of additively manufactured cellular materials using simulations and validation experiments

36 MATERIALS SCIENCE↗

DOE OSTI2023

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗

DOE OSTI2021

Shotgun Immunoproteomic Approach for the Discovery of Linear B-Cell Epitopes in Biothreat Agents Francisella tularensis and Burkholderia pseudomallei

Peptide-based subunit vaccines are coming to the forefront of current vaccine approaches, with safety and cost-effective production among their top advantages. Peptide vaccine formulations consist of multiple synthetic linear epitopes that together trigger desired immune responses that can result in robust immune memory. The advantages of linear compared to conformational epitopes are their simple structure, ease of synthesis, and ability to stimulate immune responses by means that do not require complex 3D conformation. Prediction of linear epitopes through use of computational tools is fast and cost-effective, but typically of low accuracy, necessitating extensive experimentation to verify results. On the other hand, identification of linear epitopes through experimental screening has been an inefficient process that requires thorough characterization of previously identified full-length protein antigens, or laborious techniques involving genetic manipulation of organisms. In this study, we apply a newly developed generalizable screening method that enables efficient identification of B-cell epitopes in the proteomes of pathogenic bacteria. As a test case, we used this method to identify epitopes in the proteome of Francisella tularensis (Ft), a Select Agent with a well-characterized immunoproteome. Our screen identified many peptides that map to known antigens, including verified and predicted outer membrane proteins and extracellular proteins, validating the utility of this approach. We then used the method to identify seroreactive peptides in the less characterized immunoproteome of Select Agent Burkholderia pseudomallei (Bp). This screen revealed known Bp antigens as well as proteins that have not been previously identified as antigens. Although B-cell epitope prediction tools Bepipred 2.0 and iBCE-EL classified many of our seroreactive peptides as epitopes, they did not score them significantly higher than the non-reactive tryptic peptides in our study, nor did they assign higher scores to seroreactive peptides from known Ft or Bp antigens, highlighting the need for experimental data instead of relying on computational epitope predictions alone. The present workflow is easily adaptable to detecting peptide targets relevant to the immune systems of other mammalian species, including humans (depending upon the availability of convalescent sera from patients), and could aid in accelerating the discovery of B-cell epitopes and development of vaccines to counter emerging biological threats.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Using system identification in modeling the yaw response of tail fins for small wind turbines with bearing friction

Here, we describe three main developments of our previous study of the nonlinear yaw dynamics of the tail fins for small wind turbines [Khedr et al., J. Renewable Sustainable Energy 16, 053305 (2024)]. First, the model constants derived from archived computational and experimental studies are adjusted by employing system identification (SI) to maximize the model's agreement with wind tunnel tests. This adjustment was done for high wind speeds, where yaw bearing friction can be ignored. When starting a turbine at low wind speed, however, friction can become important. Our second development is to implement a model for the frictional resistive torque and use SI to maximize its accuracy. These developments used wind tunnel experiments on generic delta, elliptical, and rectangular planforms that were described by Khedr et al. [J. Renewable Sustainable Energy 16, 053305 (2024)]. Since the aerodynamic and friction models employ a large number of constants, we describe ways to constrain the values using linearized solutions of the response equations for small and large yaw angles. Third, we test the generality of the aerodynamic and friction modeling using the complex planform from a commercial small turbine, for which limited theoretical and computational guidance is available in selecting the model constants. The model implemented with SI is shown to provide an accurate description of the yaw response of the complex planform. Guidelines are given for the use of wind tunnel tests to determine the model constants for tail fins of any planform.

17 WIND ENERGY↗

DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

Quasi‐dynamic breathing model of the lung incorporating viscoelasticity of the lung tissue

Abstract We advanced a novel model to calculate viscoelastic lung compliance and airflow resistance in presence of mucus, accounting for the quasi‐linear viscoelastic stress–strain response of the parenchyma (alveoli) tissue. We adapted a continuum‐based numerical modeling approach for the lung, integrating the fluid mechanics of the airflow within individual generations of the bronchi and alveoli. The model accounts for elasticity of the deformable bronchioles, resistance to airflow due to the presence of mucus within the bronchioles, and subsequent mucus flow. Simulated quasi‐dynamic inhalation and expiration cycles were used to characterize the net compliance and resistance of the lung, considering the rheology of the mucus and viscoelastic properties of the parenchyma tissue. The structure and material properties of the lung were identified to have an important contribution to the lung compliance and airflow resistance. The secondary objective of this work was to assess whether a higher frequency and smaller volume of harmonic air flow rate compared to a normal ventilator breathing cycle enhanced mucus outflow. Results predict, lower mucus viscosity and higher excitation frequency of breathing are favorable for the flow of mucus up the bronchi tree, towards the trachea.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

High-performance organic photodetectors enabled by a refined fibrillar multiphase morphology

Film morphology of donor:acceptor blend layers plays a critical role in photon-to-current efficiency and dark/noise current of organic photodetectors (OPDs). One effective approach to manipulate crystallization and mesoscale phase separation of such blend layers is sequential casting (SC). However, the guiding strategies to control the morphology and the impacts on OPD performance of the SC films remain elusive, as the film structural evolution during SC is different from conventional blend casting (BC). In this report a refined fibrillar multiphase morphology is demonstrated by SC in an NT812:IEICO-4F blend film, where the high-quality IEICO-4F crystallites infiltrate the robust NT812 fibril network from the surface to the bulk. Such a morphology leads to improved charge generation and collection, reduced trap states, and enhanced charge block capability of resultant device, enabling simultaneous achievement of high external quantum efficiency and low dark/noise current. A maximum special detectivity of 5 × 10 13 Jones is achieved at 860 nm under –0.1 V, which is among the highest detectivities for vis-to-NIR OPDs. The linear dynamic range and response speed are also improved. Such enhancements are parallelly observed from OPD devices based on other blend systems with similar fibrillar refinement, which provides guidelines for film structure manipulation towards high OPD performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Mapping the proteogenomic landscape enables prediction of drug response in acute myeloid leukemia

Acute myeloid leukemia is a poor prognosis cancer commonly stratified by genetic aberrations, but these mutations are often heterogeneous and don’t always predict therapeutic response. Here we combine transcriptomic, proteomic, and phosphoproteomic datasets with ex vivo drug sensitivity data to help understand the underlying pathophysiology of AML beyond mutations. We measured the proteome and phosphoproteome of 210 patients and combined them with genomics and transcriptomic measurements to identify four proteogenomic subtypes that complemented existing genetic subtypes. We then built a predictor to classify samples into subtypes based on 147 molecular features and mapped them to a ‘landscape’. Each region of this landscape corresponded to specific drug response patterns. We then built a drug response prediction model to identify drugs that target distinct subtypes. We can ultimately use these models to predict drug treatment response and prioritize treatments. Finally, we extended our models and mapped a series of cell lines representing various stages of quizartinib resistance into our subtype landscape, predicting and experimentally validating a switch in sensitivity to venetoclax to panobinostat, two drugs with very different mechanisms than quizartinib. Our results show how multi-omics data together with drug sensitivity data can inform therapy stratification and drug combinations in AML.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Theory of an Automatic Seepage Meter and Ramifications for Applications

A new approach for measuring fluxes across surface water—groundwater interfaces was recently proposed. The Automatic Seepage Meter (ASM) is equipped with a precise water level sensor and digital memory that analyzes water level time series in a vertical tube inserted into a streambed. The ability to infer flux values with high temporal resolution relies on an accurate interpretation of water level dynamics inside the tube. Here, we reduce the three-dimensional hydrodynamic problem that describes the ASM water level in a variety of field conditions to a single ordinary differential equation. This novel general analytical solution for estimating ASM responses is more comprehensive and flexible than previous approaches and is applicable to the entire range of field conditions, including steady or transient stream stages, evaporation, rainfall, and noise. For example, our analysis determines the timing of the nonmonotonic ASM response to a monotonic linear stream stage variation and explains previously used empirical parabolic approximation for estimating fluxes. We present algorithms for simultaneous inference of vertical interface flux and hydraulic conductivity values together with an example code. We quantify how the accuracy of parameter estimation depends on test duration and noise amplitude and propose how our analysis can be used to optimize field test protocols. On this basis, changing the ASM geometry by increasing the radius and decreasing tube insertion depth may enable ASM field test protocols that estimate interface flux and hydraulic conductivity faster while maintaining desired accuracy. Potential applications of joint parameter estimation are suggested.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Origin of Rashba Spin-Orbit Coupling in 2D and 3D Lead Iodide Perovskites

isaWhittaker-Brooks2 & Tho D. Nguyen1* We studied spin dynamics of charge carriers in the superlattice-like Ruddlesden-Popper hybrid lead iodide perovskite semiconductors, 2D (BA)2(MA)Pb2I7 (with MA= CH3NH3, and BA= CH3(CH2)3NH3), and 3D MAPbI3 using the magnetic field effect (MFE) on conductivity and electroluminescence in their light emitting diodes (LEDs) at cryogenic temperatures. The semiconductors with distinct structural/ bulk inversion symmetry breaking, when combined with colossal intrinsic spin–orbit coupling (SOC), theoretically give rise to giant Rashba-type SOC. We found that the magneto-conductance (MC) magnitude increases monotonically with the emission intensity and saturates at ≈0.05% and 0.11% for the MAPbI3 and (BA)2(MA)Pb2I7, respectively. The magneto-electroluminescence (MEL) response with similar line shapes as the MC response has a significantly larger magnitude, and essentially stays constant at ≈0.22% and ≈0.20% for MAPbI3 and (BA)2(MA)Pb2I7, respectively. The sign and magnitude of the MC and MEL responses can be quantitatively explained in the framework of the Δg-based excitonic model using rate equations. Remarkably, the width of the MEL response in those materials linearly increases with increasing the applied electric field, where the Rashba coefficient in (BA)2(MA) Pb2I7 is estimated to be about 7 times larger than that in MAPbI3. Our studies might have significant impact on future development of electrically-controlled spin logic devices via Rashba-like effects

36 MATERIALS SCIENCE↗

DOE OSTI2023

Sensitivity measurements for a 250 MHz quartz shear-horizontal surface acoustic wave biosensor under liquid viscous loading

Surface acoustic wave (SAW) devices have been used in biochemical assays due to their high sensitivity. The device sensitivity is a function of changes in the density and viscosity of the liquid. Here, we studied the effect of fluid viscosity using a 250 MHz quartz shear-horizontal (SH)-SAW biosensor by monitoring different concentrations of binary aqueous/glycerol solutions. In this study, the sensitivity of the biosensor was determined by fitting the data to models derived from perturbation theory. Measurements in water were used as the reference. For a 0% to 50% glycerol solution, an 87°–204° separation in the phase shift was observed. The slope of the plot of the phase shift vs (ηρ)0.5 was used to indicate the sensor’s sensitivity. The sensitivity for our 250 MHz quartz SH-SAW sensors was calculated to be 3.7×10−3m2sKg. The corresponding mass sensitivity was determined to be 9.25 × 105m2Kg. The limit of detection was calculated to be 36 picograms (pg), while the limit of quantification or LOQ was calculated to be 109 pg. Traditionally, liquid phase measurements have been challenging for SAW devices because liquids dampen the vibrating sensors severely. This problem has been largely solved using a transverse (shear) wave instead of the more popular longitudinal or Rayleigh waves. Liquid measurements are now possible using transverse waves, also known as shear waves, because transverse waves are only minimally attenuated by liquids. Shear-horizontal SAW sensors (SH-SAW) show great promise as label-free biosensors because of their ability to handle liquid samples. However, the viscosity of the liquid still induces loading effects and can be measured when the liquid is loaded onto the SH-SAW propagating surface (delay line). When the liquid above the delay line is perturbed by physical or chemical changes, such as binding to a receptor, it alters the propagating acoustic wave. The SH-SAW device can measure these changes in liquid properties as a change in the wave’s phase compared to the original wave. The device’s phase shift was recorded as a function of the changes in the density and viscosity of the binary glycerol solution and used to determine the sensitivity in the linear dynamic range of responses.

Materials Science↗

DOE OSTI2021

Divertor detachment in the pre-fusion power operation phase in ITER during application of resonant magnetic perturbations

Detachment of the divertor plasma during application of resonant magnetic perturbation (RMP) fields is evaluated for hydrogen H-mode plasma during the first pre-fusion power operation (PFPO-1) phase in ITER by 3D plasma boundary modelling with EMC3–EIRENE. Plasma response effects from a linearized, resistive, single fluid MHD model are discussed, which includes partial screening of the externally applied field—but also field amplification near the separatrix. Furthermore, this field amplification is found to play a pivotal role for the magnetic footprint on the divertor targets, but is sensitive to model parameters. Extensions of the footprint beyond the straight portions of the ITER vertical divertor targets, optimized for high stationary heat flux handling, may be possible depending on the level of toroidal rotation in the plasma. Exhaust from the bulk plasma is guided by the helical corrugations (lobes) of the perturbed separatrix, and this results in an upstream heat flux that is distributed over these lobes with lower peak values than in the typical radial heat flux profiles seen in the absence of magnetic perturbations. As a consequence, an earlier onset (with respect to the upstream density) of detachment is found in the traditional strike zone when RMPs are applied, but secondary, non-axisymmetric strike locations appear—and those remain attached at temperatures above 10 eV. Neon seeding can mitigate these non-axisymmetric heat loads, but this becomes less efficient for large magnetic footprints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Probing the structure of the initial state of heavy-ion collisions with p T -dependent flow fluctuations

The connection between initial-state geometry and anisotropic flow can be quantified through a well-established mapping between p T -integrated flow harmonics and cumulants of the initial transverse energy distribution. In this paper we successfully extend this mapping to also include p T -differential flow. In doing so, we find that subleading principal components of anisotropic flow can reveal previously unobserved details of the hydrodynamic response, in both the linear and the nonlinear regimes. Most importantly, we show that they provide novel information on the small-scale structures present in the initial stage of relativistic heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Prediction of Power Measurements Using Adaptive Filters

With the advent of smart grid concept, Internet of Things (IoT) and the deployment of smart meters, the cyberattack threats on power networks have increased due to the use of communication systems that can be accessed by adversaries. Attackers will have the ability to manipulate the outcomes of smart meters which in turn influence the core application of Energy Management System 9EMS): State Estimation (SE). Bad data analytic tools may fail to detect some attacks into measurements. Meanwhile, Machine Learning (ML) solutions have been proposed for detecting False Data Injection (FDI) attacks. However, there is a lack of ML time-series solutions presented in the state-of-the-art that is yet to be complex. In signal processing, time-series solutions do not only consider the signal, but also the statistics of the signal over time. Therefore, in this paper, a machine learning for time-series solutions is presented as an application to model the measurements of the power grid that are used in SE. The presented model takes into account adaptive linear and non-linear filters: Finite Impulse Response (FIR), and Infinite Impulse Response (IIR). The presented models are implemented and performed on the IEEE-118 bus system. The results indicate the advantage of applying those filters over the state-of-the-art machine learning solutions.

Hamad, Khaled↗

DOE OSTI2020

Sensitivity of surface temperature to oceanic forcing via q-Flux Green’s function experiments. Part III: asymmetric response to warming and cooling

Climate response is often assumed to be linear in climate sensitivity studies. However, by examining the surface temperature (TS) response to pairs of oceanic forcings of equal amplitude but opposite sign in a large set of local q-flux perturbation experiments with CAM5 coupled to a slab, we find strong asymmetry in TS responses to the heating and cooling forcings, indicating a strong nonlinearity intrinsic to the climate system examined. Regardless of where the symmetric forcing is placed, the cooling response to the negative forcing always exceeds the warming to the positive forcing, implying an intrinsic inclination towards cooling of our current climate. Thus, the ongoing global warming induced by increasing greenhouse gases may have already been alleviated by the asymmetric component of the response. The common asymmetry in TS response peaks in high latitudes, especially along sea ice edges, with notable seasonal dependence. Decomposition into different radiative feedbacks through a radiative kernel indicates that the asymmetry in the TS response is realized largely through lapse rate and albedo feedbacks. Further process interference experiments disabling the seasonal cycle and/or sea ice reveal that the asymmetry originates ultimately from the presence of the sea ice component and is further amplified by the seasonal cycle. The fact that a pair of opposite tropical q-flux forcings can excite very similar asymmetric response as a pair placed at 55°S strongly suggests the asymmetric response is a manifestation of an internal mode of the climate model system.

Liu, Fukai↗