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

Spiner vs EOSPAC6: capabilities, performance, and accuracy considerations

The Cross-Cutting Capabilities Project (XCAP) has significant interest in collecting and implementing a common set of libraries to be used across the various production hydro-codes at LANL. A common set of libraries will help facilitate comparisons between codes. In addition to minimizing variables for comparisons and physics validation, efforts toward optimization and porting of packages for future architectures will be more efficient. As part of this push Singularity is being considered as an inclusive materials interface library. One crucial part of this library is the importing, inverting, and interpolation of equation of state (EOS) data. Currently the Lagrangian Applications Project (LAP) and Safety Applications Project (SAP) are using EOSPAC6 for delivery EOS data. The Eulerian Applications Project (EAP) has historically used TEOS and more recently an implementation of Singularly via a package called Spiner for EOS interpolation. Singularity currently has the option of using several analytical EOS models, directly employing SESAME via Spiner, or directly using SESAME via EOSPAC6 (not optimized yet). Assuming Singularity moves forward as a common platform for the implementation of materials models, a decision will need to be made on which EOS interpolation package XCAP should move forward with, given resources and people are finite. We will attempt to address pros and cons of Spiner and EOSPAC6 and the trade-offs that should be considered during the decision making process. This report is intended as an ASC-PEM-EOS perspective on what is needed in an EOS interpolation package. Performance and Accuracy sections will mostly address data from comparison studies in LA-UR-22-22699.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2021

Real-World Driving Features for Identifying Intelligent Driver Model Parameters

Driver behavior models play a significant role in representing different driving styles and the associated relationships with traffic patterns and vehicle energy consumption in simulation studies. The models often serve as a proxy for baseline human driving when assessing energy-saving strategies that alter vehicle velocity. Such models are especially important in connectivity-enabled energy-saving strategy research because they can easily adapt to changing driving conditions like posted speed limits or change in traffic light state. While numerous driver models exist, parametric driver models provide the flexibility required to represent variability in real-world driving through different combinations of model parameters. These model parameters must be informed by a representative set of parameter values for the driver model to adequately represent a real-world driver. It stands to reason that determining the parameter values from real-world driving data would serve the purpose of representing a real-world driver. Although the published literature is replete with techniques and consequences of using real-world driving data to determine parameter values for parametric driver models, none have explored them in the context of using shorter driving features where the parameter values may change over the course of a single trip for the same driver. In this study we consider the “intelligent driver model” (IDM) as our driver behavior model and use real-world driving data from the Transportation Secure Data Center (TSDC) maintained at the National Renewable Energy Laboratory (NREL). The TSDC includes real-world travel data from across the United States, from which NREL has created a wide range of driving routes consisting of road features such as speed limits, stop locations, and turn locations. The real-world driving data are categorized into different driving regimes and extracted into driving features. The driving features are then used to calibrate the parameter values for the IDM. The distribution of the parameters and the relationships among them are reported. The insights obtained from this study enable judicious usage of IDM or similar parametric driver models to represent baseline human driver behavior in simulations.

27 ARPA - Advanced Research Projects Agency-Energy↗
DOE OSTI2021

Real-World Driving Features for Identifying Intelligent Driver Model Parameters: Preprint

Driver behavior models play a significant role in representing different driving styles and the associated relationships with traffic patterns and vehicle energy consumption in simulation studies. The models often serve as a proxy for baseline human driving when assessing energy-saving strategies that alter vehicle velocity. Such models are especially important in connectivity-enabled energy-saving strategy research because they can easily adapt to changing driving conditions like posted speed limits or change in traffic light state. While numerous driver models exist, parametric driver models provide the flexibility required to represent variability in real-world driving through different combinations of model parameters. These model parameters must be informed by a representative set of parameter values for the driver model to adequately represent a real-world driver. It stands to reason that determining the parameter values from real-world driving data would serve the purpose of representing a real-world driver. While the published literature is replete with techniques and consequences of using real-world driving data to determine parameter values for parametric driver models, none have explored them in the context of using shorter driving features where the parameter values may change over the course of a single trip for the same driver. In this study we consider the “Intelligent Driver Model” (IDM) as the driver behavior model to explore and real-world driving data from the Transportation Secure Data Center (TSDC) maintained at the National Renewable Energy Laboratory (NREL). The TSDC includes real-world travel data from across the United States, from which NREL has created a wide range of driving routes consisting of road features such as speed limits, stop locations and turn locations. The real-world driving data are categorized into different driving regimes and extracted into driving features. The driving features are then used to calibrate the parameter values for the IDM. The distribution of the parameters and the relationships among them are reported. The insights obtained from this study enable judicious usage of IDM or similar parametric driver models to represent baseline human driver behavior in simulations.

27 ARPA - Advanced Research Projects Agency-Energy↗
DOE OSTI2020

Exploring the Use of DSCOVR/EPIC Satellite Observations to Monitor Vegetation Phenology

Vegetation phenology plays a pivotal role in regulating several ecological processes and has profound impacts on global carbon exchange. Large-scale vegetation phenology monitoring mostly relies on Low-Earth-Orbit satellite observations with low temporal resolutions, leaving gaps in data that are important for monitoring seasonal vegetation phenology. High temporal resolution satellite observations have the potential to fill this gap by frequently collecting observations on a global scale, making it easier to study change over time. This study explored the potential of using the Earth Polychromatic Imaging Camera (EPIC) onboard the Deep Space Climate Observatory (DSCOVR) satellite, which captures images of the entire sunlit face of the Earth at a temporal resolution of once every 1–2 h, to observe vegetation phenology cycles in North America. We assessed the strengths and shortcomings of EPIC-based phenology information in comparison with the Moderate-resolution Imaging Spectroradiometer (MODIS), Enhanced Thematic Mapper (ETM+) onboard Landsat 7, and PhenoCam ground-based observations across six different plant functional types. Our results indicated that EPIC could capture and characterize seasonal changes of vegetation across different plant functional types and is particularly consistent in the estimated growing season length. Our results also provided new insights into the complementary features and benefits of the four datasets, which is valuable for improving our understanding of the complex response of vegetation to global climate variability and other disturbances and the impact of phenology changes on ecosystem productivity and global carbon exchange.

Weber, Maridee A.↗
DOE OSTI2023

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

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

Bioenergy Feedstock Library Annual Summary Report 2024

The Bioenergy Feedstock Library (BFL), part of the Biomass Feedstock National User Facility (BFNUF) located at Idaho National Laboratory (INL), is a physical sample repository and a web-accessible electronic database. The BFL stores physical and chemical characteristics of biomass and waste carbon sources for energy use, as well as samples generated from U.S. Department of Energy (DOE) Bioenergy Technologies Office (BETO) and U.S. Department of Agriculture-funded projects. The objective of this Bioenergy Feedstock Library Annual Summary Report for 2024, similar to the 2023 Annual Summary Report , is to focus on the updates to: (1) publicly available analytical data and equipment tracked through the BFNUF, (2) significant increases in the physical samples available for request, (3) sample and data archival progress from recent BETO-funded projects, and (4) publicly available data sets created upon request from BETO, INL projects, or outside entities compared to the previous annual summary reports. This report highlights key statistics and available data and information important for INL, BFL users, academics, and industry.

09 BIOMASS FUELS↗
DOE OSTI2022

Identifications of RR Lyrae Stars and Quasars from the Simulated Data of Mephisto-W Survey

We have investigated the feasibilities and accuracies of the identifications of RR Lyrae stars and quasars from the simulated data of the Multi-channel Photometric Survey Telescope (Mephisto) W Survey. Based on the variable sources light curve libraries from the Sloan Digital Sky Survey (SDSS) Stripe 82 data and the observation history simulation from the Mephisto-W Survey Scheduler, we have simulated the uvgriz multi-band light curves of RR Lyrae stars, quasars and other variable sources for the first-year observation of Mephisto W Survey. We have applied the ensemble machine learning algorithm Random Forest Classifier (RFC) to identify RR Lyrae stars and quasars, respectively. We build training and test samples and extract ~150 features from the simulated light curves and train two RFCs respectively for the RR Lyrae star and quasar classification. We find that, our RFCs are able to select the RR Lyrae stars and quasars with remarkably high precision and completeness, with purity = 95.4% and completeness = 96.9% for the RR Lyrae RFC and purity = 91.4% and completeness = 90.2% for the quasar RFC. In conclusion, we have also derived relative importances of the extracted features utilized to classify RR Lyrae stars and quasars.

(galaxies:) quasars: general↗
DOE OSTI2023

Experimental and computational design tools for industrial drying processes: Challenges in process‐limit prediction

Abstract The coating and drying of inks and slurries are important steps to manufacture a plethora of products. Drying processes, particularly, comprise energy‐intensive steps that affect product cost and quality. Prior work has highlighted failures of various multicomponent diffusivity models to conserve mass in dryer modeling and challenges in predicting process limits given variability in published values of key thermodynamic parameters. Herein, we develop a computational model and benchtop drying experiments to investigate these concerns for drying polymer‐laden coatings. Model predictions of process limits in a single‐zone drying oven demonstrate that published variability in Flory–Huggins parameter yields large variations in predicted operating temperatures above which blistering occurs. This indicates that caution should be exercised when choosing approaches to obtain or predict the Flory–Huggins parameter, and that both benchtop drying experiments and a set of additional experiments, such as sorption experiments, are needed to fully characterize and optimize a given drying process.

Parrish, Chance↗
DOE OSTI2025

Precision DIS thrust predictions for HERA and EIC

We present predictions for the DIS 1-jettiness event shape $τ^b_1$, or DIS thrust, using the framework of Soft Collinear Effective Theory (SCET) for factorization, resummation of large logarithms, and rigorous treatment of nonperturbative power corrections, matched to fixed-order QCD away from the resummation region. Our predictions reach next-to-next-to-next-to-leading-logarithmic (N 3 LL) accuracy in resummed perturbation theory, matched to $\mathcal{O}(α^2_s)$ fixed-order QCD calculations obtained using the program NLOJet++. We include a rigorous treatment of hadronization corrections, which are universal across different event shapes and kinematic variables x and Q at leading power, and supplement them with a systematic scheme to remove $\mathcal{O}$(Λ QCD ) renormalon ambiguities in their definition. The framework of SCET allows us to connect smoothly the nonperturbative, resummation, and fixed-order regions, whose relative importance varies with x and Q, and to rigorously estimate theoretical uncertainties, across a broad range of x and Q covering existing experimental results from HERA as well as expected new measurements from the upcoming Electron- Ion-Collider (EIC). Our predictions will serve as an important benchmark for the EIC program, enabling the precise determination of the QCD strong coupling α s and the universal nonperturbative first moment parameter Ω 1 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Southern Himalayas rainfall as a key driver of interannual variation of pre-monsoon aerosols over the Tibetan Plateau

The Tibetan Plateau (TP) is one of the most climate-sensitive regions around the world. Aerosols imported from adjacent regions reach their peak during the pre-monsoon season and play a vital role in the TP environment. However, the strong interannual variation in aerosols transported to the TP has not been fully understood. Here, we show that the interannual variability of pre-monsoon aerosols transported to the TP is influenced more by rainfall over the southern Himalayas than near-surface wind. Rainfall modulates fire events and biomass burning emissions and reduces aerosols over the TP by wet scavenging. Contrary to the role of wind in increasing aerosol transport, the positive correlation between wind and aerosols in the TP reported in previous studies is contributed by the negative interannual correlations between wind and rainfall and between rainfall and fire events over the southern Himalayas. This study highlights the co-variability of wind and rainfall and their confounding impacts on aerosols in the southern Himalayas and over the TP. With pre-monsoon rainfall projected to increase in adjacent regions of southern TP, aerosol transport to the TP may be mitigated in the future.

Liu, Weichien↗
DOE OSTI2020

Structural basis for divergent and convergent evolution of catalytic machineries in plant aromatic amino acid decarboxylase proteins

Radiation of the plant pyridoxal 5'-phosphate (PLP)-dependent aromatic L -amino acid decarboxylase (AAAD) family has yielded an array of paralogous enzymes exhibiting divergent substrate preferences and catalytic mechanisms. Plant AAADs catalyze either the decarboxylation or decarboxylation-dependent oxidative deamination of aromatic L -amino acids to produce aromatic monoamines or aromatic acetaldehydes, respectively. These compounds serve as key precursors for the biosynthesis of several important classes of plant natural products, including indole alkaloids, benzylisoquinoline alkaloids, hydroxycinnamic acid amides, phenylacetaldehyde-derived floral volatiles, and tyrosol derivatives. In this paper, we present the crystal structures of four functionally distinct plant AAAD paralogs. Through structural and functional analyses, we identify variable structural features of the substrate-binding pocket that underlie the divergent evolution of substrate selectivity toward indole, phenyl, or hydroxyphenyl amino acids in plant AAADs. Moreover, we describe two mechanistic classes of independently arising mutations in AAAD paralogs leading to the convergent evolution of the derived aldehyde synthase activity. Applying knowledge learned from this study, we successfully engineered a shortened benzylisoquinoline alkaloid pathway to produce (S)-norcoclaurine in yeast. This work highlights the pliability of the AAAD fold that allows change of substrate selectivity and access to alternative catalytic mechanisms with only a few mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Climate models show colorado drying sooner and with greater certainty east of the continental divide

Many studies have examined the aridity of the Colorado River Basin and the possible impacts of climate change which could further strain already over-allocated water resources in the region. Fewer studies have examined the multiple Colorado Rocky Mountain headwater regions specifically. This is especially true of areas East of the Continental Divide, despite water originating there being critical to cities and agriculture in Eastern Colorado and further downstream. This paper explores and compares drying trends in the Eastern and Western Colorado Rocky Mountains using single-model initial-condition large ensembles from ten global climate models. The use of multiple models allows us to identify signals that are consistent across different physics parameterizations, model grids, and other model intricacies. The large ensembles also allow us to quantify the time of emergence of these climate change signals--that is, when did (or when will) the long term change due to anthropogenic greenhouse gasses exceed the internal variability of the climate system. Consistent with previous studies, we find evidence of drying on both sides of the Continental Divide. That drying is more pronounced, occurs sooner, and is more consistent across global climate models in the East, however, highlighting the region’s importance despite generally receiving less attention than the West.

Rugg, Allyson↗
DOE OSTI2022

Marine Renewable Energy Applications for Restorative Ocean Farming: Kelp

Kelp farming and kelp forest restoration have both been proposed as a solution to locally decrease the impacts of ocean acidification and eutrophication, often with co-benefits to other forms of aquaculture and mariculture. Compared to global markets, the kelp industry in the United States is still in its early phases, with the first commercial kelp farm founded in Casco Bay, Maine in 2010. Since then, interest and effort in kelp production has been increasing, with farms now present in Maine, New Hampshire, Connecticut, Rhode Island, Massachusetts, New York, Washington, and Alaska. Many research projects are underway in the United States to explore benefits of 3D ocean farming, tackle logistical problems of working in the ocean, autonomous farming, and explore viable end uses for kelp products. In seaweed farming, to remove the stored carbon or excess nutrients from the system, the biomass needs to be harvested at the optimal time to avoid the release of CO 2 that comes with decomposition. Timing of the harvest is also important for maximum crop yield, which can vary based on the final product. Additional monitoring needs can include a variety of water quality metrics, growth measurements, and visuals to ensure the health of the farm, comply with permits, support operations and maintenance functions. The variables measured may vary by desired end use of the product, location of farm, and operational design. Monitoring all of these parameters requires specialized devices that can be costly and challenging to maintain. Monitoring devices often face power and logistical constraints that could prevent kelp farmers from adopting these technologies or receiving accurate, efficient monitoring to assess ecosystem benefits and valuation. Marine energy has been identified as a possible power source for these devices. This project investigates the power needs for conducting kelp farm environmental monitoring compared with the available marine energy resource to evaluate if locally generated ocean energy could provide a solution to these monitoring challenges and benefit kelp farmers. This process was structured as follows: 1. Define what data is needed for farmers and their communities through desk research and interviews with end users. 2. Identify sensors and power requirements currently in use or available for commercial purchase. 3. Analyze current kelp and other mariculture farm locations for the potential marine energy resource. 4. Analyze farm designs and associated structures to make recommendations for marine energy design. 5. Quantify value that investment in sensors could provide in terms of carbon credit possibilities.

09 BIOMASS FUELS↗
DOE OSTI2023

Experimental Evaluation of Interference in 2.4 GHz Wireless Network

To attain automation across different applications, nuclear power plants are beginning to leverage advancements in wireless communication technologies. A “one-size-fits-all” solution cannot be applied since wireless technologies are selected according to application needs, quality of service requirements, and economic restrictions. To balance the trade-off between technical and economic requirements, a multi-band heterogeneous wireless network architecture is needed. Numerous wireless technologies including Wi-Fi, Zigbee, and Bluetooth share the 2.4 GHz industrial, scientific, and medical band. However, due to different channel access mechanisms and transmit power levels, and very importantly, uncoordinated use, coexistence of these devices in the same vicinity can cause interference and degradation in performance. This report provides the technical basis for understanding the coexistence of these wireless technologies through an experimental evaluation of their performance. This report investigates interactions encompassing variables such as transmission power level, distance between the devices, data rates, and the utilization of co-channel or adjacent channels. The results show that the operation of both Zigbee and Bluetooth is severely compromised when coexisting with Wi-Fi within the same frequency spectrum. On the other hand, the performance of Bluetooth is not impaired by Zigbee and vice versa unless there exists any external interference from Wi-Fi.

2.4 GHz↗
DOE OSTI2025

Selection of Global Climate Model Data for Downscaling With Generative Machine Learning and Use in the Power Planning for Alignment of Climate and Energy Systems Project

The range of results from climate models and scenarios is important to the understanding of uncertainty in power planning analysis. A U.S. Department of Energy-funded analytic project called Power Planning for Alignment of Climate and Energy Systems is developing data and analytic methods to reflect the effects of climate change on key variables for power system planning, as part of the Grid Modernization Lab Consortium. This project will select and prepare global climate model results for use in power system planning models. A related report (Evaluation of Global Climate Models for Use in Energy Analysis) assesses the performance of various global climate models from the Coupled Model Intercomparison Project Phase 6 data archive for their historical skill with respect to energy system performance and for their future projections under multiple climate change scenarios. Building from that report, we describe the selection of a climate scenario (Shared Socioeconomic Pathway [SSP] 2-4.5) and five climate models: TaiESM1, EC-Earth3-CC, GFDL-CM4, EC-Earth3-Veg, and MPI-ESM1-2-HR. We describe the model selection criteria, which were based on the quality of the match between model results under historical conditions and on the representation of the range of future values for several variables. These results will be downscaled via an open-source generative machine learning method called Super-Resolution for Renewable Energy Resource Data with Climate Change Impacts.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗
DOE OSTI2021

High-Resolution Spatio-Temporal Estimation of Net Ecosystem Exchange in Ice-Wedge Polygon Tundra Using In Situ Sensors and Remote Sensing Data

Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes; i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Trends in Spitzer Secondary Eclipses

It is well established that the magnitude of the incident stellar flux is the single most important factor in determining the day–night temperature gradients and atmospheric chemistries of short-period gas giant planets. However, it is likely that other factors, such as planet-to-planet variations in atmospheric metallicity, C/O ratio, and cloud properties, also contribute to the observed diversity of infrared spectra for this population of planets. In this study, we present new 3.6 and 4.5 μm secondary eclipse measurements for five transiting gas giant planets: HAT-P-5b, HAT-P-38b, WASP-7b, WASP-72b, and WASP-127b. We detect eclipses in at least one bandpass for all five planets and confirm circular orbits for all planets except WASP-7b, which shows evidence for a nonzero eccentricity. Building on the work of Garhart et al., we place these new planets into a broader context by comparing them with the sample of all planets with measured Spitzer secondary eclipses. We find that incident flux is the single most important factor for determining the atmospheric chemistry and circulation patterns of short-period gas giant planets. Although we might also expect surface gravity and host star metallicity to play a secondary role, we find no evidence for correlations with either of these two variables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗