Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice core”

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 37 records · Page 2

A Mountain Glacier Perspective on the Bipolar Seesaw

A global record of mountain glacier terminations during the last deglaciation (∼19–11 ka) dated by a large, uncurated data set of cosmogenic-nuclide exposure ages highlights a statistically significant asynchrony in termination ages between the Northern and Southern Hemispheres. This interhemispheric offset in the timing of glacier terminations is consistent with previously correlated ice core records that show a systematic interhemispheric lag in the timing of abrupt climate events, with the Southern Hemisphere leading the Northern Hemisphere by ∼300–3,500 years. Our analysis (a) aggregates cosmogenic-nuclide exposure ages from a global data set of moraines to discern climatically driven peaks in moraine emplacement events, and (b) utilizes a Monte Carlo simulation based on a null hypothesis that moraine emplacement is interhemispherically synchronous to estimate the statistical significance of the observed offset. The observed lag of Northern Hemisphere emplacement events compared to the Southern Hemisphere is statistically significant and is consistent with the “bipolar seesaw” pattern observed in ice core records.

58 GEOSCIENCES↗

CT Scans of Cores Metadata, Utqiagvik (Barrow), Alaska, 2015

Individual ice cores were collected from Barrow Environmental Observatory in Barrow, Alaska, throughout 2013 and 2014. Cores were drilled along different transects to sample polygonal features (i.e. the trough, center and rim of high, transitional and low center polygons). Most cores were drilled around 1 meter in depth and a few deep cores were drilled around 3 meters in depth. Three-dimensional images of the frozen cores were constructed using a medical X-ray computed tomography (CT) scanner. TIFF files can be uploaded to ImageJ (an open-source imaging software) to examine soil structure and soil densities within each core.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Hydraulic Conductivity Measurements, Utqiagvik (Barrow), Alaska, 2014

Six individual ice cores were collected from the Barrow Environmental Observatory in Barrow, Alaska, in May of 2013 as part of the Next Generation Ecosystem Experiment (NGEE). Each core was drilled at a different location to varying depths. After drilling, the cores were stored in coolers packed with dry ice and flown to Lawrence Berkeley National Laboratory (LBNL) in Berkeley, CA. 3-dimensional images of the cores were constructed using medical X-ray computed tomography (CT) scanner at 120kV. Hydraulic conductivity samples were extracted from these cores at LBNL Richmond Field Station in Richmond, CA, in February 2014 by cutting 5 to 8 inch segments using a chop saw. Samples were packed individually and stored at -20C freezing temperatures to minimize any changes in structure or loss of ice content prior to analysis. Hydraulic conductivity was determined through falling head tests using a permeameter [ELE International, Model #: K-770B] (Appendix A). Samples were placed in a latex membrane via a membrane stretcher while frozen. Use of a membrane stretcher made the membranes easier to secure and minimized contact with the sample. A clear polycarbonate sleeve, fabricated with a stainless steel ring at the bottom to keep the sleeve from floating, was placed around the sample inside the permeameter to minimize deformation during analysis. The permeameter was filled with water and 1.0 PSI of air was applied for confining pressure during sample defrost. Outflow valves were left open to allow for incremental thawing and samples were left to thaw for approximately 12 hours. After approximately 12 hours of thaw, initial falling head tests were performed. When the flow was significantly too fast or too slow, the analysis was stopped and the burette size was adjusted accordingly (i.e. a larger diameter burette was used for flows that were faster than desired or a smaller diameter burette was used for flows that were slower than desired). Two to four measurements were collected on each sample and collection stopped when the applied head load exceeded 25% change from the original load. Analyses were performed between 2 to 3 times for each sample. The final hydraulic conductivity calculations were computed using methodology of Das et al., 1985.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Deglaciation of northwestern Greenland during Marine Isotope Stage 11

Past interglacial climates with smaller ice sheets offer analogs for ice sheet response to future warming and contributions to sea level rise; however, well-dated geologic records from formerly ice-free areas are rare. Here we report that subglacial sediment from the Camp Century ice core preserves direct evidence that northwestern Greenland was ice free during the Marine Isotope Stage (MIS) 11 interglacial. Luminescence dating shows that sediment just beneath the ice sheet was deposited by flowing water in an ice-free environment 416 ± 38 thousand years ago. Provenance analyses and cosmogenic nuclide data and calculations suggest the sediment was reworked from local materials and exposed at the surface <16 thousand years before deposition. Ice sheet modeling indicates that ice-free conditions at Camp Century require at least 1.4 meters of sea level equivalent contribution from the Greenland Ice Sheet.

58 GEOSCIENCES↗

The Extraordinary March 2022 East Antarctica “Heat” Wave. Part II: Impacts on the Antarctic Ice Sheet

Abstract Between 15 and 19 March 2022, East Antarctica experienced an exceptional heat wave with widespread 30°–40°C temperature anomalies across the ice sheet. In Part I, we assessed the meteorological drivers that generated an intense atmospheric river (AR) that caused these record-shattering temperature anomalies. Here, we continue our large collaborative study by analyzing the widespread and diverse impacts driven by the AR landfall. These impacts included widespread rain and surface melt that was recorded along coastal areas, but this was outweighed by widespread high snowfall accumulations resulting in a largely positive surface mass balance contribution to the East Antarctic region. An analysis of the surface energy budget indicated that widespread downward longwave radiation anomalies caused by large cloud-liquid water contents along with some scattered solar radiation produced intense surface warming. Isotope measurements of the moisture were highly elevated, likely imprinting a strong signal for past climate reconstructions. The AR event attenuated cosmic ray measurements at Concordia, something previously never observed. Last, an extratropical cyclone west of the AR landfall likely triggered the final collapse of the critically unstable Conger Ice Shelf while further reducing an already record low sea ice extent. Significance Statement Using our diverse collective expertise, we explored the impacts from the March 2022 heat wave and atmospheric river across East Antarctica. One key takeaway is that the Antarctic cryosphere is highly sensitive to meteorological extremes originating from the midlatitudes and subtropics. Despite the large positive temperature anomalies driven from strong downward longwave radiation, this event led to huge amounts of snowfall across the Antarctic interior desert. The isotopes in this snow of warm airmass origin will likely be detectable in future ice cores and potentially distort past climate reconstructions. Even measurements of space activity were affected. Also, the swells generated from this storm helped to trigger the final collapse of an already critically unstable Conger Ice Shelf while further degrading sea ice coverage.

Wille, Jonathan D.↗

Ice nucleating particle concentrations in the Arctic during MOSAiC 2020

This data set contains ice nucleating particle (INP) measurements, using Colorado State University’s (CSU) Ice Spectrometer (IS) of samples collected during the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign. Samples were collected from five different media: aerosol filters, meltwater, ice cores, snow pits, and bulk sea water (BSW).

54 ENVIRONMENTAL SCIENCES↗

Global Ocean Cooling of 2.3°C During the Last Glacial Maximum

Abstract Quantitative constraints on past mean ocean temperature (MOT) critically inform our historical understanding of Earth's energy balance. A recently developed MOT proxy based on paleoatmospheric Xe, Kr, and N 2 ratios in ice core air bubbles is a promising tool rooted in the temperature dependences of gas solubilities. However, these inert gases are systematically undersaturated in the modern ocean interior, and it remains unclear how air‐sea disequilibrium may have changed in the past. Here, we carry out 30 tracer‐enabled model simulations under varying circulation, sea ice cover, and wind stress regimes to evaluate air‐sea disequilibrium in the Last Glacial Maximum (LGM) ocean. We find that undersaturation of all three gases was likely reduced, primarily due to strengthened high‐latitude winds, biasing reconstructed MOT by −0.38 ± 0.37°C (1σ). Accounting for air‐sea disequilibrium, paleoatmospheric inert gases indicate that LGM MOT was 2.27 ± 0.46°C (1σ) colder than the pre‐industrial era.

Seltzer, A. M.↗

Probing the angular and polarization reconstruction of the ARIANNA detector at the South Pole

The sources of ultra-high energy (UHE) cosmic rays, which can have energies up to 10 20 eV, remain a mystery. UHE neutrinos may provide important clues to understanding the nature of cosmic-ray sources. ARIANNA aims to detect UHE neutrinos via radio (Askaryan) emission from particle showers when a neutrino interacts with ice, which is an efficient method for neutrinos with energies between 10 16 eV and 10 20 eV. The ARIANNA radio detectors are located in Antarctic ice just beneath the surface. Neutrino observation requires that radio pulses propagate to the antennas at the surface with minimum distortion by the ice and firn medium. Using the residual hole from the South Pole Ice Core Project, radio pulses were emitted from a transmitter located up to 1.7 km below the snow surface. By measuring these signals with an ARIANNA surface station, the angular and polarization reconstruction abilities are quantified, which are required to measure the direction of the neutrino. After deconvolving the raw signals for the detector response and attenuation from propagation through the ice, the signal pulses show no significant distortion and agree with a reference measurement of the emitter made in an anechoic chamber. Furthermore, the signal pulses reveal no significant birefringence for our tested geometry of mostly vertical ice propagation. The origin of the transmitted radio pulse was measured with an angular resolution of 0.37° indicating that the neutrino direction can be determined with good precision if the polarization of the radio-pulse can be well determined. In this work, we obtained a resolution of the polarization vector of 2.7°. Neither measurement show a significant offset relative to expectation.

47 OTHER INSTRUMENTATION↗

Can We Better Constrain the Timing of GNAIW/UNADW Variability in the Western Equatorial Atlantic and Its Relationship to Climate Change During the Last Deglaciation?

We have revisited the well-trod VM28-122 core retrieved from the deep Colombian Basin, which includes sediments that reflect modern Upper North Atlantic Deep Water and extends through the last deglaciation into the last glacial period when the site was bathed in Glacial North Atlantic Intermediate Water. Here, we leverage the nearby Cariaco Basin's surface water radiocarbon reconstruction (reservoir age, and ΔR) on the IntCal20 timescale to recast the period of the last deglaciation with a newly constrained age model. Based on the revised age model, we observe that the multimillennial decrease in benthic δ 13 C and B/Ca (which record δ 13 C of dissolved inorganic carbon and Δ[CO 3 –2 ], respectively) began at 18,100 ± 240 (95% CI) calibrated years BP, synchronous with Termination 1, as identified by changes in the Antarctic Ice Sheet composite and by the onset of rapid glacier recession in the Southern Hemisphere (Denton et al., 2010, https://doi.org/10.1126/science.1184119; Denton et al., 2021, https://10.1016/j.quascirev.2020.106771). The beginning of the decrease in benthic δ 18 O is concurrent with the changes in carbon chemistry or at most, a few hundred years later. It is no later than 17,700 ± 300 (95% CI) yrs BP in our record, at the putative start of Heinrich Stadial 1. With sufficient data density (more than 2–3 control points per kyr) and an independent record of past surface water radiocarbon variations, it is possible to achieve late glacial and deglacial chronologies with fidelities similar to those of ice cores. Furthermore, doing so in more oceanographic locations should shed light more broadly on the mechanisms instrumental to abrupt climate change.

58 GEOSCIENCES↗

Lacustrine leaf wax hydrogen isotopes indicate strong regional climate feedbacks in Beringia since the last ice age

The Late-Quaternary climate of Beringia remains unresolved despite the region's role in modulating glacial-interglacial climate and as the likely conduit for human dispersal into the Americas. In this work, we investigate Beringian temperature change using an ~32,000-year lacustrine record of leaf wax hydrogen isotope ratios (δ2Hwax) from Arctic Alaska. Based on Monte Carlo iterations accounting for multiple sources of uncertainty, the reconstructed summertime temperatures were ~3 °C colder (range: –8 to +3 °C) during the Last Glacial Maximum (LGM; 21-25 ka) than the pre-industrial era (PI; 2–0.1 ka). This ice-age summer cooling is substantially smaller than in other parts of the Arctic, reflecting altered atmospheric circulation and increased continentality which weakened glacial cooling in the region. Deglacial warming was punctuated by abrupt events that are largely synchronous with events seen in Greenland ice cores that originate in the North Atlantic but which are also controlled locally, such as by the opening of the Bering Strait between 13.4 and 11 ka. Our reconstruction, together with climate modeling experiments, indicates that Beringia responds more strongly to North Atlantic freshwater forcing under modern-day, open-Bering Strait conditions than under glacial conditions. Furthermore, a 2 °C increase (Monte Carlo range: –1 to +5 °C) over the anthropogenic era reverses a 6 °C decline (Monte Carlo range: –10 to 0 °C) through the Holocene, indicating that recent warming in Arctic Alaska has not surpassed peak Holocene summer warmth.

54 ENVIRONMENTAL SCIENCES↗

Radiocarbon Fingerprinting Black Carbon Source History in the Himalayas

Abstract Black carbon (BC) is considered as an important contributor to the Himalayan glaciers melt in the past few decades. However, the long‐term source apportionment of BC remains unclear. Here we present the first radiocarbon ( 14 C)‐based annual variation of BC source apportionment in an ice core spanning the period of 1959–2012 drilled from the Southeastern Tibetan Plateau, a receptor site of South Asia outflow. We find fossil fuel combustion is a major contribution (73% ± 5%), yet the biomass burning fraction ( ƒ biomass ) has grown from 24% ± 4% to 30% ± 4% since 1990. Intriguingly, we further find the ƒ biomass demonstrating a robust correlation with South Asian wildfires linked to climate oscillations. Thus, for mitigating BC impacts on Himalayan glaciers, South Asia's transition from fossil fuels to clean energy is a more efficient and urgent strategy than reducing residential biomass burning.

Wang, Mo [State Key Laboratory of Tibetan Plateau ↗

Why Firn Quakes

Snow dampens sounds, but anecdotal reports concisely describe audible propagating collapse events—firnquakes—in Antarctic and Arctic snowfields. We propose combining granular and continuum mechanics to form a testable theory for conditioning, triggering, and propagation of firnquakes consistent with scarce data. A central condition for collapse events is unconsolidated firn at depth. As firn grains compact, stresses are transmitted along force chains which carry the overburden and transition into a continuous medium by pressure sintering. This granular legacy creates solid-like supports of denser layers that keep the material below unconsolidated. Dynamic amplification triggers local brittle failure of the supports, which induces a cascade of collapse propagation. Using bulk density from ice cores as proxy for stiffness, we find the flexural wave speed by collapsing supports matches the recorded firnquake velocities on the order of 100 m/s. Our theory is to be tested in firn sheets and other compacting granular materials.

Voigtländer, A. [Lawrence Berkeley National Labora↗

Stratospheric ozone depletion in the Antarctic region triggers intense changes in sea salt aerosol geochemistry

Abstract Since the early 1980s, the Antarctic environment has served as a natural field laboratory for researchers to investigate the effects of stratospheric ozone depletion, which has resulted in increased surface ultraviolet radiation levels. However, its effective threats still present gaps. We report new pieces of evidence of increased ultraviolet radiation impacting West Antarctica sea salt aerosols. Salt aerosols, particularly in the Southern Ocean Sea, play an important role in the radiative earth balance. To disclose the molecular details of sea salt aerosols, we used a synchrotron-based multi-element microscopic speciation of individual microparticles (Scanning Transmission X-ray Microscopy with Near-Edge X-ray Absorption Fine Structure Spectroscopy combined with Computer-Controlled Scanning Electron Microscopy). Here we identified substantial abundances of chlorine-enriched aerosols in sea salt generated by photolytic products, whereas ice core records revealed increased chlorine depletion from the onset of ozone depletion. Our findings reveal that modern sea salt modification has no Holocene precedent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating proxy influence in assimilated paleoclimate reconstructions – Testing the exchangeability of two ensembles of spatial processes

Climate field reconstructions (CFR) attempt to estimate spatiotemporal fields of climate variables in the past using climate proxies such as tree rings, ice cores, and corals. Data Assimilation (DA) methods are a recent and promising new means of deriving CFRs that optimally fuse climate proxies with climate model output. Despite the growing application of DA-based CFRs, little is understood about how much the assimilated proxies change the statistical properties of the climate model data. To address this question, we propose a robust and computationally efficient method, based on functional data depth, to evaluate differences in the distributions of two spatiotemporal processes. In this study, we apply our test to study global and regional proxy influence in DA-based CFRs by comparing the background and analysis states, which are treated as two samples of spatiotemporal fields. We find that the analysis states are significantly altered from the climate-model-based background states due to the assimilation of proxies. Moreover, the difference between the analysis and background states increases with the number of proxies, even in regions far beyond proxy collection sites. Our approach allows us to characterize the added value of proxies, indicating where and when the analysis states are distinct from the background states.

54 ENVIRONMENTAL SCIENCES↗

Water isotopes, climate variability, and the hydrological cycle: recent advances and new frontiers

Abstract The hydrologic cycle is a fundamental component of the climate system with critical societal and ecological relevance. Yet gaps persist in our understanding of water fluxes and their response to increased greenhouse gas forcing. The stable isotope ratios of oxygen and hydrogen in water provide a unique opportunity to evaluate hydrological processes and investigate their role in the variability of the climate system and its sensitivity to change. Water isotopes also form the basis of many paleoclimate proxies in a variety of archives, including ice cores, lake and marine sediments, corals, and speleothems. These records hold most of the available information about past hydrologic variability prior to instrumental observations. Water isotopes thus provide a ‘common currency’ that links paleoclimate archives to modern observations, allowing us to evaluate hydrologic processes and their effects on climate variability on a wide range of time and length scales. Building on previous literature summarizing advancements in water isotopic measurements and modeling and describe water isotopic applications for understanding hydrological processes, this topical review reflects on new insights about climate variability from isotopic studies. We highlight new work and opportunities to enhance our understanding and predictive skill and offer a set of recommendations to advance observational and model-based tools for climate research. Finally, we highlight opportunities to better constrain climate sensitivity and identify anthropogenically-driven hydrologic changes within the inherently noisy background of natural climate variability.

Dee, Sylvia (ORCID:000000022140785X)↗

The effects of rising CO 2 concentrations on terrestrial systems: scaling it up

Since the Industrial Revolution, atmospheric CO 2 concentrations have increased by c. 50%, from 280 ppm to a current level of 415 ppm and rising (Ciais et al., 2013). But CO 2 concentrations would be even higher if the terrestrial biosphere was not acting as a carbon sink, absorbing c. 30% of the CO 2 we emit every year (Le Quéré et al., 2016). Understanding how increasing CO 2 concentrations will alter the ability of vegetation and soils to sequester carbon is therefore critical for predicting the trajectory of future climate change, since a reduction in this carbon sink would cause a more rapid accumulation of CO2 in the atmosphere (Dusenge et al., 2019). In this issue of New Phytologist, Walker et al. (2021; pp. 2413–2445) synthesize data from an incredibly broad range of sources, including herbaria, free air CO 2 enrichment (FACE) studies, ice cores, eddy covariance sites, and remote sensing, and examine an enormous diversity of measurements (such as soil respiration rates, glucose isotopomers from leaves, tree ring width data, stream-gauges for runoff, and direct atmospheric CO 2 measurements) to address the question of how increasing CO2 concentrations are affecting the carbon uptake capacity of our planet.

54 ENVIRONMENTAL SCIENCES↗

Near-constant retreat rate of a terrestrial margin of the Laurentide Ice Sheet during the last deglaciation

The Laurentide Ice Sheet (LIS) was the largest ice sheet during the last glacial period. An accurate representation of its behavior during the last deglaciation is critical to understanding its influence on and response to a changing climate. Here we use 10 Be dating and Bayesian modeling to track the recession of the southwest sector of the Labrador Dome of the LIS along an ~500-km-long transect west of Lake Superior during the last deglaciation. This transect reflects terrestrial ice-margin retreat and crosses multiple moraine sets, with the southwestern part of the transect deglaciated by ca. 19 ka and the northeastern part deglaciated by ca. 10 ka. The predominant behavior of the ice margin during this interval is near-constant retreat with retreat rates varying between ~59 m/a and 38 m/a. The moraine sets mark standstills and/or readvances that in total constitute only ~17% of the retreat interval. The spatial and temporal pattern of ice-margin retreat tracked here differs from existing reconstructions that are based on using isochrons to define ice-margin positions. Acknowledging the uncertainties associated with the modeled ages of ice-margin retreat, we suggest that the overall retreat pattern is consistent with forcing by a gradual increase in Northern Hemisphere, high-latitude summer insolation. The pattern of ice-margin retreat is inconsistent with Greenland ice-core temperature records, and thus these records may not be suitable to drive models of the LIS.

58 GEOSCIENCES↗

Quantifying CanESM5 and EAMv1 sensitivities to Mt. Pinatubo volcanic forcing for the CMIP6 historical experiment

Large volcanic eruptions reaching the stratosphere have caused marked perturbations to the global climate including cooling at the Earth's surface, changes in large-scale circulation and precipitation patterns and marked temporary reductions in global ocean heat content. Many studies have investigated these effects using climate models; however, uncertainties remain in the modelled response to these eruptions. This is due in part to the diversity of forcing datasets that are used to prescribe the distribution of stratospheric aerosols resulting from these volcanic eruptions, as well as uncertainties in optical property derivations from these datasets. To improve this situation for the sixth phase of the Coupled Model Intercomparison Project (CMIP6), a two-step process was undertaken. First, a combined stratospheric aerosol dataset, the Global Space-based Stratospheric Aerosol Climatology (GloSSAC; 1979–2016), was constructed. Next, GloSSAC, along with information from ice cores and Sun photometers, was used to generate aerosol distributions, characteristics and optical properties to construct a more consistent stratospheric aerosol forcing dataset for models participating in CMIP6. This “version 3” of the stratospheric aerosol forcing has been endorsed for use in all contributing CMIP6 simulations. Recent updates to the underlying GloSSAC from version 1 to version 1.1 affected the 1991–1994 period and necessitated an update to the stratospheric aerosol forcing from version 3 to version 4. As version 3 remains the official CMIP6 input, quantification of the impact on radiative forcing and climate is both relevant and timely for interpreting results from experiments such as the CMIP6 historical simulations. This study uses two models, the Canadian Earth System Model version 5 (CanESM5) and the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1), to estimate the difference in instantaneous radiative forcing in simulated post-Pinatubo climate response when using version 4 instead of version 3. Differences in temperature, precipitation and radiative forcings are generally found to be small compared to internal variability. An exception to this is differences in monthly temperature anomalies near 24 km altitude in the tropics, which can be as large as 3 °C following the eruption of Mt. Pinatubo.

58 GEOSCIENCES↗