Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Project seismic”

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 217 records · Page 12

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects

Researchers at New Mexico Tech, Los Alamos National Laboratory (LANL), Sandia National Laboratory (SNL), and other collaborators are developing a methodology to improve characterization of stress in the subsurface by means of a model-based inversion of reservoir engineering data, time lapse seismic measurements, and microseismicity. Historically, various geophysical techniques have been used in efforts to understand the state of stress in the subsurface through direct and indirect imaging of stress sensitive features (faults and fractures) and observations of transient stress related observations (time variant elastic moduli and microseismicity). Standard direct and indirect techniques for seismic fault and fracture imaging suffer from detectability limits due to reliance on high data quality and multiplicity. Variations in effective elastic properties from inversion of high quality seismic (Vertical Seismic Profile (VSP)) data may be used to infer stress near the wellbore through integration with independent experimental characterization of the stress-velocity relationship. Given a sufficiently robust observation network, microseismic emissions may be inverted to characterize focal mechanisms which, together with supporting assumptions and constraints, inform estimates of in-situ stress. While each of these techniques contributes in part to the characterization of stress within limited spatial and temporal domains, no one method provides an unambiguous stress measurement or a predictive capability over a site scale spatial extent. Challenges associated with solution non-uniqueness, measurement ambiguity, and irregular sampling may be greatly minimized through combination of one or more of independent measurements within a common framework in which realistic geological, hydrodynamic, geomechanical, and seismological constituent process models may act as constraints.

02 PETROLEUM↗

DOE Hanford Tank Side Cesium Removal System (TSCR) Update: Fabrication and Factory Acceptance Testing - 20300

The Tank Side Cesium Removal (TSCR) System is a modular, quickly-deployable, and passively safe system for efficiently removing solids and cesium from alkaline liquid waste raffinate created by fuel processing on the Hanford site. This paper provides an overview of the project with particular emphasis on completion of fabrication and factory acceptance testing (FAT). All TSCR System components are inside three enclosures known as the Process Enclosure, Ancillary Enclosure, and Control Enclosure. The Process Enclosure is seismic design category II (SDC-II) and built to withstand the rigors of transportation and lifting, high winds, volcanic ash fall, and heavy precipitation (snow) over a temperature range of -25 to 115 deg. F. The Process Enclosure and Control Enclosure were fabricated in AVANTech's Columbia, SC facility, while the Ancillary Enclosure was fabricated in AVANTech's Richland, WA facility. The design and preliminary documented safety analysis were completed in December 2019, with all work being performed in accordance with ASME NQA-1 quality assurance requirements. The Process Enclosure, Ion Exchange Columns (IXCs), and several other TSCR components have a nuclear safety function and are designated as safety significant structures, systems, and components (SS-SSCs). The NQA-1 Commercial Grade Dedication (CGD) process was used to identify critical characteristics in order to accept materials and items used for fabrication of several safety significant TSCR components. This paper also addresses unique aspects and challenges of implementing CGD in DOE facilities and systems. After fabrication and preliminary testing, the Process and Control Enclosures were shipped to AVANTech's Richland facility where the complete TSCR System was setup and connected in a represented layout and configuration that will be used on the Hanford site. Test tanks and pumps were put in place and connected to the TSCR Process Enclosure in a manner that mimicked interfacing Double Shell Tanks in the 241-AP Tank Farm. The TSCR System was then tested/operated under conditions equivalent to those planned for onsite operation. Testing was separated into five separate phases as defined below: - Phase 1: Confirmation of TSCR System readiness for the FAT; - Phase 2: Verification of TSCR System operation under standard conditions; - Phase 3: Demonstration of TSCR System operation under non-standard conditions; - Phase 4: Demonstration of ion exchange media removal from an IXC; and - Phase 5: Completion of TSCR System lift plans and site mobilization. The Hanford Tank Operations Contractor, Washington River Protection Solutions (WRPS), participated in all phases of testing. Testing at AVANTech's Richland Facility allowed WRPS personnel to familiarize themselves with TSCR operations in a radiologically clean facility prior to onsite mobilization. Engineering and operations personnel gained valuable experience that allowed them to better refine onsite operational protocols and procedures - including handling and movement of IXCs with a forklift. This paper describes the technologies used by the TSCR System as well as how these technologies are deployed and operated to advance the Hanford mission of safely, efficiently, and effectively treating tank waste. It also provides detailed information on fabrication and FAT, including operational data from simulated tank waste runs. A project schedule for WRPS' planned construction acceptance testing, operational acceptance testing, and actual tank waste treatment is provided. This presentation should be of great interest to parties responsible for the design, testing, and mobilization of modular NQA-1 technologies for the treatment of raffinate and similar liquid wastes containing high concentrations of radioactive cesium. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Design and Analyses of the NSTX-U PF1A Poloidal Field Coil Support

New poloidal field (PF) coils and supports have been designed for the National Spherical Torus Experiment Upgrade (NSTX-U) Recovery project. The PF1A coil structural requirements demand a large preload to pre-compress the coil insulation. Coil support design accounts for all operational conditions and loads of NSTX-U including bakeout, seismic events, and electromagnetic loads. Background magnetic fields produce coil Lorentz forces that react on the structural supports. These Lorentz forces include interactions from the coil current and all plasma events: normal operations, disruptions, and vertical displacement events. This paper will detail the design and design loads and highlight some of the analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncertainty Quantification of Geophysical Inversion Using Stochastic Partial Differential Equations (LDRD #218329)

This report summarizes work completed under the Laboratory Directed Research and Development (LDRD) project "Uncertainty Quantification of Geophysical Inversion Using Stochastic Differential Equations." Geophysical inversions often require computationally expensive algorithms to find even one solution, let alone propagating uncertainties through to the solution domain. The primary purpose of this project was to find more computationally efficient means to approximate solution uncertainty in geophysical inversions. We found multiple computationally efficient methods of propagating Earth model uncertainty into uncertainties in solutions of full waveform seismic moment tensor inversions. However, the optimum method of approximating the uncertainty in these seismic source solutions was to use the Karhunen-Love theorem with data misfit residuals. This method was orders of magnitude more computationally efficient than traditional Monte Carlo methods and yielded estimates of uncertainty that closely approximated those of Monte Carlo. We will summarize the various methods we evaluated for estimating uncertainty in seismic source inversions as well as work toward this goal in the realm of 3-D seismic tomographic inversion uncertainty.

58 GEOSCIENCES↗

Grid Resilience Plans: State Requirements, Utility Practices, and Utility Plan Template

As of June 2024, 14 states and one city require jurisdictional electric utilities to file resilience plans. Drawing on these requirements and filed plans, this report offers a standard template that states and utilities can consider to improve utility filings for grid resilience plans, either as part of a distribution system plan or as a separate filing. Key elements include a vulnerability assessment, description of proposed resilience programs, and projected costs and rate impacts. While many of these requirements and plans focus on extreme weather hazards, the template also can be used to address additional threats, including cyber and physical attacks and seismic events.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SCALE 6.3.1 Radiation Source Terms and Shielding Analysis for a Postulated Sodium-Cooled Fast Reactor Accident Scenario

In support of the US Nuclear Regulatory Commission non–light-water reactor fuel cycle demonstration project, SCALE 6.3.1 capabilities for radiation source term and shielding calculations are demonstrated for scenarios in the sodium-cooled fast reactor (SFR) fuel cycle. A postulated accident scenario, which consists of a seismic event causing the refueling machine to fall and release a spent fuel assembly inside the containment building (CB), is analyzed in this paper. Radiation source terms were generated for a U/TRU-10Zr metal fuel assembly with a 16.5% initial transuranic waste content and a discharge burnup of approximately 95 GWd/tHM; source terms were also generated for a high-assay low-enriched uranium metal fuel assembly (U-10Zr) with a 16.5% initial enrichment and a discharge burnup of 149.74 GWd/tHM. These radiation source terms were then used to determine the dose rate inside the CB and near the outer surface of the CB for a range of spent fuel assembly cooling times. The dose rate produced by the analyzed SFR assemblies is similar to that produced by a typical pressurized water reactor assembly with a discharge burnup of 50 GWd/MTU. Ultimately, the validation of the source terms predicted for SFRs with SCALE will need to be demonstrated via the use of assay measurements.

Radulescu, Georgeta↗

Integrating HPC simulations and physical experiments to characterize the effects of gamma radiation on seismic protective devices

Seismic protective systems, composed of seismic isolators and dampers, can substantially reduce the effects of earthquake shaking on nuclear power plants and components therein. To enable the use of these devices to protect equipment inside a plant and close to a source of radiation, the U.S. Department of Energy (DOE) funded a project at the Idaho National Laboratory (INL) and the University at Buffalo to characterize the effects of absorbed gamma dose on their mechanical properties. An early task in the project was to determine the exposure time required in the INL Foss Therapy Services (FTS) 60 Co gamma irradiator to achieve a target absorbed dose in the materials used to construct isolators and dampers, including fluids, polymers, composites, and metals. This task required the novel integration of high-performance computing (HPC), Monte Carlo N-Particle (MCNP) simulations, and irradiation experiments using Fricke dosimetry. An MCNP model of the FTS irradiator at INL was developed and validated using Fricke dosimetry. Simulations of three experiments in the irradiator, two with Fricke vials only and one with Fricke vials and a large-size isolator, predicted the Fricke-measured absorbed dose rate to within 15% in all three cases, providing high confidence in the calculation of the gamma dose absorbed in the materials comprising the seismic protective devices. The simulations demonstrated that the effects of photon scattering on absorbed dose rate in the FTS irradiator are negligible for test articles installed close to the cobalt sources and near the rear of the irradiator. The validated MCNP model of the FTS irradiator is being used to support ongoing DOE-funded experiments on seismic protective devices and could be applied to future, non-seismic-related experiments. In conclusion, the novel validation process successfully deployed for the FTS irradiator at INL could be applied to other irradiators, requiring new MCNP models and simulations, and irradiation experiments using dosimeters.

42 - ENGINEERING↗

Utilizing 3D Mechanical Earth Models for Calibration and Validation in a Large-Scale Carbon Capture and Storage Project in North Dakota

Conference paper presented at 17th International Conference on Greenhouse Gas Control Technologies (GHGT-17), Calgary, Alberta, Canada, October 20–24, 2024. Three-dimensional (3D) mechanical earth models (MEMs) are pivotal in assessing geological sites for carbon dioxide (CO 2 ) storage and mitigating potential risks associated with storage and injection. The Energy & Environmental Research Center is exploring innovative carbon storage monitoring techniques at a CO 2 storage site near Beulah, North Dakota. These methods aim to provide lower-impact, faster feedback for commercial carbon capture and storage (CCS) projects. The research includes monitoring CO 2 injection using various approaches: 1) an automated, integrated modular monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance; 4) time-lapse seismic techniques; and 5) advanced wellbore monitoring.

02 PETROLEUM↗

Activating a Natural Fault Zone in the Swiss Alps

One major hurdle for understanding earthquake mechanics are observational limitations. Important phenomena like strain localisation, fault dilation, and fault healing are readily studied in rock mechanical laboratory experiments and with numerical models. At the scale of natural earthquakes, however, these phenomena are often unresolvable, even by state-of-the-art observatories. To overcome this limitation, we are currently building the Earthquake Physics Testbed at the Bedretto Underground Laboratory for Geosciences and Geoenergies (BedrettoLab), an experimental testbed where we can activate an extensively instrumented natural fault zone via hydraulic stimulation. The goal of the Fault Activation and Earthquake Rupture (FEAR) project is to induce earthquakes of up to Mw~1.0 on this exceptionally well characterised and instrumented fault zone. Here we summarize the main scientific goals and current FEAR project status, and present first results from conducted experiments. We discuss how this large-scale experimental approach may allow us to tackle both fundamental science as well as practical questions on earthquake physics, induced seismicity and seismic hazard.

Meier, Men-Andrin↗

Tidally Modulated Glacial Slip and Tremor at Helheim Glacier, Greenland

Numerical modeling of ice sheet motion and hence projections of global sea level rise require information about the evolving subglacial environment, which unfortunately remains largely unknown due to its difficulty of access. Here we advance such subglacial observations by reporting multi–year observations of seismic tremor likely associated with glacier sliding at Helheim Glacier. This association is confirmed by correlation analysis between tremor power and multiple environmental forcings on different timescales. Variations of the observed tremor power indicate that different factors affect glacial sliding on different timescales. Effective pressure may control glacial sliding on long (seasonal/annual) timescales, while tidal forcing modulates the sliding rate and tremor power on short (hourly/daily) timescales. Polarization results suggest that the tremor source comes from an upstream subglacial ridge. This observation provides insights on how different factors should be included in ice sheet modeling and how their timescales of variability play an essential role.

54 ENVIRONMENTAL SCIENCES↗

A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 5 Report Field Scale Stress Modeling)

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686, for the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage.” The purpose of this DOE funding is to develop non-invasive methods to obtain important information about subsurface stresses. The overall research goal of the project was to develop and improve methods for determining the spatial distribution of stresses, including magnitude and orientation of the three principal stress components in the subsurface, based on new methods that extract stress information from seismic data combined with well measurements, extended with well tests and logs, and unified in a numerical model that permits computation of the full stress tensor throughout the domains under investigation.

58 GEOSCIENCES↗

Quantitative Risk Assessment (Topical Report)

Carbon capture, utilization, and storage (CCUS) in geological formations plays a key role in mitigating anthropogenic CO 2 emissions and achieving the aggressive goal of net-zero greenhouse gas emissions. Quantitative risk assessment is crucial for ensuring the safety and reliability of geologic carbon storage (GCS) by evaluating CO 2 migration in subsurface, forecasting potential leakage and induced seismicity risks, and optimizing operational and monitoring plans. We present the use, progress, and research trends of risk assessment from the Southwest Regional Partnership on Carbon Sequestration (SWP) as an example of large-scale CCUS projects in North America. The information provided in this report can help readers understand the significance of risk and uncertainty assessment and apply them effectively in large-scale CCUS projects.

42 ENGINEERING↗

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗

Detection of high-frequency gravitational waves using SRF cavities

Today, apart from some isolated R\&D efforts, there are no GW experiments, yet which explore a large part of the vast frequency range above the LIGO/Virgo band. It is planned to establish an experiment at DESY and Fermilab to search for high-frequency GWs in the frequency range of 10\,kHz to 100\,MHz. The basic idea is to use superconducting radiofrequency (SRF) cavities to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. This paper summarizes the challenging environmental boundary requirements, and the R\&D to operate a cavity using a LLRF system which pushes beyond state-of-the-art accuracy and resolutions and a seismic noise mitigated cryostat at 1.8\,K. The focus of this paper is the warm and cold commissioning of a prototype cavity, built 20 years ago during the MAGO collaboration, and its first measurement in our collaborative research project.

Giaccone, Bianca [Fermilab]↗

Detection of high-frequency gravitational waves using SRF cavities

Today, apart from some isolated R&D efforts, there are no gravitational wave (GW) experiments, yet which explore a large part of the vast frequency range above the LIGO/Virgo band. It is planned to establish an experiment at Deutsches Elektronen-Synchrotron (DESY) and at the Superconducting Quantum Materials and Systems (SQMS) Center at Fermi National Accelerator Laboratory (Fermilab) to search for high-frequency GWs in the frequency range of 10 kHz to 100 MHz. The basic idea is to use superconducting radiofrequency (SRF) cavities to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. This paper summarizes the challenging environmental boundary requirements, and the R&D to operate a cavity using a low level RF (LLRF) system which pushes beyond state-of-the-art accuracy and resolutions and a seismic noise mitigated cryostat at 1.8 K. The focus of this paper is the warm and cold commissioning of a prototype cavity, built 20 years ago during the MAGO collaboration, and its first measurement in our collaborative research project.

Wenskat, M. [DESY; Hamburg U.] (ORCID:000000016546↗

SECARB-USA: Needs Assessment Framework for Storage Complexes (Task 2.1.b)

A team of SECARB storage experts examined 63 formations at 39 sites that were selected to represent the range of diversity of newly assessed, as well as well-advanced, storage prospects in the SECARB region. We inventoried the data needs triggered by the requirements to obtain a Class VI UIC storage permit, the data needs that results from the requirement to create geocellular fluid flow models to support that permit, and by pragmatic and best practice inputs such as public acceptance, regulatory readiness, and pore space leasing. We anchored both the needs inventory and the processes for and cost of meeting the needs with data from 9 sites in in the SECARB area that have advanced far in characterization. Results show that the total cost of characterization prior to obtaining a permit is convergent, because the permit and modeling requirements drive projects to obtain the same types of data for all cases. The high cost data that control cost are 1) drilling, coring, core-testing, logging, sampling and testing a characterization well and 2) collection of a 3-D seismic survey to map reservoir and confining system properties over the area of the plume or the area of elevated pressure. In 5 of our case study sites we determined that one or both of these costs could be avoided because the needs are met by available data.

54 ENVIRONMENTAL SCIENCES↗

Discovering Blind Geothermal Systems in the Great Basin Region: An Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways: All Phases

Most geothermal resources in the Great Basin region of the western USA are blind, and thus the discovery of new commercial-grade systems requires synthesis of favorable characteristics for geothermal activity. The geothermal play fairway concept involves integration of multiple parameters indicative of geothermal activity to identify promising areas for new development. This project integrated multiple datasets to apply the play fairway concept and assess geothermal potential in a large region of the Great Basin in Nevada. It is therefore referred to as the Nevada play fairway project. This project was a strong collaborative effort between several organizations, led by the Nevada Bureau of Mines and Geology at the University of Nevada, Reno, but with key support from the U.S. Geological Survey, ATLAS Geosciences, Inc,, Hi-Q Geophysical, Inc., Lawrence Berkeley National Laboratory, Utah Geological Survey, and Innovative Geothermal Ltd. In Budget Period 1 of this project, available data for nine geologic, geochemical, and geophysical parameters were initially synthesized to produce a new detailed geothermal potential map of 96,000 km2 from west-central to eastern Nevada. These parameters were grouped into subsets and individually weighted to delineate rankings for local permeability, intermediate permeability, regional permeability, and thermal potential, which collectively defined geothermal play fairways (i.e., most likely locations for significant geothermal fluid flow). This initial work was aimed at reducing the risks in regional exploration and therefore facilitating discovery of new commercial-grade systems in blind settings, as well as in areas with surface expressions of geothermal activity. Budget Period 2 of the project involved detailed analysis of some of the most promising areas identified in Phase 1. Twenty-four highly prospective areas, including both known undeveloped systems and previously undiscovered potential blind systems, were identified for further analysis. After reconnaissance of these areas, five of the most promising sites were selected for detailed studies. Multiple techniques were employed in the detailed studies, including geologic mapping, shallow temperature surveys, gravity surveys, Lidar, geochemical studies, seismic reflection analysis, and 3D modeling. The goal of the detailed studies was to identify specific areas with the highest likelihood for high permeability and thermal fluids, such that drill sites could be targeted. Three main sets of predictive maps were generated for each detailed study area: 1) play fairway maps, 2) play fairway error maps, and 3) direct evidence maps. Local- and intermediate-scale permeability models were revised to reflect results of the detailed geologic, geophysical, and geochemical analyses. Budget Period 3 of the project involved more detailed geophysical analyses and temperature-gradient (TG) drilling in southeastern Gabbs Valley and northern Granite Springs Valley, deemed the two most promising sites, with the goal of providing preliminary validation of the play fairway methodology. In southeastern Gabbs Valley, the collocation of a favorable structural setting (displacement transfer zone and fault intersections), Quaternary faults, intersecting and terminating gravity gradients, magnetic low, shallow (2 m) temperature anomaly, low resistivity anomaly, and promising geothermometry from nearby water wells provided evidence for a blind system. Drilling of six TG holes defines an apparent geothermal system at this locality with temperatures as high as 124°C at 152 m. This system is blind, with no surface hot springs, fumaroles, or paleo-geothermal deposits. For northern Granite Springs Valley, a favorable structural setting (termination of a major Quaternary normal fault), terminating gravity gradient, magnetic gradient, newly discovered sinter deposits, nearby warm water wells, previously drilled TG holes in the vicinity, and promising geothermometry suggest a hidden system. Drilling of six new TG holes yields temperatures of ~96°C at ~250 m, suggesting the presence of a geothermal system. Major lessons learned in the course of this project include: 1) initially identified sites commonly include multiple favorable structural settings at a finer scale; 2) promising sites in Cenozoic basins cannot be recognized without detailed geophysical surveys; and 3) play fairway analysis should be refined as the exploration program vectors into the most promising sites and finer-scale data are acquired. In addition to producing copious amounts of data, this project resulted in 16 published papers, 10 abstracts, more than 40 presentations across the U.S. and abroad (including several keynote addresses), 2 Masters theses, and 7 media reports.

15 GEOTHERMAL ENERGY↗

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗