Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “available energy”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.↗
NASA NTRS1983

Baroclinic instability in the stratosphere and mesosphere

The energy source of stratospheric planetary waves is largely due to vertical propagation from the troposphere, but the stratospheric zonal mean state which changes with time is not stable against small perturbations, whereby transient perturbation waves are generated at the expense of the stratospheric zonal mean available potential energy. This theoretical conclusion is obtained based on the quasi-geostrophic model in spherical coordinates which calculates the characteristic solutions for the basic zonal mean state changing with time.

Zhang, K. S.↗
NASA NTRS1976

Studies of the inner shelf and coastal sedimentation environment of the Beaufort Sea from ERTS-A

The author has identified the following significant results. Shearing periodically occurs between the westward moving pack ice (3 to 10 km/d) within the Pacific Gyre and the fast ice along the coast, forming major grounded shear and pressure ridges between the 10 to 40 m isobaths. Ridges occur in patterns conforming to known shoals. The zone of grounded ridges, called stamukhi zone, protects the inner shelf and coast from marine energy and pack ice forces. Relatively undeformed fast ice grows inshore of the stamukhi zone. The boundary is explained in terms of pack ice drift and major promontories and shoals. Intense ice gaging, highly disrupted sediments, and landward migration of shoals suggest that much of the available marine energy is expended on the sea floor within the stamukhi zone. Naleds (products of river icings) on the North Slope are more abundant east than west of the Colville River. Their location, growth, and decay were studied from LANDSAT imagery.

Reimnitz, E.↗
DOE OSTI2024

Cumulonimbus Clouds Convert a Smaller Fraction of CAPE into Kinetic Energy in a Warmer Atmosphere

Abstract This study investigates how entrainment’s diluting effect on cumulonimbus updraft buoyancy is affected by the temperature of the troposphere, which is expected to increase by the end of the century. A parcel model framework is constructed that allows for independent variations in the temperature ( T ), the entrainment rate ε , the free-tropospheric relative humidity (RH), and the convective available potential energy (CAPE). Using this framework, dilution of buoyancy is evaluated with T and RH independently varied and with CAPE either held constant or increased with temperature. When CAPE is held constant, buoyancy decreases as T increases, with parcels in warmer environments realizing substantially smaller fractions of their CAPE as kinetic energy (KE). This occurs because the increased moisture difference between an updraft and its surroundings at warmer temperatures drives greater updraft dilution. Similar results are found in midlatitude and tropical conditions when CAPE is increased with temperature. With the expected 6%–7% increase in CAPE per kelvin of warming, KE only increases at 2%–4% K −1 in narrow updrafts but tracks more closely with CAPE at 4%–6% in wider updrafts. Interestingly, the rate of increase in the KE with T becomes larger than that of CAPE when the later quantity increases at more than 10% K −1 . These findings emphasize the importance of considering entrainment in studies of moist convection’s response to climate change, as the entrainment-driven dilution of buoyancy may partially counteract the influence of increases in CAPE on updraft intensity. Significance Statement Cumulonimbus clouds mix air with their surrounding environment through a process called entrainment, which controls how efficiently environmental energy is converted into upward speed in thunderstorm updrafts. Our research shows that warmer temperatures will exacerbate the moisture difference between cumulonimbus updrafts and their surroundings, leading to greater mixing and less efficient conversion of environmental energy into updraft speeds. This effect should be considered in future research that investigates how climate change will affect cumulonimbus clouds.

Peters, John M.↗
DOE OSTI2021

Measurement of $d^2 \sigma/d|\vec{q}|d E_{avail}$ and 2p2h contribution using charged current $\nu_\mu$ interactions in the NOvA Near Detector

This Thesis presents the analysis methods and measurement of the $d^2 \sigma/d|\vec{q}|d E_{avail}$-CC inclusive cross section that describes neutrino scattering in a predominantly hydrocarbon medium. Studies that are prerequisite for the analysis including, variable determinations, event selections, efficiency, purity determinations, and systematic uncertainty estimations, are summarized. Tests of the analysis are also detailed, and shown to reproduce the input distributions. The double-differential cross-section measurement provides the foundation for determination of the inclusive rate for the excess induced by 2p2h processes together with nuclear medium effects that are not described by the Fermi gas model. The methodology for the extraction of this rate, defined as the excess observed relative to the rate estimated for known single-nucleon interactions, is described in detail. The following new measurements are reported in this Thesis: (1) The $\nu_\mu-$CC inclusive double-differential cross section as a function of three-momentum transfer and available hadronic energy, for an average $E_\nu$ of 1.8 GeV, is obtained. The double-differential cross section value that obtains when the bin widths are taken to be the dimensions of the entire analysis domain, e.g. $0.2 \le |\vec{q}| \le 2.0 GeV/c$ and $0.0 \le E_{avail} \le 2.0$ GeV, is inclusive cross section is $(6.05 \pm 0.75)\times10^{-39}$ cm$^2$/GeV/GeV/$c$/nucleon. An event excess attributed to 2p2h-MEC processes is observed; most of the rate occurs in a contiguous phase space region $0.3 \le |\vec{q}| \le 1.0$ GeV$/c$ $0.0 \le E_{avail} \le 0.35$ GeV. This rate represents $12.0\pm 6.5\%$ of the observed CC inclusive cross section; its cross section ratio relative to CCQE scattering is estimated to be $44.3 \pm 23.9\%$.

Olson, Travis Grant↗
NASA NTRS1981

Integrated and spectral energy flows of the GLAS GCM

Methods to analyze the generation, transport, and dissipation of energy to study geophysical fluid flows are discussed. Energetics analyses are pursued in several directions: (1) the longitudinal and time dependence on the energy flow to the stratosphere was examined as a function of geographical sector; (2) strong and weak energy flows were correlated by medium range forecasts; (3) the one dimensional spectral results (Fourier services around latitude circles) were extended to spherical harmonics over a global domain; (4) the validity of vertical velocities derived from mass convergence was examined for their effect on the conversion of eddy available potential energy to eddy kinetic energy.

Tennebaum, J.↗
NASA NTRS1986

The role of latent heat in kinetic energy conversions of South Pacific cyclones

The four-dimensional behavior of cyclone systems in the South Pacific Convergence Zone (SPCZ) is analyzed. Three cyclone systems, which occurred during the period from January 10-16, 1979, are examined using the data collected during the first special observing period of the FGGE. The effects of latent heating on the life cycles of the cyclones are investigated. Particular attention is given to the conversions of eddy available potential energy to eddy kinetic energy and of mean kinetic energy to eddy kinetic energy. The net radiation profile, sensible heat flux, total field of vertical motion, and latent heat component were computed. The life cycles of the cyclones are described. It is observed that the latent heating component accounts for nearly all the conversion in the three cyclones, and latent heating within the SPCZ is the major source of eddy kinetic energy for the cyclones.

Kann, Deirdre M.↗
DOE OSTI2022

Convective Parameters Derived from Radiosonde Data (SONDEPARAM) Value-Added Product Report

Radiosondes provide fundamental observations of the vertical profile of atmospheric state (pressure, temperature, humidity, and winds), with important implications for subsequent studies on environmental controls on cloud conditions. Within convective cloud environments, there is an increasing demand for additional value-added products (VAPs) to facilitate the use of U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility radiosonde data sets. Such VAPs should provide quick and reliable estimates for several standard radiosonde parameters or quantities of interest using common assumptions, as well as open, flexible code for visualization and user interaction. The Convective Parameters Derived from Radiosonde Data (SONDEPARAM) VAP will apply several robust algorithms used in Wang et al. (2020) for the calculation of useful radiosonde convective cloud parameters, including the convective available potential energy (CAPE), convective inhibition (CIN), and other convective parameters, for several different assumptions regarding the initial parcel characteristics (i.e., surface-based, most unstable, mixed layer). These ARM VAP codes are developed in open, flexible Python formats, with the intention that these parameters/calculations will be incorporated into traditional ARM quick-look radiosonde plotting, yet associated with user-available codes for ease in user reproduction and assumption modification.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

[Role of membrane lipids in myocardial cytoprotection]

The cardiomyocyte capacity to regulate ATP production to face any change in energy demand is a major determinant of cardiac function. This process is based on a balanced fatty acid (FA) metabolism, because FA is the main fuel of the heart, although the most expensive one in oxygen. The pathway is, however, weakly controlled by the cardiac myocyte which can well regulate FA mitochondrial entry but not cell FA uptake. For this reason, several pathological situations often result from either harmful accumulation of FA and derivatives or excess FA-oxidation. Control of the FA/glucose balance by decreased energy production from FA would thus offer an alternative strategy in the treatment of ischaemia, providing the cardiomyocytes weak ability in handling the non-metabolised FA is controlled. The initiation and the regulation of cardiac contraction both result from membrane activity; the other major role of lipids in the heart is their contribution to membrane homeostasis through phospholipid synthesis pathways and phospholipases. The anti-anginal activity of Trimetazidine, reported as a cytoprotective effect without a haemo-dynamic component; is associated with reduced use of FA for energy. However, accumulation of FA and derivatives has never been observed. Trimetazidine is reported to increase significantly the synthesis of phospholipids without influencing the other lipid classes, thus increasing the incorporation of FA in membrane structures. This cytoprotection appears to be based on the redirection of the use of FA to phospholipid synthesis, which would decrease their availability for energy production. This class of compounds, with the same properties as Trimetazidine, offers a metabolic approach to the treatment of ischaemia.

Cytoprotection/physiology↗
DOE OSTI2024

Hydropower capacity factors trending down in the United States

Abstract The United States hydropower fleet has faced increasing environmental and regulatory pressures over the last half century, potentially constraining total generation. Here we show that annual capacity factor has declined at four fifths of United States hydropower plants since 1980, with two thirds of decreasing trends significant at p < 0.05. Results are based on an analysis of annual energy generation totals and nameplate capacities for 610 plants (>5 megawatt), representing 87% of total conventional hydropower capacity in the United States. On aggregate, changes in capacity factor imply a fleetwide, cumulative generation decrease of 23% since 1980 before factoring in capacity upgrades—akin to retiring a Hoover Dam once every two to three years. Changes in water availability explain energy decline in only 21% of plants, highlighting the importance of non-climatic drivers of generation, including deterioration of plant equipment as well as changes to dam operations in support of nonpower objectives.

13 HYDRO ENERGY↗
NASA NTRS1999

Land-Atmosphere Interactions: Successes, Problems and Prospects

After two decades of active research, a much better understanding of the broader role of biospheric processes on the local climate has emerged. A surface-albedo increase, particularly in desert border regions of the subtropics (as well as the deforested tropical regions), leads to a net surface energy deficit, which in turn leads to a relative sinking and reduced rainfall. On the other hand, studies of the influence of altered ratios of evapotranspiration and sensible fluxes, in situations where the net solar income is unchanged, show that evapotranspiration is a more desirable flux for increased precipitation and vitality of the biosphere. Besides providing water vapor and convective available potential energy (CAPE) to the lower troposphere, evapotranspiration helps in building larger CAPE before "turning on" the moist-convection. Larger CAPE in the lower troposphere enables convection to reach into the deeper atmosphere thereby heating the upper troposphere; indeed, moist-convection is also accompanied by the evaporation of falling precipitation that cools and moistens the lower atmosphere. While convective, as opposed to stratiform, precipitation reduces the fractional cloud cover; it also allows more solar radiation to reach the surface thereby invigorating surface fluxes. These, together with moist convection and associated downdrafts help to maintain the characteristic upper temperature limit(s) of the moist-land as well as oceanic regions. Regardless of the above understanding, several important problems continue to hinder the accurate simulation of a realistic land atmosphere interaction in a numerical model (both GCM and/or Meso-scale models). Among the unsolved problems are parameterization of sub-grid scale land processes that include small-scale variability of soil moisture, snow-cover and snow-physics, the biodiversity of the biosphere, orography, local drainage characteristics under natural conditions, and surface flow over the natural terrain. A well-known non-linear response of surface fluxes to these variations makes the problem of parameterizing land-atmosphere interaction processes hard-to-address and simulate, particularly in a GCM. In our presentation, we will discuss how orographic, snow-cover, and water table interactions can be included into a Simple Biosphere Model such as SiB/SSiB. Figure I shows how, in the Russian region, spring snowmelt affects the soil moisture profile. Corresponding figure 2 shows how interaction with the water table decreases the natural evapotranspiration in the Sahel region simulation. While these simulations need better validation with data, the simulations reveal that surface processes are sensitive to these parameterizations. With these developments, we continue to advance our understanding of the interaction of land with the atmosphere aloft, but the intrinsic variability of the newer parameters, e. g., hydraulic properties of the soil, diminish the positive influences of these advances on the improved climate simulation with GCMs.

Sud, Y. C.↗
DOE OSTI2025

Laboratory Measurements of n >= 3 K-shell Transition Energies of Sulfur Ions from F-like S viii to Li-like S xiv

Inner-shell transitions are ubiquitous in nonequilibrium collisionally ionized plasmas, such as supernova remnants, and in photoionized plasmas, such as outflows from active galactic nuclei and X-ray binaries. Inner-shell X-ray emission can help determine key parameters of these systems, such as ionization time, τ, and ionization parameter, ξ. Despite their importance, only theoretical inner-shell transition energies are available for many ions. To provide experimental benchmarks, we have measured the dominant n → 1 K-shell transitions of sulfur ions where n ≥ 3 from Li-like S xiv to F-like S viii using LLNL’s SuperEBIT electron beam ion trap and the NASA/GSFC EBIT Calorimeter Spectrometer (ECS). We identify over 30 spectral features and measure their energies with uncertainties in the ∼0.1–1 eV range. We compare these results to Flexible Atomic Code (FAC) and multireference Møller–Plesset (MR-MP) calculations and find differences between theory and experiment of ∼1 eV for FAC and <0.5 eV for most MR-MP calculations. We also compare these results to two widely used atomic databases, AtomDB and CHIANTI, and find discrepancies as high as 7 eV. Furthermore, many transitions are missing from these databases despite being prominent in our data.

Atomic data benchmarking↗
NASA NTRS1984

SU(2) x U(1) vacuum and the Centauro events

It is proposed that the fireballs invoked to explain the Centauro events are bubbles of a metastable superdense state of nuclear matter, created in high energy (E is approximately 10 to the 15th power eV) cosmic ray collisions at the top of the atmosphere. If these bubbles are created with a Lorentz factor gamma approximately = 10 at their CM frame, the objections against the origin of these events in cosmic ray interactions are overcome. Assuming further, that the Centauro events are to the explosive decay of these metastable bubbles, a relationship between their lifetime, tau, and the threshold energy for bubble formation, E sub th, is derived. The minimum lifetime consistent with such an interpretation in tau is approximately 10 to the -8th power sec, while the E sub th appears to be insensitive to the value of tau and always close to E sub th is approximately 10 to the 15th power eV. Finally it is speculated that if the available CM energy is thermalized in such collisions, these bubbles might be manifestations of excitations of the SU(2) x U(1) false vacuum. The absence of neutral pions in the Centauro events is then explained by the decay of these excitations.

Kazanas, D.↗
NASA NTRS1984

Was the early Earth completely molten?

The nature of the interior of the primitive Earth was examined. The question is posted: was the Earth a cold solid or was it a hot liquid, much like a lava lakes seen in today's volcanic cauldrons. The various energy sources available to heat the primitive Earth to see if they are sufficient to cause melting were analyzed. The two largest contributors to the Earth's early heat appear to be the heat due to accretion and the heat of core formation. The Earth formed by the accretion of particles ranging in size from millimeters to hundreds of kilometers and each impact into the protoearth provided more energy to heat the body. It is found that early in the Earth's history, the sinking of iron to the center of the Earth to form the core released a substantial amount of energy; enough to heat the entire Earth an average 2000 deg C. Mechanisms for the removal of such a large amount of heat appear inadequate to prevent substantial melting, and it is assumed that the Earth was completely molten, i.e., a magma ocean at one time early in its history.

Turcotte, D. L.↗
DOE OSTI2023

Seasonal controls on isolated convective storm drafts, precipitation intensity, and life cycle as observed during GoAmazon2014/5

Abstract. Isolated deep convective cloud life cycle and seasonal changes in storm properties are observed for daytime events during the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Green Ocean Amazon Experiment (GoAmazon2014/5) campaign to understand controls on storm behavior. Storm life cycles are documented using surveillance radar from initiation through maturity and dissipation. Vertical air velocity estimates are obtained from radar wind profiler overpasses, with the storm environment informed by radiosondes. Dry-season storm conditions favored reduced morning shallow cloud coverage and larger low-level convective available potential energy (CAPE) than wet-season counterparts. The typical dry-season storm reached its peak intensity and size earlier in its life cycle compared with wet-season cells. These cells exhibited updrafts in core precipitation regions (Z>35 dBZ) to above the melting level as well as persistent downdrafts aloft within precipitation adjacent to their cores. Moreover, dry-season cells recorded more intense updrafts to earlier life cycle stages as well as a higher incidence of strong updrafts (i.e., >5 m s−1) at low levels. In contrast, wet-season storms were longer-lived and featured a higher incidence of moderate (i.e., 2–5 m s−1) updrafts aloft. These storms also favored a shift in their most intense properties to later life cycle stages. Strong downdrafts were less frequent within wet-season cells aloft, indicating a potential systematic difference in draft behaviors, as linked to graupel loading and other factors between the seasons. Results from a stochastic parcel model suggest that dry-season cells may expect stronger updrafts at low levels because of larger low-level CAPE in the dry season. Wet-season cells anticipate strong updrafts aloft because of larger free-tropospheric relative humidity and reduced entrainment-driven dilution. Enhanced dry-season downdrafts are partially attributed to increased evaporation, dry-air entrainment mixing, and negative buoyancy in regions adjacent to sampled dry-season cores.

54 ENVIRONMENTAL SCIENCES↗