Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency shift”

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,531 records · Page 85

Explaining Polarization Reversals in STEREO Wave Data

Recently Breneman et al. reported observations of large amplitude lightning and transmitter whistler mode waves from two STEREO passes through the inner radiation belt (L<2). Hodograms of the electric field in the plane transverse to the magnetic field showed that the transmitter waves underwent periodic polarization reversals. Specifically, their polarization would cycle through a pattern of right-hand to linear to left-hand polarization at a rate of roughly 200 Hz. The lightning whistlers were observed to be left-hand polarized at frequencies greater than the lower hybrid frequency and less than the transmitter frequency (21.4 kHz) and right-hand polarized otherwise. Only righthand polarized waves in the inner radiation belt should exist in the frequency range of the whistler mode and these reversals were not explained in the previous paper. We show, with a combination of observations and simulated wave superposition, that these polarization reversals are due to the beating of an incident electromagnetic whistler mode wave at 21.4 kHz and linearly polarized, symmetric lower hybrid sidebands Doppler-shifted from the incident wave by +/-200 Hz. The existence of the lower hybrid waves is consistent with the parametric decay mechanism of Lee and Kuo whereby an incident whistler mode wave decays into symmetric, short wavelength lower hybrid waves and a purely growing (zero-frequency) mode. Like the lower hybrid waves, the purely growing mode is Doppler-shifted by 200 Hz as observed on STEREO. This decay mechanism in the upper ionosphere has been previously reported at equatorial latitudes and is thought to have a direct connection with explosive spread F enhancements. As such it may represent another dissipation mechanism of VLF wave energy in the ionosphere and may help to explain a deficit of observed lightning and transmitter energy in the inner radiation belts as reported by Starks et al.

Breneman, A.↗

Timing the Pulsations of The Accreting Millisecond Pulsar SAX J1808.4-3658 During Its 2019 Outburst

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

Peter Bult↗

Photogalvanic Effects in Surface States of Topological Insulators under Perpendicular Magnetic Fields

We present a theoretical study of the nonlinear magneto-optical shift conductivity in the surface states of the prototypical topological insulator Bi$_2$Se$_3$ under a perpendicular quantizing magnetic field. By describing the electronic states as Landau levels and using a perturbative approach, we derive the microscopic expression for the shift conductivity $\sigma^{(2);\alpha\beta\gamma}(-\omega,\omega)$, where $\alpha,\beta,\gamma=\pm$ stand for the circular polarization of light and $\omega$ is the light frequency; the spectra are further decomposed into contributions from the interband and intraband optical transitions, for which the selection rules are identified. Considering that the system possesses $C_3$ point group of symmetry, the nonzero components of the conductivity tensor are $\sigma^{(2);-++}=[\sigma^{(2);+--}]^\ast$. Therefore, a pure circularly polarized light generates zero shift current. In the clean limit, the conductivities are nonzero only for discrete photon energies because of the discrete Landau levels and energy conservation, and they become Lorentzian lineshapes with the inclusion of damping, which relaxes the condition of energy conservation. The dependence of the spectra on the damping parameters, the magnetic fields, and the chemical potentials is investigated in detail. Our results reveal that the shift current is highly tunable by the chemical potential and the magnetic field. Furthermore, these results underscore the potential of topological insulators for tunable, strong nonlinear magneto-optical applications.

Landau level↗

One-to-one mapping between the electromagnetic modes of cylindrical and coaxial half-wave cavities

Design of radio frequency (RF) couplers and diagnostics require a good understanding of the electromagnetic mode patterns of RF cavities. This study investigates the adiabatic transformation of transverse magnetic (TM) modes in a cylindrical cavity into transverse electromagnetic (TEM) modes of a coaxial cavity by gradually introducing an inner conductor. Using CST Studio Suite, we simulate the eigenmode evolution as the geometry transforms from a pure cylindrical to a coaxial configuration. We track the behavior of TM010 through TM014 modes to observe the continuous evolution into the corresponding TEM0 through TEM4 modes of the coaxial cavity. The process is governed by the evolution of the electric field orientation as the geometry shifts, enabling the axial TM fields to reorient into the radial electric field configuration of TEM modes. Field patterns, eigen-frequencies, and mode indentities are analyzed throughtout the transition. The results provide simulation-based evidence that TM to TEM conversion occurs without generation of newer eigenmodes, offering a valuable insight into the design of transition regions in superconducting RF (SRF) systems and provides a foundation for experimental validation.

Accelerator Physics↗

Tangential mode combustion instability of solid propellant rocket motors. 2: Pressure, frequency, and burning rate histories

The relationships among excursions of mean chamber pressure, high frequency pressure oscillations and accelerations associated with tangential mode combustion instability are considered. A data reduction procedure was devised to utilize the precise pressure time data as input to a numerical solution of the chamber mass balance relationships. The results indicate that both elevated burning rate induced by the tangential wave motions and variable nozzle blockage due to passage of vortical flow through the nozzle throat must act together to generate the large mean pressure peaks observed in the tests. There is strong correlation between the amplitude of wave motions and the level of mean pressure rise. Evidence for the presence of traveling tangential waves at the time of the D.C. shift was found by use of dual Kistler pressure transducers mounted at the head end of the burner. A very unusual frequency history was found to accompany the high amplitude disturbances. This was consistent from test to test and may represent an additional source of information regarding the complicated flow combustion interactions within the propellant burning zone.

Flandro, G. A.↗

Fabrication and Test of an Optical Magnetic Mirror

Traditional mirrors at optical wavelengths use thin metalized or dielectric layers of uniform thickness to approximate a perfect electric field boundary condition. The electron gas in such a mirror configuration oscillates in response to the incident photons and subsequently re-emits fields where the propagation and electric field vectors have been inverted and the phase of the incident magnetic field is preserved. We proposed fabrication of sub-wavelength-scale conductive structures that could be used to interact with light at a nano-scale and enable synthesis of the desired perfect magnetic-field boundary condition. In a magnetic mirror, the interaction of light with the nanowires, dielectric layer and ground plate, inverts the magnetic field vector resulting in a zero degree phase shift upon reflection. Geometries such as split ring resonators and sinusoidal conductive strips were shown to demonstrate magnetic mirror behavior in the microwave and then in the visible. Work to design, fabricate and test a magnetic mirror began in 2007 at the NASA Goddard Space Flight Center (GSFC) under an Internal Research and Development (IRAD) award Our initial nanowire geometry was sinusoidal but orthogonally asymmetric in spatial frequency, which allowed clear indications of its behavior by polarization. We report on the fabrication steps and testing of magnetic mirrors using a phase shifting interferometer and the first far-field imaging of an optical magnetic mirror.

Hagopian, John G.↗

Light-induced shift current vortex crystals in moiré heterobilayers

Transition metal dichalcogenide (TMD) moiré superlattices provide an emerging platform to explore various light-induced phenomena. Recently, the discoveries of novel moiré excitons have attracted great interest. The nonlinear optical responses of these systems are however still underexplored. Here, we report investigation of light-induced shift currents (a second-order response generating DC current from optical illumination) in the WSe 2 /WS 2 moiré superlattice. We identify a striking phenomenon of the formation of shift current vortex crystals—i.e., two-dimensional periodic arrays of moiré-scale current vortices and associated magnetic fields with remarkable intensity under laboratory laser setup. Furthermore, we demonstrate high optical tunability of these current vortices—their location, shape, chirality, and magnitude can be tuned by the frequency, polarization, and intensity of the incident light. Electron-hole interactions (excitonic effects) are found to play a crucial role in the generation and nature of the shift current intensity and distribution. Our findings provide a promising all-optical control route to manipulate nanoscale shift current density distributions and magnetic field patterns, as well as shed light on nonlinear optical responses in moiré quantum matter and their possible applications.

36 MATERIALS SCIENCE↗

Spread Spectrum Time Domain Reflectometry (SSTDR) and Frequency Domain Reflectometry (FDR) for Detection of Cable Anomalies Using Machine Learning

Cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, continued use of these cables must shift to a performance-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. A variety of cable tests are available and are commonly applied during outages when the cables can be taken out of service. Frequency domain reflectometry (FDR) is one of these test methods that is being more broadly accepted and used because it not only detects anomalies along the cable with a low-voltage signal that does not stress the cable insulation, but the technique also locates the anomalies. This supports follow-up local inspection and local repair or partial replacement of a damaged cable segment. Currently, FDR testing is only applied to cables that are taken out of service since the test instrument would be damaged by operational voltages.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Spread Spectrum Time Domain Reflectometry (SSTDR) and Frequency Domain Reflectometry (FDR) for Detection of Cable Anomalies Using Machine Learning

Cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, continued use of these cables must shift to a performance-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. A variety of cable tests are available and are commonly applied during outages when the cables can be taken out of service. Frequency domain reflectometry (FDR) is one of these test methods that is being more broadly accepted and used because it not only detects anomalies along the cable with a low-voltage signal that does not stress the cable insulation, but the technique also locates the anomalies. This supports follow-up local inspection and local repair or partial replacement of a damaged cable segment. Currently, FDR testing is only applied to cables that are taken out of service since the test instrument would be damaged by operational voltages. A related technology that has found some acceptance in the aircraft and rail industry is spread spectrum time domain reflectometry (SSTDR). This technology has been implemented with a custom commercial instrument by LiveWire Innovation Inc. that is designed to operate on live cables up to 1000 volts. One of the main conclusions of a previous effort was that cable reflectometry plots can be difficult for humans to analyze due to baseline noise, low or noisy anomaly response peaks, or large responses from cable ends. Detection of cable anomalies for many of these frequencies and test conditions was challenging for manual analysis. This presented an ideal opportunity for ML analysis to distinguish undamaged cable indications from anomalous cable indications. This research discusses application of machine learning (ML) to reflectometry cable test methods. The goal was to assess feasibility to distinguish undamaged cable reflectometry responses from damaged or anomalous cable reflectometry responses. The assessment considered the 3 instruments, multiple frequency bandwidths from each instrument, multiple cable anomalies and test conditions, and both supervised and unsupervised ML approaches. Although approaches and analysis methods were not identical or directly comparable, both outputs were encouraging. The unsupervised prediction weighted accuracy was assessed by instrument and by frequency. It performed better at high frequencies with the highest prediction accuracy of 0.84 for the higher frequency FDR, 0.79 for the 48-MHz LiveWire SSTDR, and 0.77 for 300-MHz PNNL SSTDR. The initial weighted accuracy average across all frequencies for using supervised ML was 0.56 to 0.68. The supervised analysis was repeated with noisier training data removed resulting in weighted accuracies of 0.69 to 0.87. These weighted accuracies are not directly comparable due to differences in the supervised and unsupervised analysis details but do indicate an encouraging trend. Even with limited and unbalanced data, strong prediction accuracies seem encouraging for further work including more data under a wider range of conditions.

42 ENGINEERING↗

Shape-Selective Synthesis of Alkylcyclopentenyl Cations in Zeolites and Spectroscopic Distinction of Constitutional Isomers

Alkylcyclopentenyl cations belong to the long-lived intermediates that make up the “hydrocarbon pool” during the catalytic conversion of methanol on zeolites, and recent works show that such cations contribute to olefin and aromatics formation. From liquid phase chemistry, two types of alkylcyclopentenyl cations are known and distinguished by the substituent at the central carbon (C-2) of the allylic system: a more stable type with a methyl group at the C-2 and a less stable type with hydrogen at the C-2. The following three linked objectives are pursued: (i) IR spectroscopic distinction of the different substitution patterns at the allylic system of alkylcyclopentenyl cations, (ii) the role of the zeolite framework in determining the substitution pattern, and (iii) identification of alkylcyclopentenyl cations from precursors relevant to methanol-to-olefins conversion. UV–vis and IR spectroscopy are applied in situ to characterize alkylcyclopentenyl cations produced by adsorption of a pentaalkylcyclopentadiene or by adsorption and thermally induced cyclization of 2,6-dimethyl-2,4,6-octatriene. Prior knowledge of electronic spectra is combined with DFT-computed and experimental IR spectra to establish the frequency of the C–H vibration of the C-2 hydrogen-substituted type, and a characteristic red shift of the asymmetric allylic stretching vibration (Δν ≈ -20 to -30 cm –1 ) after replacing hydrogen by methyl at the central carbon of the allylic system. Although DFT demonstrates that both types of ions fit into medium- and large-pore zeolites and that the C-2 methyl-substituted type is thermodynamically favored even in the pores of the considered zeolites, formation of alkylcyclopentenyl cations is found by UV–vis and IR spectra to be shape-selective. The bulkier C-2 methyl-substituted type is detected in large-pore zeolites (MOR, BEA) and in the intersections of medium-pore zeolites (MFI), whereas in channels of medium size (TON), the less bulky C-2 hydrogen-substituted type is exclusively formed. Finally, the type of ion formed and its overall size are dictated by the zeolite framework and are independent of the precursor; the same type of alkylcyclopentenyl cation as found through cyclization of dimethyloctatriene could be generated from ethene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-fidelity qutrit entangling gates for superconducting circuits

Ternary quantum information processing in superconducting devices poses a promising alternative to its more popular binary counterpart through larger, more connected computational spaces and proposed advantages in quantum simulation and error correction. Although generally operated as qubits, transmons have readily addressable higher levels, making them natural candidates for operation as quantum three-level systems (qutrits). Recent works in transmon devices have realized high fidelity single qutrit operation. Nonetheless, effectively engineering a high-fidelity two-qutrit entanglement remains a central challenge for realizing qutrit processing in a transmon device. In this work, we apply the differential AC Stark shift to implement a flexible, microwave-activated, and dynamic cross-Kerr entanglement between two fixed-frequency transmon qutrits, expanding on work performed for the ZZ interaction with transmon qubits. We then use this interaction to engineer efficient, high-fidelity qutrit CZ † and CZ gates, with estimated process fidelities of 97.3(1)% and 95.2(3)% respectively, a significant step forward for operating qutrits on a multi-transmon device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Insights into the broadband emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Aims.We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE)γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). Methods.We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC,Fermi-LAT,NuSTAR,Swift(XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario. Results.During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2–1 TeV flux of 20%–50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHEγ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in theγ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies.NuSTARdata allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree.

Astronomy & Astrophysics↗

Integrating Radio-Over-Fiber Communication System and BOTDR Sensor System

In this paper, we propose and experimentally demonstrate for the first time, the integration of a radio-over-fiber (RoF) communication system and a Brillouin optical time-domain reflectometry (BOTDR) distributed sensor system using a single optical fiber link. In this proof-of-concept integrated system, the communication system is composed of three modulation formats of quadrature phase-shift keying (QPSK), 16-quadrature amplitude modulation (16-QAM) and 64-QAM, which are modulated onto an orthogonal frequency division multiplexing (OFDM) signal. Whereas, the BOTDR sensor system is used for strain and/or temperature monitoring over the fiber distance with a spatial resolution of 5 m using a 25 km single-mode silica fiber. The error vector magnitude (EVM) is analyzed in three modulation formats in the presence of various BOTDR input pump powers. Using QPSK modulation, optimized 18 dBm sensing and 10 dBm data power, the measured EVM values with and without bandpass filter are 3.5% and 14.5%, respectively. The proposed system demonstrates a low temperature measurement error (±0.49 °C at the end of 25 km) and acceptable EVM values, which were within the 3GPP requirements. The proposed integrated system can be effectively applied for practical applications, which significantly reduces the fiber infrastructure cost by effective usage of a single optical fiber link.

42 ENGINEERING↗

Computing Bandwidths Of PCM/PSK/PM Signals

Analysis of spectral properties of class of pulse-code-modulated, phase-shift-keyed (PSK) radio signals shows one can obtain approximate, simple, closed-form equations for widths of frequency bands occupied by such signals. These equations reduce amount of computation and time needed to determine allocations of frequency and power among competing radio-communication systems or to determine whether two different radio-communication systems compatible.

Nguyen, Tien M.↗

On phosphorus chemical shift anisotropy in metal (IV) phosphates and phosphonates: A 31 P NMR experimental study

The phosphorus chemical shift anisotropies, 31 PΔcs, and asymmetry parameters η were measured by the 31 P{ 1 H} NMR experiments in static and low-frequency spinning samples of the zirconium phosphates and phosphonates and also in the mixed Zr (IV)/Sn (IV) phosphate/phosphonate material. The data obtained have shown a 111 connectivity in the HPO 4 and PO 3 groups, which does not change at modification and intercalation of the materials. The 31 PΔcs values of the phosphonate groups (43–49 ppm) significantly surpass the values characterizing the HPO 4 groups (23–37 ppm). The 31 P Δcs values obtained for the metal (IV) phosphates were discussed in terms of P-O distances. Finally, the 31 P chemical shift anisotropy parameters can help at elucidation of local structures in phosphate and phosphonate materials.

1H MAS NMR↗

dc electrical conductivity in strongly magnetized plasmas

A generalized Ohm's law is derived to treat strongly magnetized plasmas in which the electron gyrofrequency significantly exceeds the electron plasma frequency. Here, strong magnetization of electrons causes the frictional drag between electrons and ions due to Coulomb collisions to shift, producing an additional transverse resistivity term in the generalized Ohm's law that is perpendicular to both the current ($J$) and the Hall ($JxB$) direction. In the limit of very strong magnetization the parallel resistivity is found to increase by a factor of 3/2, and the perpendicular resistivity by a factor of $\frac{3}{4} \ln (m_i/m_e)$, where $m_i$ and $m_e$ are the ion and electron masses. These results suggest that strong magnetization significantly changes the magnetohydrodynamic evolution of a plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bulk Fermi surfaces of the Dirac type-II semimetallic candidate NiTe 2

Here, we present a study on the Fermi surface of the Dirac type-II semimetallic candidate NiTe 2 via the temperature and angular dependence of the de Haas–van Alphen (dHvA) effect measured in single crystals grown through Te flux. In contrast to its isostructural compounds like PtSe 2 , band-structure calculations predict NiTe 2 to display a tilted Dirac node very close to its Fermi level that is located along the Γ to A high-symmetry direction within its first Brillouin zone. The angular dependence of the dHvA frequencies is found to be in agreement with the first-principles calculations when the electronic bands are slightly shifted with respect to the Fermi level (ε F ), and therefore provide support for the existence of a Dirac type-II node in NiTe 2 . Nevertheless, we observed mild disagreements between experimental observations and density functional theory (DFT) calculations as, for example, nearly isotropic and light experimental effective masses. This indicates that the dispersion of the bands is not well captured by DFT. Despite the coexistence of Dirac-like fermions with topologically trivial carriers, samples of the highest quality display an anomalous and large either linear or sublinear magnetoresistivity. Finally, this suggests that Lorentz invariance breaking Dirac-like quasiparticles dominate the carrier transport in this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Well Stimulation Treatment in California: Evaluation of Disclosure Data, May 2015 – October 2019

In this document, disclosure data from 1,228 well stimulations occurring May 2015 to October 2019 in California are evaluated. This evaluation updates a previous study that was based on 618 well stimulations occurring May 2015 to June 2016 (Stringfellow, W.T., Camarillo, M.K., and Jordan, P. 2017, Status of Well Stimulation in California Since Implementation of SB-4 Regulations, Berkeley National Laboratory, Berkeley, CA). While the goal of the previous study was to provide a summary of information obtained since passage of California Senate Bill No. 4—Oil and Gas, the goal of the current evaluation is to provide a basis for studying chemical indicators that could potentially be used to detect aquifer contamination. We evaluate chemical indicators is a separate document using the data contained herein as well as data from monitoring wells located within oil fields. Compared with the data used by Stringfellow et al. 2017, the current evaluation is based on a more diverse data set. While the previous study contained data for only two oil fields with more than 30 well stimulations, the current data set contains four oil fields with more than 30 well stimulations, making these data more ideal for comparisons using statistical tests. More producers and a new oil field, Buena Vista Nose, are represented in the current data set. Some well stimulation practices remain relatively unchanged since the previous study. Well stimulation is mostly occurring in Kern County with the exception of a single well stimulation in Orange County. Almost all well stimulations are hydraulic fracturing. The current data set contains two acid fracturing treatments and one matrix acidizing treatment. The median number of chemicals added per well stimulation—excluding water and proppant—was previously 21; the number in the expanded data set is 20. Median water added per well stimulation was previously 89,000 gallons and is now 98,000 gallons. While this water use represents a 10% increase, water use for hydraulic fracturing in California is still lower than water use in other oil and gas fields. Chemical formulations appear to have shifted since the previous study. This shift may be related to a change in the predominant service company. We identified 26 chemicals for which the frequency of use has changed by more than 30%. The total number of chemicals used has also expanded. Previously, 178 unique chemicals were identified as being added to well stimulation fluids. The total number of unique chemicals is now 205. Analytical data available for base and recovered fluids has expanded. The number of base fluid samples has increased from 12 to 35 and the number of recovered fluid samples has increased from 1,078 to 2,166. The expanded data set is more ideal for an evaluation of indicator chemicals and comparisons across different formations. The analytical data for recovered fluid samples indicates that many ions, radioactive constituents, and organics are consistently observed, making these chemicals ideal for consideration as potential indicators of aquifer contamination.

54 ENVIRONMENTAL SCIENCES↗