Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial 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 1,459 records · Page 81

A guide to the BRAIN Initiative Cell Census Network data ecosystem

Characterizing cellular diversity at different levels of biological organization and across data modalities is a prerequisite to understanding the function of cell types in the brain. Classification of neurons is also essential to manipulate cell types in controlled ways and to understand their variation and vulnerability in brain disorders. The BRAIN Initiative Cell Census Network (BICCN) is an integrated network of data-generating centers, data archives, and data standards developers, with the goal of systematic multimodal brain cell type profiling and characterization. Emphasis of the BICCN is on the whole mouse brain with demonstration of prototype feasibility for human and nonhuman primate (NHP) brains. Here, we provide a guide to the cellular and spatial approaches employed by the BICCN, and to accessing and using these data and extensive resources, including the BRAIN Cell Data Center (BCDC), which serves to manage and integrate data across the ecosystem. We illustrate the power of the BICCN data ecosystem through vignettes highlighting several BICCN analysis and visualization tools. Finally, we present emerging standards that have been developed or adopted toward Findable, Accessible, Interoperable, and Reusable (FAIR) neuroscience. The combined BICCN ecosystem provides a comprehensive resource for the exploration and analysis of cell types in the brain.

59 BASIC BIOLOGICAL SCIENCES↗

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the equivalent Dancoff factor cell method performing slowing down calculations on the fly for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for pin-cell, unit-cell, and fuel assembly problems with various resonance self-shielding conditions based on very high temperature reactor, high temperature test reactor, and Empire benchmark cores, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also showed that the methodology worked well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

TRACER-MAP Field Campaign Report

During July-August 2022, the U.S. Department of Energy (DOE) Atmospheric System Research (ASR)-funded Tracking Aerosol Convection Interactions Experiment (TRACER)-MAP campaign completed an Atmospheric Radiation Measurement (ARM) user facility field study to measure and map aerosol, volatile organic compounds (VOC), trace gas, and select meteorological observations across Houston, Texas. That ASR project and ARM field campaign, TRACER-MAP, was designed to complement measurements made at the first ARM Mobile Facility (AMF1) by the Aerosol Observing System (AOS) during the TRACER campaign. In this manner, TRACER-MAP effectively extended aerosol measurements from the TRACER AMF deployment, increasing the spatial coverage across the Houston metropolitan area, characterizing a greater diversity of source mixtures (e.g., industrial, traffic, residential, and biogenic) and capturing air masses with differential aging of urban source emissions (e.g., downtown versus downwind). That rich data set now offers many possibilities for in-depth analysis that align with the ASR mission to “research aerosol processes that affect Earth’s radiative balance and hydrological cycle.”

54 ENVIRONMENTAL SCIENCES↗

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the on-the-fly slowing down method for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Among various on-the-fly slowing down methods, the equivalent Dancoff factor cell method was employed. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for the very high temperature reactor (VHTR), high temperature test reactor (HTTR), and Empire benchmark problems with various resonance self-shielding conditions, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also show that the methodology works well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Discovery top-down proteomics in symbiotic soybean root nodules

Proteomic methods have been widely used to study proteins in complex biological samples to understand biological molecular mechanisms. Most well-established methods (known as bottom-up proteomics, BUP) employ an enzymatic digestion step to cleave intact proteins into smaller peptides for liquid chromatography (LC) mass spectrometry (MS) detection. In contrast, top-down proteomics (TDP) directly characterizes intact proteins including all possible post-translational modifications (PTMs), thus offering unique insights into proteoform biology where combinations of individual PTMs may play important roles. We performed TDP on soybean root nodules infected by the symbiotic Bradyrhizobium japonicum in both the wildtype bacterium and a nifH- mutant, which lacks the ability to fix nitrogen in the soybean root nodule. TDP captured 1648 proteoforms derived from 313 bacterial genes and 178 soybean genes. Leghemoglobin, the most abundant protein in the sample, existed in many truncated proteoforms. Interestingly, these truncated proteoforms were considerably more abundant in the wildtype relative to the nifH- mutant, implicating protease activity as an important factor in nitrogen fixation. Proteoforms with various PTMs and combinations thereof were identified using an unrestricted open modification search. This included less common PTMs such as myristoylation, palmitoylation, cyanylation, and sulfation. In parallel, we collected high resolution MS imaging (MSI) data of intact proteins and biopolymers (<20 kDa due to current technical limitations) from sections of the soybean root nodules using matrix-assisted laser desorption/ionization (MALDI) coupled to high resolution Orbitrap. Several detected proteoforms exhibited unique spatial distributions inside the infection zone and cortex, suggesting functional compartmentalization in these regions. A subset of peaks from the MALDI-MSI were assigned to proteoforms detected in TDP LCMS data based on matching accurate masses. Many of the proteins detected in both LCMS and MALDI-MSI are currently uncharacterized in UniProt: the PTM and spatial information presented here will be valuable in understanding their biological functions. Taken together, our study demonstrates how untargeted TDP approach can provide unique insights into plant proteoform biology. On-going technology developments are expected to further improve TDP coverage for more comprehensive high-throughput analysis of proteoforms.

59 BASIC BIOLOGICAL SCIENCES↗

Broadband Multi-wavelength Properties of M87 during the 2017 Event Horizon Telescope Campaign

In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass ∼6.5 × 109 M ⊙. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87’s spectrum. We can exclude that the simultaneous γ-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the γ-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded.

79 ASTRONOMY AND ASTROPHYSICS↗

Improvements to the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.

97 MATHEMATICS AND COMPUTING↗

Robust analysis of space-, time-, and energy-resolved soft x-ray measurements of magnetically confined fusion plasmas (invited)

A novel compact multi-energy soft x-ray (ME-SXR) diagnostic based on the PILATUS3 100K x-ray detector has been developed in collaboration between the Princeton Plasma Physics Laboratory and the University of Wisconsin-Madison and tested on the Madison Symmetric Torus (MST) reversed-field pinch. This solid-state photon-counting detector consists of a two-dimensional array of ~100 000 pixels for which the lower photon absorption cutoff energy can be independently set, allowing it to be configured for a unique combination of simultaneous spatial, spectral, and temporal resolution of ~1 cm, 100 eV, and 500 Hz, respectively. The diagnostic is highly versatile and can be readily adapted to diverse plasma operating conditions and scientific needs without any required downtime. New results from improved-confinement and quasi-single helicity plasmas in the MST demonstrate how the detector can be applied to study multiple aspects of the evolution of magnetically confined fusion-grade plasmas. These include observing the evolution of thermal emissivity, characterizing the energy of mid-Z excitation lines, extracting the T e profile, and observing the evolution of non-thermal populations. A technique for integrating the ME-SXR diagnostic into an integrated data analysis framework based on Bayesian inference is also presented. As a result, this allows ME-SXR measurements to be combined with data for complementary diagnostics in order to simultaneously infer Te and nZ from all available information.

47 OTHER INSTRUMENTATION↗

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗

A TCN-Based Hybrid Forecasting Framework for Hours-Ahead Utility-Scale PV Forecasting

This paper presents a Temporal Convolutional Network (TCN) based hybrid PV forecasting framework for enhancing hours-ahead utility-scale PV forecasting. The hybrid framework consists of two forecasting models: a physics-based trend forecasting (TF) model and a data-driven fluctuation forecasting (FF) model. Three TCNs are integrated in the framework for: i) blending the inputs from different Numerical Weather Prediction sources for the TF model to achieve superior performance on forecasting hourly PV profiles, ii) capturing spatial-temporal correlations between detector sites and the target site in the FF model to achieve more accurate forecast of intra- hour PV power drops, and iii) reconciling TF and FF results to obtain coherent hours-ahead PV forecast with both hourly trends and intra-hour fluctuations well preserved. To automatically identify the most contributive neighboring sites for forming a detector network, a scenario-based correlation analysis method is developed, which significantly improves the capability of the FF model on capturing large power fluctuations caused by cloud movements. Here, the framework is developed, tested, and validated using actual PV data collected from 95 PV farms in North Carolina. Simulation results show that the performance of 6 hours ahead PV power forecasting is improved by 20% - 30% compared with state-of-the-art methods.

42 ENGINEERING↗

Analysis of a Prototype Multi-Detector Fast-Neutron Radiography Panel

A multi-detector fast neutron radiography panel was built using the previous work on scalable neutron radiography using the IDEAS ROSSPAD readout module. A new aluminum housing was built to accommodate a large number of detectors tiled together. Additional changes to startup and processing code were made to operate the detector as one cohesive unit. Spatial resolution of the full panel using Cs-137 gammas was reported to be 0.42 line pairs per centimeter at 90% MTF and 2.09 line pairs per centimeter at 10% MTF. Here, three neutron radiographs generated using a Cf-252 fission neutron source were used to determine the spatial resolution of the panel for neutrons. The experiments had 90% MTF values of 0.24, 0.3, and 0.27 line pairs per centimeter and 10% MTF values of 1.30, 1.46, and 1.40 line pairs per centimeter. An example neutron radiograph was also used to prove that the radiography panel can perform true neutron radiography.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Meteoric 10Be Flux Calibration Data for the East River Watershed, Colorado, USA

This data package contains tabular and geospatial data used to quantify and model meteoric beryllium-10 fluxes in the East River watershed, Colorado, USA. The tabular component includes calibration-site data from five glacial moraine sites and includes environmental variables used to evaluate spatial controls on meteoric 10Be delivery, including elevation, mean annual precipitation (MAP), mean snow depth, and mean snow water equivalent (SWE). These site-level data were used to compare observed fluxes with environmental gradients across the watershed and to evaluate the effects of erosion correction on flux estimates. The package also includes supporting slope and curvature values used to assess topographic inputs to the erosion analysis. A second component of the data package contains updated manuscript tables and regression outputs used to summarize the relationships between meteoric 10Be flux and environmental predictors. These tables include meteoric 10Be sample information and AMS results, site-level environmental values, site-level meteoric 10Be inventory and flux values, watershed-averaged predicted fluxes, soil bulk density measurements, fine-fraction values, soil pH measurements, and regression statistics including slope, intercept, coefficient of determination, and p-value. The regression products include both standard linear regressions and regressions in which the intercept is constrained to pass through zero, and they support the analyses presented in the companion manuscript. Together, these tabular files provide the numerical basis for the manuscript tables and the regression-based interpretation of meteoric 10Be flux variability in a snow-dominated mountain watershed. The geospatial component of the package consists of GeoTIFF raster files used to generate the map products presented in Figures 2 and 6 of the companion manuscript. These rasters represent watershed-scale spatial layers for environmental variables and regression-based predictions of meteoric 10Be flux. This dataset contains comma-separated values files (.csv), Microsoft Excel files (.xlsx), GeoTIFF raster files (.tif), and upporting metadata files, including CSV data dictionaries and readme text files (.csv, .txt). The tabular files can be opened with standard spreadsheet software, and the raster files can be viewed and analyzed in GIS software such as ArcGIS Pro or QGIS. Together, these files document the numerical and spatial datasets used to calibrate and predict meteoric 10Be delivery in the East River watershed.

East River↗

Multi‐year hydroclimatic droughts and pluvials across the conterminous United States

Abstract Time series of water‐year runoff for 2,109 hydrologic units (HUs) across the conterminous United States (CONUS) for the 1900 through 2014 period were used to identify drought and pluvial (i.e., wet) periods. Characteristics of the drought and pluvial events including frequency, duration, and severity were examined and compared. Additionally, a similar analysis was performed using gridded tree‐ring reconstructions of the Palmer Drought Severity Index (PDSI) for the period 1475 through 2005 to place the drought and pluvial characteristics determined using water‐year runoff for 1900 through 2014 in the context of multi‐century climate variability. The temporal and spatial variability of droughts and pluvials determined using runoff for the 1900 through 2014 period indicated that most drought events in the CONUS occurred before about 1970, whereas most pluvial periods occurred after about 1970. This change in the frequencies of drought and pluvial events around 1970 was largely related to an increase in fall (October through December) precipitation across much of the central United States. Also, the duration and severity of droughts and pluvials identified using runoff for the 1900 through 2014 period generally were not significantly different from the drought and pluvial characteristics identified using the PDSI for the 1475 through 2005 period.

McCabe, Gregory J.↗

Non-singlet quark helicity $\mathrm{PDF}$s of the nucleon from pseudo-distributions

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional z 2 contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional z 2 contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale μ 2 . The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned z 2 contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A contactless in situ EFISH method for measuring electrostatic potential profile of semiconductor/electrolyte junctions

In photoelectrochemical cells, promising devices for directly converting solar energy into storable chemical fuels, the spatial variation of the electrostatic potential across the semiconductor–electrolyte junction is the key parameter that determines the cell performance. In principle, electric field induced second harmonic generation (EFISH) provides a contactless in situ spectroscopic tool to measure the spatial variation of electrostatic potential. However, the total second harmonic generation (SHG) signal contains the contributions of the EFISH signals of semiconductor space charge layer and the electric double layer, in addition to the SHG signal of the electrode surface. The interference of these complex quantities hinders their analysis. In this work, to understand and deconvolute their contributions to the total SHG signals, bias-dependent SHG measurements are performed on the rutile TiO2(100)–electrolyte junction as a function of light polarization and crystal azimuthal angle (angle of the incident plane relative to the crystal [001] axis). A quadratic response between SHG intensity and the applied potential is observed in both the accumulation and depletion regions of TiO2. The relative phase difference and amplitude ratio are extracted at selected azimuthal angles and light polarizations. At 0° azimuthal angle and s-in–p-out polarization, the SHG intensity minimum has the best match with the TiO2 flatband potential due to the orthogonal relative phase difference between bias-dependent and bias-independent SHG terms. We further measure the pH-dependent flatband potential and probe the photovoltage under open circuit conditions using the EFISH technique, demonstrating the capability of this contactless method for measuring electrostatic potential at semiconductor–electrolyte junctions.

Chemistry↗

The distribution of dust in edge-on galaxies: I. The global structure

ABSTRACT In this first paper in a series we present a study of the global dust emission distribution in nearby edge-on spiral galaxies. Our sample consists of 16 angularly large and 13 less spatially resolved galaxies selected from the DustPedia sample. To explore the dust emission distribution, we exploit the Herschel photometry in the range 100–500 $\mu $m. We employ Sérsic and 3D disc models to fit the observed 2D profiles of the galaxies. Both approaches give similar results. Our analysis unequivocally states the case for the presence of extraplanar dust in between 6 and 10 large galaxies. The results reveal that both the disc scale length and height increase as a function of wavelength between 100 and 500 $\mu $m. The dust disc scale height positively correlates with the dust disc scale length, similar to what is observed for the stellar discs. We also find correlations between the scale lengths and scale heights in the near- and far-infrared which suggest that the stellar discs and their dust counterparts are tightly connected. Furthermore, the intrinsic flattening of the dust disc is inversely proportional to the maximum rotation velocity and the dust mass of the galaxy: more massive spiral galaxies host, on average, relatively thinner dust discs. Also, there is a tendency for the dust-to-stellar scale height ratio to decrease with the dust mass and rotation velocity. We conclude that low-mass spiral galaxies host a diffuse, puffed-up dust disc with a thickness similar to that of the stellar disc.

79 ASTRONOMY AND ASTROPHYSICS↗

A study on the impact of mid-gap defects on vertical GaN diodes

GaN is a favorable martial for future efficient high voltage power switches. GaN has not dominated the power electronics market due to immature substrate, homoepitaxial growth, and immature processing technology. Understanding the impact of the substrate and homoepitaxial growth on the device performance is crucial for boosting the performance of GaN. Here, we studied vertical GaN PiN diodes that were fabricated on non-homogenous Hydride Vapor Phase Epitaxy (HVPE) substrates from two different vendors. We show that defects which stemmed from growth techniques manifest themselves as leakage hubs. Different non-homogenous substrates showed different distribution of those defects spatially with the lesser quality substrates clustering those defects in clusters that causes pre-mature breakdown. Energetically these defects are mostly mid-gap around 1.8Ev with light emission spans from 450nm to 700nm. Photon emission spectrometry and hyperspectral electroluminescence were used to locate these defects spatially and energetically.

42 ENGINEERING↗

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗