Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cross correlation”

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 145 records · Page 8

Baryon Acoustic Oscillations from the C IV Forest with DESI DR2

We present a measurement of Baryon Acoustic Oscillations (BAO) in the cross-correlation of triply ionized carbon C IV absorption with the positions of quasars (QSO) and Emission Line Galaxies (ELG). We use quasars and ELGs from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our data sample consists of 2.5 million quasars, 3.1 million ELGs, and the C IV absorption is measured along the line of sight of 1.5 million high redshift quasars with $z > 1.3$. We measure the isotropic BAO signal at 4.2$σ$ for the CIV$\times$QSO cross-correlation. This translates into a 3.0% precision measurement of the ratio of the isotropic distance scale, $D_{\rm V}$, and the sound horizon at the drag epoch, $r_{\rm d}$, with $D_{\rm V}/r_{\rm d}(z_{\rm eff} = 1.92) = 30.3 \pm 0.9$. We make the first detection of the BAO feature in the CIV$\times$ELG cross-correlation at a significance of 2.5$σ$ and find $D_{\rm V}/r_{\rm d}(z_{\rm eff} = 1.47) = 24.6 \pm 1.0$.

Bault, Abby [LBL, Berkeley] (ORCID:000000029964100↗

Measurements of cosmic expansion and growth rate of structure from voids in the Sloan Digital Sky Survey between redshift 0.07 and 1.0

ABSTRACT We present measurements of the anisotropic cross-correlation of galaxies and cosmic voids in data from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample, Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS luminous red galaxy catalogues from SDSS Data Releases 7, 12, and 16, covering the redshift range 0.07 < z < 1.0. As in our previous work analysing voids in subsets of these data, we use a reconstruction method applied to the galaxy data before void finding in order to remove selection biases when constructing the void samples. We report results of a joint fit to the multipole moments of the measured cross-correlation for the growth rate of structure, fσ8(z), and the ratio DM(z)/DH(z) of the comoving angular diameter distance to the Hubble distance, in six redshift bins. For DM/DH, we are able to achieve a significantly higher precision than that obtained from analyses of the baryon acoustic oscillations and galaxy clustering in the same data sets. Our growth rate measurements are of lower precision but still comparable with galaxy clustering results. For both quantities, the results agree well with the expectations for a Λ cold dark matter model. Assuming a flat Universe, our results correspond to a measurement of the matter density parameter $\Omega _\mathrm{m}=0.337^{+0.026}_{-0.029}$. For more general models, the degeneracy directions obtained are consistent with and complementary to those from other cosmological probes. These results consolidate void–galaxy cross-correlation measurements as a pillar of modern observational cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Angular correlations of cosmic microwave background spectrum distortions from photon diffusion

ABSTRACT During cosmic recombination, charged particles bind into neutral atoms and the mean free path of photons rapidly increases, resulting in the familiar diffusion damping of primordial radiation temperature variations. An additional effect is a small photon spectrum distortion, because photons arriving from a particular sky direction were originally in thermal equilibrium at various spatial locations with different temperatures; the combination of these different blackbody temperature distributions results in a spectrum with a Compton y-distortion. Using the approximation that photons had zero mean free path prior to their second-to-last scattering, we derive an expression for the resulting y-distortion, and compute the angular correlation function of the diffusion y-distortion and its cross-correlation with the square of the photon temperature fluctuation. Detection of the cross-correlation is within reach of existing arcminute-resolution microwave background experiments such as the Atacama Cosmology Telescope and the South Pole Telescope.

Astronomy & Astrophysics↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: cosmological implications from multitracer BAO analysis with galaxies and voids

We construct cosmic void catalogues with the dive void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy–galaxy, galaxy–void, and void–void correlation functions of all data sets, including the ones cross-correlating LRG and ELG samples. Two multitracer BAO fitting schemes are tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5–15 per cent, yielding 0.9 per cent, 0.8 per cent, 1.1 per cent, 2.3 per cent, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Lyα samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\, {\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$, Ω m = 0.290 ± 0.015, and $\Omega _\Lambda h^2 = 0.3241 \pm 0.0079$ in the flat-ΛCDM framework, where the 1σ uncertainties are around 6 per cent, 6 per cent, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG–ELG cross-correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Accretion disc sizes from continuum reverberation mapping of AGN selected from the ZTF survey

ABSTRACT We present the accretion disc-size estimates for a sample of 19 active galactic nuclei (AGNs) using the optical g-, r-, and i-band light curves obtained from the Zwicky Transient Facility survey. All the AGNs have reliable supermassive black hole (SMBH) mass estimates based on previous reverberation mapping measurements. The multiband light curves are cross-correlated, and the reverberation lag is estimated using the Interpolated Cross-Correlation Function method and the Bayesian method using the javelin code. As expected from the disc-reprocessing arguments, the g − r band lags are shorter than the g − i band lags for this sample. The interband lags for all, but five sources, are larger than the sizes predicted from the standard Shakura Sunyaev (SS) analytical model. We fit the light curves directly using a thin disc model implemented through the javelin code to get the accretion disc sizes. The disc sizes obtained using this model are on an average 3.9 times larger than the prediction based on the SS disc model. We find a weak correlation between the disc sizes and the known physical parameters, namely the luminosity and the SMBH mass. In the near future, a large sample of AGNs covering broader ranges of luminosity and SMBH mass from large photometric surveys would be helpful in a better understanding of the structure and physics of the accretion disc.

Jha, Vivek Kumar (ORCID:0000000232776335)↗

The Squeezed Bispectrum from CHIME H I Emission and Planck Cosmic Microwave Background Lensing: Current Sensitivity and Forecasts

Line intensity mapping using atomic hydrogen (H I) has the potential to efficiently map large volumes of the Universe if the signal can be successfully separated from overwhelmingly bright radio foreground emission. This motivates cross correlations, to ascertain the cosmological nature of measured H I fluctuations, and to study their connections with galaxies and the underlying matter density field. However, these same foregrounds render the cross correlation with projected fields such as the lensing of the cosmic microwave background (CMB) difficult. Indeed, the correlated Fourier modes vary slowly along the line of sight and are thus most contaminated by the smooth-spectrum radio continuum foregrounds. In this paper, we implement a method that avoids this issue by attempting to measure the nonlinear gravitational coupling of the small-scale 21 cm power from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) with large-scale Planck CMB lensing. This measurement is a position-dependent power spectrum, i.e., a squeezed integrated bispectrum. Using 94 nights of CHIME data between 1.0 < z < 1.3 and aggressive foreground filtering, we find that the expected signal is 5 times smaller than the current noise. We forecast that incorporating the additional nights of CHIME data already collected would enable a signal-to-noise ratio of 3, without any further improvements in filtering for foreground cleaning.

21cm↗

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization Universe as a means of measuring large-scale cosmic structure. Using nine nights of observations targeting the North Celestial Pole field, totaling approximately 107 hr of integration time, we analyze data in the frequency range 700–800 MHz (corresponding to redshift z ∼ 0.9). We do the data format conversion, radio frequency interference flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition. The spherically averaged power spectrum Δ 2 (k) is obtained. Furthermore, this work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

cosmology: large-scale structure of universe↗

New full-sky studies of the distribution of ultra-high-energy cosmic-ray arrival directions

Ground-based full-sky studies of the angular distribution of arrival directions of ultra-high-energy cosmic rays require combining data from different observatories, such as the Pierre Auger Observatory (Auger) and the Telescope Array (TA), because no single array can cover all declinations. A working group comprising members from the Auger and TA collaborations has been tasked with performing such studies for more than a decade and has found several indications of full-sky anisotropies. Here, we update the results for the large- and medium-scale anisotropy analyses using the latest data from TA, which include corrections for daily and yearly atmospheric effects in data for large-scale anisotropies and looser selection criteria in data for medium-scale anisotropies. We extend the latter one by considering two more galaxy catalogues, consisting of jetted or all AGNs. Finally we also introduce a new angular harmonic space analysis that allows us to measure both the auto-correlation and cross-correlation with all catalogues for all multipoles independently ($\ell_\text{max} = 20$ in this work) and scanning the energy threshold.

Abdul Halim, A. [Adelaide U.]↗

Wide-angle effects for peculiar velocities

ABSTRACT The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as ‘wide-angle’ effects in redshift-space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work, we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in autocorrelation and in cross-correlation with the density field. Using both linear theory and the Zel'dovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low-redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration space, Fourier space, and spherical analyses; providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

79 ASTRONOMY AND ASTROPHYSICS↗

3D correlations in the Lyman-α forest from early DESI data

We present the first measurements of Lyman-α (Lyα) forest correlations using early data from the Dark Energy Spectroscopic Instrument (DESI). We measure the auto-correlation of Lyα absorption using 88 509 quasars at z > 2, and its cross-correlation with quasars using a further 147 899 tracer quasars at z ≳ 1.77. Then, we fit these correlations using a 13-parameter model based on linear perturbation theory and find that it provides a good description of the data across a broad range of scales. We detect the BAO peak with a signal-to-noise ratio of 3.8σ, and show that our measurements of the auto- and cross-correlations are fully-consistent with previous measurements by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). Even though we only use here a small fraction of the final DESI dataset, our uncertainties are only a factor of 1.7 larger than those from the final eBOSS measurement. As a result, we validate the existing analysis methods of Lyα correlations in preparation for making a robust measurement of the BAO scale with the first year of DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗

Correcting correlation functions for redshift-dependent interloper contamination

ABSTRACT The construction of catalogues of a particular type of galaxy can be complicated by interlopers contaminating the sample. In spectroscopic galaxy surveys this can be due to the misclassification of an emission line; for example in the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift [O ii] emitters may make up a few per cent of the observed Ly α emitter (LAE) sample. The presence of contaminants affects the measured correlation functions and power spectra. Previous attempts to deal with this using the cross-correlation function have assumed sources at a fixed redshift, or not modelled evolution within the adopted redshift bins. However, in spectroscopic surveys like HETDEX, where the contamination fraction is likely to be redshift dependent, the observed clustering of misclassified sources will appear to evolve strongly due to projection effects, even if their true clustering does not. We present a practical method for accounting for the presence of contaminants with redshift-dependent contamination fractions and projected clustering. We show using mock catalogues that our method, unlike existing approaches, yields unbiased clustering measurements from the upcoming HETDEX survey in scenarios with redshift-dependent contamination fractions within the redshift bins used. We show our method returns autocorrelation functions with systematic biases much smaller than the statistical noise for samples with at least as high as 7 per cent contamination. We also present and test a method for fitting for the redshift-dependent interloper fraction using the LAE–[O ii] galaxy cross-correlation function, which gives less biased results than assuming a single interloper fraction for the whole sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Atacama Cosmology Telescope: Multiprobe cosmology with unWISE galaxies and ACT DR6 CMB lensing

We present a joint analysis of the cosmic microwave background (CMB) lensing power spectra measured from the Data Release 6 of the Atacama Cosmology Telescope (ACT) and Planck PR4, cross-correlations between the ACT and Planck lensing reconstruction and galaxy clustering from unWISE, and the unWISE clustering auto-spectrum. We obtain 1.5% constraints on the matter density fluctuations at late times parametrized by the best constrained parameter combination 𝑆$^{3⁢x⁢2⁢pt}_{8}$ ≡ 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.4 = 0.815 ± 0.012. The commonly used 𝑆 8 ≡ 𝜎 8⁢ (Ω 𝑚 /0.3) 0.5 parameter is constrained to 𝑆 8 = 0.816 ± 0.015. In combination with baryon acoustic oscillation (BAO) measurements we find 𝜎 8 = 0.815 ± 0.012. We also present sound-horizon-independent estimates of the present day Hubble rate of 𝐻 0 = 66.4$^{+3.2}_{−3.7}$ km s −1 Mpc −1 from our large scale structure data alone and 𝐻 0 = 64.3$^{+2.1}_{−2.4}$ km s −1 Mpc −1 in combination with uncalibrated supernovae from Pantheon+. Using parametric estimates of the evolution of matter density fluctuations, we place constraints on cosmic structure in a range of high redshifts typically inaccessible with cross-correlation analyses. Combining lensing cross- and autocorrelations, we derive a 3.3% constraint on the integrated matter density fluctuations above 𝑧 = 2.4, one of the tightest constraints in this redshift range and fully consistent with a Λ cold dark matter (Λ⁢CDM) model fit to the primary CMB from Planck. Finally, combining with primary CMB observations and using the extended low redshift coverage of these combined datasets we derive constraints on a variety of extensions to the Λ⁢CDM model including massive neutrinos, spatial curvature, and dark energy. We find in flat Λ⁢CDM⁢ ∑𝑚 𝜈 < 0.12 eV at 95% confidence using the large scale structure data, BAO measurements from Sloan Digital Sky Survey, and primary CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Graphical Optimization of Spectral Shift Reconstructions for Optical Backscatter Reflectometry

Optical backscatter reflectometry (OBR) is an interferometric technique that can be used to measure local changes in temperature and mechanical strain based on spectral analyses of backscattered light from a singlemode optical fiber. The technique uses Fourier analyses to resolve spectra resulting from reflections occurring over a discrete region along the fiber. These spectra are cross-correlated with reference spectra to calculate the relative spectral shifts between measurements. The maximum of the cross-correlated spectra—termed quality—is a metric that quantifies the degree of correlation between the two measurements. Recently, this quality metric was incorporated into an adaptive algorithm to (1) selectively vary the reference measurement until the quality exceeds a predefined threshold and (2) calculate incremental spectral shifts that can be summed to determine the spectral shift relative to the initial reference. Using a graphical (network) framework, this effort demonstrated the optimal reconstruction of distributed OBR measurements for all sensing locations using a maximum spanning tree (MST). By allowing the reference to vary as a function of both time and sensing location, the MST and other adaptive algorithms could resolve spectral shifts at some locations, even if others can no longer be resolved.

47 OTHER INSTRUMENTATION↗

Mallat Scattering Transformation based surrogate for Magnetohydrodynamics

Abstract A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x , y , and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

97 MATHEMATICS AND COMPUTING↗

Lower bias, lower noise CMB lensing with foreground-hardened estimators

Extragalactic foregrounds in temperature maps of the cosmic microwave background (CMB) severely limit the ability of standard estimators to reconstruct the weak lensing potential. These foregrounds are not fully removable by multifrequency cleaning or masking and can lead to large biases if not properly accounted for. For foregrounds made of a number of unclustered point sources, an estimator for the source amplitude can be derived and deprojected, removing any bias to the lensing reconstruction. We show with simulations that all of the extragalactic foregrounds in temperature can be approximated by a collection of sources with identical profiles, and that a simple bias hardening technique is effective at reducing any bias to lensing, at a minimal noise cost. In this work, we compare the performance and bias to other methods such as "shear-only" reconstruction, and discuss how to jointly deproject any arbitrary number of foregrounds, each with an arbitrary profile. In particular, for a Simons Observatory-like experiment foreground-hardened estimators allow us to extend the maximum multipole used in the reconstruction, increasing the overall statistical power by ~50% over the standard quadratic estimator, both in auto and cross-correlation. We conclude that source hardening outperforms the standard lensing quadratic estimator both in auto- A nd cross-correlation, and in terms of lensing signal-to-noise and foreground bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DESI DR2 Baryon Acoustic Oscillations from the Lyman Alpha Forest Multipoles

We present an alternative measurement of the Baryon Acoustic Oscillation (BAO) using the Legendre multipole representation of the Ly$α$ forest correlation functions from the second data release (DR2) of the Dark Energy Spectroscopic Instrument survey. Compressing the auto- and cross-correlation functions into Legendre multipoles yields a positive-definite covariance matrix without any smoothing -- unlike the baseline DR2 analysis -- thanks to a significantly reduced data vector size. We introduce the statistical corrections required to debias the finite-sample covariance matrix estimate and demonstrate that monopole and quadrupole terms for both auto- and cross-correlations can be used even when the correlation functions are distorted by continuum errors and contaminated by metals. This formalism has slightly diminished the constraining power of the BAO scale, while considerably weakening constraints on nuisance parameters. We measure the isotropic BAO scale with $0.93\%$ precision at $z_\mathrm{eff}=2.35$, the Hubble parameter $H(z_\mathrm{eff})=(239.5\pm3.4)~(147.09~\mathrm{Mpc}/r_d) ~\mathrm{km~s}^{-1}~\text{Mpc}^{-1}$, and the transverse comoving distance $D_M(z_\mathrm{eff})=(5.80 \pm 0.10)~(r_d/147.09~\mathrm{Mpc})$~Gpc for a given value of the sound horizon ($r_d$). Our BAO results are entirely consistent with the baseline DR2 analysis.

Karaçaylı, Naim Göksel [Chicago U., KICP; Ohio Sta↗

Unraveling emission line galaxy conformity at z ∼ 1 with DESI early data

Emission line galaxies (ELGs) are now the preeminent tracers of large-scale structure at z > 0.8 due to their high density and strong emission lines, which enable accurate redshift measurements. However, relatively little is known about ELG evolution and the ELG–halo connection, exposing us to potential modelling systematics in cosmology inference using these sources. In this paper, we use a variety of observations and simulated galaxy models to propose a physical picture of ELGs and improve ELG–halo connection modelling in a halo occupation distribution framework. We investigate Dark Energy Spectroscopic Instrument (DESI)-selected ELGs in COSMOS data, and infer that ELGs are rapidly star-forming galaxies with a large fraction exhibiting disturbed morphology, implying that many of them are likely to be merger-driven starbursts. We further postulate that the tidal interactions from mergers lead to correlated star formation in central–satellite ELG pairs, a phenomenon dubbed ‘conformity’. We argue for the need to include conformity in the ELG–halo connection using galaxy models such as IllustrisTNG, and by combining observations such as the DESI ELG autocorrelation, ELG cross-correlation with luminous red galaxies, and ELG–cluster cross-correlation. We also explore the origin of conformity using the UniverseMachine model and elucidate the difference between conformity and the well-known galaxy assembly bias effect.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗