Search NASA⌕ Search

SEARCH · Search NASA

Results for “Waveform Parameters”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Full-Wave and Circuit-Based Simulations of Cable Insulation Aging/Damage Using Time-Frequency Domain Reflectometry

This paper presents a combined full-wave and circuit-based simulation study to investigate the important topic of cable insulation aging and physical damage. The analysis employs the time-frequency domain reflectometry (TFDR) technique that combines the advantages of both time domain reflectometry (TDR) and frequency domain reflectometry (FDR) for improved location identification and feature resolution. Coaxial cable sections with a location containing insulation aging or physical damage are simulated using Ansys High Frequency Structure Simulator (HFSS). Simulated scattering (S) parameters are then imported into the Advanced Design System (ADS) solver for time domain simulation where a Gaussian windowed chirp signal or ‘chirplet’ is launched. Finally, continuous wavelet transform (CWT) is applied in MATLAB to perform time-frequency analysis of the reflected waveforms. The results reveal that a small change in the relative permittivity of the insulation or the presence of a small point of physical damage on the insulation is reflected as a considerable change in the CWT magnitude demonstrating the efficacy of the method.

cable aging, Non Destructive Examination, time-fre↗
DOE OSTI2025

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full Waveform Inversion

Concrete is a major construction material worldwide and plays a crucial role in the nuclear industry. The elastic properties of concrete are prone to change and degrade while in service, as it is often subjected to extreme operational and environmental conditions. An accurate evaluation of concrete's elastic properties is thus essential to ensure structural integrity and safety. This is especially true for concrete in nuclear power plants, where irradiation effects significantly impact concrete mechanical properties. There are various methods to assess these properties, with ultrasound-based techniques showing high potential due to their nondestructive nature, cost-effectiveness, and safety. While several nondestructive evaluation methods exist, most rely on idealizations such as assuming homogeneous material and plane wavefronts. In this work, we address these issues by introducing an ultrasound-based nondestructive method aimed at reconstructing spatially varying images of concrete mechanical properties. By accurately modeling wave physics, including scattering and reflection, we overcome several of the aforementioned idealizations and aim to utilize the full waveform for imaging material properties through depth, resulting in more reliable images. Full waveform inversion (FWI) was first introduced by geophysicists to reconstruct subsurface elastic property images. The goal is to minimize the difference between simulated and recorded wavefield signals, often through gradient-based optimization algorithms. While FWI is primarily conducted using the acoustic approximation of the wave equation, few works focus on elastic FWI, where the goal is to reconstruct images of not only the pressure wave speed but also the shear wave speed and density (or their equivalents). This work explores the potential of using elastic FWI to predict concrete mechanical properties as an initial effort for a more accurate monitoring of concrete conditions in service. Reconstructing images of different elastic parameters enables more specificity and accurate condition diagnosis. This paper will detail this approach and provide examples demonstrating the effectiveness of elastic FWI in reconstructing comprehensive maps of concrete mechanical properties.

42 - ENGINEERING↗
DOE OSTI2016

A study on crustal shear wave splitting in the western part of the Banda arc-continent collision

We analyzed shear wave splitting parameters from local shallow (< 30 km) earthquakes recorded at six seismic stations in the western part of the Banda arc-continent collision. We determined fast polarization and delay time for 195 event-stations pairs calculated from good signal-to-noise ratio waveforms. We observed that there is evidence for shear wave splitting at all stations with dominant fast polarization directions oriented about NE-SW, which are parallel to the collision direction of the Australian plate. However, minor fast polarization directions are oriented around NW-SE being perpendicular to the strike of Timor through. Furthermore, the changes in fast azimuths with the earthquake-station back azimuth suggest that the crustal anisotropy in the study area is not uniform. Splitting delay times are within the range of 0.05 s to 0.8 s, with a mean value of 0.29±0.18 s. Major seismic stations exhibit a weak tendency increasing of delay times with increasing hypocentral distance suggesting the main anisotropy contribution of the shallow crust. In addition, these variations in fast azimuths and delay times indicate that the crustal anisotropy in this region might not only be caused by extensive dilatancy anisotropy (EDA), but also by heterogeneity shallow structure such as the presence of foliations in the rock fabric and the fracture zones associated with active faults.

58 GEOSCIENCES↗
DOE OSTI2022

Semi-Supervised Disaggregation of Load Profiles at Transmission Buses with Significant Behind-the-Meter Solar Generations

It is of imperative interests for regional transmission organizations (RTOs) to effectively extract daily load profiles at transmission buses, which remains a gap in existing technology paradigm. This digest proposes an explicit yet efficient linear estimator, to disaggregate metered load profiles at buses with significant behind-the-meter (BTM) solar generations in a data driven manner. The proposed estimator is based on utility zonal load profiles and proxy solar irradiance profiles, which in reality is the aggregated waveform at each transmission bus and equivalent to the mix of summed load profiles minus actual BTM solar generation. To overcome technical challenges in the lack of “ground truth” and validate the performance of supervised learning algorithms, we propose semi-supervised mechanisms with parameter tuning, and leverage the unique characteristics of zero-crossing points in BTM solar peaking behaviors.

machine leaning, behind-the-meter, power grid anal↗
DOE OSTI2015

Re-appraisal and extension of the Gratton-Vargas two-dimensional analytical snowplow model of plasma focus. II. Looking at the singularity

The Gratton-Vargas snowplow model, recently revisited and expanded [S. K. H. Auluck, Phys. Plasmas 20, 112501 (2013)], has given rise to significant new insights into some aspects of the Dense Plasma Focus (DPF), in spite of being a purely kinematic description having no reference to plasma phenomena. It is able to provide a good fit to the experimental current waveforms in at least 4 large facilities. It has been used for construction of a local curvilinear frame of reference, in which conservation laws for mass, momentum, and energy can be reduced to effectively-one-dimensional hyperbolic conservation law equations. Its utility in global parameter optimization of device parameters has been demonstrated. These features suggest that the Gratton-Vargas model deserves a closer look at its supposed limitations near the singular phase of the DPF. This paper presents a discussion of its development near the device axis, based on the original work of Gratton and Vargas, with some differences. It is shown that the Gratton-Vargas partial differential equation has solutions for times after the current singularity, which exhibit an expanding bounded volume (which can serve as model of an expanding plasma column) and decreasing dynamic inductance of the discharge, in spite of having no built-in hydrodynamics. This enables the model to qualitatively reproduce the characteristic shape of the current derivative in DPF experiments without reference to any plasma phenomena, such as instabilities, anomalous resistance, or reflection of hydrodynamic shock wave from the axis. The axial propagation of the solution exhibits a power-law dependence on the dimensionless time starting from the time of singularity, which is similar to the power-law relations predicted by theory of point explosions in ideal gases and which has also been observed experimentally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

The 26 December 2004 tsunami source estimated from satellite radar altimetry and seismic waves

The 26 December 2004 Indian Ocean tsunami was the first earthquake tsunami of its magnitude to occur since the advent of both digital seismometry and satellite radar altimetry. Both have independently recorded the event from different physical aspects. The seismic data has then been used to estimate the earthquake fault parameters, and a three-dimensional ocean-general-circulation-model (OGCM) coupled with the fault information has been used to simulate the satellite-observed tsunami waves. Here we show that these two datasets consistently provide the tsunami source using independent methodologies of seismic waveform inversion and ocean modeling. Cross-examining the two independent results confirms that the slip function is the most important condition controlling the tsunami strength, while the geometry and the rupture velocity of the tectonic plane determine the spatial patterns of the tsunami.

ocean general circulation model (OGCM)↗
DOE OSTI2023

Fourier-DeepONet: Fourier-enhanced deep operator networks for full waveform inversion with improved accuracy, generalizability, and robustness

In this article, full waveform inversion (FWI) infers the subsurface structure information from seismic waveform data by solving a non-convex optimization problem. Data-driven FWI has been increasingly studied with various neural network architectures to improve accuracy and computational efficiency. Nevertheless, the applicability of pre-trained neural networks is severely restricted by potential discrepancies between the source function used in the field survey and the one utilized during training. Here, we develop a Fourier-enhanced deep operator network (Fourier-DeepONet) for FWI with the generalization of seismic sources, including the frequencies and locations of sources. Specifically, we employ the Fourier neural operator as the decoder of DeepONet, and we utilize source parameters as one input of Fourier-DeepONet, facilitating the resolution of FWI with variable sources. To test Fourier-DeepONet, we develop three new and realistic FWI benchmark datasets (FWI-F, FWI-L, and FWI-FL) with varying source frequencies, locations, or both. Our experiments demonstrate that compared with existing data-driven FWI methods, Fourier-DeepONet obtains more accurate predictions of subsurface structures in a wide range of source parameters. Moreover, the proposed Fourier-DeepONet exhibits superior robustness when handling data with Gaussian noise or missing traces and sources with Gaussian noise, paving the way for more reliable and accurate subsurface imaging across diverse real conditions.

42 ENGINEERING↗
DOE OSTI2019

Detection of random noise and anatomy of continuous seismic waveforms in dense array data near Anza California

We develop a methodology to separate continuous seismic waveforms into random noise (RN), not random noise (NRN) produced by earthquakes, wind, traffic and other sources of ground motions, and an undetermined mixture of signals. The analysis is applied to continuous records from a dense seismic array on the San Jacinto fault zone. To detect RN signals, we cut hourly waveforms into non-overlapping 1 s time windows and apply cross-correlations to separate RN candidates from outliers. The cross-correlation coefficients between different RN candidates fall into a tight range (i.e. 0.09–0.35), while cross-correlation coefficients of RN candidates with NRN signals (e.g. seismic or air-traffic events) are lower. The amplitude spectra of RN candidates have a well-defined level, while the amplitude spectra of other signals deviate from that level. Using these properties, we examine the amplitude spectra of moving time windows and cross-correlation coefficients with RN templates in each hour. The hourly RN is quasi-stationary and the results cluster tightly in the parameter space of cross-correlation coefficients and L2 norm deviations from the mean spectra of RN candidates. Time windows with parameters in this tight cluster are identified as RN, windows that deviate significantly from the RN cluster are identified as NRN and windows with values in between are identified as mixed signals. Several iterations on each hourly data are used to update and stabilize the selection of RN templates and mean noise spectra. For the days examined, the relative fractions of RN, NRN and mixed signals in local day (night) times are about 26 (42 percent), 40 (33 percent) and 34 percent (25 percent), respectively.

Geochemistry & Geophysics↗
DOE OSTI2014

Performance Bounds on Micro-Doppler Estimation and Adaptive Waveform Design Using OFDM Signals

We analyze the performance of a wideband orthogonal frequency division multiplexing (OFDM) signal in estimating the micro-Doppler frequency of a target having multiple rotating scatterers (e.g., rotor blades of a helicopter, propellers of a submarine). The presence of rotating scatterers introduces Doppler frequency modulation in the received signal by generating sidebands about the transmitted frequencies. This is called the micro-Doppler effects. The use of a frequency-diverse OFDM signal in this context enables us to independently analyze the micro-Doppler characteristics with respect to a set of orthogonal subcarrier frequencies. Therefore, to characterize the accuracy of micro-Doppler frequency estimation, we compute the Cram er-Rao Bound (CRB) on the angular-velocity estimate of the target while considering the scatterer responses as deterministic but unknown nuisance parameters. Additionally, to improve the accuracy of the estimation procedure, we formulate and solve an optimization problem by minimizing the CRB on the angular-velocity estimate with respect to the transmitting OFDM spectral coefficients. We present several numerical examples to demonstrate the CRB variations at different values of the signal-to-noise ratio (SNR) and the number of OFDM subcarriers. The CRB values not only decrease with the increase in the SNR values, but also reduce as we increase the number of subcarriers implying the significance of frequency-diverse OFDM waveforms. The improvement in estimation accuracy due to the adaptive waveform design is also numerically analyzed. Interestingly, we find that the relative decrease in the CRBs on the angular-velocity estimate is more pronounced for larger number of OFDM subcarriers.

Sen, Satyabrata↗
DOE OSTI2021

Void-free Copper Electrodeposition in High Aspect Ratio, Full Wafer Thickness Through-Silicon Vias with Endpoint Detection

High density interconnects are required for increased input/output for microelectronics applications, incentivizing the development of Cu electrochemical deposition (ECD) processes for high aspect ratio through-silicon vias (TSVs). This work outlines Cu ECD processes for 62.5 μ m diameter TSVs, etched into a 625 μ m thick silicon substrate, a 10:1 aspect ratio. Cu ECD in high aspect ratio features relies on a delicate balance of electrolyte composition, solution replenishment, and applied voltage. Implementing a CuSO 4 -H 2 SO 4 electrolyte, which contains suppressor and a low chloride concentration, allows for a tunable relationship between applied voltage and localized deposition in the vias. A stepped potential waveform was applied to move the Cu growth front from the bottom of the via to the top. Sample characterization was performed through mechanical cross-sections and X-ray computed tomography (CT) scans. The CT scans revealed small seam voids in the Cu electrodeposit, and process parameters were tuned accordingly to produce void-free Cu features. During the voltage-controlled experiments, measured current data showed a characteristic current minimum, which was identified as an endpoint detection method for Cu deposition in these vias. We believe this is the first report of this novel endpoint detection method for TSV filling.

42 ENGINEERING↗
NASA NTRS1986

Intraplate earthquakes and the state of stress in oceanic lithosphere

The dominant sources of stress relieved in oceanic intraplate earthquakes are investigated to examine the usefulness of earthquakes as indicators of stress orientation. The primary data for this investigation are the detailed source studies of 58 of the largest of these events, performed with a body-waveform inversion technique of Nabelek (1984). The relationship between the earthquakes and the intraplate stress fields was investigated by studying, the rate of seismic moment release as a function of age, the source mechanisms and tectonic associations of larger events, and the depth-dependence of various source parameters. The results indicate that the earthquake focal mechanisms are empirically reliable indicators of stress, probably reflecting the fact that an earthquake will occur most readily on a fault plane oriented in such a way that the resolved shear stress is maximized while the normal stress across the fault, is minimized.

Bergman, Eric A.↗
NASA NTRS1995

Geometric adjustments to account for eye eccentricity in processing horizontal and vertical eye and head movement data

Neglecting the eccentric position of the eyes in the head can lead to erroneous interpretation of ocular motor data, particularly for near targets. We discuss the geometric effects that eye eccentricity has on the processing of target-directed eye and head movement data, and we highlight two approaches to processing and interpreting such data. The first approach involves determining the true position of the target with respect to the location of the eyes in space for evaluating the efficacy of gaze, and it allows calculation of retinal error directly from measured eye, head, and target data. The second approach effectively eliminates eye eccentricity effects by adjusting measured eye movement data to yield equivalent responses relative to a specified reference location (such as the center of head rotation). This latter technique can be used to standardize measured eye movement signals, enabling waveforms collected under different experimental conditions to be directly compared, both with the measured target signals and with each other. Mathematical relationships describing these approaches are presented for horizontal and vertical rotations, for both tangential and circumferential display screens, and efforts are made to describe the sensitivity of parameter variations on the calculated results.

NASA Discipline Neuroscience↗
NASA NTRS1999

Development of a flow feedback pulse duplicator system with rhesus monkey arterial input impedance characteristics

An in vitro pulsatile pump flow system that is capable of producing physiologic pressures and flows in a mock circulatory system tuned to reproduce the first nine harmonics of the input impedance of a rhesus monkey was developed and tested. The system was created as a research tool for evaluating cardiovascular function and for the design, testing, and evaluation of electrical-mechanical cardiovascular models and chronically implanted sensors. The system possesses a computerized user interface for controlling a linear displacement pulsatile pump in a controlled flow loop format to emulate in vivo cardiovascular characteristics. Evaluation of the pump system consisted of comparing its aortic pressure and flow profiles with in vivo rhesus hemodynamic waveforms in the time and frequency domains. Comparison of aortic pressure and flow data between the pump system and in vivo data showed good agreement in the time and frequency domains, however, the pump system produced a larger pulse pressure. The pump system can be used for comparing cardiovascular parameters with predicted cardiovascular model values and for evaluating such items as vascular grafts, heart valves, biomaterials, and sensors. This article describes the development and evaluation of this feedback controlled cardiovascular dynamics simulation modeling system.

Non-NASA Center↗
NASA NTRS2006

Electrostatic Screen for Transport of Martian and Lunar Regolith

The Martian and Lunar Regolith contain fine particulate including those in the size range from 0.5 to 200 micron [1-2]. Martian dust can be transported and deposited by Aeolian processes, including "Dust Devils". Due to the ultra high vacuum (10e-12 Torr), transport of dust on the Moon is solely a result of collision/ballistic motion. Dust obscuration of solar cells is one of the primary factors limiting the duration of Martian missions, including the Mars Exploration Rovers. Dust contamination in vacuum seals is one of the primarily factors that limited lunar excursions during the Apollo missions. Controlled transportation of dust on Mars and the Moon is important for many reasons, including both contamination mitigation and in situ resource utilization (ISRU). Since both the monopole and dipole electrostatic moments result in non-trivial forces on particles in an electrostatic field, dust particles, whether charged or not, can be transported by electrostatic fields. In the electrostatic screen, alternating waveforms of voltage applied to patterned grids of electrodes will transport dust. The authors will show that the canonical methods for transporting dust via electrostatic screen can be readily applied to transport of Martian and Lunar regolith. Experiments have been performed in ambient, low humidity, Martian, and Lunar conditions. Screen parameters have been examined for application to each regolith, such as grid spacing, trace width, grid voltage, pulse pattern, pulse frequency, and coating type. The authors have also developed an electrostatic screen based on optically transparent conductors that can be placed over solar arrays, windows, visors, lenses, etc.

Immer, C.↗
NASA NTRS2004

Hollow Cathode and Keeper-region Plasma Measurements Using Ultra-fast Miniature Scanning Probes

In order to support the development of comprehensive performance and life models for future deep space missions that will utilize ion thrusters, we have undertaken a study of the plasma structure in hollow cathodes using an new pneumatic scanning probe diagnostic. This device is designed to insert a miniature probe directly into the hollow cathode orifice from either the upstream insert region in the interior of the hollow cathode, or from the downstream keeper-plasma region at the exit of the hollow cathode, to provide complete axial profiles of the discharge plasma parameters. Previous attempts to diagnose this region with probes was Limited by the melting of small probes in the intense discharge near the orifice, or caused significant perturbation of the plasma by probes large enough to survive. Our new probe is extremely compact, and when configured as a single Langmuir probe, the ceramic tube insulator is only 0.5mm in diameter and the current collecting conductor has a total area of 0.002 cm2. A series of current-voltage characteristics are obtained by applying a rapid sawtooth voltage waveform to the probe as it is scanned by the pneumatic actuator into and out of the plasma region, The bellow-sealed pneumatic drive scans the probe 4 cm in the cathode insert region and 10 cm in the anode/keeper plasmas region at average speeds of about 1 mm/msec, and the residence time at the end of the insertion stroke in the densest part of the plasma near the orifice is measured to be only 10 msec. Since the voltage sweep time is fast compared to the motion of the probe, axial profiles of the plasma density, temperature and potential with reasonable spatial resolution are obtained. Measurements of the internal cathode pressures and the axial plasma-parameter profiles for a hollow cathode operating at discharge currents of up to 35 A in xenon will be presented.

plasma generation↗
DOE OSTI2026

The role of source geometry and atmospheric propagation in global bolide infrasound detectability

Global infrasound monitoring provides a persistent means of detecting energetic bolide atmospheric entries, complementing optical observations and extending coverage over remote regions. We present a global assessment of the physical factors governing bolide infrasound detectability by correlating 623 bolide events reported by the Center for Near-Earth Object Studies between 2007 and 2025 with waveform data from the International Monitoring System. We identify 311 events with confirmed infrasound detections, corresponding to a detection rate of approximately 50%, substantially higher than inferred from earlier surveys, reflecting both the maturation of the global infrasound network and advances in automated, multi-frequency array processing. Analysis of flight parameters shows that infrasound detectability is selective rather than uniform across the bolide population. Detected events are preferentially associated with steeper entry angles and lower-altitude energy deposition, while shallow, high-altitude trajectories are less consistently observed. Very high-energy events remain detectable regardless of geometry, but for the more common lower-energy regime, observability depends on specific combinations of entry parameters and propagation conditions. This geometric dependence persists across comparable energy ranges and atmospheric conditions, indicating that entry angle exerts a primary control on detectability, with energy and propagation acting as secondary modulating factors. Furthermore, these results provide new physical constraints on bolide-atmosphere interactions and improve interpretation of global infrasound observations for planetary defense and atmospheric-entry studies.

Bolides↗
NASA NTRS2017

Chorus Whistler Wave Source Scales As Determined From Multipoint Van Allen Probe Measurements

Whistler mode chorus waves are particularly important in outer radiation belt dynamics due to their key role in controlling the acceleration and scattering of electrons over a very wide energy range. The key parameters for both nonlinear and quasi-linear treatment of wave-particle interactions are the temporal and spatial scales of the wave source region and coherence of the wave field perturbations. Neither the source scale nor the coherence scale is well established experimentally, mostly because of a lack of multipoint VLF waveform measurements. We present an unprecedentedly long interval of coordinated VLF waveform measurements (sampled at 16384 s(exp -1)) aboard the two Van Allen Probes spacecraft-9 h (0800-1200 UT and 1700-2200 UT) during two consecutive apogees on 15 July 2014. The spacecraft separations varied from about 100 to 5000 km (mostly radially); measurements covered an L shell range from 3 to 6; magnetic local time 0430-0900, and magnetic latitudes were approximately 15 and approximately 5 deg during the two orbits. Using time-domain correlation techniques, the single chorus source spatial extent transverse to the background magnetic field has been determined to be about 550-650 km for upper band chorus waves with amplitudes less than 100 pT and up to 800 km for larger amplitude, lower band chorus waves. The ratio between wave amplitudes measured on the two spacecraft is also examined to reveal that the wave amplitude distribution within a single chorus element generation area can be well approximated by a Gaussian exp(-0.5 x r (exp 2)/r(sub 0)(exp 2)), with the characteristic scale r(sub 0) around 300 km. Waves detected by the two spacecraft were found to be coherent in phase at distances up to 400 km.

Agapitov, O.↗
DOE OSTI2022

On the measurement of S diff splitting caused by lowermost mantle anisotropy

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (S diff ) are potentially strongly influenced by lowermost mantle anisotropy. S diff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. S diff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SH diff –SV diff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of S diff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which S diff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on S diff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SH diff waves. In particular, we investigate how SH diff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to S diff splitting that can be clearly distinguished from other waveform effects. Here, these results enable us to lay out a strategy for the analysis of S diff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of S diff and the influence of source-side anisotropy. We demonstrate our S diff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained S diff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

54 ENVIRONMENTAL SCIENCES↗