Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable importance”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Representing surface heterogeneity in land–atmosphere coupling in E3SMv1 single-column model over ARM SGP during summertime

Abstract. The Earth's land surface features spatial and temporal heterogeneity over a wide range of scales below those resolved by current Earth system models (ESMs). State-of-the-art land and atmosphere models employ parameterizations to represent their subgrid heterogeneity, but the land–atmosphere coupling in ESMs typically operates on the grid scale. Communicating the information on the land surface heterogeneity with the overlying atmospheric boundary layer (ABL) remains a challenge in modeling land–atmosphere interactions. In order to account for the subgrid-scale heterogeneity in land–atmosphere coupling, we implement a new coupling scheme in the Energy Exascale Earth system model version 1 (E3SMv1) that uses adjusted surface variances and covariance of potential temperature and specific water content as the lower boundary condition for the atmosphere model. The new lower boundary condition accounts for both the variability of individual subgrid land surface patches and the inter-patch variability. The E3SMv1 single-column model (SCM) simulations over the Atmospheric Radiation Measurement (ARM) Southern Great Plain (SGP) site were performed to assess the impacts. We find that the new coupling parameterization increases the magnitude and diurnal cycle of the temperature variance and humidity variance in the lower ABL on non-precipitating days. The impacts are primarily attributed to subgrid inter-patch variability rather than the variability of individual patches. These effects extend vertically from the surface to several levels in the lower ABL on clear days. We also find that accounting for surface heterogeneity increases low cloud cover and liquid water path (LWP). These cloud changes are associated with the change in cloud regime indicated by the skewness of the probability density function (PDF) of the subgrid vertical velocity. In precipitating days, the inter-patch variability reduces significantly so that the impact of accounting for surface heterogeneity vanishes. These results highlight the importance of accounting for subgrid heterogeneity in land–atmosphere coupling in next-generation ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kinetics Modeling for Design of Continuous Enzymatic Hydrolysis

Enzymatic hydrolysis of cellulose to monomeric sugars continues to be a limiting step in cost-effectively producing sugar and fermentation-based biofuels from biomass. In particular, the high cost of enzymes coupled with the long time-scale of reaction pose challenges to economic viability. Performing enzymatic hydrolysis in a continuous mode with enzyme recycle may provide a path towards substantially reduced costs for sugar production, but design and analysis are complicated by a lack of suitable kinetics models. Computationally attainable models, such as fractal-based models, require knowledge of the reaction-history of the biomass, and are thus only suitable for describing the batch reactions from which they are derived. Fundamental models, while potentially more generalizable, are often too computationally intensive to use in process or reactor modeling. In this work, a phenomenological rate model is proposed based on a two-phase substrate representation. Good agreement is seen between batch and continuous enzymatic hydrolysis (CEH) experiment data, which validates the model and enables us to solve for reactor design parameters, such as CEH reactor size and stream flow rates, based on process variables like yield. This model is integrated with techno-economic analysis software to explore economic sensitivities. Important design optimizations and tradeoffs are identified and quantified, including the relative cost imposed by rate slowdown from sugar inhibition versus the cost to remove and concentrate sugars at a lower concentration. It also identifies, high-leverage avenues for further exploration, such as increasing the maximum feasible solids concentration, and sustaining high membrane flux and reliability.

09 BIOMASS FUELS↗
DOE OSTI2025

Wind and Weather Variability within the Californian Offshore Wind Energy Areas

Weather variability over the Northeast Pacific (NEP) region and its influence on wind resources within the Californian offshore wind energy areas (WEAs) at Humboldt and Morro Bay are characterized using 20-years reanalysis model and satellite data. The hub-height (180 m) winds at both locations are predominantly northwesterly driven by the NEP high pressure system, with strong coastal gradients in surface pressure, fluxes, planetary boundary layer (PBL) depths and cloudiness. These sharp coastal gradients and strong annual cycles of temperature and moisture advections pose potential challenges in accurately modeling the local wind resource. Hub-height wind speeds and power capacity factors significantly vary for different regimes of PBL depths, surface fluxes and rain area fractions. This highlights the importance of studying the physical mechanisms driving these weather regimes, hence our analysis of how large-scale NEP weather variability drives the local meteorology at the WEAs. Furthermore, at both WEAs, PBL tops and cloud boundaries intersect the rotor layer (80-280 m) more than 30% and 20% of the time, respectively. While PBL depths significantly modulates hub-height winds and power, cloud boundaries do not have a similar impact, likely due to reanalysis errors in simulating cloud boundaries accurately. These findings underscore the challenges in deploying tall wind turbines in shallow cloudy boundary layers, where the interaction between clouds, precipitation, and atmospheric layers can impact turbine efficiency. As turbines grow taller and are deployed in more complex meteorological conditions, understanding these interactions is crucial for improving wind power forecasting and optimizing energy production in coastal regions.

17 WIND ENERGY↗
DOE OSTI2025

Streamflow Response to Glacier Mass Loss Varies With Basin Precipitation Across Alaska

Abstract Diminishing glaciers affect streamflow, and given the extent of glaciers in Alaska and adjacent Canada, continued glacier mass loss is likely to have profound effects on ecosystems sensitive to runoff. The effects of glacier mass loss on streamflow are likely to vary across the wide ranges of basin size, glacier cover, and precipitation in this region. In this study, we use U.S. Geological Survey (USGS) streamflow data with satellite‐based glacier volume change estimates to quantify how glacier mass loss subsidized streamflow over the 2000–2019 period for 116 glacierized basins. We examine interannual variability in that subsidy at three USGS‐monitored glaciers to explore the ability of the subsidy to buffer streamflow derived solely from precipitation. We found the relative importance of percent glacier cover on streamflow magnitude increases in drier basins. In the driest basins, glaciers produced 40 times greater percent glacier mass loss subsidies to streamflow for the percent glacier cover compared to the wettest basins. While the subsidy from glacier mass loss buffers interannual variability in streamflow to varying degrees, it can also increase streamflow variability. Smaller amounts of percent glacier cover are needed to produce summer‐melt‐dominated seasonal flow regimes in drier basins than in wetter basins. Decreasing glacier cover will eventually decrease summer streamflow, increasing spring streamflow in drier basins, and attenuating seasonality with increasing spring and autumnal streamflow in wetter basins. Quantifying the downstream effects of continued glacier mass loss without the computational expense of a hydrological model is broadly applicable in this changing climate.

Curran, Janet H. [U.S. Geological Survey Alaska Sc↗
DOE OSTI2022

Rain on snow (ROS) understudied in sea ice remote sensing: a multi-sensor analysis of ROS during MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate)

Abstract. Arctic rain on snow (ROS) deposits liquid water onto existing snowpacks. Upon refreezing, this can form icy crusts at the surface or within the snowpack. By altering radar backscatter and microwave emissivity, ROS over sea ice can influence the accuracy of sea ice variables retrieved from satellite radar altimetry, scatterometers, and passive microwave radiometers. During the Arctic Ocean MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition, there was an unprecedented opportunity to observe a ROS event using in situ active and passive microwave instruments similar to those deployed on satellite platforms. During liquid water accumulation in the snowpack from rain and increased melt, there was a 4-fold decrease in radar energy returned at Ku- and Ka-bands. After the snowpack refroze and ice layers formed, this decrease was followed by a 6-fold increase in returned energy. Besides altering the radar backscatter, analysis of the returned waveforms shows the waveform shape changed in response to rain and refreezing. Microwave emissivity at 19 and 89 GHz increased with increasing liquid water content and decreased as the snowpack refroze, yet subsequent ice layers altered the polarization difference. Corresponding analysis of the CryoSat-2 waveform shape and backscatter as well as AMSR2 brightness temperatures further shows that the rain and refreeze were significant enough to impact satellite returns. Our analysis provides the first detailed in situ analysis of the impacts of ROS and subsequent refreezing on both active and passive microwave observations, providing important baseline knowledge for detecting ROS over sea ice and assessing their impacts on satellite-derived sea ice variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Dual-Axis Radiographic Hydrodynamic Test Facility Capability Expansion (DCX) Strategy

The Dual-Axis Radiographic Hydrodynamic Test (DARHT) facility is a vital and important part of the Nation’s nuclear security enterprise. With more than two decades of operations, DARHT has a storied history. Conceived in the 1970s, constructed in the 1990s, and operational since 2000, DARHT has advanced from open-air hydrodynamic experiments (hydros) to foam-confined hydros, to vessel-confined hydros, and in 2022, the 75th hydro was successfully completed. Radiography has advanced from a single-axis, single-pulse system to a dual-axis, multi-pulse capability to variable fields of view (VFV) on both accelerators. The culmination of these experiences, accomplishments, and advancements has brought us to a very important question: What do the next two decades at DARHT look like? The world is not the same place it was in the 1990s when construction at DARHT was in progress. Evolving threats, an expanding mission, and technology changes necessitate adaptation. To adapt, the aging facility, accelerators, vessels, and detector systems require improvements to ensure DARHT remains the Nation’s hydrodynamic data foundation for stockpile certification, safety, surety, and global security threats.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

A data-driven sensor placement approach for detecting voltage violations in distribution systems

Stochastic fluctuations in power injections from distributed energy resources (DERs) combined with load variability can cause constraint violations (e.g., exceeded voltage limits) in electric distribution systems. To monitor grid operations, sensors are placed to measure important quantities such as the voltage magnitudes. Here, in this paper, we consider a sensor placement problem which seeks to identify locations for installing sensors that can capture all possible violations of voltage magnitude limits. We formulate a bilevel optimization problem that minimizes the number of sensors and avoids false sensor alarms in the upper level while ensuring detection of any voltage violations in the lower level. This problem is challenging due to the nonlinearity of the power flow equations and the presence of binary variables. Accordingly, we employ recently developed conservative linear approximations of the power flow equations that overestimate or underestimate the voltage magnitudes. By replacing the nonlinear power flow equations with conservative linear approximations, we can ensure that the resulting sensor locations and thresholds are sufficient to identify any constraint violations. Additionally, we apply various problem reformulations to significantly improve computational tractability while simultaneously ensuring an appropriate placement of sensors. Lastly, we improve the quality of the results via an approximate gradient descent method that adjusts the sensor thresholds. We demonstrate the effectiveness of our proposed method for several test cases, including a system with multiple switching configurations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

pH dependence of brown-carbon optical properties in cloud water

Abstract. Light-absorbing organic species present in aerosols, collectively called brown carbon (BrC), have important but highly uncertain effects on climate. Clouds likely represent a significant medium for secondary BrC production and for bleaching reactions, though the relative importance of the formation and loss processes in clouds is unknown at present. The acidity (or pH) of atmospheric particles and clouds affects the optical properties of BrC and bleaching rates. Given the wide variability of pH in the atmosphere (pH in particles and clouds ranges from −1 to 8), the optical properties of BrC and its bleaching behavior are expected to vary significantly, and the link between pH and BrC is yet another uncertainty in attempts to constrain its climate forcing effects. In this work, we characterize the pH dependence of BrC optical properties – including light absorption at 365 nm (Abs365), the mass absorption coefficient (MAC365), and the absorption Ångström exponent (AAE) – in bulk cloud water sampled from the summit of Whiteface Mountain, NY. In all samples (n=17), Abs365 and MAC365 increased linearly with increasing pH, highlighting the importance of reporting pH in studies of BrC in aqueous media. There was strong variability in the sensitivity of Abs365 to pH, with normalized slopes that ranged from 5.1 % to 17.2 % per pH unit. The normalized slope decreased strongly with increasing cloud water [K+], suggesting that the non-biomass-burning BrC has optical properties that are more sensitive to pH than BrC associated with biomass burning. AAE also showed a distinct pH dependence as it was relatively flat between pH 1.5–5 and then decreased significantly above pH 5. The cloud water composition was used to inform thermodynamic predictions of aerosol pH upwind and/or downwind of Whiteface Mountain and the subsequent changes in BrC optical properties. Overall, these results show that, in addition to secondary BrC production, photobleaching, and the altitudinal distribution, the climate forcing of BrC is quite strongly affected by its pH-dependent absorption.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2024

The trailing edge is everywhere: tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species’ geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data—such that individuals at the cool edge of a species’ distribution should benefit from warming (the “leading edge”), whereas individuals at the warm edge should suffer (the “trailing edge”). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species’ distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming—the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species’ distribution with respect to temperature and the majority of the species’ distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species’ geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data—such that individuals at the cool edge of a species’ distribution should benefit from warming (the “leading edge”), whereas individuals at the warm edge should suffer (the “trailing edge”). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species’ distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming—the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species’ distribution with respect to temperature and the majority of the species’ distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Deep Learning Coordinate-Free Quantum Chemistry

Computing quantum chemical properties of small molecules and polymers can provide insights valuable to physicists, chemists and biologists when designing new materials, catalysts, biological probes and drugs. Deep learning can compute quantum chemical properties accurately in a fraction of the time required by commonly used methods such as density functional theory (DFT). However, many of these deep learning architectures require energy minimized molecular geometries as input, which is also computationally expensive, and decreasing the reproducibility and throughput of these methods. In this study, we demonstrate that accurate quantum chemical computations can be performed without optimized geometries by operating in the coordinate-free domain using deep learning on graph encodings. Furthermore, we also find that the choice of graph-encoding architecture substantially affects the performance of these methods. The Wave architecture outperforms graph convolution architectures, particularly on complex molecules. Furthermore, the structures of these graph encoding architectures provide an opportunity to probe an important, outstanding question in quantum mechanics: What types of quantum chemical properties can be represented by local-variable models? We find that Wave, a local-variable model, is more accurately calculates quantum chemical properties. Graph convolutional architectures require global variables, and are not as effective as as Wave. We anticipate that coordinate-free, deep-learning models of quantum chemistry will become valuable tools in chemistry and biology, enabling researchers to rapidly screen chemical databases or identify new molecules using automated, de-novo design algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Microphysics Dominates Sub‐Cloud Rain Evaporation in Trade Cumuli Over Barbados

This study investigates the microphysical and thermodynamic influences on North Atlantic trade-cumulus rain evaporation in the sub-cloud layer over Barbados from January to November 2020. Using radar observations and a 1-D model, the results reveal that microphysical properties, namely geometric mean diameter (D g ) and raindrop concentration, primarily control evaporation processes under conditions of limited thermodynamic variability. Thermodynamic factors like cloud base height and relative humidity exert a smaller influence due to their narrow observed ranges. Small D g leads to rapid increases in rain evaporation fraction (REF), as small drops evaporate completely, while larger D g yields slower, more gradual REF changes. Frequent weak-rain cases reduce the mean evaporation flux, whereas infrequent intense-rain events substantially enhance it. These findings highlight the pivotal role of microphysical variability in shaping sub-cloud evaporation in this trade cumulus regime, while thermodynamic effects may be more important in other environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Intramolecular hydrogen bonding analysis

The quasi-atomic orbital (QUAO) bonding analysis is used to study intramolecular hydrogen bonding (IMHB) in salicylic acid and an intermediate that is crucial to the synthesis of aspirin. The bonding analysis rigorously explores IMHB through directly accessing information that is intrinsic to the molecular wave function, thereby bypassing the need for intrinsically biased methods. The variables that affect the strength of IMHB are determined using kinetic bond orders, QUAO populations, and QUAO hybridizations. Important properties include both the interatomic distance between hydrogen and oxygen participating in the IMHB and the hybridization on the oxygen. The bonding analysis further shows that each intramolecular hydrogen bond is a four-electron three-center bond. Here, the bonding analysis is used to understand how aromatic reactivity changes due to the effect of functional groups on the aromatic ring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Uncertainty Quantification for Capacity Expansion Planning

This report quantifies the uncertainty in output decisions from a Capacity Expansion Planning (CEP) model. The need to understand how uncertainties within CEP models and modeling assumptions affect Quantities of Interest (QoIs) such as expansion and operating costs, as well as expansion decisions remains an ongoing challenge in scientific research and industrial operations. This area of research is particularly important for models which seek to capture how large networks will evolve and operate under increased sources of variable generation, i.e., higher penetration of renewable technologies such as solar and wind generators. Uncertainty quantification (UQ) of CEP models which estimate expansion costs and decisions, and production cost models which estimate operating costs and dispatch decisions, is a key focus of research at NREL. The Regional Energy Deployment System (ReEDS) represents a state-of-the-art CEP model and considers a range of possible grid evolutions in an attempt to identify key drivers, ramifications, and decisions which contribute to better informed investment and policy decisions. However, research to quantify how uncertainties and model assumptions, such as unit commitment (UC), within ReEDS may be affecting its outputs remains challenging due to to size and complexity of the model

24 POWER TRANSMISSION AND DISTRIBUTION↗