Search NASA⌕ Search

SEARCH · Search NASA

Results for “pressurized balloons”

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 73 records · Page 4

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identifying microstructures susceptible to pulverization in commercially irradiated high burnup UO2 under LOCA conditions

High burnup fuel fragmentation (HBFF) has been a concern in the nuclear industry for many years. When UO2 reaches and exceeds pellet average burnups around 55 GWd/tU, microstructural changes in the fuel pellet begin to occur that render the pellet susceptible to fine fragmentation under loss-of-coolant accident (LOCA) conditions. During a LOCA event, the cladding balloons and bursts, potentially releasing part of these fine fuel fragments into the reactor pressure vessel. This process is known as fuel fragmentation, relocation, and dispersal (FFRD). Presently, the US nuclear industry is developing a safety basis for FFRD with the goal of increasing burnups and pressurized water reactor cycle lengths. Increased cycle lengths would push the burnup limits of the fuel rods past the known threshold for HBFF susceptibility. In this work, the microstructural impact on HBFF was investigated by comparing the microstructures of as-irradiated rod segments to their post-LOCA tested counterparts. This investigation found that varied operational histories result in different microstructural evolutions across the pellet radius, which impacts the fragmentation behavior and radial location of the fragmentation. It is noted that regions with a high bubble density and a high density of grain boundaries fragment under LOCA-relevant conditions. In addition to the fragmentation that has been previously reported at the periphery of the fuel, fragmentation was also noted in the dark zone toward the fuel center. Analysis of the power histories of the fuel samples suggests that the dark zone forms in the central regions of the pellet between temperatures of approximately 740 and 960 °C.

McKinney, Casey [ORNL] (ORCID:0000000335383614)↗

Type-I ELM mitigation by continuous lithium granule gravitational injection into the upper tungsten divertor in EAST

Large edge-localized modes (ELMs) were mitigated by gravitational injection of lithium granules into the upper X-point region of the experimental advanced superconducting tokamak (EAST) device with tungsten plasma-facing components. The maximum ELM size was reduced by ~70% in high β N H-mode plasmas. Large ELM stabilization was sustained for up to about 40 energy confinement times, with constant core radiated power and no evidence of high-Z or low-Z impurity accumulation. The lithium granules injection reduced the edge plasma pedestal density and temperature and their gradients, due to increased edge radiation and reduced recycling from the plasma-facing components. Ideal stability calculations using the ELITE code indicate that the stabilization of large ELMs correlates with improved stability of intermediate-n peeling-ballooning modes, due to reduced edge current resulting from the profile changes. The pedestal pressure reduction was partially offset by a core density increase, which resulted in a modest ~7% drop in core stored energy and normalized energy confinement time. We surmise that the remnant small ELMs are triggered by the penetration of multiple Li granules just past the separatrix, similar to small ELMs triggered by deuterium pellet Futatani et al (2014 Nucl. Fusion 54 073008). This study extends previous ELM elimination with Li powder injection Maingi et al (2018 Nucl. Fusion 58 024003) in EAST because (1) use of small, dust-like powder and the related potential health hazards were eliminated, and (2) use of macroscopic granules should be more applicable to future devices, due to deeper penetration than dust particles, e.g. inside the separatrix with velocities ~10 m s -1 in EAST.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Type-I ELM mitigation by continuous lithium granule gravitational injection into the upper tungsten divertor in EAST

Large edge-localized modes (ELMs) were mitigated by gravitational injection of lithium granules into the upper X-point region of the EAST device with tungsten plasma-facing components. The maximum ELM size was reduced by ~ 70% in high βN H-mode plasmas. Large ELM stabilization was sustained for up to about 40 energy confinement times, with constant core radiated power and no evidence of high-Z or low-Z impurity accumulation. The lithium granules injection reduced the edge plasma pedestal density and temperature and their gradients, due to increased edge radiation and reduced recycling from the plasma-facing components. Ideal stability calculations using the ELITE code indicate that the stabilization of large ELMs correlates with improved stability of intermediate-n peeling-ballooning modes, due to reduced edge current resulting from the profile changes. The pedestal pressure reduction was partially offset by a core density increase, which resulted in a modest ~ 7% drop in core stored energy and normalized energy confinement time. We surmise that the remnant small ELMs are triggered by the penetration of multiple Li granules just past the separatrix, similar to small ELMs triggered by deuterium pellet [S. Futatani et al., Nucl. Fusion 54 (2014) 073008]. This study extends previous ELM elimination with Li powder injection [R. Maingi et al., Nucl. Fusion 58 (2018) 024003] in EAST because 1) use of small, dust-like powder and the related potential health hazards were eliminated, and 2) use of macroscopic granules should be more applicable to future devices, due to deeper penetration than dust particles, e.g. inside the separatrix with velocities ~ 10 m/s in EAST.

EAST↗

Tethered Balloon System Ozone Profiles during CoURAGE Summer Intensive Operational Period Field Campaign Report

During the summer IOP, the small ARM field campaign CRGTBSO3 collected measurements to gain an improved understanding of differences between the atmospheric composition in and just above the marine layer at the CoURAGE TBS site near the eastern shore of Chesapeake Bay. CRGTBSO3 included guest instrumentation with ozone (O 3 ) profile measurements on the TBS and surface O3 and meteorological measurements at the TBS site. An En-Sci 2Z electrochemical cell (ECC) ozonesonde (Komhyr 1969, 1986, Witte et al. 2018) was included on the TBS. The ozonesonde was connected to an InterMet iMet-4RSB radiosonde, and the overall data collected included vertical profiles of ozone, relative humidity, temperature, pressure, and altitude. The CRGTBSO3 iMet-4 radiosonde data is identical to that of the iMet in the Tethered Balloon System Merged Data Product (TBSMERGED; Gaustad and Dexheimer 2025). The TBSMERGED data product also includes meteorological data from a different sensor, the iMet XQ2. In some cases, the iMet XQ2 relative humidity (RH) data may be more accurate than the iMet-4, such as for some instances when the iMet-4 RH data stays at 100% for an extended period throughout a profile.

54 ENVIRONMENTAL SCIENCES↗

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗

Observation of a new pedestal stability regime in MAST Upgrade H-mode plasmas

Abstract The first pedestal stability and structure analysis on the new MAST Upgrade (MAST-U) spherical tokamak H-mode plasmas is presented. Our results indicate that MAST-U pedestals are close to the low toroidal mode number ( n ) peeling branch of the peeling-ballooning instability, in contrast with MAST H-mode pedestals which were deeply in the high- n ballooning branch. This offers the possibility of reaching the ELM-free quiescent H-mode (Burrell et al 2005 Plasma Phys. Control. Fusion 47 B37–B52) or high-performance super H-mode (Snyder et al 2015 Nucl. Fusion 55 083026; Snyder et al 2019 Nucl. Fusion 59 086017) regimes. In addition, the coupling between the peeling and ballooning branches is weak in MAST-U, suggesting that a path to very high pedestal pressure gradient at high density may exist with sufficient heating power. A possible explanation for the differences between MAST and MAST-U pedestal stability is given in terms of plasma shaping parameters, in particular squareness and elongation, as well as the pedestal top temperature and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fuel performance analysis of Cr-coated Zircaloy-4 cladding during a prototypical LOCA event using BISON

Here, deformation and failure of chromium (Cr) coated Zircaloy-4 (Zry-4) were studied in loss-of-coolant accident (LOCA) conditions using the BISON fuel performance code. The BISON validation model simulating Halden research reactor experiments was extended to include Cr coatings and higher rod internal pressures to simulate high-burnup fuel. The transient simulations show Cr coatings help relieve stress in the Zry-4 substrate during the transient, delaying the onset of high-temperature creep, which leads to ballooning and bursting of the cladding. A nominal Cr coating thickness of 30 µm delays clad failure by 26 seconds and increases the clad burst temperature by 40 K. A parametric study showed that time to failure and burst temperature both increase with coating thickness, and Cr-coated cladding offers burst resistance under a wide range of rod internal pressures simulating high-burnup fuel. Results indicate that a thin Cr coating provides resistance against ballooning and bursting of cladding during LOCA events.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

In-situ determination of strain during transient burst testing and the temperature dependence of Zircaloy-4 claddings

Understanding fuel system behavior during postulated loss-of-coolant accidents is pertinent for continued safe and efficient operation of light water reactors, particularly as higher burnups are being pursued and safety margins re-evaluated. Conventional mechanical models for the incumbent Zr alloys typically rely on the assumption that steady-state creep is the dominant fuel cladding response during transient accident conditions. To investigate this assumption, simulated accident burst testing was performed on Zircaloy-4 claddings with balloon behavior measured in-situ. Here, two distinct loading conditions were utilized during burst testing: (1) constant-gas-inventory where pressure was allowed to increase with temperature and (2) constant pressure. In-situ strains and strain rates were measured via 2-dimensional digital image correlation techniques and synchronized with temperature to determine deformation dependencies. The temperature dependence of strain rate was characterized by a two segment Arrhenius relationship, with a distinct transition between the high and low temperature/strain regimes. The average activation energy of the lower temperature/strain regime was 328 ± 25 kJ/mol, in agreement with the ~320 kJ/mol used for conventional LOCA models. However, the higher temperature/strain segment, which encompassed most of ballooning, showed increased activation energies as well as a dependence on whether the burst region was in view. For tests that burst away from the camera view, the average high temperature/strain segment activation energy was 635 ± 150 kJ/mol. For samples where the rupture opening formed in view, the average activation energy was 1015 ± 179 kJ/mol. This observed shift in temperature dependence indicates a transition in deformation mechanism at the end of life, possibly to time independent failure mechanisms, which has not yet been visualized in the literature for Zr alloys. Parameters at the transition points were analyzed to determine thresholds for this change in behavior, which occurred at an average hoop strain of 6.9 ± 2.1 %.

36 MATERIALS SCIENCE↗

Impact of LWR assembly structural features on cladding burst behavior under LOCA conditions

This provides an initial scoping study on clad balloon and burst behavior for burnup extension of reactor fuel. The associated issues with burnup extension are fuel fragmentation, relocation, and dispersal in the event of cladding failure. The general finding of this work is that the structural features, spacer grids and mixing vanes, locally suppress cladding deformation but have little impact on the overall clad performance during loss-of-coolant accidents. The work detailed in a previous report by Capps et al. focused on core optimization via neutronics, thermal hydraulics and thermomechanical analysis for burnups beyond 62 GWD/tU and enrichments above 5%. Uncertainty of fuel fragmentation relocation and dispersal in high-burnup rods during accident conditions was also investigated. The dispersal aspect of fuel fragmentation depends on cladding rupture. Thus, assessing uncertainties in the rupture behavior is helpful in estimating the dispersal of high-burnup fuel. This study builds on the previous work by assessing the impact of assembly structural features on cladding balloon and burst behavior in a full-length fuel rod. In this work, the BISON fuel performance code was used to generate 2D radial and height meshes containing structural features commonly used in nuclear fuel assemblies. First, meshes were generated with spacer grids. Results were then compared to the cladding burst temperature and balloon strain results from the previous work. A mesh sensitivity study was performed to ensure that mixing vanes and spacer grid effects were appropriately considered, resulting in a more refined mesh than the previous study. The balloon deformation and burst times of the cladding were compared to the original case. Consideration was also given to the effect of rod initial pressure. In conclusion, 3D quarter rod simulations were also performed and found good agreement with the 2D simulations in clad deformation and reasonable agreement in burst times.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hydrogen Effects on Cladding Performance during Simulated Loss of Coolant Accident Testing

Fresh Zircaloy cladding segments were charged with hydrogen and subjected to simulated Loss of Coolant Accident (LOCA) conditions to assess the impact of hydrogen on the cladding mechanical performance and rupture resistance. The cladding rupture temperature decreased by 10 °C for every 100 wppm H added into the system for two different internal pressurization cases. Post-test analysis showed a loss of ductility with increasing hydrogen content, evident through decreased balloon sizes, smaller rupture openings, and sharper rupture interfaces. Digital image correlation was used to map the cladding strain in situ for two cladding segments. An elevated intermediate temperature uniform strain was visible for the hydride cladding, occurring prior to the decreased ductility during balloon and rupture determined from posttest fractography and microscopy. The thermal stability of Zr alloys as a function of hydrogen concentration and the thermal stability of Zr hydride as a function of temperature are likely responsible for the uniform deformation behavior at intermediate temperatures. The present framework for rapid assessment of hydrogen effects on cladding rupture can be utilized as an initial comparison to high burnup cladding thermomechanical response to a LOCA, where local regions may contain high concentrations of precipitated Zr hydride phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hydrogen Effects on Cladding Performance during Simulated Loss of Coolant Accident Testing

Fresh Zircaloy cladding segments were charged with hydrogen and subjected to simulated Loss of Coolant Accident (LOCA) conditions to assess the impact of hydrogen on the cladding mechanical performance and rupture resistance. The cladding rupture temperature decreased by 10 °C for every 100 wppm H added into the system for two different internal pressurization cases. Post-test analysis showed a loss of ductility with increasing hydrogen content, evident through decreased balloon sizes, smaller rupture openings, and sharper rupture interfaces. Digital image correlation was used to map the cladding strain in situ for two cladding segments. An elevated intermediate temperature uniform strain was visible for the hydride cladding, occurring prior to the decreased ductility during balloon and rupture determined from posttest fractography and microscopy. The thermal stability of Zr alloys as a function of hydrogen concentration and the thermal stability of Zr hydride as a function of temperature are likely responsible for the uniform deformation behavior at intermediate temperatures. The present framework for rapid assessment of hydrogen effects on cladding rupture can be utilized as an initial comparison to high burnup cladding thermomechanical response to a LOCA, where local regions may contain high concentrations of precipitated Zr hydride phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Site G - Surface Meteorological Station / Reviewed Data

This dataset contains meteorological data both from the ground and tether termination near the balloon on tethered balloon system at AWAKEN site G. The met data include time stamp, wind speed, temperature, relative humidity, and pressure.

17 WIND ENERGY↗

Site A1 - Surface Meteorological Station / Reviewed Data

This dataset contains meteorological data both from the ground and tether termination near the balloon on tethered balloon system at AWAKEN site A1. The met data include time stamp, wind speed, temperature, relative humidity, and pressure.

17 WIND ENERGY↗

Access to stable, high pressure tokamak pedestals using local electron cyclotron current drive

Abstract Experiments on the DIII-D tokamak demonstrate that edge localized mode (ELM) stability can be manipulated using localized electron cyclotron current drive (ECCD) in conjunction with resonant magnetic perturbations (RMPs). The injection of counter-plasma-current edge ECCD reduces the RMP amplitude required to suppress ELMs and bifurcates the pedestal into a high-confinement regime with 7 kPA pedestal pressure. This is the first time such a high confinement regime has been accessed through the bifurcation from the ballooning stability branch predicted by existing models. These observations are consistent with modeled ECCD manipulation of magnetic islands induced by the RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Profile observations of the Arctic atmospheric boundary layer with the BELUGA tethered balloon during MOSAiC

During the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, the Balloon-bornE moduLar Utility for profilinG the lower Atmosphere (BELUGA) was deployed from an ice floe drifting in the Fram Strait from 29 June to 27 July 2020. The BELUGA observations aimed to characterize the cloudy Arctic atmospheric boundary layer above the sea ice using a modular setup of five instrument packages. The in situ measurements included atmospheric thermodynamic and dynamic state parameters (air temperature, humidity, pressure, and three-dimensional wind), broadband solar and terrestrial irradiance, aerosol particle microphysical properties, and cloud particle images. In total, 66 profile observations were collected during 33 balloon flights from the surface to maximum altitudes of 0.3 to 1.5 km. The profiles feature a high vertical resolution of 0.01 m to 1 m, including measurements below, inside, and above frequently occurring low-level clouds. This publication describes the balloon operations, instruments, and the obtained data set. We invite the scientific community for joint analysis and model application of the freely available data on PANGAEA.

54 ENVIRONMENTAL SCIENCES↗