SEARCH · Search NASA

Results for “Separation of variables”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Nemchinov-Dyson Solutions of the Two-Dimensional Axisymmetric Inviscid Compressible Flow Equations [Review]

We investigate the two-dimensional (2D) inviscid compressible flow equations in axisymmetric coordinates, constrained by an ideal gas equation of state (EOS). Beginning with the assumption that the 2D velocity field is space-time separable and linearly variable in each corresponding spatial coordinate, we proceed to derive an in nite family of elliptic or hyperbolic, uniformly expanding or contracting \gas cloud" solutions. Construction of specific example solutions belonging to this family is dependent on the solution of a system of nonlinear, coupled, second-order ordinary differential equations, and the prescription of an additional physical process of interest (e.g., uniform temperature or uniform entropy flow). The physical and computational implications of these solutions as pertaining to quantitative code verification or model qualification studies are discussed in some detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

A Probabilistic Autoencoder for Type Ia Supernova Spectral Time Series

We construct a physically parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of Type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an autoencoder that is interpreted probabilistically after training using a normalizing flow. We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multistage training setup alongside our physically parameterized network, we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including the automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an rms of 0.091 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Timing based clustering in the Northern Finland Birth Cohorts 1966 and 1986 suggests two new patterns for childhood BMI curve [Poster]

Childhood body mass index (BMI) is a widely used measure of adiposity in children (<18 years of age). Children grow with individual tempo and individuals of the same age, or of the same BMI, might be in different phases in their individual growth curves. Variability between different childhood BMI curves can be separated in two components: phase variability (x-axis; time) and amplitude variability (y-axis; BMI). Phase variability can be thought of arising from differences in maturational age between individuals. This is related to the timing of peaks and valleys in a child’s BMI curve.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

In situ X-ray spatial profiling reveals uneven compression of electrode assemblies and steep lateral gradients in lithium-ion coin cells

Coin cells are used extensively as test devices in battery research for evaluation of new materials and optimization of cycling protocols. Here, in situ X-ray diffraction profilometry is used to characterize spatial distribution of the active materials, lithiation, and phase distribution in electrodes of NCM523/graphite coin cells. The X-ray data indicate uneven areal compression of the electrode assembly in such cells, which we trace to a specific design feature that leads to elastic deformation of a metal spacer. Steep lithiation gradients observed in the electrodes imply radially-dependent resistivity, for which uneven compression of the separator is a likely cause. Electrochemical model calculations suggest that variable porosity of the polymer separator would account for the salient features of spatial profiles observed in these coin cells.

25 ENERGY STORAGE↗

DOE OSTI2024

Evaluation of partial volume correction and analysis of longitudinal [ 18 F]GTP1 tau PET imaging in Alzheimer's disease using linear mixed-effects models

Purpose: We evaluated the impact of partial volume correction (PVC) methods on the quantification of longitudinal [ 18 F]GTP1 tau positron-emission tomography (PET) in Alzheimer's disease and the suitability of describing the tau pathology burden temporal trajectories using linear mixed-effects models (LMEM). Methods: We applied van Cittert iterative deconvolution (VC), 2-compartment, and 3-compartment, and the geometric transfer matrix plus region-based voxelwise methods to data acquired in an Alzheimer's disease natural history study over 18 months at a single imaging site. We determined the optimal PVC method by comparing the standardized uptake value ratio change (%ΔSUVR) between diagnostic and tau burden–level groups and the longitudinal repeatability derived from the LMEM. The performance of LMEM analysis for calculating %ΔSUVR was evaluated in a natural history study and in a multisite clinical trial of semorinemab in prodromal to mild Alzheimer's disease by comparing results to traditional per-visit estimates. Results: The VC, 2-compartment, and 3-compartment PVC methods had similar performance, whereas region-based voxelwise overcorrected regions with a higher tau burden. The lowest within-subject variability and acceptable group separation scores were observed without PVC. The LMEM-derived %ΔSUVR values were similar to the per-visit estimates with lower variability. Conclusion: The results indicate that the tested PVC methods do not offer a clear advantage or improvement over non-PVC images for the quantification of longitudinal [ 18 F]GTP1 PET data. LMEM offers a robust framework for the longitudinal tau PET quantification with low longitudinal test–retest variability. Clinical trial registration: NCT02640092 and NCT03289143.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

A functional global sensitivity measure and efficient reliability sensitivity analysis with respect to statistical parameters

Sensitivity analysis and reliability assessment are two important aspects of structural and system safety. Epistemic uncertainty with respect to probabilistic model of input parameters due to lack of knowledge is present in many scarce-data applications and complicates the characterization of uncertainty in model response. In this article, we present two importance measures to evaluate the impact of distribution parameters on the probability distribution function (PDF) of the output and the failure probability. The epistemic uncertainty associated with the distribution parameters is modeled as random variables. Additionally, a modified extended polynomial chaos expansion (MEPCE) approach is introduced in which aleatory and epistemic random variables are modeled and propagated simultaneously while allowing the separate assessment for any single epistemic variable. A MEPCE-based kernel density estimation (KDE) construction provides a composite map from each epistemic variable to the response PDF. The functional global sensitivity index of the PDF with respect to the distribution parameters is thus derived, as a function of output, which is both more informative and more efficient than standard scalar sensitivity measures. Reliability sensitivity indices can be readily evaluated by integrating the global sensitivity index function over the failure zone. Three illustrative examples are used to demonstrate the proposed methodology.

42 ENGINEERING↗

DOE OSTI2020

Understanding How Ligand Functionalization Influences CO2 and N2 Adsorption in a Sodalite Metal–Organic Framework

In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal–organic frameworks (MOFs) known as Cu–BTT (BTT–3 = 1,3,5-benzenetristetrazolate) and Cu–BTTri (BTTri–3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu–BTTri compared to its tetrazole counterpart, Cu–BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu–BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu–BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn–Teller effect in Cu–BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials’ potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu–BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.

adsorption capacities↗

DOE OSTI2025

COMPASS-FME Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) Experiment Level 1 Sensor Data v2-1

This is the version 2-1 Level 1 (L1) data release for COMPASS-FME environmental sensors located at our Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experimental site. This manipulative, ecosystem-scale TEMPEST experiment addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in MD, USA. The experiment uses a large-unit (2000 m2), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments. L1 data are close to raw, but are units-transformed and have out-of-instrument-bounds, out-of-service, and outlier flags added. Duplicates and missing data are removed but otherwise these data are not filtered, and have not been subject to any additional algorithmic or human QA/QC. Any scientific analyses of L1 data should be performed with care. **This dataset will be updated quarterly with new data for the duration of the project** This dataset includes: - An overall dataset README file that describes the current version, gives citation and contact information, etc. - Site- and year-specific folders, each holding variable-specific CSV (comma separated value) data files for each site and plot in that year. - Metadata files within each site-year folder provide full information on data units, expected ranges, contact information, detailed flood times, as well as a general description of the site. - Environmental sensor types that appear in the data files include weather (ClimaVUE50, CS, RM Young, and LI instruments in the graphs below); soil conditions (TEROS12); soil redox state (Redox); groundwater variables (AquaTROLL200 and AquaTROLL600); open water sondes (Exo); tree sap velocity (Sapflow); and system voltage and state (Datalogger). Data are normally logged every 15 minutes. Please see v2-1 TEMPEST L1 Sensor Package Quick Start.pdf for detailed information on data package structure, temporal coverage, and versioning. The TEMPEST flood events occurred on the following dates. They lasted for ~10 hours each day and delivered ~80,000 gallons to each plot; many data streams are available at 1 or 5 minute frequency during these periods. * Tests: Aug 25 (fresh plot) and Sep 9 (salt plot), 2021 * TEMPEST 1: June 22, 2022 * TEMPEST 2: June 6-7, 2023 * TEMPEST 3: June 11-13, 2024 This dataset was updated 2026-03-12: (i) data now go through 2025-12-31 (previous end was 2025-06-30) and (ii) dataset and file names updated to “…v2-1” (previously was “v2-0”).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

The Time Domain Spectroscopic Survey: Changing-look Quasar Candidates from Multi-epoch Spectroscopy in SDSS-IV

Active galactic nuclei (AGN) can vary significantly in their rest-frame optical/UV continuum emission, and with strong associated changes in broad line emission, on much shorter timescales than predicted by standard models of accretion disks around supermassive black holes. Most such changing-look or changing-state AGN—and at higher luminosities, changing-look quasars (CLQs)—have been found via spectroscopic follow-up of known quasars showing strong photometric variability. The Time Domain Spectroscopic Survey of the Sloan Digital Sky Survey IV (SDSS-IV) includes repeat spectroscopy of large numbers of previously known quasars, many selected irrespective of photometric variability, and with spectral epochs separated by months to decades. Our visual examination of these repeat spectra for strong broad line variability yielded 61 newly discovered CLQ candidates. We quantitatively compare spectral epochs to measure changes in continuum and Hβ broad line emission, finding 19 CLQs, of which 15 are newly recognized. The parent sample includes only broad line quasars, so our study tends to find objects that have dimmed, i.e., turn-off CLQs. However, we nevertheless find four turn-on CLQs that meet our criteria, albeit with broad lines in both dim and bright states. We study the response of Hβ and Mg ii emission lines to continuum changes. The Eddington ratios of CLQs are low, and/or their Hβ broad line width is large relative to the overall quasar population. Repeat quasar spectroscopy in the upcoming SDSS-V black hole Mapper program will reveal significant numbers of CLQs, enhancing our understanding of the frequency and duty cycle of such strong variability, and the physics and dynamics of the phenomenon.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Impacts of stochastic coalescence variability on warm rain initiation using Lagrangian microphysics in box and large-eddy simulations

Various coalescence methods for Lagrangian microphysics schemes are tested in box and large-eddy simulation (LES) models, including the stochastic all-or-nothing super-droplet method (SDM) and a version of SDM (dSDM) that applies a fractional approach similar to the average impact method. In LES, variabilities driven by microphysics and by flow realizations are separated using the “piggybacking” technique. Rain initiation averaged over many realizations of the box model is delayed and rain variability increases as the number of super-drops per collision volume ( N SD ) is decreased using SDM. In contrast, rain initiation time using SDM in LES is insensitive to N SD for 32 ≤ N SD ≤ 512. This is explained through the interaction between LES grid boxes, each acting as a separate collision volume. Variability across the ensemble of LES collision volumes using SDM results in rain quickly initiating in some of the LES grid cells at low N SD and leading to a similar overall timing of rain initiation from the cloud compared to simulations with high N SD . There is a ~20% decrease in the total rain mass and mean rain flux as N SD is increased from 32 to 256, with little additional change as N SD is increased from 256 to 512. The fractional coalescence approach in dSDM leads to reduced microphysical variability and a 15-18 min delay in rain initiation compared to SDM. Finally, an additional LES ensemble with microphysical variability feeding back to the dynamics shows that flow variability dominates the impact of microphysical variability on rain properties. Thus, flow variability must be constrained to isolate impacts of microphysical variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Physical, socio-psychological, and behavioural determinants of household energy consumption in the UK

Determining which attitudes and behaviours predict household energy consumption can help accelerate the low-carbon energy transition. Conventional approaches in this domain are limited, often relying on survey methods that produce data on individuals’ motivations and self-reported activities without pairing these with actual energy consumption records, which are particularly hard to collect for large, nationally representative samples. This challenge precludes the development of empirical evidence on which attitudes and behaviours influence patterns of energy consumption, thus limiting the extent to which these can inform energy interventions or conservation programs. This study demonstrates a novel methodology for estimating energy consumption in the absence of actual energy records by using a large, publicly available data set of energy consumption in the UK. We develop a predictive model using the Smart Energy Research Laboratory (SERL) data portal (with records from nearly 13,000 UK households) and then use this model to predict energy consumption (both electric and gas) for a sample of 1,000 UK householders for which we separately collect over 200 variables relating to climate change attitudes and practices. Our approach uses a set of over 50 independent variables that are shared between the data sets, allowing us to train a model on the SERL data and use it to analyse the relationship between energy consumption and the opinions, motivations, and daily practices of survey respondents. Results show that electricity consumption is influenced by a broader range of factors compared to gas. Household energy use is best explained by physical dwelling characteristics, socio-demographic variables, and certain behavioural and attitudinal measures. Notably, pro-environmental attitudes, frugality, and conscientiousness correlate with lower energy use, while income and consumerism are linked to higher consumption. We discuss how these findings can inform efforts to decarbonise home energy use in the UK.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

DOE OSTI2026

Constraining the Aerosol Effects on Deep Convective Clouds by Considering the Coupling Between Clouds and the Planetary Boundary Layer

Several mechanisms have been proposed for the aerosol invigoration effect. Although their principles are well established, their actual magnitudes and roles in cloud development remain uncertain and debatable. This uncertainty partly stems from observational-based studies, in which it has been challenging to separate the co-variability between aerosols and meteorology. Addressing this problem requires large data samples. To this end, this study employs the Atmospheric Radiation Measurement data set expanding to 16 years (some up to 17 years, compared to 10 years in previous work) in the U.S. Southern Great Plains. It also conducts a more careful and rigorous analysis to isolate the influences of convective available potential energy (CAPE) and synoptic patterns to address a previously raised concern. We incorporated a new key process affecting aerosol-cloud interaction: cloud-surface coupling. The state/degree of the coupling relationship turns out to play an important role in the invigoration effect. Our analysis reinforces earlier findings of a robust positive relationship between cloud thickness and aerosol loading across CAPE percentiles—but only under cloud-surface coupled conditions. The increase in cloud thickness with aerosol loading is most pronounced in coupled clouds with high CAPE and bases below 1 km. Coupled clouds with bases below 1 km thicken between 1 and 4 km, depending on the CAPE percentile. Decoupled clouds show no such systematic changes. Synoptic patterns also lead to different strengths of the invigoration effect. Clean and polluted air masses are predominantly associated with northerly and southerly winds, respectively, with a stronger invigoration effect in cleaner air masses.

Geosciences↗

DOE OSTI2025

Attribution of the record-high 2023 SST using a deep-learning framework

Abstract The global-mean sea surface temperature (SST) reached a record high in 2023, exceeding the 2016 record by 0.14 °C. This unprecedented change in global-mean SST has major implications for our understanding of internal variability and the forced response in our changing climate. In this work, we use neural networks trained on simulated climate data to separate the contributions of internal variability and the forced response within observations. Performing attribution reveals that internal variability was responsible for +0.07 °C of the 2023 global mean SST, due to anomalously warm conditions in the Pacific, Atlantic, and Indian Ocean basins. Furthermore, these results provide a line of evidence for accelerated forced warming in recent years. Continued monitoring of the climate will be critical for understanding the drivers behind this unprecedented SST record.

Rader, Jamin K. (ORCID:0000000222045977)↗

DOE OSTI2022

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Model reduction for a power grid model

We examine the complexity of constructing reduced order models for subsets of the variables needed to represent the state of the power grid. In particular, we apply model reduction techniques to the DeMarco-Zheng power grid model. We show that due to the oscillating nature of the solutions and the absence of timescale separation between resolved and unresolved variables, the construction of accurate reduced models becomes highly non-trivial because one has to account for long memory effects. In addition, we show that a reduced model that includes even a short memory is drastically better than a memoryless model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Improved Air-Conditioning Demand Response of Connected Communities over Individually Optimized Buildings

Connected communities potentially offer much greater demand response capabilities over singular building energy management systems (BEMS) through an increase of connectivity. The potential increase in benefits from this next step in connectivity is still under investigation, especially when applied to existing buildings. This work utilizes EnergyPlus simulation results on eight different commercial prototype buildings to estimate the potential savings on peak demand and energy costs using a mixed-integer linear programming model. This model is used in two cases: a fully connected community and eight separate buildings with BEMS. The connected community is optimized using all zones as variables, while the individual buildings are optimized separately and then aggregated. These optimization problems are run for a range of individual zone flexibility values. The results indicate that a connected community offered 60.0% and 24.8% more peak demand savings for low and high flexibility scenarios, relative to individually optimized buildings. Energy cost optimization results show only marginally better savings of 2.9% and 6.1% for low and high flexibility, respectively.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2020

Emulating climate extreme indices

We use simple pattern scaling and time-shift to emulate changes in a set of climate extreme indices under future scenarios, and evaluate the emulators' accuracy. We propose a metric for the error in emulation in the context of initial condition ensembles, to specifically characterize the role of internal variability in the emulation performance. Our metric separates systematic emulation errors from unavoidable discrepancies between emulated and target values due to internal variability. We compute the metricis at grid-point scale, and we show geographically resolved results, or aggregate them at global scale. We demonstrate the use of our error metric in the emulation of a suite of temperature and precipitation extreme indices. We test and compare simple pattern scaling and time-shift using a range of trajectories spanning targets inspired by the Paris agreement -- warming to 1.5C and 2.0C from the pre-industrial baseline -- and two of the longer-established trajectories, RCP4.5 and RCP8.5. With this suite of scenarios we can test the effects on the performance of the size of the temperature gap between emulation origin and target. We find that for most indices emulation the dominant source of discrepancy is internal variability. For at least one index, however, counting exceedances of a high temperature threshold, significant portions of the globally aggregated discrepancy and its regional pattern originate from the systematic emulation error. This error exceeds internal variability of both the target and the emulated quantities in large coherent regions at low latitudes, and the explanation can be found in the differential behavior of temperature distributions across latitudes. The metric also highlights a fundamental difference in the two methods related to the simulation of internal variability, which is dampened significantly by simple pattern scaling. This aspect is of consequence when using these methods for specific applications, where preserving variability for uncertainty quantification is deemed important. With this study we offer our metric as a diagnostic tool, facilitating the formulation of scientific hypotheses on the reasons for the error. In the meantime, we show that for many impact relevant indices by now traditional emulation techniques can be accurate within the variations unavoidably introduced by internal variability, establishing the fundamental condition for using their emulation in impact modeling.

54 ENVIRONMENTAL SCIENCES↗