Search NASA⌕ Search

SEARCH · Search NASA

Results for “geothermal systems”

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 307 records · Page 17

Final Report for Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity

Characterizing in situ stress is essential in Enhanced Geothermal System (EGS) for both risk assessment and operation design. Standard stress measurement techniques such as diagnostic fracture injection tests (DFITs) require creating hydraulic fractures and apply analytical methods such as the G -function analysis to estimate the stress. These classical approaches rely on several assumptions, including but not limited to: (i) the induced fracture is planar and persistent, and (ii) the fracture plane is normal to the minimum horizontal stress. However, these assumptions may not hold in EGS environments like Utah FORGE, where thermo-hydro-mechanical (THM) coupling, wellbore deviation, and rock heterogeneity can alter the local stress condition and consequently lead to fracture complexities in both near-well and far-field regions.

15 GEOTHERMAL ENERGY↗

Complementary Subsurface Characterization Methods to Develop a Geologic Model for the EGS Collab Experiment, Sanford Underground Research Facility

The EGS (Enhanced Geothermal Systems) Collab project was performed within the Sanford Underground Research Facility (SURF) with a goal of understanding processes and evaluation of models related to hydraulic stimulation of rock at depth. The present work deals with the development of Testbed 2 where experiments were conducted at a depth of 1.25 km and were located within a well-characterized testbed in a metamorphic, amphibolite host rock. A total of eleven boreholes varying in length between 10.6 m and 81.2 m were continuously cored to develop the testbed. In addition to the continuous coring of the amphibolite host rock, geophysical instrumentation supporting electrical resistivity tomography (ERT), microearthquake (MEQ) detection, and optical fiber providing distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) were installed in the monitoring boreholes, all of which produced a comprehensive complementary suite of characterization and monitoring technologies. Pre-stimulation characterization of the groundwater conditions identified only two hydraulically significant fractures, neither of which transects the central portion of the testbed. Flow and pressure monitoring indicated low preexisting pore pressure conditions likely affected by the mine openings.

15 GEOTHERMAL ENERGY↗

Utah FORGE Project 3-2417: Simulations for Distributed Acoustic Sensing Strain Signatures as an Indicator of Fracture Connectivity

This dataset encompasses simulations of strain signatures from both hydraulically connected and "near-miss" fractures in enhanced geothermal systems (EGS). The files and results are presented from the perspective of digital acoustic sensing's (DAS) potential to differentiate the two fracture types. This dataset was acquired by the FOGMORE R&D project (Fiber Optic Geophysical MOnitoring of Reservoir Evolution), under Utah FORGE R&D Project 3-2417. Included are simulation and results via MatLab and COMSOL files, as well as a thesis and paper summarizing the results. Some stimulated fractures may be incomplete, approaching but not intersecting the production well. These "near-miss" fractures can be addressed in future stimulation stages or re-stimulated to complete the connection. We propose the use of fiber optic distributed acoustic sensing (DAS) as a method by which near-miss stimulated fractures may be identified and distinguished from hydraulically connected fractures. The low-frequency sub-nanostrain signatures of both complete and near-miss fractures in DAS data are simulated in this study using a hydrogeomechanical discrete fracture network model. The spatial distribution of strain was found to be an accurate indicator. However, this indicator must be evaluated in the context of DAS gauge length and spatial sampling. These simulations are a precursor to tests conducted at FORGE in 2023.

15 GEOTHERMAL ENERGY↗

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗

Project Phase 1 Report: Reducing Data Center Peak Cooling Demand and Energy Costs With Cold Underground Thermal Energy Storage (Cold UTES)

Cold Underground Thermal Energy Storage (Cold UTES) is an ultra-long duration grid energy storage technology. With Cold-UTES, low-cost grid power is converted to cold thermal energy and stored in the native subsurface rock at the point of use. Cold UTES is one approach within the general category of engineered geothermal systems. Cold UTES for peak-hour cooling of data centers (DCs) was studied for deployment in Maricopa County, Arizona and Loudoun County, Virgina using thermal storage capacities from 4 GWh-th to over 1,000 GWh-th (>1 Terawatt-hour). The two sites have different power and transmission systems, available grid energy resources, daily and seasonal load profiles, weather conditions, and grid regulatory requirements. The study results indicate high value for both locations and because of this, likely indicates value across most of the US and the world. The basis of the study was a 1,000 MW-e hourly electric use of DC computing and auxiliary loads, which was modeled as 1 GW-th of thermal load to a dry cooled heat rejection system - i.e. a cooling system that does not consume water. The electric power required for cooling the DC varies as the air temperature changes. In cold weather only the dry-coolers are used, with an electrical load for cooling load as low as 10 MW-e. In hot summer hours, chillers and dry- coolers are required, which raises the electrical load for cooling load to as much as 300 MW-e. The continuous and peak cooling electrical loads result in a grid interconnection requirement of no less than 1,300 MW-e. From both a grid and thermal design modeling perspective the 1.3 GW-e could either be a single facility or result from the total load at multiple sites.

15 GEOTHERMAL ENERGY↗

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗

Seismic monitoring of EGS fracture stimulations at Utah FORGE (Part 1): Time-lapse variations of b-values and Shear-wave Splitting Rates of induced microearthquakes

Characterizing geologic formations near geothermal wells is crucial for efficient extraction of subsurface geothermal energy resources. To monitor effective fracture stimulation of enhanced geothermal systems (EGS), we analyze a catalog of induced microearthquakes to characterize seismicity induced by Phase 2C fracture stimulations of the Utah FORGE project and derive the b-values for 424 induced microearthquakes in three stages of fracture stimulation. Here, we find the time-lapse variations of these b-values. We apply eigenvector rotation and waveform correlation to 3C geophone data to separate fast S 1 and slow S 2 waves, determine their directions, measure their lag time, and calculate the S-wave splitting rate (SSR) for six induced microearthquakes in Stage 1 and 21 in Stage 3. Both SSR values and b-values increase with time. The increase in these values may be related to heavily fractured zones created by fracture stimulations during Stage 2 injections. Large and highly variable SSR values in Stage 3 may reflect heterogeneities in the fracture zones. Time-lapse variations of SSR values and b-values may be one of the most useful approaches to monitoring and characterizing the effectiveness and efficiency of various types of fracture stimulations.

58 GEOSCIENCES↗

Ecology of methyl-coenzyme M reductase encoding Thermoproteota

The recent demonstration that members of at least three classes of archaea affiliated with the Thermoproteota superphylum are involved in the production of the climate-active gas methane has sparked discussions about how well we understand the diversity of methanogens. Here, in this work, we show that members of all three of these lineages, as well as several other, yet uncultured and physiologically uncharacterized groups within the Thermoproteota that encode the key enzyme of anaerobic methane cycling, methyl-coenzyme M reductase (MCR), are widely distributed in anoxic ecosystems. We postulate that the taxonomic, metabolic, and ecological diversity of methanogenic and MCR-encoding Thermoproteota are poorly understood, and that the contribution of methylotrophic and thermoproteotal methanogenesis to methane production is largely unknown. We hypothesize that thermoproteotal methanogens could contribute, potentially substantially, to methane emissions in many anoxic environments that harbor methylated precursors, including wetlands, sediments, peat, rice paddies, wastewater sludge, and geothermal systems. We highlight the necessity to experimentally test the (eco)physiology of these widely distributed archaea using both culture-dependent ( in vitro ) and culture-independent ( in situ ) approaches to assess their potential contribution to methane emissions. Last, we stress the importance of remaining agnostic about the physiology of MCR-encoding Thermoproteota in the absence of experimental data because most of these archaea also carry the genetic potential to grow non-methanogenically.

MCR↗

Modeling aqueous association constants and mineral solubilities at subcritical and supercritical temperatures

The need for sustainable power generation has increased interest in the use of hydrothermal fluids for industrial applications. New high-enthalpy geothermal systems and biowaste-to-fuel processes are two relevant examples that employ supercritical fluids which require an in-depth understanding of complex chemical reactions occurring near the supercritical temperature of water (374 °C). As these processes operate in thermodynamic regimes that are not currently covered by a standard molar Gibbs energy of formation model, only empirical fits for single reaction systems are available which limit the use of multi-component phase equilibria calculations that are standard practice for less extreme environments. Here, in this paper, we advance a standard molar Gibbs energy of formation model able to operate in these otherwise inaccessible thermodynamic states to include species needed for key mineral solubility systems and ion association reactions. This work extends a model based on molecular statistical thermodynamics (MST) into four new systems (Na 3 PO 4 -H 2 O, LiOH-H 2 O, KOH-H 2 O, and BaSO 4 -H 2 O) by extending the model to cover 10 new species. For each of these systems, model predictions were consistently within the experimental uncertainties for the new systems covered. A breakdown of MST contributions to the model revealed that electrostatic and hard sphere contributions were key to reproducing density dependencies of standard molar Gibbs energy of formation values around the critical point of water.

42 ENGINEERING↗

Nitrogen fixation in Pampinifervens , a new species-rich genus of Aquificaceae that inhabits a wide pH range in terrestrial hot springs

Members of the Aquificota play important roles in hydrogen, sulfur, and carbon cycles in geothermal systems as ubiquitous and abundant chemolithoautotrophs; however, their roles in the nitrogen cycle are poorly defined. Here, we show that Aquificaceae isolate T-2 T from Tengchong, China, fixes nitrogen based on growth with dinitrogen as the sole nitrogen source and incorporation of 15 N atoms when grown with 15 N 2 . We further show evidence that suggests its nitrogen fixation (nif) genes are transcribed in situ in a nearby spring. We propose the name Pampinifervens diazotrophicum gen. et sp. nov., with strain T-2 T (= JCM 35475 T = DSM 116324 T ) as the nomenclatural type. Its close relative designated Pampinifervens florentissimum sp. nov. T-8 T (= CGMC 1.5214 T = JCM 33569 T ) does not fix nitrogen but respires nitrate. Analysis of 104 metagenome-assembled genomes from 84 metagenomes allowed us to circumscribe the genus Pampinifervens and identify at least 16 sympatric Pampinifervens species that dominate the Aquificaceae populations in diverse hot springs in southwest China, with distinct species globally. Analysis of these Pampinifervens genomes revealed variable presence of nif and denitrification genes and frequent gene loss and horizontal transfer. Genes for pH homeostasis were also present in Pampinifervens genomes, with a proliferation in species from Tengchong, consistent with the wide pH growth range of both Pampinifervens isolates and the prevalence of Pampinifervens in both alkaline and acidic springs. The dynamic evolution of nif and denitrification genes, along with pH homeostasis genes, may facilitate the diversification of Pampinifervens into diverse springs with different nitrogen availability.

Aquificae↗

Microbial Community Composition in Deep‐Subsurface Reservoir Fluids Reveals Natural Interwell Connectivity

Abstract The identification of natural fractures and the wells they connect is crucial for the development of geological reservoirs because it may have an important impact on reservoir model construction and hydraulic fracture propagation. In this study we investigated the use of a novel data source, the microbial community composition in the reservoir formation fluids, for identification of interwell connectivity caused by natural fractures. We verified this concept at a newly developed mesoscale enhanced geothermal system (EGS) field testbed located 4,850 ft (1,478 m) beneath ground surface at the Sanford Underground Research Facility in Lead, SD. Fluids produced at or near the EGS testbed were sampled and subjected to high‐throughput 16S rRNA gene amplicon sequencing to analyze the microbial community profile therein. Despite the typically substantial heterogeneity across the community profiles of samples spatially distributed (10 m to 1.9 km apart) throughout the site, samples from two wells at the EGS testbed showed highly similar microbial community composition, suggesting the two wells intersected the same natural fracture. This evidence of natural connectivity between the two wells at the EGS testbed was later corroborated by core log analysis and sewer camera surveys into the boreholes. Besides the field case described in this study, microbial community composition as a reservoir diagnostic tool would be applicable in a much broader context such as unconventional hydrocarbon exploitation, groundwater reservoir characterization, and environmental remediation, adding valuable “hard” data capable of pinpointing the origins of fluids unambiguously.

Zhang, Yuran↗

Creation of a Mixed-Mode Fracture Network at Mesoscale Through Hydraulic Fracturing and Shear Stimulation

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. In this study we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5-km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100-kHz, continuous seismic monitoring system deployed in six 60-m length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear-reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

58 GEOSCIENCES↗

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗

Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah

Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.

Geochemistry & Geophysics↗

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

Causal mechanism of injection-induced earthquakes through the M w 5.5 Pohang earthquake case study

Causal mechanisms for fluid injection-induced earthquakes remain a challenge to identify. Past studies largely established spatiotemporal correlations. Here, we propose a multi-process causal mechanism for injection-induced earthquakes through a case study of the 2017 M w 5.5 induced earthquake near Pohang Enhanced Geothermal System, Korea, where detailed hydraulic stimulation and on-site seismicity monitoring data provide an unprecedented opportunity. Pore pressure modeling reveals that pore pressure changes initiate seismicity on critically stressed faults and Coulomb static stress transfer modeling reveals that earthquake interactions promote continued seismicity, leading to larger events. On the basis of these results, we propose the following causal mechanism for induced seismicity: pore pressure increase and earthquake interactions lead to fault weakening and ultimately triggering larger earthquakes later in the process. We suggest that it is prudent that pore pressure change, initial seismicity locations, and Coulomb static stress transfer from seismicity earlier in the sequence are assessed in real-time.

58 GEOSCIENCES↗

Branched-chain amino acid specialization drove diversification within Calditenuaceae ( Caldarchaeia ) and enables their cultivation

Many thermophiles that are abundant in high-temperature geothermal systems have never been cultivated and are poorly understood, including deeply branching members of the archaeal phylum Thermoproteota. Here, we describe the genome-guided cultivation of one such organism, Calditenuis ramacidaminiphagus, and show that it has evolved a heterotrophic metabolism focused on branched-chain amino acids (BCAAs). Initially, fluorescence in situ hybridization and nanoscale secondary ion mass spectrometry (FISH-nanoSIMS) showed that Cal. ramacidaminiphagus assimilated amino acids rapidly in casamino acid-amended enrichment cultures. Metagenome and metaproteome analyses showed a high abundance and expression of BCAA transporter genes, suggesting a BCAA-focused metabolism. This inference was supported by the subsequent enrichment of Cal. ramacidaminiphagus in BCAA-fed cultures, reaching 2.66×10 6 cells/mL and 48.7% of the community, whereas it was outcompeted when polar amino acids were included. Metabolic reconstruction and metaproteomics suggest that BCAAs are channeled into the mevalonate pathway for lipid biosynthesis and fuel ATP production through the TCA cycle coupled with aerobic respiration and through production of branched-chain organic acids by overflow metabolism. Ancestral state reconstructions and phylogenetic analyses of 62 Caldarchaeales genomes revealed multiple horizontal transfers of BCAA transporters to the ancestor of the genus Calditenuis. Our study highlights the crucial role of BCAAs in the early evolution and niche of this genus, and suggests a high degree of resource partitioning even within low-diversity thermophilic communities.

Biological and medical sciences↗