Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MT data”

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 55 records · Page 3

2025 TEM Workshop

The TEM Data Management Workshop will take place on August 26 from 9 a.m. to 12 p.m. MT, and will be held virtually on TEAMS. The primary goal of this workshop is to engage NSUF users and stakeholders in discussions about the data needs for the utilization of AI and ML in the analysis of TEM data. Key topics to be covered include data storage, data sharing, data tagging, metadata inclusion, standardized data formats, data augmentation, and annotated training datasets. Additionally, the workshop will provide valuable insights into resources such as the Nuclear Research Data System (NRDS) for data storage and sharing, as well as open-source codes for data analysis.

Bachhav, Mukesh↗

Kilauea Magnetotelluric Dataset

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the USGS Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilaueas plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone (ERZ), where electrical contact resistances are on the order of 100 kOhm. This data supports the continuing efforts to increase geothermal power on the island of Hawaii. Each of the 70 EDI files are the MT impedance tensors for 1 site. There is also a description of the processing of the data and a site map showing the locations of each site.

15 GEOTHERMAL ENERGY↗

Wabash CarbonSAFE Static and Dynamic Modeling: Task 9.0 (Technical Report)

The objective of the Wabash CarbonSAFE project’s static and dynamic modeling task is to assess the feasibility of storing 50 million tonnes (1.67 million metric tonnes annually; MMTA) of industrially-sourced carbon dioxide (CO 2 ) in a commercial-scale geological storage complex at Wabash Valley Resources LLC (WVR) gasification facility near Terre Haute, Indiana over a period of 30 years. The targeted formations for storing CO 2 are: 1) Mt. Simon Sandstone (MSS) and the 2) Potosi Dolomite (Knox Group). All of the available data from the recently drilled Wabash #1 stratigraphic test well (now plugged and abandoned) were used in the construction of both the static and dynamic models. Geologic models were constructed to characterize both the Mt. Simon Sandstone and Potosi Dolomite storage complexes. Dynamic simulation models were constructed and used to assess the feasibility of injecting CO 2 into the Mt. Simon and Potosi formations. The geocellular models for the Potosi Dolomite and Mt. Simon Sandstone were built using Petrel™, Schlumberger’s reservoir modeling software. The dynamic simulations were run using Landmark’s Nexus ® reservoir simulation software.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Tracing the impacts of Mount Pinatubo eruption on regional climate using spatially-varying changepoint detection

Significant events, such as volcanic eruptions, can have global and long-lasting impacts on climate. These global impacts, however, are not uniform across space and time. Understanding how the Mt. Pinatubo eruption affects global and regional climate is of great interest for predicting the impact on climate due to similar events as well as understanding the possible effect of the stratospheric aerosol injections proposed to combat climate change. While many studies illustrated the impact of the Pinatubo eruption on a global scale, studies at a fine regional scale are scarce. Here, we propose a Bayesian spatially-varying changepoint detection and estimation method to trace the impact of Mt. Pinatubo eruption on regional climate. Our approach takes into account the diffusing nature and spatial correlation of the climate changes attributed to the volcanic eruption. We illustrate our method and demonstrate its advantages over an existing changepoint detection method through simulations. Finally, we apply our method to monthly stratospheric aerosol optical depth and surface temperature data from 1985 to 1995 to detect and estimate changepoints following the 1991 Mt. Pinatubo eruption. Our results quantitatively characterize the spatial pattern of the eruption’s impact on regional climate, complementing the previous studies on the global impact of the Pinatubo eruption.

Aerosol optical depth↗

Data-Efficient Strategies for Probabilistic Voltage Envelopes under Network Contingencies

This work presents an efficient data-driven method to construct probabilistic voltage envelopes (PVE) using power flow learning in grids with network contingencies. First, a network-aware Gaussian process (GP) termed Vertex-Degree Kernel (VDK-GP), developed in prior work, is used to estimate voltage–power functions for a few network configurations. The paper introduces a novel multi-task vertex degree kernel (MT-VDK) that amalgamates the learned VDK-GPs to determine power flows for unseen networks, with a significant reduction in the computational complexity and hyperparameter requirements compared to alternate approaches. Simulations on the IEEE 30-Bus network demonstrate the retention and transfer of power flow knowledge in both N-1 and N-2 contingency scenarios. The MT-VDK-GP approach achieves over 50 % reduction in mean prediction error for novel N-1 contingency network configurations in low training data regimes (50–250 samples) over VDK-GP. Additionally, MT-VDK-GP outperforms a hyper-parameter based transfer learning approach in over 75 % of N-2 contingency network structures, even without historical N-2 outage data. Furthermore, the proposed method demonstrates the ability to achieve PVEs using sixteen times fewer power flow solutions compared to Monte-Carlo sampling-based methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗

One Earth Energy Static and Dynamic Reservoir Modeling

This report presents the static and dynamic reservoir modeling conducted for the CarbonSAFE Phase III Illinois Storage Corridor project to assess the feasibility of commercial-scale CO 2 storage in the Mt. Simon Sandstone at the One Earth Energy (OEE) site in McLean County, Illinois. Three-dimensional geocellular models of the Mt. Simon storage complex were developed in Petrel ® by integrating petrophysical log data, core analyses, and seismic surveys from the OEE #1 stratigraphic test well and two nearby wells, with multiple model versions created as new data became available. Dynamic reservoir simulations, performed using Landmark's Nexus software, progressed through three phases (preliminary, sensitivity, and UIC Class VI permit studies) evaluating injection scenarios across varying rates, well orientations, permeability models, and multi-well configurations. Results demonstrate that commercial-scale storage is feasible: three injection wells spaced approximately one mile apart can store a total of 90 million tonnes of CO 2 over 20 years, producing a combined plume with an equivalent radius of 3.2 miles and a maximum pressure-front-defined Area of Review of 178 mi 2 at the end of injection that diminishes to 34 mi 2 after 50 years of post-injection monitoring. Sensitivity analyses indicate that a 20% change in porosity or permeability yields approximately a 7% change in AoR radius, and that perforating the high-permeability arkosic zone minimizes the pressure front compared to injection in the upper Mt. Simon Sandstone.

09 BIOMASS FUELS↗

Aerosol Thermodynamics (Aerosol Liquid Water Concentration and Aerosol pH) at S2 during CoURAGE

These data report aerosol liquid water content (ALWC) and aerosol pH for the S2 (Mt. Airy) site during CoURAGE for April - June 2025. The ISORROPIA-II aerosol thermodynamic equilibrium model was used to compute ALWC and pH. The model was run in forward mode, with solids formation disabled (metastable mode), according to Pye et al. (Atmospheric Chemistry and Physics, 2020). Information about the model can be found at https://www.epfl.ch/labs/lapi/models-and-software/isorropia/, and in Fountoukis and Nenes (Atmospheric Chemistry and Physics, 2007). Inputs to the model were: (1) measured Temperature and relative humidity, both from the DoE ARM “metwxt” data products at S2, (2) aerosol chemical composition from the ARM ACSM, and (3) gas-phase ammonia concentrations. For ISORROPIA-II, ambient RH values above 0.995 were input to the model as 0.995; therefore, caution should be exercised interpreting ALWC when RH was above 0.995. Aerosol composition or NH3(g) concentrations below the LOD were input into the model as 0.5*LOD. Model runs were only conducted for data points at 30-min resolution with all three inputs (meteorology, ACSM, and NH3). If any of the three model inputs were missing, aerosol thermodynamic outputs are flagged as -999. Similarly, aerosol thermodynamic outputs are flagged as -999 if measured NH3(g) and aerosol ammonium (NH4+) were simultaneously below the respective LODs. The ALWC includes liquid water from inorganic aerosol components, calculated directly by ISORROPIA-II, as well as liquid water from organics, estimated from the ACSM organics assuming a kappa value of 0.12 according to Guo et al. (Atmospheric Chemistry and Physics, 2015).

Aerosol Liquid Water Content↗

Application of Electromagnetic Methods for Reservoir Monitoring with Emphasis on Carbon Capture, Utilization, and Storage

The Controlled-Source ElectroMagnetic (CSEM) method provides crucial information about reservoir fluids and their spatial distribution. Carbon dioxide (CO2) storage, enhanced oil recovery (EOR), geothermal exploration, and lithium exploration are ideal applications for the CSEM method. The versatility of CSEM permits its customization to specific reservoir objectives by selecting the appropriate components of a multi-component system. To effectively tailor the CSEM approach, it is essential to determine whether the primary target reservoir is resistive or conductive. This task is relatively straightforward in CO2 monitoring, where the injected fluid is resistive. However, for scenarios involving brine-saturated (water-wet) or oil-wet (carbon capture, utilization, and storage—CCUS) reservoirs, consideration must also be given to conductive reservoir components. The optimization of data acquisition before the survey involves analyzing target parameters and the sensitivity of multi-component CSEM. This optimization process typically includes on-site noise measurements and 3D anisotropic modeling. Based on our experience, subsequent surveys tend to proceed smoothly, yielding robust measurements that align with scientific objectives. Other critical aspects to be considered are using magnetotelluric (MT) measurements to define the overall background resistivities and integrating real-time quality assurance during data acquisition with 3D modeling. This integration allows the fine tuning of acquisition parameters such as acquisition time and necessary repeats. As a result, data can be examined in real-time to assess subsurface information content while the acquisition is ongoing. Consequently, high-quality data sets are usually obtained for subsequent processing and initial interpretation with minimal user intervention. The implementation of sensitivity analysis during the inversion process plays a pivotal role in ensuring that the acquired data accurately respond to the target reservoirs’ expected depth range. To elucidate these concepts, we present an illustrative example from a CO2 storage site in North Dakota, USA, wherein the long-offset transient electromagnetic method (LOTEM), a variation of the CSEM method, and the MT method were utilized. This example showcases how surface measurements attain appropriately upscaled log-scale sensitivity. Furthermore, the sensitivity of the CSEM and MT methods was examined in other case histories, where the target reservoirs exhibited conductive properties, such as those encountered in enhanced oil recovery (EOR), geothermal, and lithium exploration applications. The same equipment specifications were utilized for CSEM and MT surveys across all case studies.

Barajas-Olalde, César↗

Computed Tomography Scanning and Geophysical Measurements of the One Earth Energy Well #1 Core

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize the Lower Mt. Simon Sandstone and Eau Claire Formation core from the One Earth Energy Well #1 (OEE Well #1) in the Illinois Basin. The primary impetus of this work is a collaboration between the U.S. Department of Energy (DOE), the Illinois State Geological Survey (ISGS), the University of Illinois, and One Earth Energy to characterize and make publicly available core information of interest to carbon sequestration efforts in the Illinois Basin. This stratigraphic well and the core data produced in this report will aid in understanding the carbon sequestration potential of the Lower Mt. Simon Sandstone and the sealing capacity of the overlying units. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/illinois-storage-corridor-one-earth-energy-1-core.

47 OTHER INSTRUMENTATION↗

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗

Multi-Scale 3D Imaging for Machine Learning Property Upscaling: Mt. Simon Sandstone Case Study

Petrographic properties of principal target reservoirs for carbon sequestration, such as the Mt. Simon Sandstone, are relevant to broad interest groups. The Mt. Simon Sandstone is a deep, saline, regionally extensive Cambrian sandstone, overlain by low permeability sealing formations, making it one of the viable geologic carbon storage reservoirs in the Midwestern US. Its thickness (exceeding 2400 ft in some localities), depth, and lateral extent, combined with high porosity and permeability make it a high-priority target of multiple ongoing geologic carbon sequestration efforts in the United States of America. The National Energy Technology Laboratory in Morgantown, West Virginia, has been engaged in characterization efforts of the Mt. Simon for over a decade, with a strong focus on Computed Tomographic data acquisition. Data generated during this period has been hitherto not accessible to the public. This archival effort focused on preservation of historical CT data and associated metadata, and facilitating their accessibility, culminating with the publication of the entire dataset on NETL’s Energy Data eXchange (EDX) and the associated Gill et. al (2024) paper.

Gill, Magdalena K.↗

MicroFisher: Fungal taxonomic classification for metatranscriptomic and metagenomic data using multiple short hypervariable markers

AbstractProfiling the taxonomic and functional composition of microbes using metagenomic (MG) and metatranscriptomic (MT) sequencing is advancing our understanding of microbial functions. However, the sensitivity and accuracy of microbial classification using genome– or core protein-based approaches, especially the classification of eukaryotic organisms, is limited by the availability of genomes and the resolution of sequence databases. To address this, we propose the MicroFisher, a novel approach that applies multiple hypervariable marker genes to profile fungal communities from MGs and MTs. This approach utilizes the hypervariable regions of ITS and large subunit (LSU) rRNA genes for fungal identification with high sensitivity and resolution. Simultaneously, we propose a computational pipeline (MicroFisher) to optimize and integrate the results from classifications using multiple hypervariable markers. To test the performance of our method, we applied MicroFisher to the synthetic community profiling and found high performance in fungal prediction and abundance estimation. In addition, we also used MGs from forest soil and MTs of root eukaryotic microbes to test our method and the results showed that MicroFisher provided more accurate profiling of environmental microbiomes compared to other classification tools. Overall, MicroFisher serves as a novel pipeline for classification of fungal communities from MGs and MTs.

Wang, Haihua↗

Design and demonstration of a low‐field magnetic resonance imaging rhizotron for in‐field imaging of energy sorghum roots

Abstract Root phenotyping provides critical information to plant breeders for developing varieties with improved drought tolerance, greater root biomass, and greater nutrient use efficiency. Phenotyping roots in the natural environment is important for understanding the effect of the soil environment on root genotypic expressions. The goal of this work was to design and test a field‐scale mobile low‐field magnetic resonance imaging (LF‐MRI) Rhizotron that produces actionable root phenotyping data. We demonstrated this novel technology for root visualization and quantification using a LF‐MRI Rhizotron operating at 47 mT with two soil types. The LF‐MRI Rhizotron weights 453 kg, with a height of 90 cm, a diameter of 28 cm and an imaging field of view of 28 cm × 28 cm. The unit was operated in a Belk clay (Entic Hapluderts) and Weswood silt loam (Udifluventic Halustepts) generating 2‐D and 3‐D image data sets. The 2‐D image data had a collection time of 16.5 min per image at an image resolution of 2.2 mm per pixel. The 3‐D data had a collection time of 13 h per image with a 2.2 × 2.2 × 2.2 mm voxel resolution. Low‐field magnetic resonance imaging worked well for visualizing roots in moderate to high clay soils, demonstrating the potential for this technology; however, the broad application of this platform is hampered due to the prohibitively long scanning time to obtain 3‐D images. By increasing the field strength, and therefore the signal‐to‐noise ratio, faster scan times can enable a more useful system for root phenotyping.

Bagnall, G. Cody↗

Elastic functional changepoint detection of climate impacts from localized sources

Detecting changepoints in functional data has become an important problem as interest in monitoring of climate phenomenon has increased, where the data is functional in nature. Here, the observed data often contains both amplitude (y-axis) and phase (x-axis) variability. If not accounted for properly, true changepoints may be undetected, and the estimated underlying mean change functions will be incorrect. In this article, an elastic functional changepoint method is developed which properly accounts for these types of variability. The method can detect amplitude and phase changepoints which current methods in the literature do not, as they focus solely on the amplitude changepoint. This method can easily be implemented using the functions directly or can be computed via functional principal component analysis to ease the computational burden. We apply the method and its nonelastic competitors to both simulated data and observed data to show its efficiency in handling data with phase variation with both amplitude and phase changepoints. We use the method to evaluate potential changes in stratospheric temperature due to the eruption of Mt. Pinatubo in the Philippines in June 1991. Using an epidemic changepoint model, we find evidence of a increase in stratospheric temperature during a period that contains the immediate aftermath of Mt. Pinatubo, with most detected changepoints occurring in the tropics as expected.

54 ENVIRONMENTAL SCIENCES↗

Diurnal Trends and Meteorological Factors Influencing the Variability of Fluorescent Bioaerosol in Mt. Crested Butte, Colorado During SAIL

Abstract Despite the crucial role of bioaerosol particles (BAP) in our climate system, local ecosystems, and human health, our grasp on their atmospheric interactions is hampered by a lack of high‐resolution and long‐term data, which is essential for understanding their abundance and variability in response to meteorological conditions. To discern these relationships in a high‐altitude mountainous terrain is even less well‐studied. Therefore, we deployed a Wideband Integrated Bioaerosol Sensor (WIBS‐Neo) for three months during the first biologically active season of the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in summer 2022 on Mt. Crested Butte in Colorado (elevation 3,140 m ASL). Here, we report real‐time fluorescent BAP (FBAP) data collected from 15 June to 13 September 2022, sampled from within the ARM Aerosol Observing System (AOS). To our knowledge, these are the first high altitude (>3,000 m ASL) continuous measurements of FBAP made in North America. During this deployment, we observed on average 21% and as many as 48% (hourly maximum) of particles within the detection size range ( to ) of the WIBS as FBAP. Our analysis presents the diurnal cycles for seven distinct types of FBAP, showing unique patterns and, in some cases, correlations with temperature and solar radiation cycles. Abundance and composition of FBAP varied with relative humidity and precipitation. Precipitation events appeared to both cause emission and removal of FBAP, whereas dust and smoke events had no significant effect, highlighting the critical role of meteorology on FBAP at high altitudes.

54 ENVIRONMENTAL SCIENCES↗

Wabash CarbonSAFE. Final Report

This document summarizes work detailed in separate Wabash CarbonSAFE reports; the report describes the data collection efforts of the project and consolidates the geologic characterization, well testing, and storage complex modeling results for the Mt. Simon Sandstone and Potosi Dolomite, two distinct reservoirs characterized at the Wabash CarbonSAFE project site. Also presented are summaries of work to characterize the CO 2 source and infrastructure network, as well as summaries of reports analyzing stakeholder engagement options, policy, regulatory and legal considerations, and risk assessment associated with the Wabash CarbonSAFE project site. The report then presents recommendations for the next steps for site characterization, identifies data gaps for future activities, and provides an overall assessment of site potential.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗

Charge Transfer and Spin Dynamics in a Zinc Porphyrin Donor Covalently Linked to One or Two Naphthalenediimide Acceptors

We report quantum coherence effects on charge transfer and spin dynamics in a system having two degenerate electron acceptors are studied using a zinc 5,10,15-tri(n-pentyI)-20-phenylporphyrin (ZnP) electron donor covalently linked to either one or two naphthalene-1,8:4,5-bis(dicarboximide) (NDI) electron acceptors using an anthracene (An) spacer, ZnP-An-NDI (1) and ZnP-An-NDI 2 (2), respectively. Following photoexcitation of 1 and 2 in toluene at 295 K, femtosecond transient absorption spectroscopy shows that the electron transfer (ET) rate constant for 2 is about three times larger than that of 1, which can be accounted for by the statistical nature of incoherent ET as well as the electron couplings for the charge separation reactions. In contrast, the rate constant for charge recombination (CR) of 1 is about 25% faster than that of 2. Using femtosecond transient infrared spectroscopy and theoretical analysis, we find that the electron on NDI 2 ∙- in 2 localizes onto one of the two NDIs prior to CR, thus precluding electronically coherent CR from NDI 2 ∙- . Conversely, CR in both 1 and 2 is spin coherent as indicated by the observation of a resonance in the 3* ZnP yield following CR as a function of applied magnetic field, giving spin-spin exchange interaction energies of 2J = 210 and 236 mT, respectively, where the line width of the resonance for 2 is greater than 1. These data show that while CR is a spin-coherent process, incoherent hopping of the electron between the two NDIs in 2, consistent with the lack of delocalization noted above, results in greater spin decoherence in 2 relative to 1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗