Search NASA⌕ Search

SEARCH · Search NASA

Results for “intersectionality”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

COMPASS-FME Synoptic Sites Level 2 Sensor Data v2-1

This is the version 2-1 Level 2 (L2) data release for COMPASS-FME environmental sensors located at our synoptic field sites. COMPASS-FME is studying sites in two distinct regions, the Chesapeake Bay and the Western Lake Erie Basin. We established the network at seven "synoptic" (observational) sites along the Chesapeake Bay and Lake Erie coastlines, collectively generating over three million observations per month, to track and comprehend environmental changes where land and water intersect. Additionally, the two regions provide an interesting contrast of saltwater and freshwater coasts that allow us to differentiate the impacts of inundation and coastal water chemistries in two nationally important coastal systems. Level 2 (L2) data consist of sensor observations from the COMPASS-FME synoptic sites, TEMPEST, and DELUGE. Compared to the L1 data, these are more consistent (always 15-minute timestamps for the entire year); better QA/QC’d (out of bounds, out of service, and extreme outlier values are removed); and more complete, with a gap-filled time series available alongside the main observations, and additional derived (calculated) variables. L2 data are intended to be rapidly and easily usable in analyses and simulations. However, algorithmic outlier identification always carries the risk of removing valid data, and Level 1 data may be more suitable for analyses that focus on variability or extreme events. This dataset includes: - An overall dataset README file that describes the current version, gives citation and contact information, etc. - Site- and year-specific folders, each holding variable-specific Parquet (a high performance, space efficient format; see https://parquet.apache.org) data files for each site and plot in that year. - Metadata files within each site-year folder provide full information on data units, expected ranges, contact information, detailed flood times, as well as a general description of the site. - Environmental sensor types that appear in the data files include weather (ClimaVUE50, CS, RM Young, and LI instruments in the graphs below); soil conditions (TEROS12); soil redox state (Redox); groundwater variables (AquaTROLL200 and AquaTROLL600); open water sondes (Exo); tree sap velocity (Sapflow); and system voltage and state (Datalogger). Data are reported every 15 minutes. Data files are in Apache Parquet, a high performance, space efficient format for tabular data. These files can be read using R's `arrow` package (https://arrow.apache.org/docs/r/), with similar tools available in other languages. Please see v2-1 L2 Sensor Package QStart.pdf for detailed information on data package structure, temporal coverage, and versioning.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗
DOE OSTI2022

Integration of energy systems

This article in MRS Bulletin and the framework set out in the introductory article articulate a scenario of renewable electrons and electrification of end use appliances and industrial processes as a plausible paradigm to realize a carbon-free energy economy. The subsequent articles cover specific sectoral or chemical applications of those renewable electrons (e.g., for hydrogen, transportation, building use, electrochemical storage, and within the chemical industry). This article addresses the intersections among and across those sectors. We describe the importance of considering integrated systems and systems of systems as we consider pathways to a decarbonized energy economy. Further, we review and summarize key insights into the innovation challenges that reside at the particular integration interfaces among sectors, and highlight the opportunity for advances in materials and processes that will be critical to successful achievement of economy-wide, low-carbon energy systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Mechano-chemical understanding of NaSICON for aqueous redox-flow batteries

Redox-flow batteries utilizing a sodium superionic conductor (NaSICON) can be cost-effective systems for grid energy storage by combining high sodium selectivity with reliance on abundant, low-cost elements. However, improving membrane toughness while maintaining a sufficiently thin membrane for ion conduction is needed. Addressing this issue requires deeper insight into the mechanical properties of NaSICON, its interactions with aqueous chemistries, and the chemo-mechanical degradation mechanisms that arise at the intersection of these phenomena. Here, we provide a framework for understanding these problems, strategies to address them, and highlight the potential of unconventional sintering and thin-film fabrication to optimize the performance of NaSICON in practical flow cells.

Grid energy storage↗
DOE OSTI2018

AmeriFlux US-A74 ARM SGP milo field

This is the AmeriFlux version of the carbon flux data for the site US-A74 ARM SGP milo field. Site Description - This site is located near the ARM SGP Central Facility at the intersection of Oklahoma highways 11 and 74, 16.5 km East of Medford, OK

Billesbach, Dave↗
DOE OSTI2023

AmeriFlux FLUXNET-1F US-A74 ARM SGP milo field

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-A74 ARM SGP milo field. This is the FLUXNET version of the carbon flux data for the site US-A74 ARM SGP milo field produced by applying the standard ONEFlux (1F) software. Site Description - This site is located near the ARM SGP Central Facility at the intersection of Oklahoma highways 11 and 74, 16.5 km East of Medford, OK

Billesbach, Dave↗
DOE OSTI2025

AmeriFlux US-DUF Denver Urban Field Station

This is the AmeriFlux version of the carbon flux data for the site US-DUF Denver Urban Field Station. Site Description - This site was established on a 25m freestanding former radio tower at the Forest Service's Rocky Mountain Research Station headquarters in Fort Collins, CO. The tower is situated next to a parking lot and several one and two story buildings. To the north and west is the campus of Colorado State University, and to the north and east is low density residential development. To the south and east are a busy intersection and rapid transit bus line, and to the south is green space surrounding the City of Fort Collins' Spring Creek Trail.

Frank, John [US Forest Service, Rocky Mountain Res↗
DOE OSTI2017

Properties of Noise Cross-Correlation Functions Obtained from a Distributed Acoustic Sensing Array at Garner Valley, California

A field test that was conducted at Garner Valley, California, on 11 and 12 September 2013 using distributed acoustic sensing (DAS) to sense ground vibrations provided a continuous overnight record of ambient noise. Furthermore, the energy of ambient noise was concentrated between 5 and 25 Hz, which falls into the typical traffic noise frequency band. A standard procedure (Bensen et al., 2007) was adopted to calculate noise cross-correlation functions (NCFs) for 1-min intervals. The 1-min-long NCFs were stacked using the time–frequency domain phase-weighted-stacking method, which significantly improves signal quality. The obtained NCFs were asymmetrical, which was a result of the nonuniform distributed noise sources. A precursor appeared on NCFs along one segment, which was traced to a strong localized noise source or a scatterer at a nearby road intersection. NCF for the radial component of two surface accelerometers along a DAS profile gave similar results to those from DAS channels. Here, we calculated the phase velocity dispersion from DAS NCFs using the multichannel analysis of surface waves technique, and the result agrees with active-source results. We then conclude that ambient noise sources and the high spatial sampling of DAS can provide the same subsurface information as traditional active-source methods.

58 GEOSCIENCES↗
DOE OSTI2018

Seismicity Rate Surge on Faults after Shut–in: Poroelastic Response to Fluid Injection

Injection of large amounts of fluid into the subsurface alters the states of pore pressure and stress in the formation, potentially inducing earthquakes. Increase in the seismicity rate after shut–in is often observed at fluid–injection operation sites, but mechanistic study of the rate surge has not been investigated thoroughly. Considering full poroelastic coupling of pore pressure and stress, the earthquake occurrence after shut–in can be driven by two mechanisms: (1) post shut–in diffusion of pore pressure into distant faults and (2) poroelastic stressing caused by fluid injection. Interactions of these mechanisms can depend on fault geometry, hydraulic and mechanical properties of the formation, and injection operation. In this work, a 2D aerial view of the target reservoir intersected by strike–slip basement faults is used to evaluate the impact of injection–induced pressure buildup on seismicity rate surge. A series of sensitivity tests are performed by considering the variation in (1) permeability of the fault zone, (2) locations and the number of faults with respect to the injector, and (3) well operations with time–dependent injection rates. Lower permeability faults have higher seismicity rates than more permeable faults after shut–in due to delayed diffusion and poroelastic stressing. Hydraulic barriers, depending on their relative location to injection, can either stabilize or weaken a conductive fault via poroelastic stresses. In conclusion, gradual reduction of the injection rate minimizes the coulomb stress change and the least seismicity rates are predicted due to slower relaxation of coupling–induced compression as well as pore–pressure dissipation.

58 GEOSCIENCES↗
DOE OSTI2024

Carbon Transport and Storage Planning and Viability Support Tools

The EDX disCO2ver Carbon Transport and Storage Planning and Viability Support Tools are made up of the Carbon Storage Planning Inquiry Tool (CS PlanIT, Justman et al. 2024) and the Carbon Storage Technical Viability Approach Support Tool (CS TVA). Together, these tools support data access to support understanding data availability to support planning efforts for carbon transport and storage. The Carbon Storage Planning Inquiry Tool (CS PlanIT) is an online web mapping application designed to help users explore, query, and evaluate multiple data layers to support and accelerate carbon storage resource and feasibility assessments and planning efforts. CS PlanIT currently contains a range of datasets associated with geologic, technical, and infrastructure factors. The data sets can be filtered geographically for an area of interest to update statistics and charts within the dashboard. The dashboard is divided into different sections called widgets, relating to different steps in the carbon storage planning process. The resources in this submission include a link to PlanIT, as well as a data catalog and link to user documentation. The original citation for the CS PlanIT tool, which has now been integrated into the toolset here, was: - Devin Justman, Scott Pantaleone, Maneesh Sharma, Lucy Romeo, Paige Morkner, CS PlanIT (Carbon Storage Planning Inquiry Tool) , 6/28/2024, https://edx.netl.doe.gov/dataset/cs-planit-carbon-storage-planning-inquiry-tool, DOI: 10.18141/2377953 The Carbon Storage Technical Viability Approach Support (CS TVA) Tool displays spatial data availability for the many components of Geologic Carbon Storage (GCS) technical viability assessment (Creason al 2025). Identifying sites suitable for GCS requires evaluating the intersection of myriad factors, including reservoir conditions, subsurface and surface hazards, infrastructure requirements, and energy community metrics. The technical viability of a site can only be confirmed for instances where all these factors have data available, and where those data support viability. Additional Resources related to the Technical Viability Assessment Tool: - Julia Mulhern, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zachary Jackson, Jacob Shay, Gabriel Creason, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Database, 3/26/2025, https://edx.netl.doe.gov/dataset/edx4ccs-carbon-storage-technical-viability-approach-database , DOI:10.18141/1984655 - Julia Mulhern, MacKenzie Mark-Moser, Gabriel Creason, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zach Jackson, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Matrix, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-matrix , DOI: 10.18141/2539979 - Gabriel Creason, Zach Jackson, Neyda Cordero Rodriguez, Julia Mulhern, Casey White, Araceli Lara, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Data Availability Results Database, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-data-availability-results-database, DOI:10.18141/2538557

Carbon storage↗
DOE OSTI2025

CT and Geophysical Data of Clinton Sandstone Cores from Ohio

Collection of computed tomography and multi-sensor core logger data of 12 wells in Ohio that intersect the Clinton Sandstone formation. This data is described along with well information in a technical report series document: Paronish, T.; Holleran, A.; Pohl, M.; Crandall, D.; Jarvis, K.; Workman, S.; Drosche, J.; McKisic, T.; Collins, C.; Thomas, M.; McDonald, J. Computed Tomography Scanning and Geophysical Measurements of the Clinton Sandstone in Ohio; DOE.NETL-2025.4949; NETL Technical Report Series; U.S. Department of Energy, National Energy Technology Laboratory, Morgantown, WV, 2025; p 68. https://doi.org/10.2172/2589262

AS↗
DOE OSTI2017

Update on the JLEIC Electron Collider Ring Design

We present an update on the lattice design of the electron ring of the Jefferson Lab Electron-Ion Collider (JLEIC). The electron and ion collider rings feature a unique figure-8 layout providing optimal conditions for preservation of beam polarization. The rings include two arcs and two intersecting long straight sections containing a low-beta interaction region (IR) with special optics for detector polarimetry, electron beam spin rotator sections, ion beam cooling sections, and RF-cavity sections. Recent development of the electron ring lattice has been focused on minimizing the beam emittance while providing an efficient non-linear chromaticity correction and large dynamic aperture. We describe and compare three lattice designs, from which we determine the best option.

Nosochkov, Y.↗
DOE OSTI2025

ESTIMATION OF THE WAKEFIELD RESONANT FREQUENCY USING DIFFERENT SIMULATION TOOLS

Reliable and fast wakefield calculations are important for beam dynamics and THz generation in accelerators. We compare approaches and limitations with different available simulation tools and compare results. As an example, we analyse a cylindrical corrugated waveguide with narrower and wider radii of 5 and 6 mm, and aperture width and periodicity of 1 and 2 mm using ABCI by KEK and ANSYS HFSS simulation software. First, ABCI was used to determine the resonant frequency of the cylindri-cal corrugated waveguide with different lengths. The corresponding results were taken as a reference to exam-ine the simulation method in the ANSYS HFSS environ-ment. The eigenmode solver of ANSYS HFSS was used to determine different resonant frequencies. This was fol-lowed by examining a cylindrical corrugated waveguide that was implemented having a length of 17 mm for dif-ferent conductor settings. The corresponding waveguide was fed by a plane wave having a resonant frequency that satisfies the intersection of the light line and dispersion curve of the corresponding mode. The analysis showed that the maximum loss factor was achieved at 39.20 and 40.07 GHz using ANSYS HFSS simulation data for differ-ent conductors, whereas ABCI resulted in 38.50 GHz. The talk will present the results gathered by different simula-tion setups implemented in ANSYS HFSS.

Unutmaz, Muhammed [SLAC National Accelerator Labor↗
DOE OSTI2025

ICAT: The Interactive Corpus Analysis Tool

The Interactive Corpus Analysis Tool (ICAT) is a Python library for creating dashboards to explore textual datasets and build simple binary classification models to help filter through them and focus on entries of interest. This tool uses a form of interactive machine learning (IML), a paradigm of “machine teaching” (Simard et al., 2017) that sits at the intersection of the fields of human computer interaction (HCI), visual analytics, and machine learning. The intent of ICAT is to allow subject matter experts (SME) with limited to no experience in machine learning to benefit from an iterative human-in-the-loop (HITL) approach to building their own model without needing to understand the details of the underlying algorithm. This interactivity is achieved by allowing the user to create features, label data points, and visually manipulate a representation of the features to manually cluster and investigate data, while a model is trained on the fly based on these actions. ICAT is built on top of the Panel (Holoviz, 2018) library, using a combination of Vega, a custom IPyWidget using D3, and ipyvuetify, and is intended to be used inside of a Jupyter environment.

Martindale, Nathan [Oak Ridge National Laboratory ↗
DOE OSTI2022

Classification and Localization of Fracture-Hit Events in Low-Frequency Distributed Acoustic Sensing Strain Rate with Convolutional Neural Networks

Summary Distributed acoustic sensing (DAS) has been used in the oil and gas industry as an advanced technology for surveillance and diagnostics. Operators use DAS to monitor hydraulic fracturing activities, examine well stimulation efficacy, and estimate complex fracture system geometries. Particularly, low-frequency DAS can detect geomechanical events such as fracture hits because hydraulic fractures propagate and create strain rate variations in the rock. Analysis of DAS data today is mostly done post-job and subject to interpretation methods. However, the continuous and dense data stream generated live by DAS poses the opportunity for more efficient and accurate real-time data-driven analysis. The objective of this study is to develop a machine learning-based workflow that can identify and locate fracture-hit events in simulated strain rate responses correlated with low-frequency DAS data. In this paper, “fracture hit” refers to a hydraulic fracture originating from a stimulated well intersecting an offset well. We start with building a single fracture propagation model to produce strain rate patterns observed at a hypothetical monitoring well. This model is used to generate two sets of strain rate responses with one set containing fracture-hit events. The labeled synthetic data are then used to train a custom convolutional neural network (CNN) model for identifying the presence of fracture-hit events. The same model is trained again for locating the event with the output layer of the model replaced with linear units. We achieved near-perfect predictions for both event classification and localization. These promising results prove the feasibility of using CNN for real-time event detection from fiber-optic sensing data. Additionally, we use edge detection techniques to recognize fracture-hit event patterns in strain rate images. The fracture-hit location can be identified using recognized pixels in the image. The accuracy of edge detection-based location identification is also plausible, but edge detection is dependent on the assumption of pattern shape and image quality, hence it is less robust compared to CNN models. This comparison further supports the need for CNN applications in image-based real-time fiber-optic sensing event detection.

Engineering↗
DOE OSTI2022

Experimental Investigation of Low-Frequency Distributed Acoustic Strain-Rate Responses to Propagating Fractures

Summary Low-frequency distributed acoustic sensing (LF-DAS) responds to strain changes due to far-field fracture propagation. To better understand the LF-DAS response to fracture propagation, we performed laboratory-scale hydraulic fracture experiments with embedded optical strain sensors. The objectives of this research are to generate hydraulic fractures of known geometry, measure the strain response along the embedded fiber-optic cable comparable to LF-DAS measurements, and use the results to inform the interpretation of field-derived LF-DAS data. The experiments were conducted in unconfined transparent cubic blocks with a dimension of 8 in. on each side. The block was made of transparent epoxy to visualize the fracture propagation. Fiber-optic sensing cables were embedded in the block at different distances to the source of injection. We injected dyed water through an injection tubing to generate a transverse, radial fracture along an initial flaw. An optical interrogator recorded the response of offset fiber Bragg grating (FBG) strain sensors normal to the plane of the fracture. The strain data were visualized on a waterfall plot, akin to visualizations of field-derived LF-DAS data. Dimensional analysis was used to scale the laboratory results to field conditions. We compared the evolution of the strain response at the fiber-optic cable, injection pressure, and rate with known fracture geometry. The measured strains were compared with Sneddon’s (1946) linear elastic solution for a penny-shaped crack and found to follow this behavior. The generated radial fractures in transparent media can be modeled with Sneddon’s linear elastic radial fracture model and a Mode I critical stress intensity factor. The LF-DAS characteristic response of a narrowing region of extension surrounded by compression was exhibited as a fracture that approached and intersected the fiber-optic cable. The experimentally derived strain and strain-rate waterfall plots with known fracture geometry, injection rate, and pressure response provide insight in understanding LF-DAS responses in the field. Furthermore, we developed a method to estimate fracture geometry evolution from the fiber-optic strain data and validated the method against the experimental data.

Engineering↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗
DOE OSTI2025

Using multimodal X-ray computed tomography to advance 3D petrography: A non-destructive investigation of olivine inside a carbonaceous chondrite

Rocks form in three dimensions through time and studying them provides information from inside dynamic systems we cannot otherwise observe. Yet how we typically access the interior of the rocks themselves to gain that information may limit our understanding and influence how we reconstruct the processes that formed them. Here, we demonstrate combined non-destructive 3D X-ray imaging techniques that produce quantitative densitometric and crystallographic maps of entire individual grains inside a rock. Olivine grains throughout a sample of the carbonaceous chondrite Northwest Africa (NWA) 11346 were each characterized by size, shape, composition, zoning intensity, and crystallographic orientation. The addition of 3D crystallographic mapping to calibrated 3D densitometric analysis—used to calculate chemical composition—demonstrates a fully non-destructive petrographic method and provides unique insight. For instance, in our case, using crystallographic data to delineate individual grains and then measuring the 3D size, shape, and composition of each distinguishes variably reset relict grains from those later crystallized after a melting event. Intersection in a 2D slice could not have led to this interpretation because the integration of three-dimensional size, rounding, composition, location, and crystallographic orientation measured from each grain forms the key patterns. Multimodal laboratory X-ray imaging has strong potential to advance 3D petrography.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gender in Mineral Names

Minerals are the fundamental constituents of Earth, and mineral names appear in scientific literature for disciplines including geology, chemistry, materials science, biology, and medicine, among others. Choosing a name is the full responsibility of the authors of new mineral proposals submitted to the International Mineralogical Association (IMA). Scientific nomenclature and its traditions have evolved over time and, consequently, mineral names track changes in the landscape of mineralogy with respect to language, technology, and culture. To evaluate these changes, the namesake information for all 5896 minerals approved by the IMA or ‘grandfathered’ into use as of December 2022 was recorded and categorized within a workable database. The compiled information yields diverse insights into the intersection of science and culture and could also be used to project future trends. In this study, we used the name database to investigate gender diversity among mineral eponyms. More than half (c. 54%) of all mineral species are named after people, the identities of whom are largely a reflection of the people that have historically been involved, in one way or another, in the geosciences and in the mining industry. Of the 2738 people with minerals named for them, approximately 6.1% are (interpreted to be) women. Nearly all minerals named for women were named during the last sixty years, although the rate of growth in the year-on-year percentage of women among new mineral namesakes has slowed since about 1985. If current and historical trends hold, our model predicts that women will not comprise more than about 10.35% of newly established mineral namesakes in future years. The representation of women among mineral namesakes also differs starkly among countries. For example, Russians comprise 43.11% of women with minerals named for them, but account for only 15.12% of all eponyms. However, there are additional disparities beyond the proportions of namesakes. For scientists who were alive when a mineral was named for them, women were an average of 3.74 years older than men when evaluated over the same timespan (1954–2022). These results demonstrate that gender-based disparities are imprinted into current mineral nomenclature and indicate that gender parity among new mineral namesakes is impossible without unprecedented changes in the upstream demographics that are most likely to affect naming trends.

58 GEOSCIENCES↗