Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “underground”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A novel shallow bore ground heat exchanger for ground source heat pump applications—Model development and validation

A novel ground heat exchanger design, the Underground Thermal Battery (UTB), was proposed as a low-cost alternative to the conventional vertical bore ground heat exchanger (VBGHE) for the application of ground source heat pumps (GSHPs). The UTB is designed to be installed in a shallow borehole (less than 6 m deep) and thus could cost much less than the conventional VBGHE, which usually is installed in vertical bores 60 m deep. By utilizing natural convection of water and phase change materials, the UTB tempers its temperature change in response to thermal loads, which helps improve the efficiency of GSHPs. A one-dimensional (1D) model of the UTB has been developed and validated with the measured performance data of a small-scale UTB, as well as the simulation results of a more detailed three dimensional (3D) numerical model, which accounts for both heat transfer and fluid dynamics in the UTB. Further, this 1D model is computationally much more efficient than the 3D model and thus can simulate the annual performance of the UTB with acceptable computation time. The 1D model has been used to evaluate the performance of the UTB so that it can be compared with the conventional VBGHE.

42 ENGINEERING↗
DOE OSTI2023

Development of ultralow-background cryogenic calorimeters for the measurement of surface α contamination

Next-generation experiments searching for rare events must satisfy increasingly stringent requirements on the bulk and surface radioactive contamination of their active and structural materials. The measurement of surface contamination is particularly challenging, as no existing technology is capable of separately measuring parts of the 232 Th and 238 U decay chains that are commonly found to be out of secular equilibrium. We will present the results obtained with a detector prototype consisting of 8 silicon wafers of 150mm diameter instrumented as bolometers and operated in a low-background dilution refrigerator at the Gran Sasso Underground Laboratory of INFN, Italy. Furthermore, the prototype was characterized by a baseline energy resolution of few keV and a background <100nBq/cm 2 in the full range of α energies, obtained with simple procedures for cleaning of all employed materials and no specific measures to prevent recontamination. Such performance, together with the modularity of the detector design, demonstrate the possibility to realize an alpha detector capable of separately measuring all alpha emitters of the 232 Th and 238 U chains, possibly reaching a sensitivity of few nBq/cm 2 .

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

Increases in Cell-Wall Homogalacturonan but Decreases in Xylogalacturonan Accompany Transition from Dormant to Vegetative Stages in Chrysolaena Obovata Rhizophores

The rhizophores of Chrysolaena obovata, a perennial herb native to the Brazilian Cerrado, are underground storage organs that accumulate fructans and support vegetative propagation. Here, we investigated how the content of fructans and the composition of rhizophore cell wall polysaccharides change across developmental stages and contribute to seasonal adaptation. In addition to the accumulation of fructans in vacuoles (~26-55% of the dry mass according to the developmental stage), cell-wall polysaccharide composition changed throughout the growth cycle, reflecting structural adjustments linked to developmental and environmental cues. These cell walls have a type I architecture, typical of eudicots. Although total pectin content remained stable, homogalacturonans increased in proportion during the transition from dormant to vegetative stages, whereas amounts of xylogalacturonans decreased. Xylogalacturonans accumulate within the middle lamella during the dormant phase and are potentially associated with flexibility and water retention. Hemicelluloses consisted mainly of low-substituted arabinoxylans and xyloglucan, consistent with a secondary cell wall profile. Xylans decreased slightly in the transition from the dormant to vegetative stage, while xyloglucan proportions increased. Overall, our findings show that seasonal metabolic changes in C. obovata involve dynamic modifications of rhizophore cell wall polysaccharides, which contribute to organ development and drought adaptation throughout the seasonal cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Well–based monitoring of CO 2 geological sequestration operations in saline aquifers: Critical insights into key questions

Geological carbon sequestration in saline aquifers is one of the most promising strategies to help mitigate emissions of CO 2 to the atmosphere. Significant challenges in ensuring the security of the sequestration process rest in the evolution and expansion of the CO 2 plume in the subsurface. The ability to track the movement of the injected CO 2 poses another challenge. Critical questions related to the integrity of the sequestration operations in saline aquifers relate to plume characteristics that can we monitor using well-based variables. We addressed this and related questions using an integrated modeling framework through a numerical investigation of carbon sequestration in saline aquifers during long-term and post-injection periods. This modeling paradigm incorporates the effect of structural, geological, and petrophysical characteristics. That way, we can account for critical physicochemical processes, rock-fluid interactions, and lithology dependencies. The well fluid variables investigated include fluid composition, pH, fluid density, and ion activity. We learned that fluid property analytics and diagnostics can be powerful tools to estimate the movement of CO 2 and its storage in different trapping mechanisms. These analytics can help optimize operational aspects and simplify reservoir-scale models while still reflecting the complex nature of the CO 2 interactions underground and offering insights into plume evolution.

58 GEOSCIENCES↗
DOE OSTI2020

Mineral-specific heterogeneous neptunium sorption onto geological repository rocks in oxic and anoxic conditions and various temperatures

Neptunium poses challenges for radioactive wastes management due to its long half-life and ingrowth from americium decay, as well as its high mobility in water in its oxidized forms. Np(V)-bearing aqueous solutions were contacted with polished and crushed gneiss and dolerite from the Yeniseisky site, Russia, which is under consideration as a location for an underground repository, both in oxic and anoxic conditions at 20 and 90 °C. Neptunium sorption on crushed gneiss and dolerite in oxic conditions at 20 °C was quite low, but increased substantially in anoxic conditions. Autoradiography on samples treated at 20 °C revealed that sorbed Np was associated mostly with dark colored iron-bearing minerals, and solvent extraction indicated partial reduction of Np(V) to Np(IV). Sorption of Np from aqueous solutions contacted with dolerite and gneiss at 90 °C was in all cases higher than at the lower temperature. Autoradiography of polished rock samples treated in anoxic conditions at 90 °C revealed that Fe-bearing minerals had 3–4 orders of magnitude more Np sorbed on them than the quartz and feldspar grains. X-ray photoelectron spectroscopy indicated this sorbed Np was present as NpO 2 . This study supports the concept that the geologic conditions present at the planned repository site will retard Np migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Mobilization of trace metals from caprock and formation rocks at the Illinois Basin – Decatur Project demonstration site under geological carbon dioxide sequestration conditions

One concern for geologic CO 2 sequestration is the potential leakage of CO 2 or CO 2 -saturated brine containing trace metals into overlying aquifers, which poses the risk of adversely affecting underground sources of drinking water. In this work, rock and brine samples collected in the Illinois Basin – Decatur Project (IBDP), a large-scale geological CO 2 sequestration demonstration project in Decatur, Illinois, USA, as well as synthetic brine samples were used to understand the mobilization of trace metals as a result of CO 2 –rock–brine interactions under IBDP-specific conditions (50 °C and 20.6 MPa). Rock sample characterization indicated that, at the IBDP site, trace metal concentrations were greatest in cap rock samples from the Eau Claire Formation, compared with those from formations above and below this Formation. The natural brine samples collected from the Mt. Simon Sandstone reservoir were highly reducing and saline, with trace metal concentrations up to 680 times greater than the U.S. Environmental Protection Agency-prescribed drinking water standards. Batch leaching experiments indicated that both trace metal mobilization from rock and immobilization from the brine occurred when high-pressure CO 2 was introduced into the rock-brine system. The amounts of metals mobilized from the rock generally accounted for <5% of the total metals in the rock, but for some metals, including Ni, Pb, and Tl, up to 63% of the metals in the rock were mobilized to the brine. Leaching of trace metals into synthetic brines was different from that into the natural brine. The results of this study provided current information on the trace element sources and will help in risk simulations and experimental design to further evaluate potential impact of leakage on groundwater quality, which is important for future GCS projects to consider for monitoring and containment assurance purposes.

58 GEOSCIENCES↗
DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗
DOE OSTI2025

From clutter to clarity: Emergent neural operators via questionnaire metrics

Real-world datasets in chemical engineering and bioengineering processes—such as those from catalytic reactors, multiphase flows, polymerization reactors, bioreactors, and clinical trials—can often be unlabeled or disorganized, rendering the training of existing supervised learning models ineffective at learning the underlying dynamics. To salvage these datasets for decision-making, we first seek to obtain clarity from the cluttered data. Here, we present a framework for developing “structural” generative models, discovering emergent equations, and constructing efficient emulators from scrambled datasets by integrating unsupervised organizational learning techniques (Questionnaires) with advanced deep learning architectures (Deep Hidden Physics Models and Deep Operator Networks). Our approach is demonstrated on two illustrative model systems: (a) a 1D advection–diffusion partial differential equation representing a winding underground pipe and (b) an ensemble of Stuart–Landau oscillators, an agent-based system of coupled ordinary differential equations. In both cases, we successfully reconstruct meaningful spatial, temporal, and parameter embeddings from scrambled data, enabling good predictions of system dynamics. As a result, we highlight the framework’s potential for broader applications, enabling data-driven system identification in fields with inherently disorganized or hidden parameter spaces.

42 ENGINEERING↗
DOE OSTI2020

Stimuli-responsive petroleum cement composite with giant expansion and enhanced mechanical properties

Shrinkage, as an inherent property of cement during the curing process, has been a concern in the oil & gas industry for centuries. Therefore, expansive cement has been studied for many decades, primarily through incorporation of expansive additives. Nevertheless, traditional expansive additives are not able to achieve the required expansive ratio or maintain a good mechanical property for the cement composite at extreme underground condition. Accordingly, a new generation of expansive additive, which not only has the required expansion, but also work at the downhole high temperature and high pressure environment, and maintain or even enhance the mechanical property, is highly desired. In this work, a high enthalpy storage shape memory polymer (SMP) as a new class of expansive additive is investigated. It is found that the cement composite achieves 1.4% circumferential expansion by only 6% by weight of SMP additives. The compressive strength and flexural strength are also enhanced at the same time, which is hardly achievable by other expansive additives with the same concertation. Good pumpability is also proved by rheological study. The mechanism controlling the enhanced properties is reveled through morphological study and element mapping analysis at the SMP/cement interface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling↗
DOE OSTI2020

Reactive chemical transport simulations of geologic carbon sequestration: Methods and applications

Chemical reaction simulations are considerably used to quantitatively assess the long-term geologic carbon sequestration (GCS), such as CO 2 sequestration capacity estimations, leakage pathway analyses, enhanced oil recovery (EOR) efficiency studies, and risk assessments of sealing formations (caprocks), wellbores, and overlying underground water resources. All these require a deep understanding of the CO 2 - associated chemical reactions. To ensure long-term, safe CO 2 sequestration in the intended formations, modeling is the only way to plausibly assess the CO 2 flow, reaction, and transport over thousands of years. This review summarizes the multiple methodologies for describing homogeneous and heterogeneous chemical reaction patterns and multiscale application examples, the recent progress and current status of chemical reaction simulations for GCS, and the impact of such simulations on geological CO2 sequestration performance. Technical gaps and future challenges are also discussed for further study. The trends and challenges of such studies include: (1) the combination of coupled chemical, mechanical, and transport processes with calibrated experiments and associated uncertainty/risk assessments; (2) enhancement of the ability to simulate detailed geophysical and geochemical equations to mimic in situ conditions; and (3) characterization of multiscale subsurface systems with detailed conceptual models and assignment of suitable boundary conditions for field-scale sequestration fields. One major gap remaining is the current lack of accurate (and scale-justified) kinetic and equilibrium chemical reaction parameters under reservoir conditions. Advanced models that couple chemical, mechanical, and transport processes with scale-justified parameters, from lab to field-scale experiments, are required for quantitative assessments of sequestration capacity and the long-term safety of GCS projects.

58 GEOSCIENCES↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2021

Distribution of rare earth elements and other critical elements in beneficiated Pennsylvania anthracites

The Pennsylvania Anthracite Fields are in a complex tectonic and metamorphic terrain, historically hosting one of the largest concentrations of coal mining in the USA. Anthracite mining now largely consists of the surface mining of the pillars remaining from the prime years of underground mining. The geochemical study of the sized coal products and the refuse (largely rock) from three preparation plants (breakers) demonstrates that Principal components analysis (PCA) of select major oxide, minor element, and rare earth elements illustrates some differentiation among the products from the individual plants. The rock samples, with abundant quartz and metamorphic Al-Si minerals and with a lower ash-basis REE concentration than the coals, were distinctly separated from the coal samples on the PCA plots. Plots of Gd N /Gd N * vs. Eu N /Eu N * and Ce N /Ce N * vs. Eu N / Eu N * showed differentiation between the individual suites of coals showed that the refuse samples had distinct REE distributions compared to the associated coals. Further, several minor and trace elements show enrichments in the coal samples. Lithium, with concentrations of up to 314 ppm on an ash basis, is among the most promising of the critical elements, exceeding the enrichment of the REY and Sc.

58 GEOSCIENCES↗
DOE OSTI2022

The 2017 Pohang, South Korea, Mw 5.4 main shock was either natural or triggered, but not induced

Understanding the cause of the November 2017 Pohang main shock is of considerable importance to the future of the geothermal industry because of its large magnitude compared to prior expectations based on case histories of other projects involving underground fluid injection. Of the three possibilities – induced, triggered or natural, “induced” can be ruled out based on the disproportionately large seismic moment of the main shock. Furthermore, whether natural or triggered, the source of seismic hazard at Pohang was tectonic strain accumulation, not fluid injection. Arguably, the most timely indicator of seismic hazard and risk in the environs of Pohang was the September 2016 Mw 5.4 Gyeongju earthquake, which was natural and located about 40 km south of Pohang along the same active fault system.

58 GEOSCIENCES↗
DOE OSTI2022

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗
DOE OSTI2020

Multi-physics modeling of injected nanoparticles effect on remediation of CO 2 leakage through cracks

One potential risk in CO 2 sequestration is the leakage of carbon dioxide, which can result in contamination of underground water, creating potential threats to existing ecosystems. The common leakage pathway is through the pre-existing fractures or discontinuities within cement in the wellbore, incurred by the environmental conditions imposed on the cement. Injecting nanoparticles into pre-existing cracks is one of the most recently proposed ideas for mitigating fracture propagation in cement CO 2 sequestration. To demonstrate the feasibility of this new technology, a numerical approach was taken in this work, as it is challenging to investigate it in a laboratory setting. We proposed a coupled ALE (Arbitrary Lagrangian–Eulerian)–DEM (Discrete Element Method)–peridynamic modeling strategy within the LS-Dyna package to investigate the intertwined interaction among the CO 2 fluid flow, native fluid (e.g., brine), particle clusters, and cracks within the cement. The numerical results demonstrate that injected nanoparticles can effectively reduce the pressure exerted on the crack surface. Accordingly, the potential fracture propagation at the crack tip would be reduced compared to corresponding cases without nanoparticels as pressurized by fluid flow. This result verifies the effectiveness of proposed nanoparticle injection technology. Finally, using this established modeling strategy, the effect of filling particles on the fracture mitigation for different crack geometries (e.g. particle cluster patterns, aspect ratio of crack aperture and length) and CO 2 reservoir pressure are examined. The result shows that injected particles successfully reduce the fracture propagation in these scenarios.

42 ENGINEERING↗
DOE OSTI2020

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗
DOE OSTI2020

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗