SEARCH · Engineering Papers

Results for “k mean”

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 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗

DOE OSTI2022

EdgeAI: Machine learning via direct attached accelerator for streaming data processing at high shot rate x-ray free-electron lasers

We present a case for low batch-size inference with the potential for adaptive training of a lean encoder model. We do so in the context of a paradigmatic example of machine learning as applied in data acquisition at high data velocity scientific user facilities such as the Linac Coherent Light Source-II x-ray Free-Electron Laser. We discuss how a low-latency inference model operating at the data acquisition edge can capitalize on the naturally stochastic nature of such sources. We simulate the method of attosecond angular streaking to produce representative results whereby simulated input data reproduce high-resolution ground truth probability distributions. By minimizing the mean-squared error between the decoded output of the latent representation and the ground truth distributions, we ensure that the encoding layers and resulting latent representation maintains full fidelity for any downstream task, be it classification or regression. We present throughput results for data-parallel inference of various batch sizes, some with throughput exceeding 100 k images per second. We also show in situ training below 10 s per epoch for the full encoder–decoder model as would be relevant for streaming and adaptive real-time data production at our nation’s scientific light sources.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Dynamic Wind Loading on CSP Collectors Caused by Turbulent Wind Fluctuations: Insights from a 2-Year Field Campaign

Concentrating Solar Power (CSP) is a promising solar technology for electricity generation with thermal energy storage and with the additional benefit of industrial heat production. Wind loading on CSP collector structures, such as parabolic troughs or heliostats, is one of the primary drivers of their structural design costs. In particular, dynamic wind loading is a major source of uncertainty in the collector design process, which heavily relies on wind tunnel testing. In the field, the turbulent nature of the incoming wind creates fluctuating loads (support structure loads and resulting mirror deflections) on the collectors, with impacts on fatigue lifetime and optical performance. As is well known, wind tunnel tests cannot entirely reproduce the complex turbulent wind conditions typically observed at full-scale plants. To shed light on this topic, NREL initiated a field campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs. We use the published 2-year dataset of high-resolution combined wind and structural loads measurements [1] to characterize the dynamic structural wind response. For quantifying dynamic wind loading, we apply the concept of admittance functions, which are spectral transfer functions that couple the turbulent wind to resulting structural loads (aerodynamic admittance), and to the structural response (mechanical admittance). In practice, aerodynamic admittance describes which turbulent eddy sizes are effective in creating structural loads. The mechanical admittance describes in which frequency ranges these loads are reinforced or dampened by the structure. While these functions are an established concept in civil engineering, their recent application to a single full-scale heliostat [2] proved their broader applicability to CSP collectors. Here, we present a characterization of admittance functions for full-scale parabolic trough collectors and show how wind characteristics (mean wind speed and direction, turbulent kinetic energy, turbulent length scales), the sun-tracking trough angle, and row position alter the admittance functions. Further, we study to which extent the admittance functions are universal for a specific trough geometry and how our findings compare to reported heliostat results. References [1] https://data.openei.org/submissions/5938. [2] Blume, K., Roger, M., and Pitz-Paal, R. 2023b. "Simplified analytical model to describe wind loads and wind-induced tracking deviations of heliostats." Solar Energy, 256, 96-109. https://doi.org/10.1016/j.solener.2023.03.055.

admittance functions↗

DOE OSTI2021

Rotational Dynamics of an Amphidynamic Zirconium Metal–Organic Framework Determined by Dielectric Spectroscopy

A zirconium metal–organic framework with a difluorophenylene rotator bearing a permanent electric dipole of ~3.2 D was synthesized, and its rotational motion was analyzed by temperature- and frequency-dependent broadband dielectric spectroscopy. While solid-state NMR confirms fast rotation qualitatively, the dissipation factors measured between 113 and 153 K suggested an activation energy E a = 2.6 kcal/mol, but deviations from a single Debye relaxation suggested a dynamic process that cannot be accounted for by a well-defined potential with a single activation barrier. The dynamic heterogeneity of the dipolar rotor was confirmed by analysis in terms of a Cole–Cole relaxation, which suggested a mean barrier of ~1.9 kcal/mol, with a heterogeneity that decreases as temperature increases. Based on the single-crystal structure, we propose that the kinetic heterogeneity results from a temperature-dependent potential where rotation motion is mediated by the escape of the rotator from an energy well created by a double Ph–H···F–Ph hydrogen bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Evaluating Climate Models’ Cloud Feedbacks Against Expert Judgment

Abstract The persistent and growing spread in effective climate sensitivity (ECS) across global climate models necessitates rigorous evaluation of their cloud feedbacks. Here we evaluate several cloud feedback components simulated in 19 climate models against benchmark values determined via an expert synthesis of observational, theoretical, and high‐resolution modeling studies. We find that models with smallest feedback errors relative to these benchmark values generally have moderate total cloud feedbacks (0.4–0.6 W m −2 K −1 ) and ECS (3–4 K). Those with largest errors generally have total cloud feedback and ECS values that are too large or too small. Models tend to achieve large positive total cloud feedbacks by having several cloud feedback components that are systematically biased high rather than by having a single anomalously large component, and vice versa. In general, better simulation of mean‐state cloud properties leads to stronger but not necessarily better cloud feedbacks. The Python code base provided herein could be applied to developmental versions of models to assess cloud feedbacks and cloud errors and place them in the context of other models and of expert judgment in real‐time during model development.

58 GEOSCIENCES↗

DOE OSTI2024

Air‐Ice‐Ocean Coupling During a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions Across Scales

Abstract Arctic cyclones are key drivers of sea ice and ocean variability. During the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, joint observations of the coupled air‐ice‐ocean system were collected at multiple spatial scales. Here, we present observations of a strong mid‐winter cyclone that impacted the MOSAiC site as it drifted in the central Arctic pack ice. The sea ice dynamical response showed spatial structure at the scale of the evolving and translating cyclonic wind field. Internal ice stress and ocean stress play significant roles, resulting in timing offsets between the atmospheric forcing and the ice response and post‐cyclone inertial ringing in the ice and ocean. Ice motion in response to the wind field then forces the upper ocean currents through frictional drag. The strongest impacts to the sea ice and ocean from the passing cyclone occur as a result of the surface impacts of a strong atmospheric low‐level jet (LLJ) behind the trailing cold front and changing wind directions between the warm‐sector LLJ and post cold‐frontal LLJ. Impacts of the cyclone are prolonged through the coupled ice‐ocean inertial response. Local impacts of the approximately 120 km wide LLJ occur over a 12 hr period or less and at scales of a kilometer to a few tens of kilometers, meaning that these impacts occur at combined smaller spatial scales and faster time scales than most satellite observations and coupled Earth system models can resolve.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Field, frequency, and temperature dependencies of the surface resistance of nitrogen diffused niobium superconducting radio frequency cavities

We investigate the rf performance of several single-cell superconducting radio-frequency cavities subjected to low temperature heat treatment in nitrogen environment. The cavities were treated at temperature 120 – 165 ° C for an extended period of time (24–48 h) either in high vacuum or in a low partial pressure of ultrapure nitrogen. The improvement in Q 0 with a Q rise was observed when nitrogen gas was injected at ∼ 300 ° C during the cavity cooldown from 800 ° C and held at 165 ° C , without any degradation in accelerating gradient over the baseline performance. The treatment was applied to several elliptical cavities with frequency ranging from 0.75 to 3.0 GHz, showing an improved quality factor as a result of low temperature nitrogen treatments. The Q rise feature is similar to that achieved by nitrogen alloying Nb cavities at higher temperature, followed by material removal by electropolishing. The surface modification was confirmed by the change in electronic mean free path and tuned with the temperature and duration of heat treatment. The decrease of the temperature-dependent surface resistance with increasing rf field, resulting in a Q rise, becomes stronger with increasing frequency and decreasing temperature. The data suggest a crossover frequency of ∼ 0.95 GHz above that the Q rise phenomenon occurs at 2 K. Some of these results can be explained qualitatively with an existing model of intrinsic field-dependence of the surface resistance with both equilibrium and nonequilibrium quasiparticle distribution functions. The change in the Q slope below 0.95 GHz may result from masking contribution of trapped magnetic flux to the residual surface resistance. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Optical phonon dominated heat transport: A first-principles thermal conductivity study of BaSn 2

Acoustic phonons with long mean free paths have long been believed to control the lattice thermal conductivity κ L in solids dominantly. In this study, however, we demonstrate an optical phonon dominated κ L in BaSnS 2 . By solving the Peierls-Boltzmann transport equation, we predict a low diagonal lattice thermal conductivity κ L (D) of 0.34 W m -1 K -1 at 850 K, which is less than half the κ L (D) of SnS at the same temperature. Further calculations following the Allen-Feldman model suggest the additional off-diagonal lattice thermal conductivity κ L (OD) contributed by wavelike tunneling phonons. The κ L (OD) becomes pronounced at the high temperature ( 0.17 W m -1 K -1 at 850 K) and leads to a deviation of the temperature dependence of κ L from T -1 to T -0.76 , suggesting the potential lattice anharmonicity in BaSnS 2 . Further analyses indicate BaSnS 2 has over 68% of κ L contributed by optical phonons. We show this uncommon optical phonon dominated κ L is due to the relatively high group velocities of optical phonons in BaSnS 2 . The phonon mode visualization suggests these relatively high-velocity optical phonons correspond to the antiphase vibrations in BaSnS 2 monolayers, which is originated from the unique permutation of SnS 3 tetrahedra. Finally, by investigating the mode-resolved group velocity, relaxation time, and Grüneisen parameter, we attribute the intrinsic low κ L of BaSnS 2 to the soft lattice and the relatively high lattice anharmonicity induced by the Ba-S weak bonding and Sn(II) lone-pair electrons. Our study explicitly analyzes the microscopic mechanism of optical phonon dominated heat transport in BaSnS 2 and suggests it worthy of further experimental studies as an intrinsic low-κ L material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Hanford low-activity waste glass composition-temperature-melt viscosity relationships

This study developed a model for predicting viscosity of alkali-alumino-borosilicate glass melts as functions of composition and temperature. The model is based on a total of 3935 viscosity-temperature data from 574 glasses with viscosity values ranging from 2.53 to 7260 Poise (P) in the temperature range of 900–1260°C. Several different model forms were surveyed, including those based on Arrhenius, Vogel-Fulcher-Tammann, Avramov-Milchev, and Mauro-Yue-Ellison-Gupta-Allen. For each of these models, combinations of the temperature-independent parameters were fitted to composition. It was found that generally fitting more than one temperature-independent parameter as functions of composition resulted in overfitting. The Avramov-Milchev-based model was found to best represent the Hanford low-activity waste glass melt viscosity data based on model fit and validation statistics. A 21-term partial quadratic mixture model was recommended for use. This model predicts melt viscosity with a root-mean square error of .1736 ln(P), which is similar to the error in viscosity measurements from replicate glass analyses of .1383 ln(P). Viscosity was found to be most increased by SiO 2 > Al 2 O 3 > ZrO 2 > SnO 2 and most decreased by Li 2 O > Na 2 O > B 2 O 3 > CaO > K 2 O > MgO, at temperatures from 900 to 1260°C.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

An Assessment of Persistent Acoustic Monitoring of a Nuclear Reactor during Full Power Generation

Persistent low-frequency (<180 Hz) acoustic detection took place within the boundaries of Oak Ridge National Laboratory to monitor full power operations of the High Flux Isotope Reactor. Three acoustic sensors were installed at distances of 69, 101, and 914 m from the northeast corner of the cooling towers to monitor and assess four reactor power generation cycles. Features were extracted from power spectral density calculations where data were collected during reactor on and off operations. Diverse spectral features were present during full reactor power, including a 21.4 Hz fundamental frequency and ascending harmonics. Using bandpass filters, these related frequencies were isolated and summed, and the root mean square energy was calculated. The method of isolating and summing characteristic features provided a significant improvement in identifying acoustic behavior related to reactor power when the raw signals were obscured by noise.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

No Evidence for Orbital Clustering in the Extreme Trans-Neptunian Objects

The apparent clustering in longitude of perihelion piv and ascending node Ω of extreme trans-Neptunian objects (ETNOs) has been attributed to the gravitational effects of an unseen 5–10 Earth-mass planet in the outer solar system. To investigate how selection bias may contribute to this clustering, we consider 14 ETNOs discovered by the Dark Energy Survey, the Outer Solar System Origins Survey, and the survey of Sheppard and Trujillo. Using each survey's published pointing history, depth, and TNO tracking selections, we calculate the joint probability that these objects are consistent with an underlying parent population with uniform distributions in $\bar{w}$ and Ω. We find that the mean scaled longitude of perihelion and orbital poles of the detected ETNOs are consistent with a uniform population at a level between 17% and 94% and thus conclude that this sample provides no evidence for angular clustering.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Correlating and Simulating Socio-Demographically Driven Residential End-Use Activity Schedules

Incorporating socio-demographic and behavioral considerations into decision-support tools is crucial for identifying gaps and addressing consumer needs to ensure reliable and affordable energy solutions. In energy simulation models, the correlation between socio-demographics and time-use behavior is not well-captured. Thus, we developed a large-scale simulation workflow to generate schedules for 10 residential activities across 24 population segments defined by age, income, and employment status. Using pre-pandemic 2015-2019 American Time Use Survey (ATUS) data, we used ANOVA to confirm the correlation between demographic factors and time use. We explored three k-modes clustering methods-backward, forward, and a new hybrid approach-to delineate the occupancy patterns based on demographics. Using the probability of cluster membership for each population segment and a time inhomogeneous Markov chain to generate activity transition probabilities for each cluster, we simulated 50,000 schedules per segment and validated them against the ATUS data. The hybrid method produced the most socio-demographically differentiated clusters while demonstrating comparable performance to other approaches, with an overall root mean square error of 0.12 for both weekday and weekend schedules. Thus, the hybrid method, where each cluster is dominated by certain demographic segments and occupancy patterns, offers more modeling versatility in terms of scenario analysis. The new workflow improves the socio demographic differentiation of energy consumption by considering differences in time use. This approach enables future research on demographically segmented time of use (TOU) energy consumption, including impacts of TOU utility bills and rate analysis, long-run marginal emissions, and energy retrofits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Boron coordination change in barium borate melts and glasses and its contribution to configurational heat capacity, entropy, and fragility

High-energy x-ray diffraction from molten and glassy BaB 2 O 4 and BaB 4 O 7 has been performed using aerodynamic levitation and laser heating over a wide range of temperatures. Remarkably, even in the presence of a heavy metal modifier dominating x-ray scattering, it was possible to extract accurate values for the tetrahedral, sp 3 , boron fraction, N 4 , which declines with increasing temperature, using bond valence-based mapping from the measured mean B–O bond lengths while accounting for vibrational thermal expansion. These are used within a boron-coordination-change model to extract enthalpies, ΔH, and entropies, ΔS, of isomerization between sp 2 and sp 3 boron. The results for BaB 4 O 7 , ΔH = 22(3) kJ mol -1 boron, ΔS = 19(2) J mol -1 boron K -1 , agree quantitatively with those found previously for Na 2 B 4 O 7 . Analytical expressions for N 4 (J, T) and associated configurational heat capacity, $C_{P}^{conf}$(J, T), and entropy, S conf (J, T), contributions are extended to cover a wide composition range 0 ≤ J = BaO/B 2 O 3 ≤ 3 using a model for ΔH(J) and ΔS(J) derived empirically for lithium borates. Maxima in the $C_{P}^{conf}$(J, T g ) and fragility index contributions are thereby predicted for J ≲ 1, higher than the maximum observed and predicted in N 4 (J, T g ) at J ≃ 0.6. We discuss the utility of the boron-coordination-change isomerization model in the context of borate liquids containing other modifiers and the prospect of neutron diffraction to aid in empirical determinations of modifier-dependent effects, illustrated by new neutron diffraction data on Ba 11 B 4 O 7 glass, its well-known α-polymorph, and lesser-known δ-phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Analysis of a Heat Balance Experiment in the TREAT Reactor with SPARTA and OpenMC

This article presents the interpretation of a heat balance experiment in TREAT using the SPARTA methodology and an OpenMC model. The study examines the effects of control rods and temperature variations on the effective multiplication factor ($k_{\rm eff}$) and correction factors. Results indicate that both core temperature and rod positions have a linear effect on $k_{\rm eff}$ and correction factors in the upper-west channel. A critical hypothesis of SPARTA was validated: reactivities and correction factors can be obtained by linearly combining individual effects. The previously considered negligible effect of temperature on correction factors was found to be significant, necessitating updates to the SPARTA methodology. By fitting an exponential curve to mean temperature data, SPARTA's accuracy and convergence have been improved. Comparisons between SPARTA and OpenMC showed consistent reactivity removal due to temperature increase. While SPARTA effectively interprets complex transients, simpler methods may suffice when correction factors have minimal influence on reactivity calculations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Asteroseismological analysis of the polluted ZZ Ceti star G 29 – 38 with TESS

G29 – 38 (TIC 422526868) is one of the brightest (V = 13.1) and closest (d = 17.51 pc) pulsating white dwarfs with a hydrogen-rich atmosphere (DAV/ZZ Ceti class). It was observed by the TESS spacecraft in sectors 42 and 56. The atmosphere of G 29 – 38 is polluted by heavy elements that are expected to sink out of visible layers on short time-scales. Here, the photometric TESS data set spans ~51 d in total, and from this, we identified 56 significant pulsation frequencies, that include rotational frequency multiplets. In addition, we identified 30 combination frequencies in each sector. The oscillation frequencies that we found are associated with g-mode pulsations, with periods spanning from ~ 260 to ~ 1400 s. We identified rotational frequency triplets with a mean separation δν ℓ = 1 of 4.67 μHz and a quintuplet with a mean separation δν ℓ = 2 of 6.67 μHz, from which we estimated a rotation period of about 1.35 ± 0.1 d. We determined a constant period spacing of 41.20 s for ℓ = 1 modes and 22.58 s for ℓ = 2 modes. We performed period-to-period fit analyses and found an asteroseismological model with M * /M ⊙ = 0.632 ± 0.03, K, and log g = 8.048 ± 0.005 (with a hydrogen envelope mass of M H ~ 5.6 × 10 –5 M * ), in good agreement with the values derived from spectroscopy. We obtained an asteroseismic distance of 17.54 pc, which is in excellent agreement with that provided by Gaia (17.51 pc).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Effective Missing Value Imputation Methods for Building Monitoring Data

To understand behaviors of natural and man-made events, such as energy consumption of buildings, which accounts for 40% of energy uses in the US, we deploy automated monitoring devices to record periodic observations. However, such experimental and observation data often contains problems and irregularities that have to be cleaned up before analyses. Due to various conditions affecting sensor operations, the communication channels, recording steps, or the recording media, the recorded data might have missing values, errors, or anomalous values. An effective way to clean up these problems is to replace these missing values, errors and anomalous values with expected values, a process generally known as imputation. In this work, we survey commonly used missing value imputation techniques and compare their performance on a set of building monitoring data. To compare the different types of sensor measurements with widely varying characteristics, we use normalized root mean squared error (NRMSE) as the key metric for the effectiveness of the imputation methods. We additionally consider periodicity and run time when considering comparing methods. Through extensive testing, we find that for small gap sizes, up to 8 consecutive missing values, linear interpolation performs the best; for larger gaps stretching up to 48 consecutive missing values, K-nearest neighbors provides the most accurate imputations; for even larger gaps, more computational intensive methods, such as matrix factorization, achieve the smallest NRMSE. Additionally, we observe that these computationally intensive algorithms not only provide accurate imputations for large gaps, but are also more robust across all types of sensors.

Cho, B↗

DOE OSTI2025

Work In Progress: Using Internships as Means for Indirect Assessment of ABET Criteria 3 “1-7” Student Outcomes

The Electrical & Computer Engineering (ECE) Department at Portland State University (PSU) has developed a Power Engineering Internship (PEI) program that provides engineering career development pathways within the electric utility industry. The PEI is supported by several U.S. Department of Energy grants that aim to develop quality career opportunities and develop a future electric utility workforce that represents the nation’s diverse populations. The PSU ECE Department intends to use surveys of internship participants as assessment tools for its ABET accreditation process, in particular, the ABET Criteria 3 Student Outcomes (SOs). SOs relate to the knowledge, skills, and behaviors that students acquire as they progress through an engineering program [1]. They describe what engineering students are expected to know and should be able to do by the time they graduate. Engineering education literature presents multiple methods for using internships as means for assessing SOs. Methods include using self-assessments that derive from internship experiences as well as mapping survey responses from employers and interns to SOs. Such methods gain value as they are assessed over multiple years. However, the majority of this literature discusses assessment tools that were developed for the ABET 2000 Criteria 3 “a through k” SOs. One of our motivations is to provide assessment tools that can be used to assess modern Criteria 3 SOs based on students’ internship experiences.

Bass, Robert [Portland State University]↗

DOE OSTI2023

Explaining Forcing Efficacy With Pattern Effect and State Dependence

The magnitude of global surface temperature change in response to unit radiative forcing depends on the type and magnitude of forcing agent—a concept known as a “forcing efficacy.” However, the mechanisms behind the forcing efficacy are still unclear. In this study, we perform a set of simulations using CESM1 to calculate the efficacy of 10 different forcing agents defined in terms of fixed-SST effective radiative forcing, and then use a Green's function approach to show that each forcing efficacy can be largely understood in terms of the radiative feedbacks associated with the different surface temperature patterns induced by the forcing agents (a pattern effect). We also quantify how the state dependence of feedbacks on global mean surface temperature anomalies impacts forcing efficacies. The results show that the forcing efficacy can be well reconstructed with a combination of pattern effect and state dependence.

58 GEOSCIENCES↗