Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “favorability”

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 1,243 records · Page 69

Optimization under uncertainty of a hybrid waste tire and natural gas feedstock flexible polygeneration system using a decomposition algorithm

Market uncertainties motivate the development of flexible polygeneration systems that are able to adjust operating conditions to favor production of the most profitable product portfolio. However, this operational flexibility comes at the cost of higher capital expenditure. A scenario-based two-stage stochastic nonconvex Mixed-Integer Nonlinear Programming (MINLP) approach lends itself naturally to optimizing these trade-offs. This work studies the optimal design and operation under uncertainty of a hybrid feedstock flexible polygeneration system producing electricity, methanol, dimethyl ether, olefins or liquefied (synthetic) natural gas. A recently developed C++ based software framework (named GOSSIP) is used for modeling the optimization problem as well as its efficient solution using the Nonconvex Generalized Benders Decomposition (NGBD) algorithm. Two different cases are studied: The first uses estimates of the means and variances of the uncertain parameters from historical data, whereas the second assesses the impact of increased uncertain parameter volatility. The value of implementing flexible designs characterized by the value of the stochastic solution (VSS) is in the range of 260–405 M$ for a scale of approximately 893 MW of thermal input. Increased price volatility around the same mean results in higher expected net present value and VSS as operational flexibility allows for asymmetric exploitation of price peaks.

42 ENGINEERING↗

Thermally induced cracking on the massive concrete structure of the NSLS II synchrotron and its engineering remediation

Synchrotron accelerator facilities such as the NSLS II require extreme stability, both transient (short-term) and quasi-static (long-term) to achieve the desired resolution performance. Consequently, even μm-level movements, particularly differential movements between locations in the concrete structure supporting the accelerator electron beam lattice (storage ring) or high sensitivity experiments (experimental floor) will lead to serious degradation of its performance. Differential settlement in the overall structure or structural movement exceeding anticipated levels will inevitably degrade the performance and will require intervention. Presented in this paper are the design philosophy of the NSLS II ring structure favoring a monolithic ring, the observed cracking behavior of the young NSLS II concrete following casting and in combination with extreme ambient temperature fluctuation, the results of a non-linear, high-fidelity numerical analysis used to emulate the observed cracking and establish the driving mechanism, the numerical analysis-based identification of the crack-arresting solution and finally the implementation of the remediation solution and the long-term performance of the adopted engineering solution. The multi-stage process revealed that computational methods such as non-linear finite element methods have the potential of providing engineering guidance even when complex structures and in combination with non-linear materials, such as steel reinforcement and concrete are involved.

42 ENGINEERING↗

A review of composite polymer-ceramic electrolytes for lithium batteries

All solid-state lithium batteries are garnering attention in both academia and industry. Lithium-ion conductive polymers and lithium-ion conductive ceramics are the two major classes of solid electrolytes that have prevalently been pursued for many years. However, each of them has its own advantages and disadvantages. One approach to overcome the disadvantages and get the best out of each of those materials is a solid composite electrolyte that combines the advantages of inorganic ceramic electrolytes and solid polymer electrolytes. Such composite electrolytes can offer acceptable ionic conductivity, high mechanical strength, and favorable interfacial contact with electrodes, which can greatly improve the electrochemical performance of all-solid-state batteries compared to cells based on a polymer electrolyte alone or a ceramic electrolyte alone. We present in this review the state-of-the-art composite polymer-ceramic electrolytes in view of their electrochemical and physical properties for the applications in lithium batteries. The review mainly encompasses the polymer matrices, various ceramic filler materials, and the polymer/ceramics composite systems. Specifically, the structures, ionic conductivities, electrochemical/chemical stabilities, and fabrications of solid composite electrolytes are discussed in-depth. On the basis of previous work, a perspective on future research directions is highlighted for developing high-performance composite polymer-ceramic electrolytes.

25 ENERGY STORAGE↗

Anomalously high elastic modulus of a poly(ethylene oxide)-based composite electrolyte

The practical use of lithium metal anodes in solid-state batteries requires a polymer membrane with high lithium-ion conductivity, thermal/electrochemical stability, and mechanical strength. The primary challenge is to effectively decouple the ionic conductivity and mechanical strength of the polymer electrolytes. We report a remarkably facile single step synthetic strategy based on in-situ crosslinking of poly(ethylene oxide) (xPEO) in the presence of a woven glass fiber (GF). Such a simple method yields composite polymer electrolytes (CPE) of anomalously high elastic modulus up to 2.5 GPa over a broad temperature range (20 °C – 245 °C) that has never been previously documented. An unsupervised machine learning algorithm, K-mean clustering analysis, was implemented on the hyperspectral Raman mapping at the xPEO/GF interface. Using such a unique means, we show for the first time that the promoted mechanical strength originates from xPEO and GF interactions through dynamic hydrogen and ionic bonding. High ionic conductivity is achieved by the addition plasticizer (e.g. tetraglyme), where trifluoromethanesulfonate anions are tethered to the xPEO matrix and Li + cations are favorably transported through coordination with the plasticizer. Further, stringent galvanostatic cycling tests indicates the CPE can be stably cycled for >3000 h in a Li-metal symmetric cell at a moderate temperature (nearly 1500 Coulombs/cm 2 Li equivalents), outperforming most of the PEO-based electrolytes. The GF reinforced CPE reported here has multifunctional uses, such as solid electrolytes for all solid-state batteries and membranes for redox-flow batteries. Although the focus of this study is on lithium-based batteries, the results are equally promising for other alkali metal based batteries such as sodium and potassium.

25 ENERGY STORAGE↗

Unlocking the self-supported thermal runaway of high-energy lithium-ion batteries

Layered Ni-rich LiNi x Mn y Co 1-x-y O 2 (NMC) materials are the most promising cathode materials for Li-ion batteries due to their favorable energy densities. However, the low thermal stability typically caused by detrimental oxygen release leads to significant safety concerns. Determining the pathways of oxygen evolution reaction is essential, as the ideal safety countermeasure is to break the reaction chain of thermal runaway. In this study, we demonstrate that two endogenous pathways of oxygen involved in strong exothermic reactions lead the NMC811|graphite pouch cell to an uncontrollable state, and we quantify the individual contribution of the pathways to thermal runaway. Approximately 41% of thermal-induced oxygen reacts aggressively with ethylene carbonate (EC) at the cathode/electrolyte interface with 16% heat generation, accelerating the self-heating rate and thereby further triggering thermal runaway. The residual oxygen that survives the reaction with carbonate spreads to the lithiated anode with major heat generation (65%), bringing the battery to the maximum destructive temperature during thermal runaway. By confirming the significant roles of EC and anode, a deeper understanding on battery fire was achieved. Finally, the revealed mechanism can help guide studies on stopping the two reaction pathways, allowing for the safer use of high-energy lithium-ion batteries in the future.

25 ENERGY STORAGE↗

Spatial and temporal characterization of municipal solid waste based on resource recovery pathways

This study presents a two-year, quarterly assessment of MSW across four source sectors (residential, schools, restaurants, and grocery stores) from sixteen sites across five U.S. states. MSW was manually sorted into 27 categories and aggregated into pathway fractions: high-moisture (HM) organics, low-moisture (LM) organics, recyclable (RC) materials, and residuals for disposal. Organics represented 89 % of the MSW stream. The largest fraction was HM organics consisting of food waste (31 %) and yard waste (3 %), with large coefficient of variations (CV), 79 and 278 %, respectively, reflecting high seasonal and site variability that varied significantly (p < 0.01) across sampling periods. The HM fraction showed properties favorable for anaerobic digestion, with moisture content ranging from 56 to 95 % and volatile solids ranges of 86-95 %. In contrast, the LM and RC fractions remained more stable (plastics CV = 41 %; paper CV = 53 %) with heating values up to 26.9 MJ/kg across sources, reflecting suitability for gasification. Microstructural analysis revealed less porosity in residential waste sampled at the landfill, which can influence preprocessing efficiency and microbial accessibility. Pathway informed allocations showed that 35 % of MSW is suitable for anaerobic digestion, 36 % for gasification, and 18 % for recycling, leaving 11 % requiring landfill disposal. These results provide quantitative evidence to determine feedstock allocation, waste-to-energy system design, and the development of data-driven sustainability and resource recovery strategies within a circular bioeconomy.

09 BIOMASS FUELS↗

Manganese: The overlooked contaminant in the world largest mine tailings dam collapse

Manganese (Mn) is an abundant element in terrestrial and coastal ecosystems and an essential micronutrient in the metabolic processes of plants and animals. Mn is generally not considered a potentially toxic element due to its low content in both soil and water. However, in coastal ecosystems, the Mn dynamic (commonly associated with the Fe cycle) is mostly controlled by redox processes. Here, we assessed the potential contamination of the Rio Doce estuary (SE Brazil) by Mn after the world’s largest mine tailings dam collapse, potentially resulting in chronic exposure to local wildlife and humans. Estuarine soils, water, and fish were collected and analyzed seven days after the arrival of the tailings in 2015 and again two years after the dam collapse in 2017. Using a suite of solid-phase analyses including X-ray absorption spectroscopy and sequential extractions, our results indicated that a large quantity of Mn II arrived in the estuary in 2015 bound to Fe oxyhydroxides. Over time, dissolved Mn and Fe were released from soils when Fe III oxyhydroxides underwent reductive dissolution. Due to seasonal redox oscillations, both Fe and Mn were then re-oxidized to Fe III , Mn III , and Mn IV and re-precipitated as poorly crystalline Fe oxyhydroxides and poorly crystalline Mn oxides. In 2017, redox conditions (Eh: -47 ± 83 mV; pH: 6.7 ± 0.5) favorable to both Fe and Mn reduction led to an increase (~880%) of dissolved Mn (average for 2015: 66 ± 130 μg L -1 ; 2017: 582 ± 626 μg L -1 ) in water and a decrease (~75%, 2015: 547 ± 498 mg kg -1 ; 2017: 135 ± 80 mg kg -1 ) in the total Mn content in soils. The crystalline Fe oxyhydroxides content significantly decreased while the fraction of poorly ordered Fe oxides increased in the soils limiting the role of Fe in Mn retention. The high concentration of dissolved Mn found within the estuary two years after the arrival of mine tailings indicates a possible chronic contamination scenario, which is supported by the high levels of Mn in two species of fish living in the estuary. Our work suggests a high risk to estuarine biota and human health due to the rapid Fe and Mn biogeochemical dynamic within the impacted estuary.

54 ENVIRONMENTAL SCIENCES↗

Unique impacts of strong and westward-extended western Pacific subtropical high on ozone pollution over eastern China

As a subtropical anticyclonic high-pressure system that typically forms over the northwestern Pacific Ocean in summer, the Western Pacific subtropical high (WPSH) affects meteorological conditions and ozone pollution in China. The relationship between maximum daily 8-h average ozone (MDA8 O 3 ) concentrations and the extremely strong and westward-extended WPSH occurred in 2022 is investigated using observations, reanalysis data and GEOS-Chem model simulations. The intensity of WPSH has a significant positive correlation with MDA8 O3 over southern China during July-August in 2022, with a correlation coefficient of +0.44, but the correlation is negative (–0.40) in northern China. During the strong WPSH days, MDA8 O 3 increased by 16.5µgm -3 (16.4% relative to July-August average) over southern China and decreased by 19.0µgm -3 (14.5%) in northern China compared to the weak WPSH days. The unique dipole pattern in the relationship between ozone levels and the WPSH in 2022 exhibited a contrast to that during 2015–2021. The difference is primarily due to the extremely strong WPSH intensity and its unusual westward expansion in 2022. In this case, an anomalous anticyclone at 500 hPa dominates over southern China, which creates conditions conducive for ozone formation and accumulation. The anticyclone weakened horizontal winds and reduced the dispersion of ozone, alongside a high temperature and low relative humidity, which favored the chemical production of ozone. In contrast, abnormal northerly winds enhanced ozone diffusion in northern China and the low temperature reduced ozone chemical production. Here, this study reveals the mechanism for the significant impact of strong and westward-extended WPSH on ozone concentrations over China, emphasizing the role of the WPSH location in modulating meteorology and ozone levels.

54 ENVIRONMENTAL SCIENCES↗

Scaling open-weight large language models for hydropower regulatory information extraction: A systematic analysis

Information extraction from regulatory and technical documents using large language models (LLMs) involves practical trade-offs between extraction quality and computational cost. We evaluate eight open-weight LLMs spanning 0.6B–70B parameters on hydropower licensing documents and report deployment-oriented evidence under a unified extraction schema and evaluation protocol. Across the model set, we observe clear scale-dependent trends in both baseline extraction quality and the effectiveness of reflective reasoning (self-checking) under our fixed-prompt, no-augmentation setting. Mid-scale models often provide a favorable balance of accuracy and efficiency, whereas the smallest models show limited or inconsistent gains from the reasoning variants tested. Larger models achieve the highest overall F1 scores but incur substantially greater compute and infrastructure requirements. We further find that reliability failure modes can distort conventional metrics in this domain: in particular, high recall can coincide with systematic extraction errors when models fabricate values for fields that are absent from the source text, underscoring the importance of conservative null handling and evidence-grounded evaluation. Overall, our study provides a reproducible resource–performance comparison for open-weight LLM-based extraction in hydropower regulatory documentation and offers practical guidance for model selection under different deployment constraints.

Evaluation protocol↗

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗

A tri-level distribution locational marginal price-based demand response framework

Here, in this paper, we propose a tri-level, nested, two-stage price-based demand response (PBDR) framework that considers distribution locational marginal price (DLMP) as DR enabler between load-serving entities (LSE), demand response providers (DRPs), and customers in the day-ahead distribution market. It enables LSE and customer interactions by using multiple DRPs, positioned in-between, and independently optimizes their objectives. The problem is formulated using linear power flow with approximated power losses and its application in DLMP as DR pricing. The tri-level problem is solved using a nested reformulation & decomposition (R&D) method and tested on the real Indian-108 bus distribution system under various dynamic pricings. Further, the temporal–spatial variations in DLMPs are assessed using fairness criteria. Numerical analyses demonstrate that DLMP applications can effectively improve economic efficiency, and transparency in DR programs valuation with a favorable fairness margin. The results show that DLMP as DR pricing signal induces (0-2) % variation in DLMP for DR participation up to 10 %. Further, it gives over 90 % fairness over temporal–spatial variation for all the customers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Dispatch optimization of electric thermal energy storage within System Advisor Model

A stand-alone electric thermal energy storage (ETES) system converts low-value electricity into heat using resistance heating elements. During periods of high-value electricity, an ETES system uses a thermodynamic power cycle to convert stored thermal energy back to electricity. These dispatchable systems derive value from their ability to store energy when prices are low and generate electricity when prices are favorable, i.e., energy arbitrage. Consequently, dispatch optimization of system operations, through maximizing revenue subject to system constraints, is essential to evaluate the economic value of a particular system design. While stand-alone ETES systems offer potential advantages as dispatchable grid storage technologies, there is a lack of a neutral third-party, publicly available, open-source model to evaluate the performance, dispatch, and financial viability of these systems. To address this problem, we have developed a techno-economic model for stand-alone ETES systems, within National Renewable Energy Laboratory's (NREL's) System Advisor Model (SAM). We implement a mixed-integer linear program to determine an ETES optimal operating schedule that maximizes electricity sales less maintenance costs caused by operation and cycling given temporal-varying grid electricity prices. Our contributions include a mixed-integer linear program for energy arbitrage of an ETES system, an ETES performance model through a publicly-available software (i.e., SAM), and an exercise of our model through case studies that compare ETES operational strategies and annual financial metrics. With our dispatch optimization model, we were able to improve revenue by 20% compared to a myopic heuristic while reducing the operational cost of the ETES system through decreases in cycle starts, cycles per day, and heater starts.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Energy storage planning for enhanced resilience of power systems against wildfires and heatwaves

Extreme weather events pose significant risks to power grid stability due to their severe consequences and potential for widespread failures. Energy storage systems hold great potential for enhancing grid resilience against such events by providing reliable power during peak demand periods. However, accurately quantifying the size, location, and investment costs of new energy storage assets is a complex task, as energy storage planning decisions depend on the investment choices of other generation technologies and the integration of new transmission projects. Here, this paper presents a novel capacity expansion planning framework that simultaneously optimizes investments in energy storage, generation, and transmission, determining their optimal size, location, and type, while incorporating extreme weather events into long-term planning. More specifically, our stress-event-informed planning framework integrates the impact of heatwaves and wildfires into the planning process, identifying least-cost investment solutions that comply with policy goals and enhance grid resilience. The proposed framework employs machine-learning-based modeling to project heatwave-induced loads and performance-based risk assessment to evaluate wildfire-driven transmission line derates. Using industry-standard datasets to accurately represent the transmission topology of the Western Interconnection (WI) system, the proposed framework is applied to the WI 40-zone system, with investment decisions reported for the years 2030, 2035, and 2040. Simulation results reveal that with just a 10% increase in investment costs, resilience against extreme events can be significantly improved, with investment decisions heavily favoring energy storage, particularly 4-hour energy storage systems.

25 ENERGY STORAGE↗

Results of a fishery-independent longline survey targeting coastal sharks in the eastern Bahamas between 1979 and 2013

Long-term trends in shark abundance offer important insights for fisheries management. Few fisheries-independent, extended time-series data exist for coastal shark species in the northwest Atlantic Ocean or Caribbean outside of the United States. A decades old dataset comprised of standardized longline surveys provided us with the opportunity to characterize diversity and evaluate changes in the relative abundance of sharks in the eastern Bahamas between 1979–1990 and 2011–2013. Caribbean reef sharks Carcharhinus perezi and tiger sharks Galeocerdo cuvier were the two most commonly captured species throughout the survey and were the only two species analyzed in detail. For both species, sex ratios favored females. Using Bayesian generalized linear models, we found that Caribbean reef shark catch-per-unit-effort (CPUE) exhibited a positive relationship with year and temperature, while size exhibited no relationship with year and a negative relationship with temperature. There was no relationship between tiger shark CPUE and year but there was a positive relationship between size and year. These results are presented in the context of Bahamian shark conservation efforts and the results of other surveys in the region.

Fisheries-independent survey, Caribbean reef shark↗

Combustion, criteria pollutant and soot property assessment of mixing-controlled high-load engine operation with gasoline and diesel

Gasoline compression ignition (GCI) is an emerging technology that can achieve high fuel efficiency and low criteria emissions in diesel engine platforms using gasoline and gasoline-like fuels. While low- and mid-range loads are enabled with low temperature combustion (LTC) approaches, high-loads are more efficient with high temperature combustion (HTC) requiring mixing-controlled combustion (MCC). Lower soot emissions are achieved with gasoline even in high-load MCC than with diesel, but detailed soot properties from gasoline and diesel have been rarely discussed. In order to better understand how high-load operations affect soot properties, in this study gasoline and diesel were tested in a six-cylinder heavy duty engine at 14 bar brake mean effective pressure (BMEP) at similar combustion phasing. Effects of pilot mass ratio and EGR on heat release rates, cylinder temperatures, ignition delay (ID) and combustion duration (CD) were examined, relating with emissions variations with applications of the two fuels. Soot mass and particle number emissions were slightly lower with gasoline than with diesel, resulting in a more favorable NOx-soot trade-off with gasoline like other reports. Contrary to other findings in the literature, gasoline and diesel soot samples presented similar particle distribution shapes, carbon crystalline structures and primary particle sizes, regardless of fuel type and EGR rate. Consequently, the current result proposes that soot particles have undergone similar soot formation and oxidation processes despite reactivity difference and oxygen availability under the examined high-load MCC with similar combustion phases.

02 PETROLEUM↗

Distribution of rare earth elements in the pilot-scale processing of fly ashes derived from eastern Kentucky coals: Comparisons of the feed and processed ashes

Beneficiated fly ash from the combustion of Central Appalachian high volatile bituminous coals was extracted with HNO 3 in a pilot-scale processing plant. Several major oxides (notably CaO and SO 3 , but also including Fe 2 O 3 , MgO, K 2 O, and P 2 O 5 ) and minor elements (Mn, As, Sr, Ba, and Pb) are depleted in the post-HNO 3 -extraction spent ash. The total lanthanides, Y, and Sc concentration is reduced by about 20% in the spent ash, with Gd showing the greatest decrease. Along with Gd, Nd and Dy are also well differentiated between the feed and spent ashes, with La and Sm showing minimal partitioning. The Gd decrease is correlated with the depletion of Fe 2 O 3 . The heavy rare earth elements (REE heavier than Eu) and Y are disproportionately concentrated in the HNO 3 -leachate compared to the light REE. For the ashes studied, Sc did not partition between the feed and spent ashes. Pozzolanicity tests show that the compressive strength and strength activity indices of the spent ash + ordinary Portland cement (OPC) mixes are comparable to 100% OPC, indicating that the spent ashes produced in the pilot-scale runs have the potential to be sold as a Class F fly ash. Finally, the beneficiated ash chemistry influences the chemistry of the post-HNO 3 -extraction spent ash and the HNO 3 -leachate. A 500-ppm-REE fly ash will presumably be a more economically favorable feedstock than an ash with a significantly lesser concentration.

01 COAL, LIGNITE, AND PEAT↗